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

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

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

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

Consensus molecular subtype classification of colorectal adenomas

Consensus molecular subtyping is an RNA expression-based classification system for colorectal cancer (CRC). Genomic alterations accumulate during CRC pathogenesis, including the premalignant adenoma stage, leading to changes in RNA expression. Only a minority of adenomas progress to malignancies, a transition that is associated with specific DNA copy number aberrations or microsatellite instability (MSI). We aimed to investigate whether colorectal adenomas can already be stratified into consensus molecular subtype (CMS) classes, and whether specific CMS classes are related to the presence of specific DNA copy number aberrations associated with progression to malignancy. RNA sequencing was performed on 62 adenomas and 59 CRCs. MSI status was determined with polymerase chain reaction-based methodology. DNA copy number was assessed by low-coverage DNA sequencing (n = 30) or array-comparative genomic hybridisation (n = 32). Adenomas were classified into CMS classes together with CRCs from the study cohort and from The Cancer Genome Atlas (n = 556), by use of the established CMS classifier. As a result, 54 of 62 (87%) adenomas were classified according to the CMS. The CMS3 ‘metabolic subtype’, which was least common among CRCs, was most prevalent among adenomas (n = 45; 73%). One of the two adenomas showing MSI was classified as CMS1 (2%), the ‘MSI immune’ subtype. Eight adenomas (13%) were classified as the ‘canonical’ CMS2. No adenomas were classified as the ‘mesenchymal’ CMS4, consistent with the fact that adenomas lack invasion-associated stroma. The distribution of the CMS classes among adenomas was confirmed in an independent series. CMS3 was enriched with adenomas at low risk of progressing to CRC, whereas relatively more high-risk adenomas were observed in CMS2. We conclude that adenomas can be stratified into the CMS classes. Considering that CMS1 and CMS2 expression signatures may mark adenomas at increased risk of progression, the distribution of the CMS classes among adenomas is consistent with the proportion of adenomas expected to progress to CRC

Trans-ancestry genome-wide association meta-analysis of prostate cancer identifies new susceptibility loci and informs genetic risk prediction.

Prostate cancer is a highly heritable disease with large disparities in incidence rates across ancestry populations. We conducted a multiancestry meta-analysis of prostate cancer genome-wide association studies (107,247 cases and 127,006 controls) and identified 86 new genetic risk variants independently associated with prostate cancer risk, bringing the total to 269 known risk variants. The top genetic risk score (GRS) decile was associated with odds ratios that ranged from 5.06 (95% confidence interval (CI), 4.84-5.29) for men of European ancestry to 3.74 (95% CI, 3.36-4.17) for men of African ancestry. Men of African ancestry were estimated to have a mean GRS that was 2.18-times higher (95% CI, 2.14-2.22), and men of East Asian ancestry 0.73-times lower (95% CI, 0.71-0.76), than men of European ancestry. These findings support the role of germline variation contributing to population differences in prostate cancer risk, with the GRS offering an approach for personalized risk prediction

Minimal information for studies of extracellular vesicles 2018 (MISEV2018):a position statement of the International Society for Extracellular Vesicles and update of the MISEV2014 guidelines

The last decade has seen a sharp increase in the number of scientific publications describing physiological and pathological functions of extracellular vesicles (EVs), a collective term covering various subtypes of cell-released, membranous structures, called exosomes, microvesicles, microparticles, ectosomes, oncosomes, apoptotic bodies, and many other names. However, specific issues arise when working with these entities, whose size and amount often make them difficult to obtain as relatively pure preparations, and to characterize properly. The International Society for Extracellular Vesicles (ISEV) proposed Minimal Information for Studies of Extracellular Vesicles (“MISEV”) guidelines for the field in 2014. We now update these “MISEV2014” guidelines based on evolution of the collective knowledge in the last four years. An important point to consider is that ascribing a specific function to EVs in general, or to subtypes of EVs, requires reporting of specific information beyond mere description of function in a crude, potentially contaminated, and heterogeneous preparation. For example, claims that exosomes are endowed with exquisite and specific activities remain difficult to support experimentally, given our still limited knowledge of their specific molecular machineries of biogenesis and release, as compared with other biophysically similar EVs. The MISEV2018 guidelines include tables and outlines of suggested protocols and steps to follow to document specific EV-associated functional activities. Finally, a checklist is provided with summaries of key points