Identification of cell type-specific correlations between ERK activity and cell viability upon treatment with ERK1/2 inhibitors.

Increased MAPK signaling is a hallmark of various cancers and is a central regulator of cell survival. Direct ERK1/2 inhibition is considered a promising approach to avoid ERK1/2 reactivation caused by upstream kinases BRAF, MEK1/2, and KRAS, as well as by receptor tyrosine kinase inhibitors, but the dynamics and selectivity of ERK1/2 inhibitors are much less studied compared with BRAF or MEK inhibitors. Using ERK1/2 and downstream kinase ELK1 reporter cell lines of lung cancer (H1299; NRASQ61K), colon cancer (HCT-116; KRASG13D), neuroblastoma (SH-SY5Y), and leukemia (U937), we examined the relationship between ERK inhibition and drug-induced toxicity for five ERK inhibitors: SCH772984, ravoxertinib, LY3214996, ulixertinib, and VX-11e, as well as one MEK inhibitor, PD0325901. Comparing cell viability and ERK inhibition revealed different ERK dependencies for these cell lines. We identify several drugs, such as SCH772984 and VX-11e, which induce excessive toxicity not directly related to ERK1/2 inhibition in specific cell lines. We also show that PD0325901, LY3214996, and ulixertinib are prone to ERK1/2 reactivation over time. We distinguished two types of ERK1/2 reactivation: the first could be reversed by adding a fresh dose of inhibitors, while the second persists even after additional treatments. We also showed that cells that became resistant to the MEK1/2 inhibitor PD0325901 due to ERK1/2 reactivation remained sensitive to ERK1/2 inhibitor ulixertinib. Our data indicate that correlation of ERK inhibition with drug-induced toxicity in multiple cell lines may help to find more selective and effective ERK1/2 inhibitors.

Comprehensive Genetic Study of Malignant Cervical Paraganglioma.

Malignant middle ear paraganglioma (MEPGL) is an exceedingly rare tumor of the neuroendocrine system. In general, MEPGLs represent as slow growing and hypervascularized benign neoplasms. The genetic basis of MEPGL tumorigenesis has been poorly investigated. We report a case of malignant MEPGL accompanied by the comprehensive genetic analysis of the primary tumor and metastasis. Based on whole-exome sequencing data, the germline pathogenic mutation p.R230H in the SDHB gene, encoding for subunit B of mitochondrial complex II, was found in a patient. Analysis of somatic mutation spectra revealed five novel variants in different genes, including a potentially deleterious variant in UNC13C that was common for the tumor and metastasis. Identified somatic variants clustered into SBS1 and SBS5 mutational signatures. Of note, the primary tumor was characterized by Ki-67 4% and had an elevated mutational load (1.4/Mb); the metastasis' mutational load was about 4.5 times higher (6.4/Mb). In addition, we revealed somatic loss of the wild-type SDHB allele, as well as loss of heterozygosity (LOH) at the 11p locus. Thus, germline mutation in SDHB combined with somatic LOH seem to be drivers that lead to the tumor's initiation and progression. Other somatic changes identified can be additional disease-causing factors. Obtained results expand our understanding of molecular genetic mechanisms associated with the development of this rare tumor.

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.

Complete and Prolonged Inhibition of Herpes Simplex Virus Type 1 Infection In Vitro by CRISPR/Cas9 and CRISPR/CasX Systems.

Almost all people become infected with herpes viruses, including herpes simplex virus type 1 (HSV-1), during their lifetime. Typically, these viruses persist in a latent form that is resistant to all available antiviral medications. Under certain conditions, such as immunosuppression, the latent forms reactivate and cause disease. Moreover, strains of herpesviruses that are drug-resistant have rapidly emerged. Therefore, it is important to develop alternative methods capable of eradicating herpesvirus infections. One promising direction is the development of CRISPR/Cas systems for the therapy of herpesvirus infections. We aimed to design a CRISPR/Cas system for relatively effective long-term and safe control of HSV-1 infection. Here, we show that plasmids encoding the CRISPR/Cas9 system from Streptococcus pyogenes with a single sgRNA targeting the UL30 gene can completely suppress HSV-1 infection of the Vero cell line within 6 days and provide substantial protection within 9 days. For the first time, we show that CRISPR/CasX from Deltaproteobacteria with a single guide RNA against UL30 almost completely suppresses HSV-1 infection of the Vero cell line for 3 days and provides substantial protection for 6 days. We also found that the Cas9 protein without sgRNAs attenuates HSV-1 infection. Our results show that the developed CRISPR/Cas systems are promising therapeutic approaches to control HSV-1 infections.

Docetaxel Resistance in Castration-Resistant Prostate Cancer: Transcriptomic Determinants and the Effect of Inhibiting Wnt/ß-Catenin Signaling by XAV939.

Castration-resistant prostate cancer (CRPC) is a common form of prostate cancer in which docetaxel-based chemotherapy is used as the first line. The present study is devoted to the analysis of transcriptome profiles of tumor cells in the development of resistance to docetaxel as well as to the assessment of the combined effect with the XAV939 tankyrase inhibitor on maintaining the sensitivity of tumor cells to chemotherapy. RNA-Seq was performed for experimental PC3 cell lines as well as for plasma exosome samples from patients with CRPC. We have identified key biological processes and identified a signature based on the expression of 17 mRNA isoforms associated with the development of docetaxel resistance in PC3 cells. Transcripts were found in exosome samples, the increased expression of which was associated with the onset of progression of CRPC during therapy. The suppression of pathways associated with the participation of cellular microtubules has also been shown when cells are treated with docetaxel in the presence of XAV939. These results highlight the importance of further research into XAV939 as a therapeutic agent in the treatment of CRPC; moreover, we have proposed a number of mRNA isoforms with high predictive potential, which can be considered as promising markers of response to docetaxel.

The SDHD:p.H102R Variant Is Frequent in Russian Patients with Head and Neck Paragangliomas and Associated with Loss of 11p15.5 Region and Hypermethylation of H19-DMR.

Head and neck paragangliomas (HNPGLs) are rare neuroendocrine neoplasms derived from the parasympathetic paraganglia of the head and neck. At least 30% of HNPGLs are linked to germline mutations, predominantly in SDHx genes. In this study, we analyzed an extended cohort of Russian patients with HNPGLs using whole-exome sequencing and found a highly frequent missense variant p.H102R in the SDHD gene. We determined this variant in 34% of the SDHD mutation carriers. This variant was associated with somatic loss of the gene wild-type allele. Data from the B allele frequency method and microsatellite and microdeletion analysis indicated evident LOH at the 11p15.5 region and potential loss of the whole of chromosome 11. We found hypermethylation of H19-DMR in all tumors, whereas differential methylation of KvDMR was mostly retained. These findings do not support the paternal transmission of SDHD:p.H102R but are in agreement with the Hensen model. Using targeted sequencing, we also studied the variant frequency in a control cohort; we found SDHD:p.H102R in 1.9% of cases, allowing us to classify this variant as pathogenic. The immunohistochemistry of SDHB showed that the SDHD:p.H102R mutation, even in combination with wild-type allele loss, does not always lead to SDH deficiency. The obtained results demonstrate the frequent variant associated with HNPGLs in a Russian population and support its pathogenicity. Our findings help with understanding the mechanism of tumorigenesis and are also important for the development of cost-effective genetic screening programs.

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.

Case Report: Genetic Alterations Associated with the Progression of Carotid Paraganglioma.

Paragangliomas (PGLs) are rare neuroendocrine tumors that can develop from any paraganglion across the body. The carotid body is the most often location of PGLs in the head and neck region. Carotid PGLs (CPGLs) are characterized by predominantly non-aggressive behavior; however, all tumors have the potential to metastasize. To date, molecular mechanisms of paraganglioma progression remain elusive. We report a case of a 38-year-old woman with metastatic CPGL manifesting as a recurrent tumor with lymph node metastasis. The tumor was fast-growing and had a high Ki-67 proliferation index. Immunohistochemical (IHC) examination and whole-exome sequencing were performed for both recurrent tumor and metastasis. A germline pathogenic splice acceptor variant in the SDHB gene was found in the patient. Immunoreactivity of the SDHB subunit was weak diffuse in both samples, indicating deficiency of the succinate dehydrogenase. Moreover, the recurrent tumor exhibited loss of heterozygosity (LOH) at the SDHB locus, that is according to Knudson's "two-hit" hypothesis of cancer causation. We also identified a rare somatic promotor mutation in the TERT gene associated with the tumor progression. Obtained results confirmed the indicative role of the germline SDHB mutation for metastatic CPGLs, as well as the potential prognostic value of the TERT promoter mutation.

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.

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.

Transcriptome analysis reveals mechanisms of geroprotective effects of fucoxanthin in Drosophila.

BACKGROUND: We have previously showed that the carotenoid fucoxanthin can increase the lifespan in Drosophila melanogaster and Caenorhabditis elegans. However, the molecular mechanisms of the geroprotective effect of fucoxanthin have not been studied so far. RESULTS: Here, we studied the effects of fucoxanthin on the Drosophila aging process at the molecular and the whole organism levels. At the organismal level, fucoxanthin increased the median lifespan and had a positive effect on fecundity, fertility, intestinal barrier function, and nighttime sleep. Transcriptome analysis revealed 57 differentially expressed genes involved in 17 KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways. Among the most represented molecular pathways induced by fucoxanthin, a significant portion is related to longevity, including MAPK, mTOR, Wnt, Notch, and Hippo signaling pathways, autophagy, translation, glycolysis, oxidative phosphorylation, apoptosis, immune response, neurogenesis, sleep, and response to DNA damage. CONCLUSIONS: Life-extending effects of fucoxanthin are associated with differential expression of longevity-associated genes.

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.

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.