Parvovirus minute virus of mice interacts with sites of cellular DNA damage to establish and amplify its lytic infection
Abstract
We have developed a generally adaptable, novel high-throughput Viral Chromosome Conformation Capture assay (V3C-seq) for use in trans that allows genome-wide identification of the direct interactions of a lytic virus genome with distinct regions of the cellular chromosome. Upon infection, we found that the parvovirus Minute Virus of Mice (MVM) genome initially associated with sites of cellular DNA damage that in mock-infected cells also exhibited DNA damage as cells progressed through S-phase. As infection proceeded, new DNA damage sites were induced, and virus subsequently also associated with these. Sites of association identified biochemically were confirmed microscopically and MVM could be targeted specifically to artificially induced sites of DNA damage. Thus, MVM established replication at cellular DNA damage sites, which provide replication and expression machinery, and as cellular DNA damage accrued, virus spread additionally to newly damaged sites to amplify infection. MVM-associated sites overlap significantly with previously identified topologically-associated domains (TADs).
https://doi.org/10.7554/eLife.37750.001eLife digest
Viruses are small infectious particles that can only reproduce with the help of a host. Once they are inside their victim, they hijack the cells’ genetic material and reprogram it to become a virus factory that produces more virus particles. Parvoviruses, for example, are among the simplest of viruses and need all resources a cell has to offer to successfully replicate.
This process often takes place at so-called replication centers that contain these necessary factors. It was previously thought that parvoviruses set up such centers randomly, and gather the required molecules such as proteins to these sites. However, it was not well understood how they do this.
Now, Majumder et al. have developed a new method that enabled them to study in detail how parvoviruses gain access to the resources of the cell they need to initiate and amplify replication. The results show that parvoviruses set up their replication centers at sites on the host DNA that are already rich in proteins needed to repair and then replicate damaged DNA. Some of these sites already exist in the cell’s genetic material as a consequence of naturally occurring processes, but others are created during infection by the virus. These findings may have important implications for how other viruses may establish their replication.
Viruses, including parvoviruses, are important pathogens. Like many microbes, viruses can be beneficial for our health and environment. Others, however, can be harmful. A clearer understanding of how viruses establish and amplify an infection may provide new treatment opportunities.
https://doi.org/10.7554/eLife.37750.002Introduction
DNA viruses that replicate in the nucleus depend on host cellular functions for transcriptional and replication machinery to express and amplify their genomes. Accessing these functions is critical to productive infection, yet successful establishment of replication must also overcome cellular antiviral activity, which for larger DNA viruses includes innate immune responses, epigenetic silencing, the cellular DNA-damage response (DDR), and antiviral activity found associated with PML bodies (Weitzman et al., 2010). Replication of many DNA viruses takes place in distinct micro-nuclear compartments termed replication centers that are rich in factors viruses must interact with - either positively or negatively - to productively replicate (Schmid et al., 2014). It is not fully clear, however, how nuclear-replicating viruses initiate replication centers in order to optimize access to factors and functions they need to either utilize or inactivate.
Parvoviruses are small non-enveloped icosahedral viruses that are important pathogens in many animal species including humans. Minute Virus of Mice (MVM) is an autonomously replicating parvovirus that is lytic in murine cells and transformed human cells (Cotmore and Tattersall, 2014). The viral genome is approximately 5 kb and possesses inverted terminal repeats at each end that serve as origins of replication (Cotmore and Tattersall, 2014). MVM encodes two non-structural proteins: the larger non-structural phosphoprotein NS1 performs a number of functions required for viral replication, while NS2 plays important, currently undefined, roles during infection of the normal murine host (Cotmore and Tattersall, 2014).
Parvoviruses are the only known viruses of vertebrates that contain single-stranded linear DNA genomes, and thus, they present novel replicative DNA structures to cells during infection. They depend heavily on cellular functions for replication, and unlike the DNA tumor viruses, do not drive quiescent cells into S-phase. However, following S-phase entry, cellular DNA polymerase δ converts the single stranded viral DNA genome into a double stranded molecule that serves as a template for transcription of the viral genes (Cotmore and Tattersall, 2013). As MVM infection progresses through S-phase, it induces substantial cellular DNA damage and evokes a robust, ATM-dependent DNA damage response (DDR, [Adeyemi et al., 2010]). Infection is characterized by a pre-mitotic cell cycle arrest that is both p21 and CHK1 independent (Adeyemi and Pintel, 2012; 2014). During this block virus replication proceeds for many hours, and ATM inhibitors reduce ongoing viral replication (Adeyemi et al., 2010).
Parvoviruses establish replication factories in the nucleus (termed Autonomous Parvovirus-Associated Replication, or APAR, bodies) where active transcription of viral genes and viral replication takes place (Bashir et al., 2001). Low resolution confocal microscopy originally showed that DDR sensor and response proteins, cell cycle regulators, DNA polymerases and RNA polymerase II accumulate in MVM APAR bodies where they co-localize with replicating viral DNA and NS1 (Adeyemi et al., 2010; Bashir et al., 2001; Ruiz et al., 2011). More recently, higher magnification microscopic studies revealed that phosphorylated histone H2A variant γ-H2AX, as well as other DDR proteins ([Ruiz et al., 2011] see Figure 1, below), seemed to reside adjacent to, rather than within, early APAR bodies. This gave rise to the suggestion that these DDR factors may reside on cellular DNA in the vicinity of viral replication centers (Ruiz et al., 2011). The induction of cellular DNA damage leads to almost instantaneous recruitment of DDR factors to the break site which coordinate complex signaling cascades that recruit DNA damage sensors, repair mediators and effectors (Hashiguchi et al., 2007; Polo and Jackson, 2011). Endogenous sites of damage are also often associated with interference between replication and transcription polymerases (Durkin and Glover, 2007). DNA breaks, therefore, serve as cellular depots of DDR proteins, and factors involved in DNA replication and expression (Hashiguchi et al., 2007; Polo and Jackson, 2011). While it is plausible that during infection damaged cellular DNA is relocated to sites of MVM replication, it seemed possible, alternatively, that the virus initially established its replication centers at damaged cellular sites where these DDR, replication, and expression factors were already present.
It has been suggested that a number of viruses, including hepatitis B virus (HBV) and human papillomavirus (HPV), associate with sites of DNA damage, including early-replicating or common, fragile sites (ERFs, and CFSs, respectively) at early times during their infections (Durkin and Glover, 2007; Jang et al., 2014; Tubbs and Nussenzweig, 2017). The HPV genome was shown to be tethered to CFSs by the chromatin modifier BRD4, which facilitates subsequent integration into the host genome utilizing the cellular DDR machinery (Feitelson and Lee, 2007; Jang et al., 2014). These studies utilized chromatin immunoprecipitation (ChIP) assays of the HPV E2 protein to identify sites of HPV localization to cellular CFSs, which were then validated by 3D-FISH assays of the viral and cellular DNA (Jang et al., 2014). Consistent with these findings, crosslinked-ChIP assays have recently demonstrated that FANCD2, required for maintaining fragile site stability and coordinating their replication, also associates with HPV genomes at replication centers (Madireddy et al., 2016; Spriggs and Laimins, 2017). However, an unbiased way to map the interaction between the viral and cellular genomes directly has been largely unavailable until recently.
The development of Chromosome Conformation Capture (3C) technologies has enabled detailed analysis of cis-interactions between separated regions of the genome, as well as aspects of chromatin packaging (Dixon et al., 2016). When combined with high-throughput sequencing, these techniques (termed 4C, 5C, and Hi-C assays) have become valuable tools for studying the details of nuclear configuration (Denker and de Laat, 2016). Specifically, 4C and Hi-C assays provide genome-wide interaction data through unbiased deep sequencing. Hi-C assays are designed to provide information on all nuclear interactions whereas 4C assays enable higher resolving power and a deeper interrogation of genomic associations by utilizing inverse PCRs to amplify only the ‘bait’ (Lajoie et al., 2015). In this study, we report the novel adaptation of high-throughput circular chromosome conformation capture assay (4C) for use in trans which we term V3C (Viral Chromosome Conformation Capture). This assay, which should be generally adaptable, has allowed us to characterize, on a genome-wide scale, the direct association of the linear MVM genome with discreet regions of the cellular genome. These sites, termed Virus Association Domains (VADs), correlated initially with sites of cellular DNA damage that in mock-infected cells also exhibited damage as cells progressed through S-phase. As infection progressed these sites expanded, and additional sites of DNA damage were induced. MVM subsequently associated with the newly induced sites as infection amplified. Sites of association identified biochemically were confirmed microscopically, and in addition, MVM could be targeted specifically to sites of DNA damage artificially engineered into the cellular chromosome. Development of the V3C assay has allowed us to suggest the following model. Soon after nuclear entry MVM homes first to sites of pre-existing endogenous DNA damage to initiate infection at sites that provide cellular factors necessary for its replication. Subsequently, as cellular DNA damage accrues, virus spreads additionally to these sites of damage to amplify infection. The V3C-seq assay should be useful for characterizing the interaction of many DNA viruses that associate with the cellular genome, and provide a useful tool to characterize the molecular events leading to the initiation of infection.
Results
The replicating MVM genome localized adjacent to regions of the cellular genome undergoing a DDR
As MVM infection progresses through S-phase, it induces cellular DNA damage and evokes a robust, ATM-dependent DNA damage response (DDR) characterized by a pre-mitotic cell cycle arrest that is both p21 and CHK1 independent (Adeyemi and Pintel, 2012, 2014). During this block, virus replication proceeds for many hours, and ATM inhibitors reduce ongoing viral replication (Adeyemi et al., 2010). Standard confocal microscopy demonstrated previously that numerous cell cycle and DDR effector proteins, RNA polymerase II, as well as DNA polymerase-α and δ are found associated with MVM replication centers called APAR bodies (Adeyemi et al., 2010; Bashir et al., 2000; Kollek et al., 1982; Ruiz et al., 2011). Similar examples of such confocal images of DDR proteins associated with APAR bodies, but not the irrelevant transcription factor NR5A2 (Duggavathi et al., 2008), processed with deconvolution, can be seen in Figure 1A. Comparison of the median three-dimensional distance between FANCD2 and γ-H2AX - which localize to stalled replication forks where they facilitate DNA repair (Kim et al., 2018; Lossaint et al., 2013; Madireddy et al., 2016), and the center of APAR bodies - identified by MVM NS1 staining, indicated that at 16 hr post infection (hpi) and release (representing approximately 8–10 hr into S-phase in our para-synchronization protocol), they localized closely with MVM replication centers (Figure 1B). This was in contrast to the irrelevant transcription factor NR5A2 (Duggavathi et al., 2008) which exhibited a diffuse localization (Figure 1B). However, confocal super resolution imaging (Airyscan) demonstrated that γ-H2AX seemingly localized to the periphery of APAR bodies (Figure 1C). Super-resolution imaging using the GSD-STORM platform also demonstrated that γ-H2AX, and FANCD2 (Figure 1D), localized to the periphery of APAR bodies (Figure 1D). γ-H2AX characteristically amplifies on damaged DNA as megabase (Mb)-sized platforms, making the possibility of marking the 5 kb MVM genome with phosphorylated histone H2AX less likely (Rogakou et al., 1999). These results suggested that, MVM replication centers may localize to and expand adjacent to sites of cellular chromatin undergoing DNA damage, where replication, expression and DDR factors reside.
The MVM genome associated directly with discrete sites on the cellular genome
To characterize the association of the MVM genome with the cellular genome during lytic infection more directly we have developed a high-throughput chromosome conformation capture (3C) assay for use in trans. 3C assays have been typically used to identify long-range interactions between regions of a single chromosome (Dekker et al., 2013). Our analysis, which we term V3C-seq (Viral Chromosome Conformation Capture Sequencing) allowed us to identify direct interactions between the linear MVMp genome and the cellular chromosome in an infected cell population on a genome-wide scale in an unbiased manner. V3C-seq assays utilize formaldehyde-mediated crosslinking to first ‘freeze’ the locations of the viral and cellular genomes at various points during infection. Samples are then digested and ligated under conditions that favor intramolecular interaction, and the resultant novel virus-cell DNA fragments are subjected to high-throughput sequencing (Figure 2A). The assay provides a precise genomic map of the sites with which viral DNA interacts, and the frequency with which unique individual linked fragments are generated provides quantification of these interactions.
V3C-seq was performed in parasynchronized mouse A9 fibroblasts, the traditional host for MVM, at various times post-infection. A typical time-course of MVM infection is shown in Figure 2—figure supplement 1A. Assays utilized a viral viewpoint at a HindIII site at nucleotide 2651 and NlaIII site at 1899 in the MVMp genome, thereby capturing the interaction of the MVMp fragment containing both the viral P4 and P38 promoters upstream of the HindIII site. Clustering algorithms and visualization of interaction sites on the UCSC Genome Browser (Kent et al., 2002; Ramírez et al., 2016) revealed that by 12 hr post-infection and release (representing approximately 4–6 hr into S-phase) MVM genomes associated with discrete regions on most cellular chromosomes [Figure 2—figure supplements 1B and 2, (Kent et al., 2002)]. These cellular sites served as initial amplification points for MVM, and upon progression to 16 hpi, the virus both expanded at these regions and associated with new sites (Figure 2B, Figure 2—figure supplement 2). We term these sites of association Virus Associated Domains, or VADs. The clusters of interacting sites on each chromosome ranged in density and in size (Figure 2—figure supplement 2), from approximately 1–2 Mb to larger 5–15 Mb-size domains (chromosomes 17 and 19 are shown in Figure 2B). While most chromosomes contained multiple small VADS, many contained 1–3 larger VADs of 5–15 Mb size. Chromosomes 1, 13, 18, X and Y had fewer discernible interaction sites. The larger VADs in Figure 2B are boxed for comparison purposes but is not meant to restrict the designation of VADs to a particular size.
Comparison of MVM interaction sites across the entire mouse genome showed that approximately 51% of VADs identified at 12 hpi were retained at 16 hpi, indicating that approximately 84% of VADs identified at 16 hpi were newly generated (Figure 2C). Approximately 39% of MVM interaction sites detected at 16 hpi were retained at 20 hpi, while only approximately 18% of the interaction sites identified at 20 hpi were newly generated (Figure 2C). Clustering analysis of MVM interaction sites in multiple replicates over the time-course of infection showed that they were reproducible across replicates (Figure 2—figure supplement 1B). Importantly, MVM interaction sites clustered together at 16 and 20 hpi in a characteristic manner that was distinct from early (12 hpi), and late (24 hpi) infection. These results suggested that interactions of MVM with the cellular chromosome increased as infection progressed. By 24 hpi, MVM interaction with the host chromosome was extensive (Figure 2—figure supplement 2). This latter observation indicated that MVM interaction with the cellular genome at VADs during early stages of infection was not an artifact of preferential sequencing at the VAD sites, and is consistent with previously published profiling of APAR bodies by microscopy (Ruiz et al., 2011). Additionally, at these late times, in the presence of saturating amount of MVM DNA, we observed interactions of the viral genome with sites throughout the mouse genome. This further suggested that VADs identified earlier during infection featured properties that enabled viral recruitment and replication. It is noteworthy, however, that the cellular genome undergoes substantial DNA damage by late stages of infection (Figure 3—figure supplement 2A and as described below), likely precluding detection of some interactions effectively by V3C-seq. In addition, viral packaging would be expected to reduce available viral genomes by approximately 20 hpi (Cotmore and Tattersall, 2014), which would also be predicted to contribute to the decreased interaction seen at 20 hpi in Figure 2B.
We chose to validate the V3C assay by confirming the association of MVM with one of these VADs, murine 19qA, using a focused Taqman-based assay, in which the frequency of novel ligation junctions were determined by quantitative PCR (qPCR). MVM association with the VAD at 19qA was readily detectable in our standard protocol using Taqman probes complementary to the MVM genome (forward) and the cellular genomic site (reverse) (Figure 2D). This association was substantially diminished, either in the absence of intramolecular ligation, or when cross-links were reversed prior to intramolecular ligation (Figure 2D). These experiments defined the lower limit of background levels generated by our V3C assay, and suggested that MVM associations with cellular VADs were specific, and mediated by DNA-DNA and/or DNA-protein intramolecular crosslinks.
Focused 3C-qPCR in parasynchronized NIH-3T3 fibroblasts infected with MVMp (Figure 2E), and EL4 lymphocyte cells infected with the lymphotrophic variant MVMi (Figure 2F) both showed association with a subset of the VADs identified in A9 cells (Figure 2B, Figure 2—figure supplement 2), which suggested a common mechanism may exist for the establishment of MVM replication sites in these mouse cell lines. Differential rates of MVM replication and cell cycle kinetics in A9 compared to NIH-3T3 and EL4 cells precluded performance of 3C-qPCR assays simultaneously in these systems. Together, we interpret our results as establishing that V3C-seq is a valid means to map the interaction of a lytic linear DNA virus with specific sites on the host cell genome.
MVM infection induced distinct sites of cellular DNA damage as demonstrated by ChIP-seq for γ-H2AX
As described above, super-resolution microscopy suggested that MVM replication centers seemingly associated adjacent to genomic sites containing factors involved in replication, expression, and the DDR. ERFs have such characteristics in uninfected cells, and as mentioned, MVM continues to induce DNA damage as infection proceeds (Adeyemi et al., 2010; Barlow et al., 2013). Therefore, as we identified sites of viral interaction with the cellular genome, we looked for association with sites of cellular DNA damage.
To identify sites of cellular DNA damage, we initially performed chromatin immunoprecipitation coupled with high-throughput sequencing (ChIP-seq, [Landt et al., 2012]) for γ-H2AX in parasynchronized A9 cells, either mock infected, infected with MVM, or mock infected and treated with hydroxyurea (HU). Results for algorithm-called peaks for chromosomes 17 and 19 are shown in Figure 3A, and for the complete murine genome is shown in Figure 3—figure supplement 1. Mock infected A9 cells (taken 12 hr post release, hpr) showed a significant number of sites of damage, as identified by γ-H2AX, as they passed into S-phase. These were likely ERFs, which accrue damage during replication (Figure 3A). Whole-genome peak analysis of all γ-H2AX bound regions revealed that approximately 55% of the sites identified in mock infected cells overlapped with those identified at 16 hr post-infection (hpi); approximately 55% of sites identified at 16 hpi were newly generated (Figure 3B, peak calling using EPIC and intersection using BEDtools, (Quinlan and Hall, 2010), as described in Materials and methods). By 16 hpi, sites of damage concentrated in distribution, and expanded in number (Figure 3A). At this point in infection MVM had begun to induce additional sites of DNA damage, as evident both by ChIP-seq and increased tail moments in Comet assays (Figure 3—figure supplement 2A). The majority of γ-H2AX-containing sites at 20 hpi were newly generated, coinciding with only 10% of the sites identified at 16 hpi, indicating the widespread induction of DNA damage by this point of infection. This can be seen more clearly in the magnified view of chromosomes 17 and 19 shown in Figure 3—figure supplement 2C. Interestingly, the γ-H2AX ChIP-seq regions identified at 16 hpi correlated well with γ-H2AX ChIP-seq performed following 12 hr treatment with HU (Figures 3A and 16 hpi vs HU). Approximately 51% of γ-H2AX peaks detected 16 hpi were shared with those induced after 12 hr treatment with HU, and conversely, approximately 26% of the peaks identified following treatment with HU were shared at 16 hpi. In order to confirm the statistical significance of the intersection analyses, the γ-H2AX peaks at indicated time points were intersected with randomly permuted peaks across the mouse genome and visualized as a Jaccard Plot (Figure 3B, far right).
The MVM genome initiated infection at sites of cellular DNA damage that in mock infected cells also exhibited DNA damage as the cells cycled through S-phase, and as infection progressed, localized to additional sites of induced damage.
Comparisons of the ChIP-seq results with V3C-seq assays showed that MVM associated directly with sites of cellular DNA damage, as identified by the presence of γ-H2AX at the same region, in a manner that increased as infection progressed. Figure 3A compares MVM VADs at 16 hpi, to sites of DNA damage (as determined by γ-H2AX ChIP-seq) for chromosomes 17 and 19 as infection progressed. Large VAD regions in Figure 3A are boxed for comparison purposes, but are not meant to restrict overlap only to VADs of that size. Comparisons for the full mouse genome are shown in Figure 3—figure supplement 1 and while there is significant variation, the overlap between VADs and sites positive for γ-H2AX ChIP-seq was strikingly consistent.
Figure 3C summarizes the genome-wide correlation at the nucleotide level of VADs and γ-H2AX ChIP-seq data presented in Figure 3—figure supplement 1. For the composite comparisons at 12, 16 and 20 hr post-infection, data was taken from the same experiment (comparing VADs at various time points as shown in Figure 2B to ChIP-seq sites at those times as shown in Figure 3A). At 12 hpi, MVM associated with approximately 55% of sites that in mock infected cells exhibited DNA damage upon progression into S-phase. By 16 hpi, this association rose to close to 80%.
By 16 hpi, close to 90% of γ-H2AX occupied sites overlapped with VADs (Figure 3C), which included γ-H2AX sites present in uninfected cells (Figure 3A). By the late time point of 20 hpi, approximately 70% of the γ-H2AX sites co-localized with VADs. Visualization of MVM association in the vicinity of γ-H2AX-positive sites using hierarchical clustering further revealed that MVM association with damaged sites increased from 16 hpi to 20 hpi; however, increased incidence of cellular DNA breaks after 20 hpi led to a decrease in the proportion of MVM-associated damaged sites (Figure 3D). Notably, approximately 25% of the VADs identified at 12 hpi, and approximately 95% of VADs identified at 16 hpi, associated with the γ-H2AX sites identified following 12 hr of treatment with HU. Randomly generated peaks showed less than 1% overlap with γ-H2AX peaks identified 16 hpi (Figure 3C). A magnified view of the large VADs at 19qA and 17qA/B outlined in Figure 3A is also provided in Figure 3—figure supplement 2C (left), while further magnifications of VAD regions demarcated by red rectangles in Figure 3—figure supplement 2C (left panel) at Narfl, Vwa7, Ehd1 and Slc29a2 genes are shown on the right panel.
The strong correlation of MVM interaction sites with sites that in uninfected cells exhibit DNA damage upon replication is consistent with the notion that MVM may have initially established replication at cellular fragile sites that are susceptible to DNA damage as cells cycle through S-phase, although it cannot be formally ruled out that these sites were virally-induced at the earliest times in infection. It is also important to note that VADs also correlated strongly with γ-H2AX ChIP-seq sites identified on cellular chromosomes of non-infected, HU-treated A9 cells (Figure 3A, red rectangles). This strongly implied, although does not prove, that the γ-H2AX identified associated with MVM during infection resides on cellular DNA. It is striking that sites of damage in uninfected cells, sites of damage induced by virus, and sites of damage induced by treatment of the DNA-damaging agent HU overlap so significantly.
Cellular sites of DNA damage also often contain BRCA1, which binds DNA and can co-localize with γ-H2AX in DNA double-strand break repair foci, although typically in a more narrow pattern (Barlow et al., 2013). As expected, VADs also strongly associated with sites identified by BRCA1 ChIP-seq at 16 hr post-infection (Figure 4A, row 2). Furthermore, VADs also overlapped with BRCA1 and γ-H2AX sites in primary mouse cells induced with replication stress agents (Figure 4—figure supplement 1A and B), and characterized as ERFs, in previously published studies (Barlow et al., 2013). As MVM can infect transformed human cells, we also performed focused V3C-qPCR at 16 hpi in parasynchronized SV40-transformed human NB324K cells. As shown in Figure 4—figure supplement 1C, MVM localized to the previously characterized human fragile site FRA5H, but not FRA11F. MVM also associated weakly at this time point with the prototypical human fragile site, FRA3B.
As an indirect confirmation of MVM recruitment to the cellular genome, we performed ChIP-seq for the MVM-NS1 protein, which binds covalently to the viral genome, and non-covalently to additional ACCAACCA consensus sequences throughout the MVM genome (Christensen et al., 1995). We reasoned that ChIP-seq assays for NS1 would confirm cellular sites associated with the viral DNA by secondary crosslinking of NS1-bound MVM DNA to cellular DNA. Figure 4A shows NS1 binding profiles to cellular chromosome 17 and 19 that are concordant with the VADs at 16 hpi, further validating our findings from V3C-seq assays. A genome-wide analysis of called peaks indicated that approximately 90% of the peaks identified by NS1-, BRCA1 -, and γ-H2AX ChIP-seq overlapped with VADs identified by V3C (Figure 4B and C), while overlap was undetectable when intersected with a randomly generated library of ChIP-seq peaks of equivalent size (Figure 4B). In concordance with these findings, the binding profile of NS1, BRCA1 and γ-H2AX around a VAD site at 16 hpi centered within 1 Mb of the MVM associated cellular site (Figure 4B). Taken together, our V3C-seq and ChIP-seq experiments are consistent with a model that upon infection, MVM first localized to cellular sites susceptible to DNA damage as cells progressed into S-phase, and as infection progressed, localized to additional sites of damage that were virally induced, to amplify its replication.
FISH assays confirmed that MVM localized with cellular sites of DNA damage
We next sought to confirm the association of MVM replication with sites of cellular DNA damage using super-resolution (STORM) microscopy. For these assays, we designed PCR-based FISH probes complementary to the MVM genome, to a VAD regions at 19qA, and to a control, VAD-negative, site at 6 pA (Figure 5A). 3D-FISH combined with confocal imaging of multiple nuclei was performed at 16 hpi. Representative examples are shown in Figure 5B and D, which demonstrated close localization between the MVM genome and 19qA-VAD probes (represented by red and green probes respectively), in contrast to the lack of direct localization of MVM with the control probe at chromosome 6 pA (Figure 5C and D; represented by red and cyan foci respectively). The 3D distances between the VAD probe and MVM genome in multiple nuclei were calculated using confocal imaging. As shown in Figure 5E, the median distance between the MVM genome and the 19qA-VAD and 15qE-VAD probes were approximately 0.7, and 0.6 µm, respectively. This is similar to the median radius of Type II APAR bodies, suggesting that on average VAD sites coincide with APAR bodies. In contrast, non-VAD control sites on chromosome 12 and 17 (12qA3 and 17qA2) exhibited a much greater range of co-location and were separated from its nearest MVM genome by median distances of 1.1 and 1.1 µm, respectively. Taken together, our representative super resolution imaging and quantitative 3D-FISH analyses support V3C results demonstrating that MVM localized with VADs.
MVM associated with artificially-engineered sites of DNA damage
If MVM preferentially associates with cellular sites of DNA damage, one might expect that MVM could be targeted to artificially–engineered sites of cellular DNA damage. We tested this in two ways. First, we used laser micro-irradiation of MVM-infected A9 cells at 18 hpi to induce focused cellular DNA damage, which is evident as a γ-H2AX ‘stripe’ in the nucleus (Figure 6A). Anti-NS1 staining of these cells suggested that MVM distinctly co-localized with irradiation-induced damaged cellular DNA, and in doing so, viral replication centers adapted to the shape of the damaged DNA stripe (Figure 6A, top two panels), rather than the distinct foci characteristic of APAR bodies (representative example shown in Figure 6A, bottom panel). Localization of the cellular transcription factor NR5A2, which is not found in MVM replication centers (Figure 1B), was not affected by micro-irradiation (Figure 6A, third panel).
These experiments demonstrated that at least NS1 localized to induced sites of DNA damage. The V3C-seq/NS1 ChIP-seq experiments described above suggested that NS1 was a useful surrogate for virus replication in our assays; however, we chose to further assess the direct interaction of the MVM genome with sites of artificially induced cellular DNA damage using directed cleavage of the genome by CRISPR/Cas9 (Sanjana et al., 2014). Guide RNAs were designed that targeted a gene desert in chromosome 9 (Figure 6B, cytogenetic location at 9qE1). These guides were transfected into A9 fibroblasts stably expressing CRISPR/Cas9. When these cells were then infected with MVM and assayed by focused 3C-qPCR, we detected substantially more amplicons between MVM and the DNA break site in cells transfected with 9qE1 desert-specific guide RNAs, compared to scrambled control guides (Figure 6C). Interaction was further confirmed using complementary TaqMan probes that recognize the chromosome 9 CRISPR cleavage site (Figure 6D). As expected, MVM interaction with a previously identified VAD on Chromosome 19 (19qA) was not significantly affected in this assay (Figure 6E). As an independent verification of the localization of MVM to this site, we also detected NS1 binding to the induced damage site using ChIP, suggesting that NS1 bound to the MVM genome at the break site is secondarily crosslinked and detected by ChIP-qPCR (Figure 6F). NS1 binding to the chromosome 19 VAD at 19qA was unaffected in these experiments (Figure 6G), consistent with our V3C findings in Figure 6C and D.
Discussion
When DNA viruses enter the nucleus they must locate to sites suitable to sustain replication. Small DNA viruses, such as parvoviruses, require multiple cellular factors for the expression and replication of their genomes. It may be that DNA viruses set up replication centers essentially randomly, and factors necessary for replication are recruited to these sites. An alternative model, suggested by the present work, is that incoming DNA viruses can initially locate to cellular sites that maintain factors necessary for virus replication. Sites of cellular DNA damage present such an opportunity (Hashiguchi et al., 2007; Polo and Jackson, 2011). Using a high-throughput conformational capture assay developed here for use in trans, we show that at early times post infection MVM interacted directly with sites of cellular DNA damage that in mock infected cells also exhibited DNA damage upon entry into S-phase. As infection progressed, interaction of the MVM genome, as well as the presence of the viral replication protein NS1, increased at these sites, suggesting they were sites of ongoing viral replication. ChIP-seq analysis for γ-H2AX and BRCA1 demonstrated that DNA damage increased during infection, and MVM subsequently associated also with newly induced sites. This model is supported by our findings that MVM NS1 and the MVM genome could be re-localized to artificially induced sites of cellular DNA damage. These observations support a model consistent with the notion that MVM initially establishes replication at cellular DNA damage sites that provide replication and expression machinery, as well as other DDR factors, and as infection progresses induces additional sites of cellular DNA damage, using these to amplify infection. It should be noted, however, that even at the earliest time points examined, multiple virus genomes have accumulated at sites of replication. Thus, specific localization of input genomes remains circumstantial. Consistent with our overall model, we have previously demonstrated that an ATM inhibitor applied during infection specifically reduced virus replication (Adeyemi et al., 2010). Also, as might be expected, we find that HU pretreatment of permissive rat F111 cells, which have lower levels of endogenous DNA damage as indicated by lower levels of γ-H2AX, resulted in increased MVM replication by 20 hpi (Figure 3—figure supplement 2B).
Because VADs correlated strongly with γ-H2AX ChIP-seq sites identified on mock-infected cells as they progressed through S-phase, and on non-infected, HU-treated A9 cells, the predominantly identified γ-H2AX signal associated with MVM during infection very likely resides on cellular DNA. However it remains possible that γ-H2AX and/or other DDR signaling proteins are directly associated with the viral genome during infection. In this regard, while there is a clear potential role in MVM replication for DNA polymerase-δ and gene expression factors potentially present at DNA damage sites, the possible roles of other DDR proteins in parvovirus replication warrants additional study.
Recently Shah and O’Shea have elegantly demonstrated a bipartite cellular DDR to adenovirus (Ad) infection that initially targets the virus while sparing the cell. The Mre11-Rad50-Nbs1 (MRN) complex first inhibits adenovirus replication without inducing a global response that could interfere with cellular proliferation or viability. As Ad overcomes this block, utilizing the Ad E4Orf6/E1B 55 kDa complex, and replicates to high levels, a global DDR is subsequently induced, and cellular DNA breaks were found to sequester DDR proteins from adenovirus preventing its replication (Shah and O'Shea, 2015). While MVM localizes to sites of cellular DNA damage as replication ensues, the model that we propose for parvovirus infection is different. It is consistent with the parvovirus life cycle, which depends upon the induction of a cell cycle arrest, is susceptible to ATM inhibitors, and which continually induces cellular DNA damage during its infection. Additionally, in contrast to adenovirus, MVM is much more dependent on host cell factors for its replication and expression (Cotmore and Tattersall, 2014), and it does not have an extensive genetic capacity to encode functions that inactivate the cellular DDR as does adenovirus (Ou et al., 2012; Querido et al., 2001; Stracker et al., 2002).
Chromosome conformation assays have been traditionally used to analyze the 3D folding principles of the cellular genome, including the regulation of promoter-enhancer loops and structural folding loops (Dixon et al., 2016). However, these assays can also serve as valuable tools to study the interaction between the host genome and invading virus. In the context of virus-infected cells, Chromatin Interaction Analysis by Paired-End Tag Sequencing (ChIA-PET) assays have enabled the comprehensive mapping of EBV interaction with the cellular genome (Jiang et al., 2017). These studies, which combine ChIP-seq with chromosome conformation capture techniques, utilize the proximal interaction of distally located DNA regions bound by shared protein elements. They have demonstrated that Epstein Barr Virus enhancers can regulate the expression of the cellular Myc oncogene in lymphoblastoid cells via long-range promoter-enhancer looping, thereby contributing the EBV-mediated cellular transformation (Jiang et al., 2017). However, mapping of the virus-host interactome by ChIA-PET experiments can be limited by the necessity of having a priori knowledge of the proteins mediating this interaction. The EBV interactome has also been characterized using in-situ Hi-C assays, which generate chromosome conformation maps that provide a snapshot of how every restriction enzyme site associates with every other site throughout the genome. These studies showed that the latent EBV episome associates with gene poor regions, but relocalizes to gene-rich regions of the genome upon reactivation (Moquin et al., 2017). While highly informative, these studies did not have the resolving power of the assays described here for MVM.
For MVM, a single inverse-PCR viewpoint was sufficient to map its interactome. Moreover, using inverse PCR using primers complementary to the viral genome to generate the sequencing library ensured that the V3C-seq assay detects only MVM-host hybrid junctions. This enabled higher resolving power and a deeper interrogation of genomic associations, allowing the detection of both frequent and infrequent MVM interaction sites. The drawback of V3C-seq, however, is that it does not assay changes in nuclear architecture in response to viral infection, which invariably is altered during parvovirus replication, particularly in late stages.
Looping interactions have previously been observed at distinct loci in genomic regions proximal to some integrated viral genomes. For example, the integrated Murine Leukemia Virus (MuLV) genome has beRen, BRen, Ben shown to influence the folding properties of its proximally located Myc promoter, which in turn contributes to cellular oncogenesis (Zhang et al., 2012). In other studies, it has been shown in cell line models that the integrated HIV genome associates distally with an uncharacterized chromosomal region to promote reactivation of latent HIV (Dieudonné et al., 2009). Chromosome conformation capture analyses focused on elucidating the topological conformation of the viral genome have also revealed the conformational structure adopted by the gammaherpesvirus Kaposi’s Sarcoma Associated Herpesvirus (KSHV), which forms distinct promoter-enhancer loops mediated by the cellular architectural proteins CTCF and Cohesin (Kang et al., 2011). In this study we have utilized chromosome conformation capture technology for the first time to map the trans interaction of a lytic virus with the cellular genome in a non-biased way. However, in addition to assaying the inter-chromosomal interaction cellular sites, long-range interactions of integrated viruses such as MuLV and HIV with distal elements, or TADs, may also be assayed using V3C-seq.
Inspection of VADs show that regions of MVM interaction vary greatly in size. The size of the peaks correlates with the number of interactions, and so indicate that some regions are more densely populated by virus than others. As can be seen in the whole-genome analysis, there is quite a variation in VADs across different chromosomes (Figure 2—figure supplement 2). While most chromosomes have multiple VADS, most have only a few larger VADs of 5–15 Mb. Whether the size of the VADs correlates to the success of infection at that site is not known. There are a number of chromosomes [chromosomes 1, 13, 14, 18, X and Y (Figure 2—figure supplement 2 and Figure 3—figure supplement 1)] that have few sites of damage, as assessed by γ-H2AX ChIP-seq, as well as few VADs, reinforcing their correlation.
VADs identified by V3C-seq showed striking overlap with sites of damage in infected cells, at sites of damage incurred in uninfected cells as they progressed into S-phase, but surprisingly, also at sites in uninfected cells treated with HU. This overlap suggested that there may be a predilection in the cellular chromosome for sites sensitive to the induction of damage for which MVM has an affinity. In this vein, we find that most VADs (and by implication sites of endogenous and infection-induced damage and damage induced by HU) also overlap strikingly with Topologically Associated Domains (TADs), determined for murine CH12 B-cells, by Hi-C assays (Figure 4—figure supplement 2A–D, [Rao et al., 2014]). Such studies, in combination with high-resolution imaging, have shown that the mammalian genome is folded into megabase-sized cis-interacting regions (TADs). TADs have been divided into two main compartments: A (primarily euchromatin-like), and B (primarily heterochromatin-like) (Rao et al., 2014). The A compartment is characterized by being gene dense, containing highly expressed genes, and active chromatin marks. The A1 subcompartment was seen to finish replicating at the beginning of S-phase, while the A2 subcompartment was described as continuing to replicate into the middle of S-phase. Overlaying published chromatin marks from murine 3T3 cells induced with the DNA damaging agent aphidicolin (Kraushaar et al., 2013) suggests that MVM associations may be with the A compartment (Figure 4—figure supplement 2E). This compartment of cellular DNA would provide the necessary machinery for the expression and replication of MVM within S-phase, and would be consistent with MVM initially associating with sites that emerge as ERFs during cellular replication, which would be predicted to replicate in compartment A.
The VADs identified in our study occupy multiple adjacent TADs. This suggests that several adjacent TADs may contain the necessary environment to support MVM expression and replication. The mechanisms that maintain TAD borders are elusive, but structural proteins such as CTCF and cohesin have been implicated in these processes (Dixon et al., 2016; Rao et al., 2014, 2017). Indeed, Canela et al., have shown that CTCF and RAD21 form the anchors to cis-interacting cellular loops that are susceptible to topoisomerase TOP2B mediated DNA breaks (Canela et al., 2017). The MVM NS1 protein, which also interacts at TADs, is a known DNA binding protein with double-stranded nickase activity; however, the mechanisms by which the viral genome associates with the cellular sites remains to be determined.
Our results suggest that host chromatin states may play a significant role in permissiveness to DNA damage, and thereby influence MVM localization for replication. Consistent with this model, the HPV E2 protein has been found associated with actively transcribing genes and active chromatin, presumably to facilitate expression of HPV through hijacking the host transcriptional machinery at these sites (Jang et al., 2009). Indeed, profiling of sites on the cellular genomic that undergo replication stress have shown that these regions are encased in a protective chromatin environment to facilitate efficient repair (Kim et al., 2018). Such regions, either pre-existing prior to infection, or induced by virus, may provide the supportive environment necessary for successful infection. The generally adaptable, high-throughput V3C-seq assay allows genome-wide identification of the direct associations of viral genomes with distinct regions of the cellular chromosome. It should be useful for characterizing the interaction of many DNA viruses that associate with the cellular genome, and thus provide a useful tool to characterize the molecular events leading to the initiation of their infections.
Materials and methods
Contact for reagent and resource sharing
Request a detailed protocolFurther information and requests for resources and reagents should be directed to and will be fulfilled by David Pintel (pinteld@missouri.edu).
Experimental models and subject details
Request a detailed protocolCell lines were cultured in 5 percent FBS-containing DMEM media (5 percent CO2 and 37 degrees Celsius). Murine EL4 cells were cultured in RPMI media with 5 percent FBS. Cell lines are routinely authenticated for mycoplasma contamination, and background levels of DNA damage detected by γ-H2AX. Further information on cell line authentication and parvovirus replication are available at ATCC and published studies (Tattersall and Bratton, 1983) respectively.
Cell lines, viruses and viral infection
Request a detailed protocolMale Murine A9, NIH-3T3, EL4 and human NB324K cells were propagated and wild-type MVMp and MVMi were produced as previously described (Adeyemi et al., 2010). Infection was carried out at a Multiplicity Of Infection (MOI) of 5 unless otherwise stated, leading to infection rates of 70–80% as detected by NS1 staining.
LentiCRISPRv2 A9 cells
Request a detailed protocolLentiCRISPRv2 plasmid was obtained from Addgene (plasmid# 52961, [Sanjana et al., 2014]) and pseudotyped viruses were generated in 1 × 106 293 T cells transfected with 1 µg of LentiCRISPRv2, 1 µg of HIV Gag/Pol and 1 µg VSV-G proteins using Lipo293D (SignaGen). Supernatant containing the lentivirus was collected at 48 hr post-transfection. Independent preparations of LentiCRISPRv2 lentivirus were used to transduce A9 cells with 750 µl of lentiviral supernatant for 48 hr before selecting cells in 1 µg/ml of puromycin for 10 days. The resulting polyclonal puromycin-selected LentiCRISPRv2 A9 cells were validated for Cas9 expression by western blot, and were utilized for induced DNA break assays (described below).
Methods details
Cell synchronization and drug treatments
Request a detailed protocolMurine A9 and human NB324K cells were parasynchronized in G0 phase of cell cycle by growing them in 5% FBS containing DMEM without isoleucine for 36–42 hr (as described previously, (Adeyemi et al., 2010). NIH-3T3 cells were parasynchronized in G0 phase of cell cycle by growing them in 0.5% FBS containing DMEM for 36 hr. All cells were released into complete media containing 5% FBS in DMEM, and infected with MVMp at the time of release. Entry into S phase of cell cycle occurs approximately 8–10 hr after release into complete media. 16 hr post infection thus represents approximately 8–10 hr of transit into S-phase. Virally infected cells were harvested at the indicated timepoints and processed for experiments.
Chromosome conformation capture (3C) assay
Request a detailed protocolChromosome Conformation Capture assays were performed using 107 cultured A9, EL4, NIH-3T3 and NB-324K cells. Briefly, samples were cross-linked in 2 percent formaldehyde for 10 min, before quenching them in 0.125 M glycine. Cells were lysed in NP40 lysis buffer (0.1% NP40, NaCl, Tris-HCl) and the resulting nuclei were resuspended in restriction enzyme buffer (NEB Buffer 2.1). The nuclei were permeabilized in 0.3% SDS for an hour, followed by sequestration of SDS in 2% Triton X-100. The samples were digested in 400U of Hind III restriction enzyme overnight. Digestion was continued with a further 300U of Hind III on the next day, before inactivating the enzyme with 1% SDS at 65°C. SDS was sequestered with 1% Triton X-100, and 3C chromatin was resuspended in 1.15X T4 DNA Ligase Reaction Buffer. 50U of T4 DNA Ligase was added to the samples. Intramolecular ligation was carried out at room temperature for 4 hr, before reversing the crosslinks and digesting protein at 65 degrees C overnight with Proteinase K. 3C DNA was purified by phenol:chloroform:isoamyl alcohol extraction, isopropanol precipitation and finally using a PCR purification kit. The 3C-DNA was eluted in 200 microliters of Buffer EB (Qiagen). Cross-linking efficiencies were measured using Taqman-qPCR assays with primers and probes shown in Supplementary file 1. Relative crosslinking between two distally located HindIII fragments was determined by the ratio of the novel ligation junction to that of nearest neighbor interaction on the Ercc3 locus, as described previously (Hagège et al., 2007).
Viral chromosome conformation capture sequencing (V3C-seq) assay
View detailed protocolV3C-seq assays were performed with Hind III as the primary restriction enzyme to digest cross-linked MVM infected A9 fibroblast chromatin. The Hind III-digested DNA was intramolecular-ligated using the 3C procedure, before resuspending in Buffer EB (100 µl, Qiagen). 3C-DNA was secondary-digested with Nla III (100U, overnight at 37°C), before being heat inactivated and circularized with 100U of T4 DNA Ligase at room temperature overnight in 6 ml of ligation reaction. The V3C samples were precipitated by phenol:chloroform extraction, precipitated in isopropanol, resuspended in Qiagen Buffer EB (100 µl), and. Inverse PCR was performed on the circularized DNA using primers within the Hind III - Nla III fragments on the MVM genome using inverse PCR primers described in Supplementary file 1. Inverse PCR products were diluted 1:100 in TE buffer and used as templates for nested inverse PCRs (described in Supplementary file 1), yielding V3C-seq DNA libraries. Sequencing libraries were prepared using the NEB Ultra Kit, and twelve samples were pooled per run for 75 base-pair single end sequencing using an Illumina Next Seq 500 sequencer.
V3C-seq analysis
View detailed protocolV3C-seq samples were trimmed and aligned to the mouse reference genome (mm10 build) using Bowtie2 (Langmead and Salzberg, 2012). The Biostrings package in RStudio was used to generate a genome-wide map of HindIII restriction fragments for the assignment of reads (Pagès et al., 2017). To compare between different timepoints, reads for each fragment were averaged and quantile normalized using preprocessCore package on RStudio (Bolstad, 2013). For visualization of the V3C-seq data, a running mean was calculated using a window size of five contiguous HindIII fragments (Medvedovic et al., 2013). Bioinformatic codes provided in Table 1.
Chromatin immunoprecipitation (ChIP) assay
Request a detailed protocolThe indicated cells were cross-linked with 1% formaldehyde for 10 mins at room temperature and then quenched with 0.125 M glycine. The cells were collected and lysed using a ChIP lysis buffer (1% SDS, 10 mM EDTA, 50 mM Tris-HCl, pH 8, protease inhibitors) for 20 min on ice. The lysates were sonicated using a Diagenode Bioruptor for 75 cycles (30 s on and 30 s off per cycle), before being incubated overnight at 4°C with the indicated antibodies bound to Protein A Dynabeads (Invitrogen), in ChIP dilution buffer (0.01% SDS, 1.1% Triton X-100, 1.2 mM EDTA, 16.7 mM Tris-HCl pH8, 167 mM NaCl). Samples were washed for 3 min each at 4 degrees Celsius with low salt wash (0.01% SDS, 1% Triton X-100, 2 mM EDTA, 20 mM Tris-HCl pH8, 150 mM NaCl), high salt wash (0.01% SDS, 1% Triton X-100, 2 mM EDTA, 20 mM Tris-HCl pH8, 500 mM NaCl), lithium chloride wash (0.25M LiCl, 1% NP40, 1% DOC, 1 mM EDTA, 10 mM Tris-HCl pH8) and twice with TE buffer before being eluted with SDS- elution buffer (1% SDS, 0.1M Sodium bicarbonate). Following elution, the chromatin-antibody-DNA complexes, and the input chromatin were subjected to proteinase K treatment at 65°C overnight. The ChIP DNA was purified using a PCR purification kit (Qiagen), and eluted in 100 ul of Buffer EB (Qiagen). ChIP assays were analyzed by quantitative PCR (qPCR) with iTaq universal SYBR green mastermix (Bio-Rad), using primer sets described in Supplementary file 1, or sequenced as described below. Percent input was calculated as described previously (Fuller et al., 2017).
Sequencing libraries were generated from ChIP DNA using the NEBNext Ultra II Library Prep Kit for Illumina, and the sonication quality was determined using Agilent Bioanalyser. For ChIP-seq, twelve samples were pooled and sequenced on an Illumina Next Seq 500 using 75 base-pair Single End sequencing.
ChIP-seq analysis
Request a detailed protocolChIP-seq samples were aligned to the mouse genome (build mm10) using Bowtie2 (Langmead and Salzberg, 2012). Peaks were called with EPIC analysis software (using the SICER algorithm (Zang et al., 2009) according to default parameters. Called-peaks that were shared between replicates were identified using BEDtools software (Quinlan and Hall, 2010). Comparison between ChIP-seq and V3C-seq peaks were performed using Deeptools package (Ramírez et al., 2016). In order to compare the magnitudes of ChIP-seq peaks between different timepoints of MVM infection and mock versus Hydroxyurea treatment, rpm values were calculated (using Galaxy, [Afgan et al., 2016]) on the bedgraph files generated from EPIC, and were quantile normalized using preprocessCore package on RStudio (Bolstad, 2013). Bioinformatic codes provided in Table 1.
Laser Micro-Irradiation assays
Request a detailed protocolLaser micro-irradiation was performed on 1 million A9 cells cultured on glass bottom dishes (MatTek Corp.) infected with MVMp at an MOI of 10 for 18 hr. Cells were sensitized with 2 microliters of Hoechst dye (ThermoFisher Scientific) 5 min prior to irradiation. Samples were irradiated using a Leica TCP SP8 confocal microscope with a 405 nm laser using 25% power at 40 Hz frequency for 2 consecutive frames per field-of-view. Regions of interest (ROIs) were selected within the nucleus without traversing the nuclear membrane. Samples were processed for immunofluorescense imaging without CSK pre-extraction immediately after micro-irradiation.
CRISPR-Induced DNA break assays
Request a detailed protocolStable A9 cells expressing LentiCRISPRv2 were co-transfected with guide RNAs targeting chromosome 9 at 9qE1 (labelled as TGT), or scrambled control guide RNAs (labelled as CTRL) and human CD4 expressing vector during parasynchronization. CD4-positive cells were purified using an EasySep CD4+ T Cell Enrichment Kit (StemCell) prior to release into complete DMEM media and MVM infection. Infected cells were harvested and processed for ChIP and 3C assays at the indicated timepoints.
3D-FISH assays
Request a detailed protocolThe MVMp genome and indicated cellular regions were labelled with the DNA FISH-Tag Multicolor Kit (ThermoFisher). Briefly, 1 µg of DNA was labelled with aminoallyl-modified dNTP by nick-translation using the manufacturer’s instructions before being labelled with amine-modified Alexa-Fluor dyes (AlexaFluor 488 and AlexaFluor 555). The dye combinations were resuspended at equimolar amounts in hybridization buffer (50% formamide, 2X SSC, 40% dextran sulfate, 10% Denhardt’s solution) prior to hybridizing to the sample.
Parasynchronized MVMp-infected A9 cells were harvested at the indicated timepoints by pre-extracting with CSK Buffer (10 mM PIPES pH 6.8, 100 mM Sodium Chloride, 300 mM Sucrose, 1 mM EGTA, 1 mM Magnesium Chloride) for 3 min followed by CSK Buffer with 0.5% Triton for 3 min. Cells were crosslinked with 4% paraformaldehyde for 10 min at room temperature, before being washed with PBS. The nuclei were dehydrated by sequential treatments with 50% ethanol, 70% ethanol and 100% ethanol for 3 min each. Nuclei were subsequently rehydrated with sequential treatment with 70% ethanol, 50% ethanol and PBS for 3 min each. Cells and nuclei were permeabilized with 0.5% Triton X-100 in PBS for 15 min, before being washed with PBS. The samples were treated with 2 μg of RNAse A (Roche) in PBS for 1 hr at 37°C. Samples were denatured in 50% formamide (Ambion) before being hybridized to the suspended fluorescently labelled probes overnight at 37°C. Samples were washed in 2X SSC with 0.1% Tritox X-100 at 37°C followed by three times for 5 min each, followed by 2 washes in 2X SSC at 37°C for 5 min each. Samples were then mounted on slides and imaged on the indicated microscope.
For GSD/dSTORM super-resolution imaging, coverslips with adherent immunostained cells were mounted on cavity microscope slides with PBS (0.01 M, pH 7.4) imaging buffer containing 100 mM beta-mercaptoethylamine, 10% w/v glucose, 0.5 mg/ml glucose oxidase and 40 µg/ml catalase. GSD super-resolution imaging was performed on a Leica SR GSD 3D microscope (Leica Microsystems, Inc.) using a 560 nm (AlexaFluor 555) or a 632 nm (AlexaFluor 647) excitation lasers and a 405 nm back-pumping (activation) laser. A 160 × 1.43 NA oil-immersion objective lens was used for imaging. Two-color GSD images were acquired sequentially with an Andor iXon Ultra 897 EMCCD camera at exposure times 7–8 ms using a QGSD 561 quad filter cube and emission bandpasses 605/45 nm (AlexaFluor 555) and 695/85 nm (AlexaFluor 647). Approximately 8000 images per channel of a 18 × 18 µm field-of-view were acquired. The coordinates of single molecules were localized in all recorded raw images and high-resolution GSD images were constructed using Leica LAS X software (version 1.9).
For confocal imaging, samples were mounted on slides using Pro-Long Diamond anti-fade media with DAPI (Invitrogen). Confocal z-stacks were acquired using a Leica TCP SP8 confocal microscope with 488 nm (AlexaFluor 488) and 552 nm (AlexaFluor 555) excitation lasers and a 100 × 1.4 NA objective lens.
3D-FISH images were analyzed using ImageJ. Background noise was filtered out using the Kalman Stack Filter plugin to determine the coordinates (x,y,z) of the centers of the foci. The coordinates of the viral and cellular foci were measured using Sync Measure 3D. The 3D-distance was calculated by computing the displacement vector between the two locations as described previously (Shih and Krangel, 2010).
Immunofluorescense assays
Request a detailed protocolParasynchronized MVMp-infected A9 cells were harvested at the indicated time points processed as described above till permeabilization with 0.5% Triton X-100 in PBS. Samples were blocked with 3% BSA in PBS for 1 hr, incubated with the indicated antibodies for 1 hr, and incubated with the indicated secondary antibodies (tagged with Alexa Fluor fluorophores) for 1 hr. Samples were washed and mounted on slides with ProLong Diamond Antifade Mountant with DAPI (Invitrogen).
Alkaline comet assay
Request a detailed protocolAlkaline Comet Assays were performed using Trevigen Comet Assay kits. Murine A9 fibroblasts were grown on 10 centimeter dishes and mock infected, induced with Doxorubicin (200 nM) for 9 hr, or infected with MVMp at an MOI of 10 for 20 hr, before detaching them from the flask by scraping. Cells were washed with ice-cold PBS and resuspended at a density of 105 cells/ml in ice cold PBS. Cells were combined with molten LM Agarose at 37°C at a ratio of 1:10 and pipetted onto Comet Slides. Slides were placed at 4°C in the dark for 10 mins. Slides were immersed in 4°C Lysis solution for 30–60 min, before placing in Alkaline Unwinding Solution for 1 hr at 4°C in the dark. 850 ml Alkaline Electrophoresis Solution was added to the slide tray, and 21 Volts were applied for 30 min. Slides were immersed in water twice for 5 min each, followed by immersion in 70% ethanol for 5 min. Samples were dried at 37 degrees Celsius for 15 min, and subsequently stained with 100 ul of SYBR Gold for 30 min in the dark. Slides were briefly rinsed in water and completely dried at 37°C. Slides were imaged on a Leica widefield microscope.
Immunoblot analysis
Request a detailed protocolCells grown and infected in 60 mm dishes were harvested at the indicated timepoints, followed by lysis in modified RIPA buffer (20 mM Tris-HCl pH 7.5, 150 mM NaCl, 10% glycerol, 1% NP-40, 1% Sodium Deoxycholate, 0.1% SDS, 1 mM EDTA, 10 mM trisodium pyrophosphate, 20 mM Sodium Fluoride, 2 mM Sodium Orthovanadate and 1X Protease Inhibitor cocktail (Sigma). Protein concentrations were quantified using Bradford assay and equal amounts of lysates were loaded per well for Western blot analysis.
Southern blot analysis
Request a detailed protocolCells were grown on 25 mm plates and infected at an MOI of 5. Cells were harvested at the indicated timepoints, pelleted and resuspended in Southern Lysis Buffer. Cells were proteinase K treated overnight at 37°C, and sheared using 25 G X 5/8 inch 1 mL needle-syringe (BD Biosciences). Total DNA content was quantified using Nanodrop, equal amount of DNA loaded per well and electrophoresed on a 1 percent agarose gel. Samples were transferred to a nitrocellulose membrane and hybridized with completely homologous genomic clones.
Antibodies
Commercially available antibodies were used for ChIP assays and Immunofluorescence, and are described in the Antibody Table (Table 2) and Key Resources Table.
Plasmids
Lenti-CRISPRv2 plasmid was produced by Feng Zhang (Addgene plasmid 52961, [Sanjana et al., 2014]). pgRNA-humanized plasmid was produced by Stanley Qi (Addgene plasmid 44248, [Qi et al., 2013]). pCMV-CD4 was a gift from Dr. Marc Johnson (University of Missouri). Plasmids and reagents are available upon request.
Quantification and statistical analysis
Request a detailed protocolImaging studies (3D-FISH and Immunofluorescense) were quantified using ImageJ. Background noise was filtered out using the Kalman Stack Filter plugin, and the 3D distance between viral and cellular genome probes were calculated using Sync Measure 3D plugin to calculate the location of the center of mass between the imaged foci. The 3D-distance was calculated by computing the displacement vector between the two locations. The distances between foci were measured for multiple cells in preparations of viral infections, and were statistically analyzed using GraphPad Prism software. Statistical tests were performed using GraphPad Prism for imaging studies and chromosome conformation capture assays. The relevant statistical tests have been indicated in the respective figure legends. The code for bioinformatics analyses used to process V3C-seq and ChIP-seq data have been tabulated below:
Data availability
Request a detailed protocolThe V3C-seq and ChIP-seq data generated have been deposited in the Gene Expression Omnibus (GEO) under the accession codes GSE112957.
Data availability
Sequencing data have been deposited in GEO under accession code GSE112957.
-
Parvovirus Minute Virus of Mice Localizes to Sites of Cellular DNA Damage to Establish and Amplify its Lytic InfectionPublicly available at the NCBI Gene Expression Omnibus (accession no: GSE112957).
-
A three-dimensional map of the human genome at kilobase resolution reveals prinicples of chromatin loopingPublicly available at the NCBI Gene Expression Omnibus (accession no: GSE63525).
-
Genome-wide incorporation dynamics reveal distinct categories of turnover for the histone variant H3.3Publicly available at the NCBI Gene Expression Omnibus (accession no: GSE51505).
-
Genome-wide mapping of early replication fragile sites (ERFs)Publicly available at the NCBI Gene Expression Omnibus (accession no: GSE43504).
References
-
The ATR signaling pathway is disabled during infection with the parvovirus minute virus of miceJournal of Virology 88:10189–10199.https://doi.org/10.1128/JVI.01412-14
-
SoftwarepreprocessCore: A Collection of Pre-Processing Functions, version R package version, 1(0)preprocessCore: A Collection of Pre-Processing Functions.
-
Minute virus of mice transcriptional activator protein NS1 binds directly to the transactivation region of the viral P38 promoter in a strictly ATP-dependent mannerJournal of virology 69:5422–5430.
-
Parvovirus diversity and DNA damage responsesCold Spring Harbor Perspectives in Biology 5:a012989.https://doi.org/10.1101/cshperspect.a012989
-
Parvoviruses: small does not mean simpleAnnual Review of Virology 1:517–537.https://doi.org/10.1146/annurev-virology-031413-085444
-
Exploring the three-dimensional organization of genomes: interpreting chromatin interaction dataNature Reviews Genetics 14:390–403.https://doi.org/10.1038/nrg3454
-
The second decade of 3C technologies: detailed insights into nuclear organizationGenes & Development 30:1357–1382.https://doi.org/10.1101/gad.281964.116
-
Chromatin domains: the unit of chromosome organizationMolecular Cell 62:668–680.https://doi.org/10.1016/j.molcel.2016.05.018
-
Liver receptor homolog 1 is essential for ovulationGenes & Development 22:1871–1876.https://doi.org/10.1101/gad.472008
-
Chromosome fragile sitesAnnual Review of Genetics 41:169–192.https://doi.org/10.1146/annurev.genet.41.042007.165900
-
Minute virus of mice inhibits transcription of the cyclin B1 gene during infectionJournal of Virology 91:e00428–17.https://doi.org/10.1128/JVI.00428-17
-
Quantitative analysis of chromosome conformation capture assays (3C-qPCR)Nature Protocols 2:1722–1733.https://doi.org/10.1038/nprot.2007.243
-
Recruitment of DNA repair synthesis machinery to sites of DNA damage/repair in living human cellsNucleic Acids Research 35:2913–2923.https://doi.org/10.1093/nar/gkm115
-
The Epstein-Barr Virus Regulome in Lymphoblastoid CellsCell Host & Microbe 22:561–573.https://doi.org/10.1016/j.chom.2017.09.001
-
DNA polymerase requirements for Parvovirus H-1 DNA replication in vitroJournal of Virology 41:982–989.
-
ChIP-seq guidelines and practices of the ENCODE and modENCODE consortiaGenome Research 22:1813–1831.https://doi.org/10.1101/gr.136184.111
-
The Sequence Alignment/Map format and SAMtoolsBioinformatics 25:2078–2079.https://doi.org/10.1093/bioinformatics/btp352
-
The Epstein-Barr virus episome maneuvers between nuclear chromatin compartments during reactivationJournal of Virology 92:e01413–01417.https://doi.org/10.1128/JVI.01413-17
-
SoftwareBiostrings: Efficient manipulation of biological strings, version R Package Version, 2(0)Biostrings: Efficient manipulation of biological strings.
-
deepTools2: a next generation web server for deep-sequencing data analysisNucleic Acids Research 44:W160–W165.https://doi.org/10.1093/nar/gkw257
-
Megabase chromatin domains involved in DNA double-strand breaks in vivoThe Journal of Cell Biology 146:905–916.https://doi.org/10.1083/jcb.146.5.905
-
Improved vectors and genome-wide libraries for CRISPR screeningNature Methods 11:783–784.https://doi.org/10.1038/nmeth.3047
-
DNA virus replication compartmentsJournal of Virology 88:1404–1420.https://doi.org/10.1128/JVI.02046-13
-
NIH image to ImageJ: 25 years of image analysisNature Methods 9:671–675.https://doi.org/10.1038/nmeth.2089
-
Distinct contracted conformations of the Tcra/Tcrd locus during Tcra and Tcrd recombinationThe Journal of Experimental Medicine 207:1835–1841.https://doi.org/10.1084/jem.20100772
-
Reciprocal productive and restrictive virus-cell interactions of immunosuppressive and prototype strains of minute virus of miceJournal of Virology 46:944–955.
-
Genomes in conflict: maintaining genome integrity during virus infectionAnnual Review of Microbiology 64:61–81.https://doi.org/10.1146/annurev.micro.112408.134016
-
In Stem Cell Transcriptional Networks97–111, Spatial clustering for identification of ChIP-enriched regions (SICER) to map regions of histone methylation patterns in embryonic stem cells, In Stem Cell Transcriptional Networks, New York, Humana Press, 10.1007/978-1-4939-0512-6_5.
Article and author information
Author details
Funding
National Institute of Allergy and Infectious Diseases (AI046458)
- David J Pintel
National Institute of Allergy and Infectious Diseases (AI131468)
- Kinjal Majumder
National Institute of Allergy and Infectious Diseases (AI116595)
- David J Pintel
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank the members of the Pintel Lab for valuable discussions, Richard Adeyemi (Harvard University School of Medicine) for valuable suggestions, Andrew Huber (St. Jude Children’s Research Hospital) for helping to purify monoclonal anti-NS1 antibody from hybridoma supernatants, Lisa Burger for expert technical assistance and Alexander Jurkevich (University of Missouri Molecular Cytology Core) for assistance in GSD/dSTORM imaging. GSD/dSTORM imaging was supported by a research award from the University of Missouri Molecular Cytology Core to KM. High-throughput sequencing services were performed at the University of Missouri DNA Core Facility. Computational analysis was performed on the high performance computing infrastructure provided by Research Computing Support Services and in part by the National Science Foundation under grant number CNS-1429294 at the University of Missouri, Columbia MO. Work in the Pintel Laboratory is supported by NIH grants AI046458 and AI116595 to DJP. KM is supported by a Ruth L Kirschstein Postdoctoral Individual National Research Service Award AI131468.
Copyright
© 2018, Majumder et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 2,566
- views
-
- 312
- downloads
-
- 32
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Genetics and Genomics
- Microbiology and Infectious Disease
The sustained success of Mycobacterium tuberculosis as a pathogen arises from its ability to persist within macrophages for extended periods and its limited responsiveness to antibiotics. Furthermore, the high incidence of resistance to the few available antituberculosis drugs is a significant concern, especially since the driving forces of the emergence of drug resistance are not clear. Drug-resistant strains of Mycobacterium tuberculosis can emerge through de novo mutations, however, mycobacterial mutation rates are low. To unravel the effects of antibiotic pressure on genome stability, we determined the genetic variability, phenotypic tolerance, DNA repair system activation, and dNTP pool upon treatment with current antibiotics using Mycobacterium smegmatis. Whole-genome sequencing revealed no significant increase in mutation rates after prolonged exposure to first-line antibiotics. However, the phenotypic fluctuation assay indicated rapid adaptation to antibiotics mediated by non-genetic factors. The upregulation of DNA repair genes, measured using qPCR, suggests that genomic integrity may be maintained through the activation of specific DNA repair pathways. Our results, indicating that antibiotic exposure does not result in de novo adaptive mutagenesis under laboratory conditions, do not lend support to the model suggesting antibiotic resistance development through drug pressure-induced microevolution.
-
- Microbiology and Infectious Disease
Bacillus velezensis is a species of Bacillus that has been widely investigated because of its broad-spectrum antimicrobial activity. However, most studies on B. velezensis have focused on the biocontrol of plant diseases, with few reports on antagonizing Salmonella Typhimurium infections. In this investigation, it was discovered that B. velezensis HBXN2020, which was isolated from healthy black pigs, possessed strong anti-stress and broad-spectrum antibacterial activity. Importantly, B. velezensis HBXN2020 did not cause any adverse side effects in mice when administered at various doses (1×107, 1×108, and 1×109 CFU) for 14 days. Supplementing B. velezensis HBXN2020 spores, either as a curative or preventive measure, dramatically reduced the levels of S. Typhimurium ATCC14028 in the mice’s feces, ileum, cecum, and colon, as well as the disease activity index (DAI), in a model of infection caused by this pathogen in mice. Additionally, supplementing B. velezensis HBXN2020 spores significantly regulated cytokine levels (Tnfa, Il1b, Il6, and Il10) and maintained the expression of tight junction proteins and mucin protein. Most importantly, adding B. velezensis HBXN2020 spores to the colonic microbiota improved its stability and increased the amount of beneficial bacteria (Lactobacillus and Akkermansia). All together, B. velezensis HBXN2020 can improve intestinal microbiota stability and barrier integrity and reduce inflammation to help treat infection by S. Typhimurium.