Mapping vaccinia virus DNA replication origins at nucleotide level by deep sequencing

Parapoxvirus species revisited by whole genome sequencing: A retrospective analysis of bovine virus isolates

Parapoxviruses (PPV) of animals are spread worldwide. While the Orf virus (ORFV) species is a molecularly well-characterized prototype pathogen of small ruminants, the genomes of virus species affecting large ruminants, namely Bovine papular stomatitis virus (BPSV) and Pseudocowpox virus (PCPV), are less well known. Using Nanopore sequencing we retrospectively show the whole genome sequences (WGS) of six BPSV, three PCPV isolates and an attenuated ORFV strain, originating from different geographic locations. A phylogenetic tree shows that the de novo assembled genomes belong to PPV species including WGS of reference PPV. Remarkably, Nanopore sequencing allowed the molecular resolution of inverted terminal repeats (ITR) and the hairpin loop within the de novo assembled WGS. Additionally, peculiarities regarding map location of two genes and the heterogeneity of a genomic region were noted. Details for the molecular variability of an interferon response modulatory gene (ORF116) and the PCPV specificity of gene 073.5 are reported. In summary, WGS gained by Nanopore sequencing allowed analysis of complete PPV genomes and confident virus species attribution within a phylogenetic tree avoiding uncertainty of limited gene-based diagnostics. Nanopore-based WGS provides robust comparison of PPV genomes and reliable identity determination of new Poxviruses.

Structure and flexibility of the DNA polymerase holoenzyme of vaccinia virus

The year 2022 was marked by the mpox outbreak caused by the human monkeypox virus (MPXV), which is approximately 98% identical to the vaccinia virus (VACV) at the sequence level with regard to the proteins involved in DNA replication. We present the production in the baculovirus-insect cell system of the VACV DNA polymerase holoenzyme, which consists of the E9 polymerase in combination with its co-factor, the A20-D4 heterodimer. This led to the 3.8 Å cryo-electron microscopy (cryo-EM) structure of the DNA-free form of the holoenzyme. The model of the holoenzyme was constructed from high-resolution structures of the components of the complex and the A20 structure predicted by AlphaFold 2. The structures do not change in the context of the holoenzyme compared to the previously determined crystal and NMR structures, but the E9 thumb domain became disordered. The E9-A20-D4 structure shows the same compact arrangement with D4 folded back on E9 as observed for the recently solved MPXV holoenzyme structures in the presence and the absence of bound DNA. A conserved interface between E9 and D4 is formed by a cluster of hydrophobic residues. Small-angle X-ray scattering data show that other, more open conformations of E9-A20-D4 without the E9-D4 contact exist in solution using the flexibility of two hinge regions in A20. Biolayer interferometry (BLI) showed that the E9-D4 interaction is indeed weak and transient in the absence of DNA although it is very important, as it has not been possible to obtain viable viruses carrying mutations of key residues within the E9-D4 interface.

Biology of Variola Virus

Variola virus is an anthroponotic agent that belongs to the orthopoxvirus family. It is an etiological agent of smallpox, an ancient disease that caused massive mortality of human populations. Twentieth century has witnessed the death of about 300 million people due to the unavailability of an effective vaccine. Early detection is the primary strategy to prevent an outbreak of smallpox. Variola virus forms the characteristic pus-filled pustules and centrifugal rash distribution in the infected patients while transmission occurs mainly through respiratory droplets during the early stage of infection. No antiviral drugs are approved for variola virus till date. Generation of first-generation vaccines helped in the eradication of smallpox which was declared by the World Health Organization.

Structural insights into the assembly and mechanism of mpox virus DNA polymerase complex F8-A22-E4-H5

The DNA replication of mpox virus is performed by the viral polymerase F8 and also requires other viral factors, including processivity factor A22, uracil DNA glycosylase E4, and phosphoprotein H5. However, the molecular roles of these viral factors remain unclear. Here, we characterize the structures of F8-A22-E4 and F8-A22-E4-H5 complexes in the presence of different primer-template DNA substrates. E4 is located upstream of F8 on the template single-stranded DNA (ssDNA) and is catalytically active, highlighting a functional coupling between DNA base-excision repair and DNA synthesis. Moreover, H5, in the form of tetramer, binds to the double-stranded DNA (dsDNA) region downstream of F8 in a similar position as PCNA (proliferating cell nuclear antigen) does in eukaryotic polymerase complexes. Omission of H5 or disruption of its DNA interaction showed a reduced synthesis of full-length DNA products. These structures provide snapshots for the working cycle of the polymerase and generate insights into the mechanisms of these essential factors in viral DNA replication.

An APOBEC3 Mutational Signature in the Genomes of Human-Infecting Orthopoxviruses

The ongoing worldwide monkeypox outbreak is caused by viral lineages (globally referred to as hMPXV1) that are related to but distinct from clade IIb MPXV viruses transmitted within Nigeria. Analysis of the genetic differences has indicated that APOBEC-mediated editing might be responsible for the unexpectedly high number of mutations observed in hMPXV1 genomes. Here, using 1,624 publicly available hMPXV1 sequences, we analyzed the mutations that accrued between 2017 and the emergence of the current predominant variant (B.1), as well as those that that have been accumulating during the 2022 outbreak. We confirmed an overwhelming prevalence of C-to-T and G-to-A mutations, with a sequence context (5'-TC-3') consistent with the preferences of several human APOBEC3 enzymes. We also found that mutations preferentially occur in highly expressed viral genes, although no transcriptional asymmetry was observed. A comparison of the mutation spectrum and context was also performed against the human-specific variola virus (VARV) and the zoonotic cowpox virus (CPXV), as well as fowlpox virus (FWPV). The results indicated that in VARV genomes, C-to-T and G-to-A changes were more common than the opposite substitutions, although the effect was less marked than for hMPXV1. Conversely, no preference toward C-to-T and G-to-A changes was observed in CPXV and FWPV. Consistently, the sequence context of C-to-T changes confirmed a preference for a T in the -1 position for VARV, but not for CPXV or FWPV. Overall, our results strongly support the view that, irrespective of the transmission route, orthopoxviruses infecting humans are edited by the host APOBEC3 enzymes. IMPORTANCE Analysis of the viral lineages responsible for the 2022 monkeypox outbreak suggested that APOBEC enzymes are driving hMPXV1 evolution. Using 1,624 public sequences, we analyzed the mutations that accumulated between 2017 and the emergence of the predominant variant and those that characterize the last outbreak. We found that the mutation spectrum of hMPXV1 has been dominated by TC-to-TT and GA-to-AA changes, consistent with the editing activity of human APOBEC3 proteins. We also found that mutations preferentially affect highly expressed viral genes, possibly because transcription exposes single-stranded DNA (ssDNA), a target of APOBEC3 editing. Notably, analysis of the human-specific variola virus (VARV) and the zoonotic cowpox virus (CPXV) indicated that in VARV genomes, TC-to-TT and GA-to-AA changes are likewise extremely frequent. Conversely, no preference toward TC-to-TT and GA-to-AA changes is observed in CPXV. These results suggest that APOBEC3 proteins have an impact on the evolution of different human-infecting orthopoxviruses.

Evolution of giant pandoravirus revealed by CRISPR/Cas9

Giant viruses (GVs) are a hotspot of unresolved controversies since their discovery, including the definition of "Virus" and their origin. While increasing knowledge of genome diversity has accumulated, GV functional genomics was largely neglected. Here, we describe an experimental framework to genetically modify nuclear GVs and their host Acanthamoeba castellanii using CRISPR/Cas9, shedding light on the evolution from small icosahedral viruses to amphora-shaped GVs. Ablation of the icosahedral major capsid protein in the phylogenetically-related mollivirus highlights a transition in virion shape and size. We additionally demonstrate the existence of a reduced core essential genome in pandoravirus, reminiscent of their proposed smaller ancestors. This proposed genetic expansion led to increased genome robustness, indicating selective pressures for adaptation to uncertain environments. Overall, we introduce new tools for manipulation of the unexplored genome of nuclear GVs and provide experimental evidence suggesting that viral gigantism has aroused as an emerging trait.

The Life Cycle of the Vaccinia Virus Genome

Poxviruses, of which vaccinia virus is the prototype, are a large family of double-stranded DNA viruses that replicate exclusively in the cytoplasm of infected cells. This physical and genetic autonomy from the host cell nucleus necessitates that these viruses encode most, if not all, of the proteins required for replication in the cytoplasm. In this review, we follow the life of the viral genome through space and time to address some of the unique challenges that arise from replicating a 195-kb DNA genome in the cytoplasm. We focus on how the genome is released from the incoming virion and deposited into the cytoplasm; how the endoplasmic reticulum is reorganized to form a replication factory, thereby compartmentalizing and helping to protect the replicating genome from immune sensors; how the cellular milieu is tailored to support high-fidelity replication of the genome; and finally, how newly synthesized genomes are faithfully and specifically encapsidated into new virions. Expected final online publication date for the Annual Review of Virology, Volume 9 is September 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

The mismatched nucleotides encoded in vaccinia virus flip-and-flop hairpin telomeres serve an essential role in virion maturation

Poxvirus genomes consist of a linear duplex DNA that ends in short inverted and complementary hairpin structures. These elements also encode loops and mismatches that likely serve a role in genome packaging and perhaps replication. We constructed mutant vaccinia viruses (VACV) where the native hairpins were replaced by altered forms and tested effects on replication, assembly, and virulence. Our studies showed that structure, not sequence, likely determines function as one can replace an Orthopoxvirus (VACV) hairpin with one copied from a Leporipoxvirus with no effect on growth. Some loops can be deleted from VACV hairpins with little effect, but VACV bearing too few mismatches grew poorly and we couldn’t recover viruses lacking all mismatches. Further studies were conducted using a mutant bearing only one of six mismatches found in wild-type hairpins (SΔ1Δ3–6). This virus grew to ~20-fold lower titers, but neither DNA synthesis nor telomere resolution was affected. However, the mutant exhibited a particle-to-PFU ratio 10-20-fold higher than wild-type viruses and p4b/4b core protein processing was compromised, indicating an assembly defect. Electron microscopy showed that SΔ1Δ3–6 mutant development was blocked at the immature virus (IV) stage, which phenocopies known effects of I1L mutants. Competitive DNA binding assays showed that recombinant I1 protein had less affinity for the SΔ1Δ3–6 hairpin than the wild-type hairpin. The SΔ1Δ3–6 mutant was also attenuated when administered to SCID-NCR mice by tail scarification. Mice inoculated with viruses bearing wild-type hairpins exhibited a median survival of 30–37 days, while mice infected with SΔ1Δ3–6 virus survived >70 days. Persistent infections favor genetic reversion and genome sequencing detected one example where a small duplication near the hairpin tip likely created a new loop. These observations show that mismatches serve a critical role in genome packaging and provide new insights into how VACV “flip and flop” telomeres are arranged.

Poxviruses employ linear double-stranded DNA genomes that end in incompletely base-paired hairpin termini. These mismatched ends are thought to serve some role in virus assembly, and perhaps replication, but have not been amenable to genetic analysis. In this study we used a synthetic virology approach to alter the sequence and structure of these elements. Our research shows that although the encoded structures are of critical importance for function, the sequences are not because one can swap the ends of viruses from different poxviruses without affecting growth. When one tries to progressively delete the mismatches that are found at these ends (the telomeres) of wild-type genomes, it creates an assembly defect which shows up as an increase in the number of virus particles per infectious unit and an accumulation of incompletely assembled viruses. Electron microscopy showed that the development of mutant viruses is blocked at a stage after DNA is packaged but before the particles fully mature. This investigation supports earlier studies that had identified the telomeres as being sites where virus proteins bind and promote packaging. Viruses bearing these mutant telomeres are also less virulent but can still serve as vaccines to protect mice from a lethal virus challenge.

Nanopore sequencing and de novo assembly of a misidentified Camelpox vaccine reveals putative epigenetic modifications and alternate protein signal peptides

DNA viruses can exploit host cellular epigenetic processes to their advantage; however, the epigenome status of most DNA viruses remains undetermined. Third generation sequencing technologies allow for the identification of modified nucleotides from sequencing experiments without specialized sample preparation, permitting the detection of non-canonical epigenetic modifications that may distinguish viral nucleic acid from that of their host, thus identifying attractive targets for advanced therapeutics and diagnostics. We present a novel nanopore de novo assembly pipeline used to assemble a misidentified Camelpox vaccine. Two confirmed deletions of this vaccine strain in comparison to the closely related Vaccinia virus strain modified vaccinia Ankara make it one of the smallest non-vector derived orthopoxvirus genomes to be reported. Annotation of the assembly revealed a previously unreported signal peptide at the start of protein A38 and several predicted signal peptides that were found to differ from those previously described. Putative epigenetic modifications around various motifs have been identified and the assembly confirmed previous work showing the vaccine genome to most closely resemble that of Vaccinia virus strain Modified Vaccinia Ankara. The pipeline may be used for other DNA viruses, increasing the understanding of DNA virus evolution, virulence, host preference, and epigenomics.

Paleovirology of the DNA viruses of eukaryotes

Paleovirology is the study of ancient viruses and how they have coevolved with their hosts. An increasingly detailed understanding of the diversity, origins, and evolution of the DNA viruses of eukaryotes has been obtained through the lens of paleovirology in recent years. Members of multiple viral families have been found integrated in the genomes of eukaryotes, providing a rich fossil record to study. These elements have extended our knowledge of exogenous viral diversity, host ranges, and the timing of viral evolution, and are revealing the existence of entire new families of eukaryotic integrating dsDNA viruses and transposons. Future work in paleovirology will continue to provide insights into antiviral immunity, viral diversity, and potential applications, and reveal other secrets of the viral world.

Ancient Gene Capture and Recent Gene Loss Shape the Evolution of Orthopoxvirus-Host Interaction Genes

The survival of viruses depends on their ability to resist host defenses and, of all animal virus families, the poxviruses have the most antidefense genes. Orthopoxviruses (ORPV), a genus within the subfamily Chordopoxvirinae, infect diverse mammals and include one of the most devastating human pathogens, the now eradicated smallpox virus. ORPV encode ∼200 genes, of which roughly half are directly involved in virus genome replication and expression as well as virion morphogenesis. The remaining ∼100 "accessory" genes are responsible for virus-host interactions, particularly counter-defense of innate immunity. Complete sequences are currently available for several hundred ORPV genomes isolated from a variety of mammalian hosts, providing a rich resource for comparative genomics and reconstruction of ORPV evolution. To identify the provenance and evolutionary trends of the ORPV accessory genes, we constructed clusters including the orthologs of these genes from all chordopoxviruses. Most of the accessory genes were captured in three major waves early in chordopoxvirus evolution, prior to the divergence of ORPV and the sister genus Centapoxvirus from their common ancestor. The capture of these genes from the host was followed by extensive gene duplication, yielding several paralogous gene families. In addition, nine genes were gained during the evolution of ORPV themselves. In contrast, nearly every accessory gene was lost, some on multiple, independent occasions in numerous lineages of ORPV, so that no ORPV retains them all. A variety of functional interactions could be inferred from examination of pairs of ORPV accessory genes that were either often or rarely lost concurrently. IMPORTANCE Orthopoxviruses (ORPV) include smallpox (variola) virus, one of the most devastating human pathogens, and vaccinia virus, comprising the vaccine used for smallpox eradication. Among roughly 200 ORPV genes, about half are essential for genome replication and expression as well as virion morphogenesis, whereas the remaining half consists of accessory genes counteracting the host immune response. We reannotated the accessory genes of ORPV, predicting the functions of uncharacterized genes, and reconstructed the history of their gain and loss during the evolution of ORPV. Most of the accessory genes were acquired in three major waves antedating the origin of ORPV from chordopoxviruses. The evolution of ORPV themselves was dominated by gene loss, with numerous genes lost at the base of each major group of ORPV. Examination of pairs of ORPV accessory genes that were either often or rarely lost concurrently during ORPV evolution allows prediction of different types of functional interactions.

The Kaumoebavirus LCC10 Genome Reveals a Unique Gene Strand Bias among "Extended Asfarviridae"

Kaumoebavirus infects the amoeba Vermamoeba vermiformis and has recently been described as a distant relative of the African swine fever virus. To characterize the diversity and evolution of this novel viral genus, we report here on the isolation and genome sequencing of a second strain of Kaumoebavirus, namely LCC10. Detailed analysis of the sequencing data suggested that its 362-Kb genome is linear with covalently closed hairpin termini, so that DNA forms a single continuous polynucleotide chain. Comparative genomic analysis indicated that although the two sequenced Kaumoebavirus strains share extensive gene collinearity, 180 predicted genes were either gained or lost in only one genome. As already observed in another distant relative, i.e., Faustovirus, which infects the same host, the center and extremities of the Kaumoebavirus genome exhibited a higher rate of sequence divergence and the major capsid protein gene was colonized by type-I introns. A possible role of the Vermamoeba host in the genesis of these evolutionary traits is hypothesized. The Kaumoebavirus genome exhibited a significant gene strand bias over the two-third of genome length, a feature not seen in the other members of the “extended Asfarviridae” clade. We suggest that this gene strand bias was induced by a putative single origin of DNA replication located near the genome extremity that imparted a selective force favoring the genes positioned on the leading strand.

Comparative Analysis of the Circular and Highly Asymmetrical Marseilleviridae Genomes

Marseilleviridae members are large dsDNA viruses with icosahedral particles 250 nm in diameter infecting Acanthamoeba. Their 340 to 390 kb genomes encode 450 to 550 protein-coding genes. Since the discovery of marseillevirus (the prototype of the family) in 2009, several strains were isolated from various locations, among which 13 are now fully sequenced. This allows the organization of their genomes to be deciphered through comparative genomics. Here, we first experimentally demonstrate that the Marseilleviridae genomes are circular. We then acknowledge a strong bias in sequence conservation, revealing two distinct genomic regions. One gathers most Marseilleviridae paralogs and has undergone genomic rearrangements, while the other, enriched in core genes, exhibits the opposite pattern. Most of the genes whose protein products compose the viral particles are located in the conserved region. They are also strongly biased toward a late gene expression pattern. We finally discuss the potential advantages of Marseilleviridae having a circular genome, and the possible link between the biased distribution of their genes and the transcription as well as DNA replication mechanisms that remain to be characterized.

Long-read assays shed new light on the transcriptome complexity of a viral pathogen

Characterization of global transcriptomes using conventional short-read sequencing is challenging due to the insensitivity of these platforms to transcripts isoforms, multigenic RNA molecules, and transcriptional overlaps. Long-read sequencing (LRS) can overcome these limitations by reading full-length transcripts. Employment of these technologies has led to the redefinition of transcriptional complexities in reported organisms. In this study, we applied LRS platforms from Pacific Biosciences and Oxford Nanopore Technologies to profile the vaccinia virus (VACV) transcriptome. We performed cDNA and direct RNA sequencing analyses and revealed an extremely complex transcriptional landscape of this virus. In particular, VACV genes produce large numbers of transcript isoforms that vary in their start and termination sites. A significant fraction of VACV transcripts start or end within coding regions of neighbouring genes. This study provides new insights into the transcriptomic profile of this viral pathogen.

Novel classes of replication-associated transcripts discovered in viruses

The role of RNA molecules in the priming of DNA replication and in providing a template for telomerase extension has been known for decades. Since then, several transcripts have been discovered, which play diverse roles in governing replication, including regulation of RNA primer formation, the recruitment of replication origin (Ori) recognition complex, and the assembly of replication fork. Recent studies on viral transcriptomes have revealed novel classes of replication-associated (ra)RNAs, which are expressed from the genomic locations in close vicinity to the Ori. Many of them overlap the Ori, whereas others are terminated close to the replication origin. These novel transcripts can be both protein-coding and non-coding RNAs. The Ori-overlapping part of the mRNAs is generally either the 5ʹ-untranslated regions (UTRs), or the 3ʹ-UTRs of the longer isoforms. Several raRNAs have been identified in various viral families using primarily third-generation long-read sequencing. Hyper-editing of these transcripts has also been described.

Construction of an infectious horsepox virus vaccine from chemically synthesized DNA fragments

Edward Jenner and his contemporaries believed that his variolae vaccinae originated in horses and molecular analyses show that modern vaccinia virus (VACV) strains share common ancestry with horsepox virus (HPXV). Given concerns relating to the toxicity of modern VACV vaccines, we asked whether an HPXV-based vaccine might provide a superior alternative. Since HPXV may be extinct and the only specimen of HPXV that has been identified is unavailable for investigation, we explored whether HPXV could be obtained by large-scale gene synthesis. Ten large (10-30 kb) fragments of DNA were synthesized based on the HPXV sequence along with two 157 nt VACV terminal sequences, and were recombined into a live synthetic chimeric HPXV (scHPXV) in cells infected with Shope fibroma virus (SFV). Sequencing of the 212 kbp scHPXV confirmed it encoded a faithful copy of the input DNA. We believe this is the first complete synthesis of a poxvirus using synthetic biology approaches. This scHPXV produced smaller plaques, produced less extracellular virus and exhibited less virulence in mice than VACV, but still provided vaccine protection against a lethal VACV challenge. Collectively, these findings support further development of scHPXV as a novel replication-proficient smallpox vaccine.

Identification of Vaccinia Virus Replisome and Transcriptome Proteins by Isolation of Proteins on Nascent DNA Coupled with Mass Spectrometry

Poxviruses replicate within the cytoplasm and encode proteins for DNA and mRNA synthesis. To investigate poxvirus replication and transcription from a new perspective, we incorporated 5-ethynyl-2'-deoxyuridine (EdU) into nascent DNA in cells infected with vaccinia virus (VACV). The EdU-labeled DNA was conjugated to fluor- or biotin-azide and visualized by confocal, superresolution, and transmission electron microscopy. Nuclear labeling decreased dramatically after infection, accompanied by intense labeling of cytoplasmic foci. The nascent DNA colocalized with the VACV single-stranded DNA binding protein I3 in multiple puncta throughout the interior of factories, which were surrounded by endoplasmic reticulum. Complexes containing EdU-biotin-labeled DNA cross-linked to proteins were captured on streptavidin beads. After elution and proteolysis, the peptides were analyzed by mass spectrometry to identify proteins associated with nascent DNA. The known viral replication proteins, a telomere binding protein, and a protein kinase were associated with nascent DNA, as were the DNA-dependent RNA polymerase and intermediate- and late-stage transcription initiation and elongation factors, plus the capping and methylating enzymes. These results suggested that the replicating pool of DNA is transcribed and that few if any additional viral proteins directly engaged in replication and transcription remain to be discovered. Among the host proteins identified by mass spectrometry, topoisomerases IIα and IIβ and PCNA were noteworthy. The association of the topoisomerases with nascent DNA was dependent on expression of the viral DNA ligase, in accord with previous proteomic studies. Further investigations are needed to determine possible roles for PCNA and other host proteins detected.IMPORTANCE Poxviruses, unlike many well-characterized animal DNA viruses, replicate entirely within the cytoplasm of animal cells, raising questions regarding the relative roles of viral and host proteins. We adapted newly developed procedures for click chemistry and iPOND (Isolation of proteins on nascent DNA) to investigate vaccinia virus (VACV), the prototype poxvirus. Nuclear DNA synthesis ceased almost immediately following VACV infection, followed swiftly by the synthesis of viral DNA within discrete cytoplasmic foci. All viral proteins known from genetic and proteomic studies to be required for poxvirus DNA replication were identified in the complexes containing nascent DNA. The additional detection of the viral DNA-dependent RNA polymerase and intermediate and late transcription factors provided evidence for a temporal coupling of replication and transcription. Further studies are needed to assess the potential roles of host proteins, including topoisomerases IIα and IIβ and PCNA, which were found associated with nascent DNA.

Cytoplasmic ATR Activation Promotes Vaccinia Virus Genome Replication

In contrast to most DNA viruses, poxviruses replicate their genomes in the cytoplasm without host involvement. We find that vaccinia virus induces cytoplasmic activation of ATR early during infection, before genome uncoating, which is unexpected because ATR plays a fundamental nuclear role in maintaining host genome integrity. ATR, RPA, INTS7, and Chk1 are recruited to cytoplasmic DNA viral factories, suggesting canonical ATR pathway activation. Consistent with this, pharmacological and RNAi-mediated inhibition of canonical ATR signaling suppresses genome replication. RPA and the sliding clamp PCNA interact with the viral polymerase E9 and are required for DNA replication. Moreover, the ATR activator TOPBP1 promotes genome replication and associates with the viral replisome component H5. Our study suggests that, in contrast to long-held beliefs, vaccinia recruits conserved components of the eukaryote DNA replication and repair machinery to amplify its genome in the host cytoplasm.

Isolation and sequencing of active origins of DNA replication by nascent strand capture and release (NSCR)

Nascent strand capture and release (NSCR) is a method for isolation of short nascent strands to identify origins of DNA replication. The protocol provided involves isolation of total DNA, denaturation, size fractionation on a sucrose gradient, 5'-biotinylation of the appropriate size nucleic acids, binding to a streptavidin coated magnetic beads, intensive washing, and specific release of only the RNA-containing chimeric nascent strand DNA using ribonuclease I (RNase I). The method has been applied to mammalian cells derived from proliferative tissues and cell culture but could be used for any system where DNA replication is primed by a small RNA resulting in chimeric RNA-DNA molecules.