Tumor Suppressor Properties of Small C-Terminal Domain Phosphatases in Clear Cell Renal Cell Carcinoma.

Clear cell renal cell carcinoma (ccRCC) accounts for 80-90% of kidney cancers worldwide. Small C-terminal domain phosphatases CTDSP1, CTDSP2, and CTDSPL (also known as SCP1, 2, 3) are involved in the regulation of several important pathways associated with carcinogenesis. In various cancer types, these phosphatases may demonstrate either antitumor or oncogenic activity. Tumor-suppressive activity of these phosphatases in kidney cancer has been shown previously, but in general case, the antitumor activity may be dependent on the choice of cell line. In the present work, transfection of the Caki-1 cell line (ccRCC morphologic phenotype) with expression constructs containing the coding regions of these genes resulted in inhibition of cell growth in vitro in the case of CTDSP1 (p < 0.001) and CTDSPL (p < 0.05) but not CTDSP2. The analysis of The Cancer Genome Atlas (TCGA) data showed differential expression of some of CTDSP genes and of their target, RB1. These results were confirmed by quantitative RT-PCR using an independent sample of primary ccRCC tumors (n = 52). We observed CTDSPL downregulation and found a positive correlation of expression for two gene pairs: CTDSP1 and CTDSP2 (rs = 0.76; p < 0.001) and CTDSPL and RB1 (rs = 0.38; p < 0.05). Survival analysis based on TCGA data demonstrated a strong association of lower expression of CTDSP1, CTDSP2, CTDSPL, and RB1 with poor survival of ccRCC patients (p < 0.001). In addition, according to TCGA, CTDSP1, CTDSP2, and RB1 were differently expressed in two subtypes of ccRCC-ccA and ccB, characterized by different survival rates. These results confirm that CTDSP1 and CTDSPL have tumor suppressor properties in ccRCC and reflect their association with the more aggressive ccRCC phenotype.

Transcriptome Profiling of Prostate Cancer, Considering Risk Groups and the TMPRSS2-ERG Molecular Subtype.

Molecular heterogeneity in prostate cancer (PCa) is one of the key reasons underlying the differing likelihoods of recurrence after surgical treatment in individual patients of the same clinical category. In this study, we performed RNA-Seq profiling of 58 localized PCa and 43 locally advanced PCa tissue samples obtained as a result of radical prostatectomy on a cohort of Russian patients. Based on bioinformatics analysis, we examined features of the transcriptome profiles within the high-risk group, including within the most commonly represented molecular subtype, TMPRSS2-ERG. The most significantly affected biological processes in the samples were also identified, so that they may be further studied in the search for new potential therapeutic targets for the categories of PCa under consideration. The highest predictive potential was found with the EEF1A1P5, RPLP0P6, ZNF483, CIBAR1, HECTD2, OGN, and CLIC4 genes. We also reviewed the main transcriptome changes in the groups at intermediate risk of PCa-Gleason Score 7 (groups 2 and 3 according to the ISUP classification)-on the basis of which the LPL, MYC, and TWIST1 genes were identified as promising additional prognostic markers, the statistical significance of which was confirmed using qPCR validation.

Lymphatic Dissemination in Prostate Cancer: Features of the Transcriptomic Profile and Prognostic Models.

Radical prostatectomy is the gold standard treatment for prostate cancer (PCa); however, it does not always completely cure PCa, and patients often experience a recurrence of the disease. In addition, the clinical and pathological parameters used to assess the prognosis and choose further tactics for treating a patient are insufficiently informative and need to be supplemented with new markers. In this study, we performed RNA-Seq of PCa tissue samples, aimed at identifying potential prognostic markers at the level of gene expression and miRNAs associated with one of the key signs of cancer aggressiveness-lymphatic dissemination. The relative expression of candidate markers was validated by quantitative PCR, including an independent sample of patients based on archival material. Statistically significant results, derived from an independent set of samples, were confirmed for miR-148a-3p and miR-615-3p, as well as for the CST2, OCLN, and PCAT4 genes. Considering the obtained validation data, we also analyzed the predictive value of models based on various combinations of identified markers using algorithms based on machine learning. The highest predictive potential was shown for the "CST2 + OCLN + pT" model (AUC = 0.863) based on the CatBoost Classifier algorithm.

Transcriptomic Analysis of the Effect of Torin-2 on the Central Nervous System of Drosophila melanogaster.

Torin-2, a synthetic compound, is a highly selective inhibitor of both TORC1 and TORC2 (target of rapamycin) complexes as an alternative to the well-known immunosuppressor, geroprotector, and potential anti-cancer natural compound rapamycin. Torin-2 is effective at hundreds of times lower concentrations and prevents some negative side effects of rapamycin. Moreover, it inhibits the rapamycin-resistant TORC2 complex. In this work, we evaluated transcriptomic changes in D. melanogaster heads induced with lifetime diets containing Torin-2 and suggested possible neuroprotective mechanisms of Torin-2. The analysis included D. melanogaster of three ages (2, 4, and 6 weeks old), separately for males and females. Torin-2, taken at the lowest concentration being tested (0.5 µM per 1 L of nutrient paste), had a slight positive effect on the lifespan of D. melanogaster males (+4% on the average) and no positive effect on females. At the same time, RNA-Seq analysis revealed interesting and previously undiscussed effects of Torin-2, which differed between sexes as well as in flies of different ages. Among the cellular pathways mostly altered by Torin-2 at the gene expression level, we identified immune response, protein folding (heat shock proteins), histone modification, actin cytoskeleton organization, phototransduction and sexual behavior. Additionally, we revealed that Torin-2 predominantly reduced the expression of Srr gene responsible for the conversion of L-serine to D-serine and thus regulating activity of NMDA receptor. Via western blot analysis, we showed than in old males Torin-2 tends to increase the ratio of the active phosphorylated form of ERK, the lowest node of the MAPK cascade, which may play a significant role in neuroprotection. Thus, the complex effect of Torin-2 may be due to the interplay of the immune system, hormonal background, and metabolism. Our work is of interest for further research in the field of NMDA-mediated neurodegeneration.

Key FAD2, FAD3, and SAD Genes Involved in the Fatty Acid Synthesis in Flax Identified Based on Genomic and Transcriptomic Data.

FAD (fatty acid desaturase) and SAD (stearoyl-ACP desaturase) genes play key roles in the synthesis of fatty acids (FA) and determination of oil composition in flax (Linum usitatissimum L.). We searched for FAD and SAD genes in the most widely used flax genome of the variety CDC Bethune and three available long-read assembled flax genomes-YY5, 3896, and Atlant. We identified fifteen FAD2, six FAD3, and four SAD genes. Of all the identified genes, 24 were present in duplicated pairs. In most cases, two genes from a pair differed by a significant number of gene-specific SNPs (single nucleotide polymorphisms) or even InDels (insertions/deletions), except for FAD2a-1 and FAD2a-2, where only seven SNPs distinguished these genes. Errors were detected in the FAD2a-1, FAD2a-2, FAD3c-1, and FAD3d-2 sequences in the CDC Bethune genome assembly but not in the long-read genome assemblies. Expression analysis of the available transcriptomic data for different flax organs/tissues revealed that FAD2a-1, FAD2a-2, FAD3a, FAD3b, SAD3-1, and SAD3-2 were specifically expressed in embryos/seeds/capsules and could play a crucial role in the synthesis of FA in flax seeds. In contrast, FAD2b-1, FAD2b-2, SAD2-1, and SAD2-2 were highly expressed in all analyzed organs/tissues and could be involved in FA synthesis in whole flax plants. FAD2c-2, FAD2d-1, FAD3c-1, FAD3c-2, FAD3d-1, FAD3d-2, SAD3-1, and SAD3-2 showed differential expression under stress conditions-Fusarium oxysporum infection and drought. The obtained results are essential for research on molecular mechanisms of fatty acid synthesis, FAD and SAD editing, and marker-assisted and genomic selection for breeding flax varieties with a determined fatty acid composition of oil.

ALDH3A2, ODF2, QSOX2, and MicroRNA-503-5p Expression to Forecast Recurrence in TMPRSS2-ERG-Positive Prostate Cancer.

Following radical surgery, patients may suffer a relapse. It is important to identify such patients so that therapy tactics can be modified appropriately. Existing stratification schemes do not display the probability of recurrence with enough precision since locally advanced prostate cancer (PCa) is classified as high-risk but is not ranked in greater detail. Between 40 and 50% of PCa cases belong to the TMPRSS2-ERG subtype that is a sufficiently homogeneous group for high-precision prognostic marker search to be possible. This study includes two independent cohorts and is based on high throughput sequencing and qPCR data. As a result, we have been able to suggest a perspective-trained model involving a deep neural network based on both qPCR data for mRNA and miRNA and clinicopathological criteria that can be used for recurrence risk forecasts in patients with TMPRSS2-ERG-positive, locally advanced PCa (the model uses ALDH3A2 + ODF2 + QSOX2 + hsa-miR-503-5p + ISUP + pT, with an AUC = 0.944). In addition to the prognostic model's use of identified differentially expressed genes and miRNAs, miRNA-target pairs were found that correlate with the prognosis and can be presented as an interactome network.

The Effect of Meclofenoxate on the Transcriptome of Aging Brain of Nothobranchius guentheri Annual Killifish.

Annual fish of the genus Nothobranchius are promising models for aging research. Nothobranchius reproduces typical aspects of vertebrate aging, including hallmarks of brain aging. Meclofenoxate (MF) is a well-known compound that can enhance cognitive performance. The drug is prescribed for asthenic conditions, trauma, and vascular diseases of the brain. It is believed that MF is able to delay age-dependent changes in the human brain. However, until now, there has been no study of the MF effect on the brain transcriptome. In the present work, we performed an RNA-Seq study of brain tissues from aged Nothobranchius guentheri, which were almost lifetime administered with MF, as well as young and aged control fish. As expected, in response to MF, we revealed significant overexpression of neuron-specific genes including genes involved in synaptic activity and plasticity, neurotransmitter secretion, and neuron projection. The effect was more pronounced in female fish. In this aspect, MF alleviated age-dependent decreased expression of genes involved in neuronal activity. In both treated and untreated animals, we observed strong aging-associated overexpression of immune and inflammatory response genes. MF treatment did not prevent this effect, and moreover, some of these genes tended to be slightly upregulated under MF treatment. Additionally, we noticed upregulation of some genes associated with aging and cellular senescence, including isoforms of putative vascular cell adhesion molecule 1 (VCAM1), protein O-GlcNAcase (OGA), protein kinase C alpha type (KPCA), prolow-density lipoprotein receptor-related protein 1 (LRP1). Noteworthy, MF treatment was also associated with the elevated transcription of transposons, which are highly abundant in the N. guentheri genome. In conclusion, MF compensates for the age-dependent downregulation of neuronal activity genes, but its effect on aging brain transcriptome still cannot be considered unambiguously positive.

Spontaneous gain of susceptibility suggests a novel mechanism of resistance to hybrid dysgenesis in Drosophila virilis.

Syndromes of hybrid dysgenesis (HD) have been critical for our understanding of the transgenerational maintenance of genome stability by piRNA. HD in D. virilis represents a special case of HD since it includes simultaneous mobilization of a set of TEs that belong to different classes. The standard explanation for HD is that eggs of the responder strains lack an abundant pool of piRNAs corresponding to the asymmetric TE families transmitted solely by sperm. However, there are several strains of D. virilis that lack asymmetric TEs, but exhibit a "neutral" cytotype that confers resistance to HD. To characterize the mechanism of resistance to HD, we performed a comparative analysis of the landscape of ovarian small RNAs in strains that vary in their resistance to HD mediated sterility. We demonstrate that resistance to HD cannot be solely explained by a maternal piRNA pool that matches the assemblage of TEs that likely cause HD. In support of this, we have witnessed a cytotype shift from neutral (N) to susceptible (M) in a strain devoid of all major TEs implicated in HD. This shift occurred in the absence of significant change in TE copy number and expression of piRNAs homologous to asymmetric TEs. Instead, this shift is associated with a change in the chromatin profile of repeat sequences unlikely to be causative of paternal induction. Overall, our data suggest that resistance to TE-mediated sterility during HD may be achieved by mechanisms that are distinct from the canonical syndromes of HD.

Does fresh farmyard manure introduce surviving microbes into soil or activate soil-borne microbiota?

Manure inputs into soil strongly affect soil microbial communities leading to shifts in microbial diversity and activity. It is still not clear whether these effects are caused mainly by the survival of microbes introduced with manure or by activation of the soil-borne microbiome. Here, we investigated how the soil microbiome was changed after the introduction of fresh farmyard cattle manure, and which microorganisms originating from manure survived in soil. Manure addition led to a strong increase in soil microbial biomass, gene copies abundances, respiration activity, and diversity. High-throughput sequencing analysis showed that higher microbial diversity in manured soil was caused mainly by activation of 113 soil-borne microbial genera which were mostly minor taxa in not-fertilized soil. Two weeks after manure input, 78% of the manure-associated genera were not detected anymore. Only 15 of 237 prokaryotic genera that originated from manure survived for 144 days in soil, and only 8 of them (primarily representatives of Clostridia class) were found in manured soil after winter. Thus, an increase in microbial biomass and diversity after manure input is caused mainly by activation of soil-borne microbial communities, while most exogenous microbes from manure do not survive in soil conditions after few months.

Genetic diversity of SAD and FAD genes responsible for the fatty acid composition in flax cultivars and lines.

BACKGROUND: Flax (Linum usitatissimum L.) is grown for fiber and seed in many countries. Flax cultivars differ in the oil composition and, depending on the ratio of fatty acids, are used in pharmaceutical, food, or paint industries. It is known that genes of SAD (stearoyl-ACP desaturase) and FAD (fatty acid desaturase) families play a key role in the synthesis of fatty acids, and some alleles of these genes are associated with a certain composition of flax oil. However, data on genetic polymorphism of these genes are still insufficient. RESULTS: On the basis of the collection of the Institute for Flax (Torzhok, Russia), we formed a representative set of 84 cultivars and lines reflecting the diversity of fatty acid composition of flax oil. An approach for the determination of full-length sequences of SAD1, SAD2, FAD2A, FAD2B, FAD3A, and FAD3B genes using the Illumina platform was developed and deep sequencing of the 6 genes in 84 flax samples was performed on MiSeq. The obtained high coverage (about 400x on average) enabled accurate assessment of polymorphisms in SAD1, SAD2, FAD2A, FAD2B, FAD3A, and FAD3B genes and evaluation of cultivar/line heterogeneity. The highest level of genetic diversity was observed for FAD3A and FAD3B genes - 91 and 62 polymorphisms respectively. Correlation analysis revealed associations between particular variants in SAD and FAD genes and predominantly those fatty acids whose conversion they catalyze: SAD - stearic and oleic acids, FAD2 - oleic and linoleic acids, FAD3 - linoleic and linolenic acids. All except one low-linolenic flax cultivars/lines contained both the substitution of tryptophan to stop codon in the FAD3A gene and histidine to tyrosine substitution in the FAD3B gene, while samples with only one of these polymorphisms had medium content of linolenic acid and cultivars/lines without them were high-linolenic. CONCLUSIONS: Genetic polymorphism of SAD and FAD genes was evaluated in the collection of flax cultivars and lines with diverse oil composition, and associations between particular polymorphisms and the ratio of fatty acids were revealed. The achieved results are the basis for the development of marker-assisted selection and DNA-based certification of flax cultivars.

Immunohistochemistry and Mutation Analysis of SDHx Genes in Carotid Paragangliomas.

Carotid paragangliomas (CPGLs) are rare neuroendocrine tumors often associated with mutations in SDHx genes. The immunohistochemistry of succinate dehydrogenase (SDH) subunits has been considered a useful instrument for the prediction of SDHx mutations in paragangliomas/pheochromocytomas. We compared the mutation status of SDHx genes with the immunohistochemical (IHC) staining of SDH subunits in CPGLs. To identify pathogenic/likely pathogenic variants in SDHx genes, exome sequencing data analysis among 42 CPGL patients was performed. IHC staining of SDH subunits was carried out for all CPGLs studied. We encountered SDHx variants in 38% (16/42) of the cases in SDHx genes. IHC showed negative (5/15) or weak diffuse (10/15) SDHB staining in most tumors with variants in any of SDHx (94%, 15/16). In SDHA-mutated CPGL, SDHA expression was completely absent and weak diffuse SDHB staining was detected. Positive immunoreactivity for all SDH subunits was found in one case with a variant in SDHD. Notably, CPGL samples without variants in SDHx also demonstrated negative (2/11) or weak diffuse (9/11) SDHB staining (42%, 11/26). Obtained results indicate that SDH immunohistochemistry does not fully reflect the presence of mutations in the genes; diagnostic effectiveness of this method was 71%. However, given the high sensitivity of SDHB immunohistochemistry, it could be used for initial identifications of patients potentially carrying SDHx mutations for recommendation of genetic testing.

LINC00973 Induces Proliferation Arrest of Drug-Treated Cancer Cells by Preventing p21 Degradation.

Overcoming drug resistance of cancer cells is the major challenge in molecular oncology. Here, we demonstrate that long non-coding RNA LINC00973 is up-regulated in normal and cancer cells of different origins upon treatment with different chemotherapeutics. Bioinformatics analysis shows that this is a consequence of DNA damage response pathway activation or mitotic arrest. Knockdown of LINC0973 decreases p21 levels, activates cellular proliferation of cancer cells, and suppresses apoptosis of drug-treated cells. We have found that LINC00973 strongly increases p21 protein content, possibly by blocking its degradation. Besides, we have found that ectopic over-expression of LINC00973 inhibits formation of the pro-survival p53-Ser15-P isoform, which preserves chromosome integrity. These results might open a new approach to the development of more efficient anti-cancer drugs.

Characterization of repeated DNA sequences in genomes of blue-flowered flax.

BACKGROUND: Members of different sections of the genus Linum are characterized by wide variability in size, morphology and number of chromosomes in karyotypes. Since such variability is determined mainly by the amount and composition of repeated sequences, we conducted a comparative study of the repeatomes of species from four sections forming a clade of blue-flowered flax. Based on the results of high-throughput genome sequencing performed in this study as well as available WGS data, bioinformatic analyses of repeated sequences from 12 flax samples were carried out using a graph-based clustering method. RESULTS: It was found that the genomes of closely related species, which have a similar karyotype structure, are also similar in the repeatome composition. In contrast, the repeatomes of karyologically distinct species differed significantly, and no similar tandem-organized repeats have been identified in their genomes. At the same time, many common mobile element families have been identified in genomes of all species, among them, Athila Ty3/gypsy LTR retrotransposon was the most abundant. The 30-chromosome members of the sect. Linum (including the cultivated species L. usitatissimum) differed significantly from other studied species by a great number of satellite DNA families as well as their relative content in genomes. CONCLUSIONS: The evolution of studied flax species was accompanied by waves of amplification of satellite DNAs and LTR retrotransposons. The observed inverse correlation between the total contents of dispersed repeats and satellite DNAs allowed to suggest a relationship between both classes of repeating sequences. Significant interspecific differences in satellite DNA sets indicated a high rate of evolution of this genomic fraction. The phylogenetic relationships between the investigated flax species, obtained by comparison of the repeatomes, agreed with the results of previous molecular phylogenetic studies.

Flax (Linum usitatissimum L.) response to non-optimal soil acidity and zinc deficiency.

BACKGROUND: Flax (Linum usitatissimum L.) is grown for fiber and seed production. Unfavorable environments, such as nutrient deficiency and non-optimal soil acidity, decrease the quantity and quality of yield. Cultivation of tolerant to stress varieties can significantly reduce the crop losses. Understanding the mechanisms of flax response to the stresses and identification of resistance gene candidates will help in breeding of improved cultivars. In the present work, the response of flax plants to increased pH level and zinc (Zn) deficiency was studied. RESULTS: We performed high-throughput transcriptome sequencing of two flax cultivars with diverse tolerance to increased pH level and Zn deficiency: Norlin (tolerant) and Mogilevsky (sensitive). Sixteen cDNA libraries were created from flax plants grown under control conditions, increased pH level, Zn deficiency, and both stresses simultaneously, and about 35 million reads were obtained for each experiment type. Unfavorable pH resulted in significantly stronger gene expression alterations compared to Zn deficiency. Ion homeostasis, oxidoreductase activity, cell wall, and response to stress Gene Ontology terms were the most affected by unfavorable pH and Zn deficiency both in tolerant and sensitive flax cultivars. Upregulation of genes encoding metal transporters was identified under increased pH level, Zn deficiency, and both stresses simultaneously. Under Zn deficiency, only in tolerant cultivar Norlin, we revealed the induction of several photosynthesis-related genes and, in this way, this tolerant genotype could overcome unfavorable effects of reduced Zn content. CONCLUSIONS: We identified genes with expression alterations in flax under non-optimal soil acidity and Zn deficiency based on high-throughput sequencing data. These genes are involved in diverse processes, including ion transport, cell wall biogenesis, and photosynthesis, and could play an important role in flax response to the studied stresses. Moreover, genes with distinct expression changes between examined tolerant and sensitive genotypes could determine the mechanisms of flax tolerance to non-optimal soil acidity and Zn deficiency.

Novel potential causative genes in carotid paragangliomas.

BACKGROUND: Carotid paragangliomas (CPGLs) are rare neuroendocrine tumors that arise from the paraganglion at the bifurcation of the carotid artery and are responsible for approximately 65% of all head and neck paragangliomas. CPGLs can occur sporadically or along with different hereditary tumor syndromes. Approximately 30 genes are known to be associated with CPGLs. However, the genetic basis behind the development of these tumors is not fully elucidated, and the molecular mechanisms underlying CPGL pathogenesis remain unclear. METHODS: Whole exome and transcriptome high-throughput sequencing of CPGLs was performed on an Illumina platform. Exome libraries were prepared using a Nextera Rapid Capture Exome Kit (Illumina) and were sequenced under 75 bp paired-end model. For cDNA library preparation, a TruSeq Stranded Total RNA Library Prep Kit with Ribo-Zero Gold (Illumina) was used; transcriptome sequencing was carried out with 100 bp paired-end read length. Obtained data were analyzed using xseq which estimates the influence of mutations on gene expression profiles allowing to identify potential causative genes. RESULTS: We identified a total of 16 candidate genes (MYH15, CSP1, MYH3, PTGES3L, CSGALNACT2, NMD3, IFI44, GMCL1, LSP1, PPFIBP2, RBL2, MAGED1, CNIH3, STRA6, SLC6A13, and ATM) whose variants potentially influence their expression (cis-effect). The strongest cis-effect of loss-of-function variants was found in MYH15, CSP1, and MYH3, and several likely pathogenic variants in these genes associated with CPGLs were predicted. CONCLUSIONS: Using the xseq probabilistic model, three novel potential causative genes, namely MYH15, CSP1, and MYH3, were identified in carotid paragangliomas.

The CIMP-high phenotype is associated with energy metabolism alterations in colon adenocarcinoma.

BACKGROUND: CpG island methylator phenotype (CIMP) is found in 15-20% of malignant colorectal tumors and is characterized by strong CpG hypermethylation over the genome. The molecular mechanisms of this phenomenon are not still fully understood. The development of CIMP is followed by global gene expression alterations and metabolic changes. In particular, CIMP-low colon adenocarcinoma (COAD), predominantly corresponded to consensus molecular subtype 3 (CMS3, "Metabolic") subgroup according to COAD molecular classification, is associated with elevated expression of genes participating in metabolic pathways. METHODS: We performed bioinformatics analysis of RNA-Seq data from The Cancer Genome Atlas (TCGA) project for CIMP-high and non-CIMP COAD samples with DESeq2, clusterProfiler, and topGO R packages. Obtained results were validated on a set of fourteen COAD samples with matched morphologically normal tissues using quantitative PCR (qPCR). RESULTS: Upregulation of multiple genes involved in glycolysis and related processes (ENO2, PFKP, HK3, PKM, ENO1, HK2, PGAM1, GAPDH, ALDOA, GPI, TPI1, and HK1) was revealed in CIMP-high tumors compared to non-CIMP ones. Most remarkably, the expression of the PKLR gene, encoding for pyruvate kinase participating in gluconeogenesis, was decreased approximately 20-fold. Up to 8-fold decrease in the expression of OGDHL gene involved in tricarboxylic acid (TCA) cycle was observed in CIMP-high tumors. Using qPCR, we confirmed the increase (4-fold) in the ENO2 expression and decrease (2-fold) in the OGDHL mRNA level on a set of COAD samples. CONCLUSIONS: We demonstrated the association between CIMP-high status and the energy metabolism changes at the transcriptomic level in colorectal adenocarcinoma against the background of immune pathway activation. Differential methylation of at least nine CpG sites in OGDHL promoter region as well as decreased OGDHL mRNA level can potentially serve as an additional biomarker of the CIMP-high status in COAD.

HK3 overexpression associated with epithelial-mesenchymal transition in colorectal cancer.

BACKGROUND: Colorectal cancer (CRC) is a common cancer worldwide. The main cause of death in CRC includes tumor progression and metastasis. At molecular level, these processes may be triggered by epithelial-mesenchymal transition (EMT) and necessitates specific alterations in cell metabolism. Although several EMT-related metabolic changes have been described in CRC, the mechanism is still poorly understood. RESULTS: Using CrossHub software, we analyzed RNA-Seq expression profile data of CRC derived from The Cancer Genome Atlas (TCGA) project. Correlation analysis between the change in the expression of genes involved in glycolysis and EMT was performed. We obtained the set of genes with significant correlation coefficients, which included 21 EMT-related genes and a single glycolytic gene, HK3. The mRNA level of these genes was measured in 78 paired colorectal cancer samples by quantitative polymerase chain reaction (qPCR). Upregulation of HK3 and deregulation of 11 genes (COL1A1, TWIST1, NFATC1, GLIPR2, SFPR1, FLNA, GREM1, SFRP2, ZEB2, SPP1, and RARRES1) involved in EMT were found. The results of correlation study showed that the expression of HK3 demonstrated a strong correlation with 7 of the 21 examined genes (ZEB2, GREM1, TGFB3, TGFB1, SNAI2, TWIST1, and COL1A1) in CRC. CONCLUSIONS: Upregulation of HK3 is associated with EMT in CRC and may be a crucial metabolic adaptation for rapid proliferation, survival, and metastases of CRC cells.

CrossHub: a tool for multi-way analysis of The Cancer Genome Atlas (TCGA) in the context of gene expression regulation mechanisms.

The contribution of different mechanisms to the regulation of gene expression varies for different tissues and tumors. Complementation of predicted mRNA-miRNA and gene-transcription factor (TF) relationships with the results of expression correlation analyses derived for specific tumor types outlines the interactions with functional impact in the current biomaterial. We developed CrossHub software, which enables two-way identification of most possible TF-gene interactions: on the basis of ENCODE ChIP-Seq binding evidence or Jaspar prediction and co-expression according to the data of The Cancer Genome Atlas (TCGA) project, the largest cancer omics resource. Similarly, CrossHub identifies mRNA-miRNA pairs with predicted or validated binding sites (TargetScan, mirSVR, PicTar, DIANA microT, miRTarBase) and strong negative expression correlations. We observed partial consistency between ChIP-Seq or miRNA target predictions and gene-TF/miRNA co-expression, demonstrating a link between these indicators. Additionally, CrossHub expression-methylation correlation analysis can be used to identify hypermethylated CpG sites or regions with the greatest potential impact on gene expression. Thus, CrossHub is capable of outlining molecular portraits of a specific gene and determining the three most common sources of expression regulation: promoter/enhancer methylation, miRNA interference and TF-mediated activation or repression. CrossHub generates formatted Excel workbooks with the detailed results. CrossHub is freely available athttps://sourceforge.net/projects/crosshub/.

Evolution of blue-flowered species of genus Linum based on high-throughput sequencing of ribosomal RNA genes.

BACKGROUND: The species relationships within the genus Linum have already been studied several times by means of different molecular and phylogenetic approaches. Nevertheless, a number of ambiguities in phylogeny of Linum still remain unresolved. In particular, the species relationships within the sections Stellerolinum and Dasylinum need further clarification. Also, the question of independence of the species of the section Adenolinum still remains unanswered. Moreover, the relationships of L. narbonense and other species of the section Linum require further clarification. Additionally, the origin of tetraploid species of the section Linum (2n = 30) including the cultivated species L. usitatissimum has not been explored. The present study examines the phylogeny of blue-flowered species of Linum by comparisons of 5S rRNA gene sequences as well as ITS1 and ITS2 sequences of 35S rRNA genes. RESULTS: High-throughput sequencing has been used for analysis of multicopy rRNA gene families. In addition to the molecular phylogenetic analysis, the number and chromosomal localization of 5S and 35S rDNA sites has been determined by FISH. Our findings confirm that L. stelleroides forms a basal branch from the clade of blue-flowered flaxes which is independent of the branch formed by species of the sect. Dasylinum. The current molecular phylogenetic approaches, the cytogenetic analysis as well as different genomic DNA fingerprinting methods applied previously did not discriminate certain species within the sect. Adenolinum. The allotetraploid cultivated species L. usitatissimum and its wild ancestor L. angustifolium (2n = 30) could originate either as the result of hybridization of two diploid species (2n = 16) related to the modern L. gandiflorum and L. decumbens, or hybridization of a diploid species (2n = 16) and a diploid ancestor of modern L. narbonense (2n = 14). CONCLUSIONS: High-throughput sequencing of multicopy rRNA gene families allowed us to make several adjustments to the phylogeny of blue-flowered flax species and also reveal intra- and interspecific divergence of the rRNA gene sequences.

Differential gene expression in response to Fusarium oxysporum infection in resistant and susceptible genotypes of flax (Linum usitatissimum L.).

BACKGROUND: Flax (Linum usitatissimum L.) is a crop plant used for fiber and oil production. Although potentially high-yielding flax varieties have been developed, environmental stresses markedly decrease flax production. Among biotic stresses, Fusarium oxysporum f. sp. lini is recognized as one of the most devastating flax pathogens. It causes wilt disease that is one of the major limiting factors for flax production worldwide. Breeding and cultivation of flax varieties resistant to F. oxysporum is the most effective method for controlling wilt disease. Although the mechanisms of flax response to Fusarium have been actively studied, data on the plant response to infection and resistance gene candidates are currently very limited. RESULTS: The transcriptomes of two resistant and two susceptible flax cultivars with respect to Fusarium wilt, as well as two resistant BC2F5 populations, which were grown under control conditions or inoculated with F. oxysporum, were sequenced using the Illumina platform. Genes showing changes in expression under F. oxysporum infection were identified in both resistant and susceptible flax genotypes. We observed the predominant overexpression of numerous genes that are involved in defense response. This was more pronounced in resistant cultivars. In susceptible cultivars, significant downregulation of genes involved in cell wall organization or biogenesis was observed in response to F. oxysporum. In the resistant genotypes, upregulation of genes related to NAD(P)H oxidase activity was detected. Upregulation of a number of genes, including that encoding beta-1,3-glucanase, was significantly greater in the cultivars and BC2F5 populations resistant to Fusarium wilt than in susceptible cultivars in response to F. oxysporum infection. CONCLUSIONS: Using high-throughput sequencing, we identified genes involved in the early defense response of L. usitatissimum against the fungus F. oxysporum. In response to F. oxysporum infection, we detected changes in the expression of pathogenesis-related protein-encoding genes and genes involved in ROS production or related to cell wall biogenesis. Furthermore, we identified genes that were upregulated specifically in flax genotypes resistant to Fusarium wilt. We suggest that the identified genes in resistant cultivars and BC2F5 populations showing induced expression in response to F. oxysporum infection are the most promising resistance gene candidates.