Quantitative RNAseq analysis of Ugandan KS tumors reveals KSHV gene expression dominated by transcription from the LTd downstream latency promoter.

KSHV is endemic in Uganda and the HIV epidemic has dramatically increased the incidence of Kaposi sarcoma (KS). To investigate the role of KSHV in the development of KS, we obtained KS biopsies from ART-naïve, HIV-positive individuals in Uganda and analyzed the tumors using RNAseq to globally characterize the KSHV transcriptome. Phylogenetic analysis of ORF75 sequences from 23 tumors revealed 6 distinct genetic clusters with KSHV strains exhibiting M, N or P alleles. RNA reads mapping to specific unique coding sequence (UCDS) features were quantitated using a gene feature file previously developed to globally analyze and quantitate KSHV transcription in infected endothelial cells. A pattern of high level expression was detected in the KSHV latency region that was common to all KS tumors. The clear majority of transcription was derived from the downstream latency transcript promoter P3(LTd) flanking ORF72, with little evidence of transcription from the P1(LTc) latency promoter, which is constitutive in KSHV-infected lymphomas and tissue-culture cells. RNAseq data provided evidence of alternate P3(LTd) transcript editing, splicing and termination resulting in multiple gene products, with 90% of the P3(LTd) transcripts spliced to release the intronic source of the microRNAs K1-9 and 11. The spliced transcripts encode a regulatory uORF upstream of Kaposin A with alterations in intervening repeat sequences yielding novel or deleted Kaposin B/C-like sequences. Hierarchical clustering and PCA analysis of KSHV transcripts revealed three clusters of tumors with different latent and lytic gene expression profiles. Paradoxically, tumors with a latent phenotype had high levels of total KSHV transcription, while tumors with a lytic phenotype had low levels of total KSHV transcription. Morphologically distinct KS tumors from the same individual showed similar KSHV gene expression profiles suggesting that the tumor microenvironment and host response play important roles in the activation level of KSHV within the infected tumor cells.

Complete genome sequence of Pig-tailed macaque rhadinovirus 2 and its evolutionary relationship with rhesus macaque rhadinovirus and human herpesvirus 8/Kaposi's sarcoma-associated herpesvirus.

UNLABELLED: Two rhadinovirus lineages have been identified in Old World primates. The rhadinovirus 1 (RV1) lineage consists of human herpesvirus 8, Kaposi's sarcoma-associated herpesvirus (KSHV), and closely related rhadinoviruses of chimpanzees, gorillas, macaques and other Old World primates. The RV2 rhadinovirus lineage is distinct and consists of closely related viruses from the same Old World primate species. Rhesus macaque rhadinovirus (RRV) is the RV2 prototype, and two RRV isolates, 26-95 and 17577, were sequenced. We determined that the pig-tailed macaque RV2 rhadinovirus, MneRV2, is highly associated with lymphomas in macaques with simian AIDS. To further study the role of rhadinoviruses in the development of lymphoma, we sequenced the complete genome of MneRV2 and identified 87 protein coding genes and 17 candidate microRNAs (miRNAs). A strong genome colinearity and sequence homology were observed between MneRV2 and RRV26-95, although the open reading frame (ORF) encoding the KSHV ORFK15 homolog was disrupted in RRV26-95. Comparison with MneRV2 revealed several genomic anomalies in RRV17577 that were not present in other rhadinovirus genomes, including an N-terminal duplication in ORF4 and a recombinative exchange of more distantly related homologs of the ORF22/ORF47 interacting glycoprotein genes. The comparison with MneRV2 has revealed novel genes and important conservation of protein coding domains and transcription initiation, termination, and splicing signals, which have added to our knowledge of RV2 rhadinovirus genetics. Further comparisons with KSHV and other RV1 rhadinoviruses will provide important avenues for dissecting the biology, evolution, and pathology of these closely related tumor-inducing viruses in humans and other Old World primates. IMPORTANCE: This work provides the sequence characterization of MneRV2, the pig-tailed macaque homolog of rhesus rhadinovirus (RRV). MneRV2 and RRV belong to the rhadinovirus 2 (RV2) rhadinovirus lineage of Old World primates and are distinct but related to Kaposi's sarcoma-associated herpesvirus (KSHV), the etiologic agent of Kaposi's sarcoma. Pig-tailed macaques provide important models of human disease, and our previous studies have indicated that MneRV2 plays a causal role in AIDS-related lymphomas in macaques. Delineation of the MneRV2 sequence has allowed a detailed characterization of the genome structure, and evolutionary comparisons with RRV and KSHV have identified conserved promoters, splice junctions, and novel genes. This comparison provides insight into RV2 rhadinovirus biology and sets the groundwork for more intensive next-generation (Next-Gen) transcript and genetic analysis of this class of tumor-inducing herpesvirus. This study supports the use of MneRV2 in pig-tailed macaques as an important model for studying rhadinovirus biology, transmission and pathology.

Next-generation sequence analysis of the genome of RFHVMn, the macaque homolog of Kaposi's sarcoma (KS)-associated herpesvirus, from a KS-like tumor of a pig-tailed macaque.

The complete sequence of retroperitoneal fibromatosis-associated herpesvirus Macaca nemestrina (RFHVMn), the pig-tailed macaque homolog of Kaposi's sarcoma-associated herpesvirus (KSHV), was determined by next-generation sequence analysis of a Kaposi's sarcoma (KS)-like macaque tumor. Colinearity of genes was observed with the KSHV genome, and the core herpesvirus genes had strong sequence homology to the corresponding KSHV genes. RFHVMn lacked homologs of open reading frame 11 (ORF11) and KSHV ORFs K5 and K6, which appear to have been generated by duplication of ORFs K3 and K4 after the divergence of KSHV and RFHV. RFHVMn contained positional homologs of all other unique KSHV genes, although some showed limited sequence similarity. RFHVMn contained a number of candidate microRNA genes. Although there was little sequence similarity with KSHV microRNAs, one candidate contained the same seed sequence as the positional homolog, kshv-miR-K12-10a, suggesting functional overlap. RNA transcript splicing was highly conserved between RFHVMn and KSHV, and strong sequence conservation was noted in specific promoters and putative origins of replication, predicting important functional similarities. Sequence comparisons indicated that RFHVMn and KSHV developed in long-term synchrony with the evolution of their hosts, and both viruses phylogenetically group within the RV1 lineage of Old World primate rhadinoviruses. RFHVMn is the closest homolog of KSHV to be completely sequenced and the first sequenced RV1 rhadinovirus homolog of KSHV from a nonhuman Old World primate. The strong genetic and sequence similarity between RFHVMn and KSHV, coupled with similarities in biology and pathology, demonstrate that RFHVMn infection in macaques offers an important and relevant model for the study of KSHV in humans.

Macaque homologs of EBV and KSHV show uniquely different associations with simian AIDS-related lymphomas.

Two gammaherpesviruses, Epstein-Barr virus (EBV) (Lymphocryptovirus genus) and Kaposi's sarcoma-associated herpesvirus (KSHV) (Rhadinovirus genus) have been implicated in the etiology of AIDS-associated lymphomas. Homologs of these viruses have been identified in macaques and other non-human primates. In order to assess the association of these viruses with non-human primate disease, archived lymphoma samples were screened for the presence of macaque lymphocryptovirus (LCV) homologs of EBV, and macaque rhadinoviruses belonging to the RV1 lineage of KSHV homologs or the more distant RV2 lineage of Old World primate rhadinoviruses. Viral loads were determined by QPCR and infected cells were identified by immunolabeling for different viral proteins. The lymphomas segregated into three groups. The first group (n = 6) was associated with SIV/SHIV infections, contained high levels of LCV (1-25 genomes/cell) and expressed the B-cell antigens CD20 or BLA.36. A strong EBNA-2 signal was detected in the nuclei of the neoplastic cells in one of the LCV-high lymphomas, indicative of a type III latency stage. None of the lymphomas in this group stained for the LCV viral capsid antigen (VCA) lytic marker. The second group (n = 5) was associated with D-type simian retrovirus-2 (SRV-2) infections, contained high levels of RV2 rhadinovirus (9-790 genomes/cell) and expressed the CD3 T-cell marker. The third group (n = 3) was associated with SIV/SHIV infections, contained high levels of RV2 rhadinovirus (2-260 genomes/cell) and was negative for both CD20 and CD3. In both the CD3-positive and CD3/CD20-negative lymphomas, the neoplastic cells stained strongly for markers of RV2 lytic replication. None of the lymphomas had detectable levels of retroperitoneal fibromatosis herpesvirus (RFHV), the macaque RV1 homolog of KSHV. Our data suggest etiological roles for both lymphocryptoviruses and RV2 rhadinoviruses in the development of simian AIDS-associated lymphomas and indicate that the virus-infected neoplastic lymphoid cells are derived from different lymphocyte lineages and differentiation stages.

Consensus-degenerate hybrid oligonucleotide primers (CODEHOPs) for the detection of novel viruses in non-human primates.

Consensus-degenerate hybrid oligonucleotide primers (CODEHOPs) have proven to be a powerful tool for the identification of novel genes. CODEHOPs are designed from highly-conserved regions of multiply-aligned protein sequences from members of a gene family and are used in PCR amplification to identify distantly-related genes. The CODEHOP approach has been used to identify novel pathogens by targeting amino acid motifs conserved in specific pathogen families. We initiated a program utilizing the CODEHOP approach to develop PCR-based assays targeting a variety of viral families that are pathogens in non-human primates. We have also developed and further improved a computer program and website to facilitate the design of CODEHOP PCR primers. Here, we detail the method for the development of pathogen-specific CODEHOP PCR assays using the papillomavirus family as a target. Papillomaviruses constitute a diverse virus family infecting a wide variety of mammalian species, including humans and non-human primates. We demonstrate that our pan-papillomavirus CODEHOP assay is broadly reactive with all major branches of the virus family and show its utility in identifying a novel non-human primate papillomavirus in cynomolgus macaques.

The ORF59 DNA polymerase processivity factor homologs of Old World primate RV2 rhadinoviruses are highly conserved nuclear antigens expressed in differentiated epithelium in infected macaques.

BACKGROUND: ORF59 DNA polymerase processivity factor of the human rhadinovirus, Kaposi's sarcoma-associated herpesvirus (KSHV), is required for efficient copying of the genome during virus replication. KSHV ORF59 is antigenic in the infected host and is used as a marker for virus activation and replication. RESULTS: We cloned, sequenced and expressed the genes encoding related ORF59 proteins from the RV1 rhadinovirus homologs of KSHV from chimpanzee (PtrRV1) and three species of macaques (RFHVMm, RFHVMn and RFHVMf), and have compared them with ORF59 proteins obtained from members of the more distantly-related RV2 rhadinovirus lineage infecting the same non-human primate species (PtrRV2, RRV, MneRV2, and MfaRV2, respectively). We found that ORF59 homologs of the RV1 and RV2 Old World primate rhadinoviruses are highly conserved with distinct phylogenetic clustering of the two rhadinovirus lineages. RV1 and RV2 ORF59 C-terminal domains exhibit a strong lineage-specific conservation. Rabbit antiserum was developed against a C-terminal polypeptide that is highly conserved between the macaque RV2 ORF59 sequences. This anti-serum showed strong reactivity towards ORF59 encoded by the macaque RV2 rhadinoviruses, RRV (rhesus) and MneRV2 (pig-tail), with no cross reaction to human or macaque RV1 ORF59 proteins. Using this antiserum and RT-qPCR, we determined that RRV ORF59 is expressed early after permissive infection of both rhesus primary fetal fibroblasts and African green monkey kidney epithelial cells (Vero) in vitro. RRV- and MneRV2-infected foci showed strong nuclear expression of ORF59 that correlated with production of infectious progeny virus. Immunohistochemical studies of an MneRV2-infected macaque revealed strong nuclear expression of ORF59 in infected cells within the differentiating layer of epidermis corroborating previous observations that differentiated epithelial cells are permissive for replication of KSHV-like rhadinoviruses. CONCLUSION: The ORF59 DNA polymerase processivity factor homologs of the Old World primate RV1 and RV2 rhadinovirus lineages are phylogenetically distinct yet demonstrate similar expression and localization characteristics that correlate with their use as lineage-specific markers for permissive infection and virus replication. These studies will aid in the characterization of virus activation from latency to the replicative state, an important step for understanding the biology and transmission of rhadinoviruses, such as KSHV.

Experimental co-transmission of Simian Immunodeficiency Virus (SIV) and the macaque homologs of the Kaposi Sarcoma-Associated Herpesvirus (KSHV) and Epstein-Barr Virus (EBV).

Macaque RFHV and LCV are close homologs of human KSHV and EBV, respectively. No experimental model of RFHV has been developed due to the lack of a source of culturable infectious virus. Screening of macaques at the Washington National Primate Research Center detected RFHV in saliva of SIV-infected macaques from previous vaccine studies. A pilot experimental infection of two naïve juvenile pig-tailed macaques was initiated by inoculation of saliva from SIV-infected pig-tailed and cynomolgus macaque donors, which contained high levels of DNA (> 10(6) genomes/ml) of the respective species-specific RFHV strain. Both juvenile recipients developed SIV and RFHV infections with RFHV DNA detected transiently in saliva and/or PBMC around week 16 post-infection. One juvenile macaque was infected with the homologous RFHVMn from whole saliva of a pig-tailed donor, which had been inoculated into the cheek pouch. This animal became immunosuppressed, developing simian AIDS and was euthanized 23 weeks after inoculation. The levels of RFHV DNA in saliva and PBMC remained below the level of detection after week 17, showing no reactivation of the RFHVMn infection during the rapid development of AIDS. The other juvenile macaque was infected with the heterologous RFHVMf from i.v. inoculation of purified virions from saliva of a cynomolgus donor. The juvenile recipient remained immunocompetent, developing high levels of persistent anti-RFHV and -SIV antibodies. After the initial presence of RFHVMf DNA in saliva and PBMC decreased to undetectable levels by week 19, all attempts to reactivate the infection through additional inoculations, experimental infection with purified SRV-2 or SIV, or immunosuppressive treatments with cyclosporine or dexamethasone were unsuccessful. An heterologous LCV transmission was also detected in this recipient, characterized by continual high levels of LCVMf DNA from the cynomolgus donor in both saliva (> 10(6) genomes/ml) and PBMC (> 10(4) genomes/million cells), coupled with high levels of anti-LCV antibodies. The macaque was sacrificed 209 weeks after the initial inoculation. Low levels of LCVMf DNA were detected in salivary glands, tonsils and other lymphoid organs, while RFHVMf DNA was below the level of detection. These results show successful co-transmission of RFHV and LCV from saliva and demonstrate differential lytic activation of the different gammaherpesvirus lineages due to presumed differences in biology and tropism and control by the host immune system. Although this initial pilot transmission study utilized only two macaques, it provides the first evidence for experimental transmission of the macaque homolog of KSHV, setting the stage for larger transmission studies to examine the differential activation of rhadinovirus and lymphocryptovirus infections and the pathological effects of immunosuppression.

Macaque homologs of Kaposi's sarcoma-associated herpesvirus (KSHV) infect germinal center lymphoid cells, epithelial cells in skin and gastrointestinal tract and gonadal germ cells in naturally infected macaques.

We developed a set of rabbit antisera to characterize infections by the macaque RV2 rhadinovirus homologs of KSHV. We analyzed tissues from rhesus and pig-tailed macaques naturally infected with rhesus rhadinovirus (RRV) or Macaca nemestrina rhadinovirus 2 (MneRV2). Our study demonstrates that RV2 rhadinoviruses have a tropism for epithelial cells, lymphocytes and gonadal germ cells in vivo. We observed latent infections in both undifferentiated and differentiated epithelial cells with expression of the latency marker, LANA. Expression of the early (ORF59) and late (glycoprotein B) lytic markers were detected in highly differentiated cells in epithelial ducts in oral, renal, dermal and gastric mucosal tissue as well as differentiated germ cells in male and female gonads. Our data provides evidence that epithelial and germ cell differentiation in vivo induces rhadinovirus reactivation and suggests that infected epithelial and germ cells play a role in transmission and dissemination of RV2 rhadinovirus infections in vivo.

Quantitative Analysis of the KSHV Transcriptome Following Primary Infection of Blood and Lymphatic Endothelial Cells.

The transcriptome of the Kaposi's sarcoma-associated herpesvirus (KSHV/HHV8) after primary latent infection of human blood (BEC), lymphatic (LEC) and immortalized (TIME) endothelial cells was analyzed using RNAseq, and compared to long-term latency in BCBL-1 lymphoma cells. Naturally expressed transcripts were obtained without artificial induction, and a comprehensive annotation of the KSHV genome was determined. A set of unique coding sequence (UCDS) features and a process to resolve overlapping transcripts were developed to accurately quantitate transcript levels from specific promoters. Similar patterns of KSHV expression were detected in BCBL-1 cells undergoing long-term latent infections and in primary latent infections of both BEC and LEC cultures. High expression levels of poly-adenylated nuclear (PAN) RNA and spliced and unspliced transcripts encoding the K12 Kaposin B/C complex and associated microRNA region were detected, with an elevated expression of a large set of lytic genes in all latently infected cultures. Quantitation of non-overlapping regions of transcripts across the complete KSHV genome enabled for the first time accurate evaluation of the KSHV transcriptome associated with viral latency in different cell types. Hierarchical clustering applied to a gene correlation matrix identified modules of co-regulated genes with similar correlation profiles, which corresponded with biological and functional similarities of the encoded gene products. Gene modules were differentially upregulated during latency in specific cell types indicating a role for cellular factors associated with differentiated and/or proliferative states of the host cell to influence viral gene expression.

Genetic variability of the envelope gene of Type D simian retrovirus-2 (SRV-2) subtypes associated with SAIDS-related retroperitoneal fibromatosis in different macaque species.

BACKGROUND: D-type simian retrovirus-2 (SRV-2) causes an AIDS-like immune deficiency syndrome (SAIDS) in various macaque species. SAIDS is often accompanied by retroperitoneal fibromatosis (RF), an aggressive fibroproliferative disorder reminiscent of Kaposi's sarcoma in patients with HIV-induced AIDS. In order to determine the association of SRV-2 subtypes with SAIDS-RF, and study the evolution and transmission of SRV-2 in captive macaque populations, we have molecularly characterized the env gene of a number of SRV-2 isolates from different macaque species with and without RF. RESULTS: We sequenced the env gene from eighteen SRV-2 isolates and performed sequence comparisons and phylogenetic analyses. Our studies revealed the presence of six distinct subtypes of SRV-2, three of which were associated with SAIDS-RF cases. We found no association between SRV-2 subtypes and a particular macaque species. Little sequence variation was detected in SRV-2 isolates from the same individual, even after many years of infection, or from macaques housed together or related by descent from a common infected parent. Seventy-two amino acid changes were identified, most occurring in the larger gp70 surface protein subunit. In contrast to the lentiviruses, none of the amino acid variations involved potential N-linked glycosylation sites. Structural analysis of a domain within the gp22/gp20 transmembrane subunit that was 100% conserved between SRV-2 subtypes, revealed strong similarities to a disulfide-bonded loop that is crucial for virus-cell fusion and is found in retroviruses and filoviruses. CONCLUSION: Our study suggests that separate introductions of at least six parental SRV-2 subtypes into the captive macaque populations in the U.S. have occurred with subsequent horizontal transfer between macaque species and primate centers. No specific association of a single SRV-2 subtype with SAIDS-RF was seen. The minimal genetic variability of the env gene within a subtype over time suggests that a strong degree of adaptation to its primate host has occurred during evolution of the virus.

Conservation of the glycoprotein B homologs of the Kaposi׳s sarcoma-associated herpesvirus (KSHV/HHV8) and old world primate rhadinoviruses of chimpanzees and macaques.

The envelope-associated glycoprotein B (gB) is highly conserved within the Herpesviridae and plays a critical role in viral entry. We analyzed the evolutionary conservation of sequence and structural motifs within the Kaposi׳s sarcoma-associated herpesvirus (KSHV) gB and homologs of Old World primate rhadinoviruses belonging to the distinct RV1 and RV2 rhadinovirus lineages. In addition to gB homologs of rhadinoviruses infecting the pig-tailed and rhesus macaques, we cloned and sequenced gB homologs of RV1 and RV2 rhadinoviruses infecting chimpanzees. A structural model of the KSHV gB was determined, and functional motifs and sequence variants were mapped to the model structure. Conserved domains and motifs were identified, including an "RGD" motif that plays a critical role in KSHV binding and entry through the cellular integrin αVß3. The RGD motif was only detected in RV1 rhadinoviruses suggesting an important difference in cell tropism between the two rhadinovirus lineages.

Development of a real-time QPCR assay for the detection of RV2 lineage-specific rhadinoviruses in macaques and baboons.

BACKGROUND: Two distinct lineages of rhadinoviruses related to Kaposi's sarcoma-associated herpesvirus (KSHV/HHV8) have been identified in macaques and other Old World non-human primates. We have developed a real-time quantitative PCR (QPCR) assay using a TaqMan probe to differentially detect and quantitate members of the rhadinovirus-2 (RV2) lineage. PCR primers were derived from sequences within ORF 60 and the adjacent ORF 59/60 intergenic region which were highly conserved between the macaque RV2 rhadinoviruses, rhesus rhadinovirus (RRV) and Macaca nemestrina rhadinovirus-2 (MneRV2). These primers showed little similarity to the corresponding sequences of the macaque RV1 rhadinoviruses, retroperitoneal fibromatosis herpesvirus Macaca nemestrina (RFHVMn) and Macaca mulatta (RFHVMm). To determine viral loads per cell, an additional TaqMan QPCR assay was developed to detect the single copy cellular oncostatin M gene. RESULTS: We show that the RV2 QPCR assay is linear from less than 2 to more than 300,000 copies using MneRV2 DNA, and is non-reactive with RFHVMn DNA up to 1 billion DNA templates per reaction. RV2 loads ranging from 6 to 2,300 viral genome equivalent copies per 10(6) cells were detected in PBMC from randomly sampled macaques from the Washington National Primate Research Center. Screening tissue from other primate species, including another macaque, Macaca fascicularis, and a baboon, Papio cynocephalus, revealed the presence of novel rhadinoviruses, MfaRV2 and PcyRV2, respectively. Sequence comparison and phylogenetic analysis confirmed their inclusion within the RV2 lineage of KSHV-like rhadinoviruses. CONCLUSIONS: We describe a QPCR assay which provides a quick and sensitive method for quantitating rhadinoviruses belonging to the RV2 lineage of KSHV-like rhadinoviruses found in a variety of macaque species commonly used for biomedical research. While this assay broadly detects different RV2 rhadinovirus species, it is unreactive with RV1 rhadinovirus species which commonly co-infect the same primate hosts. We also show that this QPCR assay can be used to identify novel RV2 rhadinoviruses in different primate species.

High levels of retroperitoneal fibromatosis (RF)-associated herpesvirus in RF lesions in macaques are associated with ORF73 LANA expression in spindleoid tumour cells.

Two distinct lineages of rhadinoviruses related to Kaposi's sarcoma (KS)-associated herpesvirus (KSHV; Human herpesvirus 8), the causative agent of KS, have been identified. In macaques, the RV1 lineage is represented by retroperitoneal fibromatosis (RF) herpesvirus (RFHV), the homologue of KSHV, whilst the RV2 lineage is represented by rhesus rhadinovirus (RRV), a more distantly related virus. Real-time quantitative PCR was used to estimate the loads of RV1 and RV2 rhadinoviruses in simian acquired immunodeficiency syndrome-associated RF (SAIDS-RF), a neoplasm of macaques with similarities to AIDS-associated KS. Both RV1 and RV2 rhadinoviruses were detected in macaques with RF. The RV1 loads were 220- to 4300-fold higher in RF tumours than in spleen, showing a strong tumour association (mean loads of 1 800 000 vs 2900 copies per 10(6) cells in tumours and spleen, respectively). In contrast, RV2 loads in the RF tumours were 100-fold lower than RV1 loads and showed similar levels in tumours and spleen (mean loads of 16 000 vs 24 000 copies per 10(6) cells, respectively). Immunostaining with antibodies reactive against RFHV ORF73 latency-associated nuclear antigen (LANA) showed intense nuclear staining of the spindleoid RF tumour cells. Correlation of viral load and the number of LANA-positive cells indicated that RF tumour cells contained multiple copies of the RFHV genome per cell. This pattern of infectivity is similar to that seen in KS tumours latently infected with KSHV. Our study demonstrates similarities in the biology of KSHV and RFHV and supports a role for RFHV in the aetiology of SAIDS-RF.