RNA extraction from ten year old formalin-fixed paraffin-embedded breast cancer samples: a comparison of column purification and magnetic bead-based technologies

Abstract Background The development of protocols for RNA extraction from paraffin-embedded samples facilitates gene expression studies on archival samples with known clinical outcome. Older samples are particularly valuable because they are associated with longer clinical follow up. RNA extracted from formalin-fixed paraffin-embedded (FFPE) tissue is problematic due to chemical modifications and continued degradation over time. We compared quantity and quality of RNA extracted by four different protocols from 14 ten year old and 14 recently archived (three to ten months old) FFPE breast cancer tissues. Using three spin column purification-based protocols and one magnetic bead-based protocol, total RNA was extracted in triplicate, generating 336 RNA extraction experiments. RNA fragment size was assayed by reverse transcription-polymerase chain reaction (RT-PCR) for the housekeeping gene glucose-6-phosphate dehydrogenase (G6PD), testing primer sets designed to target RNA fragment sizes of 67 bp, 151 bp, and 242 bp. Results Biologically useful RNA (minimum RNA integrity number, RIN, 1.4) was extracted in at least one of three attempts of each protocol in 86–100% of older and 100% of recently archived ("months old") samples. Short RNA fragments up to 151 bp were assayable by RT-PCR for G6PD in all ten year old and months old tissues tested, but none of the ten year old and only 43% of months old samples showed amplification if the targeted fragment was 242 bp. Conclusion All protocols extracted RNA from ten year old FFPE samples with a minimum RIN of 1.4. Gene expression of G6PD could be measured in all samples, old and recent, using RT-PCR primers designed for RNA fragments up to 151 bp. RNA quality from ten year old FFPE samples was similar to that extracted from months old samples, but quantity and success rate were generally higher for the months old group. We preferred the magnetic bead-based protocol because of its speed and higher quantity of extracted RNA, although it produced similar quality RNA to other protocols. If a chosen protocol fails to extract biologically useful RNA from a given sample in a first attempt, another attempt and then another protocol should be tried before excluding the case from molecular analysis.

Quantification of oxidized levels of specific RNA species using an aldehyde reactive probe

Emerging evidence has shown that oxidation of RNA, including messenger RNA (mRNA), is elevated in several age-related diseases, although investigation of oxidized levels of individual RNA species has been limited. Recently we reported that an aldehyde reactive probe (ARP) quantitatively reacts with oxidatively modified depurinated/depyrimidinated (abasic) RNA. Here we report a novel method to isolate oxidized RNA using ARP and streptavidin beads. An oligo RNA containing abasic sites that were derivatized with ARP was pulled down by streptavidin beads, whereas a control oligo RNA was not. In vitro oxidized RNA, as well as total cellular RNA, isolated from oxidatively stressed cells was also pulled down, dependent on oxidation level, and concentrated in the pull-down fraction. Quantitative reverse transcription polymerase chain reaction (RT-PCR) using RNA in the pull-down fraction demonstrated that several gene transcripts were uniquely increased in the fraction by oxidative stress. Thus, our method selectively concentrates oxidized RNA by pull-down and enables the assessment of oxidation levels of individual RNA species. (C) 2011 Elsevier Inc. All rights reserved.

Comparison of amplification methods for transcriptomic analyses of low abundance prokaryotic RNA sources

Microarrays have established as instrumental for bacterial detection, identification, and genotyping as well as for transcriptomic studies. For gene expression analyses using limited numbers of bacteria (derived from in vivo or ex vivo origin, for example), RNA amplification is often required prior to labeling and hybridization onto microarrays. Evaluation of the fidelity of the amplification methods is crucial for the robustness and reproducibility of microarray results. We report here the first utilization of random primers and the highly processive Phi29 phage polymerase to amplify material for transcription profiling analyses. We compared two commercial amplification methods (GenomiPhi and MessageAmp kits) with direct reverse-transcription as the reference method, focusing on the robustness of mRNA quantification using either microarrays or quantitative RT-PCR. Both amplification methods using either poly-A tailing followed by in vitro transcription, or direct strand displacement polymerase, showed appreciable linearity. Strand displacement technique was particularly affordable compared to in vitro transcription-based (IVT) amplification methods and consisted in a single tube reaction leading to high amplification yields. Real-time measurements using low-, medium-, and highly expressed genes revealed that this simple method provided linear amplification with equivalent results in terms of relative messenger abundance as those obtained by conventional direct reverse-transcription.

Analysis of primary Wilms tumors for WT1 RNA processing alterations

Wilms tumor (WT) is an embryonal renal tumor with a heterogeneous genetic etiology that serves as a valuable model for studying tumorigenesis. Biallelic inactivation of the tumor suppressor gene WT1, a zinc-finger transcriptional regulator located at 11p13, is critical for the development of some Wilms tumors. Interestingly, WT1 genomic analysis has demonstrated mutations in less than 20% of WT cases. This suggests either other genes play a more major role in Wilms tumorigenesis or WT1 is functionally altered by mechanisms other than DNA mutation. Previous observations in rat and in WT xenograft cell lines have suggested that abnormal WT1 RNA processing (exon 6 RNA editing and aberrant exon 2 splicing, respectively) is a potential mechanism of altering WT1 function in the absence of a WT1 DNA mutation. However, the role of this abnormal RNA processing has not previously been assessed in primary Wilms tumors. ^ To test the hypothesis that abnormal WT1 RNA processing is a mechanism of WT1alteration during tumor development, WT1 RNA from 85 primary tumors was analyzed using reverse transcription and polymerase chain reaction amplification (RT-PCR). Although no evidence for WT1 RNA editing was observed, variable levels (5% to 50%) of aberrant WT1 exon 2 splicing were detected for 11 tumors in the absence of a detectable WT1 DNA mutation. Also, alteration of normal WT1 alternative splicing, observed as RNA isoform loss, was detected in five tumors with no apparent WT1 genomic alteration, although no consistent pattern of RNA isoform loss was detected. This abnormal WT1 splicing, detected by either loss of exon 2 from some of the transcripts or loss of RNA isoforms, is statistically correlated with relapse (p = 0.005). These studies demonstrate that abnormal WT1 RNA processing is not a common mechanism of abrogating normal WT1 function in primary tumors. However, in those cases in which abnormal WTI splicing is present, these data indicate that it may serve as a useful prognostic marker for relapse in WT patients. ^

Cloning of a protease gene family of Fasciola hepatica by the polymerase chain reaction

Degenerate oligonucleotide primers derived from conserved cysteine protease sequences were used in the reverse transcription polymerase chain reaction to amplify seven different cysteine protease cDNA clones, Fcp1-7, from RNA isolated from adult Fasciola hepatica. Five of the amplified F. hepatica sequences showed homology to the cathepsin L type and two were more related to the cathepsin B type. Southern blot analysis suggests that some members of this protease gene family are present in multiple copies. Northern blot analysis revealed differences in the levels of steady state mRNA expression for some of these proteases. The 5' and the 3' regions of Fcp1 were amplified using the rapid amplification of cDNA ends PCR protocol (RACE-PCR) and an additional clone was obtained by screening a lambda gt10 cDNA library using Fcp1 as a probe. The Fcp1 cDNA fragment was also subcloned in the expression vector pGEX and expressed as a glutathione-S-transferase (GST) fusion protein in Escherichia coli. Antibodies, raised in rabbits against the GST:Fcp1 fusion protein, were used in western blot analysis to examine expression in different life-cycle stages of F. hepatica. In extracts from adult and immature parasites, the immune serum recognised predominantly two proteins of 30 kDa and 38 kDa. In other parasite stages, proteins of different molecular weight were recognised by the anti-GST:Fcp1 antiserum, indicating stage-specific gene expression or processing of Fcp1. In gelatine substrate gel analysis, strong proteolytic activity could be detected at 30 kDa, but not at 38 kDa, suggesting that the 30 kDa protein represents the mature enzyme and the 38 kDa protein the proenzyme.

Cloning of a full-length cDNA encoding bovine interleukin 4 by the polymerase chain reaction

A human interleukin 4 (hIL-4)-encoding cDNA (hIL4) probe was used to screen a bovine genomic library, and three clones containing sequences with homology to the human and mouse IL4 cDNAs were isolated. Sequence information obtained from one of these genomic clones was used to design an oligodeoxyribonucleotide primer corresponding to the transcription start point region for use in the polymerase chain reaction (PCR). The PCR-RACE protocol, designed for the rapid amplification of cDNA ends, was successfully used to generate a full-length bovine IL4 (bIL4) cDNA clone from polyadenylated RNA isolated from concanavalin A-stimulated bovine lymph node cells. The bIL4 cDNA is 570 bp in length and contains an open reading frame of 405 nucleotides (nt), coding for a 15.1-kDa precursor of 135 amino acids (aa), which should be reduced to 12.6 kDa for unglycosylated bIL4 after cleavage of a putative hydrophobic leader sequence of 24 aa. The aa sequence contains one possible Asn-linked glycosylation site. Bovine IL4 is shorter than mouse (mIL4) and hIL4, because of a 51-nt deletion in the coding region. Comparison of the overall nt and deduced aa sequences shows a greater homology of bIL4 with hIL4 than with mIL4. This homology is not evenly distributed, however, with the nt sequences 5' and 3' of the coding region showing a much greater homology between all three species than the coding sequence.

RNA binding characteristics and overall topology of the vaccinia virus polyadenylation complex

mRNA 3′ polyadenylation is central to mRNA biogenesis in prokaryotes and eukaryotes, and is implicated in numerous aspects of mRNA metabolism, including efficiency of mRNA export from the nucleus, message stability, and initiation of translation. However, due to the great complexity of the eukaryotic polyadenylation apparatus, the mechanisms of RNA 3 ′ end processing have remained elusive. Although the RNA processing reactions leading to polyadenylated messenger RNA have been studied in many systems, and much progress has been made, a complete understanding of the biochemistry of the poly(A) polymerase enzyme is still lacking. My research uses Vaccinia virus as a model system to gain a better understanding of this complicated polyadenylation process, which consist of RNA binding, catalysis and polymerase translocation. ^ Vaccinia virus replicates in the cytoplasm of its host cell, so it must employ its own poly(A) polymerase (PAP), a heterodimer of two virus encoded proteins, VP55 and VP39. VP55 is the catalytic subunit, adding 30 adenylates to a non-polyadenylated RNA in a rapid processive manner before abruptly changing to a slow, non-processive mode of adenylate addition and dissociating from the RNA. VP39 is the stimulatory subunit. It has no polyadenylation catalytic activity by itself, but when associated with VP55 it facilitates the semi-processive synthesis of tails several hundred adenylates in length. ^ Oligonucleotide selection and competition studies have shown that the heterodimer binds a minimal motif of (rU)2 (N)25 U, the “heterodimer binding motif”, within an oligonucleotide, and its primer selection for polyadenylation is base-type specific. ^ Crosslinking studies using photosensitive uridylate analogs show that within a VP55-VP39-primer ternary complex, VP55 comes into contact with all three required uridylates, while VP39 only contacts the downstream uridylate. Further studies, using a backbone-anchored photosensitive crosslinker show that both PAP subunits are in close proximity to the downstream −10 to −21 region of 50mer model primers containing the heterodimer binding motif. This equal crosslinking to both subunits suggests that the dimerization of VP55 and VP39 creates either a cleft or a channel between the two subunits through which this region of RNA passes. ^ Peptide mapping studies of VP39 covalently crosslinked to the oligonucleotide have identified residue R107 as the amino acid in close proximity to the −10 uridylate. This helps us project a conceptual model onto the known physical surface of this subunit. In the absence of any tertiary structural data for VP55, we have used a series of oligonucleotide selection assays, as well as crosslinking, nucleotide transfer assays, and gel shift assays to gain insight into the requirements for binding, polyadenylation and translocation. Collectively, these data allow us to put together a comprehensive model of the structure and function of the polyadenylation ternary complex consisting of VP39, VP55 and RNA. ^

Post-transcriptional Regulation of Mammalian Gene Expression in Non-coding Region of Target RNA

Tumor Suppressor Candidate 2 (TUSC2) is a novel tumor suppressor gene located in the human chromosome 3p21.3 region. TUSC2 mRNA transcripts could be detected on Northern blots in both normal lung and some lung cancer cell lines, but no endogenous TUSC2 protein could be detected in a majority of lung cancer cell lines. Mechanisms regulating TUSC2 protein expression and its inactivation in primary lung cancer cells are largely unknown. We investigated the role of the 5’- and 3’-untranslated regions (UTRs) of the TUSC2 gene in the regulation of TUSC2 protein expression. We found that two small upstream open-reading frames (uORFs) in the 5’UTR of TUSC2 could markedly inhibit the translational initiation of TUSC2 protein by interfering with the “scanning” of the ribosome initiation complexes. Site-specific stem-loop array reverse transcription-polymerase chain reaction (SLA-RT-PCR) verified several micoRNAs (miRNAs) targeted at 3’UTR and directed TUSC2 cleavage and degradation. In addition, we used the established let-7-targeted high mobility group A2 (Hmga2) mRNA as a model system to study the mechanism of regulation of target mRNA by miRNAs in mammalian cells under physiological conditions. There have been no evidence of direct link between mRNA downregulation and mRNA cleavages mediated by miRNAs. Here we showed that the endonucleolytic cleavages on mRNAs were initiated by mammalian miRNA in seed pairing style. Let-7 directed cleavage activities among the eight predicted potential target sites have varied efficiency, which are influenced by the positional and the structural contexts in the UTR. The 5’ cleaved RNA fragments were mostly oligouridylated at their 3’-termini and accumulated for delayed 5’–3’ degradation. RNA fragment oligouridylation played important roles in marking RNA fragments for delayed bulk degradation and in converting RNA degradation mode from 3’–5’ to 5’–3’ with cooperative efforts from both endonucleolytic and non-catalytic miRNA-induced silencing complex (miRISC). Our findings point to a mammalian miRNA-mediated mechanism for the regulation of mRNA that miRNA can decrease target mRNA through target mRNA cleavage and uridine addition

A human homolog of the Saccharomyces cerevisiae REV3 gene, which encodes the catalytic subunit of DNA polymerase ζ

To get a better understanding of mutagenic mechanisms in humans, we have cloned and sequenced the human homolog of the Saccharomyces cerevisiae REV3 gene. The yeast gene encodes the catalytic subunit of DNA polymerase ζ, a nonessential enzyme that is thought to carry out translesion replication and is responsible for virtually all DNA damage-induced mutagenesis and the majority of spontaneous mutagenesis. The human gene encodes an expected protein of 3,130 residues, about twice the size of the yeast protein (1,504 aa). The two proteins are 29% identical in an amino-terminal region of ≈340 residues, 39% identical in a carboxyl-terminal region of ≈850 residues, and 29% identical in a 55-residue region in the middle of the two genes. The sequence of the expected protein strongly predicts that it is the catalytic subunit of a DNA polymerase of the pol ζ type; the carboxyl-terminal domain possesses, in the right order, the six motifs characteristic of eukaryotic DNA polymerases, most closely resembles yeast pol ζ among all polymerases in the GenBank database, and is different from the human α, δ, and ɛ enzymes. Human cells expressing high levels of an hsREV3 antisense RNA fragment grow normally, but show little or no UV-induced mutagenesis and are slightly more sensitive to killing by UV. The human gene therefore appears to carry out a function similar to that of its yeast counterpart.

Polyadenylation of stable RNA precursors in vivo

Polyadenylation at the 3′ terminus has long been considered a specific feature of mRNA and a few other unstable RNA species. Here we show that stable RNAs in Escherichia coli can be polyadenylated as well. RNA molecules with poly(A) tails are the major products that accumulate for essentially all stable RNA precursors when RNA maturation is slowed because of the absence of processing exoribonucleases; poly(A) tails vary from one to seven residues in length. The polyadenylation process depends on the presence of poly(A) polymerase I. A stochastic competition between the exoribonucleases and poly(A) polymerase is proposed to explain the accumulation of polyadenylated RNAs. These data indicate that polyadenylation is not unique to mRNA, and its widespread occurrence suggests that it serves a more general function in RNA metabolism.

Numbers and Organization of RNA Polymerases, Nascent Transcripts, and Transcription Units in HeLa Nuclei

Using HeLa cells, we have developed methods to determine 1) the number of RNA polymerases that are active at any moment, 2) the number of transcription sites, and 3) the number of polymerases associated with one transcription unit. To count engaged polymerases, cells were encapsulated in agarose, permeabilized, treated with ribonuclease, and the now-truncated transcripts extended in [32P]uridine triphosphate; then, the number of growing transcripts was calculated from the total number of nucleotides incorporated and the average increment in length of the transcripts. Approximately 15,000 transcripts were elongated by polymerase I, and ∼75,000 were elongated by polymerases II and III. Transcription sites were detected after the cells were grown in bromouridine for <2.5 min, after which the resulting bromo-RNA was labeled with gold particles; electron microscopy showed that most extranucleolar transcripts were concentrated in ∼2400 sites with diameters of ∼80 nm. The number of polymerases associated with a transcription unit was counted after templates were spread over a large area; most extranucleolar units were associated with one elongating complex. These results suggest that many templates are attached in a “cloud” of loops around a site; each site, or transcription “factory,” would contain ∼30 active polymerases and associated transcripts.

Host factor Hfq of Escherichia coli stimulates elongation of poly(A) tails by poly(A) polymerase I

Current evidence suggests that the length of poly(A) tails of bacterial mRNAs result from a competition between poly(A) polymerase and exoribonucleases that attack the 3′ ends of RNAs. Here, we show that host factor Hfq is also involved in poly(A) tail metabolism. Inactivation of the hfq gene reduces the length of poly(A) tails synthesized at the 3′ end of the rpsO mRNA by poly(A) polymerase I in vivo. In vitro, Hfq stimulates synthesis of long tails by poly(A) polymerase I. The strong binding of Hfq to oligoadenylated RNA probably explains why it stimulates elongation of primers that already harbor tails of 20–35 A. Polyadenylation becomes processive in the presence of Hfq. The similar properties of Hfq and the PABPII poly(A) binding protein, which stimulates poly(A) tail elongation in mammals, indicates that similar mechanisms control poly(A) tail synthesis in prokaryotes and eukaryotes.

Comparison of the 5′ nuclease activities of Taq DNA polymerase and its isolated nuclease domain

Many eubacterial DNA polymerases are bifunctional molecules having both polymerization (P) and 5′ nuclease (N) activities, which are contained in separable domains. We previously showed that the DNA polymerase I of Thermus aquaticus (TaqNP) endonucleolytically cleaves DNA substrates, releasing unpaired 5′ arms of bifurcated duplexes. Here, we compare the substrate specificities of TaqNP and the isolated 5′ nuclease domain of this enzyme, TaqN. Both enzymes are significantly activated by primer oligonucleotides that are hybridized to the 3′ arm of the bifurcation; optimal stimulation requires overlap of the 3′ terminal nucleotide of the primer with the terminal base pair of the duplex, but the terminal nucleotide need not hybridize to the complementary strand in the substrate. In the presence of Mn2+ ions, TaqN can cleave both RNA and circular DNA at structural bifurcations. Certain anti-TaqNP mAbs block cleavage by one or both enzymes, whereas others can stimulate cleavage of nonoptimal substrates.

Evolution of a quadripartite hybrid virus by interspecific exchange and recombination between replicase components of two related tripartite RNA viruses

Cucumber mosaic virus (CMV) and tomato aspermy virus (TAV) belong to the Cucumovirus genus. They have a tripartite genome consisting of single-stranded RNAs, designated 1, 2, and 3. Previous studies have shown that viable pseudorecombinants could be created in vitro by reciprocal exchanges between CMV and TAV RNA 3, but exchanges of RNAs 1 and 2 were replication deficient. When we coinoculated CMV RNAs 2 and 3 along with TAV RNAs 1 and 2 onto Nicotiana benthamiana, a hybrid quadripartite virus appeared that consisted of TAV RNA 1, CMV RNAs 2 and 3, and a distinctive chimeric RNA originating from a recombination between CMV RNA 2 and the 3′-terminal 320 nucleotides of TAV RNA 2. This hybrid arose by means of segment reassortment and RNA recombination to produce an interspecific hybrid with the TAV helicase subunit and the CMV polymerase subunit. To our knowledge, this is the first report demonstrating the evolution of a new plant or animal virus strain containing an interspecific hybrid replicase complex.

Kinetics of formation of hypoxanthine containing base pairs by HIV-RT: RNA template effects on the base substitution frequencies

Hypoxanthine (H), the deamination product of adenine, has been implicated in the high frequency of A to G transitions observed in retroviral and other RNA genomes. Although H·C base pairs are thermodynamically more stable than other H·N pairs, polymerase selection may be determined in part by kinetic factors. Therefore, the hypoxanthine induced substitution pattern resulting from replication by viral polymerases may be more complex than that predicted from thermodynamics. We have examined the steady-state kinetics of formation of base pairs opposite template H in RNA by HIV-RT, and for the incorporation of dITP during first- and second-strand synthesis. Hypoxanthine in an RNA template enhances the k2app for pairing with standard dNTPs by factors of 10–1000 relative to adenine at the same sequence position. The order of base pairing preferences for H in RNA was observed to be H·C >> H·T > H·A > H·G. Steady-state kinetics of insertion for all possible mispairs formed with dITP were examined on RNA and DNA templates of identical sequence. Insertion of dITP opposite all bases occurs 2–20 times more frequently on RNA templates. This bias for higher insertion frequencies on RNA relative to DNA templates is also observed for formation of mispairs at template A. This kinetic advantage afforded by RNA templates for mismatches and pairing involving H suggests a higher induction of mutations at adenines during first-strand synthesis by HIV-RT.