Nonvesicular extracellular RNAs (nv-exRNAs) constitute the majority of the extracellular RNAome, but little is known about their stability, function, and potential use as disease biomarkers. Herein, we measured the stability of several naked RNAs when incubated in human serum, urine, and cerebrospinal fluid (CSF). We identified extracellularly produced tRNA-derived small RNAs (tDRs) with half-lives of several hours in CSF. Contrary to widespread assumptions, these intrinsically stable small RNAs are full-length tRNAs containing broken phosphodiester bonds (i.e., nicked tRNAs). Standard molecular biology protocols, including phenol-based RNA extraction and heat, induce the artifactual denaturation of nicked tRNAs and the consequent in vitro production of tDRs. Broken bonds are roadblocks for reverse transcriptases, preventing amplification and/or sequencing of nicked tRNAs in their native state. To solve this, we performed enzymatic repair of nicked tRNAs purified under native conditions, harnessing the intrinsic activity of phage and bacterial tRNA repair systems. Enzymatic repair regenerated an RNase R-resistant tRNA-sized band in northern blot and enabled RT-PCR amplification of full-length tRNAs. We also separated nicked tRNAs from tDRs by chromatographic methods under native conditions, identifying nicked tRNAs inside stressed cells and in vesicle-depleted human biofluids. Dissociation of nicked tRNAs produces single-stranded tDRs that can be spontaneously taken up by human epithelial cells, positioning stable nv-exRNAs as potentially relevant players in intercellular communication pathways.

The extracellular RNA communication consortium (ERCC) is an NIH-funded program aiming to promote the development of new technologies, resources, and knowledge about exRNAs and their carriers. After Phase 1 (2013-2018), Phase 2 of the program (ERCC2, 2019-2023) aims to fill critical gaps in knowledge and technology to enable rigorous and reproducible methods for separation and characterization of both bulk populations of exRNA carriers and single EVs. ERCC2 investigators are also developing new bioinformatic pipelines to promote data integration through the exRNA atlas database. ERCC2 has established several Working Groups (Resource Sharing, Reagent Development, Data Analysis and Coordination, Technology Development, nomenclature, and Scientific Outreach) to promote collaboration between ERCC2 members and the broader scientific community. We expect that ERCC2's current and future achievements will significantly improve our understanding of exRNA biology and the development of accurate and efficient exRNA-based diagnostic, prognostic, and theranostic biomarker assays.

Purposes

Most molecular-based published studies on breast cancer do not adequately represent the unique and diverse genetic admixture of the Latin American population. Searching for similarities and differences in molecular pathways associated with these tumors and evaluating its impact on prognosis may help to select better therapeutic approaches.

Patients and methods

We collected clinical, pathological, and transcriptomic data of a multi-country Latin American cohort of 1,071 stage II-III breast cancer patients of the Molecular Profile of Breast Cancer Study (MPBCS) cohort. The 5-year prognostic ability of intrinsic (transcriptomic-based) PAM50 and immunohistochemical classifications, both at the cancer-specific (OSC) and disease-free survival (DFS) stages, was compared. Pathway analyses (GSEA, GSVA and MetaCore) were performed to explore differences among intrinsic subtypes.

Results

PAM50 classification of the MPBCS cohort defined 42·6% of tumors as LumA, 21·3% as LumB, 13·3% as HER2E and 16·6% as Basal. Both OSC and DFS for LumA tumors were significantly better than for other subtypes, while Basal tumors had the worst prognosis. While the prognostic power of traditional subtypes calculated with hormone receptors (HR), HER2 and Ki67 determinations showed an acceptable performance, PAM50-derived risk of recurrence best discriminated low, intermediate and high-risk groups. Transcriptomic pathway analysis showed high proliferation (i.e. cell cycle control and DNA damage repair) associated with LumB, HER2E and Basal tumors, and a strong dependency on the estrogen pathway for LumA. Terms related to both innate and adaptive immune responses were seen predominantly upregulated in Basal tumors, and, to a lesser extent, in HER2E, with respect to LumA and B tumors.

Conclusions

This is the first study that assesses molecular features at the transcriptomic level in a multicountry Latin American breast cancer patient cohort. Hormone-related and proliferation pathways that predominate in PAM50 and other breast cancer molecular classifications are also the main tumor-driving mechanisms in this cohort and have prognostic power. The immune-related features seen in the most aggressive subtypes may pave the way for therapeutic approaches not yet disseminated in Latin America.

Clinical trial registration

ClinicalTrials.gov (Identifier: NCT02326857).