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.

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.

Endogenous retroviruses (ERVs) are related to long terminal repeat (LTR) retrotransposons, comprising gene sequences of exogenous retroviruses integrated into the host genome and inherited according to Mendelian law. They are considered to have contributed greatly to the evolution of host genome structure and function. We previously characterized HERV-K HML-9 in the human genome. However, the biological function of this type of element in the genome of the chimpanzee, which is the closest living relative of humans, largely remains elusive. Therefore, the current study aims to characterize HML-9 in the chimpanzee genome and to compare the results with those in the human genome. Firstly, we report the distribution and genetic structural characterization of the 26 proviral elements and 38 solo LTR elements of HML-9 in the chimpanzee genome. The results showed that the distribution of these elements displayed a non-random integration pattern, and only six elements maintained a relatively complete structure. Then, we analyze their phylogeny and reveal that the identified elements all cluster together with HML-9 references and with those identified in the human genome. The HML-9 integration time was estimated based on the 2-LTR approach, and the results showed that HML-9 elements were integrated into the chimpanzee genome between 14 and 36 million years ago and into the human genome between 18 and 49 mya. In addition, conserved motifs, cis-regulatory regions, and enriched PBS sequence features in the chimpanzee genome were predicted based on bioinformatics. The results show that pathways significantly enriched for ERV LTR-regulated genes found in the chimpanzee genome are closely associated with disease development, including neurological and neurodevelopmental psychiatric disorders. In summary, the identification, characterization, and genomics of HML-9 presented here not only contribute to our understanding of the role of ERVs in primate evolution but also to our understanding of their biofunctional significance.

The HIV-1 provirus mainly consists of internal coding region flanked by the 2 same long terminal repeats (LTRs) at each terminus. The LTRs play important roles in HIV-1 reverse transcription, integration, and transcription by the association with host factors. However, despite of the significant study advances of the internal coding regions of HIV-1 by using definite reference classification, there are no systematic classifications for HIV-1 5’ LTRs, which hinders our elaboration on 5’ LTR and a better understanding of the viral origin, spread and therapy. Here, by analyzing all available resources of 5’ LTR sequences in public databases following 4 recognized principles for the reference classification, 83 representatives and 14 consensus sequences were identified as representatives of 2 groups, 6 subtypes, 6 sub-subtypes, and 9 CRFs. To test the reliability of our established classification system, the constructed references were applied to identify the 5’ LTR assignment of the 22 clinical isolates in China. The results revealed that 16 out of 22 tested strains showed a consistent subtype classification with the previous LTR-independent classification system. However, 6 strains, for which recombination events within 5’ LTR were demonstrated, unexpectedly showed a different subtype classification, leading a significant change of binding sites for important transcription factors including SP1, p53, and NF-κB. The binding change of these transcriptional factors would probably affect the transcriptional activity of 5’ LTR. This study established a unified classification system for HIV-1 5’ LTRs, which will facilitate HIV-1 characterization and be helpful for both basic and clinical research fields. IMPORTANCE Here, a scientific, reliable, and usable classification system based on the 5’ LTR sequences was established, which will allow us to effectively facilitate the precise typing of HIV-1 strains. This classification system was applied to 22 HIV-1 strains circulating in China, we found that 6 out of 22 strains analyzed, belonged to a different subtype when our results were compared to those obtained with the previous LTR-independent classification system. Thus, these data demonstrated that our classification method could greatly improve the HIV-1 subtype classification. We found that 6 5’ LTR sequences showed recombination events, leading to a significant exchange of the binding sites of transcriptional factors. Thus, this work established a comprehensive HIV-1 5’ LTR classification system, which will help the scientific community to precisely characterize HIV-1 variants, and better understand the origin and spread of HIV-1 strains, and it also may be helpful for pathogenicity and transmissibility evaluation studies.

HIV-1 provirus is flanked by one LTR at each terminal. The 5' LTR plays important roles in HIV-1 life cycle, especially, it determines HIV-1 transcription. However, there are 810 5' LTR entries exist in the HIV-1 sequence database, accounting for only 0.085% (810/949484). Here, we collected plasma samples from HIV-1-infected patients in Shenzhen province and got 219 5' LTR sequences. In addition, we found recombination in the LTR region. The recombinant (LS13145, LS11614, LS14862, LS14863) possess an insertion of CRF01_AE segment at HXB2 482-630bp (149bp) in the skeleton of 5' LTR of subtype C. At the same time, our study found that the occurrence of recombination caused changes in many TFBSs. As the increasing investigation on 5' LTRs diversity and characterization, we will get a deeper understanding of HIV-1 transmission, evolution, and the basic mechanism of transcriptional regulation.