Impact of Thawing on Reference Gene Expression Stability in Renal Cell Carcinoma Samples:

Integration of single-cell RNA-seq and bulk RNA-seq to construct liver hepatocellular carcinoma stem cell signatures to explore their impact on patient prognosis and treatment

BACKGROUND

Liver hepatocellular carcinoma (LIHC) is a prevalent form of primary liver cancer. Research has demonstrated the contribution of tumor stem cells in facilitating tumor recurrence, metastasis, and treatment resistance. Despite this, there remains a lack of established cancer stem cells (CSCs)-associated genes signatures for effectively predicting the prognosis and guiding the treatment strategies for patients diagnosed with LIHC.

METHODS

The single-cell RNA sequencing (scRNA-seq) and bulk RNA transcriptome data were obtained based on public datasets and computerized firstly using CytoTRACE package and One Class Linear Regression (OCLR) algorithm to evaluate stemness level, respectively. Then, we explored the association of stemness indicators (CytoTRACE score and stemness index, mRNAsi) with survival outcomes and clinical characteristics by combining clinical information and survival analyses. Subsequently, weighted co-expression network analysis (WGCNA) and Cox were applied to assess mRNAsi-related genes in bulk LIHC data and construct a prognostic model for LIHC patients. Single-sample gene-set enrichment analysis (ssGSEA), Cell-type Identification By Estimating Relative Subsets Of RNA Transcripts (CIBERSORT) and Tumor Immune Estimation Resource (TIMER) analysis were employed for immune infiltration assessment. Finally, the potential immunotherapeutic response was predicted by the Tumor Immune Dysfunction and Exclusion (TIDE), and the tumor mutation burden (TMB). Additionally, pRRophetic package was applied to evaluate the sensitivity of high and low-risk groups to common chemotherapeutic drugs.

RESULTS

A total of four genes (including STIP1, H2AFZ, BRIX1, and TUBB) associated with stemness score (CytoTRACE score and mRNAsi) were identified and constructed a risk model that could predict prognosis in LIHC patients. It was observed that high stemness cells occurred predominantly in the late stages of LIHC and that poor overall survival in LIHC patients was also associated with high mRNAsi scores. In addition, pathway analysis confirmed the biological uniqueness of the two risk groups. Personalized treatment predictions suggest that patients with a low risk benefited more from immunotherapy, while those with a high risk group may be conducive to chemotherapeutic drugs.

CONCLUSION

The current study developed a novel prognostic risk signature with genes related to CSCs, which provides novel ideas for the diagnosis, prognosis and treatment of LIHC.

Determination of reference genes as a quantitative standard for gene expression analysis in mouse mesangial cells stimulated with TGF-β

Reverse transcription-quantitative polymerase chain reaction (RT-PCR) is the gold standard technique for gene expression analysis, but the choice of quantitative reference genes (housekeeping genes, HKG) remains challenging. Identify the best HKG is essential for estimating the expression level of target genes. Therefore, the aim of this study was to determine the best HKG for an in vitro model with mouse mesangial cells (MMCs) stimulated with 5 ng/mL of TGF-β. Five candidates HKG were selected: Actb, Hprt, Gapdh, 18S and Ppia. After quantitative expression, the best combination of these genes was analyzed in silico using six software programs. To validate the results, the best genes were used to normalize the expression levels of fibronectin, vimentin and α-SMA. In silico analysis revealed that Ppia, Gapdh and 18S were the most stable genes between the groups. GenEX software and Spearman's correlation determined Ppia and Gapdh as the best HKG pair, and validation of the HKG by normalizing fibronectin, vimentin and α-SMA were consistent with results from the literature. Our results established the combination of Ppia and Gapdh as the best HKG pair for gene expression analysis by RT-PCR in this in vitro model using MMCs treated with TGF-β.

Bioinformatics Identification of TUBB as Potential Prognostic Biomarker for Worse Prognosis in ERα-Positive and Better Prognosis in ERα-Negative Breast Cancer

Tubulin β class I gene (TUBB) is highly expressed in various cancers and plays several roles in carcinogenesis. However, the prognostic value of TUBB in breast cancer remains to be investigated. GEPIA and Breast Cancer Gene-Expression Miner were used to explore TUBB expression in breast cancer patients. Kaplan–Meier Plotter was used to assess the relationship between TUBB expression and several prognostic indicators including overall, distant metastasis-free, and relapse-free survival in ERα-positive and ERα-negative breast cancer. The genes that correlate with TUBB in ERα-positive and ERα-negative breast cancer were explored and the pathways were investigated using GSCA. The correlation between TUBB and several gene markers of immune cells was explored using GEPIA. ERα-positive breast cancer patients with increased TUBB showed worse prognosis, possibly through the activation of the TSC/mTOR pathway, whereas ERα-negative breast cancer patients with increased TUBB mRNA showed better prognosis. Significant positive correlations were observed between TUBB and gene markers of immune cells in ERα-positive breast cancer patients, whereas significant negative correlations were observed in ERα-negative breast cancer patients. The analysis revealed that TUBB might be considered as a predictive biomarker for worse prognosis in ERα-positive and better prognosis in ERα-negative breast cancer.

Ppia is the most stable housekeeping gene for qRT-PCR normalization in kidneys of three Pkd1-deficient mouse models

Autosomal Dominant Polycystic Kidney Disease (ADPKD) is the most common inherited renal disorder, characterized by renal cyst development leading to end-stage renal disease. Although the appropriate choice of suitable reference is critical for quantitative RNA analysis, no comparison of frequently used “housekeeping” genes is available. Here, we determined the validity of 7 candidate housekeeping genes (Actb, Actg1, B2m, Gapdh, Hprt, Pgam1 and Ppia) in kidney tissues from mouse models orthologous to ADPKD, including a cystic mice (CY) 10–12 weeks old (Pkd1^flox/flox:Nestin^cre/Pkd1^flox/−:Nestin^cre, n = 10) and non-cystic (NC) controls (Pkd1^flox/^flox/Pkd1^flox/^-, n = 10), Pkd1-haploinsufficient (HT) mice (Pkd1^+/−, n = 6) and wild-type (WT) controls (Pkd1^+/+, n = 6) and a severely cystic (SC) mice 15 days old (Pkd1^V/V, n = 7) and their controls (CO, n = 5). Gene expression data were analyzed using six distinct statistical softwares. The estimation of the ideal number of genes suggested the use of Ppia alone as sufficient, although not ideal, to analyze groups altogether. Actb, Hprt and Ppia expression profiles were correlated in all samples. Ppia was identified as the most stable housekeeping gene, while Gapdh was the least stable for all kidney samples. Stat3 expression level was consistent with upregulation in SC compared to CO when normalized by Ppia expression. In conclusion, present findings identified Ppia as the best housekeeping gene for CY + NC and SC + CO groups, while Hprt was the best for the HT + WT group.

Seminal exosomal miR-210-3p as a potential marker of Sertoli cell damage in Varicocele

BACKGROUND

Varicocoele-associated stressors, such as hypoxia and heat, can damage cell function and viability, and some exosomal biomarkers released from impaired cells may reflect the cell status in testis.

OBJECTIVES

To find if seminal exosomal microRNAs can reflect the Sertoli cell function in varicocoele.

MATERIALS AND METHODS

Experimental left varicocoele rat model was established (n = 24), and patients with different grades of varicocoele (n = 104) were enrolled. Primary rat Sertoli cells were isolated with enzymatic hydrolysis. Exosomes were isolated from primary rat Sertoli cells, rat epididymis tissue, and human seminal plasma with polymer-based precipitation method. Exosomal microRNAs were quantified with qPCR. Inhibin-B was detected with enzyme immunoassay. The correlation analysis between microRNA and inhibin-B was evaluated with Spearman's correlation.

RESULTS

We screened 12 previously reported hypoxia-responsive microRNAs in the primary rat Sertoli cells and found that 4 exosomal microRNAs increased significantly in response to in vitro hypoxia treatment (P < .05). Of the 4 microRNAs, only miR-210-3p was upregulated in the rats with experimental varicocoele (P < .01). In the patients with varicocoele, we found that seminal exosomal miR-210-3p significantly increased in patients with grade II and III varicocoele (P < .01), and miR-210-3p negatively correlated with sperm count (P < .01) and seminal inhibin-B expression (r = -0.39, P < .01). For the 30 patients with microsurgical varicocelectomy, the operation notably decreased miR-210-3p (P < .01).

DISCUSSION AND CONCLUSION

Seminal exosomal miR-210-3p may be a novel, sensitive, and non-invasive biomarker of Sertoli cell damage in varicocoele.

A frequent somatic mutation in the 3'UTR of GAPDH facilitates the development of ovarian cancer by creating a miR‑125b binding site

Ovarian cancer (OVCA) is one of the most common types of cancer in women worldwide. Recent studies have focused on the presence and effect of somatic mutations in patients with OVCA; however, studies on the roles of mutations located in the untranslated regions (UTR) of genes in OVCA remain limited. In the present study, a frequent somatic mutation in the glyceraldehyde 3-phosphate dehydrogenase (GADPH) 3′UTR was identified using transcriptome sequencing of 120 pairs of OVCA tissue samples. The mutant GAPDH 3′UTR promoted tumor growth and cell motility. Furthermore, the mutation in the GAPDH 3′UTR significantly downregulated the levels of mature miR-125b by creating a new miR-125b binding site. Finally, STAT3 levels were increased in SKOV3 cells stably expressing the mutant GADPH 3′UTR, which is a critical target gene of miR-125b. In conclusion, the present study demonstrated that the mutation located in GAPDH 3′UTR promoted OVCA growth and development by sponging miR-125b and thereby affecting STAT3 expression levels.

Validation of Reference Genes for Normalization of Relative qRT-PCR Studies in Papillary Thyroid Carcinoma

Quantitative reverse transcription polymerase chain reaction (qRT-PCR) in thyroid tumors require accurate data normalization, however, there are no sufficient studies addressing the suitable reference genes for gene expression analysis in malignant and normal thyroid tissue specimens. The purpose of this study was to identify valid internal control genes for normalization of relative qRT-PCR studies in human papillary thyroid carcinoma tissue samples. The expression characteristics of 12 candidate reference genes (GAPDH, ACTB, HPRT1, TBP, B2M, PPIA, 18SrRNA, HMBS, GUSB, PGK1, RPLP0, and PGM1) were assessed by qRT-PCR in 45 thyroid tissue samples (15 papillary thyroid carcinoma, 15 paired normal tissues and 15 multinodular goiters). These twelve candidate reference genes were selected by a systematic literature search. GeNorm, NormFinder, and BestKeeper statistical algorithms were applied to determine the most stable reference genes. The three algorithms were in agreement in identifying GUSB and HPRT1 as the most stably expressed genes in all thyroid tumors investigated. According to the NormFinder software, the pair of genes including ‘GUSB and HPRT1’ or ‘GUSB and HMBS’ or ‘GUSB and PGM1’ were the best combinations for selection of pair reference genes. The optimal number of genes required for reliable normalization of qPCR data in thyroid tissues would be three according to calculations made by GeNorm algorithm. These results suggest that GUSB and HPRT1 are promising reference genes for normalization of relative qRT-PCR studies in papillary thyroid carcinoma.

Identification of suitable reference genes for real-time quantitative PCR analysis of hydrogen peroxide-treated human umbilical vein endothelial cells

Background

Oxidative stress can induce cell injury in vascular endothelial cells, which is the initial event in the development of atherosclerosis. Although quantitative real-time polymerase chain reaction (qRT-PCR) has been widely used in gene expression studies in oxidative stress injuries, using carefully validated reference genes has not received sufficient attention in related studies. The objective of this study, therefore, was to select a set of stably expressed reference genes for use in qRT-PCR normalization in oxidative stress injuries in human umbilical vein endothelial cells (HUVECs) induced by hydrogen peroxide (H₂O₂).

Results

Using geNorm analysis, we found that five stably expressed reference genes were sufficient for normalization in qRT-PCR analysis in HUVECs treated with H₂O₂. Genes with the most stable expression according to geNorm were U6, TFRC, RPLP0, GAPDH, and ACTB, and according to NormFinder were ALAS1, TFRC, U6, GAPDH, and ACTB.

Conclusion

Taken together, our study demonstrated that the expression stability of reference genes may differ according to the statistical program used. U6, TFRC, RPLP0, GAPDH, and ACTB was the optimal set of reference genes for studies on gene expression performed by qRT-PCR assays in HUVECs under oxidative stress study.

The focus on sample quality: Influence of colon tissue collection on reliability of qPCR data

Successful molecular analyses of human solid tissues require intact biological material with well-preserved nucleic acids, proteins, and other cell structures. Pre-analytical handling, comprising of the collection of material at the operating theatre, is among the first critical steps that influence sample quality. The aim of this study was to compare the experimental outcomes obtained from samples collected and stored by the conventional means of snap freezing and by PAXgene Tissue System (Qiagen). These approaches were evaluated by measuring rRNA and mRNA integrity of the samples (RNA Quality Indicator and Differential Amplification Method) and by gene expression profiling. The collection procedures of the biological material were implemented in two hospitals during colon cancer surgery in order to identify the impact of the collection method on the experimental outcome. Our study shows that the pre-analytical sample handling has a significant effect on the quality of RNA and on the variability of qPCR data. PAXgene collection mode proved to be more easily implemented in the operating room and moreover the quality of RNA obtained from human colon tissues by this method is superior to the one obtained by snap freezing.

Testosterone regulates the autophagic clearance of androgen binding protein in rat Sertoli cells

Dysregulation of androgen-binding protein (ABP) is associated with a number of endocrine and andrology diseases. However, the ABP metabolism in Sertoli cells is largely unknown. We report that autophagy degrades ABP in rat Sertoli cells, and the autophagic clearance of ABP is regulated by testosterone, which prolongs the ABP biological half-life by inhibiting autophagy. Further studies identified that the autophagic clearance of ABP might be selectively regulated by testosterone, independent of stress (hypoxia)-induced autophagic degradation. These data demonstrate that testosterone up-regulates ABP expression at least partially by suppressing the autophagic degradation. We report a novel finding with respect to the mechanisms by which ABP is cleared, and by which the process is regulated in Sertoli cells.

Identification of potential serum proteomic biomarkers for clear cell renal cell carcinoma

Objective

To investigate discriminating protein patterns and serum biomarkers between clear cell renal cell carcinoma (ccRCC) patients and healthy controls, as well as between paired pre- and post-operative ccRCC patients.

Methods

We used magnetic bead-based separation followed by matrix-assisted laser desorption ionization (MALDI) time-of-flight (TOF) mass spectrometry (MS) to identify patients with ccRCC. A total of 162 serum samples were analyzed in this study, among which there were 58 serum samples from ccRCC patients, 40 from additional paired pre- and post-operative ccRCC patients (n = 20), and 64 from healthy volunteers as healthy controls. ClinProTools software identified several distinct markers between ccRCC patients and healthy controls, as well as between pre- and post-operative patients.

Results

Patients with ccRCC could be identified with a mean sensitivity of 88.38% and a mean specificity of 91.67%. Of 67 m/z peaks that differed among the ccRCC, healthy controls, pre- and post-operative ccRCC patients, 24 were significantly different (P<0.05). Three candidate peaks, which were upregulated in ccRCC group and showed a tendency to return to healthy control values after surgery, were identified as peptide regions of RNA-binding protein 6 (RBP6), tubulin beta chain (TUBB), and zinc finger protein 3 (ZFP3) with the m/z values of 1466.98, 1618.22, and 5905.23, respectively.

Conclusion

MB-MALDI-TOF-MS method could generate serum peptidome profiles of ccRCC, and provide a new approach to identify potential biomarkers for diagnosis as well as prognosis of this malignancy.

Biphasic regulation of autophagy by miR-96 in prostate cancer cells under hypoxia

Autophagy favors cell survival under hypoxia, and increasing evidence revealed that microRNAs regulate autophagy. We report here hypoxia increased the expression of miR-96 in prostate cancer cells, and miR-96 stimulated autophagy by suppressing MTOR. We found that inhibition of miR-96 abolished hypoxia-induced autophagy. Paradoxically, ectopic over-expression of miR-96 to a certain threshold, also abolished the hypoxia-induced autophagy. Further studies have shown that high levels of miR-96 inhibited autophagy through suppressing ATG7, a key autophagy-associated gene. Importantly, the miR-96 expression level threshold was determined, and the effects of miR-96 on autophagy on either side of the threshold were opposite. These data demonstrate hypoxia-induced autophagy is at least partially regulated by miR-96; miR-96 can promote or inhibit autophagy by principally inhibiting MTOR or ATG7 depending on the expression levels of miR-96. Our observation might reveal a novel regulatory mode of autophagy by microRNAs under hypoxia.

Autophagy inhibits chemotherapy-induced apoptosis through downregulating Bad and Bim in hepatocellular carcinoma cells

The tumor microenvironment, including ischemia, has been increasingly recognized as a critical factor in the process of tumor development. Hypoxia and nutrient deficiency resulting from ischemia widely exist in solid tumors. Recent studies have shown that hypoxia and nutrient deficiency contribute to chemoresistance by inducing autophagy, but the underlying mechanism remains unknown. This study aimed to explore the role of autophagy induced by low glucose and hypoxia (LH) in the chemoresistance of hepatocellular carcinoma cells. Our results demonstrated that LH induced autophagy and downregulated Bad and Bim in hepatocellular carcinoma cells. The inhibition of autophagy reversed the reduction of these pro-apoptotic factors during the LH treatment. Furthermore, Bad and Bim were also significantly downregulated by autophagy during the process that LH promoted the chemoresistance of hepatocellular carcinoma cells. In addition, RNAi or the overexpression of Bad and Bim can significantly reduce or increase chemotherapy-induced cell death, respectively. Taken together, these data indicate that the downregulation of Bad and Bim plays a significant role in the autophagy-induced chemoresistance of hepatocellular carcinoma cells.

Renal tissue thawed for 30 minutes is still suitable for gene expression analysis

Some biosamples obtained from biobanks may go through thawing before processing. We aim to evaluate the effects of thawing at room temperature for different time periods on gene expression analysis. A time course study with four time points was conducted to investigate the expression profiling on 10 thawed normal mice renal tissue samples through Affymetrix GeneChip mouse gene 2.0 st array. Microarray results were validated by quantitative real time polymerase chain reactions (qPCR) on 6 candidate reference genes and 11 target genes. Additionally, we used geNorm plus and NormFinder to identify the most stably expressed reference genes over time. The results showed RNA degraded more after longer incubation at room temperature. However, microarray results showed only 240 genes (0.91%) altered significantly in response to thawing at room temperature. The signal of majority altered probe sets decreased with thawing time, and the crossing point (Cp) values of all candidate reference genes correlated positively with the thawing time (p<0.05). The combination of B2M, ACTB and PPIA was identified as the best choice for qPCR normalization. We found most target genes were stable by using this normalization method. However, serious gene quantification errors were resulted from improper reference genes. In conclusion, thirty minutes of thawing at room temperature has a limited impact on microarray and qPCR analysis, gene expression variations due to RNA degradation in early period after thawing can be largely reduced by proper normalization.

Identification of optimal reference genes for RT-qPCR in the rat hypothalamus and intestine for the study of obesity

BACKGROUND

Obesity has a complicated metabolic pathology, and defining the underlying mechanisms of obesity requires integrative studies with molecular end points. Real-time quantitative PCR (RT-qPCR) is a powerful tool that has been widely utilized. However, the importance of using carefully validated reference genes in RT-qPCR seems to have been overlooked in obesity-related research. The objective of this study was to select a set of reference genes with stable expressions to be used for RT-qPCR normalization in rats under fasted vs re-fed and chow vs high-fat diet (HFD) conditions.

DESIGN

Male long-Evans rats were treated under four conditions: chow/fasted, chow/re-fed, HFD/fasted and HFD/re-fed. Expression stabilities of 13 candidate reference genes were evaluated in the rat hypothalamus, duodenum, jejunum and ileum using the ReFinder software program. The optimal number of reference genes needed for RT-qPCR analyses was determined using geNorm.

RESULTS

Using geNorm analysis, we found that it was sufficient to use the two most stably expressed genes as references in RT-qPCR analyses for each tissue under specific experimental conditions. B2M and RPLP0 in the hypothalamus, RPS18 and HMBS in the duodenum, RPLP2 and RPLP0 in the jejunum and RPS18 and YWHAZ in the ileum were the most suitable pairs for a normalization study when the four aforementioned experimental conditions were considered.

CONCLUSIONS

Our study demonstrates that gene expression levels of reference genes commonly used in obesity-related studies, such as ACTB or RPS18, are altered by changes in acute or chronic energy status. These findings underline the importance of using reference genes that are stable in expression across experimental conditions when studying the rat hypothalamus and intestine, because these tissues have an integral role in the regulation of energy homeostasis. It is our hope that this study will raise awareness among obesity researchers on the essential need for reference gene validation in gene expression studies.