Identification and spontaneous immune targeting of an endogenous retrovirus K envelope protein in the Indian rhesus macaque model of human disease

Contemporary Distribution, Estimated Age, and Prehistoric Migrations of Old World Monkey Retroviruses

Old World monkeys (OWM), simians inhabiting Africa and Asia, are currently affected by at least four infectious retroviruses, namely, simian foamy virus (SFV), simian immunodeficiency virus (SIV), simian T-lymphotropic virus (STLV), and simian type D retrovirus (SRV). OWM also show chromosomal evidence of having been infected in the past with four more retroviral species, baboon endogenous virus (BaEV), Papio cynocephalus endogenous virus (PcEV), simian endogenous retrovirus (SERV), and Rhesus endogenous retrovirus-K (RhERV-K/SERV-K1). For some of the viruses, transmission to other primates still occurs, resulting, for instance, in the HIV pandemic. Retroviruses are intimately connected with their host as they are normally spread by close contact. In this review, an attempt to reconstruct the distribution and history of OWM retroviruses will be made. A literature overview of the species infected by any of the eight retroviruses as well as an age estimation of the pathogens will be given. In addition, primate genomes from databases have been re-analyzed for the presence of endogenous retrovirus integrations. Results suggest that some of the oldest retroviruses, SERV and PcEV, have travelled with their hosts to Asia during the Miocene, when a higher global temperature allowed simian expansions. In contrast, younger viruses, such as SIV and SRV, probably due to the lack of a primate continuum between the continents in later times, have been restricted to Africa and Asia, respectively.

Transcriptional analysis of multiple ovarian cancer cohorts reveals prognostic and immunomodulatory consequences of ERV expression

BACKGROUND

Endogenous retroviruses (ERVs) play a role in a variety of biological processes, including embryogenesis and cancer. DNA methyltransferase inhibitors (DNMTi)-induced ERV expression triggers interferon responses in ovarian cancer cells via the viral sensing machinery. Baseline expression of ERVs also occurs in cancer cells, though this process is poorly understood and previously unexplored in epithelial ovarian cancer (EOC). Here, the prognostic and immunomodulatory consequences of baseline ERV expression was assessed in EOC.

METHODS

ERV expression was assessed using EOC transcriptional data from The Cancer Genome Atlas (TCGA) and from an independent cohort (Hammersmith Hospital, HH), as well as from untreated or DNMTi-treated EOC cell lines. Least absolute shrinkage and selection operator (LASSO) logistic regression defined an ERV expression score to predict patient prognosis. Immunohistochemistry (IHC) was conducted on the HH cohort. Combination of DNMTi treatment with γδ T cells was tested in vitro, using EOC cell lines and patient-derived tumor cells.

RESULTS

ERV expression was found to define clinically relevant subsets of EOC patients. An ERV prognostic score was successfully generated in TCGA and validated in the independent cohort. In EOC patients from this cohort, a high ERV score was associated with better survival (log-rank p=0.0009) and correlated with infiltration of CD8+PD1+T cells (r=0.46, p=0.0001). In the TCGA dataset, a higher ERV score was found in BRCA1/2 mutant tumors, compared to wild type (p=0.015), while a lower ERV score was found in CCNE1 amplified tumors, compared to wild type (p=0.019). In vitro, baseline ERV expression dictates the level of ERV induction in response to DNMTi. Manipulation of an ERV expression threshold by DNMTi resulted in improved EOC cell killing by cytotoxic immune cells.

CONCLUSIONS

These findings uncover the potential for baseline ERV expression to robustly inform EOC patient prognosis, influence tumor immune infiltration and affect antitumor immunity.

In vitro and In vivo Susceptibility of Baboons (Papio sp.) to Infection with and Apparent Antibody Reactivity to Simian Betaretrovirus (SRV)

Despite the lack of confirmed reports of an exogenous Simian betaretrovirus (SRV) isolated from baboons (Papio sp.), reports of simian endogenous gammaretrovirus (SERV) in baboons with complete genomes suggest that such viruses may be potentially infectious. In addition, serologic tests have repeatedly demonstrated antibody reactivity to SRV in baboons from multiple colonies. These findings complicate the management and use of such animals for research. To provide further insight into this situation, we performed in vitro and in vivo studies to determine if baboons are or can be infected with SRV. In our initial experiment, we were not able to isolate SRV from 6 seropositive or sero-indeterminate baboons by coculturing their peripheral blood mononuclear cells (PBMC) with macaque PBMC or permissive cell lines. In a subsequent experiment, we found that baboon PBMC infected in vitro with high dose SRV were permissive to virus replication. To test in vivo infectibil- ity, groups of naive baboons were infused intravenously with either (i) the same SRV tissue culture virus stocks used for the in vitro studies, (ii) SRV antibody positive and PCR positive macaque blood, (iii) SRV antibody positive or indeterminate, but PCR negative baboon blood, or (iv) SRV antibody and PCR negative baboon blood. Sustained SRV infection, as defined by reproducible PCR detection and/or antibody seroconversion, was confirmed in 2 of 3 baboons receiving tissue culture virus but not in any recipients of transfused blood from seropositive macaques or baboons. In conclusion, the data indicate that even though baboon cells can be infected experimentally with high doses of tissue culture grown SRV, baboons that are repeatedly SRV antibody positive and PCR negative are unlikely to be infected with exogenous SRV and thus are unlikely to transmit a virus that would threaten the SPF status of captive baboon colonies.

Genome-wide characterization of endogenous retroviruses in snub-nosed monkeys

Background

Endogenous retroviruses (ERV) are remnants of former exogenous retroviruses that had previously invaded the germ line of the host that can be vertically transmitted across generations. While the majority of ERVs lack infectious capacity due to the accumulation of deleterious mutations, some ERVs remain active and produce potentially infectious viral particles. ERV sequences have been reported in all mammals; however, the distribution and diversity of ERVs in several primate taxa remains unclear. The aim of this study was to identify and classify the ERV sequences in the genomes of the golden snub-nosed monkey (Rhinopithecus roxellana) and the black and white snub-nosed monkey (Rhinopithecus bieti), two endangered primate species that exploit high altitude (2,500–4,500 m) temperate forests in southern and central China.

Methods

We used a TBLASTN program to search the ERV sequences of golden snub-nosed monkey genome and the black and white snub-nosed monkey genome. We retrieved all complete accession sequences from the homology search and then used the program, RetroTector, to check and identify the ERV sequences.

Results

We identified 284 and 263 endogenous retrovirus sequences in R. roxellana and R. bieti respectively. The proportion of full-length sequences of all ERV was 30% in R. roxellana and 21% in R. bieti and they were described as class I and class II or gamma-retrovirus and beta-retrovirus genera. The truncation pattern distribution in the two species was virtually identical. By analyzing and comparing ERV orthologues among 6 primate species, we identified the co-evolution of ERVs with their host. We also examined ERV-like sequences and found 48 such genes in R. roxellana and 63 in R. bieti. Some of those genes are associated with diseases, suggesting that ERVs might have involved the abnormal expression of certain genes that have contributed to deleterious consequences for the host.

Conclusions

Our results indicate that ERV sequences are widely distributed in snub-nosed monkeys, and their phylogenetic history can mirror that of their hosts over long evolutionary time scales. In addition, ERV sequences appear to have an important influence on the evolution of host pathology.

Endogenous retrovirus expression is associated with response to immune checkpoint blockade in clear cell renal cell carcinoma

Although a subset of clear cell renal cell carcinoma (ccRCC) patients respond to immune checkpoint blockade (ICB), predictors of response remain uncertain. We investigated whether abnormal expression of endogenous retroviruses (ERVs) in tumors is associated with local immune checkpoint activation (ICA) and response to ICB. Twenty potentially immunogenic ERVs (πERVs) were identified in ccRCC in The Cancer Genome Atlas data set, and tumors were stratified into 3 groups based on their expression levels. πERV-high ccRCC tumors showed increased immune infiltration, checkpoint pathway upregulation, and higher CD8+ T cell fraction in infiltrating leukocytes compared with πERV-low ccRCC tumors. Similar results were observed in ER+/HER2- breast, colon, and head and neck squamous cell cancers. ERV expression correlated with expression of genes associated with histone methylation and chromatin regulation, and πERV-high ccRCC was enriched in BAP1 mutant tumors. ERV3-2 expression correlated with ICA in 11 solid cancers, including the 4 named above. In a small retrospective cohort of 24 metastatic ccRCC patients treated with single-agent PD-1/PD-L1 blockade, ERV3-2 expression in tumors was significantly higher in responders compared with nonresponders. Thus, abnormal expression of πERVs is associated with ICA in several solid cancers, including ccRCC, and ERV3-2 expression is associated with response to ICB in ccRCC.

Emerging diagnostic challenges and characteristics of simian betaretrovirus infections in captive macaque colonies

To better understand Simian betaretrovirus (SRV) seropositivity in virus-negative macaques, we transfused blood from SRV-infected or suspect donors into immunosuppressed naive recipients. Our results do not support typical SRV1-5 infection as the cause, but provide evidence for several possibilities including serological artifact, new/different SRV, or an endogenous virus.

Identification of novel simian endogenous retroviruses that are indistinguishable from simian retrovirus (SRV) on current SRV diagnostic assays

Simian betaretroviruses include the well-known exogenous simian retroviruses (SRV-1 through SRV-8), and some closely related simian endogenous retroviruses (SERV). Here, we characterized two new viral genomes, which appear to represent novel SERVs but have characteristics of both SRV and SERV highlighting the need to develop new assays providing molecular and serologic differentiation of SERV and SRV to avoid false positives.

A computational method for prediction of matrix proteins in endogenous retroviruses

Human endogenous retroviruses (HERVs) encode active retroviral proteins, which may be involved in the progression of cancer and other diseases. Matrix protein (MA), in group-specific antigen genes (gag) of retroviruses, is associated with the virus envelope glycoproteins in most mammalian retroviruses and may be involved in virus particle assembly, transport and budding. However, the amount of annotated MAs in ERVs is still at a low level so far. No computational method to predict the exact start and end coordinates of MAs in gags has been proposed yet. In this paper, a computational method to identify MAs in ERVs is proposed. A divide and conquer technique was designed and applied to the conventional prediction model to acquire better results when dealing with gene sequences with various lengths. Initiation sites and termination sites were predicted separately and then combined according to their intervals. Three different algorithms were applied and compared: weighted support vector machine (WSVM), weighted extreme learning machine (WELM) and random forest (RF). G − mean (geometric mean of sensitivity and specificity) values of initiation sites and termination sites under 5-fold cross validation generated by random forest models are 0.9869 and 0.9755 respectively, highest among the algorithms applied. Our prediction models combine RF & WSVM algorithms to achieve the best prediction results. 98.4% of all the collected ERV sequences with complete MAs (125 in total) could be predicted exactly correct by the models. 94,671 HERV sequences from 118 families were scanned by the model, 104 new putative MAs were predicted in human chromosomes. Distributions of the putative MAs and optimizations of model parameters were also analyzed. The usage of our predicting method was also expanded to other retroviruses and satisfying results were acquired.