Structural origins of the exonuclease resistance of a zwitterionic RNA

Nuclease resistance and RNA affinity are key criteria in the search for optimal antisense nucleic acid modifications, but the origins of the various levels of resistance to nuclease degradation conferred by chemical modification of DNA and RNA are currently not understood. The 2′-O-aminopropyl (AP)-RNA modification displays the highest nuclease resistance among all phosphodiester-based analogues and its RNA binding affinity surpasses that of phosphorothioate DNA by 1°C per modified residue. We found that oligodeoxynucleotides containing AP-RNA residues at their 3′ ends competitively inhibit the degradation of single-stranded DNA by the Escherichia coli Klenow fragment (KF) 3′-5′ exonuclease and snake venom phosphodiesterase. To shed light on the origins of nuclease resistance brought about by the AP modification, we determined the crystal structure of an A-form DNA duplex with AP-RNA modifications at 1.6-Å resolution. In addition, the crystal structures of complexes between short DNA fragments carrying AP-RNA modifications and wild-type KF were determined at resolutions between 2.2 and 3.0 Å and compared with the structure of the complex between oligo(dT) and the D355A/E357A KF mutant. The structural models suggest that interference of the positively charged 2′-O-substituent with the metal ion binding site B of the exonuclease allows AP-RNA to effectively slow down degradation.

Comparative analyses of the secondary structures of synthetic and intracellular yeast MFA2 mRNAs

The overall folded (global) structure of mRNA may be critical to translation and turnover control mechanisms, but it has received little experimental attention. Presented here is a comparative analysis of the basic features of the global secondary structure of a synthetic mRNA and the same intracellular eukaryotic mRNA by dimethyl sulfate (DMS) structure probing. Synthetic MFA2 mRNA of Saccharomyces cerevisiae first was examined by using both enzymes and chemical reagents to determine single-stranded and hybridized regions; RNAs with and without a poly(A) tail were compared. A folding pattern was obtained with the aid of the mfold program package that identified the model that best satisfied the probing data. A long-range structural interaction involving the 5′ and 3′ untranslated regions and causing a juxtaposition of the ends of the RNA, was examined further by a useful technique involving oligo(dT)-cellulose chromatography and antisense oligonucleotides. DMS chemical probing of A and C nucleotides of intracellular MFA2 mRNA was then done. The modification data support a very similar intracellular structure. When low reactivity of A and C residues is found in the synthetic RNA, ≈70% of the same sites are relatively more resistant to DMS modification in vivo. A slightly higher sensitivity to DMS is found in vivo for some of the A and C nucleotides predicted to be hybridized from the synthetic structural model. With this small mRNA, the translation process and mRNA-binding proteins do not block DMS modifications, and all A and C nucleotides are modified the same or more strongly than with the synthetic RNA.

Cloning of chlorophyllase, the key enzyme in chlorophyll degradation: Finding of a lipase motif and the induction by methyl jasmonate

Chlorophyllase (Chlase) is the first enzyme involved in chlorophyll (Chl) degradation and catalyzes the hydrolysis of ester bond to yield chlorophyllide and phytol. In the present study, we isolated the Chlase cDNA. We synthesized degenerate oligo DNA probes based on the internal amino acid sequences of purified Chlase from Chenopodium album, screened the C. album cDNA library, and cloned a cDNA (CaCLH, C. album chlorophyll-chlorophyllido hydrolase). The deduced amino acid sequence (347 aa residues) had a lipase motif overlapping with an ATP/GTP-binding motif (P-loop). CaCLH possibly was localized in the extraplastidic part of the cell, because a putative signal sequence for endoplasmic reticulum is at the N terminus. The amino acid sequence shared 37% identity with a function-unknown gene whose mRNA is inducible by coronatine and methyl jasmonate (MeJA) in Arabidopsis thaliana (AtCLH1). We expressed the gene products of AtCLH1 and of CaCLH in Escherichia coli, and they similarly exhibited Chlase activity. Moreover, we isolated another full-length cDNA based on an Arabidopsis genomic fragment and expressed it in E. coli, demonstrating the presence of the second Arabidopsis CLH gene (AtCLH2). No typical feature of signal sequence was identified in AtCLH1, whereas AtCLH2 had a typical signal sequence for chloroplast. AtCLH1 mRNA was induced rapidly by a treatment of MeJA, which is known to promote senescence and Chl degradation in plants, and a high mRNA level was maintained up to 9 h. AtCLH2, however, did not respond to MeJA.

Oligomerization of activated d- and l-guanosine mononucleotides on templates containing d- and l-deoxycytidylate residues

The oligomerization of activated d- and l- and racemic guanosine-5′-phosphoro-2-methylimidazole on short templates containing d- and l-deoxycytidylate has been studied. Results obtained with d-oligo(dC)s as templates are similar to those previously reported for experiments with a poly(C) template. When one l-dC or two consecutive l-dCs are introduced into a d-template, regiospecific synthesis of 3′-5′ oligo(G)s proceeds to the end of the template, but three consecutive l-dCs block synthesis. Alternating d-,l-oligomers do not facilitate oligomerization of the d-, l-, and racemic 2-guanosine-5′-phosphoro-2-methylimidazole. We suggest that once a “predominately d-metabolism” existed, occasional l-residues in a template would not have led to the termination of self-replication.

Comprehensive Identification of Cell Cycle–regulated Genes of the Yeast Saccharomyces cerevisiae by Microarray HybridizationD⃞

We sought to create a comprehensive catalog of yeast genes whose transcript levels vary periodically within the cell cycle. To this end, we used DNA microarrays and samples from yeast cultures synchronized by three independent methods: α factor arrest, elutriation, and arrest of a cdc15 temperature-sensitive mutant. Using periodicity and correlation algorithms, we identified 800 genes that meet an objective minimum criterion for cell cycle regulation. In separate experiments, designed to examine the effects of inducing either the G1 cyclin Cln3p or the B-type cyclin Clb2p, we found that the mRNA levels of more than half of these 800 genes respond to one or both of these cyclins. Furthermore, we analyzed our set of cell cycle–regulated genes for known and new promoter elements and show that several known elements (or variations thereof) contain information predictive of cell cycle regulation. A full description and complete data sets are available at http://cellcycle-www.stanford.edu

In vitro cloning of complex mixtures of DNA on microbeads: Physical separation of differentially expressed cDNAs

We describe a method for cloning nucleic acid molecules onto the surfaces of 5-μm microbeads rather than in biological hosts. A unique tag sequence is attached to each molecule, and the tagged library is amplified. Unique tagging of the molecules is achieved by sampling a small fraction (1%) of a very large repertoire of tag sequences. The resulting library is hybridized to microbeads that each carry ≈106 strands complementary to one of the tags. About 105 copies of each molecule are collected on each microbead. Because such clones are segregated on microbeads, they can be operated on simultaneously and then assayed separately. To demonstrate the utility of this approach, we show how to label and extract microbeads bearing clones differentially expressed between two libraries by using a fluorescence-activated cell sorter (FACS). Because no prior information about the cloned molecules is required, this process is obviously useful where sequence databases are incomplete or nonexistent. More importantly, the process also permits the isolation of clones that are expressed only in given tissues or that are differentially expressed between normal and diseased states. Such clones then may be spotted on much more cost-effective, tissue- or disease-directed, low-density planar microarrays.

Generation of longer cDNA fragments from serial analysis of gene expression tags for gene identification

We have developed a technique called the generation of longer cDNA fragments from serial analysis of gene expression (SAGE) tags for gene identification (GLGI), to convert SAGE tags of 10 bases into their corresponding 3′ cDNA fragments covering hundred bases. A primer containing the 10-base SAGE tag is used as the sense primer, and a single base anchored oligo(dT) primer is used as an antisense primer in PCR, together with Pfu DNA polymerase. By using this approach, a cDNA fragment extending from the SAGE tag toward the 3′ end of the corresponding sequence can be generated. Application of the GLGI technique can solve two critical issues in applying the SAGE technique: one is that a longer fragment corresponding to a SAGE tag, which has no match in databases, can be generated for further studies; the other is that the specific fragment corresponding to a SAGE tag can be identified from multiple sequences that match the same SAGE tag. The development of the GLGI method provides several potential applications. First, it provides a strategy for even wider application of the SAGE technique for quantitative analysis of global gene expression. Second, a combined application of SAGE/GLGI can be used to complete the catalogue of the expressed genes in human and in other eukaryotic species. Third, it can be used to identify the 3′ cDNA sequence from any exon within a gene. It can also be used to confirm the reality of exons predicted by bioinformatic tools in genomic sequences. Fourth, a combined application of SAGE/GLGI can be applied to define the 3′ boundary of expressed genes in the genomic sequences in human and in other eukaryotic genomes.

Intranuclear diffusion and hybridization state of oligonucleotides measured by fluorescence correlation spectroscopy in living cells

Fluorescein-labeled oligodeoxynucleotides (oligos) were introduced into cultured rat myoblasts, and their molecular movements inside the nucleus were studied by fluorescence correlation spectroscopy (FCS) and fluorescence recovery after photobleaching (FRAP). FCS revealed that a large fraction of both intranuclear oligo(dT) (43%) and oligo(dA) (77%) moves rapidly with a diffusion coefficient of 4 × 10−7 cm2/s. Interestingly, this rate of intranuclear oligo movement is similar to their diffusion rates measured in aqueous solution. In addition, we detected a large fraction (45%) of the intranuclear oligo(dT), but not oligo(dA), diffusing at slower rates (≤1 × 10−7 cm2/s). The amount of this slower-moving oligo(dT) was greatly reduced if the oligo(dT) was prehybridized in solution with (unlabeled) oligo(dA) prior to introduction to cells, presumably because the oligo(dT) was then unavailable for subsequent hybridization to endogenous poly(A) RNA. The FCS-measured diffusion rate for much of the slower oligo(dT) population approximated the diffusion rate in aqueous solution of oligo(dT) hybridized to a large polyadenylated RNA (1.0 × 10−7 cm2/s). Moreover, this intranuclear movement rate falls within the range of calculated diffusion rates for an average-sized heterogeneous nuclear ribonucleoprotein particle in aqueous solution. A subfraction of oligo(dT) (15%) moved over 10-fold more slowly, suggesting it was bound to very large macromolecular complexes. Average diffusion coefficients obtained from FRAP experiments were in agreement with the FCS data. These results demonstrate that oligos can move about within the nucleus at rates comparable to those in aqueous solution and further suggest that this is true for large ribonucleoprotein complexes as well.

Importance of replication in microarray gene expression studies: Statistical methods and evidence from repetitive cDNA hybridizations

We present statistical methods for analyzing replicated cDNA microarray expression data and report the results of a controlled experiment. The study was conducted to investigate inherent variability in gene expression data and the extent to which replication in an experiment produces more consistent and reliable findings. We introduce a statistical model to describe the probability that mRNA is contained in the target sample tissue, converted to probe, and ultimately detected on the slide. We also introduce a method to analyze the combined data from all replicates. Of the 288 genes considered in this controlled experiment, 32 would be expected to produce strong hybridization signals because of the known presence of repetitive sequences within them. Results based on individual replicates, however, show that there are 55, 36, and 58 highly expressed genes in replicates 1, 2, and 3, respectively. On the other hand, an analysis by using the combined data from all 3 replicates reveals that only 2 of the 288 genes are incorrectly classified as expressed. Our experiment shows that any single microarray output is subject to substantial variability. By pooling data from replicates, we can provide a more reliable analysis of gene expression data. Therefore, we conclude that designing experiments with replications will greatly reduce misclassification rates. We recommend that at least three replicates be used in designing experiments by using cDNA microarrays, particularly when gene expression data from single specimens are being analyzed.

Ongoing immunoglobulin somatic mutation in germinal center B cell-like but not in activated B cell-like diffuse large cell lymphomas

B cell diffuse large cell lymphoma (B-DLCL) is a heterogeneous group of tumors, based on significant variations in morphology, clinical presentation, and response to treatment. Gene expression profiling has revealed two distinct tumor subtypes of B-DLCL: germinal center B cell-like DLCL and activated B cell-like DLCL. In a separate study, we determined that B-DLCL can also be subdivided into two groups based on the presence or absence of ongoing Ig gene hypermutation. Here, we evaluated the correlation between these B-DLCL subtypes established by the two different methods. Fourteen primary B-DLCL cases were studied by gene expression profiling using DNA microarrays and for the presence of ongoing mutations in their Ig heavy chain gene. All seven cases classified as germinal center B cell-like DLCL by gene expression showed the presence of ongoing mutations in the Ig genes. Five of the seven cases classified by gene expression as activated B cell-like DLCL had no ongoing somatic mutations, whereas, in the remaining two cases, a single point mutation was observed in only 2 of 15 and 21 examined molecular clones of variable heavy (VH) chain gene, respectively. These two cases were distantly related to the rest of the activated B cell-like DLCL tumors by gene expression. Our findings validate the concept that lymphoid malignancies are derived from cells at discrete stages of normal lymphocyte maturation and that the malignant cells retain the genetic program of those normal cells.

Hypoxia-induced gene expression profiling in the euryoxic fish Gillichthys mirabilis

Hypoxia is important in both biomedical and environmental contexts and necessitates rapid adaptive changes in metabolic organization. Mammals, as air breathers, have a limited capacity to withstand sustained exposure to hypoxia. By contrast, some aquatic animals, such as certain fishes, are routinely exposed and resistant to severe environmental hypoxia. Understanding the changes in gene expression in fishes exposed to hypoxic stress could reveal novel mechanisms of tolerance that may shed new light on hypoxia and ischemia in higher vertebrates. Using cDNA microarrays, we have studied gene expression in a hypoxia-tolerant burrow-dwelling goby fish, Gillichthys mirabilis. We show that a coherent picture of a complex transcriptional response can be generated for a nonmodel organism for which sequence data were unavailable. We demonstrate that: (i) although certain shifts in gene expression mirror changes in mammals, novel genes are differentially expressed in fish; and (ii) tissue-specific patterns of expression reflect the different metabolic roles of tissues during hypoxia.

Immobilization of oligodeoxyribonucleotides with multiple anchors to microchips

Facile modification of oligodeoxyribonucleotides is required for efficient immobilization to a pre-activated glass surface. This report presents an oligodeoxyribonucleotide which contains a hairpin stem–loop structure with multiple phosphorothioate moieties in the loop. These moieties are used to anchor the oligo to glass slides that are pre-activated with bromoacetamidopropylsilane. The efficiency of the attachment reaction was improved by increasing the number of phosphorothioates in the loop, as shown in the remarkable enhancement of template hybridization and single base extension through catalysis by DNA polymerase. The loop and stem presumably serve as lateral spacers between neighboring oligodeoxyribonucleotides and as a linker arm between the glass surface and the single-stranded sequence of interest. The oligodeoxyribonucleotides of this hairpin stem–loop architecture with multiple phosphorothioate moieties have broad application in DNA chip-based gene analysis.

Stimulation of human 8-oxoguanine-DNA glycosylase by AP-endonuclease: potential coordination of the initial steps in base excision repair

8-Oxoguanine-DNA glycosylase 1 (OGG1), with intrinsic AP lyase activity, is the major enzyme for repairing 7,8-dihydro-8-oxoguanine (8-oxoG), a critical mutagenic DNA lesion induced by reactive oxygen species. Human OGG1 excised the damaged base from an 8-oxoG·C-containing duplex oligo with a very low apparent kcat of 0.1 min–1 at 37°C and cleaved abasic (AP) sites at half the rate, thus leaving abasic sites as the major product. Excision of 8-oxoG by OGG1 alone did not follow Michaelis–Menten kinetics. However, in the presence of a comparable amount of human AP endonuclease (APE1) the specific activity of OGG1 was increased ∼5-fold and Michaelis–Menten kinetics were observed. Inactive APE1, at a higher molar ratio, and a bacterial APE (Nfo) similarly enhanced OGG1 activity. The affinity of OGG1 for its product AP·C pair (Kd ∼ 2.8 nM) was substantially higher than for its substrate 8-oxoG·C pair (Kd ∼ 23.4 nM) and the affinity for its final β-elimination product was much lower (Kd ∼ 233 nM). These data, as well as single burst kinetics studies, indicate that the enzyme remains tightly bound to its AP product following base excision and that APE1 prevents its reassociation with its product, thus enhancing OGG1 turnover. These results suggest coordinated functions of OGG1 and APE1, and possibly other enzymes, in the DNA base excision repair pathway.

The transcriptional program of a human B cell line in response to Myc

The proto-oncogene c-myc (myc) encodes a transcription factor (Myc) that promotes growth, proliferation and apoptosis. Myc has been suggested to induce these effects by induction/repression of downstream genes. Here we report the identification of potential Myc target genes in a human B cell line that grows and proliferates depending on conditional myc expression. Oligonucleotide microarrays were applied to identify downstream genes of Myc at the level of cytoplasmic mRNA. In addition, we identified potential Myc target genes in nuclear run-on experiments by changes in their transcription rate. The identified genes belong to gene classes whose products are involved in amino acid/protein synthesis, lipid metabolism, protein turnover/folding, nucleotide/DNA synthesis, transport, nucleolus function/RNA binding, transcription and splicing, oxidative stress and signal transduction. The identified targets support our current view that myc acts as a master gene for growth control and increases transcription of a large variety of genes.

Effect of 2′-O-methyl antisense ORNs on expression of thymidylate synthase in human colon cancer RKO cells

Translation of thymidylate synthase (TS) mRNA is controlled by its own protein end-product TS in a negative autoregulatory manner. Disruption of this regulation results in increased synthesis of TS and may lead to the development of cellular drug resistance to TS-directed anticancer agents. As a strategy to inhibit TS expression, antisense 2′-O-methyl RNA oligoribonucleotides (ORNs) were designed to directly target the 5′ upstream cis-acting regulatory element (nucleotides 80–109) of TS mRNA. A 30 nt ORN, HYB0432, inhibited TS expression in human colon cancer RKO cells in a dose-dependent manner but had no effect on the expression of β-actin, α-tubulin or topoisomerase I. TS expression was unaffected by treatment with control sense or mismatched ORNs. HYB0504, an 18 nt ORN targeting the same core sequence, also repressed expression of TS protein. However, further reduction in oligo size resulted in loss of antisense activity. Following HYB0432 treatment, TS protein levels were reduced by 60% within 6 h and were maximally reduced by 24 h. Expression of p53 protein was inversely related to that of TS, suggesting that p53 expression may be directly linked to intracellular levels of TS. Northern blot analysis demonstrated that TS mRNA was unaffected by HYB0432 treatment. The half-life of TS protein was unchanged after antisense treatment suggesting that the mechanism of action of antisense ORNs is mediated through a process of translational arrest. These findings demonstrate that an antisense ORN targeted at a critical cis-acting element on TS mRNA can specifically inhibit expression of TS protein in RKO cells.