Down-regulated MYH11 expression correlates with poor prognosis in stage II and III colorectal cancer

Comprehensive identification of hub mRNAs and lncRNAs in colorectal cancer using galaxy: an in silico transcriptome analysis

Colorectal cancer (CRC) is the second leading cause of cancer-related mortality. Using the Galaxy platform, the present study aimed to assess the differentially expressed genes (DEGs) in CRC patients. The expression data was obtained from the Gene Expression Omnibus database (GSE137327). DEGs were analyzed using Gene Ontology (GO) and GeneMANIA databases to detect the most critical biological pathways and processes. Protein-Protein Interaction Studies (PPIS) identified four hub genes (CCN1, CCL2, FLNC, MYH11). This article presents findings on three mRNAs (CEMIP, MMP7, and DPEP1) and also two notable lncRNAs, EVADR and DLX6-AS1, that have an impact on CRC pathogenesis and play a role in the epithelial-mesenchymal transition in tumor cells. The identified genes and lncRNAs are putative therapeutic targets and diagnostic markers. For instance, CRISPR/Cas9 editing systems can be designed in order to modulate expression of these genes, or edit them for the purpose of inducing sensitivity to conventional therapies. Besides, these genes can be incorporated into clinical prognostic models, offering panels of genes to choose appropriate personalized methods of treatment. Together, these genes represent novel markers and possible therapeutic targets for CRC.

Identification of hub genes and drug candidates for NF2-related vestibular schwannoma by bioinformatics tools

Neurofibromatosis type 2 (NF2)-related vestibular schwannoma (NF2-VS) is a rare genetic disorder that results in bilateral acoustic neuromas. However, the exact pathogenesis of the disease is still unclear. This study aims to use bioinformatics analyses to identify potential hub genes and therapeutic. We retrieved the mRNA expression profiles (GSE108524 and GSE141801) of NF2-VS from the database, and selected the leading 25% genes with the most variance across samples for weighted correlation network analysis. Subsequently, we conducted gene ontology term and Kyoto Encyclopedia of Genes and Genomes signaling network enrichment analyses. The STRING database was employed for protein-protein interaction (PPI) axis construction. The mRNA-miRNA modulatory network was generated via the miRTarBase database. Differentially expressed genes (DEGs) were identified via the R package "limma" in both datasets, and hub genes were screened via intersection of common DEGs, candidate hub genes from the PPI axis, and candidate hub genes from the key module. Finally, common DEGs were uploaded onto the connectivity map database to determine drug candidates. Based on our observations, the blue module exhibited the most significant relation to NF2-VS, and it included the NF2 gene. Using enrichment analysis, we demonstrated that the blue modules were intricately linked to modulations of cell proliferation, migration, adhesion, junction, and actin skeleton. Overall, 356 common DEGs were screened in both datasets, and 33 genes carrying a degree > 15 were chosen as candidate hub genes in the PPI axis. Subsequently, 4 genes, namely, GLUL, CAV1, MYH11, and CCND1 were recognized as real hub genes. In addition, 10 drugs with enrichment scores < -0.7 were identified as drug candidates. Our conclusions offered a novel insight into the potential underlying mechanisms behind NF2-VS. These findings may facilitate the identification of novel therapeutic targets in the future.

Prognostic and immune correlation analysis of mitochondrial autophagy and aging-related genes in lung adenocarcinoma

PURPOSE

Mitophagy and aging (MiAg) are very important pathophysiological mechanisms contributing to tumorigenesis. MiAg-related genes have prognostic value in lung adenocarcinoma (LUAD). However, prognostic, and immune correlation studies of MiAg-related genes in LUAD are lacking.

METHODS

MiAg differentially expressed genes (DEGs) in LUAD were obtained from public sequencing datasets. A prognostic model including MiAg DEGs was constructed according to patients divided into low- and high-risk groups. Gene Ontology, gene set enrichment analysis, gene set variation analysis, CIBERSORT immune infiltration analysis, and clinical characteristic correlation analyses were performed for functional annotation and correlation of MiAgs with prognosis in patients with LUAD.

RESULTS

Seven MiAg DEGs of LUAD were identified: CAV1, DSG2, DSP, MYH11, NME1, PAICS, PLOD2, and the expression levels of these genes were significantly correlated (P < 0.05). The RiskScore of the MiAg DEG prognostic model demonstrated high predictive ability of overall survival of patients diagnosed with LUAD. Patients with high and low MiAg phenotypic scores exhibited significant differences in the infiltration levels of eight types of immune cells (P < 0.05). The multi-factor DEG regression model showed higher efficacy in predicting 5-year survival than 3- and 1-year survival of patients with LUAD.

CONCLUSIONS

Seven MiAg-related genes were identified to be significantly associated with the prognosis of patients diagnosed with LUAD. Moreover, the identified MiAg DEGs might affect the immunotherapy strategy of patients with LUAD.

Proteomic Fingerprint of Lung Fibrosis Progression and Response to Therapy in Bleomycin-Induced Mouse Model

Idiopathic pulmonary fibrosis (IPF) is a chronic lung disease characterized by the aberrant accumulation of extracellular matrix in the lungs. nintedanib is one of the two FDA-approved drugs for IPF treatment; however, the exact pathophysiological mechanisms of fibrosis progression and response to therapy are still poorly understood. In this work, the molecular fingerprint of fibrosis progression and response to nintedanib treatment have been investigated by mass spectrometry-based bottom-up proteomics in paraffin-embedded lung tissues from bleomycin-induced (BLM) pulmonary fibrosis mice. Our proteomics results unveiled that (i) samples clustered depending on the tissue fibrotic grade (mild, moderate, and severe) and not on the time course after BLM treatment; (ii) the dysregulation of different pathways involved in fibrosis progression such as the complement coagulation cascades, advanced glycation end products (AGEs) and their receptors (RAGEs) signaling, the extracellular matrix-receptor interaction, the regulation of actin cytoskeleton, and ribosomes; (iii) Coronin 1A (Coro1a) as the protein with the highest correlation when evaluating the progression of fibrosis, with an increased expression from mild to severe fibrosis; and (iv) a total of 10 differentially expressed proteins (p_adj-value ≤ 0.05 and Fold change ≤-1.5 or ≥1.5), whose abundance varied in the base of the severity of fibrosis (mild and moderate), were modulated by the antifibrotic treatment with nintedanib, reverting their trend. Notably, nintedanib significantly restored lactate dehydrogenase B (Ldhb) expression but not lactate dehydrogenase A (Ldha). Notwithstanding the need for further investigations to validate the roles of both Coro1a and Ldhb, our findings provide an extensive proteomic characterization with a strong relationship with histomorphometric measurements. These results unveil some biological processes in pulmonary fibrosis and drug-mediated fibrosis therapy.

Upregulation of SPI1 in Ectopic Endometrium Contributes to an Invasive Phenotype

BACKGROUD AND AIM

Endometriosis is one of the most common gynecological diseases associated with chronic pelvic pain, infertility, and cancer. However, its molecular pathogenesis remains unclear. This study aimed to identify key genes involved in the pathogenesis of endometriosis.

METHODS

Bioinformatic analyses were perfomed to identify key differentially expressed genes (DEGs), transcription factors (TFs), and functionally enriched pathways. Effect of SPI1 on migration, invasion, expression of ADH1B, MYH11, and PLN were analyzed in human endometriotic cells.

RESULTS

By screening three transcriptome datasets from the GEO for overlapping DEGs between eutopic and ectopic endometria in patients with endometriosis, we found that the expression of ADH1B, MYH11, and PLN was markedly upregulated in the ectopic endometrium. Knockdown of ADH1B, MYH11, and PLN significantly inhibited the migration and invasion of human endometriotic 12Z cells. Additionally, gene set enrichment analysis revealed that epithelial-mesenchymal transition gene signature was positively correlated with ADH1B, MYH11, and PLN expression. Notably, the TF SPI1 was found to regulate the expression of these three genes in the endometriotic tissues and 12TZ cells. Moreover, SPI1 expression was associated with the invasion of endometriotic cells and was increased in the ectopic endometrium of patients with endometriosis.

CONCLUSION

These data suggest that SPI1 plays a key role in the progression of endometriosis by regulating ADH1B, MYH11, and PLN expression and may therefore serve as a potential prognostic and therapeutic factor for endometriosis.

Contribution of the Immune Response in the Ileum to the Development of Diarrhoea caused by Helminth Infection: Studies with the Sheep Model

Gastrointestinal helminths are a global health issue, for humans as well as domestic animals. Most studies focus on the tissues that are infected with the parasite, but here we studied the ileum, a tissue that is rarely infected by helminths. We tested whether inflammation in the ileum contributes to the development and severity of diarrhoea, by comparing sheep that are susceptible (n = 4) or resistant (n = 4) to the disease. We analyzed the ileum transcriptome using RNASeq sequencing approach and various bioinformatics tools including FastQC, STAR, featureCounts, DESeq2, DAVID, clusterProfiler, Cytoscape (ClusterONE) and EnrichR. We identified 243 differentially expressed genes (DEGs), of which 118 were up-regulated and 125 were down-regulated DEGs in the diarrhoea-susceptible animals compared to the diarrhoea-resistant animals. The resulting DEGs were functionally enriched for biological processes, pathways and gene set enrichment analysis. The up-regulated DEGs suggested that an inflammatory immune response was coupled with genes involved in 'Th2 immune response' and 'anti-inflammatory response'. The down-regulated DEGs were related to ion transport, muscle contraction and pathways preventing inflammation. We conclude that i) susceptibility to helminth-induced diarrhoea involves an inflammatory response at a non-infectious site; ii) down-regulation of pathways preventing inflammation can contribute to the severity of diarrhoea; and iii) genes involved in anti-inflammatory responses can reduce the inflammation and diarrhoea.

Integrative Multi-Omics Analysis Reveals Candidate Biomarkers for Oral Squamous Cell Carcinoma

Oral squamous cell carcinoma (OSCC) is one of the most common types of cancer worldwide. Due to the lack of early detection and treatment, the survival rate of OSCC remains poor and the incidence of OSCC has not decreased during the past decades. To explore potential biomarkers and therapeutic targets for OSCC, we analyzed differentially expressed genes (DEGs) associated with OSCC using RNA sequencing technology. Methylation-regulated and differentially expressed genes (MeDEGs) of OSCC were further identified via an integrative approach by examining publicly available methylomic datasets together with our transcriptomic data. Protein-protein interaction (PPI) networks of MeDEGs were constructed and highly connected hub MeDEGs were identified from these PPI networks. Subsequently, expression and survival analyses of hub genes were performed using The Cancer Genome Atlas (TCGA) database and the Gene Expression Profiling Interactive Analysis (GEPIA) online tool. A total of 56 upregulated MeDEGs and 170 downregulated MeDEGs were identified in OSCC. Eleven hub genes with high degree of connectivity were picked out from the PPI networks constructed by those MeDEGs. Among them, the expression level of four hub genes (CTLA4, CDSN, ACTN2, and MYH11) were found to be significantly changed in the head and neck squamous carcinoma (HNSC) patients. Three hypomethylated hub genes (CTLA4, GPR29, and TNFSF11) and one hypermethylated hub gene (ISL1) were found to be significantly associated with overall survival (OS) of HNSC patients. Therefore, these hub genes may serve as potential DNA methylation biomarkers and therapeutic targets of OSCC.

Lysine 2-hydroxyisobutyrylation proteomics reveals protein modification alteration in the actin cytoskeleton pathway of oral squamous cell carcinoma

As the most commonplace malignant carcinoma in the oral cavity, oral squamous cell carcinoma (OSCC) is highly invasive and prone to recurrence. The nosogenesis of OSCC are affected by epigenetics. Recently, a newly-found post-translational modification of lysine, 2-hydroxyisobutylation (Khib), has been proved to play a critical role in biological regulation. However, no research has evaluated the mechanism of Khib in oral cancer. Here, we performed liquid chromatography-mass spectrometry-based quantitative proteomics combined with bioinformatics analysis to reveal and evaluate Khib protein alterations in OSCC. Numerous proteins in OSCC undergo up-regulated modification of Khib. We quantified and identified 967 proteins with differential expression levels, and 617 2-hydroxyisobutylated proteins with 938 Khib sites. Among them, 125 proteins both differentially expressed and accompanied by obvious Khib modification were further identified and analyzed through KEGG-based and ingenuity pathway analysis (IPA). These proteins are enriched in the actin cytoskeleton regulatory pathway, and IPA predicted that they alter the state of actin aggregation and stability, hence impacting and regulating the actin cytoskeleton in OSCC. This is the first 2-hydroxyisobutylated modification proteomics performed for OSCC. Khib protein is significantly concentrated in the actin cytoskeleton regulatory pathway, indicating that this pathway may mediate the tumorigenesis or exacerbation of OSCC. SIGNIFICANCE: This is the first study that revealed the alterations of Khib protein in oral squamous cell carcinoma through LC-MS/MS-based modified proteomic. Our data showed that the protein in the actin cytoskeleton regulatory pathway was underwent significant Khib modification and abundance changes. We applied predictive function in IPA software to analyze and clarify that the aggregation of actin and the regulation of actin stability that mediated by the actin cytoskeleton regulatory pathway may be the potential mechanism of the occurrence and development of oral squamous cell carcinoma. Our research broadens the understanding of the pathogenesis of oral squamous cell carcinoma and provides new insights for future research.

Identification of potential biomarkers with colorectal cancer based on bioinformatics analysis and machine learning

Colorectal cancer (CRC) is one of the most common malignancies worldwide. Biomarker discovery is critical to improve CRC diagnosis, however, machine learning offers a new platform to study the etiology of CRC for this purpose. Therefore, the current study aimed to perform an integrated bioinformatics and machine learning analyses to explore novel biomarkers for CRC prognosis. In this study, we acquired gene expression microarray data from Gene Expression Omnibus (GEO) database. The microarray expressions GSE103512 dataset was downloaded and integrated. Subsequently, differentially expressed genes (DEGs) were identified and functionally analyzed via Gene Ontology (GO) and Kyoto Enrichment of Genes and Genomes (KEGG). Furthermore, protein protein interaction (PPI) network analysis was conducted using the STRING database and Cytoscape software to identify hub genes; however, the hub genes were subjected to Support Vector Machine (SVM), Receiver operating characteristic curve (ROC) and survival analyses to explore their diagnostic values. Meanwhile, TCGA transcriptomics data in Gene Expression Profiling Interactive Analysis (GEPIA) database and the pathology data presented by in the human protein atlas (HPA) database were used to verify our transcriptomic analyses. A total of 105 DEGs were identified in this study. Functional enrichment analysis showed that these genes were significantly enriched in biological processes related to cancer progression. Thereafter, PPI network explored a total of 10 significant hub genes. The ROC curve was used to predict the potential application of biomarkers in CRC diagnosis, with an area under ROC curve (AUC) of these genes exceeding 0.92 suggesting that this risk classifier can discriminate between CRC patients and normal controls. Moreover, the prognostic values of these hub genes were confirmed by survival analyses using different CRC patient cohorts. Our results demonstrated that these 10 differentially expressed hub genes could be used as potential biomarkers for CRC diagnosis.

Interaction between DNMT3B and MYH11 via hypermethylation regulates gastric cancer progression

BACKGROUND

Gastric cancer (GC) has an unwelcoming prognosis when diagnosed at an advanced stage. The purpose of this study was to examine the expression of myosin heavy chain 11 (MYH11) in GC and mechanisms related.

METHODS

The MYH11 expression in GC was investigated via the SangerBox platform. MYH11 expression in GC tissues and cell lines was examined by immunohistochemistry, RT-qPCR, and western blot. The relationship between MYH11 expression and patients' prognosis was analyzed. The effects of MYH11 on the biological behaviors of GC cells were investigated by gain-of-function experiments. Bioinformatics analysis was used to find genes with relevance to MYH11 expression in GC. The relationship was verified by luciferase and ChIP-qPCR assays, followed by rescue assay validation. The causes of MYH11 downregulation in GC were verified by quantitative methylation-specific PCR. Finally, the effect of MYH11 on tumor growth was examined.

RESULTS

MYH11 was downregulated in GC and predicted poor prognoses. MYH11 reverted the malignant phenotype of GC cells. MYH11 repressed the TNFRSF14 expression by binding to the TNFRSF14 promoter. TNFRSF14 reversed the inhibitory effect of MYH11 on the malignant phenotype of GC cells. The methylation of the MYH11 promoter was elevated in GC, which was correlated with the elevated DNMT3B in GC. Overexpression of DNMT3B repressed transcription of MYH11 by promoting its methylation. Also, MYH11 upregulation inhibited tumor growth.

CONCLUSION

DNMT3B inhibits MYH11 expression by promoting its DNA methylation, thereby attenuating the repressive effect of MYH11 on the transcriptional of TNFRSF14 and promoting the progression of GC.

Differentiation between thyroid-associated orbitopathy and Graves' disease by iTRAQ-based quantitative proteomic analysis

Graves' ophthalmopathy, also known as thyroid-associated orbitopathy (TAO), is the most common inflammatory eye disease in adults. The most common etiology for TAO is Graves' disease (GD); however, proteomic research focusing on differences between GD and TAO is limited. This study aimed to identify differentially expressed proteins between thyroid-associated orbitopathy (TAO) and GD. Furthermore, we sought to explore the pathogenesis of TAO and elucidate the differentiation process via specific markers. Serum samples of three patients with TAO, GD, and healthy controls, respectively, were collected. These samples were measured using the iTRAQ technique coupled with mass spectrometry. Differentially expressed proteins in TAO and GD were identified by proteomics; 3172 quantified proteins were identified. Compared with TAO, we identified 110 differential proteins (27 proteins were upregulated and 83 were downregulated). In addition, these differentially expressed proteins were closely associated with cellular processes, metabolic processes, macromolecular complexes, signal transduction, and the immune system. The corresponding functions were protein, calcium ion, and nucleic acid binding. Among the differential proteins, MYH11, P4HB, and C4A were markedly upregulated in TAO patients and have been reported to participate in apoptosis, autophagy, the inflammatory response, and the immune system. A protein-protein interaction network analysis was performed. Proteomics demonstrated valuable large-scale protein-related information for expounding the pathogenic mechanism underlying TAO. This research provides new insights and potential targets for studying GD with TAO.

LncRNA Sox2ot modulates the progression of thoracic aortic aneurysm by regulating miR-330-5p/Myh11

Thoracic aortic aneurysm (TAA) has been causing the death of elder people. Myosin heavy chain 11 (Myh11) has been reported associated with aortic aneurysm, but there is no specific study on its function on TAA. Here we aimed to explore the function of Myh11 on mouse aortic smooth muscle cells (SMCs) for studying the inner mechanism of TAA. H₂O₂ treatment was implemented on mouse aortic SMCs for detecting cell apoptosis. Meanwhile, functional assays were conducted to verify the function of Myh11 on mouse aortic SMCs. Also, pull-down assay, RIP assay were implemented to identify the potential RNAs for study. Quantitative real-time polymerase chain reaction (qRT-PCR) and luciferase reporter assay were implemented to identify the expression and binding relationships of RNAs. Myh11 expression was increased by treatment of H₂O₂. Myh11 could decrease proliferation and enhance apoptosis of mouse aortic SMCs. At the same time, mmu-miR-330-5p could bind to Myh11 and Sox2ot, forming a competing endogenous RNA (ceRNA) pathway to regulate the proliferation and apoptosis of mouse aortic SMCs. Moreover, both Sox2ot and Myh11 were proved to be up-regulated whereas miR-330-5p down-regulated in Fbn1^C1039G/+ mice, the in vivo model of TAA. In a word, long noncoding RNA (lncRNA) Sox2ot modulates the progression of TAA by regulating miR-330-5p/Myh11 axis.

Proteomics study on the effect of silybin on cardiomyopathy in obese mice

Due to the increase in the number of obese individuals, the incidence of obesity-related complications such as cardiovascular disease and type 2 diabetes is higher. The aim of the present study was to explore the effects of silybin on protein expression in obese mice. Firstly, serum was collected, and it was used to detect serum lipids and other serological indicators. Secondly, total protein from epididymal adipose tissue was extracted for differential expression analysis by quantitative tandem mass tag (TMT) combined with liquid chromatography-tandem mass spectrometry (LC-MS/MS), followed by bioinformatics and protein-protein interaction (PPI) network analyses of these proteins. Lastly, real-time polymerase chain reaction (RT-PCR) and parallel reaction monitoring (PRM) were used to further validate the expression of identified differentially expressed proteins (DEPs) at the mRNA and protein level, respectively. The results revealed that silybin could improve abnormal lipid metabolism caused by the high fat diet in obese mice. A total of 341, 538 and 243 DEPs were found in the high fat/control (WF/WC), silybin/high fat (WS/WF) and WS/WC groups, respectively. These DEPs mainly participated in lipid metabolism and energy metabolism. Notably, tropomyosin 1 (TPM1), myosin light chain 2 (MYL2), myosin heavy chain 11 (MYH11) and other DEPs were involved in hypertrophic cardiomyopathy, dilated cardiomyopathy and other pathways. Silybin could protect cardiac function by inducing the protein expression of TPM1, MYL2 and MYH11 in the adipose tissue of obese mice.

Integrated analysis of label-free quantitative proteomics and bioinformatics reveal insights into signaling pathways in male breast cancer

Male breast cancer (MBC) is a rare malignancy that accounts for about 1.8% of all breast cancer cases. In contrast to the high number of the “omics” studies in breast cancer in women, only recently molecular approaches have been performed in MBC research. High-throughput proteomics based methodologies are promisor strategies to characterize the MBC proteomic signatures and their association with clinico-pathological parameters. In this study, the label-free quantification-mass spectrometry and bioinformatics approaches were applied to analyze the proteomic profiling of a MBC case using the primary breast tumor and the corresponding axillary metastatic lymph nodes and adjacent non-tumor breast tissues. The differentially expressed proteins were identified in the signaling pathways of granzyme B, sirtuins, eIF2, actin cytoskeleton, eNOS, acute phase response and calcium and were connected to the upstream regulators MYC, PI3K SMARCA4 and cancer-related chemical drugs. An additional proteomic comparative analysis was performed with a primary breast tumor of a female patient and revealed an interesting set of proteins, which were mainly involved in cancer biology. Together, our data provide a relevant data source for the MBC research that can help the therapeutic strategies for its management.

Construction and prognostic analysis of miRNA-mRNA regulatory network in liver metastasis from colorectal cancer

Background

Colorectal cancer (CRC) is one of the most common cancers in the world, and liver metastasis is the leading cause of colorectal cancer-related deaths. However, the mechanism of liver metastasis in CRC has not been clearly elucidated.

Methods

Three datasets from the Gene Expression Omnibus (GEO) database were analyzed to obtain differentially expressed genes (DEGs), which were subjected to functional enrichment analysis and protein-protein interaction analysis. Subsequently, mRNA-miRNA network was constructed, and the associated DEGs and DEMs were performed for prognostic analysis. Finally, we did infiltration analysis of growth arrest specific 1 (GAS1)-associated immune cells.

Results

We obtained 325 DEGs and 9 differentially expressed miRNAs (DEMs) between primary CRC and liver metastases. Enrichment analysis and protein-protein interactions (PPI) further revealed the involvement of DEGs in the formation of the inflammatory microenvironment and epithelial-mesenchymal transition (EMT) during the liver metastases process in CRC. Survival analysis demonstrated that low-expressed GAS1 as well as low-expressed hsa-miR-33b-5p was a favorable prognostic indicator of overall survival. Further exploration of GAS1 revealed that its expression was interrelated with the infiltration of immune cells in tumor tissues.

Conclusions

In summary, DEGs, DEMs, and their interactions found in liver metastasis of CRC may provide a basis for further understanding of the mechanism of CRC metastasis.

Myosins and Disease

Myosins constitute a superfamily of actin-based molecular motor proteins that mediates a variety of cellular activities including muscle contraction, cell migration, intracellular transport, the formation of membrane projections, cell adhesion, and cell signaling. The 12 myosin classes that are expressed in humans share sequence similarities especially in the N-terminal motor domain; however, their enzymatic activities, regulation, ability to dimerize, binding partners, and cellular functions differ. It is becoming increasingly apparent that defects in myosins are associated with diseases including cardiomyopathies, colitis, glomerulosclerosis, neurological defects, cancer, blindness, and deafness. Here, we review the current state of knowledge regarding myosins and disease.

Identification of potential crucial genes associated with the pathogenesis and prognosis of prostate cancer

Aim: Prostate cancer (PCa) is the sixth leading cause of cancer-related deaths in men throughout the world. This study aimed to investigate genes associated with the pathogenesis and prognosis of PCa. Materials & methods: Data of PCa cases were obtained from public datasets and were analyzed using an integrated bioinformatics strategy. Results: A total of 969 differential expression genes were identified. Moreover, GSE16560 and The Cancer Genome Atlas (TCGA) data showed a prognostic prompt function of the nine-gene signature, as well as in PCa with Gleason 7. Finally, majority of the nine hub genes were associated with drug sensitivity, mutational landscape, immune infiltrates and clinical characteristics of PCa. Conclusion: The nine-gene signature was correlated with drug sensitivity, mutational landscape, immune infiltrates, clinical characteristics and survival from PCa.

Clinical and prognostic significance of MYH11 in lung cancer

Myosin heavy chain 11 (MYH11), encoded by the MYH11 gene, is a protein that participates in muscle contraction through the hydrolysis of adenosine triphosphate. Although previous studies have demonstrated that MYH11 gene expression levels are downregulated in several types of cancer, its expression levels have rarely been investigated in lung cancer. The present study aimed to explore the clinical significance and prognostic value of MYH11 expression levels in lung cancer and to further study the underlying molecular mechanisms of the function of this gene. The Oncomine database showed that the MYH11 expression levels were decreased in lung cancer compared with those noted in the normal lung tissue (P<0.05). Kaplan-Meier plotter results revealed that the decreased MYH11 expression levels were correlated with poor prognosis in lung cancer patients. Among the lung cancer cases with gene alteration of MYH11, mutation was the most common of all alteration types. Coexpedia and Metascape analyses revealed that the target genes were primarily enriched in ‘muscle contraction’, ‘contractile fiber part’, ‘actin cytoskeleton’ and the ‘adherens junction’. These results indicated that MYH11 is a potential novel drug target and prognostic indicator of lung cancer.

TMT-based proteomics analysis of the anti-hepatocellular carcinoma effect of combined dihydroartemisinin and sorafenib

Hepatocellular carcinoma (HCC), as the major primary liver cancer, is one of the most prevalent malignant diseases with a high mortality rate worldwide. Prior studies have demonstrated that dihydroartemisinin (DHA), the semisynthetic derivative of artemisinin, possesses anti-HCC activity. The multikinase inhibitor sorafenib has been approved for the treatment of HCC. However, the anti-HCC efficacy of DHA combined with sorafenib has not been reported. In this study, we confirmed the significantly enhanced anti-HCC efficacy of DHA in combination with sorafenib compared with that of each agent alone. Tandem Mass Tag (TMT) peptide labeling coupled with LC-MS/MS was used to quantify the proteins from the control, DHA, sorafenib, and DHA + sorafenib groups. In total, 532, 426, 628 differentially expressed proteins (fold change >1.20 or <0.83 and P-value <0.05) were determined by comparing DHA versus control, sorafenib versus control and DHA + sorafenib versus control groups, respectively. Moreover, optimized screening was performed, and 101 optimized differentially expressed proteins were identified. The results of functional analysis of the optimized differentially expressed proteins suggested that they were enriched in cell components such as membrane-bound vesicles, extracellular vesicles, and organelle lumens, and they were mainly involved in biological processes such as cellular component organization, response to stress, and response to chemicals; in addition, they were related to various molecular functions such as protein binding, chromatin binding and enzyme binding. KEGG pathway analysis showed that the optimized differentially expressed proteins were enriched in pyrimidine metabolism, RNA polymerase, base excision repair, and osteoclast differentiation. Protein-protein interaction (PPI) networks of some of the optimized upregulated proteins suggested that they might not only affect vitamin and fat digestion and absorption but may also be involved in tight junctions. In the PPI network, some of the optimized downregulated proteins were enriched in base excision repair, RNA polymerase, purine metabolism, pyrimidine metabolism and mucin type O-glycan biosynthesis. Overall, this research explored the anti-HCC efficacy of DHA combined with sorafenib by using the TMT-based quantitative proteomics technique and might facilitate the understanding of the related anti-HCC molecular mechanism.

Compound heterozygous variants in MYH11 underlie autosomal recessive megacystis-microcolon-intestinal hypoperistalsis syndrome in a Chinese family

Megacystis-microcolon-intestinal-hypoperistalsis syndrome (MMIHS) is a rare and severe disorder characterized by functional obstruction in the urinary and gastrointestinal tract. The molecular basis of this condition has been defined recently. Heterozygous variants in ACTG2, homozygous mutations in LMOD1, MYLK, and MYH9 were related to the pathogenesis of the syndrome, which encodes proteins involved in the process of smooth muscle contraction, supporting a myopathic basis for the disease. Recent studies have identified homozygous or compound heterozygous variants in MYH11 as a candidate gene of MMIHS. In this report, we described a nonconsanguineous Chinese family with three male fetuses affected with megacystis. Trio-targeted exome sequencing identified compound heterozygous variants, c.2051 G > A (p.R684H) and c.3540_3541delinsTT (p.(E1180D, Q1181Ter)), in MYH11 (NM_001040114). The variants were inherited from the parents, respectively. Western blotting showed a marked decrease in MYH11 protein in the proband's umbilical cord tissue compared with the control sample. The study's results confirmed that MYH11 is a candidate gene for MMIHS with autosomal recessive (AR) inheritance and expanded the mutation spectrum for this clinical condition. Combining clinical phenotype with molecular diagnosis may enable the identification of candidate genes for potential monogenic diseases and facilitate accurate genetic counseling, informed decision-making, and prenatal diagnosis.

lncRNA SNHG1 cooperated with miR-497/miR-195-5p to modify epithelial-mesenchymal transition underlying colorectal cancer exacerbation

Our study was intended to provide evidence for whether long noncoding RNA (lncRNA) SNHG1 would accelerate the epithelial-mesenchymal transition (EMT) course intrinsic in colorectal cancer (CRC) by sponging downstream miR-497-5p and miR-195-5p. We altogether collected 338 pairs of CRC and noncancerous tissues, and meanwhile purchased five CRC cell lines (i.e., SW480, HCT116, Lovo, CaCO-2, and HT29) and human embryo intestinal mucosal tissue-sourced cell line (i.e., CCC-HIE-2). The CRC cells as mentioned above were appraised regarding their potencies in proliferation, migration, and invasion, after being transfected with pcDNA3.1-SNHG1, si-SNHG1, miR-195-5p mimic/inhibitor, and miR-497-5p mimic/inhibitor. Eventually, we depended on reverse transcription-polymerase chain reaction to assess SNHG1, miR-497-5p, and miR-195-5p expressions, and the protein levels of EMT-specific molecules were determined on the strength of western blotting. It seemed that there was a high potential for highly expressed SNHG1 and lowly expressed miR-497/miR-195 to symbolize CRC patients' unfavorable prognosis (p < .05). Concurrently, CRC cells were detected with higher SNHG1 expression and lower miR-497/miR-195 expression than CCC-HIE-2 cells (p < .05). In addition, the EMT process of CRC cells was facilitated markedly against the contexts of overexpressed SNHG1 and underexpressed miR-497-5p/miR-195-5p. Intriguingly, the strength of miR-195-5p collaborating with miR-497-5p in affecting the activity of CRC cells seemed to overweigh that of miR-497/miR-195-5p alone. Besides, both miR-195-5p and miR-497-5p were subjected to in vivo and in vitro modification of SNHG1 (p < .05). Conclusively, application of lncRNA SNHG1 for treating CRC might be promising, given its dual modulation of miR-497 and miR-195 underlying CRC pathogenesis.

Identification of LIFR, PIK3R1, and MMP12 as Novel Prognostic Signatures in Gallbladder Cancer Using Network-Based Module Analysis

Background: Gallbladder cancer (GBC) is a rare and aggressive malignancy of the biliary tract with a dismal survival rate. Effective biomarkers and therapeutic targets are urgently needed.

Methods: We analyzed gene expression profiles of GBC to identify differentially expressed genes (DEGs) and then used these DEGs to identify functional module biomarkers based on protein functional interaction (FI) networks. We further evaluated the module-gene protein expression and clinical significance with immunohistochemistry staining (IHC) in a tissue microarray (TMA) from 80 GBC samples.

Results: Five functional modules were identified. Module 0 included classical cancer signaling pathways, such as Ras and PI3K-Akt; and modules 1–4 included genes associated with muscle cells, fibrinogen, extracellular matrix, and integrins, respectively. We validated the expression of LIFR, PIK3R1, and MMP12, which were hubs or functional nodes in modules. Compared with paired peritumoural tissues, we found that the expression of LIFR (P = 0.002) and PIK3R1 (P = 0.046) proteins were significantly downregulated, and MMP12 (P = 0.006) was significantly upregulated. Further prognostic analysis showed that patients with low expression of LIFR had shorter overall survival than those with high expression (log-rank test P = 0.028), the same trend as for PIK3R1 (P = 0.053) and MMP12 (P = 0.006). Multivariate analysis indicated that expression of MMP12 protein (hazard ratio [HR] = 0.429; 95% confidence interval [CI] 0.198, 0.930; P = 0.032) was one of the significant independent prognostic factors for overall survival.

Conclusions: We found a highly reliable FI network, which revealed LIFR, PIK3R1, and MMP12 as novel prognostic biomarker candidates for GBC. These findings could accelerate biomarker discovery and therapeutic development in this cancer.

Atractylenolide II reverses the influence of lncRNA XIST/miR-30a-5p/ROR1 axis on chemo-resistance of colorectal cancer cells

This investigation was conducted to elucidate whether atractylenolide II could reverse the role of lncRNA XIST/miR-30a-5p/ROR1 axis in modulating chemosensitivity of colorectal cancer cells. We totally collected 294 pairs of colorectal cancer tissues and adjacent normal tissues and also purchased colorectal cancer cell lines and human embryonic kidney cell line. 5-fluorouracil, cisplatin, mitomycin and adriamycin were designated as the chemotherapies for colorectal cell lines, and atractylenolides were arranged as the Chinese drug. The expressions of XIST, miR-30a-5p and ROR1 were quantified with aid of qRT-PCR or Western blot, and luciferase reporter gene assay was implemented to determine the relationships among XIST, miR-30a-5p and ROR1. Our results demonstrated that XIST and ROR1 expressions were dramatically up-regulated, yet miR-30a-5p expression was down-regulated within colorectal cancer tissues (P < 0.05). The overexpressed XIST and ROR1, as well as under-expressed miR-30a-5p, were inclined to promote viability and proliferation of colorectal cells under the influence of chemo drugs (P < 0.05). In addition, XIST could directly target miR-30a-5p, and ROR1 acted as the targeted molecule of miR-30a-5p. Interestingly, atractylenolides not only switched the expressions of XIST, miR-30a-5p and ROR1 within colorectal cancer cells but also significantly intensified the chemosensitivity of colorectal cancer cells (P < 0.05). Finally, atractylenolide II was discovered to slow down the viability and proliferation of colorectal cancer cells (P < 0.05). In conclusion, the XIST/miR-30a-5p/ROR1 axis could be deemed as pivotal markers underlying colorectal cancer, and administration of atractylenolide II might improve the chemotherapeutic efficacy for colorectal cancer.

Characterization of DNA Methylation Associated Gene Regulatory Networks During Stomach Cancer Progression

DNA methylation plays a critical role in tumorigenesis through regulating oncogene activation and tumor suppressor gene silencing. Although extensively analyzed, the implication of DNA methylation in gene regulatory network is less characterized. To address this issue, in this study we performed an integrative analysis on the alteration of DNA methylation patterns and the dynamics of gene regulatory network topology across distinct stages of stomach cancer. We found the global DNA methylation patterns in different stages are generally conserved, whereas some significantly differentially methylated genes were exclusively observed in the early stage of stomach cancer. Integrative analysis of DNA methylation and network topology alteration yielded several genes which have been reported to be involved in the progression of stomach cancer, such as IGF2, ERBB2, GSTP1, MYH11, TMEM59, and SST. Finally, we demonstrated that inhibition of SST promotes cell proliferation, suggesting that DNA methylation-associated SST suppression possibly contributes to the gastric cancer progression. Taken together, our study suggests the DNA methylation-associated regulatory network analysis could be used for identifying cancer-related genes. This strategy can facilitate the understanding of gene regulatory network in cancer biology and provide a new insight into the study of DNA methylation at system level.

Identification of Potential Key Genes and Pathways in Early-Onset Colorectal Cancer Through Bioinformatics Analysis

This study was designed to identify the potential key protein interaction networks, genes, and correlated pathways in early-onset colorectal cancer (CRC) via bioinformatics methods. We selected microarray data GSE4107 consisting 12 patient's colonic mucosa and 10 healthy control mucosa; initially, the GSE4107 were downloaded and analyzed using limma package to identify differentially expressed genes (DEGs). A total of 131 DEGs consisting of 108 upregulated genes and 23 downregulated genes of patients in early-onset CRC were selected by the criteria of adjusted P values <.01 and |log2 fold change (FC)| ≥ 2. The gene ontology functional enrichment analysis and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were accomplished to view the biological process, cellular components, molecular function, and the KEGG pathways of DEGs. Finally, protein-protein interactions (PPIs) were constructed, and the hub protein module was identified. Genes such as ACTA2, ACTG2, MYH11, CALD1, MYL9, TPM2, and LMOD1 were strongly implicated in CRC. In summary, in this study, we indicated that molecular mechanisms were involved in muscle contraction and vascular smooth muscle contraction signaling pathway, which improve our understanding of CRC and could be used as new therapeutic targets for CRC.

Pan-cancer analysis for studying cancer stage using protein and gene expression data

Pan-cancer analyses attempt to discover similar features among multiple cancers to identify fundamental patterns common to cancer development and progression. A pan-cancer analysis integrating both protein expression and transcriptomic data is important because it can identify genes that are linked to proteins potentially responsible for a patient's status. This study aims to identify differentially expressed (DE) genes between early and advanced cases of multiple cancer types through the usage of RNA sequencing data. The relevance of these genes is further investigated by developing predictive models using K-nearest neighbor and linear discriminant analysis classifiers. The use of cancer-specific and non-cancer specific features resulted in several moderately performing models. Highlighted genes were further investigated to determine if they encoded for proteins identified in a previously conducted pan-cancer analysis. The results of this study suggest that a pan-cancer analysis may be highly complementary to standard analyses of individual cancers for identifying biologically relevant DE genes and can assist in developing effective predictive models for cancer progression.

The up-regulation of Myb may help mediate EGCG inhibition effect on mouse lung adenocarcinoma

Background

Green tea polyphenol epigallocatechin-3-gallate (EGCG) has been demonstrated to inhibit cancer in experimental studies through its antioxidant activity and modulations on cellular functions by binding specific proteins. By means of computational analysis and functional genomic approaches, we previously identified a set of protein coding genes and microRNAs whose expressions were significantly modulated in response to the EGCG treatment in tobacco carcinogen-induced lung adenocarcinoma in A/J mice. However, to what degree these genes are involved in the cancer inhibition of EGCG remains unclear.

Results

In this study, we further employed statistical methods and literature research to analyze these data in combination with The Cancer Genome Atlas (TCGA) lung adenocarcinoma datasets for additional data mining. Under the assumption that, if a gene mediates EGCG’s cancer inhibition, its expression level change caused by EGCG should be opposite to what occurred in the carcinogenesis, we identified Myb and Peg3 as the primary putative genes involved in the cancer inhibitory activity. Further analysis suggested that the regulation of Myb could be mediated through an EGCG-upregulated microRNA, miR-449c-5p.

Conclusions

Although the actions of EGCG involve multiple targets/pathways, further analysis by mining the existing genomic datasets revealed that the upregulations of Myb and Peg3 are likely the key anti-cancer events of EGCG in vivo.

Electronic supplementary material

The online version of this article (doi:10.1186/s40246-016-0072-4) contains supplementary material, which is available to authorized users.

Multiple region whole-exome sequencing reveals dramatically evolving intratumor genomic heterogeneity in esophageal squamous cell carcinoma

Cancer is a disease of genome instability and genomic alterations; now, genomic heterogeneity is rapidly emerging as a defining feature of cancer, both within and between tumors. Motivation for our pilot study of tumor heterogeneity in esophageal squamous cell carcinoma (ESCC) is that it is not well studied, but the highest incidences of esophageal cancers are found in China and ESCC is the most common type. We profiled the mutations and changes in copy number that were identified by whole-exome sequencing and array-based comparative genomic hybridization in multiple regions within an ESCC from two patients. The average mutational heterogeneity rate was 90% in all regions of the individual tumors in each patient; most somatic point mutations were nonsynonymous substitutions, small Indels occurred in untranslated regions of genes, and copy number alterations varied among multiple regions of a tumor. Independent Sanger sequencing technology confirmed selected gene mutations with more than 88% concordance. Phylogenetic analysis of the somatic mutation frequency demonstrated that multiple, genomically heterogeneous divergent clones evolve and co-exist within a primary ESCC and metastatic subclones result from the dispersal and adaptation of an initially non-metastatic parental clone. Therefore, a single-region sampling will not reflect the evolving architecture of a genomically heterogeneous landscape of mutations in ESCC tumors and the divergent complexity of this genomic heterogeneity among patients will complicate any promise of a simple genetic or epigenetic diagnostic signature in ESCC. We conclude that any potential for informative biomarker discovery in ESCC and targeted personalized therapies will require a deeper understanding of the functional biology of the ontogeny and phylogeny of the tumor heterogeneity.

Characterization of the Transcriptional Complexity of the Receptive and Pre-receptive Endometria of Dairy Goats

Endometrium receptivity is essential for successful embryo implantation in mammals. However, the lack of genetic information remains an obstacle to understanding the mechanisms underlying the development of a receptive endometrium from the pre-receptive phase in dairy goats. In this study, more than 4 billion high-quality reads were generated and de novo assembled into 102,441 unigenes; these unigenes were annotated using published databases. A total of 3,255 unigenes that were differentially expressed (DEGs) between the PE and RE were discovered in this study (P-values < 0.05). In addition, 76,729–77,102 putative SNPs and 12,837 SSRs were discovered in this study. Bioinformatics analysis of the DEGs revealed a number of biological processes and pathways that are potentially involved in the establishment of the RE, notably including the GO terms proteolysis, apoptosis, and cell adhesion and the KEGG pathways Cell cycle and extracellular matrix (ECM)-receptor interaction. We speculated that ADCY8, VCAN, SPOCK1, THBS1, and THBS2 may play important roles in the development of endometrial receptivity. The de novo assembly provided a good starting point and will serve as a valuable resource for further investigations into endometrium receptivity in dairy goats and future studies on the genomes of goats and other related mammals.