Aberrant Methylation of 20 miRNA Genes Specifically Involved in Various Steps of Ovarian Carcinoma Spread: From Primary Tumors to Peritoneal Macroscopic Metastases.

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.

PPLine: An Automated Pipeline for SNP, SAP, and Splice Variant Detection in the Context of Proteogenomics.

The fundamental mission of the Chromosome-Centric Human Proteome Project (C-HPP) is the research of human proteome diversity, including rare variants. Liver tissues, HepG2 cells, and plasma were selected as one of the major objects for C-HPP studies. The proteogenomic approach, a recently introduced technique, is a powerful method for predicting and validating proteoforms coming from alternative splicing, mutations, and transcript editing. We developed PPLine, a Python-based proteogenomic pipeline providing automated single-amino-acid polymorphism (SAP), indel, and alternative-spliced-variants discovery based on raw transcriptome and exome sequence data, single-nucleotide polymorphism (SNP) annotation and filtration, and the prediction of proteotypic peptides (available at https://sourceforge.net/projects/ppline). In this work, we performed deep transcriptome sequencing of HepG2 cells and liver tissues using two platforms: Illumina HiSeq and Applied Biosystems SOLiD. Using PPLine, we revealed 7756 SAP and indels for HepG2 cells and liver (including 659 variants nonannotated in dbSNP). We found 17 indels in transcripts associated with the translation of alternate reading frames (ARF) longer than 300 bp. The ARF products of two genes, SLMO1 and TMEM8A, demonstrate signatures of caspase-binding domain and Gcn5-related N-acetyltransferase. Alternative splicing analysis predicted novel proteoforms encoded by 203 (liver) and 475 (HepG2) genes according to both Illumina and SOLiD data. The results of the present work represent a basis for subsequent proteomic studies by the C-HPP consortium.

Downregulation of Ribosomal Protein Genes Is Revealed in a Model of Rat Hippocampal Neuronal Culture Activation with GABA(A)R/GlyRa2 Antagonist Picrotoxin.

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.

Zebrafish pigment cells develop directly from persistent highly multipotent progenitors.

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.

Wide-scale identification of novel/eliminated genes responsible for evolutionary transformations.

BACKGROUND: It is generally accepted that most evolutionary transformations at the phenotype level are associated either with rearrangements of genomic regulatory elements, which control the activity of gene networks, or with changes in the amino acid contents of proteins. Recently, evidence has accumulated that significant evolutionary transformations could also be associated with the loss/emergence of whole genes. The targeted identification of such genes is a challenging problem for both bioinformatics and evo-devo research. RESULTS: To solve this problem we propose the WINEGRET method, named after the first letters of the title. Its main idea is to search for genes that satisfy two requirements: first, the desired genes were lost/emerged at the same evolutionary stage at which the phenotypic trait of interest was lost/emerged, and second, the expression of these genes changes significantly during the development of the trait of interest in the model organism. To verify the first requirement, we do not use existing databases of orthologs, but rely purely on gene homology and local synteny by using some novel quickly computable conditions. Genes satisfying the second requirement are found by deep RNA sequencing. As a proof of principle, we used our method to find genes absent in extant amniotes (reptiles, birds, mammals) but present in anamniotes (fish and amphibians), in which these genes are involved in the regeneration of large body appendages. As a result, 57 genes were identified. For three of them, c-c motif chemokine 4, eotaxin-like, and a previously unknown gene called here sod4, essential roles for tail regeneration were demonstrated. Noteworthy, we established that the latter gene belongs to a novel family of Cu/Zn-superoxide dismutases lost by amniotes, SOD4. CONCLUSIONS: We present a method for targeted identification of genes whose loss/emergence in evolution could be associated with the loss/emergence of a phenotypic trait of interest. In a proof-of-principle study, we identified genes absent in amniotes that participate in body appendage regeneration in anamniotes. Our method provides a wide range of opportunities for studying the relationship between the loss/emergence of phenotypic traits and the loss/emergence of specific genes in evolution.

Dysregulation of lncRNA-miRNA-mRNA Interactome as a Marker of Metastatic Process in Ovarian Cancer.

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.

Chromosomal Translocations in NK-Cell Lymphomas Originate from Inter-Chromosomal Contacts of Active rDNA Clusters Possessing Hot Spots of DSBs.

Endogenous hot spots of DNA double-strand breaks (DSBs) are tightly linked with transcription patterns and cancer. There are nine hot spots of DSBs (denoted Pleiades) in human rDNA units that are located exclusively inside the intergenic spacer (IGS). Profiles of Pleiades coincide with the profiles of γ-H2AX, suggesting a high level of in vivo breakage inside rDNA genes. The data were confirmed by microscopic observation of the largest γ-H2AX foci inside nucleoli in interphase chromosomes. Circular chromosome conformation capture (4C) data indicate that the rDNA units often make contact with a specific set of chromosomal regions containing genes that are involved in differentiation and cancer. Interestingly, these regions also often possess hot spots of DSBs that provide the potential for Robertsonian and oncogenic translocations. In this study, we searched for translocations in which rDNA clusters are involved. The whole genome sequence (WGS) data of normal T cells and NK-cell lymphomas from the same individuals revealed numerous translocations in which Pleiades were involved. The sites of these translocations in normal T cells and in the lymphomas were mostly different, although there were also some common sites. The genes at translocations in normal cells and in lymphomas are associated with predominantly non-overlapping lists of genes that are depleted with silenced genes. Our data indicate that rDNA-mediated translocations occur at about the same frequency in the normal T cells and NK-lymphoma cells but differ at particular sites that correspond to open chromatin. We conclude that oncogenic translocations lead to dysregulation of a specific set of genes controlling development. In normal T cells and in NK cells, there are hot spots of translocations at sites possessing strong H3K27ac marks. The data indicate that Pleiades are involved in rDNA-mediated translocation.

A Method for Identification of the Methylation Level of CpG Islands From NGS Data.

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.

Adaptive Changes in the Vestibular System of Land Snail to a 30-Day Spaceflight and Readaptation on Return to Earth.

The vestibular system receives a permanent influence from gravity and reflexively controls equilibrium. If we assume gravity has remained constant during the species' evolution, will its sensory system adapt to abrupt loss of that force? We address this question in the land snail Helix lucorum exposed to 30 days of near weightlessness aboard the Bion-M1 satellite, and studied geotactic behavior of postflight snails, differential gene expressions in statocyst transcriptome, and electrophysiological responses of mechanoreceptors to applied tilts. Each approach revealed plastic changes in the snail's vestibular system assumed in response to spaceflight. Absence of light during the mission also affected statocyst physiology, as revealed by comparison to dark-conditioned control groups. Readaptation to normal tilt responses occurred at ~20 h following return to Earth. Despite the permanence of gravity, the snail responded in a compensatory manner to its loss and readapted once gravity was restored.

Excess fertilizer responsive miRNAs revealed in Linum usitatissimum L.

Effective fertilizer application is necessary to increase crop yields and reduce risk of plant overdosing. It is known that expression level of microRNAs (miRNAs) alters in plants under different nutrient concentrations in soil. The aim of our study was to identify and characterize miRNAs with expression alterations under excessive fertilizer in agriculturally important crop - flax (Linum usitatissimum L.). We have sequenced small RNAs in flax grown under normal and excessive fertilizer using Illumina GAIIx. Over 14 million raw reads was obtained for two small RNA libraries. 84 conserved miRNAs from 20 families were identified. Differential expression was revealed for several flax miRNAs under excessive fertilizer according to high-throughput sequencing data. For 6 miRNA families (miR395, miR169, miR408, miR399, miR398 and miR168) expression level alterations were evaluated on the extended sampling using qPCR. Statistically significant up-regulation was revealed for miR395 under excessive fertilizer. It is known that target genes of miR395 are involved in sulfate uptake and assimilation. However, according to our data alterations of the expression level of miR395 could be associated not only with excess sulfur application, but also with redundancy of other macro- and micronutrients. Furthermore expression level was evaluated for miRNAs and their predicted targets. The negative correlation between miR399 expression and expression of its predicted target ubiquitin-conjugating enzyme E2 gene was shown in flax for the first time. So we suggested miR399 involvement in phosphate regulation in L. usitatissimum. Revealed in our study expression alterations contribute to miRNA role in flax response to excessive fertilizer.

Mining gene expression data for pollutants (dioxin, toluene, formaldehyde) and low dose of gamma-irradiation.

General and specific effects of molecular genetic responses to adverse environmental factors are not well understood. This study examines genome-wide gene expression profiles of Drosophila melanogaster in response to ionizing radiation, formaldehyde, toluene, and 2,3,7,8-tetrachlorodibenzo-p-dioxin. We performed RNA-seq analysis on 25,415 transcripts to measure the change in gene expression in males and females separately. An analysis of the genes unique to each treatment yielded a list of genes as a gene expression signature. In the case of radiation exposure, both sexes exhibited a reproducible increase in their expression of the transcription factors sugarbabe and tramtrack. The influence of dioxin up-regulated metabolic genes, such as anachronism, CG16727, and several genes with unknown function. Toluene activated a gene involved in the response to the toxins, Cyp12d1-p; the transcription factor Fer3's gene; the metabolic genes CG2065, CG30427, and CG34447; and the genes Spn28Da and Spn3, which are responsible for reproduction and immunity. All significantly differentially expressed genes, including those shared among the stressors, can be divided into gene groups using Gene Ontology Biological Process identifiers. These gene groups are related to defense response, biological regulation, the cell cycle, metabolic process, and circadian rhythms. KEGG molecular pathway analysis revealed alteration of the Notch signaling pathway, TGF-beta signaling pathway, proteasome, basal transcription factors, nucleotide excision repair, Jak-STAT signaling pathway, circadian rhythm, Hippo signaling pathway, mTOR signaling pathway, ribosome, mismatch repair, RNA polymerase, mRNA surveillance pathway, Hedgehog signaling pathway, and DNA replication genes. Females and, to a lesser extent, males actively metabolize xenobiotics by the action of cytochrome P450 when under the influence of dioxin and toluene. Finally, in this work we obtained gene expression signatures pollutants (dioxin, toluene), low dose of gamma-irradiation and common molecular pathways for different kind of stressors.

Retrotransposon-based molecular markers for analysis of genetic diversity within the Genus Linum.

SSAP method was used to study the genetic diversity of 22 Linum species from sections Linum, Adenolinum, Dasylinum, Stellerolinum, and 46 flax cultivars. All the studied flax varieties were distinguished using SSAP for retrotransposons FL9 and FL11. Thus, the validity of SSAP method was demonstrated for flax marking, identification of accessions in genebank collections, and control during propagation of flax varieties. Polymorphism of Fl1a, Fl1b, and Cassandra insertions were very low in flax varieties, but these retrotransposons were successfully used for the investigation of Linum species. Species clusterization based on SSAP markers was in concordance with their taxonomic division into sections Dasylinum, Stellerolinum, Adenolinum, and Linum. All species of sect. Adenolinum clustered apart from species of sect. Linum. The data confirmed the accuracy of the separation in these sections. Members of section Linum are not as closely related as members of other sections, so taxonomic revision of this section is desirable. L. usitatissimum accessions genetically distant from modern flax cultivars were revealed in our work. These accessions are of utmost interest for flax breeding and introduction of new useful traits into flax cultivars. The chromosome localization of Cassandra retrotransposon in Linum species was determined.

Empirical potentials for ion binding in proteins.

Empirical potentials for interaction of proteins with intracellular ions are presented. We derive the potentials using a training dataset of the protein 3D structure bank, PDB, based on the statistical analysis of contacts between ions and protein atoms of different types. The potentials are derived using Monte Carlo Reference State, simulating non-interacting structure elements as random 3D points in the structure space. The resulting potentials are detailed, continuous, and cover a wide range of contact distances. The obtained potentials were tested for prediction of ion-binding sites in proteins and are shown to reproduce locations and specificities of ion-binding sites with a high accuracy. A web server is created for predictions of ion-binding sites in proteins.