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.

Multiple modes of cellular activation and virus transmission in HIV infection: A role for chronically and latently infected cells in sustaining viral replication

CD4+ T cell activation, required for virus replication in these cells, occurs in local microenvironmental domains in transient bursts. Thus, although most HIV originates from short-lived virus-producing cells, it is unlikely that chronic infection is generally sustained in rapid continuous cycles of productive infection as has been proposed. Such continuity of productive infection cycles would depend on efficient long-range transmission of HIV from one set of domains to another, in turn requiring the maintenance of sufficiently high concentrations of cell-free virus across lymphoid tissues at all times. By contrast, long-lived cellular sources of HIV maintain the capacity to infect newly activated cells at close range despite the temporal and spatial discontinuities of activation events. Such proximal activation and transmission (PAT) involving chronically and latently infected cells may be responsible for sustained infection, particularly when viral loads are low. Once CD4 cells are productively infected through PAT, they can infect other activated cells in their immediate vicinity. Such events propagate locally but generally do not spread systemically, unlike in the acute phase of the infection, because of the early establishment of protective anergy. Importantly, antiretroviral drug treatment is likely to differentially impact long-range transmission and PAT.

Efficient translesion replication in the absence of Escherichia coli Umu proteins and 3′–5′ exonuclease proofreading function

Translesion replication (TR) past a cyclobutane pyrimidine dimer in Escherichia coli normally requires the UmuD′2C complex, RecA protein, and DNA polymerase III holoenzyme (pol III). However, we find that efficient TR can occur in the absence of the Umu proteins if the 3′–5′ exonuclease proofreading activity of the pol III ɛ-subunit also is disabled. TR was measured in isogenic uvrA6 ΔumuDC strains carrying the dominant negative dnaQ allele, mutD5, or ΔdnaQ spq-2 mutations by transfecting them with single-stranded M13-based vectors containing a specifically located cis-syn T–T dimer. As expected, little TR was observed in the ΔumuDC dnaQ+ strain. Surprisingly, 26% TR occurred in UV-irradiated ΔumuDC mutD5 cells, one-half the frequency found in a uvrA6 umuDC+mutD5 strain. lexA3 (Ind−) derivatives of the strains showed that this TR was contingent on two inducible functions, one LexA-dependent, responsible for ≈70% of the TR, and another LexA-independent, responsible for the remaining ≈30%. Curiously, the ΔumuDC ΔdnaQ spq-2 strain exhibited only the LexA-independent level of TR. The cause of this result appears to be the spq-2 allele, a dnaE mutation required for viability in ΔdnaQ strains, since introduction of spq-2 into the ΔumuDC mutD5 strain also reduces the frequency of TR to the LexA-independent level. The molecular mechanism responsible for the LexA-independent TR is unknown but may be related to the UVM phenomenon [Palejwala, V. A., Wang, G. E., Murphy, H. S. & Humayun, M. Z. (1995) J. Bacteriol. 177, 6041–6048]. LexA-dependent TR does not result from the induction of pol II, since TR in the ΔumuDC mutD5 strain is unchanged by introduction of a ΔpolB mutation.

Yeast mutations in multiple complementation groups inhibit brome mosaic virus RNA replication and transcription and perturb regulated expression of the viral polymerase-like gene

Brome mosaic virus (BMV), a member of the alphavirus-like superfamily of positive-strand RNA viruses, encodes two proteins, 1a and 2a, that interact with each other, with unidentified host proteins, and with host membranes to form the viral RNA replication complex. Yeast expressing 1a and 2a support replication and subgenomic mRNA synthesis by BMV RNA3 derivatives. Using a multistep selection and screening process, we have isolated yeast mutants in multiple complementation groups that inhibit BMV-directed gene expression. Three complementation groups, represented by mutants mab1–1, mab2–1, and mab3–1 (for maintenance of BMV functions), were selected for initial study. Each of these mutants has a single, recessive, chromosomal mutation that inhibits accumulation of positive- and negative-strand RNA3 and subgenomic mRNA. BMV-directed gene expression was inhibited when the RNA replication template was introduced by in vivo transcription from DNA or by transfection of yeast with in vitro transcripts, confirming that cytoplasmic RNA replication steps were defective. mab1–1, mab2–1, and mab3–1 slowed yeast growth to varying degrees and were temperature-sensitive, showing that the affected genes contribute to normal cell growth. In wild-type yeast, expression of the helicase-like 1a protein increased the accumulation of 2a mRNA and the polymerase-like 2a protein, revealing a new level of viral regulation. In association with their other effects, mab1–1 and mab2–1 blocked the ability of 1a to stimulate 2a mRNA and protein accumulation, whereas mab3–1 had elevated 2a protein accumulation. Together, these results show that BMV RNA replication in yeast depends on multiple host genes, some of which directly or indirectly affect the regulated expression and accumulation of 2a.

Respiratory syncytial virus (RSV) SH and G proteins are not essential for viral replication in vitro: Clinical evaluation and molecular characterization of a cold-passaged, attenuated RSV subgroup B mutant

A live, cold-passaged (cp) candidate vaccine virus, designated respiratory syncytial virus (RSV) B1 cp-52/2B5 (cp-52), replicated efficiently in Vero cells, but was found to be overattenuated for RSV-seronegative infants and children. Sequence analysis of reverse-transcription–PCR-amplified fragments of this mutant revealed a large deletion spanning most of the coding sequences for the small hydrophobic (SH) and attachment (G) proteins. Northern blot analysis of cp-52 detected multiple unique read-through mRNAs containing SH and G sequences, consistent with a deletion mutation spanning the SH:G gene junction. Immunological studies confirmed that an intact G glycoprotein was not produced by the cp-52 virus. Nonetheless, cp-52 was infectious and replicated to high titer in tissue culture despite the absence of the viral surface SH and G glycoproteins. Thus, our characterization of this negative-strand RNA virus identified a novel replication-competent deletion mutant lacking two of its three surface glycoproteins. The requirement of SH and G for efficient replication in vivo suggests that selective deletion of one or both of these RSV genes may provide an alternative or additive strategy for developing an optimally attenuated vaccine candidate.

Two forms of replication initiator protein: Positive and negative controls

The pir gene of plasmid R6K encodes the protein, π, a replication and transcription factor. Two translational options for the pir gene give rise to two forms of π protein: a 35.0-kDa form (π35.0) and a shortened 30.5-kDa form (π30.5). Although both proteins bind to a series of 22-bp direct repeats essential for plasmid R6K replication, only π35.0 can bind to a site in the (A⋅T)-rich segment of its γ ori and activate the γ ori in vivo and in vitro. However, unlike π35.0, π30.5can inhibit in vivo and in vitro replication (activated by π35.0). We propose that the two forms of π might have distinct functions in replication. We show that although both forms of π produce dimers, the nature of these dimers is not identical. The N-terminal 37 amino acid residues appear to control the formation of the more stable π35.0 dimers, whereas another, apparently weaker interface holds together dimers of π30.5. We speculate that the leucine zipper-like motif, absent in π30.5, controls very specific functions of π protein.

Evaluation of the role of heterogeneous nuclear ribonucleoprotein A1 as a host factor in murine coronavirus discontinuous transcription and genome replication

Viruses with RNA genomes often capture and redirect host cell components to assist in mechanisms particular to RNA-dependent RNA synthesis. The nidoviruses are an order of positive-stranded RNA viruses, comprising coronaviruses and arteriviruses, that employ a unique strategy of discontinuous transcription, producing a series of subgenomic mRNAs linking a 5′ leader to distal portions of the genome. For the prototype coronavirus mouse hepatitis virus (MHV), heterogeneous nuclear ribonucleoprotein (hnRNP) A1 has been shown to be able to bind in vitro to the negative strand of the intergenic sequence, a cis-acting element found in the leader RNA and preceding each downstream ORF in the genome. hnRNP A1 thus has been proposed as a host factor in MHV transcription. To test this hypothesis genetically, we initially constructed MHV mutants with a very high-affinity hnRNP A1 binding site inserted in place of, or adjacent to, an intergenic sequence in the MHV genome. This inserted hnRNP A1 binding site was not able to functionally replace, or enhance transcription from, the intergenic sequence. This finding led us to test more directly the role of hnRNP A1 by analysis of MHV replication and RNA synthesis in a murine cell line that does not express this protein. The cellular absence of hnRNP A1 had no detectable effect on the production of infectious virus, the synthesis of genomic RNA, or the quantity or quality of subgenomic mRNAs. These results strongly suggest that hnRNP A1 is not a required host factor for MHV discontinuous transcription or genome replication.

An initial ATP-independent step in the unwinding of a herpes simplex virus type I origin of replication by a complex of the viral origin-binding protein and single-strand DNA-binding protein

Using a spectrophotometric assay that measures the hyperchromicity that accompanies the unwinding of a DNA duplex, we have identified an ATP-independent step in the unwinding of a herpes simplex virus type 1 (HSV-1) origin of replication, Oris, by a complex of the HSV-1 origin binding protein (UL9 protein) and the HSV-1 single-strand DNA binding protein (ICP8). The sequence unwound is the 18-bp A + T-rich segment that links the two high-affinity UL9 protein binding sites, boxes I and II of Oris. P1 nuclease sensitivity of Oris and single-strand DNA-dependent ATPase measurements of the UL9 protein indicate that, at 37°C, the A + T-rich segment is sufficiently single stranded to permit the binding of ICP8. Binding of the UL9 protein to boxes I and II then results in the formation of the UL9 protein–ICP8 complex, that can, in the absence of ATP, promote unwinding of the A + T-rich segment. On addition of ATP, the helicase activity of the UL9 protein–ICP8 complex can unwind boxes I and II, permitting access of the replication machinery to the Oris sequences.

Mitochondrial Single-stranded DNA-binding Protein Is Required for Mitochondrial DNA Replication and Development in Drosophila melanogaster

The discovery that several inherited human diseases are caused by mtDNA depletion has led to an increased interest in the replication and maintenance of mtDNA. We have isolated a new mutant in the lopo (low power) gene from Drosophila melanogaster affecting the mitochondrial single-stranded DNA-binding protein (mtSSB), which is one of the key components in mtDNA replication and maintenance. lopo1 mutants die late in the third instar before completion of metamorphosis because of a failure in cell proliferation. Molecular, histochemical, and physiological experiments show a drastic decrease in mtDNA content that is coupled with the loss of respiration in these mutants. However, the number and morphology of mitochondria are not greatly affected. Immunocytochemical analysis shows that mtSSB is expressed in all tissues but is highly enriched in proliferating tissues and in the developing oocyte. lopo1 is the first mtSSB mutant in higher eukaryotes, and its analysis demonstrates the essential function of this gene in development, providing an excellent model to study mitochondrial biogenesis in animals.

Similarities between the DNA replication initiators of Gram-negative bacteria plasmids (RepA) and eukaryotes (Orc4p)/archaea (Cdc6p)

The proteins responsible for the initiation of DNA replication are thought to be essentially unrelated in bacteria and archaea/eukaryotes. Here we show that RepA, the initiator from the Pseudomonas plasmid pPS10, and the C-terminal domain of ScOrc4p, a subunit of Saccharomyces cerevisiae (Sc) origin recognition complex (ORC), share sequence similarities. Based on biochemical and spectroscopic evidence, these similarities include common structural elements, such as a winged-helix domain and a leucine-zipper dimerization motif. We have also found that ScOrc4p, as previously described for RepA-type initiators, interacts with chaperones of the Hsp70 family both in vitro and in vivo, most probably to regulate the assembly of active ORC. In evolutionary terms, our results are compatible with the recruitment of the same protein module for initiation of DNA replication by the ancestors of present-day Gram-negative bacteria plasmids, archaea, and eukaryotes.

DNA polymerase ɛ is required for coordinated and efficient chromosomal DNA replication in Xenopus egg extracts

DNA polymerase ɛ (Polɛ) is thought to be involved in DNA replication, repair, and cell-cycle checkpoint control in eukaryotic cells. Although the requirement of other replicative DNA polymerases, DNA polymerases α and δ (Polα and δ), for chromosomal DNA replication has been well documented by genetic and biochemical studies, the precise role, if any, of Polɛ in chromosomal DNA replication is still obscure. Here we show, with the use of a cell-free replication system with Xenopus egg extracts, that Xenopus Polɛ is indeed required for chromosomal DNA replication. In Polɛ-depleted extracts, the elongation step of chromosomal DNA replication is markedly impaired, resulting in significant reduction of the overall DNA synthesis as well as accumulation of small replication intermediates. Moreover, despite the decreased DNA synthesis, excess amounts of Polα are loaded onto the chromatin template in Polɛ-depleted extracts, indicative of the failure of proper assembly of DNA synthesis machinery at the fork. These findings strongly suggest that Polɛ, along with Polα and Polδ, is necessary for coordinated chromosomal DNA replication in eukaryotic cells.

UV-induced Hyperphosphorylation of Replication Protein A Depends on DNA Replication and Expression of ATM Protein

Exposure to DNA-damaging agents triggers signal transduction pathways that are thought to play a role in maintenance of genomic stability. A key protein in the cellular processes of nucleotide excision repair, DNA recombination, and DNA double-strand break repair is the single-stranded DNA binding protein, RPA. We showed previously that the p34 subunit of RPA becomes hyperphosphorylated as a delayed response (4–8 h) to UV radiation (10–30 J/m2). Here we show that UV-induced RPA-p34 hyperphosphorylation depends on expression of ATM, the product of the gene mutated in the human genetic disorder ataxia telangiectasia (A-T). UV-induced RPA-p34 hyperphosphorylation was not observed in A-T cells, but this response was restored by ATM expression. Furthermore, purified ATM kinase phosphorylates the p34 subunit of RPA complex in vitro at many of the same sites that are phosphorylated in vivo after UV radiation. Induction of this DNA damage response was also dependent on DNA replication; inhibition of DNA replication by aphidicolin prevented induction of RPA-p34 hyperphosphorylation by UV radiation. We postulate that this pathway is triggered by the accumulation of aberrant DNA replication intermediates, resulting from DNA replication fork blockage by UV photoproducts. Further, we suggest that RPA-p34 is hyperphosphorylated as a participant in the recombinational postreplication repair of these replication products. Successful resolution of these replication intermediates reduces the accumulation of chromosomal aberrations that would otherwise occur as a consequence of UV radiation.

Targeting peptide nucleic acid (PNA) oligomers to mitochondria within cells by conjugation to lipophilic cations: implications for mitochondrial DNA replication, expression and disease

The selective manipulation of mitochondrial DNA (mtDNA) replication and expression within mammalian cells has proven difficult. One promising approach is to use peptide nucleic acid (PNA) oligomers, nucleic acid analogues that bind selectively to complementary DNA or RNA sequences inhibiting replication and translation. However, the potential of PNAs is restricted by the difficulties of delivering them to mitochondria within cells. To overcome this problem we conjugated a PNA 11mer to a lipophilic phosphonium cation. Such cations are taken up by mitochondria through the lipid bilayer driven by the membrane potential across the inner membrane. As anticipated, phosphonium–PNA (ph–PNA) conjugates of 3.4–4 kDa were imported into both isolated mitochondria and mitochondria within human cells in culture. This was confirmed by using an ion-selective electrode to measure uptake of the ph–PNA conjugates; by cell fractionation in conjunction with immunoblotting; by confocal microscopy; by immunogold-electron microscopy; and by crosslinking ph–PNA conjugates to mitochondrial matrix proteins. In all cases dissipating the mitochondrial membrane potential with an uncoupler prevented ph–PNA uptake. The ph–PNA conjugate selectively inhibited the in vitro replication of DNA containing the A8344G point mutation that causes the human mtDNA disease ‘myoclonic epilepsy and ragged red fibres’ (MERRF) but not the wild-type sequence that differs at a single nucleotide position. Therefore these modified PNA oligomers retain their selective binding to DNA and the lipophilic cation delivers them to mitochondria within cells. When MERRF cells were incubated with the ph–PNA conjugate the ratio of MERRF to wild-type mtDNA was unaffected, even though the ph–PNA content of the mitochondria was sufficient to inhibit MERRF mtDNA replication in a cell-free system. This unexpected finding suggests that nucleic acid derivatives cannot bind their complementary sequences during mtDNA replication. In summary, we have developed a new strategy for targeting PNA oligomers to mitochondria and used it to determine the effects of PNA on mutated mtDNA replication in cells. This work presents new approaches for the manipulation of mtDNA replication and expression, and will assist in the development of therapies for mtDNA diseases.

Rescue of stalled replication forks by RecG: Simultaneous translocation on the leading and lagging strand templates supports an active DNA unwinding model of fork reversal and Holliday junction formation

Modification of damaged replication forks is emerging as a crucial factor for efficient chromosomal duplication and the avoidance of genetic instability. The RecG helicase of Escherichia coli, which is involved in recombination and DNA repair, has been postulated to act on stalled replication forks to promote replication restart via the formation of a four-stranded (Holliday) junction. Here we show that RecG can actively unwind the leading and lagging strand arms of model replication fork structures in vitro. Unwinding is achieved in each case by simultaneous interaction with and translocation along both the leading and lagging strand templates at a fork. Disruption of either of these interactions dramatically inhibits unwinding of the opposing duplex arm. Thus, RecG translocates simultaneously along two DNA strands, one with 5′-3′ and the other with 3′-5′ polarity. The unwinding of both nascent strands at a damaged fork, and their subsequent annealing to form a Holliday junction, may explain the ability of RecG to promote replication restart. Moreover, the preferential binding of partial forks lacking a leading strand suggests that RecG may have the ability to target stalled replication intermediates in vivo in which lagging strand synthesis has continued beyond the leading strand.