NORMSEQ: a tool for evaluation, selection and visualization of RNA-Seq normalization methods.

RNA-sequencing has become one of the most used high-throughput approaches to gain knowledge about the expression of all different RNA subpopulations. However, technical artifacts, either introduced during library preparation and/or data analysis, can influence the detected RNA expression levels. A critical step, especially in large and low input datasets or studies, is data normalization, which aims at eliminating the variability in data that is not related to biology. Many normalization methods have been developed, each of them relying on different assumptions, making the selection of the appropriate normalization strategy key to preserve biological information. To address this, we developed NormSeq, a free web-server tool to systematically assess the performance of normalization methods in a given dataset. A key feature of NormSeq is the implementation of information gain to guide the selection of the best normalization method, which is crucial to eliminate or at least reduce non-biological variability. Altogether, NormSeq provides an easy-to-use platform to explore different aspects of gene expression data with a special focus on data normalization to help researchers, even without bioinformatics expertise, to obtain reliable biological inference from their data. NormSeq is freely available at: https://arn.ugr.es/normSeq.

sRNAbench and sRNAtoolbox 2022 update: accurate miRNA and sncRNA profiling for model and non-model organisms.

The NCBI Sequence Read Archive currently hosts microRNA sequencing data for over 800 different species, evidencing the existence of a broad taxonomic distribution in the field of small RNA research. Simultaneously, the number of samples per miRNA-seq study continues to increase resulting in a vast amount of data that requires accurate, fast and user-friendly analysis methods. Since the previous release of sRNAtoolbox in 2019, 55 000 sRNAbench jobs have been submitted which has motivated many improvements in its usability and the scope of the underlying annotation database. With this update, users can upload an unlimited number of samples or import them from Google Drive, Dropbox or URLs. Micro- and small RNA profiling can now be carried out using high-confidence Metazoan and plant specific databases, MirGeneDB and PmiREN respectively, together with genome assemblies and libraries from 441 Ensembl species. The new results page includes straightforward sample annotation to allow downstream differential expression analysis with sRNAde. Unassigned reads can also be explored by means of a new tool that performs mapping to microbial references, which can reveal contamination events or biologically meaningful findings as we describe in the example. sRNAtoolbox is available at: https://arn.ugr.es/srnatoolbox/.

Epstein-Barr virus mRNA profiles and viral DNA methylation status in nasopharyngeal brushings from nasopharyngeal carcinoma patients reflect tumor origin.

Undifferentiated nasopharyngeal carcinoma (NPC) is 100% associated with Epstein-Barr virus (EBV) as oncogenic driver. NPC is often diagnosed late due to initial vague complaints and obscured location. Prior studies suggest that measurement of EBV DNA load and RNA transcripts in nasopharyngeal (NP) brushings is useful for minimally invasive NPC diagnosis. However, whether these EBV markers relate to local virus replication or reflect tumor origin remains to be demonstrated. To resolve this, we analysed EBV-DNA characteristics and quantified latent and lytic viral RNA transcripts in NP brushings and matching frozen NP-biopsy specimens from patients suspected of having NPC. We observed non-fragmented and Cp-promotor methylated EBV-DNA in both NP brushings and biopsies suggestive of tumor origin. Using quantitative RT-PCR we determined expression levels of 7 critical latent (EBER1, Qp-EBNA1, EBNA2, BART, LMP1, LMP2, BARF1) and 5 lytic (Zta, Rta, TK, PK and VCA-p18) RNA transcripts. Although latent and early lytic RNA transcripts were frequently detected in conjunction with high EBV viral load, in both brushings and biopsies the latent transcripts prevailed and reflected a typical NPC-associated latency-II transcription profile without EBNA2. Late lytic RNA transcripts were rare and detected at low levels mainly in NP brushings, suggestive of abortive viral reactivation rather than complete virus replication. EBV-IgA serology (EBNA1, VCA, Zta) did not correlate to the level of viral reactivation in situ. Overall, viral RNA profiling, DNA fragmentation and methylation analysis in NP brushings and parallel biopsies indicate that NP brush sampling provides a true and robust indicator of NPC tumor presence.

Multi-modal cell-free DNA genomic and fragmentomic patterns enhance cancer survival and recurrence analysis.

The structure of cell-free DNA (cfDNA) is altered in the blood of patients with cancer. From whole-genome sequencing, we retrieve the cfDNA fragment-end composition using a new software (FrEIA [fragment end integrated analysis]), as well as the cfDNA size and tumor fraction in three independent cohorts (n = 925 cancer from >10 types and 321 control samples). At 95% specificity, we detect 72% cancer samples using at least one cfDNA measure, including 64% early-stage cancer (n = 220). cfDNA detection correlates with a shorter overall (p = 0.0086) and recurrence-free (p = 0.017) survival in patients with resectable esophageal adenocarcinoma. Integrating cfDNA measures with machine learning in an independent test set (n = 396 cancer, 90 controls) achieve a detection accuracy of 82% and area under the receiver operating characteristic curve of 0.96. In conclusion, harnessing the biological features of cfDNA can improve, at no extra cost, the diagnostic performance of liquid biopsies.

A novel persistence associated EBV miRNA expression profile is disrupted in neoplasia.

We have performed the first extensive profiling of Epstein-Barr virus (EBV) miRNAs on in vivo derived normal and neoplastic infected tissues. We describe a unique pattern of viral miRNA expression by normal infected cells in vivo expressing restricted viral latency programs (germinal center: Latency II and memory B: Latency I/0). This includes the complete absence of 15 of the 34 miRNAs profiled. These consist of 12 BART miRNAs (including approximately half of Cluster 2) and 3 of the 4 BHRF1 miRNAs. All but 2 of these absent miRNAs become expressed during EBV driven growth (Latency III). Furthermore, EBV driven growth is accompanied by a 5-10 fold down regulation in the level of the BART miRNAs expressed in germinal center and memory B cells. Therefore, Latency III also expresses a unique pattern of viral miRNAs. We refer to the miRNAs that are specifically expressed in EBV driven growth as the Latency III associated miRNAs. In EBV associated tumors that employ Latency I or II (Burkitt's lymphoma, Hodgkin's disease, nasopharyngeal carcinoma and gastric carcinoma), the Latency III associated BART but not BHRF1 miRNAs are up regulated. Thus BART miRNA expression is deregulated in the EBV associated tumors. This is the first demonstration that Latency III specific genes (the Latency III associated BARTs) can be expressed in these tumors. The EBV associated tumors demonstrate very similar patterns of miRNA expression yet were readily distinguished when the expression data were analyzed either by heat-map/clustering or principal component analysis. Systematic analysis revealed that the information distinguishing the tumor types was redundant and distributed across all the miRNAs. This resembles "secret sharing" algorithms where information can be distributed among a large number of recipients in such a way that any combination of a small number of recipients is able to understand the message. Biologically, this may be a consequence of functional redundancy between the miRNAs.

Comprehensive profiling of Epstein-Barr virus microRNAs in nasopharyngeal carcinoma.

Epstein-Barr Virus (EBV) establishes a long-term latent infection and is associated with a number of human malignancies that are thought to arise from deregulation of different stages of the viral life cycle. Recently, a large number of microRNAs (miRNAs) have been described for EBV, and it has been suggested that their expression may vary between the different latency states found in normal and malignant tissue. To date, however, no technique has been utilized to comprehensively and quantitatively test this idea by profiling expression of the EBV miRNAs in primary infected tissues. We describe here a multiplex reverse transcription-PCR assay that allows the profiling of 39 of the 40 known mature EBV miRNAs from as little as 250 ng of RNA. With this approach, we present a comprehensive profile of EBV miRNAs in primary nasopharyngeal carcinoma (NPC) tumors including estimates of miRNA copy number per tumor cell. This is the first comprehensive profiling of EBV miRNAs in any EBV-associated tumor. In contrast to previous suggestions, we show that the BART-derived miRNAs are present in a wide range of copy numbers from < or =10(3) per cell in both primary tumors and the widely used NPC-derived C666-1 cell line. However, we confirm the hypothesis that the BHRF1 miRNAs are not expressed in NPC. Lastly, we demonstrate that EBV miRNA expression in the widely used NPC line C666-1 is, with some caveats, broadly representative of primary NPC tumors.

Glycan modification of glioblastoma-derived extracellular vesicles enhances receptor-mediated targeting of dendritic cells.

Glioblastoma is the most prevalent and aggressive primary brain tumour for which total tumour lysate-pulsed dendritic cell vaccination is currently under clinical evaluation. Glioblastoma extracellular vesicles (EVs) may represent an enriched cell-free source of tumour-associated (neo-) antigens to pulse dendritic cells (DCs) for the initiation of an anti-tumour immune response. Capture and uptake of EVs by DCs could occur in a receptor-mediated and presumably glycan-dependent way, yet the glycan composition of glioblastoma EVs is unknown. Here, we set out to characterize the glycocalyx composition of glioblastoma EVs by lectin-binding ELISA and comprehensive immunogold transmission electron microscopy (immuno-TEM). The surface glycan profile of human glioblastoma cell line-derived EVs (50-200 nm) was dominated by α-2,3- and α-2,6 linked sialic acid-capped complex N-glycans and bi-antennary N-glycans. Since sialic acids can trigger immune inhibitory sialic acid-binding Ig-like lectin (Siglec) receptors, we screened for Siglec ligands on the EVs. Glioblastoma EVs showed significant binding to Siglec-9, which is highly expressed on DCs. Surprisingly, however, glioblastoma EVs lack glycans that could bind Dendritic Cell-Specific Intercellular adhesion molecule-3-Grabbing Non-integrin (DC-SIGN, CD209), a receptor that mediates uptake and induction of CD4+ and CD8+ T cell activation. Therefore, we explored whether modification of the EV glycan surface could reduce immune inhibitory Siglec binding, while enhancing EV internalization by DCs in a DC-SIGN dependent manner. Desialylation with a pan-sialic acid hydrolase led to reduction of sialic acid expression on EVs. Moreover, insertion of a high-affinity ligand (LewisY) for DC-SIGN resulted in a four-fold increase of uptake by monocyte-derived DCs. In conclusion, we show that the glycocalyx composition of EVs is a key factor of efficient DC targeting and that modification of the EV glycocalyx potentiates EVs as anti-cancer vaccine.

Targeted proteomics in urinary extracellular vesicles identifies biomarkers for diagnosis and prognosis of prostate cancer.

Rapid and reliable diagnosis of prostate cancer (PCa) is highly desirable as current used methods lack specificity. In addition, identification of PCa biomarkers that can classify patients into high- and low-risk groups for disease progression at early stage will improve treatment decision-making. Here, we describe a set of protein-combination panels in urinary extracellular vesicles (EVs), defined by targeted proteomics and immunoblotting techniques that improve early non-invasive detection and stratification of PCa patients.We report a two-protein combination in urinary EVs that classifies benign and PCa patients (ADSV-TGM4), and a combination of five proteins able to significantly distinguish between high- and low-grade PCa patients (CD63-GLPK5-SPHM-PSA-PAPP). Proteins composing the panels were validated by immunohistochemistry assays in tissue microarrays (TMAs) confirming a strong link between the urinary EVs proteome and alterations in PCa tissues. Moreover, ADSV and TGM4 abundance yielded a high diagnostic potential in tissue and promising TGM4 prognostic power. These results suggest that the proteins identified in urinary EVs distinguishing high- and low grade PCa are a reflection of histological changes that may be a consequence of their functional involvement in PCa development. In conclusion, our study resulted in the identification of protein-combination panels present in urinary EVs that exhibit high sensitivity and specificity for PCa detection and patient stratification. Moreover, our study highlights the potential of targeted proteomic approaches-such as selected reaction monitoring (SRM)-as diagnostic assay for liquid biopsies via urinary EVs to improve diagnosis and prognosis of suspected PCa patients.

Protein Complexes in Urine Interfere with Extracellular Vesicle Biomarker Studies.

Urine exosomes (extracellular vesicles; EVs) contain (micro)RNA (miRNA) and protein biomarkers that are useful for the non-invasive diagnosis of various urological diseases. However, the urinary Tamm-Horsfall protein (THP) complex, which forms at reduced temperatures, may affect EV isolation and may also lead to contamination by other molecules including microRNAs (miRNAs). Therefore, we compared the levels of three miRNAs within the purified EV fraction and THP- protein-network. Urine was collected from healthy donors and EVs were isolated by ultracentrifugation (UC), two commercial kits or sepharose size-exclusion chromatography (SEC). SEC enables the separation of EVs from protein-complexes in urine. After UC, the isolation of EV-miRNA was compared with two commercial kits. The EV isolation efficiency was evaluated by measuring the EV protein markers, Alix and TSG101, CD63 by Western blotting, or miR-375, miR-204 and miR-21 by RT-qPCR. By using commercial kits, EV isolation resulted in either low yields or dissimilar miRNA levels. Via SEC, the EVs were separated from the protein-complex fraction. Importantly, a different ratio was observed between the three miRNAs in the protein fraction compared to the EV fraction. Thus, protein-complexes within urine may influence EV-biomarker studies. Therefore, the characterization of the isolated EV fraction is important to obtain reproducible results.

Non­invasive prostate cancer detection by measuring miRNA variants (isomiRs) in urine extracellular vesicles.

In many cancer types, the expression and function of ~22 nucleotide-long microRNAs (miRNA) is deregulated. Mature miRNAs can be stably detected in extracellular vesicles (EVs) in biofluids, therefore they are considered to have great potential as biomarkers. In the present study, we investigated whether miRNAs have a distinct expression pattern in urine-EVs of prostate cancer (PCa) patients compared to control males. By next generation sequencing, we determined the miRNA expression in a discovery cohort of 4 control men and 9 PCa patients. miRNAs were validated by using a stemloop RT-PCR in an independent cohort of 74 patients (26 control and 48 PCa-patients). Whereas standard mapping protocols identified > 10 PCa associated miRNAs in urinary EVs, miR-21, miR-375 and miR-204 failed to robustly discriminate for disease in a validation study with RT-PCR-detection of mature miRNA sequences. In contrast, we observed that miRNA isoforms (isomiRs) with 3' end modifications were highly discriminatory between samples from control men and PCa patients. Highly differentially expressed isomiRs of miR-21, miR-204 and miR-375 were subsequently validated in an independent group of 74 patients. Receiver-operating characteristic analysis was performed to evaluate the diagnostic performance of three isomiRs, resulting in a 72.9% sensitivity with a high (88%) specificity and an area under the curve (AUC) of 0.866. In comparison, prostate specific antigen had an AUC of 0.707 and measuring the mature form of these miRNAs yielded a lower 70.8% sensitivity and 72% specificity (AUC 0.766). We propose that isomiRs may carry discriminatory information which is useful to generate stronger biomarkers.

Human Salivary Micro-RNA in Patients with Parotid Salivary Gland Neoplasms.

BACKGROUND: Currently, clinical examination, ultrasound scanning (with or without fine needle aspiration cytology), preoperative CT-scan and MRI are available for the differential diagnosis of parotid gland swelling. A preliminary non-invasive salivary diagnostic tool may be helpful in the clinical decision making process. Altered salivary micro-RNA (miRNA) expression levels have been observed in saliva from patients with various cancers. Therefore, we investigated miRNA expression levels in saliva samples from patients with a parotid gland neoplasm using Human miRNA cards in comparison to controls. RESULTS: In the discovery phase, eight miRNAs were identified having different expression levels in patients compared to controls. In the validation phase, the differences in miRNA expression levels between patients and controls were confirmed for seven out of eight discovered miRNAs (p < 0.001). A combination of two miRNAs yielded a receiver-operator-characteristics curve with an AUC of 0.94 (95% CI: 0.87-1.00; sensitivity 91%; specificity 86%). Validation of discovered miRNAs in segregated collected parotid saliva revealed that expression of these miRNAs differ between whole saliva and parotid saliva. CONCLUSIONS: A two miRNA combination can predict the presence of a parotid gland neoplasm. Furthermore, this study suggested that the identified, patient-specific, salivary miRNAs were not derived from the parotid gland itself.