Exogenous Coronavirus Interacts With Endogenous Retrotransposon in Human Cells

Beyond acute infection: molecular mechanisms underpinning cardiovascular complications in long COVID

SARS-CoV-2, responsible for the global COVID-19 pandemic, has manifested significant cardiovascular implications for the infected population. These cardiovascular repercussions not only linger beyond the initial phase of illness but have also been observed in individuals who remain asymptomatic. This extended and pervasive impact is often called the post-acute COVID-19 syndrome (PACS) or "Long COVID". With the number of confirmed global cases approaching an alarming 756 million, the multifaceted challenges of Long COVID are undeniable. These challenges span from individual health complications to considerable burdens on worldwide healthcare systems. Our review comprehensively examines the complications of the persistent cardiovascular complications associated with COVID-19. Furthermore, we shed light on emerging therapeutic strategies that promise to manage and possibly mitigate these complications. We also introduce and discuss the profound concerns regarding the potential transgenerational repercussions of SARS-CoV-2, emphasizing the need for a proactive and informed approach to future research and clinical practice.

The controversy of SARS-CoV-2 integration into the human genome

Bat borne disease have attracted many researchers for years. The ability of the bat to host several exogenous viruses has been a focal point in research lately. The latest pandemic shifted the focus of scholars towards understanding the difference in response to viral infection between humans and bats. In a way to understand the basis of the interaction and behaviour between SARS-CoV-2 and the environment, a conflict between different researchers across the globe arose. This conflict asked many questions about the truth of virus-host integration, whether an interaction between RNA viruses and human genomes has ever been reported, the possible route and mechanism that could lead to genomic integration of viral sequences and the methods used to detect integration. This article highlights those questions and will discuss the diverse opinions of the controversy and provide examples on reported integration mechanisms and possible detection techniques.

Computational analysis of the interactions between Ebselen and derivatives with the active site of the main protease from SARS-CoV-2

The main protease (Mpro) of the novel coronavirus SARS-CoV-2 is a key target for developing antiviral drugs. Ebselen (EbSe) is a selenium-containing compound that has been shown to inhibit Mpro in vitro by forming a covalent bond with the cysteine (Cys) residue in the active site of the enzyme. However, EbSe can also bind to other proteins, like albumin, and low molecular weight compounds that have free thiol groups, such as Cys and glutathione (GSH), which may affect its availability and activity. In this study, we analyzed the Mpro interaction with EbSe, its analogues, and its metabolites with Cys, GSH, and albumin by molecular docking. We also simulated the electronic structure of the generated molecules by density functional theory (DFT) and explored the stability of EbSe and one of its best derivatives, EbSe-2,5-MeClPh, in the catalytic pocket of Mpro through covalent docking and molecular dynamics. Our results show that EbSe and its analogues bound to GSH/albumin have larger distance between the selenium atom of the ligands and the sulfur atom of Cys145 of Mpro than the other compounds. This suggests that EbSe and its GSH/albumin-analogues may have less affinity for the active site of Mpro. EbSe-2,5-MeClPh was found one of the best molecules, and in molecular dynamics simulations, it showed to undergo more conformational changes in the active site of Mpro, in relation to EbSe, which remained stable in the catalytic pocket. Moreover, this study also reveals that all compounds have the potential to interact closely with the active site of Mpro, providing us with a concept of which derivatives may be promising for in vitro analysis in the future. We propose that these compounds are potential covalent inhibitors of Mpro and that organoselenium compounds are molecules that should be studied for their antiviral properties.

Single-cell analysis reveals transcriptomic and epigenomic impacts on the maternal-fetal interface following SARS-CoV-2 infection

During pregnancy the maternal-fetal interface plays vital roles in fetal development. Its disruption is frequently found in pregnancy complications. Recent studies show increased incidences of adverse pregnancy outcomes in patients with COVID-19; however, the mechanism remains unclear. Here we analysed the molecular impacts of SARS-CoV-2 infection on the maternal-fetal interface. Generating bulk and single-nucleus transcriptomic and epigenomic profiles from patients with COVID-19 and control samples, we discovered aberrant immune activation and angiogenesis patterns in distinct cells from patients. Surprisingly, retrotransposons were also dysregulated in specific cell types. Notably, reduced enhancer activities of LTR8B elements were functionally linked to the downregulation of pregnancy-specific glycoprotein genes in syncytiotrophoblasts. Our findings revealed that SARS-CoV-2 infection induced substantial changes to the epigenome and transcriptome at the maternal-fetal interface, which may be associated with pregnancy complications.

Contribution of Retrotransposons to the Pathogenesis of Type 1 Diabetes and Challenges in Analysis Methods

Type 1 diabetes (T1D) is one of the most common chronic diseases of the endocrine system, associated with several life-threatening comorbidities. While the etiopathogenesis of T1D remains elusive, a combination of genetic susceptibility and environmental factors, such as microbial infections, are thought to be involved in the development of the disease. The prime model for studying the genetic component of T1D predisposition encompasses polymorphisms within the HLA (human leukocyte antigen) region responsible for the specificity of antigen presentation to lymphocytes. Apart from polymorphisms, genomic reorganization caused by repeat elements and endogenous viral elements (EVEs) might be involved in T1D predisposition. Such elements are human endogenous retroviruses (HERVs) and non-long terminal repeat (non-LTR) retrotransposons, including long and short interspersed nuclear elements (LINEs and SINEs). In line with their parasitic origin and selfish behaviour, retrotransposon-imposed gene regulation is a major source of genetic variation and instability in the human genome, and may represent the missing link between genetic susceptibility and environmental factors long thought to contribute to T1D onset. Autoreactive immune cell subtypes with differentially expressed retrotransposons can be identified with single-cell transcriptomics, and personalized assembled genomes can be constructed, which can then serve as a reference for predicting retrotransposon integration/restriction sites. Here we review what is known to date about retrotransposons, we discuss the involvement of viruses and retrotransposons in T1D predisposition, and finally we consider challenges in retrotransposons analysis methods.

Transcriptome and DNA methylome analysis of peripheral blood samples reveals incomplete restoration and transposable element activation after 3-months recovery of COVID-19

Comprehensive analyses showed that SARS-CoV-2 infection caused COVID-19 and induced strong immune responses and sometimes severe illnesses. However, cellular features of recovered patients and long-term health consequences remain largely unexplored. In this study, we collected peripheral blood samples from nine recovered COVID-19 patients (median age of 36 years old) from Hubei province, China, 3 months after discharge as well as 5 age- and gender-matched healthy controls; and carried out RNA-seq and whole-genome bisulfite sequencing to identify hallmarks of recovered COVID-19 patients. Our analyses showed significant changes both in transcript abundance and DNA methylation of genes and transposable elements (TEs) in recovered COVID-19 patients. We identified 425 upregulated genes, 214 downregulated genes, and 18,516 differentially methylated regions (DMRs) in total. Aberrantly expressed genes and DMRs were found to be associated with immune responses and other related biological processes, implicating prolonged overreaction of the immune system in response to SARS-CoV-2 infection. Notably, a significant amount of TEs was aberrantly activated and their activation was positively correlated with COVID-19 severity. Moreover, differentially methylated TEs may regulate adjacent gene expression as regulatory elements. Those identified transcriptomic and epigenomic signatures define and drive the features of recovered COVID-19 patients, helping determine the risks of long COVID-19, and guiding clinical intervention.

SARS-CoV-2-encoded inhibitors of human LINE-1 retrotransposition

The ongoing pandemic of severe acute respiratory coronavirus 2 (SARS-CoV-2) is causing a devastating impact on public health worldwide. However, details concerning the profound impact of SARS-CoV-2 on host cells remain elusive. Here, we investigated the effects of SARS-CoV-2-encoded viral proteins on the intracellular activity of long interspersed element 1 (L1) retrotransposons using well-established reporter systems. Several nonstructural or accessory proteins (Nsps) of SARS-CoV-2 (i.e., Nsp1, Nsp3, Nsp5, and Nsp14) significantly suppress human L1 mobility, and these viral L1 inhibitors generate a complex network that modulates L1 transposition. Specifically, Nsp1 and Nsp14 inhibit the intracellular accumulation of L1 open reading frame proteins (ORF1p), whereas Nsp3, Nsp5, and Nsp14 repress the reverse transcriptase activity of L1 ORF2p. Given recent findings concerning the roles of L1 in antiviral immune activation and host genome instability, the anti-L1 activities mediated by SARS-CoV-2-encoded inhibitors suggest that SARS-CoV-2 employs different strategies to optimize the host genetic environment.

SARS-CoV-2 and the Nucleus

The ongoing COVID-19 pandemic is caused by an RNA virus, SARS-CoV-2. The genome of SARS-CoV-2 lacks a nuclear phase in its life cycle and is replicated in the cytoplasm. However, interfering with nuclear trafficking using pharmacological inhibitors greatly reduces virus infection and virus replication of other coronaviruses is blocked in enucleated cells, suggesting a critical role of the nucleus in virus infection. Here, we summarize the alternations of nuclear pathways caused by SARS-CoV-2, including nuclear translocation pathways, innate immune responses, mRNA metabolism, epigenetic mechanisms, DNA damage response, cytoskeleton regulation, and nuclear rupture. We consider how these alternations contribute to virus replication and discuss therapeutic treatments that target these pathways, focusing on small molecule drugs that are being used in clinical studies.

Melatonin: Regulation of Viral Phase Separation and Epitranscriptomics in Post-Acute Sequelae of COVID-19

The relentless, protracted evolution of the SARS-CoV-2 virus imposes tremendous pressure on herd immunity and demands versatile adaptations by the human host genome to counter transcriptomic and epitranscriptomic alterations associated with a wide range of short- and long-term manifestations during acute infection and post-acute recovery, respectively. To promote viral replication during active infection and viral persistence, the SARS-CoV-2 envelope protein regulates host cell microenvironment including pH and ion concentrations to maintain a high oxidative environment that supports template switching, causing extensive mitochondrial damage and activation of pro-inflammatory cytokine signaling cascades. Oxidative stress and mitochondrial distress induce dynamic changes to both the host and viral RNA m6A methylome, and can trigger the derepression of long interspersed nuclear element 1 (LINE1), resulting in global hypomethylation, epigenetic changes, and genomic instability. The timely application of melatonin during early infection enhances host innate antiviral immune responses by preventing the formation of "viral factories" by nucleocapsid liquid-liquid phase separation that effectively blockades viral genome transcription and packaging, the disassembly of stress granules, and the sequestration of DEAD-box RNA helicases, including DDX3X, vital to immune signaling. Melatonin prevents membrane depolarization and protects cristae morphology to suppress glycolysis via antioxidant-dependent and -independent mechanisms. By restraining the derepression of LINE1 via multifaceted strategies, and maintaining the balance in m6A RNA modifications, melatonin could be the quintessential ancient molecule that significantly influences the outcome of the constant struggle between virus and host to gain transcriptomic and epitranscriptomic dominance over the host genome during acute infection and PASC.

Multi-Omics Integration Reveals Only Minor Long-Term Molecular and Functional Sequelae in Immune Cells of Individuals Recovered From COVID-19

The majority of COVID-19 patients experience mild to moderate disease course and recover within a few weeks. An increasing number of studies characterized the long-term changes in the specific anti-SARS-CoV-2 immune responses, but how COVID-19 shapes the innate and heterologous adaptive immune system after recovery is less well known. To comprehensively investigate the post-SARS-CoV-2 infection sequelae on the immune system, we performed a multi-omics study by integrating single-cell RNA-sequencing, single-cell ATAC-sequencing, genome-wide DNA methylation profiling, and functional validation experiments in 14 convalescent COVID-19 and 15 healthy individuals. We showed that immune responses generally recover without major sequelae after COVID-19. However, subtle differences persist at the transcriptomic level in monocytes, with downregulation of the interferon pathway, while DNA methylation also displays minor changes in convalescent COVID-19 individuals. However, these differences did not affect the cytokine production capacity of PBMCs upon different bacterial, viral, and fungal stimuli, although baseline release of IL-1Ra and IFN-γ was higher in convalescent individuals. In conclusion, we propose that despite minor differences in epigenetic and transcriptional programs, the immune system of convalescent COVID-19 patients largely recovers to the homeostatic level of healthy individuals.

Human endogenous retrovirus K in the respiratory tract is associated with COVID-19 physiopathology

BACKGROUND

Critically ill 2019 coronavirus disease (COVID-19) patients under invasive mechanical ventilation (IMV) are 10 to 40 times more likely to die than the general population. Although progression from mild to severe COVID-19 has been associated with hypoxia, uncontrolled inflammation, and coagulopathy, the mechanisms involved in the progression to severity are poorly understood.

METHODS

The virome of tracheal aspirates (TA) from 25 COVID-19 patients under IMV was assessed through unbiased RNA sequencing (RNA-seq), and correlation analyses were conducted using available clinical data. Unbiased sequences from nasopharyngeal swabs (NS) from mild cases and TA from non-COVID patients were included in our study for further comparisons.

RESULTS

We found higher levels and differential expression of human endogenous retrovirus K (HERV-K) genes in TA from critically ill and deceased patients when comparing nasopharyngeal swabs from mild cases to TA from non-COVID patients. In critically ill patients, higher HERV-K levels were associated with early mortality (within 14 days of diagnosis) in the intensive care unit. Increased HERV-K expression in deceased patients was associated with IL-17-related inflammation, monocyte activation, and an increased consumption of clotting/fibrinolysis factors. Moreover, increased HERV-K expression was detected in human primary monocytes from healthy donors after experimental SARS-CoV-2 infection in vitro.

CONCLUSION

Our data implicate the levels of HERV-K transcripts in the physiopathology of COVID-19 in the respiratory tract of patients under invasive mechanical ventilation. Video abstract.

mRNA Vaccines: Why Is the Biology of Retroposition Ignored?

The major advantage of mRNA vaccines over more conventional approaches is their potential for rapid development and large-scale deployment in pandemic situations. In the current COVID-19 crisis, two mRNA COVID-19 vaccines have been conditionally approved and broadly applied, while others are still in clinical trials. However, there is no previous experience with the use of mRNA vaccines on a large scale in the general population. This warrants a careful evaluation of mRNA vaccine safety properties by considering all available knowledge about mRNA molecular biology and evolution. Here, I discuss the pervasive claim that mRNA-based vaccines cannot alter genomes. Surprisingly, this notion is widely stated in the mRNA vaccine literature but never supported by referencing any primary scientific papers that would specifically address this question. This discrepancy becomes even more puzzling if one considers previous work on the molecular and evolutionary aspects of retroposition in murine and human populations that clearly documents the frequent integration of mRNA molecules into genomes, including clinical contexts. By performing basic comparisons, I show that the sequence features of mRNA vaccines meet all known requirements for retroposition using L1 elements-the most abundant autonomously active retrotransposons in the human genome. In fact, many factors associated with mRNA vaccines increase the possibility of their L1-mediated retroposition. I conclude that is unfounded to a priori assume that mRNA-based therapeutics do not impact genomes and that the route to genome integration of vaccine mRNAs via endogenous L1 retroelements is easily conceivable. This implies that we urgently need experimental studies that would rigorously test for the potential retroposition of vaccine mRNAs. At present, the insertional mutagenesis safety of mRNA-based vaccines should be considered unresolved.

Virus-Induced Membrane Fusion in Neurodegenerative Disorders

A growing body of epidemiological and research data has associated neurotropic viruses with accelerated brain aging and increased risk of neurodegenerative disorders. Many viruses replicate optimally in senescent cells, as they offer a hospitable microenvironment with persistently elevated cytosolic calcium, abundant intracellular iron, and low interferon type I. As cell-cell fusion is a major driver of cellular senescence, many viruses have developed the ability to promote this phenotype by forming syncytia. Cell-cell fusion is associated with immunosuppression mediated by phosphatidylserine externalization that enable viruses to evade host defenses. In hosts, virus-induced immune dysfunction and premature cellular senescence may predispose to neurodegenerative disorders. This concept is supported by novel studies that found postinfectious cognitive dysfunction in several viral illnesses, including human immunodeficiency virus-1, herpes simplex virus-1, and SARS-CoV-2. Virus-induced pathological syncytia may provide a unified framework for conceptualizing neuronal cell cycle reentry, aneuploidy, somatic mosaicism, viral spreading of pathological Tau and elimination of viable synapses and neurons by neurotoxic astrocytes and microglia. In this narrative review, we take a closer look at cell-cell fusion and vesicular merger in the pathogenesis of neurodegenerative disorders. We present a "decentralized" information processing model that conceptualizes neurodegeneration as a systemic illness, triggered by cytoskeletal pathology. We also discuss strategies for reversing cell-cell fusion, including, TMEM16F inhibitors, calcium channel blockers, senolytics, and tubulin stabilizing agents. Finally, going beyond neurodegeneration, we examine the potential benefit of harnessing fusion as a therapeutic strategy in regenerative medicine.

Intracellular Reverse Transcription of Pfizer BioNTech COVID-19 mRNA Vaccine BNT162b2 In Vitro in Human Liver Cell Line

Preclinical studies of COVID-19 mRNA vaccine BNT162b2, developed by Pfizer and BioNTech, showed reversible hepatic effects in animals that received the BNT162b2 injection. Furthermore, a recent study showed that SARS-CoV-2 RNA can be reverse-transcribed and integrated into the genome of human cells. In this study, we investigated the effect of BNT162b2 on the human liver cell line Huh7 in vitro. Huh7 cells were exposed to BNT162b2, and quantitative PCR was performed on RNA extracted from the cells. We detected high levels of BNT162b2 in Huh7 cells and changes in gene expression of long interspersed nuclear element-1 (LINE-1), which is an endogenous reverse transcriptase. Immunohistochemistry using antibody binding to LINE-1 open reading frame-1 RNA-binding protein (ORFp1) on Huh7 cells treated with BNT162b2 indicated increased nucleus distribution of LINE-1. PCR on genomic DNA of Huh7 cells exposed to BNT162b2 amplified the DNA sequence unique to BNT162b2. Our results indicate a fast up-take of BNT162b2 into human liver cell line Huh7, leading to changes in LINE-1 expression and distribution. We also show that BNT162b2 mRNA is reverse transcribed intracellularly into DNA in as fast as 6 h upon BNT162b2 exposure.

Transcription of the Envelope Protein by 1-L Protein-RNA Recognition Code Leads to Genes/Proteins That Are Relevant to the SARS-CoV-2 Life Cycle and Pathogenesis

The theoretical protein-RNA recognition code was used in this study to research the compatibility of the SARS-CoV-2 envelope protein (E) with mRNAs in the human transcriptome. According to a review of the literature, the spectrum of identified genes showed that the virus post-transcriptionally promotes or represses the genes involved in the SARS-CoV-2 life cycle. The identified genes/proteins are also involved in adaptive immunity, in the function of the cilia and wound healing (EMT and MET) in the pulmonary epithelial tissue, in Alzheimer's and Parkinson's disease and in type 2 diabetes. For example, the E-protein promotes BHLHE40, which switches off the IL-10 inflammatory "brake" and inhibits antiviral THαβ cells. In the viral cycle, E supports the COPII-SCAP-SREBP-HSP90α transport complex by the lowering of cholesterol in the ER and by the repression of insulin signaling, which explains the positive effect of HSP90 inhibitors in COVID-19 (geldanamycin), and E also supports importin α/β-mediated transport to the nucleus, which explains the positive effect of ivermectin, a blocker of importins α/β. In summary, transcription of the envelope protein by the 1-L protein-RNA recognition code leads to genes/proteins that are relevant to the SARS-CoV-2 life cycle and pathogenesis.

Competing endogenous RNA network mediated by circ_3205 in SARS-CoV-2 infected cells

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) is a new member of the Betacoronaviridae family, responsible for the recent pandemic outbreak of COVID-19. To start exploring the molecular events that follow host cell infection, we queried VirusCircBase and identified a circular RNA (circRNA) predicted to be synthesized by SARS-CoV-2, circ_3205, which we used to probe: (i) a training cohort comprised of two pools of cells from three nasopharyngeal swabs of SARS-CoV-2 infected (positive) or uninfected (negative, UCs) individuals; (ii) a validation cohort made up of 12 positive and 3 negative samples. The expression of circRNAs, miRNAs and miRNA targets was assayed through real-time PCR. CircRNA-miRNA interactions were predicted by TarpMiR, Analysis of Common Targets for circular RNAs (ACT), and STarMir tools. Enrichment of the biological processes and the list of predicted miRNA targets were retrieved from DIANA miRPath v3.0. Our results showed that the predicted SARS-CoV-2 circ_3205 was expressed only in positive samples and its amount positively correlated with that of SARS-CoV-2 Spike (S) mRNA and the viral load (r values = 0.80952 and 0.84867, Spearman's correlation test, respectively). Human (hsa) miR-298 was predicted to interact with circ_3205 by all three predictive tools. KCNMB4 and PRKCE were predicted as hsa-miR-298 targets. Interestingly, the function of both is correlated with blood coagulation and immune response. KCNMB4 and PRKCE mRNAs were upregulated in positive samples as compared to UCs (6 and 8.1-fold, p values = 0.049 and 0.02, Student's t test, respectively) and their expression positively correlated with that of circ_3205 (r values = 0.6 and 0.25, Spearman's correlation test, respectively). We propose that our results convincingly suggest that circ_3205 is a circRNA synthesized by SARS-CoV-2 upon host cell infection and that it may behave as a competitive endogenous RNA (ceRNA), sponging hsa-miR-298 and contributing to the upregulation of KCNMB4 and PRKCE mRNAs.

The Inhibition of RNA Viruses by Amaryllidaceae Alkaloids: Opportunities for the Development of Broad-Spectrum Anti-Coronavirus Drugs

The global COVID-19 pandemic has claimed the lives of millions and disrupted nearly every aspect of human society. Currently, vaccines remain the only widely available medical means to address the cause of the pandemic, the SARS-CoV-2 virus. Unfortunately, current scientific consensus deems the emergence of vaccine-resistant SARS-CoV-2 variants highly likely. In this context, the design and development of broad-spectrum, small-molecule based antiviral drugs has been described as a potentially effective, alternative medical strategy to address circulating and re-emerging CoVs. Small molecules are well-suited to target the least-rapidly evolving structures within CoVs such as highly conserved RNA replication enzymes, and this renders them less vulnerable to evolved drug resistance. Examination of the vast literature describing the inhibition of RNA viruses by Amaryllidaceae alkaloids suggests that future, broad-spectrum anti-CoV drugs may be derived from this family of natural products.

Enterovirus Infection Restricts Long Interspersed Element 1 Retrotransposition

Long interspersed element 1 (LINE-1 or L1) is the only active autonomous retrotransposon in the human genome that can serve as an endogenous upstream activator of cytoplasmic nucleic acid sensing pathways to elicit an antiviral immune response. In this study, we investigated the influence of enteroviral infection on L1 mobility. The results showed that infection with different enteroviruses, both EV-D68 and EV-A71, blocked L1 transposition. We screened diverse viral accessory proteins for L1 activity and identified EV-D68 2A, 3A, 3C, and EV-A71 ORF2p proteins as viral L1 inhibitors. EV-D68 2A suppressed L1 mobility by expression suppression of L1 proteins. Viral proteins 3A and 3C restricted ORF2p-mediated L1 reverse transcription in isolated L1 ribonucleoproteins. The newly identified enteroviral protein ORF2p inhibited the expression of L1 ORF1p. Altogether, our findings shed light on the strict modulation of L1 retrotransposons during enterovirus replication.

Integration of SARS-CoV-2 RNA in infected human cells by retrotransposons: an unlikely hypothesis and old viral relationships

Zhang et al. (Proc Natl Acad Sci 118:e2105968118, 2021) recently reported that SARS-CoV-2 RNA can be retrotranscribed and integrated into the DNA of human cells by the L1 retrotransposon machinery. This phenomenon could cause persistence of viral sequences in patients and may explain the prolonged PCR-positivity of SARS-CoV-2 infected patients, even long after the phase of active virus replication has ended. This commentary does critically review the available data on this topic and discusses them in the context of findings made for other exogenous viruses and ancestral endogenous retroviral elements.

COVID-19 in Children: Expressions of Type I/II/III Interferons, TRIM28, SETDB1, and Endogenous Retroviruses in Mild and Severe Cases

Children with the new coronavirus disease 2019 (COVID-19) have milder symptoms and a better prognosis than adult patients. Several investigations assessed type I, II, and III interferon (IFN) signatures in SARS-CoV-2 infected adults, however no data are available for pediatric patients. TRIM28 and SETDB1 regulate the transcription of multiple genes involved in the immune response as well as of human endogenous retroviruses (HERVs). Exogenous viral infections can trigger the activation of HERVs, which in turn can induce inflammatory and immune reactions. Despite the potential cross-talks between SARS-CoV-2 infection and TRIM28, SETDB1, and HERVs, information on their expressions in COVID-19 patients is lacking. We assessed, through a PCR real time Taqman amplification assay, the transcription levels of six IFN-I stimulated genes, IFN-II and three of its sensitive genes, three IFN-lIIs, as well as of TRIM28, SETDB1, pol genes of HERV-H, -K, and -W families, and of env genes of Syncytin (SYN)1, SYN2, and multiple sclerosis-associated retrovirus (MRSV) in peripheral blood from COVID-19 children and in control uninfected subjects. Higher expression levels of IFN-I and IFN-II inducible genes were observed in 36 COVID-19 children with mild or moderate disease as compared to uninfected controls, whereas their concentrations decreased in 17 children with severe disease and in 11 with multisystem inflammatory syndrome (MIS-C). Similar findings were found for the expression of TRIM-28, SETDB1, and every HERV gene. Positive correlations emerged between the transcriptional levels of type I and II IFNs, TRIM28, SETDB1, and HERVs in COVID-19 patients. IFN-III expressions were comparable in each group of subjects. This preserved induction of IFN-λs could contribute to the better control of the infection in children as compared to adults, in whom IFN-III deficiency has been reported. The upregulation of IFN-I, IFN-II, TRIM28, SETDB1, and HERVs in children with mild symptoms, their declines in severe cases or with MIS-C, and the positive correlations of their transcription in SARS-CoV-2-infected children suggest that they may play important roles in conditioning the evolution of the infection.

Host-Virus Chimeric Events in SARS-CoV-2-Infected Cells Are Infrequent and Artifactual

The pathogenic mechanisms underlying severe SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) infection remain largely unelucidated. High-throughput sequencing technologies that capture genome and transcriptome information are key approaches to gain detailed mechanistic insights from infected cells. These techniques readily detect both pathogen- and host-derived sequences, providing a means of studying host-pathogen interactions. Recent studies have reported the presence of host-virus chimeric (HVC) RNA in transcriptome sequencing (RNA-seq) data from SARS-CoV-2-infected cells and interpreted these findings as evidence of viral integration in the human genome as a potential pathogenic mechanism. Since SARS-CoV-2 is a positive-sense RNA virus that replicates in the cytoplasm, it does not have a nuclear phase in its life cycle. Thus, it is biologically unlikely to be in a location where splicing events could result in genome integration. Therefore, we investigated the biological authenticity of HVC events. In contrast to true biological events like mRNA splicing and genome rearrangement events, which generate reproducible chimeric sequencing fragments across different biological isolates, we found that HVC events across >100 RNA-seq libraries from patients with coronavirus disease 2019 (COVID-19) and infected cell lines were highly irreproducible. RNA-seq library preparation is inherently error prone due to random template switching during reverse transcription of RNA to cDNA. By counting chimeric events observed when constructing an RNA-seq library from human RNA and spiked-in RNA from an unrelated species, such as the fruit fly, we estimated that ∼1% of RNA-seq reads are artifactually chimeric. In SARS-CoV-2 RNA-seq, we found that the frequency of HVC events was, in fact, not greater than this background "noise." Finally, we developed a novel experimental approach to enrich SARS-CoV-2 sequences from bulk RNA of infected cells. This method enriched viral sequences but did not enrich HVC events, suggesting that the majority of HVC events are, in all likelihood, artifacts of library construction. In conclusion, our findings indicate that HVC events observed in RNA-sequencing libraries from SARS-CoV-2-infected cells are extremely rare and are likely artifacts arising from random template switching of reverse transcriptase and/or sequence alignment errors. Therefore, the observed HVC events do not support SARS-CoV-2 fusion to cellular genes and/or integration into human genomes. IMPORTANCE The pathogenic mechanisms underlying SARS-CoV-2, the virus responsible for COVID-19, are not fully understood. In particular, relatively little is known about the reasons some individuals develop life-threatening or persistent COVID-19. Recent studies identified host-virus chimeric (HVC) reads in RNA-sequencing data from SARS-CoV-2-infected cells and suggested that HVC events support potential "human genome invasion" and "integration" by SARS-CoV-2. This suggestion has fueled concerns about the long-term effects of current mRNA vaccines that incorporate elements of the viral genome. SARS-CoV-2 is a positive-sense, single-stranded RNA virus that does not encode a reverse transcriptase and does not include a nuclear phase in its life cycle, so some doubts have rightfully been expressed regarding the authenticity of HVCs and the role played by endogenous retrotransposons in this phenomenon. Thus, it is important to independently authenticate these HVC events. Here, we provide several lines of evidence suggesting that the observed HVC events are likely artifactual.

SARS-CoV-2-Host Chimeric RNA-Sequencing Reads Do Not Necessarily Arise From Virus Integration Into the Host DNA

The human genome bears evidence of extensive invasion by retroviruses and other retroelements, as well as by diverse RNA and DNA viruses. High frequency of somatic integration of the RNA virus severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) into the DNA of infected cells was recently suggested, based on a number of observations. One key observation was the presence of chimeric RNA-sequencing (RNA-seq) reads between SARS-CoV-2 RNA and RNA transcribed from human host DNA. Here, we examined the possible origin specifically of human-SARS-CoV-2 chimeric reads in RNA-seq libraries and provide alternative explanations for their origin. Chimeric reads were frequently detected also between SARS-CoV-2 RNA and RNA transcribed from mitochondrial DNA or episomal adenoviral DNA present in transfected cell lines, which was unlikely the result of SARS-CoV-2 integration. Furthermore, chimeric reads between SARS-CoV-2 RNA and RNA transcribed from nuclear DNA were highly enriched for host exonic, rather than intronic or intergenic sequences and often involved the same, highly expressed host genes. Although these findings do not rule out SARS-CoV-2 somatic integration, they nevertheless suggest that human-SARS-CoV-2 chimeric reads found in RNA-seq data may arise during library preparation and do not necessarily signify SARS-CoV-2 reverse transcription, integration in to host DNA and further transcription.

Reverse-transcribed SARS-CoV-2 RNA can integrate into the genome of cultured human cells and can be expressed in patient-derived tissues

An unresolved issue of SARS-CoV-2 disease is that patients often remain positive for viral RNA as detected by PCR many weeks after the initial infection in the absence of evidence for viral replication. We show here that SARS-CoV-2 RNA can be reverse-transcribed and integrated into the genome of the infected cell and be expressed as chimeric transcripts fusing viral with cellular sequences. Importantly, such chimeric transcripts are detected in patient-derived tissues. Our data suggest that, in some patient tissues, the majority of all viral transcripts are derived from integrated sequences. Our data provide an insight into the consequence of SARS-CoV-2 infections that may help to explain why patients can continue to produce viral RNA after recovery.

Prolonged detection of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) RNA and recurrence of PCR-positive tests have been widely reported in patients after recovery from COVID-19, but some of these patients do not appear to shed infectious virus. We investigated the possibility that SARS-CoV-2 RNAs can be reverse-transcribed and integrated into the DNA of human cells in culture and that transcription of the integrated sequences might account for some of the positive PCR tests seen in patients. In support of this hypothesis, we found that DNA copies of SARS-CoV-2 sequences can be integrated into the genome of infected human cells. We found target site duplications flanking the viral sequences and consensus LINE1 endonuclease recognition sequences at the integration sites, consistent with a LINE1 retrotransposon-mediated, target-primed reverse transcription and retroposition mechanism. We also found, in some patient-derived tissues, evidence suggesting that a large fraction of the viral sequences is transcribed from integrated DNA copies of viral sequences, generating viral–host chimeric transcripts. The integration and transcription of viral sequences may thus contribute to the detection of viral RNA by PCR in patients after infection and clinical recovery. Because we have detected only subgenomic sequences derived mainly from the 3′ end of the viral genome integrated into the DNA of the host cell, infectious virus cannot be produced from the integrated subgenomic SARS-CoV-2 sequences.