Rong-Quan He

Synergistic Effect of MiR-146a Mimic and Cetuximab on Hepatocellular Carcinoma Cells

Previously, we found that the expression of microRNA-146a (miR-146a) was downregulated in hepatocellular carcinoma (HCC) formalin-fixed paraffin-embedded (FFPE) tissues compared to the adjacent noncancerous hepatic tissues. In the current study, we have explored the in vitro effect of miR-146a on the malignant phenotypes of HCC cells. MiR-146a mimic could suppress cell growth and increase cellular apoptosis in HCC cell lines HepG2, HepB3, and SNU449, as assessed by spectrophotometry, fluorimetry, and fluorescence microscopy, respectively. Furthermore, western blot showed that miR-146a mimic downregulated EGFR, ERK1/2, and stat5 signalings. These effects were less potent compared to that of a siRNA targeting EGFR, a known target gene of miR-146a. Moreover, miR-146a mimic could enhance the cell growth inhibition and apoptosis induction impact of various EGFR targeting agents. The most potent combination was miR-146a mimic with cetuximab, presenting a synergistic effect. In conclusion, miR-146a plays a vital role in the cell growth and apoptosis of HCC cells and inducing miR-146a level might be a critical targeted molecular therapy strategy for HCC.

Clinical roles of the aberrantly expressed lncRNAs in lung squamous cell carcinoma: a study based on RNA-sequencing and microarray data mining

Lung squamous cell carcinoma (LUSC) accounts for a significant proportion of lung cancer and there have been few therapeutic alternatives for recurrent LUSC due to the lack of specific driver molecules. To investigate the prospective role of lncRNAs in the tumorigenesis and progression of LUSC, the aberrantly expressed lncRNAs were calculated based on The Cancer Genome Atlas RNA-seq data. Of 7589 lncRNAs with 504 LUSC cases, 884 lncRNAs were identified as being aberrantly expressed (|log2 fold change| >2 and adjusted P<0.05) by DESeq R. The top 10 lncRNAs with the highest diagnostic value were SFTA1P,LINC00968, LINC00961, LINC01572,RP1-78O14.1, FENDRR, LINC01314,LINC01272, GATA6-AS1, and MIR3945HG. In addition to the significant roles in the carcinogenesis of LUSC, several lncRNAs also played vital parts in the survival and progression of LUSC. SFTA1P, LINC01272, GATA6-AS1 and MIR3945HG were closely related to the survival time of LUSC. Furthermore, LINC01572 and LINC01314 could distinguish the LUSC at early stage from that at advanced stage. The prospective molecular assessment of key lncRNAs showed that a certain series of genes could be involved in the regulation network. Furthermore, the OncoPrint from cBioPortal indicated that 14% (69/501) LUSC cases with genetic alterations could be obtained, including amplification, deep deletion and mRNA upregulation. More interestingly, the cases with genetic alterations had a poorer survival as compared to those without alterations. Overall, the study propounds a potentiality for interpreting the pathogenesis and development of LUSC with lncRNAs, and provides a novel platform for searching for more capable diagnostic biomarkers for LUSC.Lung squamous cell carcinoma (LUSC) accounts for a significant proportion of lung cancer and there have been few therapeutic alternatives for recurrent LUSC due to the lack of specific driver molecules. To investigate the prospective role of lncRNAs in the tumorigenesis and progression of LUSC, the aberrantly expressed lncRNAs were calculated based on The Cancer Genome Atlas RNA-seq data. Of 7589 lncRNAs with 504 LUSC cases, 884 lncRNAs were identified as being aberrantly expressed (|log2 fold change| >2 and adjusted P<0.05) by DESeq R. The top 10 lncRNAs with the highest diagnostic value were SFTA1P,LINC00968, LINC00961, LINC01572,RP1-78O14.1, FENDRR, LINC01314,LINC01272, GATA6-AS1, and MIR3945HG. In addition to the significant roles in the carcinogenesis of LUSC, several lncRNAs also played vital parts in the survival and progression of LUSC. SFTA1P, LINC01272, GATA6-AS1 and MIR3945HG were closely related to the survival time of LUSC. Furthermore, LINC01572 and LINC01314 could distinguish the LUSC at early stage from that at advanced stage. The prospective molecular assessment of key lncRNAs showed that a certain series of genes could be involved in the regulation network. Furthermore, the OncoPrint from cBioPortal indicated that 14% (69/501) LUSC cases with genetic alterations could be obtained, including amplification, deep deletion and mRNA upregulation. More interestingly, the cases with genetic alterations had a poorer survival as compared to those without alterations. Overall, the study propounds a potentiality for interpreting the pathogenesis and development of LUSC with lncRNAs, and provides a novel platform for searching for more capable diagnostic biomarkers for LUSC.

Comprehensive investigation of a novel differentially expressed lncRNA expression profile signature to assess the survival of patients with colorectal adenocarcinoma

Growing evidence has shown that long non-coding RNAs (lncRNAs) can serve as prospective markers for survival in patients with colorectal adenocarcinoma. However, most studies have explored a limited number of lncRNAs in a small number of cases. The objective of this study is to identify a panel of lncRNA signature that could evaluate the prognosis in colorectal adenocarcinoma based on the data from The Cancer Genome Atlas (TCGA). Altogether, 371 colon adenocarcinoma (COAD) patients with complete clinical data were included in our study as the test cohort. A total of 578 differentially expressed lncRNAs (DELs) were observed, among which 20 lncRNAs closely related to overall survival (OS) in COAD patients were identified using a Cox proportional regression model. A risk score formula was developed to assess the prognostic value of the lncRNA signature in COAD with four lncRNAs (LINC01555, RP11-610P16.1, RP11-108K3.1 and LINC01207), which were identified to possess the most remarkable correlation with OS in COAD patients. COAD patients with a high-risk score had poorer OS than those with a low-risk score. The multivariate Cox regression analyses confirmed that the four-lncRNA signature could function as an independent prognostic indicator for COAD patients, which was largely mirrored in the validating cohort with rectal adenocarcinoma (READ) containing 158 cases. In addition, the correlative genes of LINC01555 and LINC01207 were enriched in the cAMP signaling and mucin type O-Glycan biosynthesis pathways. With further validation in the future, our study indicates that the four-lncRNA signature could serve as an independent biomarker for survival of colorectal adenocarcinoma.

Identification of a RNA-Seq based prognostic signature with five lncRNAs for lung squamous cell carcinoma

Long non-coding RNAs (lncRNAs) expression profile signature for survival assessment in lung squamous cell carcinoma (LUSC) are largely inconsistent due to distinct detecting approaches and small sample size. Systematic and integrative investigation of RNA-Seq based data from The Cancer Genome Atlas (TCGA) herein was performed to determine candidate lncRNAs for prognosis evaluation of LUSC. A total of 60483 genes, including 7589 lncRNAs were assessed in a cohort including 478 LUSC cases with follow-up data. Firstly, 4225 differentially expressed lncRNAs were obtained via R packages. Next, univariate and multivariate Cox proportional hazards regression revealed that 41 lncRNAs were closely related to the survival of LUSC. Finally, lncRNA based prognosis index (PI) could predict overall survival of LUSC with high accuracy (AUC = 0.652, CI: 0.598, 0.705), PI = expCYP4F26P*βCYP4F26P+expRP11-108M12.3*βRP11-108M12.3+expRP11-38M8.1*βRP11-38M8.1+expRP11-54H7.4*βRP11-54H7.4+expZNF503-AS1*βZNF503-AS1. Furthermore, it was confirmed that the five-lncRNA signature could act as an independent prognostic indicator for LUSC (HR = 2.068, p < 0.001 with univariate analysis, HR = 1.928, p = 0.038 with multivariate). Besides, we constructed a weighted gene co-expression network analysis (WGCNA) of key lncRNA RP11-54H7.4 according to the p-value of related genes’ weight. This study provides a RNA-Seq based prognostic signature with five lncRNAs for further clinical application to LUSC patients.Long non-coding RNAs (lncRNAs) expression profile signature for survival assessment in lung squamous cell carcinoma (LUSC) are largely inconsistent due to distinct detecting approaches and small sample size. Systematic and integrative investigation of RNA-Seq based data from The Cancer Genome Atlas (TCGA) herein was performed to determine candidate lncRNAs for prognosis evaluation of LUSC. A total of 60483 genes, including 7589 lncRNAs were assessed in a cohort including 478 LUSC cases with follow-up data. Firstly, 4225 differentially expressed lncRNAs were obtained via R packages. Next, univariate and multivariate Cox proportional hazards regression revealed that 41 lncRNAs were closely related to the survival of LUSC. Finally, lncRNA based prognosis index (PI) could predict overall survival of LUSC with high accuracy (AUC = 0.652, CI: 0.598, 0.705), PI = expCYP4F26P*βCYP4F26P+expRP11-108M12.3*βRP11-108M12.3+expRP11-38M8.1*βRP11-38M8.1+expRP11-54H7.4*βRP11-54H7.4+expZNF503-AS1*βZNF503-AS1. Furthermore, it was confirmed that the five-lncRNA signature could act as an independent prognostic indicator for LUSC (HR = 2.068, p < 0.001 with univariate analysis, HR = 1.928, p = 0.038 with multivariate). Besides, we constructed a weighted gene co-expression network analysis (WGCNA) of key lncRNA RP11-54H7.4 according to the p-value of related genes weight. This study provides a RNA-Seq based prognostic signature with five lncRNAs for further clinical application to LUSC patients.

Clinicopathological role of miR-30a-5p in hepatocellular carcinoma tissues and prediction of its function with bioinformatics analysis

Background It has been reported that deregulation or dysfunction of microRNAs (miRNAs) plays an essential part in the hepatocarcinogenesis. However, the contribution and mechanism of microRNA-30a-5p (miR-30a-5p) in hepatocellular carcinoma (HCC) remains largely unknown. Therefore, our aim was to investigate the clinicopathological role of miR-30a-5p in HCC tissues and explore its potential pathways in this study. Methods The expression of miR-30a-5p was measured in 95 HCC and adjacent noncancer tissues by real-time reverse transcription quantitative polymerase chain reaction. The relationship between miR-30a-5p expression levels and clinicopathological parameters was also analyzed. Furthermore, the potential target genes of miR-30a-5p were collected via online prediction and literature searching. Gene ontology and pathway enrichment analyses were used to identify the possible function of miR-30a-5p in HCC. Results Compared with adjacent noncancer tissues (2.23±0.77), expression level of miR-30a-5p was significantly lower in HCC tissues (1.26±0.66, P<0.001). MiR-30a-5p expression was evidently correlated with tumor nodes, metastasis, tumor–node–metastasis stage, portal vein tumor embolus, vascular invasion, and status of tumor capsule (all P<0.05). A total of 878 genes were finally used for the biological informatics analyses. These prospective target genes were highly enriched in various key pathways, for instance, Ubiquitin-mediated proteolysis, Axon guidance, Neurotrophin signaling pathway, Amyotrophic lateral sclerosis, and ErbB signaling pathway. Conclusion In conclusion, this study clarifies that the downregulation of miRNA-30a-5p might play a vital part in the incidence and progression of HCC via targeting various prospective genes and pathways. Future validation is required to further explore the prospective molecular mechanism of miR-30a-5p in HCC.

Human papillomavirus as a potential risk factor for gastric cancer: a meta-analysis of 1,917 cases

Background Human papillomaviruses (HPVs) are causally associated with the tumorigenesis of several classes of cancers. However, the prevalence of HPV in gastric cancer (GC) has not yet been systematically reviewed. Hence, a meta-analysis was conducted to estimate the HPV prevalence in patients with GC, and its potential etiologic significance was assessed. Methods The pooled HPV prevalence and 95% confidence intervals (CIs) were estimated among all GC patients. Heterogeneity was described by using the I2 statistic. Sources of heterogeneity were explored by meta-regression and stratified analyses. The meta-influence was applied to evaluate the influence of a single study on the pooled estimates. Odds ratios (ORs) and 95% CIs were computed for case–control studies. For research providing clinicopathological parameters of age, sex, pathological, differentiated, and clinical stages, and HPV subtypes, the corresponding pooled ORs and 95% CIs were also calculated. Results Thirty studies were included in the current meta-analysis, involving 1,917 patients with GC and 576 controls. The pooled HPV prevalence was 28.0% (95% CI: 23.2%, 32.7%) among all the patients with GC, and the I2 was 96.9% (P<0.001). A pooled OR of 7.388 (95% CI: 3.876, 14.082) was achieved based on 15 case–control studies (I2=56.7%, P=0.004). Moreover, the HPV prevalence was significantly higher in patients from China than in those from non-Chinese regions (31% vs 9%, I2=95.0%, P<0.001). The pooled prevalence of HPV16 was 21% in GC tissues, and the pooled prevalence of HPV18 was 7% with an OR of 3.314 (95% CI =1.617, 6.792). HPV16 was 3 times more frequently detected than HPV18. Conclusion HPV could play a potential role in the pathogenesis of GC. A causal relationship can be confirmed only by detecting HPV in the cells of GC precursor lesions (gastric dysplasia or adenoma). In addition, this study might be beneficial for expounding the potential etiologic significance of molecular mechanism of gastric tumorigenesis and providing opinions regarding precautionary measures.

RNA-sequencing investigation identifies an effective risk score generated by three novel lncRNAs for the survival of papillary thyroid cancer patients

Scholars are striving to apply molecular biology involving long non-coding RNA (lncRNA) in the prognostication of papillary thyroid cancer (PTC). However, the clinical role of lncRNAs in the prognostic setting of PTC is still unclear. Herein, a comprehensive inquiry was performed to screen lncRNA expression profiling with 507 PTC patients from The Cancer Genome Atlas RNA-sequencing datasets. A total of 734 lncRNAs were detected to be aberrantly expressed, among which three novel lncRNAs including AC079630.2, CRNDE and CTD-2171N6.1 were markedly related to the progression and survival of PTC. Furthermore, the aberrant expression of these lncRNAs could be verified by other cohorts from gene expression omnibus (GEO) as detected by microarrays. Next, we established a three-lncRNA signature and divided the PTC patients into two subgroups of high- and low-risk. Interestingly, patients with high-risk tended to gain obviously poorer outcome. Most importantly, this three-lncRNA signature was an independent biomarker to predict the patient survival of PTC. The accurate molecular roles of these three lncRNAs remains unclarified and warrants further characterization, but our current data propose that they might play pivotal roles in PTC tumorigenesis and more importantly, these novel lncRNAs are closely related to patients’ survival. These discoveries will have far-reaching consequences with respect to molecular prediction of PTC.Scholars are striving to apply molecular biology involving long non-coding RNA (lncRNA) in the prognostication of papillary thyroid cancer (PTC). However, the clinical role of lncRNAs in the prognostic setting of PTC is still unclear. Herein, a comprehensive inquiry was performed to screen lncRNA expression profiling with 507 PTC patients from The Cancer Genome Atlas RNA-sequencing datasets. A total of 734 lncRNAs were detected to be aberrantly expressed, among which three novel lncRNAs including AC079630.2, CRNDE and CTD-2171N6.1 were markedly related to the progression and survival of PTC. Furthermore, the aberrant expression of these lncRNAs could be verified by other cohorts from gene expression omnibus (GEO) as detected by microarrays. Next, we established a three-lncRNA signature and divided the PTC patients into two subgroups of high- and low-risk. Interestingly, patients with high-risk tended to gain obviously poorer outcome. Most importantly, this three-lncRNA signature was an independent biomarker to predict the patient survival of PTC. The accurate molecular roles of these three lncRNAs remains unclarified and warrants further characterization, but our current data propose that they might play pivotal roles in PTC tumorigenesis and more importantly, these novel lncRNAs are closely related to patients survival. These discoveries will have far-reaching consequences with respect to molecular prediction of PTC.

Low Serum Sex Hormone-Binding Globulin Associated with Insulin Resistance in Men with Nonalcoholic Fatty Liver Disease

The presence of nonalcoholic fatty liver disease (NAFLD) is a strong risk predictor for type 2 diabetes (T2D). A reduction in sex hormone-binding globulin (SHBG) is associated with NAFLD. Low SHBG is also associated with insulin resistance (IR). However, very limited data are available for the association of SHBG and IR in patients with NAFLD. The study aims to clarify the association between SHBG and IR in patients with NAFLD. In this cross-sectional study, 334 men with NAFLD were recruited. SHBG, total testosterone, free testosterone, total cholesterol, triglyceride, insulin, and glucose concentrations were measured. Homeostatic model assessment (HOMA)-IR and HOMA-β were calculated. Spearmans correlations and multiple linear regressions were used to analyze the association between SHBG and IR. Men with moderate-severe NAFLD had higher waist circumference, BMI, total cholesterol, triglyceride, insulin, HOMA-IR, HOMA-β, and free testosterone, but lower SHBG than the mild NAFLD. The moderate-severe NAFLD group exhibited higher HOMA-IR (2.38±1.35 vs. 4.16±2.84, p<0.001) and lower SHBG (25.89±11.89 vs. 30.13±12.97u2009nmol/l, p<0.001) than the other group. SHBG value was negatively correlated with insulin, and HOMA-IR, but was not significantly correlated with glucose and testosterone. The multiple linear regression analysis showed that SHBG was significantly associated with insulin (β=-u20090.241, p<0.001), and HOMA-IR (β=-u20090.229, p<0.001), even adjusting for potential confounders. In conclusion, low serum SHBG is associated with IR in men with NAFLD.

Comprehensive investigation of aberrant microRNA profiling in bladder cancer tissues.

There has been accumulative evidence that microRNAs (miRNAs) play essential roles in the tumorigenesis and progression of bladder cancer. However, individual studies and small sample size caused discrepant outcomes. Thus, the current study focused on a comprehensive profiling of all differentially expressed miRNAs in a total of 519 bladder cancer tissue samples, based on miRNA microarray data. Altogether, 11 prioritized miRNAs stated by 21 published microarray datasets, including five down-regulated (miR-133a-3p, miR-1-3p, miR-99a-5p, miR-490-5p, and miR-133b) and six up-regulated candidate miRNAs (miR-182-5p, miR-935, miR-518e-3p, miR-573, miR-100-3p, and miR-3171) were analyzed with vote-counting strategy and a Robust Rank Aggregation method. Subsequently, miRNA in silico target prediction and potential pathway enrichment analysis were performed to investigate the prospective molecular mechanism of miRNAs in the tumorigenesis of bladder cancer. We found that most of the relative pathways of the aberrantly expressed miRNAs found in the current study were closely correlated with different biological processes, cellular components, molecular functions, cancer pathogeneses, and some cell signalings, such as Wnt signaling, insulin/IGF, PI3 kinase, and FGF signaling pathways. Hence, a comprehensive overview on the miRNA expression pattern in bladder cancer tissues was gained by the current study. These miRNAs might be involved in the tumorigenesis and deterioration of bladder cancer.

Down-regulation of miR-26a-5p in hepatocellular carcinoma: a qRT-PCR and bioinformatics study

BACKGROUNDnTo practically verify the clinical value of miR-26a-5p and thoroughly explore its target genes as well as its potential functions in hepatocellular carcinoma (HCC).nnnMETHODSnHCC and adjacent non-cancerous hepatic tissues of 95 HCC patients were collected for analysis using reverse transcription quantitative real-time PCR (qRT-PCR). For the bioinformatics analysis, we identified potential target genes for miR-26a-5p from differentially expressed genes (DEGs) in Gene Expression Omnibus (GEO) data sets and miRWalk predicted database. Gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway and protein-protein interaction (PPI) analyses were applied to analyze the prospective mechanisms of the predicted target genes.nnnRESULTSnMiR-26a-5p showed a significantly lower expression level in HCC tissues (1.56±1.07) than adjacent benign liver tissues (2.28±1.06, P<0.001). The area under the curve (AUC) of the receiver operating characteristic (ROC) was 0.665 (95% CI: 0.588-0.743, P<0.001). Significant correlations between miR-26a-5p expression and clinicopathological features such as gender (r=0.275, P<0.01), clinical TNM stage (r=-0.306, P<0.01), and metastasis (r=-0.321, P<0.01) were observed. To examine potential target genes, we obtained 175 genes for further function analysis, by attaining the intersection of 2062 up-regulated DEGs and 1390 online-predicted target genes. The GO and KEGG pathway annotation indicated focal adhesion, regulation of actin cytoskeleton and the PI3K-Akt signaling pathway as significant prospective mechanisms. The PPI network indicated that NRAS was the most essential hub gene in the whole network.nnnCONCLUSIONnDown-regulated miR-26a-5p was closely correlated with the status of metastasis and the progression of HCC. MiR-26a-5p might play protective roles by targeting diverse genes and pathways.