Identification of the Causative Bacteria in Musculoskeletal Infections Using 16s rDNA - Denaturing Gradient Gel Electrophoresis Analysis

Musculoskeletal infections are infections of the bone and surrounding tissues. They are currently diagnosed based on culture analysis, which is the gold standard for pathogen identification. However, these clinical laboratory methods are frequently inadequate for the identification of the causative agents, because a large percentage (25-50%) of confirmed musculoskeletal infections are false negatives in which no pathogen is identified in culture. My data supports these results. The goal of this project was to use PCR amplification of a portion of the 16S rRNA gene to test an alternative approach for the identification of these pathogens and to assess the diversity of the bacteria involved. The advantages of this alternative method are that it should increase sample sensitivity and the speed of detection. In addition, bacteria that are non-culturable or in low abundance can be detected using this molecular technique. However, a complication of this approach is that the majority of musculoskeletal infections are polymicrobial, which prohibits direct identification from the infected tissue by DNA sequencing of the initial 16S rDNA amplification products. One way to solve this problem is to use denaturing gradient gel electrophoresis (DGGE) to separate the PCR products before DNA sequencing. Denaturing gradient gel electrophoresis (DGGE) separates DNA molecules based on their melting point, which is determined by their DNA sequence. This analytical technique allows a mixture of PCR products of the same length that electrophoreses through agarose gels as one band, to be separated into different bands and then used for DNA sequence analysis. In this way, the DGGE allows for the identification of individual bacterial species in polymicrobial-infected tissue, which is critical for improving clinical outcomes. By combining the 16S rDNA amplification and the DGGE techniques together, an alternative approach for identification has been used. The 16S rRNA gene PCR-DGGE method includes several critical steps: DNA extraction from tissue biopsies, amplification of the bacterial DNA, PCR product separation by DGGE, amplification of the gel-extracted DNA, and DNA sequencing and analysis. Each step of the method was optimized to increase its sensitivity and for rapid detection of the bacteria present in human tissue samples. The limit of detection for the DNA extraction from tissue was at least 20 Staphylococcus aureus cells and the limit of detection for PCR was at least 0.05 pg of template DNA. The conditions for DGGE electrophoreses were optimized by using a double gradient of acrylamide (6 – 10%) and denaturant (30-70%), which increased the separation between distinct PCR products. The use of GelRed (Biotium) improved the DNA visualization in the DGGE gel. To recover the DNA from the DGGE gels the gel slices were excised, shredded in a bead beater, and the DNA was allowed to diffuse into sterile water overnight. The use of primers containing specific linkers allowed the entire amplified PCR product to be sequenced and then analyzed. The optimized 16S rRNA gene PCR-DGGE method was used to analyze 50 tissue biopsy samples chosen randomly from our collection. The results were compared to those of the Memorial Hermann Hospital Clinical Microbiology Laboratory for the same samples. The molecular method was congruent for 10 of the 17 (59%) culture negative tissue samples. In 7 of the 17 (41%) culture negative the molecular method identified a bacterium. The molecular method was congruent with the culture identification for 7 of the 33 (21%) positive cultured tissue samples. However, in 8 of the 33 (24%) the molecular method identified more organisms. In 13 of the 15 (87%) polymicrobial cultured tissue samples the molecular method identified at least one organism that was also identified by culture techniques. Overall, the DGGE analysis of 16S rDNA is an effective method to identify bacteria not identified by culture analysis.

Transcriptional regulation and functional analysis of the human Wilms' tumor 1 gene

The Wilms' tumor 1 gene (WT1) encodes a zinc-finger transcription factor and is expressed in urogenital, hematopoietic and other tissues. It is expressed in a temporal and spatial manner in both embryonic and adult stages. To obtain a better understanding of the biological function of WT1, we studied two aspects of WT1 regulation: one is the identification of tissue-specific cis-regulatory elements that regulate its expression, the other is the downstream genes which are modulated by WT1.^ My studies indicate that in addition to the promoter, other regulatory elements are required for the tissue specific expression of this gene. A 259-bp hematopoietic specific enhancer in intron 3 of the WT1 gene increased the transcriptional activity of the WT1 promoter by 8- to 10-fold in K562 and HL60 cells. Sequence analysis revealed both GATA and c-Myb motifs in the enhancer fragment. Mutation of the GATA motif decreased the enhancer activity by 60% in K562 cells. Electrophoretic mobility shift assays showed that both GATA-1 and GATA-2 proteins in K562 nuclear extracts bind to this motif. Cotransfection of the enhancer containing reporter construct with a GATA-1 or GATA-2 expression vector showed that both GATA-1 and GATA-2 transactivated this enhancer, increasing the CAT reporter activity 10-15 fold and 5-fold respectively. Similar analysis of the c-Myb motif by cotransfection with the enhancer CAT reporter construct and a c-Myb expression vector showed that c-Myb transactivated the enhancer by 5-fold. A DNase I-hypersensitive site has been identified in the 258 bp enhancer region. These data suggest that GATA-1 and c-Myb are responsible for the activity of this enhancer in hematopoietic cells and may bind to the enhancer in vivo. In the process of searching for cis-regulatory elements in transgenic mice, we have identified a 1.0 kb fragment that is 50 kb downstream from the promoter and is required for the central nervous system expression of WT1.^ In the search for downstream target genes of WT1, we noted that the proto-oncogene N-myc is coexpressed with the tumor suppressor gene WT1 in the developing kidney and is overexpressed in many Wilms' tumors. Sequence analysis revealed eleven consensus WT1 binding sites located in the 1 kb mouse N-myc promoter. We further showed that the N-myc promoter was down-regulated by WT1 in transient transfection assays. Electrophoretic mobility shift assays showed that oligonucleotides containing the WT1 motifs could bind WT1 protein. Furthermore, a Denys-Drash syndrome mutant of WT1, R394W, that has a mutation in the DNA binding domain, failed to repress the N-myc promoter. This suggests that the repression of the N-myc promoter is mediated by DNA binding of WT1. This finding helps to elucidate the relationship of WT1 and N-myc in tumorigenesis and renal development. ^

COMPARISON OF DNA TRANSFECTION AND MICROCELL-MEDIATED CHROMOSOME TRANSFER FOR DETECTING TRANSFORMING GENES IN HUMAN TUMOR CELL LINES

DNA mediated gene transfection is an important tool for moving and isolating genes from one cell type and putting them into a foreign genetic background. DNA transfection studies have been done routinely in many laboratories to identify and isolate transforming sequences in human tumors and tumor cell lines. A second technique, microcell-mediated chromosome transfer, allows the transfer of small numbers of intact human chromosome from one cell to another. This work was done to compare the efficiency of these two techniques in the transformation of NIH 3T3 mouse fibroblast cells.^ My intent in comparing these two techniques was to see if there was a difference in the transforming capability of DNA which has been purified of all associated protein and RNAs, and that of DNA which is introduced into a cell in its native form, the chromosome. If chromosomal sequences were capable of transforming the 3T3 cells in culture, the method could then be used as a way to isolate the relevant tumorigenic chromosomes from human tumors.^ The study shows, however, that even for those cell lines that contain transforming sequences identified by DNA-mediated gene transfer, those same sequences were unable to transform 3T3 cells when introduced to the cells by somatic fusion of human tumor microcells. I believe that the human transforming sequences in their original genetic conformation are not recognized by the mouse cell as genes which should be expressed; therefore, no noticeable transformation event was selected by this technique. ^

Roles of mismatch repair proteins in the recognition and processing of DNA interstrand crosslinks in human cells

DNA interstrand crosslinks (ICLs) are among the most toxic type of damage to a cell. Many ICL-inducing agents are widely used as therapeutic agents, e.g. cisplatin, psoralen. A bettor understanding of the cellular mechanism that eliminates ICLs is important for the improvement of human health. However, ICL repair is still poorly understood in mammals. Using a triplex-directed site-specific ICL model, we studied the roles of mismatch repair (MMR) proteins in ICL repair in human cells. We are also interested in using psoralen-conjugated triplex-forming oligonucleotides (TFOs) to direct ICLs to a specific site in targeted DNA and in the mammalian genomes. ^ MSH2 protein is the common subunit of two MMR recognition complexes, and MutSα and MutSβ. We showed that MSH2 deficiency renders human cell hypersensitive to psoralen ICLs. MMR recognition complexes bind specifically to triplex-directed psoralen ICLs in vitro. Together with the fact that psoralen ICL-induced repair synthesis is dramatically decreased in MSH2 deficient cell extracts, we demonstrated that MSH2 function is critical for the recognition and processing of psoralen ICLs in human cells. Interestingly, lack of MSH2 does not reduce the level of psoralen ICL-induced mutagenesis in human cells, suggesting that MSH2 does not contribute to error-generating repair of psoralen ICLs, and therefore, may represent a novel error-free mechanism for repairing ICLs. We also studied the role of MLH1, anther key protein in MMR, in the processing of psoralen ICLs. MLH1-deficient human cells are more resistant to psoralen plus UVA treatment. Importantly, MLH1 function is not required for the mutagenic repair of psoralen ICLs, suggesting that it is not involved in the error-generating repair of this type of DNA damage in human cells. ^ These are the first data indicating mismatch repair proteins may participate in a relatively error-free mechanism for processing psoralen ICL in human cells. Enhancement of MMR protein function relative to nucleotide excision repair proteins may reduce the mutagenesis caused by DNA ICLs in humans. ^ In order to specifically target ICLs to mammalian genes, we identified novel TFO target sequences in mouse and human genomes. Using this information, many critical mammalian genes can now be targeted by TFOs.^

Human DNA ligase and DNA polymerase as molecular targets for heavy metals and anticancer drugs

DNA ligase and DNA polymerase play important roles in DNA replication, repair, and recombination. Frequencies of spontaneous and chemical- and physical-induced mutations are correlated to the fidelity of DNA replication. This dissertation elucidates the mechanisms of the DNA ligation reaction by DNA ligases and demonstrates that human DNA ligase I and DNA polymerase $\alpha$ are the molecular targets for two metal ions, Zn$\sp{2+}$ and Cd$\sp{2+},$ and an anticancer drug, F-ara-ATP.^ Human DNA ligases were purified to homogeneity and their AMP binding domains were mapped. Although their AMP-binding domains are similar, there could be difference between the two ligases in their DNA binding domains.^ The formation of the AMP-DNA intermediate and the successive ligation reaction by human DNA ligases were analyzed. Both reactions showed their substrate specificity for ligases I and II, required Mg2+, and were inhibited by ATP.^ A protein inhibitor from HeLa cells and specific for human DNA ligase I but not ligase II and T4 ligase was discovered. It reversibly inhibited DNA ligation activity but not the AMP-binding activity due to the formation of a reversible ligase I-inhibitor complex.^ F-ara-ATP inhibited human DNA ligase I activity by competing with ATP for the AMP-binding site of DNA ligase I, forming a ligase I-F-ara-AMP complex, as well as when it was incorporated at 3$\sp\prime$-terminus of DNA nick by DNA polymerase $\alpha.$^ All steps of the DNA ligation reaction were inhibited by Zn$\sp{2+}$ and Cd$\sp{2+}$ in a concentration-dependent manner. Both ions did not show the ability to change the fidelity of DNA ligation reaction catalyzed by human DNA ligase I. However, Zn$\sp{2+}$ and Cd$\sp{2+}$ showed their contradictory effects on the fidelity of the reaction by human DNA polymerase $\alpha.$ Zn$\sp{2+}$ decreased the frequency of misinsertion but less affected that of mispair extension. On the contrary, Cd$\sp{2+}$ increased the frequencies of both misinsertion and mispair extension at very low concentration. Our data provided strong evidence in the molecular mechanisms for the mutagenicity of zinc and cadmium, and were comparable with the results previously reported. ^

Analysis of a human placental lactogen gene promoter

Human placental lactogen (hPL) and human growth hormone (hGH) comprise a multigene family that share $>$90% nucleic acid sequence homology including 500 bp of 5$\sp\prime$ flanking sequence. Despite these similarities, hGH is produced in the anterior pituitary while hPL is expressed in the placenta. For most genes studied to date, regulation of expression occurs by alterations at the level of transcriptional initiation. Nuclear proteins bind specific DNA sequences in the promoter to regulate gene expression. In this study, the hPL$\sb3$ promoter was analyzed for DNA sequences that contribute to its expression. The interaction between the hPL$\sb3$ promoter and nuclear proteins was examined using nuclear extracts from placental and non-placental cells.^ To identify regulatory elements in the promoter of the hPL$\sb3$ gene, 5$\sp\prime$ deletion mutants were constructed by cleaving 1200 bp of upstream sequence with various restriction enzymes. These DNA fragments were ligated 5$\sp\prime$ to a promoterless bacterial gene chloramphenicol acetyltransferase (CAT) and transfected into JEG-3 cells, a human placental choriocarcinoma cell line. The level of CAT activity reflects the ability of the promoter mutants to activate transcription. Deletion of the sequence between $-$142 bp and $-$129 bp, relative to the start of transcription, resulted in an 8-fold decrease in CAT activity. Nuclear proteins from JEG-3, HeLa, and HepG2 (human liver cells), formed specific binding complexes with this region of the hPL$\sb3$ promoter, as shown by gel mobility shift assay. The $-$142 bp to $-$129 bp region contains a sequence similar to that of a variant binding site for the transcription factor Sp1. Sp1-like proteins were identified by DNA binding assay, in the nuclear extracts of the three cell lines. A series of G nucleotides in the hPL$\sb3$ promoter regulatory region were identified by methylation interference assay to interact with the DNA-binding proteins and the pattern obtained is similar to that for other Sp1 binding sites that have been studied. This suggests that hPL$\sb3$ may be transcriptionally regulated by Sp1 or a Sp1-like transacting factor. ^

Mapping and cloning of genes from portions of the human genome using somatic cell hybrids

Human x rodent somatic cell hybrids have played an important role in human genetics research. They have been especially useful for assigning genes to chromosomes and isolating DNA markers from specific regions of the human genome.^ By employing a combination of somatic cell genetic, recombinant DNA, and cytogenetic techniques, human DNA excision repair gene ERCC4 was mapped regionally to human 16p13.13-13.2, even though the gene has not been cloned. Human x Chinese hamster ovary (CHO) cell hybrids selected for human ERCC4 activity and containing 16p13.1-p13.3 as the only human genetic material were identified. These hybrids were used to order DNA markers located in 16p13.1-p13.3. New DNA markers physically close to ERCC4 were isolated from such hybrids. Using amplified human DNA from the hybrids as probe in fluorescent in situ hybridization, the short arm breakpoint in the chromosome 16 inversion associated with acute myelomonocytic leukemia (AMML) was found to be physically close to the ERCC4 gene. The physical mapping and eventually, the cloning of the ERCC4 gene, will benefit the understanding of the DNA repair system and the study of other important biomedical problems such as tumorigenesis.^ To facilitate the cloning of ERCC4 gene and, in general, the cloning of genes from any defined regions of the human genome, a method was developed for the direct isolation of human transcribed genes ffom somatic cell hybrids. cDNA was prepared from human x rodent hybrid by using consensus 5$\sp\prime$ splice site sequences as primers. These primers were designed to select immature, unspliced messenger RNA (still retaining species specific repeat sequences) as templates. Screening of a derived cDNA library for human repeat sequences resulted in the isolation of human clones at the anticipated frequency with characteristics expected of exons of transcribed human genes. The usefulness of the splice site specific primers was analyzed and the cDNA synthesis conditions with these primers were optimized. The procedure was shown to be sensitive enough to clone weakly expressed genes. Studying the expression of the represented genes with the isolated clones was shown to be feasible. Such regional specific human gene fragments will be very valuable for many human genetic studies such as the search of inherited disease genes and the construction of a cDNA map of the human genome. ^

Characterization of the human tissue transglutaminase gene promoter

Transglutaminases are a family of calcium-dependent enzymes, that catalyze the covalent cross-linking of proteins by forming $\varepsilon(\gamma$-glutamyl)lysine isopeptide bonds. In order to investigate the molecular mechanisms regulating the expression of the tissue transglutaminase gene and to determine its biological functions, the goal of this research has been to clone and characterize the human tissue transglutaminase promoter. Thirteen clones of the tissue transglutaminase gene were obtained from the screening of a human placental genomic DNA library. A 1.74 Kb fragment derived from DNA located immediately upstream of the translation start site was subcloned and sequenced. Sequence analysis of this DNA fragment revealed that it contains a TATA box (TATAA), a CAAT box (GGACAAT), and a series of potential transcription factor binding sites and hormone response elements. Four regions of significant homology, a GC-rich region, a TG-rich region, an AG-rich region, and HR1, were identified by aligning 1.8 Kb of DNA flanking the human, mouse, and guinea pig tissue transglutaminase genes.^ To measure promoter activity, we subcloned the 1.74 Kb fragment of the tissue transglutaminase gene into a luciferase reporter vector to generate transglutaminase promoter/luciferase reporter constructs. Transfection experiments showed that this DNA segment includes a functional promoter with high constitutive activity. Deletion analysis revealed that the SP1 sites or corresponding sequences contribute to this activity. We investigated the role of DNA methylation in regulating the activity of the promoter and found that in vitro methylation of tissue transglutaminase promoter/luciferase reporter constructs suppressed their basal activity. Methylation of the promoter is inversely correlated with the expression of the tissue transglutaminase gene in vivo. These results suggest that DNA methylation may be one of the mechanisms regulating the expression of the gene. The tumor suppressor gene product p53 was also shown to inhibit the activity of the promoter, suggesting that induction of the tissue transglutaminase gene is not involved in the p53-dependent programmed cell death pathway. Although retinoids regulate the expression of the tissue transglutaminase gene in vivo, retinoid-inducible activity can not be identified in 3.7 Kb of DNA 5$\sp\prime$ to the tissue transglutaminase gene.^ The structure of the 5$\sp\prime$ end of the tissue transglutaminase gene was mapped. Alignment analysis of the human tissue transglutaminase gene with other human transglutaminases showed that tissue transglutaminase is the simplest member of transglutaminase superfamily. Transglutaminase genes show a conserved core of exons and introns but diverse N-terminuses and promoters. These observations suggest that key regulatory sequences and promoter elements have been appended upstream of the core transglutaminase gene to generate the diversity of regulated expression and regulated activity characteristic of the transglutaminase gene family. ^

The utility of N5 for gene therapy of human cancer

Cancer is a result of defects in the coordination of cell proliferation and programmed cell death. The extent of cell death is physiologically controlled by the activation of a programmed suicide pathway that results in a morphologically recognizable form of death termed apoptosis. Inducing apoptosis in tumor cells by gene therapy provides a potentially effective means to treat human cancers. The p84N5 is a novel nuclear death domain containing protein that has been shown to bind an amino terminal domain of retinoblastoma tumor suppressor gene product (pRb). Expression of N5 can induce apoptosis that is dependent upon its intact death domain and is inhibited by pRb. In many human cancer cells the functions of pRb are either lost through gene mutation or inactivated by different mechanisms. N5 based gene therapy may induce cell death preferentially in tumor cells relative to normal cells. We have demonstrated that N5 gene therapy is less toxic to normal cells than to tumor cells. To test the possibility that N5 could be used in gene therapy of cancer, we have generated a recombinant adenovirus engineered to express N5 and test the effects of viral infection on growth and tumorigenicity of human cancer cells. Adenovirus N5 infection significantly reduced the proliferation and tumorigenicity of breast, ovarian, and osteosarcoma tumor cell lines. Reduced proliferation and tumorigenicity were mediated by an induction of apoptosis as indicated by DNA fragmentation in infected cells. We also test the potential utility of N5 for gene therapy of pancreatic carcinoma that typically respond poorly to conventional treatment. Adenoviral mediated N5 gene transfer inhibits the growth of pancreatic cancer cell lines in vitro. N5 gene transfer also reduces the growth and metastasis of human pancreatic adenocarcinoma in subcutaneous and orthotopic mouse model. Interestingly, the pancreatic adenocarcinoma cells are more sensitive to N5 than they are to p53, suggesting that N5 gene therapy may be effective in tumors resistant to p53. We also test the possibilities of the use of N5 and p53 together on the inhibition of pancreatic cancer cell growth in vitro and vivo. Simultaneous use of N5 and RbΔCDK has been found to exert a greater extent on the inhibition of pancreatic cancer cell growth in vitro and in vivo. ^

Mechanism of damage sensing and cell killing following the inhibition of nucleotide excision repair by fludarabine in quiescent cells

Inhibition of DNA repair by the nucleoside of fludarabine (F-ara-A) induces toxicity in quiescent human cells. The sensing and signaling mechanisms following DNA repair inhibition by F-ara-A are unknown. The central hypothesis of this project was that the mechanistic interaction of a DNA repair initiating agent and a nucleoside analog initiates an apoptotic signal in quiescent cells. The purpose of this research was to identify the sensing and signaling mechanism(s) that respond to DNA repair inhibition by F-ara-A. Lymphocytes were treated with F-ara-A, to accumulate the active triphosphate metabolite and subsequently DNA repair was activated by UV irradiation. Pre-incubation of lymphocytes with 3 μM F-ara-A inhibited DNA repair initiated by 2 J/m2 UV and induced greater than additive apoptosis after 24 h. Blocking the incorporation of F-ara-A nucleotide into repairing DNA using 30 μM aphidicolin considerably lowered the apoptotic response. ^ Wild-type quiescent cells showed a significant loss in viability than did cells lacking functional sensor kinase DNA-PKcs or p53 as measured by colony formation assays. The functional status of ATM did not appear to affect the apoptotic outcome. Immunoprecipitation studies showed an interaction between the catalytic sub-unit of DNA-PK and p53 following DNA repair inhibition. Confocal fluorescence microscopy studies have indicated the localization pattern of p53, DNA-PK and γ-H2AX in the nucleus following DNA damage. Foci formation by γ-H2AX was seen as an early event that is followed by interaction with DNA-PKcs. p53 serine-15 phosphorylation and accumulation were detected 2 h after treatment. Fas/Fas ligand expression increased significantly after repair inhibition and was dependent on the functional status of p53. Blocking the interaction between Fas and Fas ligand by neutralizing antibodies significantly rescued the apoptotic fraction of cells. ^ Collectively, these results suggest that incorporation of the nucleoside analog into repair patches is critical for cytotoxicity and that the DNA damage, while being sensed by DNA-PK, may induce apoptosis by a p53-mediated signaling mechanism. Based on the results, a model is proposed for the sensing of F-ara-A-induced DNA damage that includes γ-H2AX, DNA-PKcs, and p53. Targeting the cellular DNA repair mechanism can be a potential means of producing cytotoxicity in a quiescent population of neoplastic cells. These results also provide mechanistic support for the success of nucleoside analogs with cyclophosphamide or other agents that initiate excision repair processes, in the clinic. ^

DNA damage-induced Fas ligand expression by epidermal dendritic cells mediates UV-induced immune suppression of contact hypersensitivity

Exposure to UVB radiation induces local and systemic immune suppression, evidenced by inhibition of the contact hypersensitivity response (CHS). Epidermal dendritic cells, the primary antigen presenting cells responsible for the induction of CHS, are profoundly altered in phenotype and function by UVB exposure and possess UV-specific DNA damage upon migrating to skin-draining lymph nodes. Expression of the proapoptotic protein FasL has been demonstrated in both skin and lymph node cells following UVB exposure. Additionally, functional FasL expression has recently been demonstrated to be required in the phenomenon of UV-induced immune suppression. To test the hypothesis that FasL expression by DNA-damaged Langerhans cells migrating to the skin-draining lymph nodes is a crucial event in the generation of this phenomenon, mice were given a single 5KJ/m2 UV-B exposure and sensitized to 0.5% FITC through the exposed area. Dendritic cells (DC) harvested from skin-draining lymph nodes (DLN) 18 hours following sensitization by magnetic CD11c-conjugated microbeads expressed high levels of Iab, CD80 and CD86, DEC-205 and bore the FITC hapten, suggesting epidermal origin. Radioimmunoassay of UV-specific DNA damage showed that DC contained the vast majority of cyclobutane pyrimidine dimers (CPDs) found in the DLN after UVB and exhibited increased FasL mRNA expression, a result which correlated with greatly increased FasL-mediated cytotoxicity. The ability of DCs to transfer sensitization to naïve hosts was lost following UVB exposure, a phenomenon which required DC FasL expression, and was completely reversed by cutaneous DNA repair. Collectively, these results demonstrate the central importance of DNA damage-induced FasL expression on migrating dendritic cells in mediating UV-induced suppression of contact hypersensitivity. ^

A comprehensive assessment of environmental ultraviolet radiation induced DNA damage in marine microorganisms

There is evidence that ultraviolet radiation (UVR) is increasing over certain locations on the Earth's surface. Of primary concern is the annual pattern of ozone depletion over Antarctica and the Southern Ocean. Reduction of ozone concentration selectively limits absorption of solar UV-B (290–320 nm), resulting in higher irradiance at the Earth's surface. The effects of ozone depletion on the human population and natural ecosystems, particularly the marine environment, are a matter of considerable concern. Indeed, marine plankton may serve as sensitive indicators of ozone depletion and UV-B fluctuations. Direct biological effects of UVR result from absorption of UV-B by DNA. Once absorbed, energy is dissipated by a variety of pathways, including covalent chemical reactions leading to the formation of photoproducts. The major types of photoproduct formed are cyclobutyl pyrimidine dimer (CPD) and pyrimidine(6-4)pyrimidone dimer [(6-4)PD]. Marine plankton repair these photoproducts using light-dependent photoenzymatic repair or nucleotide excision repair. The studies here show that fluctuations in CPD concentrations in the marine environment at Palmer Station, Antarctica correlate well with ozone concentration and UV-B irradiance at the Earth's surface. A comparison of photoproduct levels in marine plankton and DNA dosimeters show that bacterioplankton display higher resistance to solar UVR than phytoplankton in an ozone depleted environment. DNA damage in marine microorganisms was investigated during two separate latitudinal transects which covered a total range of 140°. We observed the same pattern of change in DNA damage levels in dosimeters and marine plankton as measured using two distinct quantitative techniques. Results from the transects show that differences in photosensitivity exist in marine plankton collected under varying UVR environments. Laboratory studies of Antarctic bacterial isolates confirm that marine bacterioplankton possess differences in survival, DNA damage induction, and repair following exposure to UVR. Results from DNA damage measurements during ozone season, along a latitudinal gradient, and in marine bacterial isolates suggest that changes in environmental UVR correlate with changes in UV-B induced DNA damage in marine microorganisms. Differences in the ability to tolerate UVR stress under different environmental conditions may determine the composition of the microbial communities inhabiting those environments. ^

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

Studies of the expression of DNA repair related genes

This dissertation examines the biological functions and the regulation of expression of DNA ligase I by studying its expression under different conditions.^ The gene expression of DNA ligase I was induced two- to four-fold in S-phase lymphoblastoid cells but was decreased to 15% of control after administration of a DNA damaging agent, 4-nitroquinoline-1-oxide. When cells were induced into differentiation, the expression level of DNA ligase I was decreased to less than 15% of that of the control cells. When the gene of DNA ligase I was examined for tissue specific expression in adult rats, high levels of DNA ligase I mRNA were observed in testis (8-fold), intermediate levels in ovary and brain (4-fold), and low levels were found in intestine, spleen, and liver (1- to 2-fold).^ In confluent cells of normal skin fibroblasts, UV irradiation induced the gene expression of DNA ligase I at 24 and 48 h. The induction of DNA ligase I gene expression requires active p53 protein. Introducing a vector containing the wild type p53 protein in the cells caused an induction of the DNA ligase I protein 24 h after the treatment.^ Our results indicate that, in addition to the regulation by phosphorylation/dephosphorylation, cellular DNA ligase I activity can be regulated at the gene transcription level, and the p53 tumor suppresser is one of the transcription factors for the DNA ligase I gene. Also, our results suggest that DNA ligase I is involved in DNA repair as well as in DNA replication.^ Also, as an early attempt to clone the human homolog of the yeast CDC9 gene which has been shown to be involved in DNA replication, DNA repair, and DNA recombination, we have identified a human gene with mRNA of 1.7 kb. This dissertation studies the gene regulation and the possible biological functions of this new human gene by examining its expression at different stages of the cell cycle, during cell differentiation, and in cellular response to DNA damage.^ The new gene that we recently identified from human cells is highly expressed in brain and reproductive organs (BRE). This BRE gene encodes an mRNA of 1.7-1.9 kb, with an open reading frame of 1,149 bp, and gives rise to a deduced polypeptide of 383 amino acid residues. No extensive homology was found between BRE and sequences from the EMBL-Gene Banks. BRE showed tissue-specific expression in adult rats. The steady state mRNA levels were high in testis (5-6 fold), ovary and brain (3-4 fold) compared to the spleen level, but low in intestine and liver (1-2 fold). The expression of this gene is responsive to DNA damage and/or retinoic acid (RA) treatment. Treatment of fibroblast cells with UV irradiation and 4-nitroquinoline-1-oxide caused more than 90% and 50% decreases in BRE mRNA, respectively. Similar decreases in BRE expression were observed after treatment of the brain glioma cell line U-251 and the promyelocytic cell line HL-60 with retinoic acid. (Abstract shortened by UMI). ^

RNA-SEQUENCING APPLICATIONS: GENE EXPRESSION QUANTIFICATION AND METHYLATOR PHENOTYPE IDENTIFICATION

My dissertation focuses on two aspects of RNA sequencing technology. The first is the methodology for modeling the overdispersion inherent in RNA-seq data for differential expression analysis. This aspect is addressed in three sections. The second aspect is the application of RNA-seq data to identify the CpG island methylator phenotype (CIMP) by integrating datasets of mRNA expression level and DNA methylation status. Section 1: The cost of DNA sequencing has reduced dramatically in the past decade. Consequently, genomic research increasingly depends on sequencing technology. However it remains elusive how the sequencing capacity influences the accuracy of mRNA expression measurement. We observe that accuracy improves along with the increasing sequencing depth. To model the overdispersion, we use the beta-binomial distribution with a new parameter indicating the dependency between overdispersion and sequencing depth. Our modified beta-binomial model performs better than the binomial or the pure beta-binomial model with a lower false discovery rate. Section 2: Although a number of methods have been proposed in order to accurately analyze differential RNA expression on the gene level, modeling on the base pair level is required. Here, we find that the overdispersion rate decreases as the sequencing depth increases on the base pair level. Also, we propose four models and compare them with each other. As expected, our beta binomial model with a dynamic overdispersion rate is shown to be superior. Section 3: We investigate biases in RNA-seq by exploring the measurement of the external control, spike-in RNA. This study is based on two datasets with spike-in controls obtained from a recent study. We observe an undiscovered bias in the measurement of the spike-in transcripts that arises from the influence of the sample transcripts in RNA-seq. Also, we find that this influence is related to the local sequence of the random hexamer that is used in priming. We suggest a model of the inequality between samples and to correct this type of bias. Section 4: The expression of a gene can be turned off when its promoter is highly methylated. Several studies have reported that a clear threshold effect exists in gene silencing that is mediated by DNA methylation. It is reasonable to assume the thresholds are specific for each gene. It is also intriguing to investigate genes that are largely controlled by DNA methylation. These genes are called “L-shaped” genes. We develop a method to determine the DNA methylation threshold and identify a new CIMP of BRCA. In conclusion, we provide a detailed understanding of the relationship between the overdispersion rate and sequencing depth. And we reveal a new bias in RNA-seq and provide a detailed understanding of the relationship between this new bias and the local sequence. Also we develop a powerful method to dichotomize methylation status and consequently we identify a new CIMP of breast cancer with a distinct classification of molecular characteristics and clinical features.