MPAC: a computational framework for inferring cancer pathway activities from multi-omic data.

Fully capturing cellular state requires examining genomic, epigenomic, transcriptomic, proteomic, and other assays for a biological sample and comprehensive computational modeling to reason with the complex and sometimes conflicting measurements. Modeling these so-called multi-omic data is especially beneficial in disease analysis, where observations across omic data types may reveal unexpected patient groupings and inform clinical outcomes and treatments. We present Multi-omic Pathway Analysis of Cancer (MPAC), a computational framework that interprets multi-omic data through prior knowledge from biological pathways. MPAC uses network relationships encoded in pathways using a factor graph to infer consensus activity levels for proteins and associated pathway entities from multi-omic data, runs permutation testing to eliminate spurious activity predictions, and groups biological samples by pathway activities to prioritize proteins with potential clinical relevance. Using DNA copy number alteration and RNA-seq data from head and neck squamous cell carcinoma patients from The Cancer Genome Atlas as an example, we demonstrate that MPAC predicts a patient subgroup related to immune responses not identified by analysis with either input omic data type alone. Key proteins identified via this subgroup have pathway activities related to clinical outcome as well as immune cell compositions. Our MPAC R package, available at https://bioconductor.org/packages/MPAC, enables similar multi-omic analyses on new datasets.

Key aspects of papillomavirus infection influence the host cervicovaginal microbiome in a preclinical murine papillomavirus (MmuPV1) infection model.

Human papillomaviruses (HPVs) are the most common sexually transmitted infection in the United States and are a major etiological agent of cancers in the anogenital tract and oral cavity. Growing evidence suggests changes in the host microbiome are associated with the natural history and ultimate outcome of HPV infection. We sought to define changes in the host cervicovaginal microbiome during papillomavirus infection, persistence, and pathogenesis using the murine papillomavirus (MmuPV1) cervicovaginal infection model. Cervicovaginal lavages were performed over a time course of MmuPV1 infection in immunocompetent female FVB/N mice and extracted DNA was analyzed by qPCR to track MmuPV1 viral copy number. 16S ribosomal RNA (rRNA) gene sequencing was used to determine the composition and diversity of microbial communities throughout this time course. We also sought to determine whether specific microbial communities exist across the spectrum of MmuPV1-induced neoplastic disease. We, therefore, performed laser-capture microdissection to isolate regions of disease representing all stages of neoplastic disease progression (normal, low- and high-grade dysplasia, and cancer) from female reproductive tract tissue sections from MmuPV1-infected mice and performed 16S rRNA sequencing. Consistent with other studies, we found that the natural murine cervicovaginal microbiome is highly variable across different experiments. Despite these differences in initial microbiome composition between experiments, we observed that MmuPV1 persistence, viral load, and severity of disease influenced the composition of the cervicovaginal microbiome. These studies demonstrate that papillomavirus infection can alter the cervicovaginal microbiome.IMPORTANCEHuman papillomaviruses (HPVs) are the most common sexually transmitted infection in the United States. A subset of HPVs that infect the anogenital tract (cervix, vagina, anus) and oral cavity cause at least 5% of cancers worldwide. Recent evidence indicates that the community of microbial organisms present in the human cervix and vagina, known as the cervicovaginal microbiome, plays a role in HPV-induced cervical cancer. However, the mechanisms underlying this interplay are not well-defined. In this study, we infected the female reproductive tract of mice with a murine papillomavirus (MmuPV1) and found that key aspects of papillomavirus infection and disease influence the host cervicovaginal microbiome. This is the first study to define changes in the host microbiome associated with MmuPV1 infection in a preclinical animal model of HPV-induced cervical cancer. These results pave the way for using MmuPV1 infection models to further investigate the interactions between papillomaviruses and the host microbiome.

Positive-strand RNA virus genome replication organelles: structure, assembly, control.

Positive-strand RNA [(+)RNA] viruses include pandemic SARS-CoV-2, tumor-inducing hepatitis C virus, debilitating chikungunya virus (CHIKV), lethal encephalitis viruses, and many other major pathogens. (+)RNA viruses replicate their RNA genomes in virus-induced replication organelles (ROs) that also evolve new viral species and variants by recombination and mutation and are crucial virus control targets. Recent cryo-electron microscopy (cryo-EM) reveals that viral RNA replication proteins form striking ringed 'crowns' at RO vesicle junctions with the cytosol. These crowns direct RO vesicle formation, viral (-)RNA and (+)RNA synthesis and capping, innate immune escape, and transfer of progeny (+)RNA genomes into translation and encapsidation. Ongoing studies are illuminating crown assembly, sequential functions, host factor interactions, etc., with significant implications for control and beneficial uses of viruses.

Exploiting a rodent cell block for intrinsic resistance to HIV-1 gene expression in human T cells.

IMPORTANCE: Unlike humans, mice are unable to support HIV-1 infection. This is due, in part, to a constellation of defined minor, species-specific differences in conserved host proteins needed for viral gene expression. Here, we used precision CRISPR/Cas9 gene editing to engineer a "mousified" version of one such host protein, cyclin T1 (CCNT1), in human T cells. CCNT1 is essential for efficient HIV-1 transcription, making it an intriguing target for gene-based inactivation of virus replication. We show that isogenic cell lines engineered to encode CCNT1 bearing a single mouse-informed amino acid change (tyrosine in place of cysteine at position 261) exhibit potent, durable, and broad-spectrum resistance to HIV-1 and other pathogenic lentiviruses, and with no discernible impact on host cell biology. These results provide proof of concept for targeting CCNT1 in the context of one or more functional HIV-1 cure strategies.

HIV-1 virological synapse formation enhances infection spread by dysregulating Aurora Kinase B.

HIV-1 spreads efficiently through direct cell-to-cell transmission at virological synapses (VSs) formed by interactions between HIV-1 envelope proteins (Env) on the surface of infected cells and CD4 receptors on uninfected target cells. Env-CD4 interactions bring the infected and uninfected cellular membranes into close proximity and induce transport of viral and cellular factors to the VS for efficient virion assembly and HIV-1 transmission. Using novel, cell-specific stable isotope labeling and quantitative mass spectrometric proteomics, we identified extensive changes in the levels and phosphorylation states of proteins in HIV-1 infected producer cells upon mixing with CD4+ target cells under conditions inducing VS formation. These coculture-induced alterations involved multiple cellular pathways including transcription, TCR signaling and, unexpectedly, cell cycle regulation, and were dominated by Env-dependent responses. We confirmed the proteomic results using inhibitors targeting regulatory kinases and phosphatases in selected pathways identified by our proteomic analysis. Strikingly, inhibiting the key mitotic regulator Aurora kinase B (AURKB) in HIV-1 infected cells significantly increased HIV activity in cell-to-cell fusion and transmission but had little effect on cell-free infection. Consistent with this, we found that AURKB regulates the fusogenic activity of HIV-1 Env. In the Jurkat T cell line and primary T cells, HIV-1 Env:CD4 interaction also dramatically induced cell cycle-independent AURKB relocalization to the centromere, and this signaling required the long (150 aa) cytoplasmic C-terminal domain (CTD) of Env. These results imply that cytoplasmic/plasma membrane AURKB restricts HIV-1 envelope fusion, and that this restriction is overcome by Env CTD-induced AURKB relocalization. Taken together, our data reveal a new signaling pathway regulating HIV-1 cell-to-cell transmission and potential new avenues for therapeutic intervention through targeting the Env CTD and AURKB activity.

Nodavirus RNA replication crown architecture reveals proto-crown precursor and viral protein A conformational switching.

Positive-strand RNA viruses replicate their genomes in virus-induced membrane vesicles, and the resulting RNA replication complexes are a major target for virus control. Nodavirus studies first revealed viral RNA replication proteins forming a 12-fold symmetric "crown" at the vesicle opening to the cytosol, an arrangement recently confirmed to extend to distantly related alphaviruses. Using cryoelectron microscopy (cryo-EM), we show that mature nodavirus crowns comprise two stacked 12-mer rings of multidomain viral RNA replication protein A. Each ring contains an ~19 nm circle of C-proximal polymerase domains, differentiated by strikingly diverged positions of N-proximal RNA capping/membrane binding domains. The lower ring is a "proto-crown" precursor that assembles prior to RNA template recruitment, RNA synthesis, and replication vesicle formation. In this proto-crown, the N-proximal segments interact to form a toroidal central floor, whose 3.1 Å resolution structure reveals many mechanistic details of the RNA capping/membrane binding domains. In the upper ring, cryo-EM fitting indicates that the N-proximal domains extend radially outside the polymerases, forming separated, membrane-binding "legs." The polymerase and N-proximal domains are connected by a long linker accommodating the conformational switch between the two rings and possibly also polymerase movements associated with RNA synthesis and nonsymmetric electron density in the lower center of mature crowns. The results reveal remarkable viral protein multifunctionality, conformational flexibility, and evolutionary plasticity and insights into (+)RNA virus replication and control.

Multifunctional Protein A Is the Only Viral Protein Required for Nodavirus RNA Replication Crown Formation.

Positive-strand RNA virus RNA genome replication occurs in membrane-associated RNA replication complexes (RCs). Nodavirus RCs are outer mitochondrial membrane invaginations whose necked openings to the cytosol are "crowned" by a 12-fold symmetrical proteinaceous ring that functions as the main engine of RNA replication. Similar protein crowns recently visualized at the openings of alphavirus and coronavirus RCs highlight their broad conservation and functional importance. Using cryo-EM tomography, we earlier showed that the major nodavirus crown constituent is viral protein A, whose polymerase, RNA capping, membrane interaction and multimerization domains drive RC formation and function. Other viral proteins are strong candidates for unassigned EM density in the crown. RNA-binding RNAi inhibitor protein B2 co-immunoprecipitates with protein A and could form crown subdomains that protect nascent viral RNA and dsRNA templates. Capsid protein may interact with the crown since nodavirus virion assembly has spatial and other links to RNA replication. Using cryoelectron tomography and complementary approaches, we show that, even when formed in mammalian cells, nodavirus RC crowns generated without B2 and capsid proteins are functional and structurally indistinguishable from mature crowns in infected Drosophila cells expressing all viral proteins. Thus, the only nodaviral factors essential to form functional RCs and crowns are RNA replication protein A and an RNA template. We also resolve apparent conflicts in prior results on B2 localization in infected cells, revealing at least two distinguishable pools of B2. The results have significant implications for crown structure, assembly, function and control as an antiviral target.

HIV RGB: Automated Single-Cell Analysis of HIV-1 Rev-Dependent RNA Nuclear Export and Translation Using Image Processing in KNIME.

Single-cell imaging has emerged as a powerful means to study viral replication dynamics and identify sites of virus−host interactions. Multivariate aspects of viral replication cycles yield challenges inherent to handling large, complex imaging datasets. Herein, we describe the design and implementation of an automated, imaging-based strategy, "Human Immunodeficiency Virus Red-Green-Blue" (HIV RGB), for deriving comprehensive single-cell measurements of HIV-1 unspliced (US) RNA nuclear export, translation, and bulk changes to viral RNA and protein (HIV-1 Rev and Gag) subcellular distribution over time. Differentially tagged fluorescent viral RNA and protein species are recorded using multicolor long-term (>24 h) time-lapse video microscopy, followed by image processing using a new open-source computational imaging workflow dubbed "Nuclear Ring Segmentation Analysis and Tracking" (NR-SAT) based on ImageJ plugins that have been integrated into the Konstanz Information Miner (KNIME) analytics platform. We describe a typical HIV RGB experimental setup, detail the image acquisition and NR-SAT workflow accompanied by a step-by-step tutorial, and demonstrate a use case wherein we test the effects of perturbing subcellular localization of the Rev protein, which is essential for viral US RNA nuclear export, on the kinetics of HIV-1 late-stage gene regulation. Collectively, HIV RGB represents a powerful platform for single-cell studies of HIV-1 post-transcriptional RNA regulation. Moreover, we discuss how similar NR-SAT-based design principles and open-source tools might be readily adapted to study a broad range of dynamic viral or cellular processes.

Crowning Touches in Positive-Strand RNA Virus Genome Replication Complex Structure and Function.

Positive-strand RNA viruses, the largest genetic class of eukaryotic viruses, include coronaviruses and many other established and emerging pathogens. A major target for understanding and controlling these viruses is their genome replication, which occurs in virus-induced membrane vesicles that organize replication steps and protect double-stranded RNA intermediates from innate immune recognition. The structure of these complexes has been greatly illuminated by recent cryo-electron microscope tomography studies with several viruses. One key finding in diverse systems is the organization of crucial viral RNA replication factors in multimeric rings or crowns that among other functions serve as exit channels gating release of progeny genomes to the cytosol for translation and encapsidation. Emerging results suggest that these crowns serve additional important purposes in replication complex assembly, function, and interaction with downstream processes such as encapsidation. The findings provide insights into viral function and evolution and new bases for understanding, controlling, and engineering positive-strand RNA viruses.

Human cytomegalovirus lytic infection inhibits replication-dependent histone synthesis and requires stem loop binding protein function.

Replication-dependent (RD) histones are deposited onto human cytomegalovirus (HCMV) genomes at the start of infection. We examined how HCMV affects the de novo production of RD histones and found that viral infection blocked the accumulation of RD histone mRNAs that normally occurs during the S phase. Furthermore, RD histone mRNAs present in HCMV-infected cells did not undergo the unique 3' processing required for their normal nuclear export and translation. The protein that orchestrates processing in the nucleus, stem loop­binding protein (SLBP), was found predominantly in the cytoplasm, and RD histone proteins were not de novo synthesized in HCMV-infected cells. Intriguingly, however, we found that SLBP was required for the efficient synthesis and assembly of infectious progeny virions. We conclude that HCMV infection attenuates RD histone mRNA accumulation and processing and the de novo protein synthesis of the RD histones, while utilizing SLBP for an alternative purpose to support infectious virion production.

Perturbing HIV-1 Ribosomal Frameshifting Frequency Reveals a cis Preference for Gag-Pol Incorporation into Assembling Virions.

HIV-1 virion production is driven by Gag and Gag-Pol (GP) proteins, with Gag forming the bulk of the capsid and driving budding, while GP binds Gag to deliver the essential virion enzymes protease, reverse transcriptase, and integrase. Virion GP levels are traditionally thought to reflect the relative abundances of GP and Gag in cells (∼1:20), dictated by the frequency of a -1 programmed ribosomal frameshifting (PRF) event occurring in gag-pol mRNAs. Here, we exploited a panel of PRF mutant viruses to show that mechanisms in addition to PRF regulate GP incorporation into virions. First, we show that GP is enriched ∼3-fold in virions relative to cells, with viral infectivity being better maintained at subphysiological levels of GP than when GP levels are too high. Second, we report that GP is more efficiently incorporated into virions when Gag and GP are synthesized in cis (i.e., from the same gag-pol mRNA) than in trans, suggesting that Gag/GP translation and assembly are spatially coupled processes. Third, we show that, surprisingly, virions exhibit a strong upper limit to trans-delivered GP incorporation; an adaptation that appears to allow the virus to temper defects to GP/Gag cleavage that may negatively impact reverse transcription. Taking these results together, we propose a "weighted Goldilocks" scenario for HIV-1 GP incorporation, wherein combined mechanisms of GP enrichment and exclusion buffer virion infectivity over a broad range of local GP concentrations. These results provide new insights into the HIV-1 virion assembly pathway relevant to the anticipated efficacy of PRF-targeted antiviral strategies. IMPORTANCE HIV-1 infectivity requires incorporation of the Gag-Pol (GP) precursor polyprotein into virions during the process of virus particle assembly. Mechanisms dictating GP incorporation into assembling virions are poorly defined, with GP levels in virions traditionally thought to solely reflect relative levels of Gag and GP expressed in cells, dictated by the frequency of a -1 programmed ribosomal frameshifting (PRF) event that occurs in gag-pol mRNAs. Herein, we provide experimental support for a "weighted Goldilocks" scenario for GP incorporation, wherein the virus exploits both random and nonrandom mechanisms to buffer infectivity over a wide range of GP expression levels. These mechanistic data are relevant to ongoing efforts to develop antiviral strategies targeting PRF frequency and/or HIV-1 virion maturation.

Cryo-electron microscopy of nodavirus RNA replication organelles illuminates positive-strand RNA virus genome replication.

The nodavirus flock house virus recently provided a well-characterized model for the first cryo-electron microscope tomography of membrane-bound, positive-strand RNA ((+)RNA) virus genome replication complexes (RCs). The resulting first views of RC organization and complementary biochemical results showed that the viral RNA replication vesicle is tightly packed with the dsRNA genomic RNA replication intermediate, and that (+)ssRNA replication products are released through the vesicle neck to the cytosol through a 12-fold symmetric ring or crown of multi-functional viral RNA replication proteins, which likely also contribute to viral RNA synthesis. Subsequent studies identified similar crown-like RNA replication protein complexes in alphavirus and coronavirus RCs, indicating related mechanisms across highly divergent (+)RNA viruses. As outlined in this review, these results have significant implications for viral function, evolution and control.

ESCRT components ISTL1 andLIP5 are required for tapetal function and pollen viability.

Pollen wall assembly is crucial for pollen development and plant fertility. The durable biopolymer sporopollenin and the constituents of the tryphine coat are delivered to developing pollen grains by the highly coordinated secretory activity of the surrounding tapetal cells. The role of membrane trafficking in this process, however, is largely unknown. In this study, we used Arabidopsis thaliana to characterize the role of two late-acting endosomal sorting complex required for transport (ESCRT) components, ISTL1 and LIP5, in tapetal function. Plants lacking ISTL1 and LIP5 form pollen with aberrant exine patterns, leading to partial pollen lethality. We found that ISTL1 and LIP5 are required for exocytosis of plasma membrane and secreted proteins in the tapetal cells at the free microspore stage, contributing to pollen wall development and tryphine deposition. Whereas the ESCRT machinery is well known for its role in endosomal trafficking, the function of ISTL1 and LIP5 in exocytosis is not a typical ESCRT function. The istl1 lip5 double mutants also show reduced intralumenal vesicle concatenation in multivesicular endosomes in both tapetal cells and developing pollen grains as well as morphological defects in early endosomes/trans-Golgi networks, suggesting that late ESCRT components function in the early endosomal pathway and exocytosis.

ZBTB2 represses HIV-1 transcription and is regulated by HIV-1 Vpr and cellular DNA damage responses.

Previously, we reported that cellular transcription factor ZASC1 facilitates DNA-dependent/RNA-independent recruitment of HIV-1 TAT and the cellular elongation factor P-TEFb to the HIV-1 promoter and is a critical factor in regulating HIV-1 transcriptional elongation (PLoS Path e1003712). Here we report that cellular transcription factor ZBTB2 is a novel repressor of HIV-1 gene expression. ZBTB2 strongly co-immunoprecipitated with ZASC1 and was dramatically relocalized by ZASC1 from the cytoplasm to the nucleus. Mutations abolishing ZASC1/ZBTB2 interaction prevented ZBTB2 nuclear relocalization. We show that ZBTB2-induced repression depends on interaction of cellular histone deacetylases (HDACs) with the ZBTB2 POZ domain. Further, ZASC1 interaction specifically recruited ZBTB2 to the HIV-1 promoter, resulting in histone deacetylation and transcription repression. Depleting ZBTB2 by siRNA knockdown or CRISPR/CAS9 knockout in T cell lines enhanced transcription from HIV-1 vectors lacking Vpr, but not from these vectors expressing Vpr. Since HIV-1 Vpr activates the viral LTR by inducing the ATR kinase/DNA damage response pathway, we investigated ZBTB2 response to Vpr and DNA damaging agents. Expressing Vpr or stimulating the ATR pathway with DNA damaging agents impaired ZASC1's ability to localize ZBTB2 to the nucleus. Moreover, the effects of DNA damaging agents and Vpr on ZBTB2 localization could be blocked by ATR kinase inhibitors. Critically, Vpr and DNA damaging agents decreased ZBTB2 binding to the HIV-1 promoter and increased promoter histone acetylation. Thus, ZBTB2 is recruited to the HIV-1 promoter by ZASC1 and represses transcription, but ATR pathway activation leads to ZBTB2 removal from the promoter, cytoplasmic sequestration and activation of viral transcription. Together, our data show that ZASC1/ZBTB2 integrate the functions of TAT and Vpr to maximize HIV-1 gene expression.

Transmembrane redox regulation of genome replication functions in positive-strand RNA viruses.

Positive-strand RNA virus genome replication takes place on intracellular membranes that separate the reduced cytosol from the oxidized extracellular/luminal milieu. Ongoing studies of these membrane-bounded genome replication complexes have revealed underlying common principles in their structure, assembly and functionalization, including transmembrane features and redox dependencies. Among these, members of the alphavirus, flavivirus, and picornavirus supergroups all encode membrane-permeabilizing viroporins required for efficient RNA replication. For flaviviruses and particularly alphavirus supergroup members, these viroporins are linked to activating viral RNA capping and potentially other later-stage RNA replication functions, and to local transmembrane release of oxidizing potential to trigger these changes in cytoplasmic RNA replication complexes. Further exploration of these emerging shared principles could spur development of broad-spectrum antivirals.

Subdomain cryo-EM structure of nodaviral replication protein A crown complex provides mechanistic insights into RNA genome replication.

For positive-strand RNA [(+)RNA] viruses, the major target for antiviral therapies is genomic RNA replication, which occurs at poorly understood membrane-bound viral RNA replication complexes. Recent cryoelectron microscopy (cryo-EM) of nodavirus RNA replication complexes revealed that the viral double-stranded RNA replication template is coiled inside a 30- to 90-nm invagination of the outer mitochondrial membrane, whose necked aperture to the cytoplasm is gated by a 12-fold symmetric, 35-nm diameter "crown" complex that contains multifunctional viral RNA replication protein A. Here we report optimizing cryo-EM tomography and image processing to improve crown resolution from 33 to 8.5 Å. This resolves the crown into 12 distinct vertical segments, each with 3 major subdomains: A membrane-connected basal lobe and an apical lobe that together comprise the ∼19-nm-diameter central turret, and a leg emerging from the basal lobe that connects to the membrane at ∼35-nm diameter. Despite widely varying replication vesicle diameters, the resulting two rings of membrane interaction sites constrain the vesicle neck to a highly uniform shape. Labeling protein A with a His-tag that binds 5-nm Ni-nanogold allowed cryo-EM tomography mapping of the C terminus of protein A to the apical lobe, which correlates well with the predicted structure of the C-proximal polymerase domain of protein A. These and other results indicate that the crown contains 12 copies of protein A arranged basally to apically in an N-to-C orientation. Moreover, the apical polymerase localization has significant mechanistic implications for template RNA recruitment and (-) and (+)RNA synthesis.

Augmenting subnetwork inference with information extracted from the scientific literature.

Many biological studies involve either (i) manipulating some aspect of a cell or its environment and then simultaneously measuring the effect on thousands of genes, or (ii) systematically manipulating each gene and then measuring the effect on some response of interest. A common challenge that arises in these studies is to explain how genes identified as relevant in the given experiment are organized into a subnetwork that accounts for the response of interest. The task of inferring a subnetwork is typically dependent on the information available in publicly available, structured databases, which suffer from incompleteness. However, a wealth of potentially relevant information resides in the scientific literature, such as information about genes associated with certain concepts of interest, as well as interactions that occur among various biological entities. We contend that by exploiting this information, we can improve the explanatory power and accuracy of subnetwork inference in multiple applications. Here we propose and investigate several ways in which information extracted from the scientific literature can be used to augment subnetwork inference. We show that we can use literature-extracted information to (i) augment the set of entities identified as being relevant in a subnetwork inference task, (ii) augment the set of interactions used in the process, and (iii) support targeted browsing of a large inferred subnetwork by identifying entities and interactions that are closely related to concepts of interest. We use this approach to uncover the pathways involved in interactions between a virus and a host cell, and the pathways that are regulated by a transcription factor associated with breast cancer. Our experimental results demonstrate that these approaches can provide more accurate and more interpretable subnetworks. Integer program code, background network data, and pathfinding code are available at https://github.com/Craven-Biostat-Lab/subnetwork_inference.