EV-miRome-wide profiling uncovers miR-320c for detecting metastatic colorectal cancer and monitoring the therapeutic response.

PURPOSE: Molecular composition of circulating small extracellular vesicles (EVs) does not merely reflect the cells of origin, but also is enriched in specific biomolecules directly associated with the cellular transformation. However, while most of the currently identified EV-miRs are only geared towards one-dimensional disease detection, their application for long-term tracking and treatment response monitoring has been largely elusive. METHODS: We established and optimized a rapid, sensitive and robust liquid biopsy sampling method, and further used small RNA sequencing to comprehensively catalogue EV-miRomes in association with the progression and outcome of metastatic colorectal cancer (mCRC). RESULTS: By cross-comparison of EV-miRomes (n = 290) from multi-stage and longitudinal cohorts, we uncovered a 15-EV-miR signature with dual detection and long-term monitoring of tumor size progression for mCRC. From this panel, EV-miR-320c was uncovered as a strong clinical marker - aside from its diagnostic power and a therapeutic monitoring performance superior to carcinoembryonic antigen (CEA), its high expression has also been linked to lower overall survival and a greater likelihood of disease recurrence. Further, integrative analyses of tissue transcriptomic and liquid biopsy implicated this 15-EV-miR signature in programming the mesenchymal-epithelial transition (MET) for distant localization of the metastasized cells and also in creating a tumor-favoring metastatic niche. CONCLUSION: Our clinically-oriented delineation of the mCRC-associated circulating EV-miRomes systematically revealed the functional significance of these liquid biopsy markers and further strengthen their translational potential in mCRC therapeutic monitoring.

Comprehensive transcriptome profiling of Taiwanese colorectal cancer implicates an ethnic basis for pathogenesis.

Colorectal cancer (CRC) is one of the most commonly diagnosed cancers worldwide. While both genetic and environmental factors have been linked to the incidence and mortality associated with CRC, an ethnic aspect of its etiology has also emerged. Since previous large-scale cancer genomics studies are mostly based on samples of European ancestry, the patterns of clinical events and associated mechanisms in other minority ethnic patients suffering from CRC are largely unexplored. We collected 104 paired and adjacent normal tissue and CRC tumor samples from Taiwanese patients and employed an integrated approach - paired expression profiles of mRNAs and microRNAs (miRNAs) combined with transcriptome-wide network analyses - to catalog the molecular signatures of this regional cohort. On the basis of this dataset, which is the largest ever reported for this type of systems analysis, we made the following key discoveries: (1) In comparison to the The Cancer Genome Atlas (TCGA) data, the Taiwanese CRC tumors show similar perturbations in expressed genes but a distinct enrichment in metastasis-associated pathways. (2) Recurrent as well as novel CRC-associated gene fusions were identified based on the sequencing data. (3) Cancer subtype classification using existing tools reveals a comparable distribution of tumor subtypes between Taiwanese cohort and TCGA datasets; however, this similarity in molecular attributes did not translate into the predicted subtype-related clinical outcomes (i.e., death event). (4) To further elucidate the molecular basis of CRC prognosis, we developed a new stratification strategy based on miRNA-mRNA-associated subtyping (MMAS) and consequently showed that repressed WNT signaling activity is associated with poor prognosis in Taiwanese CRC. In summary, our findings of distinct, hitherto unreported biosignatures underscore the heterogeneity of CRC tumorigenesis, support our hypothesis of an ethnic basis of disease, and provide prospects for translational medicine.

A novel antisense RNA ASPACT confers multi-level suppression of PACT and associated signalling.

The innate immune system is the frontline host protection against pathogens. Effective antiviral immunity is elicited upon recognition of viral RNAs by the host pattern recognition receptors. One of the major viral RNA sensors is retinoic acid inducible gene-1, which triggers the production of interferons (IFNs). In turn, this protective response requires another viral sensor and immunity factor interferon-inducible protein kinase RNA activator (PACT/PRKRA). Here, we report the identification and characterization of a novel antisense PACT gene that expresses a non-coding RNA in a convergent and interferon-inducible manner. Publicly available gene structure and expression data revealed that this gene, that we termed ASPACT, overlaps with the 3' -end of the PACT locus and is highly expressed during viral infection. Our results confirm the IFN-ß-inducibility of ASPACT, which is dependent on STAT-1/2. We further discovered that downregulation of ASPACT impacts both the expression and localization of the PACT transcript. At the transcription level, ChIP and ChIRP assays demonstrated that the ASPACT non-coding RNA occupies distinct chromatin regions of PACT gene and is important for promoter recruitment of the epigenetic silencer HDAC1. In parallel, ASPACT was also found to mediate nuclear retention of the PACT mRNA via direct RNA-RNA interaction, as revealed by RNA antisense purification assay. In summary, our results support the model that the non-coding RNA ASPACT acts as a negative regulator of PACT at multiple levels, and reveal a novel regulator of the viral counteractive response.

ADAR1 promotes robust hypoxia signaling via distinct regulation of multiple HIF-1α-inhibiting factors.

Adenosine deaminase acting on RNA (ADAR)-catalyzed adenosine-to-inosine RNA editing is potentially dysregulated in neoplastic progression. However, how this transcriptome recoding process is functionally correlated with tumorigenesis remains largely elusive. Our analyses of RNA editome datasets identify hypoxia-related genes as A-to-I editing targets. In particular, two negative regulators of HIF-1A-the natural antisense transcript HIF1A-AS2 and the ubiquitin ligase scaffold LIMD1-are directly but differentially modulated by ADAR1. We show that HIF1A-AS2 antagonizes the expression of HIF-1A in the immediate-early phase of hypoxic challenge, likely through a convergent transcription competition in cis ADAR1 in turn suppresses transcriptional progression of the antisense gene. In contrast, ADAR1 affects LIMD1 expression post-transcriptionally, by interfering with the cytoplasmic translocation of LIMD1 mRNA and thus protein translation. This multi-tier regulation coordinated by ADAR1 promotes robust and timely accumulation of HIF-1α upon oxygen depletion and reinforces target gene induction and downstream angiogenesis. Our results pinpoint ADAR1-HIF-1α axis as a hitherto unrecognized key regulator in hypoxia.

Tumor-associated intronic editing of HNRPLL generates a novel splicing variant linked to cell proliferation.

Processing of the eukaryotic transcriptome is a dynamic regulatory mechanism that confers genetic diversity, and splicing and adenosine to inosine (A-to-I) RNA editing are well-characterized examples of such processing. Growing evidence reveals the cross-talk between the splicing and RNA editing, but there is a paucity of substantial evidence for its mechanistic details and contribution in a physiological context. Here, our findings demonstrate that tumor-associated differential RNA editing, in conjunction with splicing machinery, regulates the expression of variants of HNRPLL, a gene encoding splicing factor. We discovered an HNRPLL transcript variant containing an additional exon 12A (E12A), which is a substrate of ADAR1 and ADAR2. Adenosine deaminases acting on RNA (ADAR) direct deaminase-dependent expression of the E12A transcript, and ADAR-mediated regulation of E12A is largely splicing-based, and does not affect the stability or nucleocytoplasmic distribution of the transcript. Furthermore, ADAR-mediated modification of exon 12A generates an enhancer for the oncogenic splicing factor SRSF1 and consequently promotes the frequency of alternative splicing. Gene expression profiling by RNA-seq revealed that E12A acts distinctly from HNRPLL and regulates a set of growth-related genes, such as cyclin CCND1 and growth factor receptor TGFBR1 Accordingly, silencing E12A expression leads to impaired clonogenic ability and enhanced sensitivity to doxorubicin, thus highlighting the significance of this alternative isoform in tumor cell survival. In summary, we present the interplay of RNA editing and splicing as a regulatory mechanism of gene expression and also its physiological relevance. These findings extend our understanding of transcriptional dynamics and provide a mechanistic explanation to the link of RNA editors to tumorigenesis.

The PPARγ-SETD8 axis constitutes an epigenetic, p53-independent checkpoint on p21-mediated cellular senescence.

Cellular senescence is a permanent proliferative arrest triggered by genome instability or aberrant growth stresses, acting as a protective or even tumor-suppressive mechanism. While several key aspects of gene regulation have been known to program this cessation of cell growth, the involvement of the epigenetic regulation has just emerged but remains largely unresolved. Using a systems approach that is based on targeted gene profiling, we uncovered known and novel chromatin modifiers with putative link to the senescent state of the cells. Among these, we identified SETD8 as a new target as well as a key regulator of the cellular senescence signaling. Knockdown of SETD8 triggered senescence induction in proliferative culture, irrespectively of the p53 status of the cells; ectopic expression of this epigenetic writer alleviated the extent doxorubicin-induced cellular senescence. This repressive effect of SETD8 in senescence was mediated by directly maintaining the silencing mark H4K20me1 at the locus of the senescence switch gene p21. Further in support of this regulatory link, depletion of p21 reversed this SETD8-mediated cellular senescence. Additionally, we found that PPARγ acts upstream and regulates SETD8 expression in proliferating cells. Downregulation of PPARγ coincided with the senescence induction, while its activation inhibited the progression of this process. Viewed together, our findings delineated a new epigenetic pathway through which the PPARγ-SETD8 axis directly silences p21 expression and consequently impinges on its senescence-inducing function. This implies that SETD8 may be part of a cell proliferation checkpoint mechanism and has important implications in antitumor therapeutics.

A 17-billion-solar-mass black hole in a group galaxy with a diffuse core.

Quasars are associated with and powered by the accretion of material onto massive black holes; the detection of highly luminous quasars with redshifts greater than z = 6 suggests that black holes of up to ten billion solar masses already existed 13 billion years ago. Two possible present-day 'dormant' descendants of this population of 'active' black holes have been found in the galaxies NGC 3842 and NGC 4889 at the centres of the Leo and Coma galaxy clusters, which together form the central region of the Great Wall--the largest local structure of galaxies. The most luminous quasars, however, are not confined to such high-density regions of the early Universe; yet dormant black holes of this high mass have not yet been found outside of modern-day rich clusters. Here we report observations of the stellar velocity distribution in the galaxy NGC 1600--a relatively isolated elliptical galaxy near the centre of a galaxy group at a distance of 64 megaparsecs from Earth. We use orbit superposition models to determine that the black hole at the centre of NGC 1600 has a mass of 17 billion solar masses. The spatial distribution of stars near the centre of NGC 1600 is rather diffuse. We find that the region of depleted stellar density in the cores of massive elliptical galaxies extends over the same radius as the gravitational sphere of influence of the central black holes, and interpret this as the dynamical imprint of the black holes.

Functional Impact of RNA editing and ADARs on regulation of gene expression: perspectives from deep sequencing studies.

Cells regulate gene expression at multiple levels leading to a balance between robustness and complexity within their proteome. One core molecular step contributing to this important balance during metazoan gene expression is RNA editing, such as the co-transcriptional recoding of RNA transcripts catalyzed by the adenosine deaminse acting on RNA (ADAR) family of enzymes. Understanding of the adenosine-to-inosine RNA editing process has been broadened considerably by the next generation sequencing (NGS) technology, which allows for in-depth demarcation of an RNA editome at nucleotide resolution. However, critical issues remain unresolved with regard to how RNA editing cooperates with other transcript-associated events to underpin regulated gene expression. Here we review the growing body of evidence, provided by recent NGS-based studies, that links RNA editing to other mechanisms of post-transcriptional RNA processing and gene expression regulation including alternative splicing, transcript stability and localization, and the biogenesis and function of microRNAs (miRNAs). We also discuss the possibility that systematic integration of NGS data may be employed to establish the rules of an "RNA editing code", which may give us new insights into the functional consequences of RNA editing.

Two ten-billion-solar-mass black holes at the centres of giant elliptical galaxies.

Observational work conducted over the past few decades indicates that all massive galaxies have supermassive black holes at their centres. Although the luminosities and brightness fluctuations of quasars in the early Universe suggest that some were powered by black holes with masses greater than 10 billion solar masses, the remnants of these objects have not been found in the nearby Universe. The giant elliptical galaxy Messier 87 hosts the hitherto most massive known black hole, which has a mass of 6.3 billion solar masses. Here we report that NGC 3842, the brightest galaxy in a cluster at a distance from Earth of 98 megaparsecs, has a central black hole with a mass of 9.7 billion solar masses, and that a black hole of comparable or greater mass is present in NGC 4889, the brightest galaxy in the Coma cluster (at a distance of 103 megaparsecs). These two black holes are significantly more massive than predicted by linearly extrapolating the widely used correlations between black-hole mass and the stellar velocity dispersion or bulge luminosity of the host galaxy. Although these correlations remain useful for predicting black-hole masses in less massive elliptical galaxies, our measurements suggest that different evolutionary processes influence the growth of the largest galaxies and their black holes.

Are halos of collisionless cold dark matter collisionless?

We study whether gravitational scattering of halo dark matter particles by subhalos can connect two seemingly independent problems: the abundance of subhalos in dark matter halos and the cuspiness of the halos' inner density profiles. Our numerical experiments indicate that subhalos can cause the collisionless dark matter particles in the centers of main halos to diffuse. Combined with tidal mass loss of the subhalos, this process introduces significant scatter in the inner density profiles and offers an explanation for the range of profiles seen in both observations and cosmological simulations.

Nonlinear kinetic Sunyaev-Zeldovich effect.

We derive fully nonlinear expressions for temperature fluctuations from the kinetic Sunyaev-Zeldovich (SZ) effect, the scattering of cosmic microwave background photons off hot electrons in bulk motion. Our result reproduces the Ostriker-Vishniac effect to second order in perturbation theory but contains nonlinear corrections to the electron velocities and densities that were neglected previously. We use the recently developed halo model for nonlinear gravitational clustering to compute the nonlinear kinetic SZ power spectrum, which dominates the primary anisotropy on small angular scales.