Predicting host-based, synthetic lethal antiviral targets from omics data.

Traditional antiviral therapies often have limited effectiveness due to toxicity and the emergence of drug resistance. Host-based antivirals are an alternative, but can cause nonspecific effects. Recent evidence shows that virus-infected cells can be selectively eliminated by targeting synthetic lethal (SL) partners of proteins disrupted by viral infection. Thus, we hypothesized that genes depleted in CRISPR knockout (KO) screens of virus-infected cells may be enriched in SL partners of proteins altered by infection. To investigate this, we established a computational pipeline predicting antiviral SL drug targets. First, we identified SARS-CoV-2-induced changes in gene products via a large compendium of omics data. Second, we identified SL partners for each altered gene product. Last, we screened CRISPR KO data for SL partners required for cell viability in infected cells. Despite differences in virus-induced alterations detected by various omics data, they share many predicted SL targets, with significant enrichment in CRISPR KO-depleted datasets. Our comparison of SARS-CoV-2 and influenza infection data revealed potential broad-spectrum, host-based antiviral SL targets. This suggests that CRISPR KO data are replete with common antiviral targets due to their SL relationship with virus-altered states and that such targets can be revealed from analysis of omics datasets and SL predictions.

Predicting host-based, synthetic lethal antiviral targets from omics data.

Traditional antiviral therapies often have limited effectiveness due to toxicity and development of drug resistance. Host-based antivirals, while an alternative, may lead to non-specific effects. Recent evidence shows that virus-infected cells can be selectively eliminated by targeting synthetic lethal (SL) partners of proteins disrupted by viral infection. Thus, we hypothesized that genes depleted in CRISPR KO screens of virus-infected cells may be enriched in SL partners of proteins altered by infection. To investigate this, we established a computational pipeline predicting SL drug targets of viral infections. First, we identified SARS-CoV-2-induced changes in gene products via a large compendium of omics data. Second, we identified SL partners for each altered gene product. Last, we screened CRISPR KO data for SL partners required for cell viability in infected cells. Despite differences in virus-induced alterations detected by various omics data, they share many predicted SL targets, with significant enrichment in CRISPR KO-depleted datasets. Comparing data from SARS-CoV-2 and influenza infections, we found possible broad-spectrum, host-based antiviral SL targets. This suggests that CRISPR KO data are replete with common antiviral targets due to their SL relationship with virus-altered states and that such targets can be revealed from analysis of omics datasets and SL predictions.

Molecular and Phylogenetic Characterization of Novel Papillomaviruses Isolated from Oral and Anogenital Neoplasms of Japanese Macaques (Macaca fuscata).

Papillomaviruses (PVs) are a diverse group of host species-specific DNA viruses, etiologically linked with various benign and malignant neoplasms of cutaneous and mucosal epithelia. Here, we describe the detection and characterization of the first two PVs naturally infecting Japanese macaques (Macaca fuscata), including the determination of their etiological association(s) with the development of original neoplasms. The molecular and phylogenetic analyses were performed on complete genome sequences of Macaca fuscata PV types 1 (MfuPV1) and 2 (MfuPV2), which were completely sequenced in samples of a malignant oral tumor and benign anogenital neoplasm of Japanese macaques, respectively. Subsequently, two type-specific quantitative real-time PCRs were developed to estimate viral loads of MfuPV1 and MfuPV2 and to evaluate their etiological roles. The in silico molecular analyses revealed that both viral genomes encode characteristic PV proteins with conserved functional domains and have a non-coding genomic region with regulatory sequences to regulate and complete the viral life cycle. However, additional experimental evidence is needed to finally confirm the presence and biological functionality of the molecular features of both novel PVs. While MfuPV1, together with PVs identified in other macaques, is classified into the Alphapapillomavirus (Alpha-PV) species 12, MfuPV2 is most likely a representative of the novel viral species within the Alpha-PV genus. Their relatively high viral loads suggest that both PVs are etiologically linked with the development of the original neoplasms.

Detection of 14 High-Risk Human Papillomaviruses Using Digital LAMP Assays on a Self-Digitization Chip.

Cervical cancer is the fourth-leading cause of cancer deaths among women worldwide and most cases occur in developing countries. Detection of high-risk (HR) HPV, the etiologic agent of cervical cancer, is a primary screening method for cervical cancer. However, the current gold standard for HPV detection, real-time PCR, is expensive, time-consuming, and instrumentation-intensive. A rapid, low-cost HPV detection method is needed for cervical cancer screening in low-resource settings. We previously developed a digital loop-mediated isothermal amplification (dLAMP) assay for rapid, quantitative detection of nucleic acids without the need for thermocycling. This assay employs a microfluidic self-digitization chip to automatically digitize a sample into an array of nanoliter wells in a simple assay format. Here we evaluate the dLAMP assay and self-digitization chip for detection of the commonly tested 14 high-risk HPVs in clinical samples. The dLAMP platform provided reliable genotyping and quantitative detection of the 14 high-risk HPVs with high sensitivity, demonstrating its potential for simple, rapid, and low-cost diagnosis of HPV infection.

Self-digitization chip for quantitative detection of human papillomavirus gene using digital LAMP.

Digital nucleic acid amplification and detection methods provide excellent sensitivity and specificity and allow absolute quantification of target nucleic acids. Isothermal methods such as digital loop-mediated isothermal amplification (digital LAMP) have potential for use in rapid disease diagnosis in low-resource settings due to their speed and lack of thermal cycling. We previously developed a self-digitization (SD) chip, a simple microfluidics device that automatically digitizes a sample into an array of nanoliter wells, for use in digital LAMP. In this work, we improve the SD chip design to increase sample loading efficiency, speed, and completeness, and test a range of well volumes and numbers. We demonstrate the diagnostic capability of this platform by applying it to quantifying human papillomavirus 18 gene.

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.

Base-By-Base Version 3: New Comparative Tools for Large Virus Genomes.

Base-By-Base is a comprehensive tool for the creation and editing of multiple sequence alignments that is coded in Java and runs on multiple platforms. It can be used with gene and protein sequences as well as with large viral genomes, which themselves can contain gene annotations. This report describes new features added to Base-By-Base over the last 7 years. The two most significant additions are: (1) The recoding and inclusion of "consensus-degenerate hybrid oligonucleotide primers" (CODEHOP), a popular tool for the design of degenerate primers from a multiple sequence alignment of proteins; and (2) the ability to perform fuzzy searches within the columns of sequence data in multiple sequence alignments to determine the distribution of sequence variants among the sequences. The intuitive interface focuses on the presentation of results in easily understood visualizations and providing the ability to annotate the sequences in a multiple alignment with analytic and user data.

Detection of novel Betapapillomaviruses and Gammapapillomaviruses in eyebrow hair follicles using a single-tube 'hanging droplet' PCR assay with modified pan-PV CODEHOP primers.

A modified pan-PV consensus-degenerate hybrid oligonucleotide primer (CODEHOP) PCR was developed for generic and sensitive detection of a broad-spectrum of human papillomaviruses (HPVs) infecting the cutaneous epithelium. To test the analytical sensitivity of the assay we examined 149 eyebrow hair follicle specimens from immunocompetent male patients. HPV DNA was detected in 60 % (89/149) of analysed eyebrow samples with a total of 48 different HPV sequences, representing 21 previously described HPVs and 27 putative novel HPV types. Evidence for ten novel HPV subtypes and seven viral variants, clustering to three out of five genera containing cutaneous HPVs, was also obtained. Thus, we have shown that the modified pan-PV CODEHOP PCR assay is able to identify multiple HPV types, even from different genera, in the same clinical sample. Overall, these results demonstrate that the pan-PV CODEHOP PCR is an excellent tool for screening and identification of novel cutaneous HPVs, even in samples with low viral loads.

Complete Unique Genome Sequence, Expression Profile, and Salivary Gland Tissue Tropism of the Herpesvirus 7 Homolog in Pigtailed Macaques.

UNLABELLED: Human herpesvirus 6A (HHV-6A), HHV-6B, and HHV-7 are classified as roseoloviruses and are highly prevalent in the human population. Roseolovirus reactivation in an immunocompromised host can cause severe pathologies. While the pathogenic potential of HHV-7 is unclear, it can reactivate HHV-6 from latency and thus contributes to severe pathological conditions associated with HHV-6. Because of the ubiquitous nature of roseoloviruses, their roles in such interactions and the resulting pathological consequences have been difficult to study. Furthermore, the lack of a relevant animal model for HHV-7 infection has hindered a better understanding of its contribution to roseolovirus-associated diseases. Using next-generation sequencing analysis, we characterized the unique genome of an uncultured novel pigtailed macaque roseolovirus. Detailed genomic analysis revealed the presence of gene homologs to all 84 known HHV-7 open reading frames. Phylogenetic analysis confirmed that the virus is a macaque homolog of HHV-7, which we have provisionally named Macaca nemestrina herpesvirus 7 (MneHV7). Using high-throughput RNA sequencing, we observed that the salivary gland tissue samples from nine different macaques had distinct MneHV7 gene expression patterns and that the overall number of viral transcripts correlated with viral loads in parotid gland tissue and saliva. Immunohistochemistry staining confirmed that, like HHV-7, MneHV7 exhibits a natural tropism for salivary gland ductal cells. We also observed staining for MneHV7 in peripheral nerve ganglia present in salivary gland tissues, suggesting that HHV-7 may also have a tropism for the peripheral nervous system. Our data demonstrate that MneHV7-infected macaques represent a relevant animal model that may help clarify the causality between roseolovirus reactivation and diseases. IMPORTANCE: Human herpesvirus 6A (HHV-6A), HHV-6B, and HHV-7 are classified as roseoloviruses. We have recently discovered that pigtailed macaques are naturally infected with viral homologs of HHV-6 and HHV-7, which we provisionally named MneHV6 and MneHV7, respectively. In this study, we confirm that MneHV7 is genetically and biologically similar to its human counterpart, HHV-7. We determined the complete unique MneHV7 genome sequence and provide a comprehensive annotation of all genes. We also characterized viral transcription profiles in salivary glands from naturally infected macaques. We show that broad transcriptional activity across most of the viral genome is associated with high viral loads in infected parotid glands and that late viral protein expression is detected in salivary duct cells and peripheral nerve ganglia. Our study provides new insights into the natural behavior of an extremely prevalent virus and establishes a basis for subsequent investigations of the mechanisms that cause HHV-7 reactivation and associated disease.

Discovery and biological characterization of two novel pig-tailed macaque homologs of HHV-6 and HHV-7.

Human herpesvirus-6 (HHV-6) and -7 (HHV-7) are Roseoloviruses within the Betaherpesvirus family, which have a high prevalence and suspected involvement in a number of diseases. Using CODEHOP-based PCR, we identified homologs of both viruses in saliva of pig-tailed macaques, provisionally named MneHV-6 and MneHV-7. This finding supports the existence of two distinct Roseolovirus lineages before the divergence of humans and macaques. Using specific qPCR assays, high levels of MneHV-6 and MneHV-7 DNA were detected in macaque saliva, although the frequency was greater for MneHV-7. A blood screen of 283 macaques revealed 10% MneHV-6 DNA positivity and 25% MneHV-7 positivity, with higher prevalences of MneHV-6 in older females and of MneHV-7 in younger males. Levels of MneHV-6 were increased in animals coinfected with MneHV-7, and both viruses were frequently detected in salivary gland and stomach tissues. Our discovery provides a unique animal model to answer unresolved questions regarding Roseolovirus pathology.

Oral squamous cell carcinoma in a pigtailed macaque (Macaca nemestrina).

An adult, gravid, female pigtailed macaque (Macaca nemestrina) presented for facial swelling centered on the left mandible that was approximately 5 cm wide. Differential diagnoses included infectious, inflammatory, and neoplastic origins. Definitive antemortem diagnosis was not possible, and the macaque's condition worsened despite supportive care. Necropsy findings included a mandibular mass that was locally invasive and expansile, encompassing approximately 80% of the left mandibular bone. The mass replaced portions of the soft palate, hard palate, sinuses, ear canal, and the caudal-rostral calvarium and masseter muscle. Histologically, the mass was a neoplasm that was poorly circumscribed, unencapsulated, and infiltrative invading regional bone and soft tissue. The mass consisted of polygonal squamous epithelial cells with intercellular bridging that breached the epithelial basement membrane and formed invasive nests, cords, and trabeculae. The mitotic rate averaged 3 per 400× field of view, with occasional bizarre mitotic figures. Epithelial cells often exhibited dyskeratosis, and the nests often contained compact lamellated keratin (keratin pearls). The neoplasm was positive via immunohistochemistry for pancytokeratin, variably positive for S100, and negative for vimentin, smooth muscle actin, and desmin. The gross, histologic, and immunohistochemical findings were consistent with an aggressive oral squamous cell carcinoma. The neoplasm was negative via PCR for papilloma virus. In general, neoplasia in macaques is rare. Although squamous cell carcinomas are one of the most common oral neoplasia in many species, to our knowledge this case represents the first reported oral squamous cell carcinoma in a pigtailed macaque.

CODEHOP PCR and CODEHOP PCR primer design.

While PCR primer design for the amplification of known sequences is usually quite straightforward, the design, and successful application of primers aimed at the detection of as yet unknown genes is often not. The search for genes that are presumed to be distantly related to a known gene sequence, such as homologous genes in different species, paralogs in the same genome, or novel pathogens in diverse hosts, often turns into the proverbial search for the needle in the haystack. PCR-based methods commonly used to address this issue involve the use of either consensus primers or degenerate primers, both of which have significant shortcomings regarding sensitivity and specificity. We have developed a novel primer design approach that diminishes these shortcomings and instead takes advantage of the strengths of both consensus and degenerate primer designs, by combining the two concepts into a Consensus-Degenerate Hybrid Oligonucleotide Primer (CODEHOP) approach. CODEHOP PCR primers contain a relatively short degenerate 3' core and a 5' nondegenerate clamp. The 3' degenerate core consists of a pool of primers containing all possible codons for a 3-4 aminoacid motif that is highly conserved in multiply aligned sequences from known members of a protein family. Each primer in the pool also contains a single 5' nondegenerate nucleotide sequence derived from a codon consensus across the aligned aminoacid sequences flanking the conserved motif. During the initial PCR amplification cycles, the degenerate core is responsible for specific binding to sequences encoding the conserved aminoacid motif. The longer consensus clamp region serves to stabilize the primer and allows the participation of all primers in the pool in the efficient amplification of products during later PCR cycles. We have developed an interactive web site and algorithm (iCODEHOP) for designing CODEHOP PCR primers from multiply aligned protein sequences, which is freely available online. Here, we describe the workflow of a typical CODEHOP PCR assay design and optimization and give a specific implementation example along with "best-practice" advice.

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.

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.

Distinct roles for p53, p27Kip1, and p21Cip1 during tumor development.

Mutations in p53 and reduced expression of the Cdk inhibitor p27Kip1 are frequently observed in a wide variety of human cancers, but it is not known whether alterations in these tumor suppressors interact to influence tumor progression. To address this, we measured tumor latency and spectrum in p53 and p27 single and compound mutant mice. p53-/- (null) mice developed T-cell lymphomas and soft-tissue sarcomas, while p27-/- mice developed adenomas of the pituitary and lung, but with much longer latency. The latency for tumor development in p53-/- p27-/- and p53-/- p27+/- compound mutant mice was significantly reduced, by 15-30%, compared to single mutant p53-/- mice. The tumor spectrum in the compound mutants was similar to that of p53-/- mice, and additional tumors of diverse histotypes. In tumors from p53-/- mice, p27 protein levels were reduced to a greater extent than were mRNA levels, indicating that p27 is downregulated in tumors at the transcriptional as well as post-transcriptional levels. In contrast, mice deficient in another Cdk inhibitor p21Cip1, which is also a transcriptional target and effector of p53, showed only a marginal increase in tumor predisposition in response to ENU treatment. Thus, downregulation of p27 is a common feature in p53-/- tumors. Germline deletion of one or both alleles of p27 accelerates tumor development and associated mortality in p53-/- mice, indicating potent synergy between loss of p27 and p53. Although p21 is functionally similar to both p53 and p27, it plays a lesser role in tumor suppression. These results further highlight the highly cooperative nature of p27 and its central role in tumor suppression.

Pathway-specific tumor suppression. Reduction of p27 accelerates gastrointestinal tumorigenesis in Apc mutant mice, but not in Smad3 mutant mice.

Expression of the cyclin-dependent kinase inhibitor p27(Kip1) (p27) is frequently reduced in human colorectal cancer, and this correlates with poor patient prognosis. To clarify the role of p27 in gastrointestinal (GI) cancer, we measured p27 expression, as well as the effect of germline deletion of p27, in 3 different mouse models of GI neoplasia. p27 expression was frequently reduced in GI tumors arising in 1,2-dimethylhydrazine (DMH) treated mice, and in Apc mutant Min/+ mice, but not in GI tumors arising in Smad3 mutant mice. Germline deletion of p27 resulted in accelerated tumor development and increased tumor cell proliferation in both DMH treated and Min/+ mice, but not in Smad3 mutant mice. p27 deficiency also led to increased adenoma to adenocarcinoma progression. These results indicate that reduction of p27 cooperates with mutations in Apc but not in Smad3 during GI tumorigenesis. Thus, tumor suppression by p27 is contingent on the specific oncogenic pathway that drives tumor development.