SARS-CoV-2 RNA reverse-transcribed and integrated into the human genome
Prolonged SARS-CoV-2 RNA shedding and recurrence of PCR-positive tests have been widely reported in patients after recovery, yet these patients most commonly are non-infectious[1][1]–[14][2]. Here we investigated the possibility that SARS-CoV-2 RNAs can be reverse-transcribed and integrated into...
Summary
Prolonged SARS-CoV-2 RNA shedding and recurrence of PCR-positive tests have been widely reported in patients after recovery, yet these patients most commonly are non-infectious1–14. Here we investigated the possibility that SARS-CoV-2 RNAs can be reverse-transcribed and integrated into the human genome and that transcription of the integrated sequences might account for PCR-positive tests. In support of this hypothesis, we found chimeric transcripts consisting of viral fused to cellular sequences in published data sets of SARS-CoV-2 infected cultured cells and primary cells of patients, consistent with the transcription of viral sequences integrated into the genome. To experimentally corroborate the possibility of viral retro-integration, we describe evidence that SARS-CoV-2 RNAs can be reverse transcribed in human cells by reverse transcriptase (RT) from LINE-1 elements or by HIV-1 RT, and that these DNA sequences can be integrated into the cell genome and subsequently be transcribed. Human endogenous LINE-1 expression was induced upon SARS-CoV-2 infection or by cytokine exposure in cultured cells, suggesting a molecular mechanism for SARS-CoV-2 retro-integration in patients. This novel feature of SARS-CoV-2 infection may explain why patients can continue to produce viral RNA after recovery and suggests a new aspect of RNA virus replication.Introduction
Continuous or recurrent positive SARS-CoV-2 PCR tests have been reported in patients weeks or months after recovery from an initial infection1–14. Although bona fide re-infection of SARS-CoV-2 after recovery has been reported lately15, cohort-based studies with strict quarantine on subjects recovered from COVID-19 suggested “re-positive” cases were not caused by re-infection16,17. Furthermore, no replication-competent virus was isolated or spread from these PCR-positive patients1–3,5,6,12. The cause for such prolonged and recurrent viral RNA production is unknown. As positive-stranded RNA viruses, SARS-CoV-2 and other beta-coronaviruses such as SARS-CoV-1 and MERS employ an RNA-dependent RNA polymerase to replicate their genomic RNA and transcribe their sub-genomic RNAs18–20. One possibility is that SARS-CoV-2 RNAs could be reverse-transcribed and integrated into the human genome, and transcription of the integrated DNA copies could be responsible for positive PCR tests.Endogenous reverse transcriptase (RT) activity has been observed in human cells, and the products of reverse transcription have been shown to become integrated into the genome21,22. For example, APP transcripts have been shown to be reverse-transcribed by endogenous RT, with resultant APP fragments integrated into the genome of neurons and transcribed22. Human LINE-1 elements (~17% of the human genome), a type of autonomous retrotransposons, are a potential source of endogenous RT, able to retro-transpose themselves and other non-autonomous elements such as Alu21,23.
Results
Expression of viral-cellular chimeric transcripts in infected cultured and in patient-derived cells is consistent with genomic integration of viral sequences
To investigate the possibility of viral integration into virus infected cells we analyzed published RNA-Seq data from SARS-CoV-2-infected cells for evidence of chimeric transcripts, which would be indicative of viral integration into the genome and expression. Examination of these data sets 24–30 (Fig. S1a-b) revealed a substantial number of host-viral chimeric reads (Fig. 1a-c, S1c). These occurred in multiple sample types, including cells and organoids from lung/heart/brain/stomach tissues, as well as BALF cells directly isolated from COVID-19 patients (Fig. 1c). Chimeric read abundance was positively correlated with viral RNA level across the sample types (Fig. 1c). Chimeric reads generally accounted for 0.004% - 0.14% of total SARS-CoV-2 reads across the samples, with a 69.24% maximal number of reads in bronchoalveolar lavage fluid cells derived from severe COVID19 patients and near no chimeric reads from patient blood buffy coat cells (corresponding to almost no total SARS-CoV-2 reads). A majority of chimeric junctions mapped to SARS-CoV-2 nucleocapsid (N) sequence (Fig. 1d-e). This is consistent with the finding that nucleocapsid (N) RNA is the most abundant SARS-CoV-2 sub-genomic RNA31, and thus is most likely to be a target for reverse transcription and integration. These analyses support the hypothesis that SARS-CoV-2 RNA may retro-integrate into the genome of infected cells resulting in the production of chimeric viral-cellular transcripts.
Supplementary Figure 1. Human – SARS-CoV-2 chimeric reads identified from published RNA-Seq data.
a) Published data used to identify human – CoV2 chimeric reads summarizing GEO accession number (data ID), sample type, infection method/type (MOI: Multiplicity Of Infection), RNA-Seq format (single or paired-end with read length), and threshold to call chimeric reads (Min overhang: minimum number of bases mapped to either human or SARS-CoV-2 genome/transcriptome to call a chimeric reads). b) Comparison of SARS-CoV-2 read fraction of total mappable reads in the published RNA-Seq datasets as shown in a). c) One chimeric read example (149 nt) from Calu3 (infected) RNA-Seq with 57 nt mapped to human Chromosome × (green) and 92 nt (magenta) mapped to the SARS-CoV-2 genome.
Figure 1. Human – SARS-CoV-2 chimeric transcripts identified in published data sets of infected cultured and patient-derived cells.
a) Pipeline to identify human-CoV2 chimeric RNA-Seq reads. b) Diagram showing human-CoV2 chimeric reads mapped to potential SARS-CoV-2 integration sites in the human genome from published RNA-Seq data. c) Scatter plot showing human-CoV2 chimeric read number (per million total mappable reads, y-axis) versus SARS-CoV-2 read fraction of total mappable reads (x-axis) in published RNA-Seq datasets (summarized in Supplementary Figure 1a) from different bio-samples with SARS-CoV-2 infection. d-e) Human-CoV2 chimeric read junctions (duplicates removed) mapped to the SARS-CoV-2 genome (d) and distribution among SARS-CoV-2 genes (e, three biological replicates; mean ± s.e.m.). RNA-Seq data is from SARS-CoV-2 infected Calu3 cells (see Supplementary Figure 1a). Chimeric read junction is defined as the “breaking point” of sequences mapped to human or SARS-CoV-2 genome/transcriptome in a given RNA-Seq read.
SARS-CoV-2 RNA can be reverse-transcribed and integrated into the human genome in cells overexpressing a reverse transcriptase
To provide experimental evidence for reverse-transcription and integration of SARS-CoV-2 RNA, we overexpressed human LINE-1 or HIV-1 reverse transcriptase (RT) in HEK293T cells and infected the transduced cells with SARS-CoV-2. The cells were tested 2 days after infection for viral sequences by PCR or fluorescence in situ hybridization (FISH) (Fig. 2a). Considering that the N RNA is the most abundant SARS-CoV-2 sub-genomic RNA31 and is most likely to be retro-integrated (Fig. 1d-e), we chose four N – targeting PCR primer sets that are used in COVID-19 tests (primer source from WHO32, Fig. 2a). PCR amplification of purified cellular DNA showed positive gel-bands in cells with human LINE-1 or HIV-1 RT overexpression (Fig. 2b) but not in non-transfected or non-infected cells. To test whether the DNA copies of N sequences were integrated into the cellular genome, we gel-purified cell genomic DNA (gDNA, >23 kb, Fig. S2a) and qPCR confirmed N sequences in gDNA of cells with expression of all three types of RT (Fig. 2c). Cells with strong expression of LINE-1 driven by a CMV promoter showed ~8-fold higher signals of N sequence detection suggesting a higher copy-number of integrated N sequences than in cells expressing LINE-1 driven by its natural promoter (5’UTR) or HIV-1 RT (Fig. 2c). We were able to clone full-length N DNA from gDNA of cells overexpressing CMV-LINE-1 and confirmed its sequence by Sanger sequencing (Fig. S2b). We did not detect the full-length N sequence from gDNA of cells transfected with 5’UTR-LINE-1 or HIV-1 RT, which may be due to lower expression of RT in these cells (Fig. S2b). We further confirmed that purified SARS-CoV-2 RNA from infected cells can be reverse-transcribed in vitro by lysates of cells expressing either LINE-1 or HIV-1 RT (Fig. S2c-d).
Supplementary Figure 2. SARS-CoV-2 RNA can be reverse-transcribed in vivo and in vitroby other sources of reverse transcriptase.
a) Gel purification of large-fragment genomic DNA (yellow boxes) from SARS-CoV-2 infected HEK293T cells carrying transfected HIV-1 RT, CMV-LINE1 or 5’UTR-LINE1. b) Cloning and Sanger sequencing of DNA copy of full-length SARS-CoV-2 N sequence (CoV2 fl N, yellow box) from gel-purified HEK293T genomic DNA as shown in a). CoV2 fl N: amplification of full-length N sequence (1.26 kb) by primers targeting the two ends of N. HSPA1A: human HSPA1A gene as a control. Note that we can only detect full-length N sequence in gDNA from cells with CMV-LINE-1 expression, corresponding to the high copy-number of integrated N sequences as shown in Figure 2c. c)In vitro reverse transcription of SARS-CoV-2 RNA by adding RNA purified from SARS-Cov-2 infected HEK293T cells to a commercial reverse transcriptase (I, SSIII, with oligo dT and random hexamer primers, positive control), or HEK293T cell lysate (II), or lysates of HEK293T cells expressing HIV-1 reverse transcriptase without (III) or with (IV) heat inactivation. Gel images showing PCR detection of SARS-CoV-2 N sequences from the in vitro reverse transcription products using primer sets (N1 – N4) as shown in Figure 2a. HSPA1A and TUBB: PCR primer sets against human HSPA1A and TUBB genes as controls. d) Same in vitro reverse transcription and PCR detection setup as in c) using lysates of HEK293T cells expressing HIV-1 reverse transcriptase or human LINE-1.
Figure 2. SARS-CoV-2 RNA can be reverse-transcribed and integrated into the host genome in cells with reverse transcriptase expression.
a) Experimental workflow (top), PCR primer sets (shown as amplicons, brown) and single-molecule RNA-FISH probes (red) to detect reverse-transcription and integration of SARS-CoV-2 nucleocapsid (N) sequence (middle, blue). b) PCR detection of SARS-CoV-2 N sequences in cellular DNA purified from mock (left) or SARS-CoV-2 (right) infected HEK293T cells without or with transfection of human LINE-1 (CMV-LINE1 or 5’UTR-LINE1) or HIV-1 RT. HSPA1A: human HSPA1A gene as control; N1 – N4: SARS-CoV-2 N sequences as shown in a). c) qPCR detection and copy-number estimation of SARS-CoV-2 N sequences on gel-purified HEK293T genomic DNA. HSPA1A: human HSPA1A gene as a reference; N2, N3: SARS-CoV-2 N sequences as shown in a). Three biological replicates; mean ± standard error of the mean (s.e.m.). d) Single-molecule RNA-FISH (red) targeting SARS-CoV-2 N sequence using probes shown in a) plus LINE-1 ORF1 protein immuno-staining (green) and merged channels with DAPI (blue) in SARS-CoV-2 infected HEK293T cells with (left) or without (right) transfected LINE-1. Insets: 2.5x enlargement of region in white-box to show nuclear signals of SARS-CoV-2 N sequence (white arrows). Images were single z-slices from 3D optical sections (0.2-μm z-steps). e) Fraction of HEK293T cells infected by SARS-CoV-2 (indicated by cytoplasmic FISH signals) showing nuclear FISH signals of N sequence with (+ LINE-1 O/E, n = 75) or without (-LINE-1 O/E, n = 57) LINE-1 overexpression (indicated by LINE-1 ORF1 protein immuno-staining). Combination of two independent cell samples; Chi-Square Test of Homogeneity.
We conducted single-molecule RNA-FISH (smRNA-FISH) using fluorophore-labeled oligo-nucleotide probes targeting N (Fig. 2a) to confirm that viral N sequences were integrated and detected their transcription in the nucleus. SARS-CoV-2 infected cells showed the expected cytoplasmic FISH signals of N RNA (Fig. S3a). N RNA FISH signals were detected in cell nuclei with cells overexpressing LINE-1 (Fig. 2d, S3b), indicating nascent transcription sites of integrated N sequences. In the same cell population, a significantly higher fraction (~35%) of infected cells overexpressing LINE-1, as indicated by LINE-1 ORF1p immunostaining, showed nuclear N signals than cells not overexpressing LINE-1 (~12%) (Fig. 2e). A significantly higher fraction of infected cells that were transfected with LINE-1 plasmid (~80% transfection efficiency) showed positive nuclear N FISH signals (~30%) as compared to non-transfected cells (13%; Fig. S3c). Infected but not transfected cells also exhibited nuclear N signals, albeit at a lower frequency (~10%; Fig. 2e, S3c), implying integration of SARS-CoV-2 N RNA by cell endogenous RT activity.
Supplementary Figure 3. SARS-CoV-2 N RNA signals detected in cell nuclei by single-molecule RNA-FISH.
a-b) Example images of single-molecule RNA-FISH (red/grey) targeting SARS-CoV-2 N sequence using probes shown in Figure 2a and merged channels with DAPI (blue) in SARS-CoV-2 infected HEK293T cells without (a) or with (b) human LINE-1 transfection. Insets in b): 4x enlargement of regions in white-boxes to show nuclear signals of SARS-CoV-2 N sequence (white arrows). c) Comparison of nuclear N RNA-FISH signals in SARS-CoV-2 infected HEK293T cells without or with human LINE-1 transfection. Left: example images as in a) and b); Right: fraction of HEK293T cells infected by SARS-CoV-2 (indicated by cytoplasmic FISH signals) showing nuclear N RNA-FISH signals in cell populations without (left bar, n = 109) or with (right bar, n = 103) CMV-LINE-1 plasmid transfection (~80% transfection efficiency). Combination of two independent cell samples; Chi-Square Test of Homogeneity. All images shown were single z-slices from 3D optical sections (0.2-μm z-steps).
Human endogenous LINE-1 expression induced by SARS-CoV-2 infection and cytokines correlates with retro-integration
Human LINE-1 elements are autonomous retro-transposons with their encoded reverse transcriptase (ORF2p) and supporting protein (ORF1p) also aiding non-autonomous elements to retro-transpose, such as Alu and other cellular RNAs21. We found that expression of LINE-1 elements was significantly up-regulated in published RNA-Seq data of cells upon infection with SARS-CoV-2 and correlated with chimeric read abundance (Fig. 3a-b, S4a-d, compare Calu3 cells that are efficiently infected versus NHBE cells that are resistant to infection). Although the upregulation in Calu3 was not higher than that in NHBE, multiple LINE-1 elements were upregulated as compared to just one in NHBE (Fig. 3a, S4b, d). Expression analysis using LINE-1 specific primers33,34 showed a ~3-4-fold up-regulation of LINE-1 in Calu3 cells when infected by SARS-CoV-2 (Fig. 3c). Moreover, PCR analysis on Calu3 cellular DNA showed retro-integration of SARS-CoV-2 N sequences after infection (Fig. 3d-e), possibly by the activated LINE-1 reverse transcriptase.
Supplementary Figure 4. LINE-1 induction in human cells correlates with SARS-CoV-2 infection.
a, c) Log2 fold-changes (x-axis) of different types of human repetitive elements (y-axis) with significant (FDR < 0.05) expression changes in SARS-CoV-2 versus mock infected Calu3 (a) or NHBE (c) cells from published RNA-Seq data (GSE147507). b, d) Fold changes (y-axis) of different human LINE-1 families (x-axis) with significant (FDR < 0.05) expression changes in SARS-CoV-2 versus mock infected Calu3 (b) or NHBE (d) cells from published RNA-Seq data (GSE147507, see Supplementary Figure 1a).
Figure 3. LINE-1 expression as an endogenous reverse-transcriptase source in human cells is induced by SARS-CoV-2 infection and cytokine-containing conditioned media treatment.
a) RNA-Seq (GSE147507, see Supplementary Figure 1a) differential expression analysis for all human repetitive elements in SARS-CoV-2 versus mock-infected Calu3 (left) or NHBE (right) cells. Volcano plots showing −Log10 p-values (y-axis) versus Log2 fold-changes (x-axis) for all human repetitive elements with (orange circle) or without (grey circle) significant expression changes (SARS-CoV-2 versus mock-infected); dots: LINE-1 families with (dark blue) or without (light blue) significant expression changes. b) Scatter plot showing human-CoV2 chimeric read number (per million total mappable reads, y-axis) versus SARS-CoV-2 read fraction of total mappable reads (x-axis) in published RNA-Seq (GSE147507, see Supplementary Figure 1a) from infected Calu3 (magenta) or NHBE (red) cells. c) Endogenous LINE-1 expression fold-changes between SARS-CoV-2 versus mock-infected Calu3 cells measured by RT-qPCR with primers probing 5’UTR, ORF1, or 3’UTR regions of LINE-1. Reference genes: GAPDH and TUBB. Three biological replicates; mean ± s.e.m. d) qPCR detection and copy-number estimation of SARS-CoV-2 N sequences in mock (green) or SARS-CoV-2 infected (magenta) Calu3 cellular DNA. HSPA1A: human HSPA1A gene as a reference; N2, N3: SARS-CoV-2 N sequences as shown in Figure 1a. Three biological replicates; mean ± s.e.m; n.d.: not detected. e) Gel purification of large-fragment genomic DNA (yellow box, top) from SARS-CoV-2 infected Calu3 cells and PCR detection of SARS-CoV-2 N sequences in the purified genomic DNA (bottom) with same primer sets as in d). f) Endogenous LINE-1 expression fold-changes in Calu3 cells comparing Myeloid conditioned versus basal media treatment measured by RT-qPCR with primers probing 5’UTR, ORF1, or 3’UTR regions of LINE-1. Reference genes: GAPDH and TUBB. Three biological replicates; mean ± s.e.m. g) LINE-1 ORF1 protein immuno-staining (magenta, same exposure and intensity scaling) plus merged channels with DAPI (blue) in Calu3 cells cultured in basal or myeloid conditioned media. Scale bar: 10 μm. h) Normalized cellular total LINE-1 ORF1p immuno-staining signals of Calu3 cells cultured in basal (n = 84, mean = 1.0, median = 0.9) or myeloid conditioned media (n = 126, mean = 1.7, median = 1.5). Combination of two independent cell samples. Box plots show median (inside line), means (red dot), interquartile range (IQR, box), and upper/lower quartile ± 1.5-times IQRs (whiskers). Welch’s t-test.
Patients infected with SARS-CoV-2 and other corona viruses show evidence of cytokine induction associated with the immune response, and in severe cases experience a cytokine storm35–37, prompting us to investigate whether cytokines alone can induce LINE-1 activation. We treated cells with cytokine-containing conditioned media from Myeloid, Microglia, or CAR-T cell cultures and found a ~2-3-fold upregulation of endogenous LINE-1 expression by PCR analysis (Fig. 3f, S5b). Expressed LINE-1 protein (ORF1p) was also confirmed by immunofluorescence staining (Fig. 3g-h, S5a). In summary, our results show induced LINE-1 expression in cells stressed by viral infection or exposed to cytokines, suggesting a molecular mechanism for SARS-CoV-2 retro-integration in human cells.
Supplementary Figure 5. Cytokine containing media treatment triggers LINE-1 expression in human cells.
a) LINE-1 ORF1 protein immuno-staining (magenta, same exposure and intensity scaling, 1st column: no primary antibody control) plus merged channels with DAPI (blue) in HEK293T cells cultured in basal (1st and 2nd columns) or microglia conditioned media (3rd column) or LPS-treated microglia conditioned media (4th column). Scale bar: 10 μm. b) Endogenous LINE-1 expression fold-changes in Calu3 cells between CAR-T conditioned (diluted with basal media at indicated percentage in volume) versus basal media treatment measured by RT-qPCR with primers probing 5’UTR, ORF1, or 3’UTR regions of LINE-1. Reference genes: GAPDH and TUBB. Three independent cell samples treated with two batches of media; mean ± s.e.m.
Discussion
In this study, we showed evidence that SARS-CoV-2 RNAs can be reverse-transcribed and integrated into the human genome by several sources of reverse transcriptase such as activated human LINE-1 or co-infected retrovirus (HIV). We found LINE-1 expression can be induced upon SARS-CoV-2 infection or cytokine exposure, suggesting a molecular mechanism responsible for SARS-CoV-2 retro-integration in patients. Moreover, our results suggest that the integrated SARS-CoV-2 sequences can be transcribed, as shown by RNA-Seq and smRNA-FISH data, providing a possible explanation for the presence of viral sequences at later times after initial virus exposure and in the absence of detectable infectious virus1–14. The retro-inserted SRAS-CoV-2 sequences are most likely sub-genomic fragments, as the integration junctions are mostly enriched at the N sequence (Fig. 1d-e), excluding the production of infectious virus. Our data may also explain that patients, after recovery from disease symptoms, may become again positive for viral sequences as detected by PCR1,8–14.An important follow-up question is whether these integrated SARS-CoV-2 sequences can express viral antigens. If so, it will be of clinical interest to assess whether viral antigens expressed from integrated virus fragments could trigger an immune response in patients that could affect the course and treatment of the disease. It is possible that the clinical consequences of the integrated viral fragments may depend on their insertion sites in the human genome, and on epigenetic regulation which has been shown in HIV patients38. Careful analysis on SARS-CoV-2 retro-integration sites in patient samples and correlation with disease severity will help to elucidate potential clinical consequences. Furthermore, immune response may vary depending on an individual’s underlying conditions. More generally, our results suggest a novel aspect of infection possibly also for other common disease-causing RNA viruses such as Dengue, Zika or Influenza virus, which could be subject to retro-integration and perhaps affect disease progression.
Human LINE-1 accounts for ~17% of the human genome, ~100 out of 500,000 copies of which are active21,23. LINE-1 – encoded reverse-transcriptase (ORF2p) and supporting protein (ORF1p) are known to retro-transpose not only LINE-1 transcripts (in Cis), but also other RNA species such as Alu (SINE) and cellular mRNA (in Trans, creating processed pseudogenes), with a “target-site – primed reverse transcription” mechanism21. LINE-1 proteins have been shown as nucleic acid chaperones with high RNA binding affinity39, therefore it is perhaps not surprising that they can retro-integrate exogenous viral RNAs. From an evolutionarily perspective, retro-integration of viral RNA by LINE-1 could be an adaptive response by the host to provide sustaining antigen expression possibly enhancing protective immunity. Conversely, retro-integration of viral RNAs could be detrimental and cause a more severe immune response in patients such as a “cytokine storm” or auto-immune reactions.
Our results may also be relevant for current clinical trials of antiviral therapies40. The reliance of PCR tests to assess the effect of treatments on viral replication and viral load may not reflect the efficacy of the treatment to suppress viral replication as the PCR assay may detect viral transcripts from viral sequences stably integrated into the genome rather than infectious virus.
Methods
Cell culture and plasmid transfection
HEK293T cells were obtained from ATCC (CRL-3216) and cultured in DMEM supplemented with 10% heat-inactivated FBS (Hyclone, SH30396.03) and 2mM L-glutamine (MP Biomedicals, IC10180683) following ATCC’s method. Calu3 cells were obtained from ATCC (HTB-55) and cultured in EMEM (ATCC 30-2003) supplemented with 10% heat-inactivated FBS (Hyclone, SH30396.03) following ATCC’s method.Plasmid for HIV-1 reverse transcriptase expression: pCMV-dR8.2 dvpr was a gift from Bob Weinberg (Addgene plasmid # 8455; Addgene: pCMV-dR8.2 dvpr; RRID:Addgene_8455)41. Plasmids for human LINE-1 expression: pBS-L1PA1-CH-mneo (CMV-LINE-1) was a gift from Astrid Roy-Engel (Addgene plasmid # 51288; Addgene: pBS-L1PA1-CH-mneo; RRID:Addgene_51288)42; EF06R (5’UTR-LINE-1) was a gift from Eline Luning Prak (Addgene plasmid # 42940; Addgene: EF06R; RRID:Addgene_42940)43. Transfection was done with Lipofectamine™ 3000 (Invitrogen L3000001) following manufacturer’s protocol.
SARS-CoV-2 infection
SARS-CoV-2 USA-WA1/2020 (Gen Bank: MN985325.1) was obtained from BEI Resources and expanded and tittered on Vero cells. Cells were infected in DMEM +2% FBS for 48 hrs using multiplicity (MOI) of 0.5 for infection of HEK293T cells and an MOI of 2 for Calu3 cells. All sample processing and harvest with infectious virus were done in the BSL3 facility at the Ragon Institute.Nucleic acids extraction, in vitro reverse transcription and PCR/qPCR
DNA extraction was following a published protocol22. For purification of genomic DNA, extracted total cellular DNA was run on 0.4% (w/v) agarose/1x TAE gel for 1.5 hrs with a 3V/cm voltage, with λ DNA-HindIII Digest (NEB N3012S) as size markers. Large fragment bands (>23.13 kb) were cut off, frozen in −80 °C and then crushed by a pipette tip. 3 times of volume (v/w) of high T-E buffer (10 mM Tris – 10 mM EDTA, pH 8.0) was added and then NaCl was added to 200 mM. Gel solution was heated at 70 °C for 15 mins with constant mixing and then extracted with Phenol:Chloroform:Isoamyl Alcohol (25:24:1, v/v) (Life Technologies 15593031) and Chloroform:Isoamyl alcohol 24:1 (Sigma C0549-1PT). DNA was then precipitated by sodium acetate and isopropyl alcohol. For small amount of DNA, glycogen (Life Technologies 10814010) was added as a carrier to aid precipitation.RNA extraction was done with either TRIzol™ LS Reagent (Invitrogen 10296010) or RNeasy Plus Micro Kit (Qiagen 74034) following manufacturers’ protocols. RNA reverse transcription was done with either SuperScript™ III First-Strand Synthesis SuperMix (oligo dT + random hexamer, Invitrogen 18080400) or qScript cDNA SuperMix (QuantaBio 95048-500), following manufacturers’ protocols. In vitro reverse transcription assay for viral RNA by cell lysates was done following a published protocl22.
PCR was done using AccuPrime Taq DNA Polymerase, high fidelity (Life Technologies 12346094). qPCR was done using SYBR™ Green PCR Master Mix (Applied Biosystems 4309155) or PowerUp™ SYBR™ Green Master Mix (Applied Biosystems A25742) in a QuantStudio™ 6 system (Applied Biosystems). See Supplementary Table 1 for primer sequences used in this study. qPCR plots were generated with Prism 8 (Prism).
Supplementary Table 1.
Primer sequences used in this study
Immuno-fluorescence staining and single-molecule RNA-FISH
Cells subject to SARS-CoV-2 infection were grown in μ-Slide 8 Well (#1.5 polymer, Ibidi 80826) and fixed with 4% paraformaldehyde/CMF-PBS at room temperature (RT) for 30 mins. Otherwise, cells were grown on 12 mm round coverslips (#1.5, Warner Instruments 64-0712) and fixed with 1.6% paraformaldehyde/CMF-PBS at room temperature (RT) for 15 mins. Cells were permeabilized with 0.5% (v/v) Triton X-100/PBS, blocked with 4% (w/v) BSA/CMF-PBS at RT for 1 hr, incubated with 1:200 diluted anti-LINE-1 ORF1p mouse monoclonal antibody (clone 4H1, Sigma MABC1152, Lot 3493991), and then with 1:400 diluted Donkey-anti-Mouse-Alexa Fluor 594 second antibody (Invitrogen 21203).Single-molecule RNA-FISH probes (Stellaris®) were ordered from LGC Biosearch Technologies with Quasar® 670 Dye labeling. See Supplementary Table 2 for probe sequences. FISH procedure combining with immuno-fluorescence staining was following previous publications44,45.
Supplementary Table 2.
Single-molecule RNA FISH probe sequences used in this study
Cells in μ-Slide were mounted with Ibidi Mounting Medium With DAPI (Ibidi 50011). Cells on coverslips were mounted with VECTASHIELD® HardSet™ Antifade Mounting Medium with DAPI (Vector Laboratories H-1500-10).
Microscopy and imaging analysis
3D optical sections were acquired with 0.2-μm z-steps using a DeltaVision Elite Imaging System microscope system with a 100 × oil objective (NA 1.4) and a pco.edge 5.5 camera and DeltaVision SoftWoRx software (GE Healthcare). Image deconvolution was done using SoftWoRx. All figure panel images were prepared using FIJI software (ImageJ, NIH) and Adobe Illustrator 2020 (Adobe), showing deconvolved single z-slices.To measure the LINE-1 ORF1p immuno-staining signal intensity, we projected cell optical sections (sum, 42 slices) with the “z projection” function in FIJI. We measured the sum of intensity of the entire cell area in the z-projected image as the signal intensity, subtracted the background intensity outside of cells and then divided by the mean of the “Basal media treatment” group to have the normalized signal intensity, as previously described44,45. All images from the same experiment were using the same exposure time and transmitted exciting light. All intensity measurements were done with non-deconvolved raw images. Box plot was done in R (version 4.0.3)46.
RNA-Seq data analysis
RNA-Seq data were downloaded from GEO with the accession numbers GSE14750724, GSE15327725, GSE15675426, GSE15785227, GSE15368428, GSE14592629, GSE15499830 (summarized in Supplementary Figure 1a).To identify human – SARS-CoV-2 chimeric reads, raw sequencing reads were aligned to concatenated human and SARS-CoV-2 genomes plus transcriptomes by STAR (version 2.7.1a)47. Human genome version hg38 with no alternative chromosomes and gene annotation version GRCh38.97 were used. SARS-CoV-2 genome version NC_045512.2 and gene annotation (Index of /goldenPath/wuhCor1/bigZips/genes) were used. The following STAR parameters31 were used to call chimeric reads unless otherwise specified (Supplementary Figure 1a): --chimOutType Junctions SeparateSAMold WithinBAM HardClip \ --chimScoreJunctionNonGTAG 0 \ --alignSJstitchMismatchNmax −1 −1 −1 −1 \ --chimSegmentMin 50 \ --chimJunctionOverhangMin 50.
To analyze human LINE-1 expression in RNA-Seq data, a published method, RepEnrich248, was used to map RNA-Seq reads to human repeat annotations, using human repeat masker (hg38). Differential expression was analyzed using EdgeR package (version 3.30.3)49,50 in R (version 4.0.3)46.
Conditioned media production and treatment
As previously described51, myeloid precursors were derived from human pluripotent stem cells. Briefly, human embryonic stem cells were cultured in StemFlex (ThermoFisher) feeder-free medium on Matrigel™-(Corning) coated tissue culture polystyrene. 24 hrs before single-cell harvesting via TrypLE Express (ThermoFisher), cells were treated with 10 μM ROCK Inhibitor (Y-27632) (Stem Cell Technologies) in Essential 8 (E8) medium (ThermoFisher). After harvesting, cells were centrifuged at 300 g for 3 mins in non-adherent U-bottom 96-well plates (Corning) at 10,000 cells per 150 μL/well of embryoid body (EB) medium consisting of 10 μM ROCK Inhibitor, 50 ng/mL BMP-4 (Peprotech), 20 ng/mL SCF (Peprotech), 50 ng/mL VEGF (Peprotech), and 100 U/mL Penn/Strep (ThermoFisher) in E8 base medium. EBs were cultured in the 96-well plates for 4 days with 150 μL/well of EB medium added at day 2. After 4 days, 16 EBs/well were plated in a 6-well tissue culture polystyrene plated coated with Matrigel™ in hematopoietic myeloid medium (HIM) consisting of 2mM GlutaMax (ThermoFisher), 55 μM beta-mercaptoethanol, 100 ng/mL M-CSF (Peprotech), and 25 ng/mL IL-3 (Peprotech) in X-VIVO 15 base medium (Lonza). HIM media was changed every 3-4 days for 2-3 weeks until floating CD14-positive myeloid precursors emerged. Myeloid conditioned media consisted of floating myeloid cells cultured in HIM media for 7 days at a concentration of 0.5 × 106 – 1 × 106 cells/mL. Cells in conditioned media were removed by centrifugation and filtration through 0.2 μM filters. Calu3 cells were cultured in the myeloid conditioned media or HIM media (basal) for two days with daily media change before harvest or fixation.Microglia were differentiated from human induced pluripotent stem cells (hiPSCs) via embryoid bodies and primitive macrophage precursors (PMPs)51. In brief, hiPSCs (cultured feeder-free on matrigel in StemFlexTM (Gibco)) were dissociated with TrypLE Express (Gibco), and 10,000 cells were plated per well in 96-well ultra-low attachment plates (Corning) in 100 μL embryoid body medium (10 μM ROCK inhibitor, 50 ng/mL BMP-4, 20 ng/mL SCF, and 50 ng/mL VEGF-121 in StemFlex), before centrifugation at 300 × g for 3 mins at 4 °C. Embryoid bodies were cultured for 4 days, with adding 100 μL embryoid body medium after 2 days. 12 to 16 embryoid bodies were plated per well of tissue culture-treated 6-well plates and cultured in 3 mL hematopoetic medium (2 mM GlutaMax, 100 U/mL penicillin, 100 μg/mL streptomycin, 55 μM β-mercaptoethanol, 100 ng/mL M-CSF, 25 ng/mL IL-3, 100 U/mL penicillin, 100 μg/mL streptomycin in X-VIVO 15 (Lonza, BW04418Q). From this point on, 2 mL medium was exchanged every 4–7 days. PMPs were harvested from suspension during medium exchange and plated in microglia differentiation media over 7-14 days to produce microglia like cell monocultures (Neurobasal (Life Technologies 21103049) supplemented with Gem21 NeuroPlex without Vitamin A (GeminiBio, 400-161), 2mM GlutaMAX (Gibco), 100 ng/mL IL-34, and 10 ng/mL GM-CSF, 100 U/mL penicillin, 100 μg/mL streptomycin). For microglia stimulation, microglia differentiation media was exchanged with HEK293T media (DMEM + 10% heat-inactivated FBS + final 2mM L-Glutamine) and supplemented with 100 hg/ml lipopolysaccharide (LPS, Sigma Aldrich L4391-1MG) or PBS. After 24 hrs, the microglia conditioned media was collected, centrifugated (1000 rpm 10min) and the supernatant was directly applied to HEK293T cells. HEK293T cells received microglia conditioned media or basal HEK293T media on three constitutive days before fixation.
Human anti-CD19 CAR-T cells were generated by transduction of primary T cells purified from human peripheral blood mononuclear cells (PBMC) with CD19-CAR expressing retrovirus52. Anti-CD19 CAR-T cells were co-cultured with CD19-expressing beta-like cells52 or WIBR3 cells with a luciferase-2A-CD19 expressing cassette integrated at the AAVS1 locus in RPMI1640 medium with 10% human serum AB. Cells in the conditioned medium were removed by filtration through 0.45 μM filters. RPMI1640 medium with 10% human serum AB was used as basal media control. Calu3 cells were cultured in the CAR-T conditioned media with indicated dilutions or in the basal media for two days before harvest.
Data Availability
The datasets generated during and/or analysed during the current study are available from the corresponding author on reasonable request.Author contributions
Project design by R.J. and R.A.Y, execution of experiments and data analysis by L.Z., A.R, R.J and R.A.Y; E.W., A.K., and H.M. generated cells and reagents; Manuscript preparation by L.Z. and R.J. with input from all authors.Competing interests
R.J. is an advisor/co-founder of Fate Therapeutics, Fulcrum Therapeutics, Omega Therapeutics, and Dewpoint Therapeutics. R.A.Y. is a founder and shareholder of Syros Pharmaceuticals, Camp4 Therapeutics, Omega Therapeutics, and Dewpoint Therapeutics. All other authors declare no competing interests.Acknowledgements
We thank members in the laboratories of Rudolf Jaenisch and Richard Young and other colleagues from Whitehead Institute and MIT for helpful discussions and resources. We thank Wendy Salmon from the Whitehead W.M. Keck Microscopy Facility and M. Inmaculada Barrasa from the Whitehead Bioinformatics and Research Computing for technical advice. This work was supported by grants from the NIH to RJ (1U19AI131135-01, 5R01MH104610-21) and by a generous gift from Dewpoint Therapeutics and from Jim Stone. ASK would like to acknowledge funding from the NIH (Grant: T32 EB016652). Finally, we thank Nathans Island for inspiration.References
- KDCA.
- Bullard, J. et al. Predicting Infectious Severe Acute Respiratory Syndrome Coronavirus 2 From Diagnostic Samples. Clinical Infectious Diseases, doi:10.1093/cid/ciaa638 (2020).
CrossRefPubMedGoogle Scholar - He, X. et al. Temporal dynamics in viral shedding and transmissibility of COVID-19. Nature Medicine 26, 672–675, doi:10.1038/s41591-020-0869-5 (2020).
CrossRefPubMedGoogle Scholar - Li, N., Wang, X. & Lv, T. Prolonged SARS-CoV-2 RNA shedding: Not a rare phenomenon. J Med Virol, doi:10.1002/jmv.25952 (2020).
CrossRefGoogle Scholar - Mina, M. J., Parker, R. & Larremore, D. B. Rethinking Covid-19 Test Sensitivity — A Strategy for Containment. New England Journal of Medicine 383, e120, doi:10.1056/NEJMp2025631 (2020).
CrossRefPubMedGoogle Scholar - Sethuraman, N., Jeremiah, S. S. & Ryo, A. Interpreting Diagnostic Tests for SARS-CoV-2. JAMA 323, 2249–2251, doi:10.1001/jama.2020.8259 (2020).
CrossRefGoogle Scholar - Yang, J.-R. et al. Persistent viral RNA positivity during the recovery period of a patient with SARS-CoV-2 infection. Journal of Medical Virology 92, 1681–1683, doi:Error - Cookies Turned Off (2020).
Google Scholar - An, J. et al. Clinical characteristics of recovered COVID-19 patients with re-detectable positive RNA test. Ann Transl Med 8, 1084, doi:10.21037/atm-20-5602 (2020).
CrossRefGoogle Scholar - Chen, D. et al. Recurrence of positive SARS-CoV-2 RNA in COVID-19: A case report. Int J Infect Dis 93, 297–299, doi:10.1016/j.ijid.2020.03.003 (2020).
CrossRefPubMedGoogle Scholar - Lan, L. et al. Positive RT-PCR Test Results in Patients Recovered From COVID-19. JAMA 323, 1502–1503, doi:10.1001/jama.2020.2783 (2020).
CrossRefPubMedGoogle Scholar - Loconsole, D. et al. Recurrence of COVID-19 after recovery: a case report from Italy. Infection, doi:10.1007/s15010-020-01444-1 (2020).
CrossRefGoogle Scholar - Lu, J. et al. Clinical, immunological and virological characterization of COVID-19 patients that test re-positive for SARS-CoV-2 by RT-PCR. EBioMedicine 59, 102960, doi:10.1016/j.ebiom.2020.102960 (2020).
CrossRefGoogle Scholar - Luo, S., Guo, Y., Zhang, X. & Xu, H. A follow-up study of recovered patients with COVID-19 in Wuhan, China. Int J Infect Dis 99, 408–409, doi:10.1016/j.ijid.2020.05.119 (2020).
CrossRefGoogle Scholar - Ye, G. et al. Clinical characteristics of severe acute respiratory syndrome coronavirus 2 reactivation. J Infect 80, e14–e17, doi:10.1016/j.jinf.2020.03.001 (2020).
CrossRefPubMedGoogle Scholar - To, K. K. et al. COVID-19 re-infection by a phylogenetically distinct SARS-coronavirus-2 strain confirmed by whole genome sequencing. Clin Infect Dis, doi:10.1093/cid/ciaa1275 (2020).
CrossRefPubMedGoogle Scholar - Huang, J. et al. Recurrence of SARS-CoV-2 PCR positivity in COVID-19 patients: a single center experience and potential implications. medRxiv, 2020.2005.2006.20089573, doi:10.1101/2020.05.06.20089573 (2020).
Abstract/FREE Full TextGoogle Scholar - Yuan, B. et al. Recurrence of positive SARS-CoV-2 viral RNA in recovered COVID-19 patients during medical isolation observation. Sci Rep 10, 11887, doi:10.1038/s41598-020-68782-w (2020).
CrossRefGoogle Scholar - Alanagreh, L., Alzoughool, F. & Atoum, M. The Human Coronavirus Disease COVID-19: Its Origin, Characteristics, and Insights into Potential Drugs and Its Mechanisms. Pathogens 9, doi:10.3390/pathogens9050331 (2020).
CrossRefGoogle Scholar - Fehr, A. R. & Perlman, S. Coronaviruses: an overview of their replication and pathogenesis. Methods Mol Biol 1282, 1–23, doi:10.1007/978-1-4939-2438-7_1 (2015).
CrossRefPubMedGoogle Scholar - de Wit, E., van Doremalen, N., Falzarano, D. & Munster, V. J. SARS and MERS: recent insights into emerging coronaviruses. Nat Rev Microbiol 14, 523–534, doi:10.1038/nrmicro.2016.81 (2016).
CrossRefPubMedGoogle Scholar - Kazazian, H. H., Jr.. & Moran, J. V. Mobile DNA in Health and Disease. N Engl J Med 377, 361–370, doi:10.1056/NEJMra1510092 (2017).
CrossRefPubMedGoogle Scholar - Lee, M. H. et al. Somatic APP gene recombination in Alzheimer’s disease and normal neurons. Nature 563, 639–645, doi:10.1038/s41586-018-0718-6 (2018).
CrossRefPubMedGoogle Scholar - Brouha, B. et al. Hot L1s account for the bulk of retrotransposition in the human population. Proceedings of the National Academy of Sciences of the United States of America 100, 5280–5285, doi:10.1073/pnas.0831042100 (2003).
Abstract/FREE Full TextGoogle Scholar - Blanco-Melo, D. et al. Imbalanced Host Response to SARS-CoV-2 Drives Development of COVID-19. Cell 181, 1036–1045 e1039, doi:10.1016/j.cell.2020.04.026 (2020).
CrossRefPubMedGoogle Scholar - Huang, J. et al. SARS-CoV-2 Infection of Pluripotent Stem Cell-Derived Human Lung Alveolar Type 2 Cells Elicits a Rapid Epithelial-Intrinsic Inflammatory Response. Cell stem cell, doi:10.1016/j.stem.2020.09.013 (2020).
CrossRefGoogle Scholar - Perez-Bermejo, J. A. et al. SARS-CoV-2 infection of human iPSC-derived cardiac cells predicts novel cytopathic features in hearts of COVID-19 patients. bioRxiv, doi:10.1101/2020.08.25.265561 (2020).
Abstract/FREE Full TextGoogle Scholar - Jacob, F. et al. Human Pluripotent Stem Cell-Derived Neural Cells and Brain Organoids Reveal SARS-CoV-2 Neurotropism Predominates in Choroid Plexus Epithelium. Cell stem cell, doi:10.1016/j.stem.2020.09.016 (2020).
CrossRefGoogle Scholar - Giobbe, G. G. et al. SARS-CoV-2 infection and replication in human fetal and pediatric gastric organoids. bioRxiv,2020.2006.2024.167049, doi:10.1101/2020.06.24.167049 (2020).
Abstract/FREE Full TextGoogle Scholar - Liao, M. et al. Single-cell landscape of bronchoalveolar immune cells in patients with COVID-19. Nat Med 26, 842–844, doi:10.1038/s41591-020-0901-9 (2020).
CrossRefPubMedGoogle Scholar - Gill, S. E. et al. Transcriptional Profiling of Leukocytes in Critically Ill COVID19 Patients: Implications for Interferon Response and Coagulation. PREPRINT (Version 2) available at Research Square, doi:10.21203/rs.3.rs-63632/v2 (2020).
CrossRefGoogle Scholar - Kim, D. et al. The Architecture of SARS-CoV-2 Transcriptome. Cell 181, 914–921 e910, doi:10.1016/j.cell.2020.04.011 (2020).
CrossRefPubMedGoogle Scholar - https://www.who.int/docs/default-source/coronaviruse/whoinhouseassays.pdf?sfvrsn=de3a76aa_2.
- Jones, R. B. et al. LINE-1 retrotransposable element DNA accumulates in HIV-1-infected cells. J Virol 87, 13307–13320, doi:10.1128/JVI.02257-13 (2013).
Abstract/FREE Full TextGoogle Scholar - Tiwari, B. et al. p53 directly represses human LINE1 transposons. Genes & development, doi:10.1101/gad.343186.120 (2020).
Abstract/FREE Full TextGoogle Scholar - Channappanavar, R. & Perlman, S. Pathogenic human coronavirus infections: causes and consequences of cytokine storm and immunopathology. Semin Immunopathol 39, 529–539, doi:10.1007/s00281-017-0629-x (2017).
CrossRefPubMedGoogle Scholar - Mehta, P. et al. COVID-19: consider cytokine storm syndromes and immunosuppression. The Lancet 395, 1033–1034, doi:10.1016/s0140-6736(20)30628-0 (2020).
CrossRefGoogle Scholar - Costela-Ruiz, V. J., Illescas-Montes, R., Puerta-Puerta, J. M., Ruiz, C. & Melguizo-Rodriguez, L. SARS-CoV-2 infection: The role of cytokines in COVID-19 disease. Cytokine Growth Factor Rev 54, 62–75, doi:10.1016/j.cytogfr.2020.06.001 (2020).
CrossRefGoogle Scholar - Jiang, C. et al. Distinct viral reservoirs in individuals with spontaneous control of HIV-1. Nature 585, 261–267, doi:10.1038/s41586-020-2651-8 (2020).
CrossRefPubMedGoogle Scholar - Naufer, M. N., Furano, A. V. & Williams, M. C. Protein-nucleic acid interactions of LINE-1 ORF1p. Semin Cell Dev Biol 86, 140–149, doi:10.1016/j.semcdb.2018.03.019 (2019).
CrossRefGoogle Scholar - Baum, A. et al. REGN-COV2 antibodies prevent and treat SARS-CoV-2 infection in rhesus macaques and hamsters. Science 370, 1110, doi:10.1126/science.abe2402 (2020).
Abstract/FREE Full TextGoogle Scholar - Stewart, S. A. et al. Lentivirus-delivered stable gene silencing by RNAi in primary cells. RNA 9, 493–501, doi:10.1261/rna.2192803 (2003).
Abstract/FREE Full TextGoogle Scholar - Wagstaff, B. J., Barnerssoi, M. & Roy-Engel, A. M. Evolutionary conservation of the functional modularity of primate and murine LINE-1 elements. PloS one 6, e19672, doi:10.1371/journal.pone.0019672 (2011).
CrossRefPubMedGoogle Scholar - Farkash, E. A., Kao, G. D., Horman, S. R. & Prak, E. T. Gamma radiation increases endonuclease-dependent L1 retrotransposition in a cultured cell assay. Nucleic acids research 34, 1196–1204, doi:10.1093/nar/gkj522 (2006).
CrossRefPubMedWeb of ScienceGoogle Scholar - Zhang, L. et al. TSA-Seq reveals a largely “hardwired” genome organization relative to nuclear speckles with small position changes tightly correlated with gene expression changes. bioRxiv, 824433, doi:10.1101/824433 (2020).
Abstract/FREE Full TextGoogle Scholar - Kim, J., Venkata, N. C., Hernandez Gonzalez, G. A., Khanna, N. & Belmont, A. S. Gene expression amplification by nuclear speckle association. The Journal of cell biology 219, doi:10.1083/jcb.201904046 (2020).
Abstract/FREE Full TextGoogle Scholar - R Core Team (2020). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. URL R: The R Project for Statistical Computing.
Google Scholar - Dobin, A. et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics 29, 15–21, doi:10.1093/bioinformatics/bts635 (2013).
CrossRefPubMedWeb of ScienceGoogle Scholar - Criscione, S. W., Zhang, Y., Thompson, W., Sedivy, J. M. & Neretti, N. Transcriptional landscape of repetitive elements in normal and cancer human cells. BMC Genomics 15, 583, doi:10.1186/1471-2164-15-583 (2014).
CrossRefPubMedGoogle Scholar - Robinson, M. D., McCarthy, D. J. & Smyth, G. K. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26, 139–140, doi:10.1093/bioinformatics/btp616 (2010).
CrossRefPubMedWeb of ScienceGoogle Scholar - McCarthy, D. J., Chen, Y. & Smyth, G. K. Differential expression analysis of multifactor RNA-Seq experiments with respect to biological variation. Nucleic acids research 40, 4288–4297, doi:10.1093/nar/gks042 (2012).
CrossRefPubMedWeb of ScienceGoogle Scholar - Brownjohn, P. W. et al. Functional Studies of Missense TREM2 Mutations in Human Stem Cell-Derived Microglia. Stem cell reports 10, 1294–1307, doi:10.1016/j.stemcr.2018.03.003 (2018).
CrossRefPubMedGoogle Scholar - Ma, H., Jeppesen, J. F. & Jaenisch, R. Human T Cells Expressing a CD19 CAR-T Receptor Provide Insights into Mechanisms of Human CD19-Positive beta Cell Destruction. Cell Rep Med 1, 100097, doi:10.1016/j.xcrm.2020.100097 (2020).
CrossRefGoogle Scholar
