Schlafen 12 restricts HIV-1 latency reversal by a codon-usage dependent post-transcriptional block in CD4+ T cells.

Latency is a major barrier towards virus elimination in HIV-1-infected individuals. Yet, the mechanisms that contribute to the maintenance of HIV-1 latency are incompletely understood. Here we describe the Schlafen 12 protein (SLFN12) as an HIV-1 restriction factor that establishes a post-transcriptional block in HIV-1-infected cells and thereby inhibits HIV-1 replication and virus reactivation from latently infected cells. The inhibitory activity is dependent on the HIV-1 codon usage and on the SLFN12 RNase active sites. Within HIV-1-infected individuals, SLFN12 expression in PBMCs correlated with HIV-1 plasma viral loads and proviral loads suggesting a link with the general activation of the immune system. Using an RNA FISH-Flow HIV-1 reactivation assay, we demonstrate that SLFN12 expression is enriched in infected cells positive for HIV-1 transcripts but negative for HIV-1 proteins. Thus, codon-usage dependent translation inhibition of HIV-1 proteins participates in HIV-1 latency and can restrict the amount of virus release after latency reversal.

CHIKV infection reprograms codon optimality to favor viral RNA translation by altering the tRNA epitranscriptome.

Ample evidence indicates that codon usage bias regulates gene expression. How viruses, such as the emerging mosquito-borne Chikungunya virus (CHIKV), express their genomes at high levels despite an enrichment in rare codons remains a puzzling question. Using ribosome footprinting, we analyze translational changes that occur upon CHIKV infection. We show that CHIKV infection induces codon-specific reprogramming of the host translation machinery to favor the translation of viral RNA genomes over host mRNAs with an otherwise optimal codon usage. This reprogramming was mostly apparent at the endoplasmic reticulum, where CHIKV RNAs show high ribosome occupancy. Mechanistically, it involves CHIKV-induced overexpression of KIAA1456, an enzyme that modifies the wobble U34 position in the anticodon of tRNAs, which is required for proper decoding of codons that are highly enriched in CHIKV RNAs. Our findings demonstrate an unprecedented interplay of viruses with the host tRNA epitranscriptome to adapt the host translation machinery to viral production.

Data-driven identification of inherent features of eukaryotic stress-responsive genes.

Living organisms are continuously challenged by changes in their environment that can propagate to stresses at the cellular level, such as rapid changes in osmolarity or oxygen tension. To survive these sudden changes, cells have developed stress-responsive mechanisms that tune cellular processes. The response of Saccharomyces cerevisiae to osmostress includes a massive reprogramming of gene expression. Identifying the inherent features of stress-responsive genes is of significant interest for understanding the basic principles underlying the rewiring of gene expression upon stress. Here, we generated a comprehensive catalog of osmostress-responsive genes from 5 independent RNA-seq experiments. We explored 30 features of yeast genes and found that 25 (83%) were distinct in osmostress-responsive genes. We then identified 13 non-redundant minimal osmostress gene traits and used statistical modeling to rank the most stress-predictive features. Intriguingly, the most relevant features of osmostress-responsive genes are the number of transcription factors targeting them and gene conservation. Using data on HeLa samples, we showed that the same features that define yeast osmostress-responsive genes can predict osmostress-responsive genes in humans, but with changes in the rank-ordering of feature-importance. Our study provides a holistic understanding of the basic principles of the regulation of stress-responsive gene expression across eukaryotes.

Fusion transcripts and their genomic breakpoints in polyadenylated and ribosomal RNA-minus RNA sequencing data.

BACKGROUND: Fusion genes are typically identified by RNA sequencing (RNA-seq) without elucidating the causal genomic breakpoints. However, non-poly(A)-enriched RNA-seq contains large proportions of intronic reads that also span genomic breakpoints. RESULTS: We have developed an algorithm, Dr. Disco, that searches for fusion transcripts by taking an entire reference genome into account as search space. This includes exons but also introns, intergenic regions, and sequences that do not meet splice junction motifs. Using 1,275 RNA-seq samples, we investigated to what extent genomic breakpoints can be extracted from RNA-seq data and their implications regarding poly(A)-enriched and ribosomal RNA-minus RNA-seq data. Comparison with whole-genome sequencing data revealed that most genomic breakpoints are not, or minimally, transcribed while, in contrast, the genomic breakpoints of all 32 TMPRSS2-ERG-positive tumours were present at RNA level. We also revealed tumours in which the ERG breakpoint was located before ERG, which co-existed with additional deletions and messenger RNA that incorporated intergenic cryptic exons. In breast cancer we identified rearrangement hot spots near CCND1 and in glioma near CDK4 and MDM2 and could directly associate this with increased expression. Furthermore, in all datasets we find fusions to intergenic regions, often spanning multiple cryptic exons that potentially encode neo-antigens. Thus, fusion transcripts other than classical gene-to-gene fusions are prominently present and can be identified using RNA-seq. CONCLUSION: By using the full potential of non-poly(A)-enriched RNA-seq data, sophisticated analysis can reliably identify expressed genomic breakpoints and their transcriptional effects.

Network neighbors of viral targets and differentially expressed genes in COVID-19 are drug target candidates.

The COVID-19 pandemic is raging. It revealed the importance of rapid scientific advancement towards understanding and treating new diseases. To address this challenge, we adapt an explainable artificial intelligence algorithm for data fusion and utilize it on new omics data on viral-host interactions, human protein interactions, and drugs to better understand SARS-CoV-2 infection mechanisms and predict new drug-target interactions for COVID-19. We discover that in the human interactome, the human proteins targeted by SARS-CoV-2 proteins and the genes that are differentially expressed after the infection have common neighbors central in the interactome that may be key to the disease mechanisms. We uncover 185 new drug-target interactions targeting 49 of these key genes and suggest re-purposing of 149 FDA-approved drugs, including drugs targeting VEGF and nitric oxide signaling, whose pathways coincide with the observed COVID-19 symptoms. Our integrative methodology is universal and can enable insight into this and other serious diseases.

LRRC8A-containing chloride channel is crucial for cell volume recovery and survival under hypertonic conditions.

Regulation of cell volume is essential for tissue homeostasis and cell viability. In response to hypertonic stress, cells need rapid electrolyte influx to compensate water loss and to prevent cell death in a process known as regulatory volume increase (RVI). However, the molecular component able to trigger such a process was unknown to date. Using a genome-wide CRISPR/Cas9 screen, we identified LRRC8A, which encodes a chloride channel subunit, as the gene most associated with cell survival under hypertonic conditions. Hypertonicity activates the p38 stress-activated protein kinase pathway and its downstream MSK1 kinase, which phosphorylates and activates LRRC8A. LRRC8A-mediated Cl- efflux facilitates activation of the with-no-lysine (WNK) kinase pathway, which in turn, promotes electrolyte influx via Na+/K+/2Cl- cotransporter (NKCC) and RVI under hypertonic stress. LRRC8A-S217A mutation impairs channel activation by MSK1, resulting in reduced RVI and cell survival. In summary, LRRC8A is key to bidirectional osmotic stress responses and cell survival under hypertonic conditions.

Probabilistic graphlets capture biological function in probabilistic molecular networks.

MOTIVATION: Molecular interactions have been successfully modeled and analyzed as networks, where nodes represent molecules and edges represent the interactions between them. These networks revealed that molecules with similar local network structure also have similar biological functions. The most sensitive measures of network structure are based on graphlets. However, graphlet-based methods thus far are only applicable to unweighted networks, whereas real-world molecular networks may have weighted edges that can represent the probability of an interaction occurring in the cell. This information is commonly discarded when applying thresholds to generate unweighted networks, which may lead to information loss. RESULTS: We introduce probabilistic graphlets as a tool for analyzing the local wiring patterns of probabilistic networks. To assess their performance compared to unweighted graphlets, we generate synthetic networks based on different well-known random network models and edge probability distributions and demonstrate that probabilistic graphlets outperform their unweighted counterparts in distinguishing network structures. Then we model different real-world molecular interaction networks as weighted graphs with probabilities as weights on edges and we analyze them with our new weighted graphlets-based methods. We show that due to their probabilistic nature, probabilistic graphlet-based methods more robustly capture biological information in these data, while simultaneously showing a higher sensitivity to identify condition-specific functions compared to their unweighted graphlet-based method counterparts. AVAILABILITYAND IMPLEMENTATION: Our implementation of probabilistic graphlets is available at https://github.com/Serdobe/Probabilistic_Graphlets. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Host-derived circular RNAs display proviral activities in Hepatitis C virus-infected cells.

Viruses subvert macromolecular pathways in infected host cells to aid in viral gene amplification or to counteract innate immune responses. Roles for host-encoded, noncoding RNAs, including microRNAs, have been found to provide pro- and anti-viral functions. Recently, circular RNAs (circRNAs), that are generated by a nuclear back-splicing mechanism of pre-mRNAs, have been implicated to have roles in DNA virus-infected cells. This study examines the circular RNA landscape in uninfected and hepatitis C virus (HCV)-infected liver cells. Results showed that the abundances of distinct classes of circRNAs were up-regulated or down-regulated in infected cells. Identified circRNAs displayed pro-viral effects. One particular up-regulated circRNA, circPSD3, displayed a very pronounced effect on viral RNA abundances in both hepatitis C virus- and Dengue virus-infected cells. Though circPSD3 has been shown to bind factor eIF4A3 that modulates the cellular nonsense-mediated decay (NMD) pathway, circPSD3 regulates RNA amplification in a pro-viral manner at a post-translational step, while eIF4A3 exhibits the anti-viral property of the NMD pathway. Findings from the global analyses of the circular RNA landscape argue that pro-, and likely, anti-viral functions are executed by circRNAs that modulate viral gene expression as well as host pathways. Because of their long half-lives, circRNAs likely play hitherto unknown, important roles in viral pathogenesis.

A genetic analysis reveals novel histone residues required for transcriptional reprogramming upon stress.

Cells have the ability to sense, respond and adapt to environmental fluctuations. Stress causes a massive reorganization of the transcriptional program. Many examples of histone post-translational modifications (PTMs) have been associated with transcriptional activation or repression under steady-state growth conditions. Comparatively less is known about the role of histone PTMs in the cellular adaptive response to stress. Here, we performed high-throughput genetic screenings that provide a novel global map of the histone residues required for transcriptional reprogramming in response to heat and osmotic stress. Of note, we observed that the histone residues needed depend on the type of gene and/or stress, thereby suggesting a 'personalized', rather than general, subset of histone requirements for each chromatin context. In addition, we identified a number of new residues that unexpectedly serve to regulate transcription. As a proof of concept, we characterized the function of the histone residues H4-S47 and H4-T30 in response to osmotic and heat stress, respectively. Our results uncover novel roles for the kinases Cla4 and Ste20, yeast homologs of the mammalian PAK2 family, and the Ste11 MAPK as regulators of H4-S47 and H4-T30, respectively. This study provides new insights into the role of histone residues in transcriptional regulation under stress conditions.

Functional Network Analysis Reveals the Relevance of SKIIP in the Regulation of Alternative Splicing by p38 SAPK.

Alternative splicing is a prevalent mechanism of gene regulation that is modulated in response to a wide range of extracellular stimuli. Stress-activated protein kinases (SAPKs) play a key role in controlling several steps of mRNA biogenesis. Here, we show that osmostress has an impact on the regulation of alternative splicing (AS), which is partly mediated through the action of p38 SAPK. Splicing network analysis revealed a functional connection between p38 and the spliceosome component SKIIP, whose depletion abolished a significant fraction of p38-mediated AS changes. Importantly, p38 interacted with and directly phosphorylated SKIIP, thereby altering its activity. SKIIP phosphorylation regulated AS of GADD45α, the upstream activator of the p38 pathway, uncovering a negative feedback loop involving AS regulation. Our data reveal mechanisms and targets of SAPK function in stress adaptation through the regulation of AS.

The exonuclease Xrn1 activates transcription and translation of mRNAs encoding membrane proteins.

The highly conserved 5'-3' exonuclease Xrn1 regulates gene expression in eukaryotes by coupling nuclear DNA transcription to cytosolic mRNA decay. By integrating transcriptome-wide analyses of translation with biochemical and functional studies, we demonstrate an unanticipated regulatory role of Xrn1 in protein synthesis. Xrn1 promotes translation of a specific group of transcripts encoding membrane proteins. Xrn1-dependence for translation is linked to poor structural RNA contexts for translation initiation, is mediated by interactions with components of the translation initiation machinery and correlates with an Xrn1-dependence for mRNA localization at the endoplasmic reticulum, the translation compartment of membrane proteins. Importantly, for this group of mRNAs, Xrn1 stimulates transcription, mRNA translation and decay. Our results uncover a crosstalk between the three major stages of gene expression coordinated by Xrn1 to maintain appropriate levels of membrane proteins.

Rapid reversible changes in compartments and local chromatin organization revealed by hyperosmotic shock.

Nuclear architecture is decisive for the assembly of transcriptional responses. However, how chromosome organization is dynamically modulated to permit rapid and transient transcriptional changes in response to environmental challenges remains unclear. Here we show that hyperosmotic stress disrupts different levels of chromosome organization, ranging from A/B compartment changes to reduction in the number and insulation of topologically associating domains (TADs). Concomitantly, transcription is greatly affected, TAD borders weaken, and RNA Polymerase II runs off from hundreds of transcription end sites. Stress alters the binding profiles of architectural proteins, which explains the disappearance of local chromatin organization. These processes are dynamic, and cells rapidly reconstitute their default chromatin conformation after stress removal, uncovering an intrinsic organization. Transcription is not required for local chromatin reorganization, while compartment recovery is partially transcription-dependent. Thus, nuclear organization in mammalian cells can be rapidly modulated by environmental changes in a reversible manner.

Epithelial-Mesenchymal Transition in Human Prostate Cancer Demonstrates Enhanced Immune Evasion Marked by IDO1 Expression.

Cancer invasion and metastasis are driven by epithelial-mesenchymal transition (EMT), yet the exact mechanisms that account for EMT in clinical prostate cancer are not fully understood. Expression of N-cadherin is considered a hallmark of EMT in clinical prostate cancer. In this study, we determined the molecular mechanisms associated with N-cadherin expression in patients with prostate cancer. We performed laser capture microdissection of matched N-cadherin-positive and -negative prostate cancer areas from patient samples (n = 8), followed by RNA sequencing. N-cadherin expression was significantly associated with an immune-regulatory signature including profound upregulation of indoleamine 2,3-dioxygenase (IDO1; log2-fold change = 5.1; P = 2.98E-04). Fluorescent immunostainings of patient samples confirmed expression of IDO1 protein and also its metabolite kynurenine in primarily N-cadherin-positive areas. N-cadherin-positive areas also exhibited a local decrease of intraepithelial cytotoxic (CD8+) T cells and an increase of immunosuppressive regulatory T cells (CD4+/FOXP3+). In conclusion, EMT in clinical prostate cancer is accompanied by upregulated expression of IDO1 and an increased number of regulatory T cells. These data indicate that EMT, which is an important step in tumor progression, can be protected from effective immune control in patients with prostate cancer.Significance: These findings demonstrate EMT is linked to an immunosuppressive environment in clinical prostate cancer, suggesting that patients with prostate cancer can potentially benefit from combinatorial drug therapy. Cancer Res; 78(16); 4671-9. ©2018 AACR.

Cribriform and intraductal prostate cancer are associated with increased genomic instability and distinct genomic alterations.

BACKGROUND: Invasive cribriform and intraductal carcinoma (CR/IDC) is associated with adverse outcome of prostate cancer patients. The aim of this study was to determine the molecular aberrations associated with CR/IDC in primary prostate cancer, focusing on genomic instability and somatic copy number alterations (CNA). METHODS: Whole-slide images of The Cancer Genome Atlas Project (TCGA, N = 260) and the Canadian Prostate Cancer Genome Network (CPC-GENE, N = 199) radical prostatectomy datasets were reviewed for Gleason score (GS) and presence of CR/IDC. Genomic instability was assessed by calculating the percentage of genome altered (PGA). Somatic copy number alterations (CNA) were determined using Fisher-Boschloo tests and logistic regression. Primary analysis were performed on TCGA (N = 260) as discovery and CPC-GENE (N = 199) as validation set. RESULTS: CR/IDC growth was present in 80/260 (31%) TCGA and 76/199 (38%) CPC-GENE cases. Patients with CR/IDC and ≥ GS 7 had significantly higher PGA than men without this pattern in both TCGA (2.2 fold; p = 0.0003) and CPC-GENE (1.7 fold; p = 0.004) cohorts. CR/IDC growth was associated with deletions of 8p, 16q, 10q23, 13q22, 17p13, 21q22, and amplification of 8q24. CNAs comprised a total of 1299 gene deletions and 369 amplifications in the TCGA dataset, of which 474 and 328 events were independently validated, respectively. Several of the affected genes were known to be associated with aggressive prostate cancer such as loss of PTEN, CDH1, BCAR1 and gain of MYC. Point mutations in TP53, SPOP and FOXA1were also associated with CR/IDC, but occurred less frequently than CNAs. CONCLUSIONS: CR/IDC growth is associated with increased genomic instability clustering to genetic regions involved in aggressive prostate cancer. Therefore, CR/IDC is a pathologic substrate for progressive molecular tumour derangement.

Regulation of transcription elongation in response to osmostress.

Cells trigger massive changes in gene expression upon environmental fluctuations. The Hog1 stress-activated protein kinase (SAPK) is an important regulator of the transcriptional activation program that maximizes cell fitness when yeast cells are exposed to osmostress. Besides being associated with transcription factors bound at target promoters to stimulate transcriptional initiation, activated Hog1 behaves as a transcriptional elongation factor that is selective for stress-responsive genes. Here, we provide insights into how this signaling kinase functions in transcription elongation. Hog1 phosphorylates the Spt4 elongation factor at Thr42 and Ser43 and such phosphorylations are essential for the overall transcriptional response upon osmostress. The phosphorylation of Spt4 by Hog1 regulates RNA polymerase II processivity at stress-responsive genes, which is critical for cell survival under high osmostress conditions. Thus, the direct regulation of Spt4 upon environmental insults serves to stimulate RNA Pol II elongation efficiency.

A mononucleotide repeat in PRRT2 is an important, frequent target of mismatch repair deficiency in cancer.

The DNA mismatch repair (MMR) system corrects DNA replication mismatches thereby contributing to the maintenance of genomic stability. MMR deficiency has been observed in prostate cancer but its impact on the genomic landscape of these tumours is not known. In order to identify MMR associated mutations in prostate cancer we have performed whole genome sequencing of the MMR deficient PC346C prostate cancer cell line. We detected a total of 1196 mutations in PC346C which was 1.5-fold higher compared to a MMR proficient prostate cancer sample (G089). Of all different mutation classes, frameshifts in mononucleotide repeat (MNR) sequences were significantly enriched in the PC346C sample. As a result, a selection of genes with frameshift mutations in MNR was further assessed regarding its mutational status in a comprehensive panel of prostate, ovarian, endometrial and colorectal cancer cell lines. We identified PRRT2 and DAB2IP to be frequently mutated in MMR deficient cell lines, colorectal and endometrial cancer patient samples. Further characterization of PRRT2 revealed an important role of this gene in cancer biology. Both normal prostate cell lines and a colorectal cancer cell line showed increased proliferation, migration and invasion when expressing the mutated form of PRRT2 (ΔPRRT2). The wild-type PRRT2 (PRRT2wt) had an inhibitory effect in proliferation, consistent with the low expression level of PRRT2 in cancer versus normal prostate samples.

Evolution of protein phosphorylation across 18 fungal species.

Living organisms have evolved protein phosphorylation, a rapid and versatile mechanism that drives signaling and regulates protein function. We report the phosphoproteomes of 18 fungal species and a phylogenetic-based approach to study phosphosite evolution. We observe rapid divergence, with only a small fraction of phosphosites conserved over hundreds of millions of years. Relative to recently acquired phosphosites, ancient sites are enriched at protein interfaces and are more likely to be functionally important, as we show for sites on H2A1 and eIF4E. We also observe a change in phosphorylation motif frequencies and kinase activities that coincides with the whole-genome duplication event. Our results provide an evolutionary history for phosphosites and suggest that rapid evolution of phosphorylation can contribute strongly to phenotypic diversity.

Human PDE4D isoform composition is deregulated in primary prostate cancer and indicative for disease progression and development of distant metastases.

Phosphodiesterase 4D7 was recently shown to be specifically over-expressed in localized prostate cancer, raising the question as to which regulatory mechanisms are involved and whether other isoforms of this gene family (PDE4D) are affected under the same conditions.We investigated PDE4D isoform composition in prostatic tissues using a total of seven independent expression datasets and also included data on DNA methylation, copy number and AR and ERG binding in PDE4D promoters to gain insight into their effect on PDE4D transcription.We show that expression of PDE4D isoforms is consistently altered in primary human prostate cancer compared to benign tissue, with PDE4D7 being up-regulated while PDE4D5 and PDE4D9 are down-regulated. Disease progression is marked by an overall down-regulation of long PDE4D isoforms, while short isoforms (PDE4D1/2) appear to be relatively unaffected. While these alterations seem to be independent of copy number alterations in the PDE4D locus and driven by AR and ERG binding, we also observed increased DNA methylation in the promoter region of PDE4D5, indicating a long lasting alteration of the isoform composition in prostate cancer tissues.We propose two independent metrics that may serve as diagnostic and prognostic markers for prostate disease: (PDE4D7 - PDE4D5) provides an effective means for distinguishing PCa from normal adjacent prostate, whereas PDE4D1/2 - (PDE4D5 + PDE4D7 + PDE4D9) offers strong prognostic potential to detect aggressive forms of PCa and is associated with metastasis free survival. Overall, our findings highlight the relevance of PDE4D as prostate cancer biomarker and potential drug target.

Loss of SLCO1B3 drives taxane resistance in prostate cancer.

BACKGROUND: Both taxanes, docetaxel and cabazitaxel, are effective treatments for metastatic castration-resistant prostate cancer (mCRPC). However, resistance to taxanes is common. Our objective was to investigate mechanisms of taxane resistance in prostate cancer. METHODS: Two docetaxel-resistant patient-derived xenografts (PDXs) of CRPC were established (PC339-DOC and PC346C-DOC) in male athymic nude mice by frequent intraperitoneal administrations of docetaxel. Next-generation sequencing was performed on PDX tissue pre- and post-docetaxel resistance and gene expression profiles were compared. [(14)C]-docetaxel and [(14)C]-cabazitaxel uptake assays in vitro and cytotoxicity assays were performed to validate direct involvement of transporter genes in taxane sensitivity. RESULTS: Organic anion-transporting polypeptide (SLCO1B3), an influx transporter of docetaxel, was significantly downregulated in PC346C-DOC tumours. In accordance with this finding, intratumoural concentrations of docetaxel and cabazitaxel were significantly decreased in PC346C-DOC as compared with levels in chemotherapy-naive PC346C tumours. In addition, silencing of SLCO1B3 in chemo-naive PC346C resulted in a two-fold decrease in intracellular concentrations of both taxanes. Overexpression of SLCO1B3 showed higher sensitivity to docetaxel and cabazitaxel. CONCLUSIONS: The SLCO1B3 determines intracellular concentrations of docetaxel and cabazitaxel and consequently influences taxane efficacy. Loss of the drug transporter SLCO1B3 may drive taxane resistance in prostate cancer.

Secreted Phospholipases A2 Are Intestinal Stem Cell Niche Factors with Distinct Roles in Homeostasis, Inflammation, and Cancer.

The intestinal stem cell niche provides cues that actively maintain gut homeostasis. Dysregulation of these cues may compromise intestinal regeneration upon tissue insult and/or promote tumor growth. Here, we identify secreted phospholipases A2 (sPLA2s) as stem cell niche factors with context-dependent functions in the digestive tract. We show that group IIA sPLA2, a known genetic modifier of mouse intestinal tumorigenesis, is expressed by Paneth cells in the small intestine, while group X sPLA2 is expressed by Paneth/goblet-like cells in the colon. During homeostasis, group IIA/X sPLA2s inhibit Wnt signaling through intracellular activation of Yap1. However, upon inflammation they are secreted into the intestinal lumen, where they promote prostaglandin synthesis and Wnt signaling. Genetic ablation of both sPLA2s improves recovery from inflammation but increases colon cancer susceptibility due to release of their homeostatic Wnt-inhibitory role. This "trade-off" effect suggests sPLA2s have important functions as genetic modifiers of inflammation and colon cancer.