International Society for Extracellular Vesicles Workshop. QuantitatEVs: multiscale analyses, from bulk to single extracellular vesicle.

The 'QuantitatEVs: multiscale analyses, from bulk to single vesicle' workshop aimed to discuss quantitative strategies and harmonized wet and computational approaches toward the comprehensive analysis of extracellular vesicles (EVs) from bulk to single vesicle analyses with a special focus on emerging technologies. The workshop covered the key issues in the quantitative analysis of different EV-associated molecular components and EV biophysical features, which are considered the core of EV-associated biomarker discovery and validation for their clinical translation. The in-person-only workshop was held in Trento, Italy, from January 31st to February 2nd, 2023, and continued in Milan on February 3rd with "Next Generation EVs", a satellite event dedicated to early career researchers (ECR). This report summarizes the main topics and outcomes of the workshop.

New advances of the androgen receptor in prostate cancer: report from the 1st International Androgen Receptor Symposium.

The androgen receptor (AR) is a crucial player in various aspects of male reproduction and has been associated with the development and progression of prostate cancer (PCa). Therefore, the protein is the linchpin of current PCa therapies. Despite great research efforts, the AR signaling pathway has still not been deciphered, and the emergence of resistance is still the biggest problem in PCa treatment. To discuss the latest developments in AR research, the "1st International Androgen Receptor Symposium" offered a forum for the exchange of clinical and scientific innovations around the role of the AR in prostate cancer (PCa) and to stimulate new collaborative interactions among leading scientists from basic, translational, and clinical research. The symposium included three sessions covering preclinical studies, prognostic and diagnostic biomarkers, and ongoing prostate cancer clinical trials. In addition, a panel discussion about the future direction of androgen deprivation therapy and anti-AR therapy in PCa was conducted. Therefore, the newest insights and developments in therapeutic strategies and biomarkers are discussed in this report.

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 < 0.0001) tRF-Lys-CTT-5c (p < 0.0327), piR-28004 (p = 0.0081) and miR-375-3p (p < 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.

Viability Analysis and High-Content Live-Cell Imaging for Drug Testing in Prostate Cancer Xenograft-Derived Organoids.

Tumor organoids have been pushed forward as advanced model systems for in vitro oncology drug testing, with the eventual goal to direct personalized cancer treatments. However, drug testing efforts suffer from a large variation in experimental conditions for organoid culturing and organoid treatment. Moreover, most drug tests are restricted to whole-well viability as the sole read-out, thereby losing important information about key biological aspects that might be impacted due to the use of administered drugs. These bulk read-outs also discard potential inter-organoid heterogeneity in drug responses. To tackle these issues, we developed a systematic approach for processing organoids from prostate cancer (PCa) patient-derived xenografts (PDXs) for viability-based drug testing and identified essential conditions and quality checks for consistent results. In addition, we generated an imaging-based drug testing procedure using high-content fluorescence microscopy in living PCa organoids to detect various modalities of cell death. Individual organoids and cell nuclei in organoids were segmented and quantified using a dye combination of Hoechst 33342, propidium iodide and Caspase 3/7 Green, allowing the identification of cytostatic and cytotoxic treatment effects. Our procedures provide important insights into the mechanistic actions of tested drugs. Moreover, these methods can be adapted for tumor organoids originating from other cancer types to increase organoid-based drug test validity, and ultimately, accelerate clinical implementation.

Site-specific analysis of ribosomal 2'O-methylation by quantitative reverse transcription PCR under low deoxynucleotide triphosphate concentrations.

Ribose 2'O-methylation (Nm, ribomethylation) is the most abundant RNA modification present in rRNA. It has been shown that alterations in ribosomal 2'O-methylation at individual Nm sites likely reflect regulated cellular processes. Although several analytical approaches for Nm detection and profiling have been developed, a simple and affordable method for the screening and measurement of individual Nm sites in large numbers of tissue samples is required to examine their potential for clinical translation. Here, we describe a new quantitative reverse transcription PCR-based method that can sensitively assess ribomethylation levels at specific rRNA sites at single-nucleotide resolution in low input amounts of total RNA.

Enhancer profiling identifies epigenetic markers of endocrine resistance and reveals therapeutic options for metastatic castration-resistant prostate cancer patients.

Androgen Receptor (AR) signaling inhibitors, including enzalutamide, are treatment options for patients with metastatic castration-resistant prostate cancer (mCRPC), but resistance inevitably develops. Using metastatic samples from a prospective phase II clinical trial, we epigenetically profiled enhancer/promoter activities with H3K27ac chromatin immunoprecipitation followed by sequencing, before and after AR-targeted therapy. We identified a distinct subset of H3K27ac-differentially marked regions that associated with treatment responsiveness. These data were successfully validated in mCRPC patient-derived xenograft models (PDX). In silico analyses revealed HDAC3 as a critical factor that can drive resistance to hormonal interventions, which we validated in vitro . Using cell lines and mCRPC PDX tumors in vitro , we identified drug-drug synergy between enzalutamide and the pan-HDAC inhibitor vorinostat, providing therapeutic proof-of-concept. These findings demonstrate rationale for new therapeutic strategies using a combination of AR and HDAC inhibitors to improve patient outcome in advanced stages of mCRPC.

Genetic Aspects and Molecular Testing in Prostate Cancer: A Report from a Dutch Multidisciplinary Consensus Meeting.

Background: Germline and tumour genetic testing in prostate cancer (PCa) is becoming more broadly accepted, but testing indications and clinical consequences for carriers in each disease stage are not yet well defined. Objective: To determine the consensus of a Dutch multidisciplinary expert panel on the indication and application of germline and tumour genetic testing in PCa. Design setting and participants: The panel consisted of 39 specialists involved in PCa management. We used a modified Delphi method consisting of two voting rounds and a virtual consensus meeting. Outcome measurements and statistical analysis: Consensus was reached if ≥75% of the panellists chose the same option. Appropriateness was assessed by the RAND/UCLA appropriateness method. Results and limitations: Of the multiple-choice questions, 44% reached consensus. For men without PCa having a relevant family history (familial PCa/BRCA-related hereditary cancer), follow-up by prostate-specific antigen was considered appropriate. For patients with low-risk localised PCa and a family history of PCa, active surveillance was considered appropriate, except in case of the patient being a BRCA2 germline pathogenic variant carrier. Germline and tumour genetic testing should not be done for nonmetastatic hormone-sensitive PCa in the absence of a relevant family history of cancer. Tumour genetic testing was deemed most appropriate for the identification of actionable variants, with uncertainty for germline testing. For tumour genetic testing in metastatic castration-resistant PCa, consensus was not reached for the timing and panel composition. The principal limitations are as follows: (1) a number of topics discussed lack scientific evidence, and therefore the recommendations are partly opinion based, and (2) there was a small number of experts per discipline. Conclusions: The outcomes of this Dutch consensus meeting may provide further guidance on genetic counselling and molecular testing related to PCa. Patient summary: A group of Dutch specialists discussed the use of germline and tumour genetic testing in prostate cancer (PCa) patients, indication of these tests (which patients and when), and impact of these tests on the management and treatment of PCa.

A Compendium of AR Splice Variants in Metastatic Castration-Resistant Prostate Cancer.

Treatment-induced AR alterations, including AR alternative splice variants (AR-Vs), have been extensively linked to harboring roles in primary and acquired resistance to conventional and next-generation hormonal therapies in prostate cancer and therefore have gained momentum. Our aim was to uniformly determine recurrent AR-Vs in metastatic castration-resistant prostate cancer (mCRPC) using whole transcriptome sequencing in order to assess which AR-Vs might hold potential diagnostic or prognostic relevance in future research. This study reports that in addition to the promising AR-V7 as a biomarker, AR45 and AR-V3 were also seen as recurrent AR-Vs and that the presence of any AR-V could be associated with higher AR expression. With future research, these AR-Vs may therefore harbor similar or complementary roles to AR-V7 as predictive and prognostic biomarkers in mCRPC or as proxies for abundant AR expression.

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

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.

Isolation-free measurement of single urinary extracellular vesicles by imaging flow cytometry.

Urinary extracellular vesicles (uEVs) are promising biomarkers for various diseases. However, many tools measuring uEVs rely on time-consuming uEV isolation methods, which could induce sample bias. This study demonstrates the detection of single uEVs without isolation using imaging flow cytometry (IFCM). Unstained urine samples contained auto-fluorescent (A-F) particles when characterized with IFCM. Centrifugation successfully removed A-F particles from the unprocessed urine. Based on the disappearance of A-F particles, a gate was defined to distinguish uEVs from A-F particles. The final readouts of IFCM were verified as single EVs based on detergent treatment and serial dilutions. When developing this protocol to measure urine samples with abnormally high protein levels, 25 mg/mL dithiothreitol (DTT) showed improved uEV recovery over 200 mg/mL DTT. This study provides an isolation-free protocol using IFCM to quantify and phenotype single uEVs, eliminating the hindrance and influence of A-F particles, protein aggregates, and coincidence events.

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.

Identifying genes targeted by disease-associated non-coding SNPs with a protein knowledge graph.

Genome-wide association studies (GWAS) have identified many single nucleotide polymorphisms (SNPs) that play important roles in the genetic heritability of traits and diseases. With most of these SNPs located on the non-coding part of the genome, it is currently assumed that these SNPs influence the expression of nearby genes on the genome. However, identifying which genes are targeted by these disease-associated SNPs remains challenging. In the past, protein knowledge graphs have often been used to identify genes that are associated with disease, also referred to as "disease genes". Here, we explore whether protein knowledge graphs can be used to identify genes that are targeted by disease-associated non-coding SNPs by testing and comparing the performance of six existing methods for a protein knowledge graph, four of which were developed for disease gene identification. We compare our performance against two baselines: (1) an existing state-of-the-art method that is based on guilt-by-association, and (2) the leading assumption that SNPs target the nearest gene on the genome. We test these methods with four reference sets, three of which were obtained by different means. Furthermore, we combine methods to investigate whether their combination improves performance. We find that protein knowledge graphs that include predicate information perform comparable to the current state of the art, achieving an area under the receiver operating characteristic curve (AUC) of 79.6% on average across all four reference sets. Protein knowledge graphs that lack predicate information perform comparable to our other baseline (genetic distance) which achieved an AUC of 75.7% across all four reference sets. Combining multiple methods improved performance to 84.9% AUC. We conclude that methods for a protein knowledge graph can be used to identify which genes are targeted by disease-associated non-coding SNPs.

Cell Line Characteristics Predict Subsequent Resistance to Androgen Receptor-Targeted Agents (ARTA) in Preclinical Models of Prostate Cancer.

Treatment of prostate cancer (PCa) has changed considerably in the last decade due to the introduction of novel androgen receptor (AR)-targeted agents (ARTAs) for patients progressing on androgen deprivation therapy (ADT). Preclinical research however still relies heavily on AR-negative cell line models. In order to investigate potential differences in castration-resistant PCa (CRPC) growth, we set out to create a comprehensive panel of ARTA-progressive models from 4 androgen-responsive AR wild-type PCa cell lines and analyzed its androgen response as opposed to its ADT-progressive counterparts. Parallel cultures of VCaP, DuCaP, PC346C, and LAPC4 were established in their respective culture media with steroid-stripped fetal calf serum (FCS) [dextran-coated charcoal-stripped FCS (DCC)] without androgen (ADT) or in DCC plus 1 µM of the ARTAs bicalutamide, OH-flutamide, or RD162 (an enzalutamide/apalutamide analog). Cell growth was monitored and compared to those of parental cell lines. Short-term androgen response was measured using cell proliferation 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) assay. qRT-PCR was performed to assess the mRNA expression of markers for AR signaling, steroidogenesis, glucocorticoid receptor (GR) signaling, epithelial-mesenchymal transition (EMT), and WNT signaling. Out of 35 parallel cultures per cell line, a total of 24, 15, 34, and 16 CRPC sublines emerged for VCaP, DuCaP, PC346C, and LAPC4, respectively. The addition of bicalutamide or OH-flutamide significantly increased CRPC growth compared to ADT or RD162. VCaP, DuCaP, and PC346C CRPC clones retained an AR-responsive phenotype. The expression of AR and subsequent androgen response were completely lost in all LAPC4 CRPC lines. Markers for EMT and WNT signaling were found to be elevated in the resilient PC346C model and CRPC derivatives of VCaP, DuCaP, and LAPC4. Although the resistant phenotype is pluriform between models, it seems consistent within models, regardless of type of ARTA. These data suggest that the progression to and the phenotype of the CRPC state might already be determined early in carcinogenesis.

CRISPRs in the human genome are differentially expressed between malignant and normal adjacent to tumor tissue.

Clustered Regularly Interspaced Short Palindromic Repeats (CRISPRs) have been identified in bacteria, archaea and mitochondria of plants, but not in eukaryotes. Here, we report the discovery of 12,572 putative CRISPRs randomly distributed across the human chromosomes, which we termed hCRISPRs. By using available transcriptome datasets, we demonstrate that hCRISPRs are distinctively expressed as small non-coding RNAs (sncRNAs) in cell lines and human tissues. Moreover, expression patterns thereof enabled us to distinguish normal from malignant tissues. In prostate cancer, we confirmed the differential hCRISPR expression between normal adjacent and malignant primary prostate tissue by RT-qPCR and demonstrate that the SHERLOCK and DETECTR dipstick tools are suitable to detect these sncRNAs. We anticipate that the discovery of CRISPRs in the human genome can be further exploited for diagnostic purposes in cancer and other medical conditions, which certainly will lead to the development of point-of-care tests based on the differential expression of the hCRISPRs.

Abiraterone switches castration-resistant prostate cancer dependency from adrenal androgens towards androgen receptor variants and glucocorticoid receptor signalling.

INTRODUCTION: Castration-resistant prostate cancer (CRPC) remains dependent on androgen receptor (AR) signalling, which is largely driven by conversion of adrenal androgen precursors lasting after castration. Abiraterone, an inhibitor of the steroidogenic enzyme CYP17A1, has been demonstrated to reduce adrenal androgen synthesis and prolong CRPC patient survival. To study mechanisms of resistance to castration and abiraterone, we created coculture models using human prostate and adrenal tumours. MATERIALS AND METHODS: Castration-naïve and CRPC clones of VCaP were incubated with steroid substrates or cocultured with human adrenal cells (H295R) and treated with abiraterone or the antiandrogen enzalutamide. Male mice bearing VCaP xenografts with and without concurrent H295R xenografts were castrated and treated with placebo or abiraterone. Response was assessed by tumour growth and PSA release. Plasma and tumour steroid levels were assessed by LC/MS-MS. Quantitative polymerase chain reaction determined steroidogenic enzyme, nuclear receptor and AR target gene expression. RESULTS: In vitro, adrenal androgens induced castration-naïve and CRPC cell growth, while precursors steroids for de novo synthesis did not. In a coculture system, abiraterone blocked H295R-induced growth of VCaP cells. In vivo, H295R promoted castration-resistant VCaP growth. Abiraterone only inhibited VCaP growth or PSA production in the presence of H295R. Plasma steroid levels demonstrated CYP17A1 inhibition by abiraterone, whilst CRPC tumour tissue steroid levels showed no evidence of de novo intratumoural androgen production. Castration-resistant and abiraterone-resistant VCaP tumours had increased levels of AR, AR variants and glucocorticoid receptor (GR) resulting in equal AR target gene expression levels compared to noncastrate tumours. CONCLUSIONS: In our model, ligand-dependent AR-regulated regrowth of CRPC was predominantly supported via adrenal androgen precursor production while there was no evidence for intratumoural androgen synthesis. Abiraterone-resistant tumours relied on AR overexpression, expression of ligand-independent AR variants and GR signalling.

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 among its most common genomic mutations. There are conflicting reports about its prognostic role and it remains unclear whether and how it differs in signaling compared with wildtype EGFR. METHODS: To better understand the oncogenic role of EGFRvIII, we leveraged 4 large datasets into 1 large glioblastoma transcriptome dataset (n = 741) alongside 81 whole-genome samples from 2 datasets. RESULTS: The EGFRvIII/EGFR expression ratios differ strongly between tumors and range 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 = .007), which may point toward crosstalk between these pathways. EGFRvIII-expressing tumors 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 toward 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 tumors.

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.

The role of OncoSnoRNAs and Ribosomal RNA 2'-O-methylation in Cancer.

Ribosomes are essential nanomachines responsible for all protein production in cells. Ribosome biogenesis and function are energy costly processes, they are tightly regulated to match cellular needs. In cancer, major pathways that control ribosome biogenesis and function are often deregulated to ensure cell survival and to accommodate the continuous proliferation of tumour cells. Ribosomal RNAs (rRNAs) are abundantly modified with 2'-O-methylation (Nm, ribomethylation) being one of the most common modifications. In eukaryotic ribosomes, ribomethylation is performed by the methyltransferase Fibrillarin guided by box C/D small nucleolar RNAs (snoRNAs). Accumulating evidences indicate that snoRNA expression and ribosome methylation profiles are altered in cancer. Here we review our current knowledge on differential snoRNA expression and rRNA 2'-O methylation in the context of human malignancies, and discuss the consequences and opportunities for cancer diagnostics, prognostics, and therapeutics.

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.