Public transcriptome database-based selection and validation of reliable reference genes for breast cancer research

Screening and validating the optimal panel of housekeeping genes for 4T1 breast carcinoma and metastasis studies in mice

The 4T1 model is extensively employed in murine studies to elucidate the mechanisms underlying the carcinogenesis of triple-negative breast cancer. Molecular biology serves as a cornerstone in these investigations. However, accurate gene expression analyses necessitate data normalization employing housekeeping genes (HKGs) to avert spurious results. Here, we initially delve into the characteristics of the tumor evolution induced by 4T1 in mice, underscoring the imperative for additional tools for tumor monitoring and assessment methods for tracking the animals, thereby facilitating prospective studies employing this methodology. Subsequently, leveraging various software platforms, we assessed ten distinct HKGs (GAPDH, 18 S, ACTB, HPRT1, B2M, GUSB, PGK1, CCSER2, SYMPK, ANKRD17) not hitherto evaluated in the 4T1 breast cancer model, across tumors and diverse tissues afflicted by metastasis. Our principal findings underscore GAPDH as the optimal HKG for gene expression analyses in tumors, while HPRT1 emerged as the most stable in the liver and CCSER2 in the lung. These genes demonstrated consistent expression and minimal variation among experimental groups. Furthermore, employing these HKGs for normalization, we assessed TNF-α and VEGF expression in tissues and discerned significant disparities among groups. We posit that this constitutes the inaugural delineation of an ideal HKG for experiments utilizing the 4T1 model, particularly in vivo settings.

Validation of Endogenous Control Genes by Real-Time Quantitative Reverse Transcriptase Polymerase Chain Reaction for Acute Leukemia Gene Expression Studies

Reference genes are used as internal reaction controls for gene expression analysis, and for this reason, they are considered reliable and must meet several important criteria. In view of the absence of studies regarding the best reference gene for the analysis of acute leukemia patients, a panel of genes commonly used as endogenous controls was selected from the literature for stability analysis: Glyceraldehyde-3-phosphate dehydrogenase (GAPDH), Abelson murine leukemia viral oncogene human homolog 1 (ABL), Hypoxanthine phosphoribosyl-transferase 1 (HPRT1), Ribosomal protein lateral stalk subunit P0 (RPLP0), β-actin (ACTB) and TATA box binding protein (TBP). The stability of candidate reference genes was analyzed according to three statistical methods of assessment, namely, NormFinder, GeNorm and R software (version 4.0.3). From this study's analysis, it was possible to identify that the endogenous set composed of ACTB, ABL, TBP and RPLP0 demonstrated good performances and stable expressions between the analyzed groups. In addition to that, the GAPDH and HPRT genes could not be classified as good reference genes, considering that they presented a high standard deviation and great variability between groups, indicating low stability. Given these findings, this study suggests the main endogenous gene set for use as a control/reference for the gene expression in peripheral blood and bone marrow samples from patients with acute leukemias is composed of the ACTB, ABL, TBP and RPLP0 genes. Researchers may choose two to three of these housekeeping genes to perform data normalization.

Selection of reference genes in liproxstatin-1-treated K562 Leukemia cells via RT-qPCR and RNA sequencing

BACKGROUND

Reverse transcription quantitative polymerase chain reaction (RT-qPCR) can accurately detect relative gene expression levels in biological samples. However, widely used reference genes exhibit unstable expression under certain conditions.

METHODS AND RESULTS

Here, we compared the expression stability of eight reference genes (RPLP0, RPS18, RPL13, EEF1A1, β-actin, GAPDH, HPRT1, and TUBB) commonly used in liproxstatin-1 (Lip-1)-treated K562 cells using RNA-sequencing and RT-qPCR. The expression of EEF1A1, ACTB, GAPDH, HPRT1, and TUBB was considerably lower in cells treated with 20 μM Lip-1 than in the control, and GAPDH also showed significant downregulation in the 10 μM Lip-1 group. Meanwhile, when we used geNorm, NormFinder, and BestKeeper to compare expression stability, we found that GAPDH and HPRT1 were the most unstable reference genes among all those tested. Stability analysis yielded very similar results when geNorm or BestKeeper was used but not when NormFinder was used. Specifically, geNorm and BestKeeper identified RPL13 and RPLP0 as the most stable genes under 20 μM Lip-1 treatment, whereas RPL13, EEF1A1, and TUBB were the most stable under 10 μM Lip-1 treatment. TUBB and EEF1A1 were the most stable genes in both treatment groups according to the results obtained using NormFinder. An assumed most stable gene was incorporated into each software to validate the accuracy. The results suggest that NormFinder is not an appropriate algorithm for this study.

CONCLUSIONS

Stable reference genes were recognized using geNorm and BestKeeper but not NormFinder. Overall, RPL13 and RPLP0 were the most stable reference genes under 20 μM Lip-1 treatment, whereas RPL13, EEF1A1, and TUBB were the most stable genes under 10 μM Lip-1 treatment.

Identification of Prognostic Biomarkers for Breast Cancer Metastasis Using Penalized Additive Hazards Regression Model

Background:

Breast cancer (BC) has been reported as one of the most common cancers diagnosed in females throughout the world. Survival rate of BC patients is affected by metastasis. So, exploring its underlying mechanisms and identifying related biomarkers to monitor BC relapse/recurrence using new statistical methods is essential. This study investigated the high-dimensional gene-expression profiles of BC patients using penalized additive hazards regression models.

Methods:

A publicly available dataset related to the time to metastasis in BC patients (GSE2034) was used. There was information of 22 283 genes expression profiles related to 286 BC patients. Penalized additive hazards regression models with different penalties, including LASSO, SCAD, SICA, MCP and Elastic net were used to identify metastasis related genes.

Results:

Five regression models with penalties were applied in the additive hazards model and jointly found 9 genes including SNU13, CLINT1, MAPK9, ABCC5, NKX3-1, NCOR2, COL2A1, and ZNF219. According the median of the prognostic index calculated using the regression coefficients of the penalized additive hazards model, the patients were labeled as high/low risk groups. A significant difference was detected in the survival curves of the identified groups. The selected genes were examined using validation data and were significantly associated with the hazard of metastasis.

Conclusion:

This study showed that MAPK9, NKX3-1, NCOR1, ABCC5, and CD44 are the potential recurrence and metastatic predictors in breast cancer and can be taken into account as candidates for further research in tumorigenesis, invasion, metastasis, and epithelial-mesenchymal transition of breast cancer.