A Pilot Study on the Whole Exome Sequencing of Prostate Cancer in the Indian Phenotype Reveals Distinct Polymorphisms

Identification of Novel Stemness-based Subtypes and Construction of a Prognostic Risk Model for Patients with Lung Squamous Cell Carcinoma

BACKGROUND

Although cancer stem cells (CSCs) contribute to tumorigenesis, progression, and drug resistance, stemness-based classification and prognostic signatures of lung squamous cell carcinoma (LUSC) remain unclarified. This study attempted to identify stemness-based subtypes and develop a prognostic risk model for LUSC.

METHODS

Based on RNA-seq data from The Cancer Genome Atlas (TCGA), Gene-Expression Omnibus (GEO) and Progenitor Cell Biology Consortium (PCBC), mRNA expression-based stemness index (mRNAsi) was calculated by one-class logistic regression (OCLR) algorithm. A weighted gene coexpression network (WGCNA) was employed to identify stemness subtypes. Differences in mutation, clinical characteristics, immune cell infiltration, and antitumor therapy responses were determined. We constructed a prognostic risk model, followed by validations in GEO cohort, pan-cancer and immunotherapy datasets.

RESULTS

LUSC patients with subtype C2 had a better prognosis, manifested by higher mRNAsi, higher tumor protein 53 (TP53) and Titin (TTN) mutation frequencies, lower immune scores and decreased immune checkpoints. Patients with subtype C2 were more sensitive to Imatinib, Pyrimethamine, and Paclitaxel therapy, whereas those with subtype C1 were more sensitive to Sunitinib, Saracatinib, and Dasatinib. Moreover, we constructed stemness-based signatures using seven genes (BMI1, CCDC51, CTNS, EIF1AX, FAM43A, THBD, and TRIM68) and found high-risk patients had a poorer prognosis in the TCGA cohort. Similar results were found in the GEO cohort. We verified the good performance of risk scores in prognosis prediction and therapy responses.

CONCLUSION

The stemness-based subtypes shed novel insights into the potential roles of LUSC-stemness in tumor heterogeneity, and our prognostic signatures offer a promising tool for prognosis prediction and guide therapeutic decisions in LUSC.

Whole-exome sequencing of Indian prostate cancer reveals a novel therapeutic target: POLQ

PURPOSE

Prostate cancer is the second most common cancer diagnosed worldwide and the third most common cancer among men in India. This study's objective was to characterise the mutational landscape of Indian prostate cancer using whole-exome sequencing to identify population-specific polymorphisms.

METHODS

Whole-exome sequencing was performed of 58 treatment-naive primary prostate tumors of Indian origin. Multiple computational and statistical analyses were used to profile the known common mutations, other deleterious mutations, driver genes, prognostic biomarkers, and gene signatures unique to each clinical parameter. Cox analysis was performed to validate survival-associated genes. McNemar test identified genes significant to recurrence and receiver-operating characteristic (ROC) analysis was conducted to determine its accuracy. OncodriveCLUSTL algorithm was used to deduce driver genes. The druggable target identified was modeled with its known inhibitor using Autodock.

RESULTS

TP53 was the most commonly mutated gene in our cohort. Three novel deleterious variants unique to the Indian prostate cancer subtype were identified: POLQ, FTHL17, and OR8G1. COX regression analysis identified ACSM5, a mitochondrial gene responsible for survival. CYLC1 gene, which encodes for sperm head cytoskeletal protein, was identified as an unfavorable prognostic biomarker indicative of recurrence. The novel POLQ mutant, also identified as a driver gene, was evaluated as the druggable target in this study. POLQ, a DNA repair enzyme implicated in various cancer types, is overexpressed and is associated with a poor prognosis. The mutant POLQ was subjected to structural analysis and modeled with its known inhibitor novobiocin resulting in decreased binding efficiency necessitating the development of a better drug.

CONCLUSION

In this pilot study, the molecular profiling using multiple computational and statistical analyses revealed distinct polymorphisms in the Indian prostate cancer cohort. The mutational signatures identified provide a valuable resource for prognostic stratification and targeted treatment strategies for Indian prostate cancer patients. The DNA repair enzyme, POLQ, was identified as the druggable target in this study.

Towards Understanding the Key Signature Pathways Associated from Differentially Expressed Gene Analysis in an Indian Prostate Cancer Cohort

Prostate cancer (PCa) is one of the most prevalent cancers among men in India. Although studies on PCa have dealt with genetics, genomics, and the environmental influence in the causality of PCa, not many studies employing the Next Generation Sequencing (NGS) approaches of PCa have been carried out. In our previous study, we identified some causal genes and mutations specific to Indian PCa using Whole Exome Sequencing (WES). In the recent past, with the help of different cancer consortiums such as The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC), along with differentially expressed genes (DEGs), many cancer-associated novel non-coding RNAs have been identified as biomarkers. In this work, we attempt to identify differentially expressed genes (DEGs) including long non-coding RNAs (lncRNAs) associated with signature pathways from an Indian PCa cohort using the RNA-sequencing (RNA-seq) approach. From a cohort of 60, we screened six patients who underwent prostatectomy; we performed whole transcriptome shotgun sequencing (WTSS)/RNA-sequencing to decipher the DEGs. We further normalized the read counts using fragments per kilobase of transcript per million mapped reads (FPKM) and analyzed the DEGs using a cohort of downstream regulatory tools, viz., GeneMANIA, Stringdb, Cytoscape-Cytohubba, and cbioportal, to map the inherent signatures associated with PCa. By comparing the RNA-seq data obtained from the pairs of normal and PCa tissue samples using our benchmarked in-house cuffdiff pipeline, we observed some important genes specific to PCa, such as STEAP2, APP, PMEPA1, PABPC1, NFE2L2, and HN1L, and some other important genes known to be involved in different cancer pathways, such as COL6A1, DOK5, STX6, BCAS1, BACE1, BACE2, LMOD1, SNX9, CTNND1, etc. We also identified a few novel lncRNAs such as LINC01440, SOX2OT, ENSG00000232855, ENSG00000287903, and ENST00000647843.1 that need to be characterized further. In comparison with publicly available datasets, we have identified characteristic DEGs and novel lncRNAs implicated in signature PCa pathways in an Indian PCa cohort which perhaps have not been reported. This has set a precedent for us to validate candidates further experimentally, and we firmly believe this will pave a way toward the discovery of biomarkers and the development of novel therapies.

Characterization of the Impact of Density Gradient Centrifugation on the Profile of the Pig Sperm Transcriptome by RNA-Seq

RNA-Seq data from human semen suggests that the study of the sperm transcriptome requires the previous elimination from the ejaculates of somatic cells carrying a larger load of RNA. Semen purification is also carried to study the sperm transcriptome in other species including swine and it is often done by density gradient centrifugation to obtain viable spermatozoa from fresh ejaculates or artificial insemination doses, thereby limiting the throughput and remoteness of the samples that can be processed in one study. The aim of this work was to evaluate the impact of purification with density gradient centrifugation by BoviPure^TM on porcine sperm. Four boar ejaculates were purified with BoviPure^TM and their transcriptome sequenced by RNA-Seq was compared with the RNA-Seq profiles of their paired non-purified sample. Seven thousand five hundred and nineteen protein coding genes were identified. Correlation, cluster, and principal component analysis indicated high—although not complete—similarity between the purified and the paired non-purified ejaculates. 372 genes displayed differentially abundant RNA levels between treatments. Most of these genes had lower abundances after purification and were mostly related to translation, transcription and metabolic processes. We detected a significant change in the proportion of genes of epididymal origin within the differentially abundant genes (1.3%) when compared with the catalog of unaltered genes (0.2%). In contrast, the proportion of testis-specific genes was higher in the group of unaltered genes (4%) when compared to the list of differentially abundant genes (0%). No proportion differences were identified for prostate, white blood, lymph node, tonsil, duodenum, skeletal muscle, liver, and mammary gland. Altogether, these results suggest that the purification impacts on the RNA levels of a small number of genes which are most likely caused by the removal of epididymal epithelial cells but also premature germinal cells, immature or abnormal spermatozoa or seminal exosomes with a distinct load of RNAs.