HervD Atlas: a curated knowledgebase of associations between human endogenous retroviruses and diseases.

Human endogenous retroviruses (HERVs), as remnants of ancient exogenous retrovirus infected and integrated into germ cells, comprise ∼8% of the human genome. These HERVs have been implicated in numerous diseases, and extensive research has been conducted to uncover their specific roles. Despite these efforts, a comprehensive source of HERV-disease association still needs to be added. To address this gap, we introduce the HervD Atlas (https://ngdc.cncb.ac.cn/hervd/), an integrated knowledgebase of HERV-disease associations manually curated from all related published literature. In the current version, HervD Atlas collects 60 726 HERV-disease associations from 254 publications (out of 4692 screened literature), covering 21 790 HERVs (21 049 HERV-Terms and 741 HERV-Elements) belonging to six types, 149 diseases and 610 related/affected genes. Notably, an interactive knowledge graph that systematically integrates all the HERV-disease associations and corresponding affected genes into a comprehensive network provides a powerful tool to uncover and deduce the complex interplay between HERVs and diseases. The HervD Atlas also features a user-friendly web interface that allows efficient browsing, searching, and downloading of all association information, research metadata, and annotation information. Overall, the HervD Atlas is an essential resource for comprehensive, up-to-date knowledge on HERV-disease research, potentially facilitating the development of novel HERV-associated diagnostic and therapeutic strategies.

Increased prevalence and stable clustering rate of HIV-1 transmitted drug resistance strains after 'treat-all' in a megacity of China.

OBJECTIVES: The 'treat-all' strategy was implemented in Shenzhen, China in 2016. The effect of this extensive treatment on transmitted drug resistance (TDR) of HIV is unclear. METHODS: TDR analysis was performed, based on the partial HIV-1 pol gene obtained from the newly reported HIV-1 positive cases from 2011 to 2019 in Shenzhen, China. The HIV-1 molecular transmission networks were inferred to analyse the spread of TDR. Logistic regression was used to identify the potential risk factors with TDR mutations (TDRMs) to cluster. RESULTS: A total of 12 320 partial pol sequences were included in this study. The prevalence of TDR was 2.95% (363/12 320), which increased from 2.57% to 3.52% after 'treat-all'. The TDR prevalence was increased in populations with the characteristics of CRF07_BC, being single, educated to junior college level and above, MSM and male. The sensitivities of viruses to six antiretroviral drugs were decreased. The clustering rate of TDRMs remained stable, and the sequences in the three drug resistance transmission clusters (DRTCs) were mainly found during 2011-16. CRF07_BC and CRF55_01B were the factors associated with TDRMs clustering in the networks. CONCLUSIONS: The 'treat-all' strategy might have contributed to a small increase in TDR, while most of the TDRMs were distributed sporadically, which implies that the 'treat-all' strategy is helpful for the control of TDR in high-risk populations.

Comprehensive identification and characterization of the HERV-K (HML-9) group in the human genome.

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.

Genetic characterization of HIV-1 epidemic in Anhui Province, China.

BACKGROUND: Anhui Province in China is facing a severe HIV epidemic with an increasing number of newly diagnosed cases. METHODS: In this study, HIV genetic characteristics in the province were investigated. Newly reported HIV-positive individuals from 15 districts of Anhui Province were enrolled and interviewed. Total viral RNA was extracted from plasma isolated from blood samples. We amplified and sequenced an HIV pol fragment of the 1062 bp. The sequences were used for determination of HIV subtypes and the presence of drug resistance mutations. Transmission networks were constructed to explore possible relationships. And all of assembled partial pol genes were submitted to the Stanford HIV Drug Resistance Database website to find the transmitted drug resistance. RESULTS: Partial pol gene sequences were obtained from 486 cases. The results showed that MSM was the most dominant transmission route (253, 52.06%), followed by heterosexual transmission (210, 43.21%) and blood-borne transmission (1, 0.21%). Many subtypes were identified, including CRF01_AE (226, 46.50%), CRF07_BC (151, 31.07%), subtype B (28, 5.76%), CRF08_BC (20, 4.12%), CRF55_01B (15, 3.09%), CRF68_01B (7, 1.44%), CRF67_01B (3, 0.62%), CRF57_BC (2, 0.41%), CRF59_01B (2, 0.41%), CRF79_0107 (2, 0.41%), subtype C (2, 0.41%), CRF64_BC (1, 0.21%), and circulating recombinant forms (URFs) (27, 5.55%). Four transmission subnetworks containing high transmission risk individuals (with degree ≥4) were identified based on CRF01_AE and CRF07_BC sequences, including two CRF01_AE transmission subnetworks constituted by elderly people with average ages of 67.9 and 61.5 years. Infection occurred most likely through heterosexual transmission, while the other two CRF07_BC transmission subnetworks consist mainly of MSMs with average ages of 31.73 and 34.15. The level of HIV-transmitted drug resistance is 3.09%. CONCLUSIONS: The simultaneous spread of multiple HIV subtypes in Anhui province underscores that close surveillance of the local HIV epidemic is necessary. Furthermore, the elderly people were frequently involved, arguing for behaviour intervention in this specific population besides the MSM risk group.

High prevalence and viremia of human pegivirus 2 in the HIV-infected population in Honghe Prefecture, Yunnan Province.

Human pegivirus 2 (HPgV-2) is a recently recognized pegivirus of the family Flaviviridae. To investigate the epidemic features of HPgV-2 circulating in the human immunodeficiency virus (HIV)-infected population, we tested for antibodies and viral RNA of HPgV-2 and hepatitis C virus (HCV) with retrospective plasma samples collected from 771 HIV infections with multiple risk behaviors in Honghe Prefecture of Yunnan Province. A total of 195 subjects (25.29%) were seroreactive to HPgV-2, and 41 (5.32%) were RNA positive. Although the positive rate of HPgV-2 antibodies in HIV/HCV-coinfected individuals (27.69%) was significantly higher than that of HIV monoinfections (20.82%) (p = 0.036), this is the first report of HPgV-2 viremia in HIV-infected individuals without HCV infection and the presence of two HPgV-2 lineages in China. Our data indicate that HPgV-2 can also be transmitted sexually, which might be facilitated when combined with HCV infection, injecting drug use, and risky sexual behavior, which appear to have a synergistic effect on HPgV-2 infection. Phylogenetic analysis of 26 near-full-length genome sequences showed that the HPgV-2 strains in China are divided into two clusters.

Natural presence of the V179D and K103R/V179D mutations associated with resistance to nonnucleoside reverse transcriptase inhibitors in HIV-1 CRF65_cpx strains.

BACKGROUND: There is increasing evidence that HIV-1 genetic diversity can have an impact on drug resistance. The aim of this study is to investigate the epidemiological situation of CRF65_cpx and the impact of natural polymorphisms of this variant on genotypic resistance. METHODS: We used the BLAST search program followed by phylogenetic analysis to identify additional CRF65_cpx pol sequences from the Los Alamos HIV Sequence Database. Maximum likelihood phylogeny was estimated to clarify the epidemiological relationship of CRF65_cpx strains. Genotypic resistance was determined by submitting sequences to the Stanford HIV Drug Resistance Database. RESULTS: A total of 32 CRF65_cpx pol sequences were obtained. The CRF65_cpx strains were detected in seven provinces with large geographic distance. Yunnan CRF65_cpx sequences were mainly derived from a heterosexual risk group, whereas the CRF65_cpx sequences in other provinces were almost exclusively derived from an MSM population. With one exception of V179E, the other 31 strains harbored V179D mutation. The combination of V179D and K103R, conferring intermediate resistance to EFV and NVP, was detected in seven treatment-naive MSM patients. CONCLUSIONS: This study confirmed the expansion CRF65_cpx in China. Furthermore, we found the natural presence of the V179D and K103R/V179D mutations associated with resistance to NNRTIs in HIV-1 CRF65_cpx. Our findings highlight the contribution of polymorphic mutations to drug resistance and underscore the challenges in treating patients harboring CRF65_cpx strains.

Two-year cross-sectional studies reveal that single, young MSMs in Shenzhen, China are at high risk for HIV infection.

BACKGROUND: Shenzhen City is a rapidly growing area with a large number of floating populations, thus making it difficult to control HIV. Serial cross-sectional studies are helpful for the prediction of epidemiological tendency. In this study, two parallel cross-sectional studies were compared to explore changes in HIV epidemiology in Shenzhen, China. METHODS: Two hundred and fifty newly reported HIV-positive cases were randomly selected in Shenzhen City in 2013 and 2015. Socio-demographical information was collected with informed consent. Full-length gag and partial pol genes were amplified using nested RT-PCR followed by sequencing and phylogenetic analysis. The genotypes of anti-HIV drug resistance were also analyzed. The characteristics of the HIV epidemics of 2013 and 2015 were compared to identify patterns. RESULTS: The proportion of single, young MSMs dramatically increased in 2015 compared to 2013. Many subtypes, including CRF07_BC (36.4%), CRF01_AE (34.1%), CRF55_01B (10.2%), B (6.4%), CRF08_BC (3.4%), CRF59_01B (0.9%), C (0.7%), D (0.2%), CRF68_01B (0.2%), CRF67_01B (0.2%), and unique recombinant forms (URFs, 7.3%), were identified. Close phylogenetic relationships between strains prevalent in Shenzhen and other areas of China was observed. No epidemic cluster confined to single, young MSMs was identified. 0.4 and 2.8% of the strains contained transmitted drug-resistant mutations in 2013 and 2015, respectively. CONCLUSION: Although the interval period is short, changes in HIV epidemiology in Shenzhen City are distinct. Frequent surveillance of HIV epidemics in Shenzhen City is thus necessary. Single, young MSMs have become a high-risk population for HIV infection and should be considered as focus population for HIV prevention and behavior intervention in Shenzhen City.

Phylogenetic analysis highlights the role of older people in the transmission of HIV-1 in Fuyang, Anhui Province, China.

BACKGROUND: The proportion of older HIV-1 infected people in China has increased rapidly in recent years. Elucidation of the transmission characteristics of this high-risk population subgroup is helpful for the development of tailored interventions. METHODS: A phylogenetic analysis was performed that uses available HIV-1 pol sequences amplified with nested RT-PCR from plasma samples of all newly diagnosed participants spanning from October 2017 to September 2018 in Fuyang, Anhui Province. Transmission clusters were identified as two or more sequences that shared a corresponding node with an aLRT-SH value ≥90 in the maximum-likelihood phylogenetic tree and had an overall mean genetic distance of ≤1.5%. A local transmission cluster was defined as a cluster that had more than 80% of its sequences from Fuyang. The role of older people in local HIV-1 transmission was determined using an integration of molecular and demographic data. RESULTS: Of 362 available sequences, 14 subtypes, and 28 local transmission clusters were identified. It was found that the proportion of older people in the local transmission cluster (69/77, 89.61%) was much higher than that of younger people (46/114, 40.35%) (χ2 test, P < 0.001). In the pretreatment drug resistance analysis, the proportion of sequences with PDRMs in the local transmission cluster was not significantly different between the older people group (57.14%, 4/7) and non-old-aged group (11.11%, 1/9) (Fisher's exact test, P > 0.05). CONCLUSION: By combining phylogenetic analyses with demographic data, more detailed information was provided about the local transmission structure in Fuyang. These findings suggested that older people play an important role in local transmission, and more tailored interventions for this population subgroup are urgently needed.

Characterization of small genomic regions of the hepatitis B virus should be performed with more caution.

BACKGROUND: Hepatitis B virus is a hepatotropic DNA virus that reproduces via an RNA intermediate. It can lead to an increased risk of serious liver diseases such as hepatocellular carcinoma and is a serious threat to public health. Currently, the HBV are designated based on greater than 8% nucleotide variation along the whole genome. The recombination of HBV is very common, a large majority of which are recombinants between 2 genotypes. The current work aims to characterize a suspected recombinant involving 3 genotypes. METHODS: Fifty-seven HBV full-genome sequences were obtained from 57 patients co-infected with HBV and HIV-1 by amplification coupled with sequencing. JpHMM and RDP4 were used to perform recombination analysis respectively. The recombination results of a suspected 3-genotypic recombinant were further confirmed by both maximum likelihood phylogenetic tree and Mrbayes tree. RESULTS: JpHMM recombination analysis clearly indicated one 3-genotypic HBV recombinant composing of B/C/D. The genotype assignments are supported by significant posterior probabilities. The subsequent phylogenetic analysis of sub-regions derived from inferred breakpoints led to a disagreement on the assignment of D segment. Investigating the conflict, further exploration by RDP4 and phylogenies revealed that the jpHMM-derived 3-genotypic recombinant is actually a B/C genotypic recombinant with C fragment spanning 1899 to 2295 (jpHMM) or 1821 to 2199 (RDP4). CONCLUSIONS: The whole analysis indicated that (i) determination of small genomic regions should be performed with more caution, (ii) combinations of various recombination detection approaches conduce to obtain impartial results, and (iii) a unified system of nomenclature of HBV genotypes is necessary.

HIV-1 Thai B strain has spread out of former plasma donors into general population through sexual contact in Henan, China.

Henan, China is characterized by the outbreak of HIV epidemic of Thai B strain in former plasma donors in 1990s. After the forbidden of paid blood donation, whether Thai B strain will spread out of former plasma donors into sexual transmitted population is unknown. To answer the question, phylogenetic analysis was used to explore relationships of HIV strains circulating in those two populations in the study. HIV-1 sero-positive drug-naïve patients infected through sexual contact were enrolled into the study. Full length gag and pol genes were amplified with nested RT-PCR followed by sequencing and phylogenetic analysis. The genotypes of anti-HIV drug resistance were also analyzed with available pol genes. HIV subtypes were determined in 249 individuals from 288 participants. Subtype B was dominant (202/249, 81.1%), followed by CRF01_AE (25/249, 10.0%), CRF07_BC (14/249, 5.6%), C (4/249, 1.6%), URF (3/249, 1.2%), and CRF08_BC (1/249, 0.4%). Most of subtype B strains belong to Thailand B lineage. All of Thai B strains identified in sexual transmitted population intermixed with those from former blood donors in phylogenetic tree, suggesting close phylogenetic relationship between strains epidemic in those two populations. TDR was identified in 9.9% individuals. Thai B strain has spread out of former blood donors in Henan province. The finding will contribute to understanding the distribution and evolution of HIV-1 in Henan province and also provide clue to behavior change intervention.

Analysis of HIV-1 intersubtype recombination breakpoints suggests region with high pairing probability may be a more fundamental factor than sequence similarity affecting HIV-1 recombination.

BACKGROUND: 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. METHODS: 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. RESULTS: 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. CONCLUSIONS: 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.

Characterization of two HIV-1 infectors during initial antiretroviral treatment, and the emergence of phenotypic resistance in reverse transcriptase-associated mutation patterns.

BACKGROUND: Highly active antiretroviral therapy (HAART) is recommended to control the infection of HIV-1. HIV-1 drug resistance becomes an obstacle to HAART due to the accumulation of specific mutations in the RT coding region. The development of resistance mutations may be more complex than previously thought. METHODS: We followed two HIV-1 infectors from a HIV-1 drug resistance surveillance cohort in Henan province and evaluated CD4+ T-cell number and viral load thereafter at ten time-periods and characterized their reverse transcriptase-associated mutation patterns at each time point. Then we constructed the recombinant virus strains with these mutation patterns to mimick the viruses and test the phenotypic resistance caused by the mutation patterns on TZM-b1 cells. RESULTS: CD4+ T-cell number initially increased and then decreased rapidly, while viral load decreased and then dropped sharply during initial antiretroviral treatment. The number of mutations and the combination patterns of mutations increased over time. According to the phenotypic resistance performed by recombinant virus strains, VirusT215Y/V179E/Y181C/H221Y exhibited high levels of resistance to EFV (5.57-fold), and T215Y/V179E-containing virus increased 20.20-fold in AZT resistance (p < 0.01). VirusT215Y/V179E/Y181C increased markedly in EFV resistance (p < 0.01). The IC50 for VirusT215Y/V179E/H221Y was similar to that for VirusT215Y/V179E/Y181C. VirusT215Y/K103N/Y181C/H221Y induced a dramatic IC50 increase of all the four agents (Efavirenz EFV, Zidovudine AZT, Lamivudine 3TC, and Stavudine d4T) (p < 0.01). As for VirusT215Y/K103N/Y181C, only the IC50 of EFV was significantly increased. T215Y/K103N resulted in a 26.36-fold increase in EFV (p < 0.01). T215Y/K103N/H221Y significantly increased the resistance to AZT and 3TC. The IC50 of EFV with T215Y/V179E was lower than with T215Y/K103N (F = 93.10, P < 0.0001). With T215Y/V179E, Y181C significantly increase in EFV resistance, while the interaction between 181 and 221 in EFV was not statistically significant (F = 1.20, P = 0.3052). With T215Y/K103N, neither H221Y nor Y181C showed a significant increase in EFV resistance, but the interaction between 181 and 221 was statistically significant (F = 38.12, P = 0.0003). CONCLUSIONS: Data in this study suggests that pathways of viral evolution toward drug resistance appear to proceed through distinct steps and at different rates. Phenotypic resistance using recombinant virus strains with different combination of mutation patterns reveals that interactions among mutations may provide information on the impact of these mutations on drug resistance. All the result provides reference to optimize clinical treatment schedule.

A considerable proportion of CRF01_AE strains in China originated from circulating intrasubtype recombinant forms (CIRF).

BACKGROUND: In this study, the prevalence of HIV-1 CRF01_AE intrasubtype recombinants in China is estimated and their contributions to the epidemic are explored. METHODS: Available HIV-1 complete genomes of CRF01_AE were retrieved from the HIV database. The two alignments were evaluated with RDP3. Recombinants were defined as cases in which the recombination signal was supported by at least 3 methods with P-values of ≤0.05 after Bonferroni correction for multiple comparisons implemented in RDP3. Phylogenetic analysis was performed to further investigate the role of intrasubtype recombinants in epidemics. RESULTS: Here, 124 out of the 339 sequences from around the world (36.6 %) showed significant evidence of recombination. Here, 84 of these recombinants were from China, accounting for 54.9 % of local total sequences (84 out of 153). The results indicated non-negligible levels of intrasubtype recombination. Subsequent phylogenetic analysis indicated that a considerable proportion of CRF01_AE strains in China originated from circulating intrasubtype recombinant forms. Three large, well-supported intrasubtype recombinants clusters were identified here. Through a survey of risk factors and sampling cities and provinces, cluster I and cluster II were found to be prevalent primarily among men who have sex with men in major northern cities. Cluster III was prevalent among heterosexuals and intravenous drug users in southern and southwestern provinces. CONCLUSIONS: The current work highlighted the remarkable prevalence of intrasubtype recombination within the CRF01_AE epidemic and emphasized the value of intrasubtype recombinants, which came to circulate in the same manner as intersubtype recombinants.

Higher Expression of Human Endogenous Retrovirus-K was Observed in Peripheral B Lymphocytes of Leukemia and Lymphoma Patients.

Hematological malignant tumors (HMTs) are serious diseases that threaten human health and life with high mortality. Therefore, it is necessary to develop novel strategies for diagnosis and treatment. Human endogenous retroviruses (HERVs) have recently attracted increasing attention as potential targets for cancer diagnosis and therapy. In this study, we explored the association between HERV-K expression levels and HMTs development. Clinical data and peripheral blood samples were collected from 236 leukemia, 384 lymphoma patients, and 69 healthy controls. Quantitative polymerase chain reaction was used to detect the expression of HERV-K gag, pol, and env genes in peripheral blood mononuclear cells or different cell subpopulations. Differently expressed HERV-K genes were further tested by using deep sequencing method, and further analyzed with gene ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment. B cell- and T cell-related cytokines in patients were also detected by enzyme-linked immunosorbent assay (ELISA). The results showed that the expression levels of the HERV-K gag, pol, and env genes in patients were significantly higher than in healthy controls. There was a correlation between the expression level of HERV-K and the clinicopathological parameters of leukemia patients. HERV-K expression was increased in the B lymphocytes of leukemia and lymphoma patients, but not in the T cells or neutrophils. The GO and KEGG analyses showed that abnormal expression of the HERV-K locus in patients affected immune regulation. The analysis of cytokines proved that the B cell-related cytokines, including interleukin (IL)-1ß, IL-2, IL-4, IL-6, IL-10, tumor necrosis factor (TNF)-α, and interferon-gamma, were significantly decreased in patients, while the T cell-related cytokines, including IL-3, IL-12, and TNF-ß, were not significantly changed. In conclusion, HERV-K genes might participate in the occurrence and development of leukemia and lymphoma, and might be biomarkers for the detection or evaluation of leukemia and lymphoma.

Sequence Note: Characterization of Two HIV-1 Strains with Novel Unique Recombinant Genome in Hebei, China.

In China, the proportion of HIV-1 infections due to men who have sex with men (MSM) has increased rapidly. More and more new subtypes are found among the MSM population besides known CRF01_AE, CRF07_BC, and B. The co-circulation of several HIV subtypes in the same population provides the opportunity to develop a new circulating recombinant form (CRF) and unique recombinant form (URF). Here we reported two new URFs from two HIV-1 positive subjects infected through homosexual contact in Hebei, China. Phylogenetic and recombinant analyses based on the near full-length genome (NFLG) of the two URFs are the second-generation recombinant strains that originated from B, CRF01_AE, and CRF07_BC. The CRF01_AE segments in the genome of two URFs originated from cluster 4 of CRF01_AE strains, while the CRF07_BC segments were clustered with 07BC_N in the phylogenetic tree. The emergence of the novel CRF01_AE/CRF07_BC and CRF01_AE/B recombinant forms indicated the importance of the continuous monitoring of the HIV-1 epidemic and new URFs among the MSM population.

Comprehensive characterization of ERV-K (HML-8) in the chimpanzee genome revealed less genomic activity than humans.

Endogenous retroviruses (ERVs) originate from ancestral germline infections caused by exogenous retroviruses. Throughout evolution, they have become fixed within the genome of the animals into which they were integrated. As ERV elements coevolve with the host, they are normally epigenetically silenced and can become upregulated in a series of physiological and pathological processes. Generally, a detailed ERV profile in the host genome is critical for understanding the evolutionary history and functional performance of the host genome. We previously characterized and cataloged all the ERV-K subtype HML-8 loci in the human genome; however, this has not been done for the chimpanzee, the nearest living relative of humans. In this study, we aimed to catalog and characterize the integration of HML-8 in the chimpanzee genome and compare it with the integration of HML-8 in the human genome. We analyzed the integration of HML-8 and found that HML-8 pervasively invaded the chimpanzee genome. A total of 76 proviral elements were characterized on 23/24 chromosomes, including detailed elements distribution, structure, phylogeny, integration time, and their potential to regulate adjacent genes. The incomplete structure of HML-8 proviral LTRs will undoubtedly affect their activity. Moreover, the results indicated that HML-8 integration occurred before the divergence between humans and chimpanzees. Furthermore, chimpanzees include more HML-8 proviral elements (76 vs. 40) and fewer solo long terminal repeats (LTR) (0 vs. 5) than humans. These results suggested that chimpanzee genome activity is less than the human genome and that humans may have a better ability to shape and screen integrated proviral elements. Our work is informative in both an evolutionary and a functional context for ERVs.

Comprehensive Identification and Characterization of HML-9 Group in Chimpanzee Genome.

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.

A unified classification system for HIV-1 5' long terminal repeats.

The HIV-1 provirus mainly consists of internal coding region flanked by 1 long terminal repeats (LTRs) at each terminus. The LTRs play important roles in HIV-1 reverse transcription, integration, and transcription. However, despite of the significant study advances of the internal coding regions of HIV-1 by using definite reference classification, there are no systematic and phylogenetic 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 the supplemented 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 supplemented 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.

Near Full-Length Genome Characterization of Two Novel Unique Recombinants (CRF01_AE/CRF07_BC) in Beijing, China.

With the prevalence of human immunodeficiency virus type 1 (HIV-1) CRF01_AE and CRF07_BC subtypes in China, the co-circulation of multiple subtypes in the HIV-1-positive population may result in dual infection or superinfection in the population, leading to the emergence of unique recombinant forms (URFs) of the HIV-1 virus. In this study, two second-generation unique recombinant strains, BI0114 and BI0116, were identified, and their near full-length genome sequences were obtained. Recombination analysis showed that both sequences were isoforms of URF_0107, and they were second-generation unique recombinant strains formed by the recombination of CRF01_AE and CRF07_BC, with the isoforms being CRF01_AE and CRF0107_BC, respectively. The continued emergence of novel CRF01_AE/CRF07_BC recombinant strains suggests that the epidemiological, preventive, and control situation of HIV-1 is complex and that the relevant health authorities urgently need to establish responses to the challenges posed by changes in the pattern of strain recombination.

Identification of Two Novel HIV-1 Second-Generation Recombinant Forms (CRF01_AE/CRF07_BC) in Men Who Have Sex with Men in Hebei, China.

In recent years, men who have sex with men (MSM) have been identified as the primary source of HIV-1 transmission in Hebei Province, China. Co-circulation of multiple subtypes in HIV-1 positive MSM populations may contribute to the emergence of the second generation of recombinant HIV-1 strains, indicating the occurrence of dual infections or superinfections in MSM populations. Thus, the discovery of new recombinant strains is important to indicate the appearance of multiple infected individuals and the prevalence caused by changes in the parent strains. In this study, we present two new unique recombinant forms (URFs) from two HIV-1-positive subjects (HB070052 and HB070056) infected through homosexual contact in Hebei Province, China. The near full-length genome of the two URFs revealed that HB070052 was divided into seven segments by six breakpoints in the gag, pol, vif, and vpr genes; HB070056 was separated into five fragments by four breakpoints, with two regions of CRF07_BC inserted into a CRF01_AE backbone's gag, pol regions. The subregion tree showed CRF01_AE segments were traced back to the cluster 4 and 6 of the CRF01_AE phylogenetic tree, which were prevalent among HIV-1 infections through MSM in China. The continued emergence of the novel CRF01_AE/CRF07_BC recombinant forms indicates the HIV-1 epidemic is complex and long-term surveillance of recombinant strains is necessary among MSM in this region.