FASTAFS:file system virtualisation of random access compressed FASTA files

Background: The FASTA file format, used to store polymeric sequence data, has become a bioinformatics file standard used for decades. The relatively large files require additional files, beyond the scope of the original format, to identify sequences and to provide random access. Multiple compressors have been developed to archive FASTA files back and forth, but these lack direct access to targeted content or metadata of the archive. Moreover, these solutions are not directly backwards compatible to FASTA files, resulting in limited software integration. Results: We designed a linux based toolkit that virtualises the content of DNA, RNA and protein FASTA archives into the filesystem by using filesystem in userspace. This guarantees in-sync virtualised metadata files and offers fast random-access decompression using bit encodings plus Zstandard (zstd). The toolkit, FASTAFS, can track all its system-wide running instances, allows file integrity verification and can provide, instantly, scriptable access to sequence files and is easy to use and deploy. The file compression ratios were comparable but not superior to other state of the art archival tools, despite the innovative random access feature implemented in FASTAFS. Conclusions: FASTAFS is a user-friendly and easy to deploy backwards compatible generic purpose solution to store and access compressed FASTA files, since it offers file system access to FASTA files as well as in-sync metadata files through file virtualisation. Using virtual filesystems as in-between layer offers format conversion without the need to rewrite code into different programming languages while preserving compatibility.</p

Validation of potential RNA biomarkers for prostate cancer diagnosis and monitoring in plasma and urinary extracellular vesicles

Introduction: Prostate cancer (PCa), one of the most prevalent malignancies affecting men worldwide, presents significant challenges in terms of early detection, risk stratification, and active surveillance. In recent years, liquid biopsies have emerged as a promising non-invasive approach to complement or even replace traditional tissue biopsies. Extracellular vesicles (EVs), nanosized membranous structures released by various cells into body fluids, have gained substantial attention as a source of cancer biomarkers due to their ability to encapsulate and transport a wide range of biological molecules, including RNA. In this study, we aimed to validate 15 potential RNA biomarkers, identified in a previous EV RNA sequencing study, using droplet digital PCR. Methods: The candidate biomarkers were tested in plasma and urinary EVs collected before and after radical prostatectomy from 30 PCa patients and their diagnostic potential was evaluated in a test cohort consisting of 20 benign prostate hyperplasia (BPH) and 20 PCa patients’ plasma and urinary EVs. Next, the results were validated in an independent cohort of plasma EVs from 31 PCa and 31 BPH patients. Results: We found that the levels of NKX3-1 (p = 0.0008) in plasma EVs, and tRF-Phe-GAA-3b (p &lt; 0.0001) tRF-Lys-CTT-5c (p &lt; 0.0327), piR-28004 (p = 0.0081) and miR-375-3p (p &lt; 0.0001) in urinary EVs significantly decreased after radical prostatectomy suggesting that the main tissue source of these RNAs is prostate and/or PCa. Two mRNA biomarkers—GLO1 and NKX3-1 showed promising diagnostic potential in distinguishing between PCa and BPH with AUC of 0.68 and 0.82, respectively, in the test cohort and AUC of 0.73 and 0.65, respectively, in the validation cohort, when tested in plasma EVs. Combining these markers in a biomarker model yielded AUC of 0.85 and 0.71 in the test and validation cohorts, respectively. Although the PSA levels in the blood could not distinguish PCa from BPH in our cohort, adding PSA to the mRNA biomarker model increased AUC from 0.71 to 0.76. Conclusion: This study identified two novel EV-enclosed RNA biomarkers–NKX3-1 and GLO1–for the detection of PCa, and highlights the complementary nature of GLO1, NKX3-1 and PSA as combined biomarkers in liquid biopsies of PCa.</p

Androgen receptor profiling predicts prostate cancer outcome

Prostate cancer is the second most prevalent malignancy in men. Biomarkers for outcome prediction are urgently needed, so that high-risk patients could be monitored more closely postoperatively. To identify prognostic markers and to determine causal players in prostate cancer progression, we assessed changes in chromatin state during tumor development and progression. Based on this, we assessed genomewide androgen receptor/chromatin binding and identified a distinct androgen receptor/chromatin binding profile between primary prostate cancers and tumors with an acquired resistance to therapy. These differential androgen receptor/chromatin interactions dictated expression of a distinct gene signature with strong prognostic potential. Further refinement of the signature provided us with a concise list of nine genes that hallmark prostate cancer outcome in multiple independent validation series. In this report, we identified a novel gene expression signature for prostate cancer outcome through generation of multilevel genomic data on chromatin accessibility and transcriptional regulation and integration with publically available transcriptomic and clinical datastreams. By combining existing technologies, we propose a novel pipeline for biomarker discovery that is easily implementable in other fields of oncology

Patient-Derived Xenografts and Organoids Recapitulate Castration-Resistant Prostate Cancer with Sustained Androgen Receptor Signaling

Castration-resistant prostate cancer (CRPC) remains an incurable and lethal malignancy. The development of new CRPC treatment strategies is strongly impeded by the scarcity of representative, scalable and transferable preclinical models of advanced, androgen receptor (AR)-driven CRPC. Here, we present contemporary patient-derived xenografts (PDXs) and matching PDX-derived organoids (PDXOs) from CRPC patients who had undergone multiple lines of treatment. These models were comprehensively profiled at the morphologic, genomic ( n = 8) and transcriptomic levels ( n = 81). All are high-grade adenocarcinomas that exhibit copy number alterations and transcriptomic features representative of CRPC patient cohorts. We identified losses of PTEN and RB1, MYC amplifications, as well as genomic alterations in TP53 and in members of clinically actionable pathways such as AR, PI3K and DNA repair pathways. Importantly, the clinically observed continued reliance of CRPC tumors on AR signaling is preserved across the entire set of models, with AR amplification identified in four PDXs. We demonstrate that PDXs and PDXOs faithfully reflect donor tumors and mimic matching patient drug responses. In particular, our models predicted patient responses to subsequent treatments and captured sensitivities to previously received therapies. Collectively, these PDX-PDXO pairs constitute a reliable new resource for in-depth studies of treatment-induced, AR-driven resistance mechanisms. Moreover, PDXOs can be leveraged for large-scale tumor-specific drug response profiling critical for accelerating therapeutic advances in CRPC. </p

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

AR splice variants in circulating tumor cells of patients with castration-resistant prostate cancer: relation with outcome to cabazitaxel

The androgen receptor splice variant (AR-V) 7 in circulating tumor cells (CTCs) is a predictor for resistance to anti-AR-targeted treatment, but not to taxane-based chemotherapy in metastatic castration-resistant prostate cancer (mCRPC). In this study, we investigated whether the presence of two constitutively active variants (AR-V3, AR-V7) and two other conditionally activated variants (AR-V1, AR-V9) vs full-length androgen receptor (AR-FL) measured in CTCs from patients with mCRPC were associated with outcome to therapy with the taxane cabazitaxel. Blood was collected at baseline and after two cycles of cabazitaxel from 118 mCRPC patients starting cabazitaxel in a prospective phase II trial. CellSearch-enriched CTCs were enumerated and in parallel characterized for the presence of the AR-Vs by reverse transcription quantitative polymerase chain reaction. Correlations with CTC and prostate-specific antigen response to cabazitaxel as well as associations with overall survival (OS) were investigated. All AR-Vs were frequently pre

Modulation of Androgen Receptor Signaling in Hormonal Therapy-Resistant Prostate Cancer Cell Lines

Background: Prostate epithelial cells depend on androgens for survival and function. In (early) prostate cancer (PCa) androgens also regulate tumor growth, which is exploited by hormonal therapies in metastatic disease. The aim of the present study was to characterize the androgen receptor (AR) response in hormonal therapy-resistant PC346 cells and identify potential disease markers. Methodology/Principal Findings: Human 19K oligoarrays were used to establish the androgen-regulated expression profile of androgen-responsive PC346C cells and its derivative therapy-resistant sublines: PC346DCC (vestigial AR levels), PC346Flu1 (AR overexpression) and PC346Flu2 (T877A AR mutation). In total, 107 transcripts were differentially-expressed in PC346C and derivatives after R1881 or hydroxyflutamide stimulations. The AR-regulated expression profiles reflected the AR modifications of respective therapy-resistant sublines: AR overexpression resulted in stronger and broader transcriptional response to R1881 stimulation, AR down-regulation correlated with deficient response of AR-target genes and the T877A mutation resulted in transcriptional response to both R1881 and hydroxyflutamide. This AR-target signature was linked to multiple publicly available cell line and tumor derived PCa databases, revealing that distinct functional clusters were differentially modulated during PCa progression. Differentiation and secretory functions were up-regulated in primary PCa but repressed i

Bypass Mechanisms of the Androgen Receptor Pathway in Therapy-Resistant Prostate Cancer Cell Models

Background: Prostate cancer is initially dependent on androgens for survival and growth, making hormonal therapy the cornerstone treatment for late-stage tumors. However, despite initial remission, the cancer will inevitably recur. The present study was designed to investigate how androgen-dependent prostate cancer cells eventually survive and resume growth under androgen-deprived and antiandrogen supplemented conditions. As model system, we used the androgen-responsive PC346C cell line and its therapy-resistant sublines: PC346DCC, PC346Flu1 and PC346Flu2. Methodology/Principal Findings: Microarray technology was used to analyze differences in gene expression between the androgen-responsive and therapy-resistant PC346 cell lines. Microarray analysis revealed 487 transcripts differentiallyexpressed between the androgen-responsive and the therapy-resistant cell lines. Most of these genes were common to all three therapy-resistant sublines and only a minority (,5%) was androgen-regulated. Pathway analysis revealed enrichment in functions involving cellular movement, cell growth and cell death, as well as association with cancer and reproductive system disease. PC346DCC expressed residual levels of androgen receptor (AR) and showed significant down-regulation of androgen-regulated genes (p-value = 10 27). Up-regulation of VAV3 and TWIST1 oncogenes and repression of the DKK3 tumor-suppressor was observed in PC346DCC, suggesting a potential AR bypass mechanism. Subsequent validation of these three genes in patient samples confirmed that expression was deregulated during prostate cancer progression

Dynamic release of neuronal extracellular vesicles containing miR-21a-5p is induced by hypoxia

© 2022 The Authors. Journal of Extracellular Vesicles published by Wiley Periodicals, LLC on behalf of the International Society for Extracellular Vesicles. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.Hypoxia induces changes in the secretion of extracellular vesicles (EVs) in several non-neuronal cells and pathological conditions. EVs are packed with biomolecules, such as microRNA(miR)-21-5p, which respond to hypoxia. However, the true EV association of miR-21-5p, and its functional or biomarker relevance, are inadequately characterised. Neurons are extremely sensitive cells, and it is not known whether the secretion of neuronal EVs and miR-21-5p are altered upon hypoxia. Here, we characterised the temporal EV secretion profile and cell viability of neurons under hypoxia. Hypoxia induced a rapid increase of miR-21a-5p secretion in the EVs, which preceded the elevation of hypoxia-induced tissue or cellular miR-21a-5p. Prolonged hypoxia induced cell death and the release of morphologically distinct EVs. The EVs protected miR-21a-5p from enzymatic degradation but a remarkable fraction of miR-21a-5p remained fragile and non-EV associated. The increase in miR-21a-5p secretion may have biomarker potential, as high blood levels of miR-21-5p in stroke patients were associated with significant disability at hospital discharge. Our data provides an understanding of the dynamic regulation of EV secretion from neurons under hypoxia and provides a candidate for the prediction of recovery from ischemic stroke.This work was supported by the University of Eastern Finland, Emil Aaltonen Foundation, Paavo Nurmi Foundation, Saastamoinen Foundation, Instrumentarium Science Foundation and Business Finland (Grant number 4399/31/2019). Work with clinical samples was supported by the European Regional Development Fund - Project INBIO (No. CZ.02.1.01/0.0/0.0/16_026/0008451). Work with EVQuant was supported by the IMMPROVE Alpe d'HuZes grant of the Dutch Cancer Society (EMCR2015-8022) and the Daniel den Hoed Foundation grant for Erasmus MC Cancer Treatment Screening Facility. L.R. is supported by a predoctoral fellowship grant (IFI17/00012) and J.M. is the principal investigator of the grant PI18/804 ‘MULTI-BIO-TARGETS: a new strategy for stroke management combining outcome biomarkers and neuroprotection’, both from the Instituto de Salud Carlos III.Peer reviewe

The EGFRvIII transcriptome in glioblastoma, a meta-omics analysis.

BACKGROUND: EGFR is among the genes most frequently altered in glioblastoma, with exons 2-7 deletions (EGFRvIII) being amongst its most common genomic mutations. There are conflicting reports about its prognostic role and it remains unclear whether and how it differs in signalling compared with wildtype EGFR. METHODS: To better understand the oncogenic role of EGFRvIII, we leveraged four large datasets into one large glioblastoma transcriptome dataset (n=741) alongside 81 whole-genome samples from two datasets. RESULTS: The EGFRvIII/EGFR expression ratios differ strongly between tumours and ranges from 1% to 95%. Interestingly, the slope of relative EGFRvIII expression is near-linear, which argues against a more positive selection pressure than EGFR wildtype. An absence of selection pressure is also suggested by the similar survival between EGFRvIII positive and negative glioblastoma patients. EGFRvIII levels are inversely correlated with pan-EGFR (all wildtype and mutant variants) expression, which indicates that EGFRvIII has a higher potency in downstream pathway activation. EGFRvIII-positive glioblastomas have a lower CDK4 or MDM2 amplification incidence than EGFRvIII-negative (p=0.007), which may point towards crosstalk between these pathways. EGFRvIII-expressing tumours have an upregulation of 'classical' subtype genes compared to those with EGFR-amplification only (p=3.873e-6). Genomic breakpoints of the EGFRvIII deletions have a preference towards the 3' end of the large intron-1. These preferred breakpoints preserve a cryptic exon resulting in a novel EGFRvIII variant and preserve an intronic enhancer. CONCLUSIONS: These data provide deeper insights into the complex EGFRvIII biology and provide new insights for targeting EGFRvIII mutated tumours