Background. Interferon is a marker of host antiviral immunity, which is disordered in COVID-19 patients. ERV can affect the secretion of interferon through the cGAS-STING pathway. In this study, we explored whether IFN-I and HERV-K (HML-2) were activated in COVID-19 patients and whether there was an interaction between them. Methods. We collected blood samples from COVID-19 patients and healthy controls. We first detected the expression of HERV-K (HML-2) gag, env, and pol genes and IFN-I-related genes between patients and healthy people by qPCR, synchronously detected VERO cells infected with SARS-CoV-2. Then, the chromosome distributions of highly expressed HERV-K (HML-2) gag, env, and pol genes were mapped by the next-generation sequencing results, and GO analysis was performed on the related genes. Results. We found that the HERV-K (HML-2) gag, env, and pol genes were highly expressed in COVID-19 patients and VERO cells infected with SARS-CoV-2. The interferon-related genes IFNB1, ISG15, and IFIT1 were also activated in COVID-19 patients, and GO analysis showed that HERV-K (HML-2) can regulate the secretion of interferon. Conclusions. The high expression of HERV-K (HML-2) might activate the increase of interferon in COVID-19 patients, proving that HERV-K does not only play a negative role in the human body.

Human endogenous retroviruses (HERVs) are derived from the infection and integration of exogenetic retroviruses. HERVs account for 8% of human genome, and the majority of HERVs are solitary LTRs (solo-LTRs) due to homologous recombination. Multiple findings have showed that solo-LTRs could provide an enormous reservoir of transcriptional regulatory sequences involved in diverse biological processes, especially carcinogenesis and cancer development. The link between solo-LTRs and human diseases still remains poorly understood. This review focuses on the regulatory modules of solo-LTRs, which contribute greatly to the diversification and evolution of human genes. More importantly, although inactivating mutations, insertions and deletions have been identified in solo-LTRs, the inherited regulatory elements of solo-LTRs initiate the expression of chimeric lncRNA transcripts, which have been reported to play crucial roles in human health and disease. These findings provide valuable insights into the evolutionary and functional mechanisms underlying the presence of HERVs in human genome. Taken together, in this review, we will present evidences showing the regulatory and encoding capacity of solo-LTRs as well as the significant impact on various aspects of human biology.

Background: Human endogenous retroviruses (HERVs) result from ancestral infections caused by exogenous retroviruses that became incorporated into the germline DNA and evolutionarily fixed in the human genome. HERVs can be transmitted vertically in a Mendelian fashion and be stably maintained in the human genome, of which they are estimated to comprise approximately 8%. HERV-K (HML1-10) transcription has been confirmed to be associated with a variety of diseases, such as breast cancer, lung cancer, prostate cancer, melanoma, rheumatoid arthritis, and amyotrophic lateral sclerosis. However, the poor characterization of HML-9 prevents a detailed understanding of the regulation of the expression of this family in humans and its impact on the host genome. In light of this, a precise and updated HERV-K HML-9 genomic map is urgently needed to better evaluate the role of these elements in human health. Results: We report a comprehensive analysis of the presence and distribution of HERV-K HML-9 elements within the human genome, with a detailed characterization of the structural and phylogenetic properties of the group. A total of 23 proviruses and 47 solo LTR elements were characterized, with a detailed description of the provirus structure, integration time, potential regulated genes, transcription factor binding sites (TFBS), and primer binding site (PBS) features. The integration time results showed that the HML-9 elements found in the human genome integrated into the primate lineage between 17.5 and 48.5 million years ago (mya). Conclusion: The results provide a clear characterization of HML-9 and a comprehensive background for subsequent functional studies.

With increasing data on HIV-1, a more relevant molecular model describing mechanism details of HIV-1 genetic recombination usually requires upgrades. Currently an incomplete structural understanding of the copy choice mechanism along with several other issues in the field that lack elucidation led us to perform an analysis of the correlation between breakpoint distributions and (1) the probability of base pairing, and (2) intersubtype genetic similarity to further explore structural mechanisms. Near full length sequences of URFs from Asia, Europe, and Africa (one sequence/patient), and representative sequences of worldwide CRFs were retrieved from the Los Alamos HIV database. Their recombination patterns were analyzed by jpHMM in detail. Then the relationships between breakpoint distributions and (1) the probability of base pairing, and (2) intersubtype genetic similarities were investigated. Pearson correlation test showed that all URF groups and the CRF group exhibit the same breakpoint distribution pattern. Additionally, the Wilcoxon two-sample test indicated a significant and inexplicable limitation of recombination in regions with high pairing probability. These regions have been found to be strongly conserved across distinct biological states (i.e., strong intersubtype similarity), and genetic similarity has been determined to be a very important factor promoting recombination. Thus, the results revealed an unexpected disagreement between intersubtype similarity and breakpoint distribution, which were further confirmed by genetic similarity analysis. Our analysis reveals a critical conflict between results from natural HIV-1 isolates and those from HIV-1-based assay vectors in which genetic similarity has been shown to be a very critical factor promoting recombination. These results indicate the region with high-pairing probabilities may be a more fundamental factor affecting HIV-1 recombination than sequence similarity in natural HIV-1 infections. Our findings will be relevant in furthering the understanding of HIV-1 recombination mechanisms.

Background: Endogenous retrovirus (ERV) elements can act as proximal regulatory elements in promoting interferon (IFN) responses. Previous relevant studies have mainly focused on IFN-stimulated genes (ISGs). However, the role of ERV elements as cis-regulatory motifs in regulating genes of the JAK–STAT pathway remains poorly understood. In our study, we analyzed the changes in ERV elements and genes under both IFN stimulation and blockade of the signaling pathway. Methods: The effects of interferon on cells under normal conditions and knockout of the receptor were compared based on the THP1_IFNAR1_KO and THP1_IFNAR2_mutant cell lines. The correlation between differentially expressed ERVs (DHERVs) and differentially expressed genes (DEGs) as DEHERV-G pairs was explored with construction of gene regulatory networks related to ERV and induced by proinflammatory cytokines. Results: A total of 430 DEHERV loci and 190 DEGs were identified in 842 DEHERV-G pairs that are common to the three groups. More than 87% of DEHERV-G pairs demonstrated a consistent expression pattern. ISGs such as AIM2, IFIT1, IFIT2, IFIT3, STAT1, and IRF were activated via the JAK–STAT pathway in response to interferon stimulation. Thus, STAT1, STAT2, and IRF1 appear to play core roles in regulatory networks and are closely associated with ERVs. Conclusions: The RNA expression of ISGs and ERV elements is correlated, indicating that ERV elements are closely linked to host innate immune responses.

Human endogenous retroviruses (HERV), a type of transposable elements (TE), play crucial roles in human placental morphogenesis, immune response, and cancer progression. Emerging evidence suggests that TEs have been a rich source of regulatory elements in the human genome, but little is known about the global impact of HERVs on transcriptional networks in cancer. Using genome-wide approaches, we show that HERVs are composed primarily of three ancient superfamilies: ERVL-MaLR, ERVL, and ERV1. This analysis suggests that the integration of exonic, intronic, and intergenic HERVs, as well as human or Hominidae gene-specific HERVs, contributes to human genomic innovation. HERVs exonized in genes are located mainly in the 3′ untranslated region (UTR) or 3′ end and participate in basic biological processes. Active HERVs are located mainly in intronic and intergenic regions and tend to function as enhancers and contribute to cancer cell type–specific gene expression. More importantly, HERVs may also define chromatin topologically associating domain (TAD) and loop boundaries in a cell type–specific manner. Taken together, these findings reveal that ancient HERV elements are a source of diverse regulatory sequences, including 3′ UTRs, 5′ UTRs, promoters, and enhancers, and they contribute to genetic innovation and cancer cell type–specific gene expression, highlighting the previously underestimated importance of these elements. Significance: Genome-wide analyses show that human endogenous retroviruses mediate cancer cell type–specific gene expression, epigenetic modification, and 3D chromatin architecture, elucidating the relationship between HERVs and diverse cancers.