Serum small RNAs in metastatic colorectal cancer predict response to chemotherapy and characterize high-risk patients.

Metastatic colorectal cancer (mCRC) presents significant challenges in clinical management due to its heterogeneity and variable response to treatment. In this study, we conducted comprehensive small RNA (sRNA) sequencing analyses to identify sRNA biomarkers associated with survival and treatment response in mCRC patients. We measured serum sRNAs before and after chemotherapy treatment in a discovery cohort of 189 mCRC patients. Our analysis revealed 25 microRNAs (miRNA) as significantly associated with overall survival at baseline. We found that 11 of the 25 significant miRNAs were also significant in an independent validation cohort of 20 mCRC patients, including the top five miRNAs from the discovery cohort. Importantly, all but four of the 25 significant miRNAs from the discovery cohort had hazard ratios in the same direction in the validation cohort. Among the 25 significant miRNAs, we identified the miR-320 family of miRNAs as the strongest independent prognostic marker, with high baseline levels correlating with poor survival outcomes. Furthermore, post-treatment levels of the same miRNAs were even more predictive of overall survival, emphasizing the prognostic value of serum changes in miRNA levels before and after treatment. Moreover, we observed significant changes in serum miRNAs and other sRNAs when comparing samples before and after chemotherapy, with distinct expression patterns between responders and non-responders. Leveraging these differential expression patterns, we established a serum sRNA signature that accurately predicts response to chemotherapy with an area under the curve (AUC) of 0.8. In summary, our study highlights the prognostic and predictive potential of sRNA biomarkers in mCRC, offering valuable insights into patient stratification and personalized treatment approaches.

Comprehensive small RNA-sequencing of primary myeloma cells identifies miR-105-5p as a predictor of patient survival.

BACKGROUND: Small RNAs (sRNAs), a heterogenous group of non-coding RNAs, are emerging as promising molecules for cancer patient risk stratification and as players in tumour pathogenesis. Here, we have studied microRNAs (miRNAs) and other sRNAs in relation to survival and disease severity in multiple myeloma. METHODS: We comprehensively characterised sRNA expression in multiple myeloma patients by performing sRNA-sequencing on myeloma cells isolated from bone marrow aspirates of 86 myeloma patients. The sRNA expression profiles were correlated with the patients' clinical data to investigate associations with survival and disease subgroups, by using cox proportional hazards (coxph) -models and limma-voom, respectively. A publicly available sRNA dataset was used as external validation (n = 151). RESULTS: We show that multiple miRNAs are differentially expressed between ISS Stage I and III. Interestingly, we observed the downregulation of seven different U2 spliceosomal RNAs, a type of small nuclear RNAs in severe disease stages. Further, by a discovery-based approach, we identified miRNA miR-105-5p as a predictor of poor overall survival (OS) in multiple myeloma. Multivariate analysis showed that miR-105-5p predict OS independently of established disease markers. CONCLUSIONS: Overexpression of miR-105-5p in myeloma cells correlates with reduced OS, potentially improving prognostic risk stratification in multiple myeloma.

The first mitochondrial 5-methylcytosine map in a non-model teleost (Oreochromis niloticus) reveals extensive strand-specific and non-CpG methylation.

DNA methylation is one of the main epigenetic mechanisms that regulate gene expression in a manner that depends on the genomic context and varies considerably across taxa. This DNA modification was first found in nuclear genomes of eukaryote several decades ago and it has also been described in mitochondrial DNA. It has recently been shown that mitochondrial DNA is extensively methylated in mammals and other vertebrates. Our current knowledge of mitochondrial DNA methylation in fish is very limited, especially in non-model teleosts. In this study, using whole-genome bisulfite sequencing, we determined methylation patterns within non-CpG (CH) and CpG (CG) contexts in the mitochondrial genome of Nile tilapia, a non-model teleost of high economic importance. Our results demonstrate the presence of mitochondrial DNA methylation in this species predominantly within a non-CpG context, similarly to mammals. We found a strand-specific distribution of methylation, in which highly methylated cytosines were located on the minus strand. The D-loop region had the highest mean methylation level among all mitochondrial loci. Our data provide new insights into the potential role of epigenetic mechanisms in regulating metabolic flexibility of mitochondria in fish, with implications in various biological processes, such as growth and development.

Epigenetic mapping of the somatotropic axis in Nile tilapia reveals differential DNA hydroxymethylation marks associated with growth.

In vertebrates, the somatotropic axis comprising the pituitary gland, liver and muscle plays a major role in myogenesis. Its output in terms of muscle growth is highly affected by nutritional and environmental cues, and thus likely epigenetically regulated. Hydroxymethylation is emerging as a DNA modification that modulates gene expression but a holistic characterization of the hydroxymethylome of the somatotropic axis has not been investigated to date. Using reduced representation 5-hydroxymethylcytosine profiling we demonstrate tissue-specific localization of 5-hydroxymethylcytosines at single nucleotide resolution. Their abundance within gene bodies and promoters of several growth-related genes supports their pertinent role in gene regulation. We propose that cytosine hydroxymethylation may contribute to the phenotypic plasticity of growth through epigenetic regulation of the somatotropic axis.

Uracil-DNA glycosylase UNG1 isoform variant supports class switch recombination and repairs nuclear genomic uracil.

UNG is the major uracil-DNA glycosylase in mammalian cells and is involved in both error-free base excision repair of genomic uracil and mutagenic uracil-processing at the antibody genes. However, the regulation of UNG in these different processes is currently not well understood. The UNG gene encodes two isoforms, UNG1 and UNG2, each possessing unique N-termini that mediate translocation to the mitochondria and the nucleus, respectively. A strict subcellular localization of each isoform has been widely accepted despite a lack of models to study them individually. To determine the roles of each isoform, we generated and characterized several UNG isoform-specific mouse and human cell lines. We identified a distinct UNG1 isoform variant that is targeted to the cell nucleus where it supports antibody class switching and repairs genomic uracil. We propose that the nuclear UNG1 variant, which in contrast to UNG2 lacks a PCNA-binding motif, may be specialized to act on ssDNA through its ability to bind RPA. RPA-coated ssDNA regions include both transcribed antibody genes that are targets for deamination by AID and regions in front of the moving replication forks. Our findings provide new insights into the function of UNG isoforms in adaptive immunity and DNA repair.

Kisspeptin Influences the Reproductive Axis and Circulating Levels of microRNAs in Senegalese Sole.

Kisspeptin regulates puberty and reproduction onset, acting upstream of the brain-pituitary-gonad (HPG) axis. This study aimed to test a kisspeptin-based hormonal therapy on cultured Senegalese sole (G1) breeders, known to have reproductive dysfunctions. A single intramuscular injection of KISS2-10 decapeptide (250 µg/kg) was tested in females and males during the reproductive season, and gonad maturation, sperm motility, plasma levels of gonadotropins (Fsh and Lh) and sex steroids (11-ketotestosterone, testosterone and estradiol), as well as changes in small non-coding RNAs (sncRNAs) in plasma, were investigated. Fsh, Lh, and testosterone levels increased after kisspeptin injection in both sexes, while sperm analysis did not show differences between groups. Let7e, miR-199a-3p and miR-100-5p were differentially expressed in females, while miR-1-3p miRNA was up-regulated in kisspeptin-treated males. In silico prediction of mRNAs targeted by miRNAs revealed that kisspeptin treatment might affect paracellular transporters, regulate structural and functional polarity of cells, neural networks and intracellular trafficking in Senegalese sole females; also, DNA methylation and sphingolipid metabolism might be altered in kisspeptin-treated males. Results demonstrated that kisspeptin stimulated gonadotropin and testosterone secretion in both sexes and induced an unanticipated alteration of plasma miRNAs, opening new research venues to understand how this neuropeptide impacts in fish HPG axis.

Small RNA expression from viruses, bacteria and human miRNAs in colon cancer tissue and its association with microsatellite instability and tumor location.

BACKGROUND: MicroRNAs (miRNA) and other small RNAs are frequently dysregulated in cancer and are promising biomarkers for colon cancer. Here we profile human, virus and bacteria small RNAs in normal and tumor tissue from early stage colon cancer and correlate the expression with clinical parameters. METHODS: Small RNAs from colon cancer tissue and adjacent normal mucosa of 48 patients were sequenced using Illumina high-throughput sequencing. Clinical parameters were correlated with the small RNA expression data using linear models. We performed a meta-analysis by comparing publicly available small RNA sequencing datasets with our original sequencing data to confirm the main findings. RESULTS: We identified 331 differentially expressed miRNAs between tumor and normal samples. We found that the major changes in miRNA expression between left and right colon are due to miRNAs located within the Hox-developmental genes, including miR-10b, miR-196b and miR-615. Further, we identified new miRNAs associated with microsatellite instability (MSI), including miR-335, miR-26 and miR-625. We performed a meta-analysis on all publicly available miRNA-seq datasets and identified 117 common miRNAs that were differentially expressed between tumor and normal tissue. The miRNAs miR-135b and miR-31 were the most significant upregulated miRNA in tumor across all datasets. The miRNA miR-133a was the most strongly downregulated miRNA in our dataset and also showed consistent downregulation in the other datasets. The miRNAs associated with MSI and tumor location in our data showed similar changes in the other datasets. Finally, we show that small RNAs from Epstein-Barr virus and Fusobacterium nucleatum are differentially expressed between tumor and normal adjacent tissue. CONCLUSIONS: Small RNA profiling in colon cancer tissue revealed novel RNAs associated with MSI and tumor location. We show that Fusobacterium nucleatum are detectable at the RNA-level in colon tissue, and that both Fusobacterium nucleatum and Epstein-Barr virus separate tumor and normal tissue.

Comprehensive transcriptomic analyses of tissue, serum, and serum exosomes from hepatocellular carcinoma patients.

BACKGROUND: The expression of microRNAs (miRNAs) is a promising prognostic and diagnostic tool in hepatocellular carcinoma (HCC). Here we performed small RNA sequencing (sRNA-seq) of tissue, serum and serum exosomes to investigate changes in miRNA expression between the different sample types and correlated the expression with clinical parameters. We also performed gene expression arrays on tumor and normal tissue. RESULTS: Paired tissue, serum and serum exosomes sequencing revealed consistent positive correlation of miR-21 between serum exosomes and tumor tissue, indicating that miR-21 could be exported from tissue to circulation via exosomes. We found that let-7 miRNAs are generally upregulated in serum exosomes compared to whole serum, indicating that these miRNAs could be preferentially loaded into exosomes. Comparing serum from HCC patients with serum from healthy individuals revealed a global increase of miRNAs in serum from HCC patients, including an almost 4-fold increase of several miRNAs, including the liver-specific miR-122. When correlating miRNA expression with clinical parameters we detected significant association between hepatitis B virus (HBV) infection and miR-122 in serum as well as several serum and tissue-miRNAs that correlated with surgery type. We found that miR-141 and miR-146 correlated with cirrhosis in tumor tissue and normal tissue, respectively. Finally, high expression of miR-21 in tumors were associated with poor survival. Focusing on gene expression we found several significant messenger RNAs (mRNAs) between tumor and normal tissue and a Gene Ontology (GO) analysis revealed that these changes were mainly related to cell cycle and metabolism. Further, we detected mRNAs that correlated with cirrhosis and HBV infection in tissue. Finally, GO analysis of predicted targets for miRNAs down-regulated in tumor found that these were enriched for functions related to collagen synthesis. CONCLUSIONS: Our combined data point to altered miRNA and mRNA expression contributing to both generally impaired lipid metabolism and increased cell proliferation and a miRNA-driven increase in collagen synthesis in HCC. Our results further indicate a correlation in miRNA expression between exosomes, serum, and tissue samples suggesting export from tumors via exosomes. This correlation could provide a basis for a more tumor-specific miRNA profile in serum.

Transcription profiling during the cell cycle shows that a subset of Polycomb-targeted genes is upregulated during DNA replication.

Genome-wide gene expression analyses of the human somatic cell cycle have indicated that the set of cycling genes differ between primary and cancer cells. By identifying genes that have cell cycle dependent expression in HaCaT human keratinocytes and comparing these with previously identified cell cycle genes, we have identified three distinct groups of cell cycle genes. First, housekeeping genes enriched for known cell cycle functions; second, cell type-specific genes enriched for HaCaT-specific functions; and third, Polycomb-regulated genes. These Polycomb-regulated genes are specifically upregulated during DNA replication, and consistent with being epigenetically silenced in other cell cycle phases, these genes have lower expression than other cell cycle genes. We also find similar patterns in foreskin fibroblasts, indicating that replication-dependent expression of Polycomb-silenced genes is a prevalent but unrecognized regulatory mechanism.

Toll-like receptor signaling in multiple myeloma cells promotes the expression of pro-survival genes B-cell lymphoma 2 and MYC and modulates the expression of B-cell maturation antigen.

Infections are common in plasma cell cancer multiple myeloma (MM) due to disease-related immune deficiencies and cancer treatment. Myeloma cells express Toll-like receptors (TLRs), and TLR activation has been shown to induce proliferative and pro-survival signals in cancer cells. MM is a complex and heterogeneous disease, and expression levels of TLRs as well as downstream signaling components are likely to differ between patients. Here, we show that in a large cohort of patients, TLR1, TLR4, TLR6, TLR9, and TLR10 are the most highly expressed in primary CD138+ cells. Using an MM cell line expressing TLR4 and TLR9 as a model, we demonstrate that TLR4 and TLR9 activation promoted the expression of well-established pro-survival and oncogenes in MM such as MYC, IRF4, NFKB, and BCL2. TLR4 and TLR9 activation inhibited the efficacy of proteasome inhibitors bortezomib and carfilzomib, drugs used in the treatment of MM. Inhibiting the autophagosome-lysosome protein degradation pathway by hydroxychloroquine (HCQ) diminished the protective effect of TLR activation on proteasome inhibitor-induced cytotoxicity. We also found that TLR signaling downregulated the expression of TNFRSF17, the gene encoding for B-cell maturation antigen (BCMA). MYC, BCL2, and BCL2L1 were upregulated in approximately 50% of primary cells, while the response to TLR signaling in terms of TNFRSF17 expression was dichotomous, as an equal fraction of patients showed upregulation and downregulation of the gene. While proteasome inhibitors are part of first-line MM treatment, several of the new anti-MM immune therapeutic drugs target BCMA. Thus, TLR activation may render MM cells less responsive to commonly used anti-myeloma drugs.

IL-32 is induced by activation of toll-like receptors in multiple myeloma cells.

Multiple myeloma (MM) is a hematological cancer characterized by accumulation of malignant plasma cells in the bone marrow. The patients are immune suppressed and suffer from recurrent and chronic infections. Interleukin-32 is a non-conventional, pro-inflammatory cytokine expressed in a subgroup of MM patients with a poor prognosis. IL-32 has also been shown to promote proliferation and survival of the cancer cells. Here we show that activation of toll-like receptors (TLRs) promotes expression of IL-32 in MM cells through NFκB activation. In patient-derived primary MM cells, IL-32 expression is positively associated with expression of TLRs. Furthermore, we found that several TLR genes are upregulated from diagnosis to relapse in individual patients, predominantly TLRs sensing bacterial components. Interestingly, upregulation of these TLRs coincides with an increase in IL-32. Taken together, these results support a role for IL-32 in microbial sensing in MM cells and suggest that infections can induce expression of this pro-tumorigenic cytokine in MM patients.

TiSA: TimeSeriesAnalysis-a pipeline for the analysis of longitudinal transcriptomics data.

Improved transcriptomic sequencing technologies now make it possible to perform longitudinal experiments, thus generating a large amount of data. Currently, there are no dedicated or comprehensive methods for the analysis of these experiments. In this article, we describe our TimeSeries Analysis pipeline (TiSA) which combines differential gene expression, clustering based on recursive thresholding, and a functional enrichment analysis. Differential gene expression is performed for both the temporal and conditional axes. Clustering is performed on the identified differentially expressed genes, with each cluster being evaluated using a functional enrichment analysis. We show that TiSA can be used to analyse longitudinal transcriptomic data from both microarrays and RNA-seq, as well as small, large, and/or datasets with missing data points. The tested datasets ranged in complexity, some originating from cell lines while another was from a longitudinal experiment of severity in COVID-19 patients. We have also included custom figures to aid with the biological interpretation of the data, these plots include Principal Component Analyses, Multi Dimensional Scaling plots, functional enrichment dotplots, trajectory plots, and complex heatmaps showing the broad overview of results. To date, TiSA is the first pipeline to provide an easy solution to the analysis of longitudinal transcriptomics experiments.

Increased levels of microRNA-320 in blood serum and plasma is associated with imminent and advanced lung cancer.

Lung cancer (LC) incidence is increasing globally and altered levels of microRNAs (miRNAs) in blood may contribute to identification of individuals with LC. We identified miRNAs differentially expressed in peripheral blood at LC diagnosis and evaluated, in pre-diagnostic blood specimens, how long before diagnosis expression changes in such candidate miRNAs could be detected. We identified upregulated candidate miRNAs in plasma specimens from a hospital-based study sample of 128 patients with confirmed LC and 62 individuals with suspected but confirmed negative LC (FalsePos). We then evaluated the expression of candidate miRNAs in pre-diagnostic plasma or serum specimens of 360 future LC cases and 375 matched controls. There were 1663 miRNAs detected in diagnostic specimens, nine of which met our criteria for candidate miRNAs. Higher expression of three candidates, miR-320b, 320c, and 320d, was associated with poor survival, independent of LC stage and subtype. Moreover, miR-320c and miR-320d expression was higher in pre-diagnostic specimens collected within 2 years of LC diagnosis. Our results indicated that elevated levels of miR-320c and miR-320d may be early indications of imminent and advanced LC.

Small RNAs are differentially expressed in autoimmune and non-autoimmune diabetes and controls.

Objective: Diabetes is a heterogeneous disease and a precise diagnosis of diabetes subgroups is necessary to initiate proper early treatment and clinical management of the disease. Circulating small RNAs (sRNAs) are potentially diagnostic biomarkers in diseases, including diabetes. Here we aimed to examine whether profiles of circulating sRNAs differed between patients with autoimmune and non-autoimmune diabetes and non-diabetic controls. Design: This cross-sectional case-control study included participants from the third survey of the HUNT study. Methods: We performed sRNA sequencing in serum from adult-onset type 1 diabetes (n = 51), type 2 diabetes (n = 50) and latent autoimmune diabetes in adult (LADA, n = 51), as well as non-diabetic HUNT3 participants as control group (n = 51). Differential expression analysis of the sRNAs was performed in R using limma-voom. Results: We identified differences in sRNA expression between autoimmune (type 1 diabetes and LADA) and non-autoimmune diabetes (type 2 diabetes) and between patients with diabetes and non-diabetic controls. Focusing on miRNA, we identified 10 differentially expressed mature miRNAs and 30 differentially expressed miRNA variants (isomiRs). We also identified significant changes within other sRNA classes, including a pronounced downregulation of a tRNA fragment in patients with diabetes compared to non-diabetic controls. We created cross-validated sRNA signatures based on the significant sRNAs that distinguished patients with diabetes from non-diabetic controls, and autoimmune from non-autoimmune diabetes, with high specificity and sensitivity. sRNA profiles did not distinguish between type 1 diabetes and LADA. Conclusions: Circulating sRNAs are differentially expressed between patients with diabetes and non-diabetic controls and between autoimmune and non-autoimmune diabetes.

Assessment of Two Commercial Comprehensive Gene Panels for Personalized Cancer Treatment.

(1) Background: Analysis of tumor DNA by next-generation sequencing (NGS) plays various roles in the classification and management of cancer. This study aimed to assess the performance of two similar and large, comprehensive gene panels with a focus on clinically relevant variant detection and tumor mutation burden (TMB) assessment; (2) Methods: DNA from 19 diagnostic small cell lung cancer biopsies and an AcroMetrix™ assessment sample with >500 mutations were sequenced using Oncomine™ Comprehensive Assay Plus (OCAP) on the Ion Torrent platform and TruSight Oncology 500 Assay (TSO500) on the Illumina platform; (3) Results: OCAP and TSO500 achieved comparable NGS quality, such as mean read coverage and mean coverage uniformity. A total of 100% of the variants in the diagnostic samples and 80% of the variants in the AcroMetrix™ assessment sample were detected by both panels, and the panels reported highly similar variant allele frequency. A proportion of 14/19 (74%) samples were classified in the same TMB category; (4) Conclusions: Comparable results were obtained using OCAP and TSO500, suggesting that both panels could be applied to screen patients for enrolment in personalized cancer treatment trials.

Paired miRNA- and messenger RNA-sequencing identifies novel miRNA-mRNA interactions in multiple myeloma.

Multiple myeloma (MM) is an incurable cancer of terminally differentiated plasma cells that proliferate in the bone marrow. miRNAs are promising biomarkers for risk stratification in MM and several miRNAs are shown to have a function in disease pathogenesis. However, to date, surprisingly few miRNA-mRNA interactions have been described for and functionally validated in MM. In this study, we performed miRNA-seq and mRNA-seq on CD138 + cells isolated from bone marrow aspirates of 86 MM patients to identify novel interactions between sRNAs and mRNAs. We detected 9.8% significantly correlated miRNA-mRNA pairs of which 5.17% were positively correlated and 4.65% were negatively correlated. We found that miRNA-mRNA pairs that were predicted by in silico target-prediction algorithms were more negatively correlated than non-target pairs, indicating functional miRNA targeting and that correlation between miRNAs and mRNAs from patients can be used to identify miRNA-targets. mRNAs for negatively correlated miRNA-mRNA target pairs were associated with gene ontology terms such as autophagy, protein degradation and endoplasmic stress response, reflecting important processes in MM. Targets for two specific miRNAs, miR-125b-5p and miR-365b-3p, were functionally validated in MM cell line transfection experiments followed by RNA-sequencing and qPCR. In summary, we identified functional miRNA-mRNA target pairs by correlating miRNA and mRNA data from primary MM cells. We identified several target pairs that are of potential interest for further studies. The data presented here may serve as a hypothesis-generating knowledge base for other researchers in the miRNA/MM field. We also provide an interactive web application that can be used to exploit the miRNA-target interactions as well as clinical parameters associated to these target-pairs.

Intracellular IL-32 regulates mitochondrial metabolism, proliferation, and differentiation of malignant plasma cells.

Interleukin-32 (IL-32) is a nonclassical cytokine expressed in cancers, inflammatory diseases, and infections. Its expression is regulated by two different oxygen sensing systems; HIF1α and cysteamine dioxygenase (ADO), indicating that IL-32 may be involved in the response to hypoxia. We here demonstrate that endogenously expressed, intracellular IL-32 interacts with components of the mitochondrial respiratory chain and promotes oxidative phosphorylation. Knocking out IL-32 in three myeloma cell lines reduced cell survival and proliferation in vitro and in vivo. High-throughput transcriptomic and MS-metabolomic profiling of IL-32 KO cells revealed that cells depleted of IL-32 had perturbations in metabolic pathways, with accumulation of lipids, pyruvate precursors, and citrate. IL-32 was expressed in a subgroup of myeloma patients with inferior survival, and primary myeloma cells expressing IL-32 had a gene signature associated with immaturity, proliferation, and oxidative phosphorylation. In conclusion, we demonstrate a previously unrecognized role of IL-32 in the regulation of plasma cell metabolism.

sMETASeq: Combined Profiling of Microbiota and Host Small RNAs.

Understanding microbial communities' roles in human health and disease requires methods that accurately characterize the microbial composition and their activity and effects within human biological samples. We present sMETASeq (small RNA Metagenomics by Sequencing), a novel method that uses sequencing of small RNAs to jointly measure host small RNA expression and create metagenomic profiles and detect small bacterial RNAs. We evaluated the performance of sMETASeq on a mock bacterial community and demonstrated its use on different human samples, including colon cancer, oral leukoplakia, cervix cancer, and a panel of human biofluids. In all datasets, the detected microbes reflected the biology of the different sample types.

Major gene expression changes and epigenetic remodelling in Nile tilapia muscle after just one generation of domestication.

The historically recent domestication of fishes has been essential to meet the protein demands of a growing human population. Selection for traits of interest during domestication is a complex process whose epigenetic basis is poorly understood. Cytosine hydroxymethylation is increasingly recognized as an important DNA modification involved in epigenetic regulation. In the present study, we investigated if hydroxymethylation plays a role in fish domestication and demonstrated for the first time at a genome-wide level and single nucleotide resolution that the muscle hydroxymethylome changes after a single generation of Nile tilapia (Oreochromis niloticus, Linnaeus) domestication. The overall decrease in hydroxymethylcytosine levels was accompanied by the downregulation of 2015 genes in fish reared in captivity compared to their wild progenitors. In contrast, several myogenic and metabolic genes that can affect growth potential were upregulated. There were 126 differentially hydroxymethylated cytosines between groups, which were not due to genetic variation; they were associated with genes involved in immune-, growth- and neuronal-related pathways. Taken together, our data unveil a new role for DNA hydroxymethylation in epigenetic regulation of fish domestication with impact in aquaculture and implications in artificial selection, environmental adaptation and genome evolution.

Identification of metastasis-associated microRNAs in serum from rectal cancer patients.

MicroRNAs (miRNAs) are promising prognostic and diagnostic biomarkers due to their high stability in blood. Here we investigate the expression of miRNAs and other noncoding (nc) RNAs in serum of rectal cancer patients. Serum from 96 rectal cancer patients was profiled using small RNA sequencing and expression of small RNAs was correlated with the clinicopathological characteristics of the patients. Multiple classes of RNAs were detected, including miRNAs and fragments of tRNAs, snoRNAs, long ncRNAs, and other classes of RNAs. Several miRNAs, miRNA variants (isomiRs) and other ncRNAs were differentially expressed between Stage IV and Stage I-III rectal cancer patients, including several members of the miR-320 family. Furthermore, we show that high expression of miR-320d as well as one tRNA fragment is associated with poor survival. We also show that several miRNAs and isomiRs are differentially expressed between patients receiving preoperative chemoradiotherapy and patients who did not receive any treatment before serum collection. In summary, our study shows that the expression of miRNAs and other small ncRNAs in serum may be used to predict distant metastasis and survival in rectal cancer.