The model student: GPT-4 performance on graduate biomedical science exams

The GPT-4 large language model (LLM) and ChatGPT chatbot have emerged as accessible and capable tools for generating English-language text in a variety of formats. GPT-4 has previously performed well when applied to questions from multiple standardized examinations. However, further evaluation of trustworthiness and accuracy of GPT-4 responses across various knowledge domains is essential before its use as a reference resource. Here, we assess GPT-4 performance on nine graduate-level examinations in the biomedical sciences (seven blinded), finding that GPT-4 scores exceed the student average in seven of nine cases and exceed all student scores for four exams. GPT-4 performed very well on fill-in-the-blank, short-answer, and essay questions, and correctly answered several questions on figures sourced from published manuscripts. Conversely, GPT-4 performed poorly on questions with figures containing simulated data and those requiring a hand-drawn answer. Two GPT-4 answer-sets were flagged as plagiarism based on answer similarity and some model responses included detailed hallucinations. In addition to assessing GPT-4 performance, we discuss patterns and limitations in GPT-4 capabilities with the goal of informing design of future academic examinations in the chatbot era.

Circulating microRNA Biomarkers of Thiazide Response in Hypertension

BACKGROUND

Thiazide diuretics are the second most frequently prescribed class of antihypertensives, but up to 50% of patients with hypertension have minimal antihypertensive response to thiazides. We explored circulating microRNAs (miRNAs) in search of predictive biomarkers of thiazide response.

METHODS AND RESULTS

We profiled 754 miRNAs in baseline plasma samples of 36 hypertensive European American adults treated with hydrochlorothiazide, categorized into responders (n=18) and nonresponders (n=18) on the basis of diastolic blood pressure response to hydrochlorothiazide. miRNAs with ≥2.5-fold differential expression between responders and nonresponders were considered for validation in 3 cohorts (n=50 each): hydrochlorothiazide-treated European Americans, chlorthalidone-treated European Americans, and hydrochlorothiazide-treated Black individuals. Different blood pressure phenotypes including categorical (responder versus nonresponder) and continuous diastolic blood pressure and systolic blood pressure were tested for association with the candidate miRNA expression using multivariate regression analyses adjusting for age, sex, and baseline blood pressure. After quality control, 74 miRNAs were available for screening, 19 of which were considered for validation. In the validation cohort, miR-193b-3p and 30d-5p showed significant associations with continuous SBP or diastolic blood pressure response or both, to hydrochlorothiazide in European Americans at Benjamini-Hochberg adjusted P<0.05. In the combined analysis of validation cohorts, let-7g (odds ratio, 0.6 [95% CI, 0.4-0.8]), miR-142-3p (odds ratio, 1.1 [95% CI, 1.0, 1.2]), and miR-423-5p (odds ratio, 0.7 [95% CI, 0.5-0.9]) associated with categorical diastolic blood pressure response at Benjamini-Hochberg adjusted P<0.05. Predicted target genes of the 5 miRNAs were mapped to key hypertension pathways: lysine degradation, fatty acid biosynthesis, and metabolism.

CONCLUSIONS

The above identified circulating miRNAs may have a potential for clinical use as biomarkers for thiazide diuretic selection in hypertension.

REGISTRATION

URL: ClinicalTrials.gov. Unique identifiers: NCT00246519, NCT01203852, NCT00005520.

Identifying locally actionable strategies to increase participant acceptability and feasibility to participate in Phase I cancer clinical trials

BACKGROUND

Recruitment of cancer clinical trial (CCT) participants, especially participants representing the diversity of the US population, is necessary to create successful medications and a continual challenge. These challenges are amplified in Phase I cancer trials that focus on evaluating the safety of new treatments and are the gateway to treatment development. In preparation for recruitment to a Phase I recurrent head and neck cancer (HNC) trial, we assessed perceived barriers to participation or referral and suggestions for recruitment among people with HNC and community physicians (oncologist, otolaryngologist or surgeon).

METHODS

Between December 2020 and February 2022, we conducted a qualitative needs assessment via semistructured interviews with a race and ethnicity-stratified sample of people with HNC (n = 30: 12 non-Hispanic White, 9 non-Hispanic African American, 8 Hispanic and 1 non-Hispanic Pacific Islander) and community physicians (n = 16) within the University of Florida Health Cancer Center catchment area. Interviews were analyzed using a qualitative content analysis approach to describe perspectives and identify relevant themes.

RESULTS

People with HNC reported thematic barriers included: concerns about side effects, safety and efficacy; lack of knowledge and systemic and environmental obstacles. Physicians identified thematic barriers of limited physician knowledge; clinic and physician barriers and structural barriers. People with HNC and physicians recommended themes included: improved patient education, dissemination of trial information and interpersonal communication between community physicians and CCT staff.

CONCLUSIONS

The themes identified by people with HNC and community physicians are consistent with research efforts and recommendations on how to increase the participation of people from minoritized populations in CCTs. This community needs assessment provides direction on the selection of strategies to increase CCT participation and referral.

PATIENT OR PUBLIC CONTRIBUTION

This study focused on people with HNC and community physicians' lived experience and their interpretations of how they would consider a future Phase I clinical trial. In addition to our qualitative data reflecting community voices, a community member reviewed the draft interview guide before data collection and both people with HNC and physicians aided interpretation of the findings.

Hybkit: a Python API and command-line toolkit for hybrid sequence data from chimeric RNA methods

SUMMARY

Experimental methods using microRNA/target ligation have recently provided significant insights into microRNA functioning through generation of chimeric (hybrid) RNA sequences. Here, we introduce Hybkit, a Python3 API, and command-line toolkit for analysis of hybrid sequence data in the "hyb" file format to enable customizable evaluation and annotation of hybrid characteristics. The Hybkit API includes a suite of python objects for developing custom analyses of hybrid data as well as miRNA-specific analysis methods, built-in plotting of analysis results, and incorporation of predicted miRNA/target interactions in Vienna format.

AVAILABILITY AND IMPLEMENTATION

Hybkit is provided free and open source under the GNU GPL license at github.com/RenneLab/hybkit and archived on Zenodo (doi.org/10.5281/zenodo.7834299). Hybkit distributions are also provided via PyPI (pypi.org/project/hybkit), Conda (bioconda.github.io/recipes/hybkit/README.html), and Docker (quay.io/repository/biocontainers/hybkit).

Sex-Biased Expression Is Associated With Chromatin State in Drosophila melanogaster and Drosophila simulans

In Drosophila melanogaster and D. simulans head tissue, 60% of orthologous genes show evidence of sex-biased expression in at least one species. Of these, ∼39% (2,192) are conserved in direction. We hypothesize enrichment of open chromatin in the sex where we see expression bias and closed chromatin in the opposite sex. Male-biased orthologs are significantly enriched for H3K4me3 marks in males of both species (∼89% of male-biased orthologs vs. ∼76% of unbiased orthologs). Similarly, female-biased orthologs are significantly enriched for H3K4me3 marks in females of both species (∼90% of female-biased orthologs vs. ∼73% of unbiased orthologs). The sex-bias ratio in female-biased orthologs was similar in magnitude between the two species, regardless of the closed chromatin (H3K27me2me3) marks in males. However, in male-biased orthologs, the presence of H3K27me2me3 in both species significantly reduced the correlation between D. melanogaster sex-bias ratio and the D. simulans sex-bias ratio. Male-biased orthologs are enriched for evidence of positive selection in the D. melanogaster group. There are more male-biased genes than female-biased genes in both species. For orthologs with gains/losses of sex-bias between the two species, there is an excess of male-bias compared to female-bias, but there is no consistent pattern in the relationship between H3K4me3 or H3K27me2me3 chromatin marks and expression. These data suggest chromatin state is a component of the maintenance of sex-biased expression and divergence of sex-bias between species is reflected in the complexity of the chromatin status.

Virology under the Microscope-a Call for Rational Discourse

Viruses have brought humanity many challenges: respiratory infection, cancer, neurological impairment and immunosuppression to name a few. Virology research over the last 60+ years has responded to reduce this disease burden with vaccines and antivirals. Despite this long history, the COVID-19 pandemic has brought unprecedented attention to the field of virology. Some of this attention is focused on concern about the safe conduct of research with human pathogens. A small but vocal group of individuals has seized upon these concerns - conflating legitimate questions about safely conducting virus-related research with uncertainties over the origins of SARS-CoV-2. The result has fueled public confusion and, in many instances, ill-informed condemnation of virology. With this article, we seek to promote a return to rational discourse. We explain the use of gain-of-function approaches in science, discuss the possible origins of SARS-CoV-2 and outline current regulatory structures that provide oversight for virological research in the United States. By offering our expertise, we - a broad group of working virologists - seek to aid policy makers in navigating these controversial issues. Balanced, evidence-based discourse is essential to addressing public concern while maintaining and expanding much-needed research in virology.

Investigating miR-9 as a mediator in laryngeal cancer health disparities

Background

For several decades, Black patients have carried a higher burden of laryngeal cancer among all races. Even when accounting for sociodemographics, a disparity remains. Differentially expressed microRNAs have been linked to racially disparate clinical outcomes in breast and prostate cancers, yet an association in laryngeal cancer has not been addressed. In this study, we present our computational analysis of differentially expressed miRNAs in Black compared with White laryngeal cancer and further validate microRNA-9-5p (miR-9-5p) as a potential mediator of cancer phenotype and chemoresistance.

Methods

Bioinformatic analysis of 111 (92 Whites, 19 Black) laryngeal squamous cell carcinoma (LSCC) specimens from the TCGA revealed miRNAs were significantly differentially expressed in Black compared with White LSCC. We focused on miR-9-5 p which had a significant 4-fold lower expression in Black compared with White LSCC (p<0.05). After transient transfection with either miR-9 mimic or inhibitor in cell lines derived from Black (UM-SCC-12) or White LSCC patients (UM-SCC-10A), cellular migration and cell proliferation was assessed. Alterations in cisplatin sensitivity was evaluated in transient transfected cells via IC50 analysis. qPCR was performed on transfected cells to evaluate miR-9 targets and chemoresistance predictors, ABCC1 and MAP1B.

Results

Northern blot analysis revealed mature miR-9-5p was inherently lower in cell line UM-SCC-12 compared with UM-SCC-10A. UM -SCC-12 had baseline increase in cellular migration (p < 0.01), proliferation (p < 0.0001) and chemosensitivity (p < 0.01) compared to UM-SCC-10A. Increasing miR-9 in UM-SCC-12 cells resulted in decreased cellular migration (p < 0.05), decreased proliferation (p < 0.0001) and increased sensitivity to cisplatin (p < 0.001). Reducing miR-9 in UM-SCC-10A cells resulted in increased cellular migration (p < 0.05), increased proliferation (p < 0.05) and decreased sensitivity to cisplatin (p < 0.01). A significant inverse relationship in ABCC1 and MAP1B gene expression was observed when miR-9 levels were transiently elevated or reduced in either UM-SCC-12 or UM-SCC-10A cell lines, respectively, suggesting modulation by miR-9.

Conclusion

Collectively, these studies introduce differential miRNA expression in LSCC cancer health disparities and propose a role for low miR-9-5p as a mediator in LSCC tumorigenesis and chemoresistance.

Virology under the Microscope-a Call for Rational Discourse

Viruses have brought humanity many challenges: respiratory infection, cancer, neurological impairment and immunosuppression to name a few. Virology research over the last 60+ years has responded to reduce this disease burden with vaccines and antivirals. Despite this long history, the COVID-19 pandemic has brought unprecedented attention to the field of virology. Some of this attention is focused on concern about the safe conduct of research with human pathogens. A small but vocal group of individuals has seized upon these concerns - conflating legitimate questions about safely conducting virus-related research with uncertainties over the origins of SARS-CoV-2. The result has fueled public confusion and, in many instances, ill-informed condemnation of virology. With this article, we seek to promote a return to rational discourse. We explain the use of gain-of-function approaches in science, discuss the possible origins of SARS-CoV-2 and outline current regulatory structures that provide oversight for virological research in the United States. By offering our expertise, we - a broad group of working virologists - seek to aid policy makers in navigating these controversial issues. Balanced, evidence-based discourse is essential to addressing public concern while maintaining and expanding much-needed research in virology.

Virology under the Microscope-a Call for Rational Discourse

Viruses have brought humanity many challenges: respiratory infection, cancer, neurological impairment and immunosuppression to name a few. Virology research over the last 60+ years has responded to reduce this disease burden with vaccines and antivirals. Despite this long history, the COVID-19 pandemic has brought unprecedented attention to the field of virology. Some of this attention is focused on concern about the safe conduct of research with human pathogens. A small but vocal group of individuals has seized upon these concerns - conflating legitimate questions about safely conducting virus-related research with uncertainties over the origins of SARS-CoV-2. The result has fueled public confusion and, in many instances, ill-informed condemnation of virology. With this article, we seek to promote a return to rational discourse. We explain the use of gain-of-function approaches in science, discuss the possible origins of SARS-CoV-2 and outline current regulatory structures that provide oversight for virological research in the United States. By offering our expertise, we - a broad group of working virologists - seek to aid policy makers in navigating these controversial issues. Balanced, evidence-based discourse is essential to addressing public concern while maintaining and expanding much-needed research in virology.

Sex-biased expression is associated with chromatin state in D. melanogaster and D. simulans

We propose a new model for the association of chromatin state and sex-bias in expression. We hypothesize enrichment of open chromatin in the sex where we see expression bias (OS) and closed chromatin in the opposite sex (CO). In this study of D. melanogaster and D. simulans head tissue, sex-bias in expression is associated with H3K4me3 (open mark) in males for male-biased genes and in females for female-biased genes in both species. Sex-bias in expression is also largely conserved in direction and magnitude between the two species on the X and autosomes. In male-biased orthologs, the sex-bias ratio is more divergent between species if both species have H3K27me2me3 marks in females compared to when either or neither species has H3K27me2me3 in females. H3K27me2me3 marks in females are associated with male-bias in expression on the autosomes in both species, but on the X only in D. melanogaster . In female-biased orthologs the relationship between the species for the sex-bias ratio is similar regardless of the H3K27me2me3 marks in males. Female-biased orthologs are more similar in the ratio of sex-bias than male-biased orthologs and there is an excess of male-bias in expression in orthologs that gain/lose sex-bias. There is an excess of male-bias in sex-limited expression in both species suggesting excess male-bias is due to rapid evolution between the species. The X chromosome has an enrichment in male-limited H3K4me3 in both species and an enrichment of sex-bias in expression compared to the autosomes.

Metabolic signatures associated with oncolytic myxoma viral infections

Oncolytic viral therapy is a recent advance in cancer treatment, demonstrating promise as a primary treatment option. To date, the secondary metabolic effects of viral infection in cancer cells has not been extensively studied. In this work, we have analyzed early-stage metabolic changes in cancer cells associated with oncolytic myxoma virus infection. Using GC–MS based metabolomics, we characterized the myxoma virus infection induced metabolic changes in three cancer cell lines—small cell (H446) and non-small cell (A549) lung cancers, and glioblastoma (SFxL). We show that even at an early stage (6 and 12 h) myxoma infection causes profound changes in cancer cell metabolism spanning several important pathways such as the citric acid cycle, fatty acid metabolism, and amino acid metabolism. In general, the metabolic effects of viral infection across cell lines are not conserved. However, we have identified several candidate metabolites that can potentially serve as biomarkers for monitoring oncolytic viral action in general.

DExD/H Box Helicases DDX24 and DDX49 Inhibit Reactivation of Kaposi’s Sarcoma Associated Herpesvirus by Interacting with Viral mRNAs

Kaposi’s sarcoma-associated herpesvirus (KSHV) is an oncogenic gammaherpesvirus that is the causative agent of primary effusion lymphoma and Kaposi’s sarcoma. In healthy carriers, KSHV remains latent, but a compromised immune system can lead to lytic viral replication that increases the probability of tumorigenesis. RIG-I-like receptors (RLRs) are members of the DExD/H box helicase family of RNA binding proteins that recognize KSHV to stimulate the immune system and prevent reactivation from latency. To determine if other DExD/H box helicases can affect KSHV lytic reactivation, we performed a knock-down screen that revealed DHX29-dependent activities appear to support viral replication but, in contrast, DDX24 and DDX49 have antiviral activity. When DDX24 or DDX49 are overexpressed in BCBL-1 cells, transcription of all lytic viral genes and genome replication were significantly reduced. RNA immunoprecipitation of tagged DDX24 and DDX49 followed by next-generation sequencing revealed that the helicases bind to mostly immediate-early and early KSHV mRNAs. Transfection of expression plasmids of candidate KSHV transcripts, identified from RNA pull-down, demonstrated that KSHV mRNAs stimulate type I interferon (alpha/beta) production and affect the expression of multiple interferon-stimulated genes. Our findings reveal that host DExD/H box helicases DDX24 and DDX49 recognize gammaherpesvirus transcripts and convey an antiviral effect in the context of lytic reactivation.

Sequencing of Argonaute-bound microRNA/mRNA hybrids reveals regulation of the unfolded protein response by microRNA-320a

MicroRNAs (miRNA) are short non-coding RNAs widely implicated in gene regulation. Most metazoan miRNAs utilize the RNase III enzymes Drosha and Dicer for biogenesis. One notable exception is the RNA polymerase II transcription start sites (TSS) miRNAs whose biogenesis does not require Drosha. The functional importance of the TSS-miRNA biogenesis is uncertain. To better understand the function of TSS-miRNAs, we applied a modified Crosslinking, Ligation, and Sequencing of Hybrids on Argonaute (AGO-qCLASH) to identify the targets for TSS-miRNAs in HCT116 colorectal cancer cells with or without DROSHA knockout. We observed that miR-320a hybrids dominate in TSS-miRNA hybrids identified by AGO-qCLASH. Targets for miR-320a are enriched for the eIF2 signaling pathway, a downstream component of the unfolded protein response. Consistently, in miR-320a mimic- and antagomir- transfected cells, differentially expressed gene products are associated with eIF2 signaling. Within the AGO-qCLASH data, we identified the endoplasmic reticulum (ER) chaperone calnexin as a direct miR-320a down-regulated target, thus connecting miR-320a to the unfolded protein response. During ER stress, but not amino acid deprivation, miR-320a up-regulates ATF4, a critical transcription factor for resolving ER stress. In summary, our study investigates the targetome of the TSS-miRNAs in colorectal cancer cells and establishes miR-320a as a regulator of unfolded protein response.

Modified Cross-Linking, Ligation, and Sequencing of Hybrids (qCLASH) to Identify MicroRNA Targets

This protocol was designed to identify microRNA (miRNA) targetomes from smaller‐input samples by performing a simplified workflow of the Cross‐Linking and Sequencing of Hybrids (CLASH) technique developed in the Tollervey group. In this ribonomics‐based technique, Cross‐Linking and Immunoprecipitation (CLIP) of Argonaute (Ago) is combined with an RNA ligase reaction that yields covalently bound “hybrids” between miRNAs and their target RNAs. While this iteration of CLIP identifies “high‐confidence” or “unambiguous” miRNA targets, the added ligation step is highly inefficient and therefore requires large numbers of cultured cells. To make this powerful approach applicable to smaller cell numbers, we created qCLASH, incorporating a workflow that performs all enzymatic reactions on bead‐bound complexes and omits gel purification of immunoprecipitated Ago complexes associated with major loss of RNA. At a sequencing depth of 100 million reads per library, which is highly feasible with rapidly decreasing sequencing costs, qCLASH, when used with three biological replicates, results in thousands of high‐confidence miRNA targets. qCLASH was first developed to identify viral miRNA targetomes of endothelial cells infected with Kaposi's sarcoma−associated herpesvirus. Since then, qCLASH has been applied to Epstein‐Barr virus− and MHV68‐infected cells, and more recently to metastatic melanoma and breast cancer cells. Currently, protocols are under development to apply qCLASH to human solid tumor specimens. © 2021 The Authors. Current Protocols published by Wiley Periodicals LLC.

Basic Protocol: Quick Cross‐Linking and Sequencing of Hybrids (qCLASH)

Support Protocol: Optimization of Ago immunoprecipitation

A noncanonical microRNA derived from the snaR-A noncoding RNA targets a metastasis inhibitor

MicroRNAs (miRNAs) are small noncoding RNAs that function as critical posttranscriptional regulators in various biological processes. While most miRNAs are generated from processing of long primary transcripts via sequential Drosha and Dicer cleavage, some miRNAs that bypass Drosha cleavage can be transcribed as part of another small noncoding RNA. Here, we develop the target-oriented miRNA discovery (TOMiD) bioinformatic analysis method to identify Drosha-independent miRNAs from Argonaute crosslinking and sequencing of hybrids (Ago-CLASH) data sets. Using this technique, we discovered a novel miRNA derived from a primate specific noncoding RNA, the small NF90 associated RNA A (snaR-A). The miRNA derived from snaR-A (miR-snaR) arises independently of Drosha processing but requires Exportin-5 and Dicer for biogenesis. We identify that miR-snaR is concurrently up-regulated with the full snaR-A transcript in cancer cells. Functionally, miR-snaR associates with Ago proteins and targets NME1, a key metastasis inhibitor, contributing to snaR-A's role in promoting cancer cell migration. Our findings suggest a functional link between a novel miRNA and its precursor noncoding RNA.

EBV miRNAs are potent effectors of tumor cell transcriptome remodeling in promoting immune escape

The Epstein Barr virus (EBV) contributes to the tumor phenotype through a limited set of primarily non-coding viral RNAs, including 31 mature miRNAs. Here we investigated the impact of EBV miRNAs on remodeling the tumor cell transcriptome. Strikingly, EBV miRNAs displayed exceptionally abundant expression in primary EBV-associated Burkitt’s Lymphomas (BLs) and Gastric Carcinomas (GCs). To investigate viral miRNA targeting, we used the high-resolution approach, CLASH in GC and BL cell models. Affinity constant calculations of targeting efficacies for CLASH hits showed that viral miRNAs bind their targets more effectively than their host counterparts, as did Kaposi’s sarcoma-associated herpesvirus (KSHV) and murine gammaherpesvirus 68 (MHV68) miRNAs. Using public BL and GC RNA-seq datasets, we found that high EBV miRNA targeting efficacies translates to enhanced reduction of target expression. Pathway analysis of high efficacy EBV miRNA targets showed enrichment for innate and adaptive immune responses. Inhibition of the immune response by EBV miRNAs was functionally validated in vivo through the finding of inverse correlations between EBV miRNAs and immune cell infiltration and T-cell diversity in BL and GC datasets. Together, this study demonstrates that EBV miRNAs are potent effectors of the tumor transcriptome that play a role in suppressing host immune response.

Burkitt’s Lymphoma and gastric cancer are both associated with EBV, a prolific DNA tumor virus that latently resides in nearly all human beings. Despite mostly restricting viral gene expression to noncoding RNAs, EBV has important influences on the fitness of infected tumor cells. Here, we show that the miRNA class of viral noncoding RNAs are a major viral contributor to remodeling the tumor cell regulatory machinery in patient tumor samples. First, an assessment of miRNA expression in clinical tumor samples showed that EBV miRNAs are expressed at unexpectedly high levels relative to cell miRNAs. Using a highly specific miRNA target identification approach, CLASH, we comprehensively identified both viral and cellular miRNA targets and the relative abundance of each miRNA-mRNA interaction. We also show that viral miRNAs bind to and alter the expression of their mRNA targets more effectively than their cellular miRNA counterparts. Pathway analysis of the most effectively targeted mRNAs revealed enrichment of immune signaling pathways and we show a corresponding inverse correlation between EBV miRNA expression and infiltrating immune cells in EBV positive primary tumors. Altogether, this study shows that EBV miRNAs are key regulators of the tumor cell phenotype and the immune cell microenvironment.

EBV miRNAs are potent effectors of tumor cell transcriptome remodeling in promoting immune escape

The Epstein Barr virus (EBV) contributes to the tumor phenotype through a limited set of primarily non-coding viral RNAs, including 31 mature miRNAs. Here we investigated the impact of EBV miRNAs on remodeling the tumor cell transcriptome. Strikingly, EBV miRNAs displayed exceptionally abundant expression in primary EBV-associated Burkitt's Lymphomas (BLs) and Gastric Carcinomas (GCs). To investigate viral miRNA targeting, we used the high-resolution approach, CLASH in GC and BL cell models. Affinity constant calculations of targeting efficacies for CLASH hits showed that viral miRNAs bind their targets more effectively than their host counterparts, as did Kaposi's sarcoma-associated herpesvirus (KSHV) and murine gammaherpesvirus 68 (MHV68) miRNAs. Using public BL and GC RNA-seq datasets, we found that high EBV miRNA targeting efficacies translates to enhanced reduction of target expression. Pathway analysis of high efficacy EBV miRNA targets showed enrichment for innate and adaptive immune responses. Inhibition of the immune response by EBV miRNAs was functionally validated in vivo through the finding of inverse correlations between EBV miRNAs and immune cell infiltration and T-cell diversity in BL and GC datasets. Together, this study demonstrates that EBV miRNAs are potent effectors of the tumor transcriptome that play a role in suppressing host immune response.

Cross-Linking Ligation and Sequencing of Hybrids (qCLASH) Reveals an Unpredicted miRNA Targetome in Melanoma Cells

MicroRNAs are key post-transcriptional gene regulators often displaying aberrant expression patterns in cancer. As microRNAs are promising disease-associated biomarkers and modulators of responsiveness to anti-cancer therapies, a solid understanding of their targetome is crucial. Despite enormous research efforts, the success rates of available tools to reliably predict microRNAs (miRNA)-target interactions remains limited. To investigate the disease-associated miRNA targetome, we have applied modified cross-linking ligation and sequencing of hybrids (qCLASH) to BRAF-mutant melanoma cells. The resulting RNA-RNA hybrid molecules provide a comprehensive and unbiased snapshot of direct miRNA-target interactions. The regulatory effects on selected miRNA target genes in predicted vs. non-predicted binding regions was validated by miRNA mimic experiments. Most miRNA-target interactions deviate from the central dogma of miRNA targeting up to 60% interactions occur via non-canonical seed pairing with a strong contribution of the 3' miRNA sequence, and over 50% display a clear bias towards the coding sequence of mRNAs. miRNAs targeting the coding sequence can directly reduce gene expression (miR-34a/CD68), while the majority of non-canonical miRNA interactions appear to have roles beyond target gene suppression (miR-100/AXL). Additionally, non-mRNA targets of miRNAs (lncRNAs) whose interactions mainly occur via non-canonical binding were identified in melanoma. This first application of CLASH sequencing to cancer cells identified over 8 K distinct miRNA-target interactions in melanoma cells. Our data highlight the importance non-canonical interactions, revealing further layers of complexity of post-transcriptional gene regulation in melanoma, thus expanding the pool of miRNA-target interactions, which have so far been omitted in the cancer field.

Age-Related Changes in miRNA Expression Influence GSTZ1 and Other Drug Metabolizing Enzymes

Previous work has shown that hepatic levels of human glutathione transferase zeta 1 (GSTZ1) protein, involved in tyrosine catabolism and responsible for metabolism of the investigational drug dichloroacetate, increase in cytosol after birth before reaching a plateau around age 7. However, the mechanism regulating this change of expression is still unknown, and previous studies showed that GSTZ1 mRNA levels did not correlate with GSTZ1 protein expression. In this study, we addressed the hypothesis that microRNAs (miRNAs) could regulate expression of GSTZ1. We obtained liver samples from donors aged less than 1 year or older than 13 years and isolated total RNA for use in a microarray to identify miRNAs that were downregulated in the livers of adults compared with children. From a total of 2578 human miRNAs tested, 63 miRNAs were more than 2-fold down-regulated in adults, of which miR-376c-3p was predicted to bind to the 3' untranslated region of GSTZ1 mRNA. There was an inverse correlation of miR-376c-3p and GSTZ1 protein expression in the liver samples. Using cell culture, we confirmed that miR-376c-3p could downregulate GSTZ1 protein expression. Our findings suggest that miR-376c-3p prevents production of GSTZ1 through inhibition of translation. These experiments further our understanding of GSTZ1 regulation. Furthermore, our array results provide a database resource for future studies on mechanisms regulating human hepatic developmental expression. SIGNIFICANCE STATEMENT: Hepatic glutathione transferase zeta 1 (GSTZ1) is responsible for metabolism of the tyrosine catabolite maleylacetoacetate as well as the investigational drug dichloroacetate. Through examination of microRNA (miRNA) expression in liver from infants and adults and studies in cells, we showed that expression of GSTZ1 is controlled by miRNA. This finding has application to the dosing regimen of the drug dichloroacetate. The miRNA expression profiles are provided and will prove useful for future studies of drug-metabolizing enzymes in infants and adults.

Human Cytomegalovirus Latency and Myelosuppression: A microRNA-Dependent Yin and Yang Regulatory Loop

In this issue of Cell Host & Microbe, Hancock et al., 2020 show that latent HCMV infection, and specifically miR-US5-2, induces TGF-β secretion, which inhibits myelopoiesis in uninfected HPCs. They also show that HCMV-infected cells become resistant to TGF-β signaling through targeting of SMAD3 by miR-UL22-A-3p and -5p.

Downregulation of the human peripheral myelin protein 22 gene by miR-29a in cellular models of Charcot-Marie-Tooth disease

The majority of hereditary neuropathies are caused by duplication of the peripheral myelin protein 22 (PMP22) gene. Therefore, mechanisms to suppress the expression of the PMP22 gene have high therapeutic significance. Here we asked whether the human PMP22 gene is a target for regulation by microRNA 29a (miR-29a). Using bioinformatics, we determined that the human PMP22 gene contains the conserved seed sequence of the miR-29a binding site and this regulatory motif is included in the duplicated region in neuropathic patients. Using luciferase reporter assays in HEK293 cells, we demonstrated that transient transfection of a miR-29a mimic is associated with reduction in PMP22 3'UTR reporter activity. Transfecting normal and humanized transgenic neuropathic mouse Schwann cells with a miR-29a expression plasmid effectively lowered both the endogenous mouse and the transgenic human PMP22 transcripts compared with control vector. In dermal fibroblasts derived from neuropathic patients, ectopic expression of miR-29a led to ~50% reduction in PMP22 mRNA, which corresponded to ~20% decrease in PMP22 protein levels. Significantly, miR-29a-mediated reduction in PMP22 mitigated the reduced mitotic capacity of the neuropathic cells. Together, these results support further testing of miR-29a and/or PMP22-targeting siRNAs as therapeutic agents for correcting the aberrant expression of PMP22 in neuropathic patients.

Computational analysis of ribonomics datasets identifies long non-coding RNA targets of γ-herpesviral miRNAs

Ribonomics experiments involving crosslinking and immuno-precipitation (CLIP) of Ago proteins have expanded the understanding of the miRNA targetome of several organisms. These techniques, collectively referred to as CLIP-seq, have been applied to identifying the mRNA targets of miRNAs expressed by Kaposi’s Sarcoma-associated herpes virus (KSHV) and Epstein–Barr virus (EBV). However, these studies focused on identifying only those RNA targets of KSHV and EBV miRNAs that are known to encode proteins. Recent studies have demonstrated that long non-coding RNAs (lncRNAs) are also targeted by miRNAs. In this study, we performed a systematic re-analysis of published datasets from KSHV- and EBV-driven cancers. We used CLIP-seq data from lymphoma cells or EBV-transformed B cells, and a crosslinking, ligation and sequencing of hybrids dataset from KSHV-infected endothelial cells, to identify novel lncRNA targets of viral miRNAs. Here, we catalog the lncRNA targetome of KSHV and EBV miRNAs, and provide a detailed in silico analysis of lncRNA–miRNA binding interactions. Viral miRNAs target several hundred lncRNAs, including a subset previously shown to be aberrantly expressed in human malignancies. In addition, we identified thousands of lncRNAs to be putative targets of human miRNAs, suggesting that miRNA–lncRNA interactions broadly contribute to the regulation of gene expression.

Identification of murine gammaherpesvirus 68 miRNA-mRNA hybrids reveals miRNA target conservation among gammaherpesviruses including host translation and protein modification machinery

Gammaherpesviruses, including the human pathogens Epstein-Barr virus (EBV) and Kaposi’s sarcoma-associated herpesvirus (KSHV), establish lifelong latent infection in B cells and are associated with a variety of tumors. In addition to protein coding genes, these viruses encode numerous microRNAs (miRNAs) within their genomes. While putative host targets of EBV and KSHV miRNAs have been previously identified, the specific functions of these miRNAs during in vivo infection are largely unknown. Murine gammaherpesvirus 68 (MHV68) is a natural pathogen of rodents that is genetically related to both EBV and KSHV, and thus serves as an excellent model for the study of EBV and KSHV genetic elements such as miRNAs in the context of infection and disease. However, the specific targets of MHV68 miRNAs remain completely unknown. Using a technique known as qCLASH (quick crosslinking, ligation, and sequencing of hybrids), we have now identified thousands of Ago-associated, direct miRNA-mRNA interactions during lytic infection, latent infection and reactivation from latency. Validating this approach, detailed molecular analyses of specific interactions demonstrated repression of numerous host mRNA targets of MHV68 miRNAs, including Arid1a, Ctsl, Ifitm3 and Phc3. Notably, of the 1,505 MHV68 miRNA-host mRNA targets identified in B cells, 86% were shared with either EBV or KSHV, and 64% were shared among all three viruses, demonstrating significant conservation of gammaherpesvirus miRNA targeting. Pathway analysis of MHV68 miRNA targets further revealed enrichment of cellular pathways involved in protein synthesis and protein modification, including eIF2 Signaling, mTOR signaling and protein ubiquitination, pathways also enriched for targets of EBV and KSHV miRNAs. These findings provide substantial new information about specific targets of MHV68 miRNAs and shed important light on likely conserved functions of gammaherpesvirus miRNAs.

Gammaherpesviruses, including the human pathogens Epstein-Barr virus (EBV) and Kaposi’s sarcoma-associated herpesvirus (KSHV), establish lifelong infections and are associated with a variety of tumors. These viruses encode numerous molecules called microRNAs (miRNAs) within their genomes, which target and suppress the products of specific genes within infected host cells. However, the function of these miRNAs during in vivo infection is largely unknown. Murine gammaherpesvirus 68 (MHV68) is a natural pathogen of rodents that is genetically related to both EBV and KSHV, and thus serves as an excellent model for the study of EBV and KSHV. Here, we describe the identification and validation of thousands of new MHV68 miRNA targets. Notably, 86% of the MHV68 miRNA targets identified were shared with either EBV or KSHV, and 64% were shared among all three viruses. Further analyses revealed enrichment of cellular pathways involved in protein synthesis and protein modification, including pathways also enriched for targets of EBV and KSHV miRNAs. These findings provide substantial new information about specific targets of MHV68 miRNAs and shed important light on likely conserved functions of gammaherpesvirus miRNAs.

Gammaherpesvirus RNAs Come Full Circle

After an adaptive immune response is mounted, gammaherpesviruses achieve persistence through the utilization of viral noncoding RNAs to craft a suitable host cell environment in an immunologically transparent manner. While gammaherpesvirus long noncoding RNAs (lncRNAs) and microRNAs have been recognized for some time and have been actively investigated, a recent spate of reports have now identified repertoires of the circular RNA (circRNA) class of noncoding RNAs in both the lymphocryptovirus and rhadinovirus genera of gammaherpesviruses. Despite the recent nature of these findings, the detection of circRNAs across viruses and viral gene expression programs, the conservation of some viral circRNAs, and their detection in the clinical setting already raises the spectrum of functional importance in gammaherpesvirus biology and associated malignancies. Here, we provide an overview of currently known gammaherpesvirus circular RNAs and discuss reported physical and contextual properties that may be germane to future functional studies. With the Epstein-Barr virus (EBV) circRNAome being the most extensively studied to date, our discussions will be weighted toward EBV circRNAs while also addressing circRNAs discovered in the rhesus macaque lymphocryptovirus (rLCV), the Kaposi's sarcoma herpesvirus (KSHV), and the murid gammaherpesvirus 68 (MHV68). We hope that this will help set the stage for future investigations into the functions and relevance of this new class of viral noncoding RNAs in infection and disease.

Comparative Analysis of Gammaherpesvirus Circular RNA Repertoires: Conserved and Unique Viral Circular RNAs

Recent studies have identified circular RNAs (circRNAs) expressed from the Epstein-Barr virus (EBV) and Kaposi's sarcoma herpesvirus (KSHV) human DNA tumor viruses. To gain initial insights into the potential relevance of EBV circRNAs in virus biology and disease, we assessed the circRNAome of the interspecies homologue rhesus macaque lymphocryptovirus (rLCV) in a naturally occurring lymphoma from a simian immunodeficiency virus (SIV)-infected rhesus macaque. This analysis revealed rLCV orthologues of the latency-associated EBV circular RNAs circRPMS1_E4_E3a and circEBNA_U. Also identified in two samples displaying unusually high lytic gene expression was a novel rLCV circRNA that contains both conserved and rLCV-specific RPMS1 exons and whose backsplice junctions flank an rLCV lytic origin of replication (OriLyt). Analysis of a lytic infection model for the murid herpesvirus 68 (MHV68) rhadinovirus identified a cluster of circRNAs near an MHV68 lytic origin of replication, with the most abundant of these, circM11_ORF69, spanning the OriLyt. Lastly, analysis of KSHV latency and reactivation models revealed the latency associated circRNA originating from the vIRF4 gene as the predominant viral circRNA. Together, the results of this study broaden our appreciation for circRNA repertoires in the Lymphocryptovirus and Rhadinovirus genera of gammaherpesviruses and provide evolutionary support for viral circRNA functions in latency and viral replication.IMPORTANCE Infection with oncogenic gammaherpesviruses leads to long-term viral persistence through a dynamic interplay between the virus and the host immune system. Critical for remodeling of the host cell environment after the immune responses are viral noncoding RNAs that modulate host signaling pathways without attracting adaptive immune recognition. Despite the importance of noncoding RNAs in persistent infection, the circRNA class of noncoding RNAs has only recently been identified in gammaherpesviruses. Accordingly, their roles in virus infection and associated oncogenesis are unknown. Here we report evolutionary conservation of EBV-encoded circRNAs determined by assessing the circRNAome in rLCV-infected lymphomas from an SIV-infected rhesus macaque, and we report latent and lytic circRNAs from KSHV and MHV68. These experiments demonstrate utilization of the circular RNA class of RNAs across 4 members of the gammaherpesvirus subfamily, and they identify orthologues and potential homoplastic circRNAs, implying conserved circRNA functions in virus biology and associated malignancies.

Connivance, Complicity, or Collusion? The Role of Noncoding RNAs in Promoting Gammaherpesvirus Tumorigenesis

EBV and KSHV are etiologic agents of multiple types of lymphomas and carcinomas. The frequency of EBV⁺ or KSHV⁺ malignancies arising in immunocompromised individuals reflects the intricate evolutionary balance established between these viruses and their immunocompetent hosts. However, the specific mechanisms by which these pathogens drive tumorigenesis remain poorly understood. In recent years an enormous array of cellular and viral noncoding RNAs (ncRNAs) have been discovered, and host ncRNAs have been revealed as contributory factors to every single cancer hallmark cellular process. As new evidence emerges that gammaherpesvirus ncRNAs also alter host processes and viral factors dysregulate host ncRNA expression, and as novel viral ncRNAs continue to be discovered, we examine the contribution of small, non-miRNA ncRNAs and long ncRNAs to gammaherpesvirus tumorigenesis.

The Epstein Barr virus circRNAome

Our appreciation for the extent of Epstein Barr virus (EBV) transcriptome complexity continues to grow through findings of EBV encoded microRNAs, new long non-coding RNAs as well as the more recent discovery of over a hundred new polyadenylated lytic transcripts. Here we report an additional layer to the EBV transcriptome through the identification of a repertoire of latent and lytic viral circular RNAs. Utilizing RNase R-sequencing with cell models representing latency types I, II, and III, we identified EBV encoded circular RNAs expressed from the latency Cp promoter involving backsplicing from the W1 and W2 exons to the C1 exon, from the EBNA BamHI U fragment exon, and from the latency long non-coding RPMS1 locus. In addition, we identified circular RNAs expressed during reactivation including backsplicing from exon 8 to exon 2 of the LMP2 gene and a highly expressed circular RNA derived from intra-exonic backsplicing within the BHLF1 gene. While expression of most of these circular RNAs was found to depend on the EBV transcriptional program utilized and the transcription levels of the associated loci, expression of LMP2 exon 8 to exon 2 circular RNA was found to be cell model specific. Altogether we identified over 30 unique EBV circRNAs candidates and we validated and determined the structural features, expression profiles and nuclear/cytoplasmic distributions of several predominant and notable viral circRNAs. Further, we show that two of the EBV circular RNAs derived from the RPMS1 locus are detected in EBV positive clinical stomach cancer specimens. This study increases the known EBV latency and lytic transcriptome repertoires to include viral circular RNAs and it provides an essential foundation and resource for investigations into the functions and roles of this new class of EBV transcripts in EBV biology and diseases.

Our understanding of the extent of EBV transcriptome complexity has recently come to light with discoveries of viral microRNAs, snoRNAs, and a diverse set of transcript isoforms expressed during reactivation. While EBV utilizes distinct transcriptional programs during it’s normal infection cycle, the latency programs are critical drivers of EBV associated oncogenesis where they activate oncogenic pathways without killing the cell via lytic viral replication. Further, pathogenic latency programs typically express minimal levels of viral protein coding genes and instead rely on non-coding RNA mechanisms to alter the host cell environment. The utilization of non-coding RNAs in these settings is presumably an effort to minimize infected cell killing by the immune system through adaptive immune responses. Using tissue culture models representing the spectrum of EBV transcriptional states that occur in vivo, we report here that EBV expresses the circular class of predominantly non-coding RNAs. This stands to reveal a previously unknown set of players in EBV associated biology and disease. The findings reported here set a key platform from which future functional and potentially clinical related studies can be carried out to bring the initial findings from discovery to a new understanding of EBV biology and disease.

A KSHV microRNA enhances viral latency and induces angiogenesis by targeting GRK2 to activate the CXCR2/AKT pathway

Kaposi's sarcoma-associated herpesvirus (KSHV) is the causative agent of Kaposi's sarcoma (KS), primary effusion lymphoma (PEL) and multicentric Castleman's disease (MCD). Most tumor cells in these malignancies are latently infected by KSHV. Thus, viral latency is critical for the development of tumor and induction of tumor-associated angiogenesis. KSHV encodes more than two dozens of miRNAs but their roles in KSHV-induced angiogenesis remains unknown. We have recently shown that miR-K12-3 (miR-K3) promoted cell migration and invasion by targeting GRK2/CXCR2/AKT signaling (PLoS Pathog, 2015;11(9):e1005171). Here, we further demonstrated a role of miR-K3 and its induced signal pathway in KSHV latency and KSHV-induced angiogenesis. We found that overexpression of miR-K3 not only promoted viral latency by inhibiting viral lytic replication, but also induced angiogenesis. Further, knockdown of GRK2 inhibited KSHV replication and enhanced KSHV-induced angiogenesis by enhancing the CXCR2/AKT signals. As a result, blockage of CXCR2 or AKT increased KSHV replication and decreased angiogenesis induced by PEL cells in vivo. Finally, deletion of miR-K3 from viral genome reduced KSHV-induced angiogenesis and increased KSHV replication. These findings indicate that the miR-K3/GRK2/CXCR2/AKT axis plays an essential role in KSHV-induced angiogenesis and promotes KSHV latency, and thus may be a potential therapeutic target of KSHV-associated malignancies.

microRNA dependent and independent deregulation of long non-coding RNAs by an oncogenic herpesvirus

Kaposi’s sarcoma (KS) is a highly prevalent cancer in AIDS patients, especially in sub-Saharan Africa. Kaposi’s sarcoma-associated herpesvirus (KSHV) is the etiological agent of KS and other cancers like Primary Effusion Lymphoma (PEL). In KS and PEL, all tumors harbor latent KSHV episomes and express latency-associated viral proteins and microRNAs (miRNAs). The exact molecular mechanisms by which latent KSHV drives tumorigenesis are not completely understood. Recent developments have highlighted the importance of aberrant long non-coding RNA (lncRNA) expression in cancer. Deregulation of lncRNAs by miRNAs is a newly described phenomenon. We hypothesized that KSHV-encoded miRNAs deregulate human lncRNAs to drive tumorigenesis. We performed lncRNA expression profiling of endothelial cells infected with wt and miRNA-deleted KSHV and identified 126 lncRNAs as putative viral miRNA targets. Here we show that KSHV deregulates host lncRNAs in both a miRNA-dependent fashion by direct interaction and in a miRNA-independent fashion through latency-associated proteins. Several lncRNAs that were previously implicated in cancer, including MEG3, ANRIL and UCA1, are deregulated by KSHV. Our results also demonstrate that KSHV-mediated UCA1 deregulation contributes to increased proliferation and migration of endothelial cells.

KS is the most prevalent cancer associated with AIDS in sub-Saharan Africa, and is also common in males not affected by AIDS. KSHV manipulates human cells by targeting protein-coding genes and cell signaling. Here we show that KSHV alters the expression of hundreds of human lncRNAs, a broad class of regulatory molecules involved in a variety of cellular pathways including cell cycle and apoptosis. KSHV uses both latency proteins and miRNAs to target lncRNAs. miRNA-mediated targeting of lncRNAs is a novel regulatory mechanism of gene expression. Given that most herpesviruses encode miRNAs, this mechanism might be a common theme during herpesvirus infections. Understanding lncRNA deregulation by KSHV will help decipher the important molecular mechanisms underlying viral pathogenesis and tumorigenesis.

A viral microRNA downregulates metastasis suppressor CD82 and induces cell invasion and angiogenesis by activating the c-Met signaling

Kaposi’s sarcoma (KS) is the most common AIDS-associated malignancy etiologically caused by Kaposi’s sarcoma-associated herpesvirus (KSHV). KS is a highly disseminated and vascularized tumor comprised of poorly differentiated spindle-shaped endothelial cells. KSHV encodes 12 pre-microRNAs (pre-miRNAs) that yield 25 mature miRNAs, but their roles in KSHV-induced tumor dissemination and angiogenesis remain largely unknown. KSHV-encoded miR-K12-6 (miR-K6) can produce two mature miRNAs, miR-K6-3p and miR-K6-5p. Recently, we have shown that miR-K6-3p promoted cell migration and angiogenesis by directly targeting SH3 domain binding glutamate-rich protein (SH3BGR) (PLoS Pathog. 2016;12(4):e1005605). Here, by using mass spectrometry, bioinformatics analysis and luciferase reporter assay, we showed that miR-K6-5p directly targeted the coding sequence (CDS) of CD82 molecule (CD82), a metastasis suppressor. Ectopic expression of miR-K6-5p specifically inhibited the expression of endogenous CD82 and strongly promoted endothelial cells invasion in vitro and angiogenesis in vivo. Overexpression of CD82 significantly inhibited cell invasion and angiogenesis induced by miR-K6-5p. Mechanistically, CD82 directly interacted with c-Met to inhibit its activation. MiR-K6-5p directly repressed CD82, relieving its inhibition on c-Met activation and inducing cell invasion and angiogenesis. Deletion of miR-K6 from KSHV genome abrogated KSHV suppression of CD82 resulting in compromised KSHV activation of c-Met pathway, and KSHV-induced invasion and angiogenesis. In conclusion, these results show that by inhibiting CD82, KSHV miR-K6-5p promotes cell invasion and angiogenesis by activating the c-Met pathway. Our findings illustrate that KSHV miRNAs may play an essential role in the dissemination and angiogenesis of KSHV-induced malignancies.

Role of heme oxygenase-1 in the pathogenesis and tumorigenicity of Kaposi's sarcoma-associated herpesvirus

Kaposi's Sarcoma-associated Herpesvirus (KSHV) is the etiologic agent of several malignancies, including Kaposi's Sarcoma (KS), which preferentially arise in immunocompromised patients such as HIV+ subpopulation and lack effective therapeutic options. Heme oxygenase-1 (HO-1) has been reported as an important regulator of endothelial cell cycle control, proliferation and angiogenesis. HO-1 has also been found to be highly expressed in KSHV-infected endothelial cells and oral AIDS-KS lesions. We previously demonstrate that the multifunctional glycoprotein CD147 is required for KSHV/LANA-induced endothelial cell invasiveness. During the identification of CD147 controlled downstream genes by microarray analysis, we found that the expression of HO-1 is significantly elevated in both CD147-overexpressing and KSHV-infected HUVEC cells when compared to control cells. In the current study, we further identify the regulation of HO-1 expression and mediated cellular functions by both CD147 and KSHV-encoded LANA proteins. Targeting HO-1 by either RNAi or the chemical inhibitor, SnPP, effectively induces cell death of KSHV-infected endothelial cells (the major cellular components of KS) through DNA damage and necrosis process. By using a KS-like nude mouse model, we found that SnPP treatment significantly suppressed KSHV-induced tumorigenesis in vivo. Taken together, our data demonstrate the important role of HO-1 in the pathogenesis and tumorigenesis of KSHV-infected endothelial cells, the underlying regulatory mechanisms for HO-1 expression and targeting HO-1 may represent a promising therapeutic strategy against KSHV-related malignancies.

The SH3BGR/STAT3 Pathway Regulates Cell Migration and Angiogenesis Induced by a Gammaherpesvirus MicroRNA

Kaposi’s sarcoma (KS)-associated herpesvirus (KSHV) is a gammaherpesvirus etiologically associated with KS, a highly disseminated angiogenic tumor of hyperproliferative spindle endothelial cells. KSHV encodes 25 mature microRNAs but their roles in KSHV-induced tumor dissemination and angiogenesis remain unknown. Here, we investigated KSHV-encoded miR-K12-6-3p (miR-K6-3p) promotion of endothelial cell migration and angiogenesis, which are the underlying mechanisms of tumor dissemination and angiogenesis. We found that ectopic expression of miR-K6-3p promoted endothelial cell migration and angiogenesis. Mass spectrometry, bioinformatics and luciferase reporter analyses revealed that miR-K6-3p directly targeted sequence in the 3’ untranslated region (UTR) of SH3 domain binding glutamate-rich protein (SH3BGR). Overexpression of SH3BGR reversed miR-K6-3p induction of cell migration and angiogenesis. Mechanistically, miR-K6-3p downregulated SH3BGR, hence relieved STAT3 from SH3BGR direct binding and inhibition, which was required for miR-K6-3p maximum activation of STAT3 and induction of cell migration and angiogenesis. Finally, deletion of miR-K6 from the KSHV genome abrogated its effect on the SH3BGR/STAT3 pathway, and KSHV-induced migration and angiogenesis. Our results illustrated that, by inhibiting SH3BGR, miR-K6-3p enhances cell migration and angiogenesis by activating the STAT3 pathway, and thus contributes to the dissemination and angiogenesis of KSHV-induced malignancies.

Kaposi’s Sarcoma (KS), caused by infection of Kaposi’s sarcoma (KS)-associated herpesvirus (KSHV), is a tumor of endothelial cells characterized by angiogenesis and invasiveness. In vitro, KSHV-infected endothelial cells display an increased invasiveness and angiogenicity. KSHV encodes twelve precursor miRNAs (pre-miRNAs), which are processed into at least 25 mature miRNAs. However, the roles of these miRNAs in KSHV-induced tumor dissemination and angiogenesis remain unknown. Here, we investigated KSHV-encoded miR-K12-6-3p (miR-K6-3p) promotion of endothelial cell migration and angiogenesis, which are the underlying mechanisms of tumor dissemination and angiogenesis. We demonstrated that miR-K6-3p promoted cell migration and angiogenesis by directly targeting SH3 domain binding glutamate-rich protein (SH3BGR). Furthermore, we found that STAT3, which was negatively regulated by SH3BGR mediated miR-K6-3p-induced cell migration and angiogenesis. MiR-K6-3p downregulation of SH3BGR, hence relieved SH3BGR direct inhibition of STAT3 resulting in the activation of STAT3 and induction of cell migration and angiogenesis. These results identify miR-K6-3p and its the downstream pathway as potential therapeutic targets for the treatment of KSHV-associated malignancies.

A unifying gene signature for adenoid cystic cancer identifies parallel MYB-dependent and MYB-independent therapeutic targets

MYB activation is proposed to underlie development of adenoid cystic cancer (ACC), an aggressive salivary gland tumor with no effective systemic treatments. To discover druggable targets for ACC, we performed global mRNA/miRNA analyses of 12 ACC with matched normal tissues, and compared these data with 14 mucoepidermoid carcinomas (MEC) and 11 salivary adenocarcinomas (ADC). We detected a unique ACC gene signature of 1160 mRNAs and 22 miRNAs. MYB was the top-scoring gene (18-fold induction), however we observed the same signature in ACC without detectable MYB gene rearrangements. We also found 4 ACC tumors (1 among our 12 cases and 3 from public databases) with negligible MYB expression that retained the same ACC mRNA signature including over-expression of extracellular matrix (ECM) genes. Integration of this signature with somatic mutational analyses suggests that NOTCH1 and RUNX1 participate with MYB to activate ECM elements including the VCAN/HAPLN1 complex. We observed that forced MYB-NFIB expression in human salivary gland cells alters cell morphology and cell adhesion in vitro and depletion of VCAN blocked tumor cell growth of a short-term ACC tumor culture. In summary, we identified a unique ACC signature with parallel MYB-dependent and independent biomarkers and identified VCAN/HAPLN1 complexes as a potential target.

Kaposi's Sarcoma-Associated Herpesvirus (KSHV) Induces the Oncogenic miR-17-92 Cluster and Down-Regulates TGF-β Signaling

KSHV is a DNA tumor virus that causes Kaposi's sarcoma. Upon KSHV infection, only a limited number of latent genes are expressed. We know that KSHV infection regulates host gene expression, and hypothesized that latent genes also modulate the expression of host miRNAs. Aberrant miRNA expression contributes to the development of many types of cancer. Array-based miRNA profiling revealed that all six miRNAs of the oncogenic miR-17-92 cluster are up-regulated in KSHV infected endothelial cells. Among candidate KSHV latent genes, we found that vFLIP and vCyclin were shown to activate the miR-17-92 promoter, using luciferase assay and western blot analysis. The miR-17-92 cluster was previously shown to target TGF-β signaling. We demonstrate that vFLIP and vCyclin induce the expression of the miR-17-92 cluster to strongly inhibit the TGF-β signaling pathway by down-regulating SMAD2. Moreover, TGF-β activity and SMAD2 expression were fully restored when antagomirs (inhibitors) of miR-17-92 cluster were transfected into cells expressing either vFLIP or vCyclin. In addition, we utilized viral genetics to produce vFLIP or vCyclin knock-out viruses, and studied the effects in infected TIVE cells. Infection with wildtype KSHV abolished expression of SMAD2 protein in these endothelial cells. While single-knockout mutants still showed a marked reduction in SMAD2 expression, TIVE cells infected by a double-knockout mutant virus were fully restored for SMAD2 expression, compared to non-infected TIVE cells. Expression of either vFLIP or vCycIin was sufficient to downregulate SMAD2. In summary, our data demonstrate that vFLIP and vCyclin induce the oncogenic miR-17-92 cluster in endothelial cells and thereby interfere with the TGF-β signaling pathway. Manipulation of the TGF-β pathway via host miRNAs represents a novel mechanism that may be important for KSHV tumorigenesis and angiogenesis, a hallmark of KS.

A KSHV microRNA Directly Targets G Protein-Coupled Receptor Kinase 2 to Promote the Migration and Invasion of Endothelial Cells by Inducing CXCR2 and Activating AKT Signaling

Kaposi's sarcoma (KS) is a highly disseminated angiogenic tumor of endothelial cells linked to infection by Kaposi's sarcoma-associated herpesvirus (KSHV). KSHV encodes more than two dozens of miRNAs but their roles in KSHV-induced tumor dissemination and metastasis remain unknown. Here, we found that ectopic expression of miR-K12-3 (miR-K3) promoted endothelial cell migration and invasion. Bioinformatics and luciferase reporter analyses showed that miR-K3 directly targeted G protein-coupled receptor (GPCR) kinase 2 (GRK2, official gene symbol ADRBK1). Importantly, overexpression of GRK2 reversed miR-K3 induction of cell migration and invasion. Furthermore, the chemokine receptor CXCR2, which was negatively regulated by GRK2, was upregulated in miR-K3-transduced endothelial cells. Knock down of CXCR2 abolished miR-K3-induced cell migration and invasion. Moreover, miR-K3 downregulation of GRK2 relieved its direct inhibitory effect on AKT. Both CXCR2 induction and the release of AKT from GRK2 were required for miR-K3 maximum activation of AKT and induction of cell migration and invasion. Finally, deletion of miR-K3 from the KSHV genome abrogated its effect on the GRK2/CXCR2/AKT pathway and KSHV-induced migration and invasion. Our data provide the first-line evidence that, by repressing GRK2, miR-K3 facilitates cell migration and invasion via activation of CXCR2/AKT signaling, which likely contribute to the dissemination of KSHV-induced tumors.

Kaposi's sarcoma-associated herpesvirus (KSHV) is the etiological agent of Kaposi's sarcoma (KS). KS is a highly disseminated tumor often involved with visceral organs. Experimentally, KSHV infection induces the invasiveness of endothelial cells. KSHV encodes twelve precursor miRNAs (pre-miRNAs), which are processed into at least 25 mature miRNAs. However, the roles of these miRNAs in KSHV-induced tumor dissemination remain unknown. Here, we investigated KSHV-encoded miR-K12-3 (miR-K3) promotion of endothelial cell migration and invasion, which are the underlying mechanisms of tumor dissemination. We demonstrated that miR-K3 promoted cell migration and invasion by directly targeting G protein-coupled receptor (GPCR) kinase 2 (GRK2). Furthermore, we found that the chemokine receptor CXCR2, which was negatively regulated by GRK2, and its downstream AKT signaling positively mediated miR-K3-induced cell migration and invasion. miR-K3 downregulation of GRK2 relieved its direct inhibitory effect on AKT, and both CXCR2 induction and the release of AKT from GRK2 were required for miR-K3 maximum activation of AKT and induction of cell migration and invasion. These results show that miR-K3 and its the downstream pathway may be potential therapeutic targets for the treatment of KSHV-associated malignancies.

HITS-CLIP and PAR-CLIP advance viral miRNA targetome analysis

MiRNAs regulate gene expression by binding predominantly to the 3' untranslated region (UTR) of target transcripts to prevent their translation and/or induce target degradation. In addition to the more than 1200 human miRNAs, human DNA tumor viruses such as Kaposi's sarcoma-associated herpesvirus (KSHV) and Epstein-Barr virus (EBV) encode miRNAs. Target predictions indicate that each miRNA targets hundreds of transcripts, many of which are regulated by multiple miRNAs. Thus, target identification is a big challenge for the field. Most methods used currently investigate single miRNA-target interactions and are not able to analyze complex miRNA-target networks. To overcome these challenges, cross-linking and immunoprecipitation (CLIP), a recently developed method to study direct RNA-protein interactions in living cells, has been successfully applied to miRNA target analysis. It utilizes Argonaute (Ago)-immunoprecipitation to isolate native Ago-miRNA-mRNA complexes. In four recent publications, two variants of the CLIP method (HITS-CLIP and PAR-CLIP) were utilized to determine the targetomes of human and viral miRNAs in cells infected with the gamma-herpesviruses KSHV and EBV, which are associated with a number of human cancers. Here, we briefly introduce herpesvirus-encoded miRNAs and then focus on how CLIP technology has largely impacted our understanding of viral miRNAs in viral biology and pathogenesis.

A systems biology approach identified different regulatory networks targeted by KSHV miR-K12-11 in B cells and endothelial cells

Background

Kaposi’s sarcoma associated herpes virus (KSHV) is associated with tumors of endothelial and lymphoid origin. During latent infection, KSHV expresses miR-K12-11, an ortholog of the human tumor gene hsa-miR-155. Both gene products are microRNAs (miRNAs), which are important post-transcriptional regulators that contribute to tissue specific gene expression. Advances in target identification technologies and molecular interaction databases have allowed a systems biology approach to unravel the gene regulatory networks (GRNs) triggered by miR-K12-11 in endothelial and lymphoid cells. Understanding the tissue specific function of miR-K12-11 will help to elucidate underlying mechanisms of KSHV pathogenesis.

Results

Ectopic expression of miR-K12-11 differentially affected gene expression in BJAB cells of lymphoid origin and TIVE cells of endothelial origin. Direct miRNA targeting accounted for a small fraction of the observed transcriptome changes: only 29 genes were identified as putative direct targets of miR-K12-11 in both cell types. However, a number of commonly affected biological pathways, such as carbohydrate metabolism and interferon response related signaling, were revealed by gene ontology analysis. Integration of transcriptome profiling, bioinformatic algorithms, and databases of protein-protein interactome from the ENCODE project identified different nodes of GRNs utilized by miR-K12-11 in a tissue-specific fashion. These effector genes, including cancer associated transcription factors and signaling proteins, amplified the regulatory potential of a single miRNA, from a small set of putative direct targets to a larger set of genes.

Conclusions

This is the first comparative analysis of miRNA-K12-11’s effects in endothelial and B cells, from tissues infected with KSHV in vivo. MiR-K12-11 was able to broadly modulate gene expression in both cell types. Using a systems biology approach, we inferred that miR-K12-11 establishes its GRN by both repressing master TFs and influencing signaling pathways, to counter the host anti-viral response and to promote proliferation and survival of infected cells. The targeted GRNs are more reproducible and informative than target gene identification, and our approach can be applied to other regulatory factors of interest.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-668) contains supplementary material, which is available to authorized users.

LANA binds to multiple active viral and cellular promoters and associates with the H3K4methyltransferase hSET1 complex

Kaposi's sarcoma-associated herpesvirus (KSHV) is a γ-herpesvirus associated with KS and two lymphoproliferative diseases. Recent studies characterized epigenetic modification of KSHV episomes during latency and determined that latency-associated genes are associated with H3K4me3 while most lytic genes are associated with the silencing mark H3K27me3. Since the latency-associated nuclear antigen (LANA) (i) is expressed very early after de novo infection, (ii) interacts with transcriptional regulators and chromatin remodelers, and (iii) regulates the LANA and RTA promoters, we hypothesized that LANA may contribute to the establishment of latency through epigenetic control. We performed a detailed ChIP-seq analysis in cells of lymphoid and endothelial origin and compared H3K4me3, H3K27me3, polII, and LANA occupancy. On viral episomes LANA binding was detected at numerous lytic and latent promoters, which were transactivated by LANA using reporter assays. LANA binding was highly enriched at H3K4me3 peaks and this co-occupancy was also detected on many host gene promoters. Bioinformatic analysis of enriched LANA binding sites in combination with biochemical binding studies revealed three distinct binding patterns. A small subset of LANA binding sites showed sequence homology to the characterized LBS1/2 sequence in the viral terminal repeat. A large number of sites contained a novel LANA binding motif (TCCAT)₃ which was confirmed by gel shift analysis. Third, some viral and cellular promoters did not contain LANA binding sites and are likely enriched through protein/protein interaction. LANA was associated with H3K4me3 marks and in PEL cells 86% of all LANA bound promoters were transcriptionally active, leading to the hypothesis that LANA interacts with the machinery that methylates H3K4. Co-immunoprecipitation demonstrated LANA association with endogenous hSET1 complexes in both lymphoid and endothelial cells suggesting that LANA may contribute to the epigenetic profile of KSHV episomes.

KSHV is a DNA tumor virus which is associated with Kaposi's sarcoma and some lymphoproliferative diseases. During latent infection, the viral genome persists as circular extrachromosomal DNA in the nucleus and expresses a very limited number of viral proteins, including LANA, a multi-functional protein. KSHV viral episomes, like host genomic DNA, are subject to chromatin formation and histone modifications which contribute to tightly controlled gene expression during latency. We determined where LANA binds on the KSHV and human genomes, and mapped activating and repressing histone marks and RNA polymerase II binding. We found that LANA bound near transcription start sites, and binding correlated with the transcription active mark H3K4me3, but not silencing mark H3K27me3. Binding sites for transcription factors including znf143, CTCF, and Stat1 are enriched at regions where LANA is bound. We identified some LANA binding sites near human gene promoters that resembled KSHV sequences known to bind LANA. We also found a novel motif that occurs frequently in the human genome and that binds LANA directly despite being different from known LANA-binding sequences. Furthermore, we demonstrate that LANA associates with the H3K4 methyltransferase hSET1 which creates activating histone marks.

Helminth infection reactivates latent γ-herpesvirus via cytokine competition at a viral promoter

Mammals are coinfected by multiple pathogens that interact through unknown mechanisms. We found that helminth infection, characterized by the induction of the cytokine interleukin-4 (IL-4) and the activation of the transcription factor Stat6, reactivated murine γ-herpesvirus infection in vivo. IL-4 promoted viral replication and blocked the antiviral effects of interferon-γ (IFNγ) by inducing Stat6 binding to the promoter for an important viral transcriptional transactivator. IL-4 also reactivated human Kaposi's sarcoma-associated herpesvirus from latency in cultured cells. Exogenous IL-4 plus blockade of IFNγ reactivated latent murine γ-herpesvirus infection in vivo, suggesting a "two-signal" model for viral reactivation. Thus, chronic herpesvirus infection, a component of the mammalian virome, is regulated by the counterpoised actions of multiple cytokines on viral promoters that have evolved to sense host immune status.

Genomic and proteomic analysis of transcription factor TFII-I reveals insight into the response to cellular stress

The ubiquitously expressed transcription factor TFII-I exerts both positive and negative effects on transcription. Using biotinylation tagging technology and high-throughput sequencing, we determined sites of chromatin interactions for TFII-I in the human erythroleukemia cell line K562. This analysis revealed that TFII-I binds upstream of the transcription start site of expressed genes, both upstream and downstream of the transcription start site of repressed genes, and downstream of RNA polymerase II peaks at the ATF3 and other stress responsive genes. At the ATF3 gene, TFII-I binds immediately downstream of a Pol II peak located 5 kb upstream of exon 1. Induction of ATF3 expression increases transcription throughout the ATF3 gene locus which requires TFII-I and correlates with increased association of Pol II and Elongin A. Pull-down assays demonstrated that TFII-I interacts with Elongin A. Partial depletion of TFII-I expression caused a reduction in the association of Elongin A with and transcription of the DNMT1 and EFR3A genes without a decrease in Pol II recruitment. The data reveal different interaction patterns of TFII-I at active, repressed, or inducible genes, identify novel TFII-I interacting proteins, implicate TFII-I in the regulation of transcription elongation and provide insight into the role of TFII-I during the response to cellular stress.

Leveraging biological replicates to improve analysis in ChIP-seq experiments

ChIP-seq experiments identify genome-wide profiles of DNA-binding molecules including transcription factors, enzymes and epigenetic marks. Biological replicates are critical for reliable site discovery and are required for the deposition of data in the ENCODE and modENCODE projects. While early reports suggested two replicates were sufficient, the widespread application of the technique has led to emerging consensus that the technique is noisy and that increasing replication may be worthwhile. Additional biological replicates also allow for quantitative assessment of differences between conditions. To date it has remained controversial about how to confirm peak identification and to determine signal strength across biological replicates, particularly when the number of replicates is greater than two. Using objective metrics, we evaluate the consistency of biological replicates in ChIP-seq experiments with more than two replicates. We compare several approaches for binding site determination, including two popular but disparate peak callers, CisGenome and MACS2. Here we propose read coverage as a quantitative measurement of signal strength for estimating sample concordance. Determining binding based on genomic features, such as promoters, is also examined. We find that increasing the number of biological replicates increases the reliability of peak identification. Critically, binding sites with strong biological evidence may be missed if researchers rely on only two biological replicates. When more than two replicates are performed, a simple majority rule (>50% of samples identify a peak) identifies peaks more reliably in all biological replicates than the absolute concordance of peak identification between any two replicates, further demonstrating the utility of increasing replicate numbers in ChIP-seq experiments.

Persistent human herpesvirus-6 infection in patients with an inherited form of the virus

Human herpesvirus-6 (HHV-6)A and 6B are ubiquitous betaherpesviruses viruses with lymphotropic and neurotropic potential. As reported earlier, these viruses establish latency by integration into the telomeres of host chromosomes. Chromosomally integrated HHV-6 (CIHHV-6) can be transmitted vertically from parent to child. Some CIHHV-6 patients are suffering from neurological symptoms, while others remain asymptomatic. Four patients with CIHHV-6 and CNS dysfunction were treated with valganciclovir or foscarnet. HHV-6 replication was detected by reverse transcriptase polymerase chain reaction amplification of a late envelope glycoprotein. In this study we also compared the inherited and persistent HHV-6 viruses by DNA sequencing. The prevalence of CIHHV-6 in this cohort of adult patients from the USA suffering from a wide range of neurological symptoms including long-term fatigue were found significantly greater than the reported 0.8% in the general population. Long-term antiviral therapy inhibited HHV-6 replication as documented by loss of viral mRNA production. Sequence comparison of the mRNA and the inherited viral genome revealed that the transcript is produced by an exogenous virus. In conclusion, the data presented here document that some individuals with CIHHV-6 are infected persistently with exogenous HHV-6 strains that lead to a wide range of neurological symptoms; the proposed name for this condition is inherited herpesvirus 6 syndrome or IHS.

Epigenetic diversity of Kaposi's sarcoma-associated herpesvirus

Spontaneous lytic reactivation of Kaposi’s sarcoma–associated herpesvirus (KSHV) occurs at a low rate in latently infected cells in disease and culture. This suggests imperfect epigenetic maintenance of viral transcription programs, perhaps due to variability in chromatin structure at specific loci across the population of KSHV episomal genomes. To characterize this locus-specific chromatin structural diversity, we used MAPit single-molecule footprinting, which simultaneously maps endogenous CG methylation and accessibility to M.CviPI at GC sites. Diverse chromatin structures were detected at the LANA, RTA and vIL6 promoters. At each locus, chromatin ranged from fully closed to fully open across the population. This diversity has not previously been reported in a virus. Phorbol ester and RTA transgene induction were used to identify chromatin conformations associated with reactivation of lytic transcription, which only a fraction of episomes had. Moreover, certain chromatin conformations correlated with CG methylation patterns at the RTA and vIL6 promoters. This indicated that some of the diverse chromatin conformations at these loci were epigenetically distinct. Finally, by comparing chromatin structures from a cell line infected with constitutively latent virus, we identified products of lytic replication. Our findings show that epigenetic drift can restrict viral propagation by chromatin compaction at latent and lytic promoters.

Viral miRNAs and immune evasion

Viral miRNAs, ~22nt RNA molecules which post-transcriptionally regulate gene expression, are emerging as important tools in immune evasion. Viral infection is a complex process that requires immune evasion in order to establish persistent life-long infection of the host. During this process viruses express both protein-coding and non-coding genes, which help to modulate the cellular environment making it more favorable for infection. In the last decade, it was uncovered that DNA viruses express a diverse and abundant pool of small non-coding RNA molecules, called microRNAs (miRNAs). These virally encoded miRNAs are non-immunogenic and therefore are important tools used to evade both innate and adaptive immune responses. This review aims to summarize our current knowledge of herpesvirus- and polyomavirus-encoded miRNAs, and how they contribute to immune evasion by targeting viral and/or host cellular genes. This article is part of a Special Issue entitled: MicroRNAs in viral gene regulation.

Viral miRNAs

Since 2004, more than 200 microRNAs (miRNAs) have been discovered in double-stranded DNA viruses, mainly herpesviruses and polyomaviruses (Nucleic Acids Res 32:D109-D111, 2004). miRNAs are short 22  ±  3 nt RNA molecules that posttranscriptionally regulate gene expression by binding to 3'-untranslated regions (3'UTR) of target mRNAs, thereby inducing translational silencing and/or transcript degradation (Nature 431:350-355, 2004; Cell 116:281-297, 2004). Since miRNAs require only limited complementarity for binding, miRNA targets are difficult to determine (Mol Cell 27:91-105, 2007). To date, targets have only been experimentally verified for relatively few viral miRNAs, which either target viral or host cellular gene expression: For example, SV40 and related polyomaviruses encode miRNAs which target viral large T antigen expression (Nature 435:682-686, 2005; J Virol 79:13094-13104, 2005; Virology 383:183-187, 2009; J Virol 82:9823-9828, 2008) and miRNAs of α-, β-, and γ-herpesviruses have been implicated in regulating the transition from latent to lytic gene expression, a key step in the herpesvirus life cycle. Viral miRNAs have also been shown to target various host cellular genes. Although this field is just beginning to unravel the multiple roles of viral miRNA in biology and pathogenesis, the current data strongly suggest that virally encoded miRNAs are able to regulate fundamental biological processes such as immune recognition, promotion of cell survival, angiogenesis, proliferation, and cell differentiation. This chapter aims to summarize our current knowledge of viral miRNAs, their targets and function, and the challenges lying ahead to decipher their role in viral biology, pathogenesis, and for γ-herepsvirus-encoded miRNAs, potentially tumorigenesis.

Viral miRNAs: tools for immune evasion

MicroRNAs (miRNAs) are noncoding RNA molecules approximately 22 nucleotides in length that post-transcriptionally regulate gene expression by complementary binding to target mRNAs. MiRNAs have been identified in a diverse range of both metazoan and plant species. Functionally, miRNAs modulate multiple cellular processes including development, hematopoiesis, immunity, and oncogenesis. More recently, DNA viruses were found to encode and express miRNAs during host infection. Although the functions of most viral miRNAs are not well understood, early analysis of target genes pointed to immune modulation suggesting that viral miRNAs are a component of the immune evasion repertoire, which facilitates viral persistence. In addition to directly targeting immune functions, viral encoded miRNAs contribute to immune evasion by targeting proapoptotic genes, and in the case of herpesviruses, by controlling viral latency. Here we summarize the recently discovered targets of viral miRNAs and discuss the complex nature of this novel emerging regulatory mechanism.