Neural crest cells are highly multipotent stem cells, but it remains unclear how their fate restriction to specific fates occurs. The direct fate restriction model hypothesises that migrating cells maintain full multipotency, whilst progressive fate restriction envisages fully multipotent cells transitioning to partially-restricted intermediates before committing to individual fates. Using zebrafish pigment cell development as a model, we show applying NanoString hybridization single cell transcriptional profiling and RNAscope in situ hybridization that neural crest cells retain broad multipotency throughout migration and even in post-migratory cells in vivo, with no evidence for partially-restricted intermediates. We find that leukocyte tyrosine kinase early expression marks a multipotent stage, with signalling driving iridophore differentiation through repression of fate-specific transcription factors for other fates. We reconcile the direct and progressive fate restriction models by proposing that pigment cell development occurs directly, but dynamically, from a highly multipotent state, consistent with our recently-proposed Cyclical Fate Restriction model.

Our work aimed to differentiate 20 aberrantly methylated miRNA genes that participate at different stages of development and metastasis of ovarian carcinoma (OvCa) using methylation-specific qPCR in a representative set of clinical samples: 102 primary tumors without and with metastases (to lymph nodes, peritoneum, or distant organs) and 30 peritoneal macroscopic metastases (PMM). Thirteen miRNA genes (MIR107, MIR124-2, MIR124-3, MIR125B-1, MIR127, MIR129-2, MIR130B, MIR132, MIR193A, MIR339, MIR34B/C, MIR9-1, and MIR9-3) were hypermethylated already at the early stages of OvCa, while hypermethylation of MIR1258, MIR137, MIR203A, and MIR375 was pronounced in metastatic tumors, and MIR148A showed high methylation levels specifically in PMM. We confirmed the significant relationship between methylation and expression levels for 11 out of 12 miRNAs analyzed by qRT-PCR. Moreover, expression levels of six miRNAs were significantly decreased in metastatic tumors in comparison with nonmetastatic ones, and downregulation of miR-203a-3p was the most significant. We revealed an inverse relationship between expression levels of miR-203a-3p and those of ZEB1 and ZEB2 genes, which are EMT drivers. We also identified three miRNA genes (MIR148A, MIR9-1, and MIR193A) that likely regulate EMT–MET reversion in the colonization of PMM. According to the Kaplan–Meier analysis, hypermethylation of several examined miRNA genes was associated with poorer overall survival of OvCa patients, and high methylation levels of MIR130B and MIR9-1 were related to the greatest relative risk of death.

Ovarian cancer (OC) is one of the most common types of cancer among malignancies of the female reproductive system. This pathology is asymptomatic until advanced stages and has a poor prognosis. Our study aimed to search for lncRNA–miRNA–mRNA competing triplets that promote ovarian tumorigenesis. For this purpose, we analyzed tumor samples from the TCGA database and verified the results experimentally in a set of 46 paired samples of tumor and matched histologically unchanged ovarian tissues from OC patients. The list of RNAs selected in silico for experimental studies included 13 mRNAs, 10 lncRNAs, and 5 miRNAs related to epithelial–mesenchymal transition and angiogenesis. We evaluated the expression of these RNAs by qRT-PCR and assessed the correlation between levels of miRNAs, mRNAs, and lncRNAs. Sixteen significant triplets were revealed, in some of which, e.g., OIP5-AS1–miR-203a–c-MET and OIP5-AS1–miR-203a–ZEB2, both lncRNA and mRNA had sites for miR-203a direct binding. Transfection of the OVCAR-3 and SKOV-3 cell lines with the miR-203a mimic was used to confirm the novel links of miR-203a with ZEB2 and c-MET in OC. These connections suggest that the interactomes have the potential for diagnostics of metastasis at early onset.

Long-read transcriptome sequencing provides us with a convenient tool for the thorough study of biological processes such as neuronal plasticity. Here, we aimed to perform transcriptional profiling of rat hippocampal primary neuron cultures after stimulation with picrotoxin (PTX) to further understand molecular mechanisms of neuronal activation. To overcome the limitations of short-read RNA-Seq approaches, we performed an Oxford Nanopore Technologies MinION-based long-read sequencing and transcriptome assembly of rat primary hippocampal culture mRNA at three time points after the PTX activation. We used a specific approach to exclude uncapped mRNAs during sample preparation. Overall, we found 23,652 novel transcripts in comparison to reference annotations, out of which ~6000 were entirely novel and mostly transposon-derived loci. Analysis of differentially expressed genes (DEG) showed that 3046 genes were differentially expressed, of which 2037 were upregulated and 1009 were downregulated at 30 min after the PTX application, with only 446 and 13 genes differentially expressed at 1 h and 5 h time points, respectively. Most notably, multiple genes encoding ribosomal proteins, with a high basal expression level, were downregulated after 30 min incubation with PTX; we suggest that this indicates redistribution of transcriptional resources towards activity-induced genes. Novel loci and isoforms observed in this study may help us further understand the functional mRNA repertoire in neuronal plasticity processes. Together with other NGS techniques, differential gene expression analysis of sequencing data obtained using MinION platform might provide a simple method to optimize further study of neuronal plasticity.

In the course of sample preparation for Next Generation Sequencing (NGS), DNA is fragmented by various methods. Fragmentation shows a persistent bias with regard to the cleavage rates of various dinucleotides. With the exception of CpG dinucleotides the previously described biases were consistent with results of the DNA cleavage in solution. Here we computed cleavage rates of all dinucleotides including the methylated CpG and unmethylated CpG dinucleotides using data of the Whole Genome Sequencing datasets of the 1000 Genomes project. We found that the cleavage rate of CpG is significantly higher for the methylated CpG dinucleotides. Using this information, we developed a classifier for distinguishing cancer and healthy tissues based on their CpG islands statuses of the fragmentation. A simple Support Vector Machine classifier based on this algorithm shows an accuracy of 84%. The proposed method allows the detection of epigenetic markers purely based on mechanochemical DNA fragmentation, which can be detected by a simple analysis of the NGS sequencing data.

Neural crest cells (NCCs) are highly multipotent stem cells. A long-standing controversy exists over the mechanism of NCC fate specification, specifically regarding the presence and potency of intermediate progenitors. The direct fate restriction (DFR) model, based on early in vivo clonal studies, hypothesised that intermediates are absent and that migrating cells maintain full multipotency ^1–6 . However, most authors favour progressive fate restriction (PFR) models, with fully multipotent early NCCs (ENCCs) transitioning to partially-restricted intermediates before committing to individual fates ^7–12 . Here, single cell transcriptional profiling of zebrafish pigment cell development leads to us proposing a Cyclical Fate Restriction mechanism of NCC development that reconciles the DFR and PFR models. Our clustering of single NCC Nanostring transcriptional profiles identifies only broadly multipotent intermediate states between ENCCs and differentiated melanocytes and iridophores. Leukocyte tyrosine kinase (Ltk) marks the multipotent progenitor and iridophores, consistent with biphasic ltk expression ^13–15 . Ltk inhibitor and constitutive activation studies support expression at an early multipotent stage, whilst lineage-tracing of ltk- expressing cells reveals their multipotency extends beyond pigment cell-types to neural fates. We conclude that pigment cell development does not involve a conventional PFR mechanism, but instead occurs directly and more dynamically from a broadly multipotent intermediate state.