TPX2/Aurora kinase A signaling as a potential therapeutic target in genomically unstable cancer cells

Aurkin-A, a TPX2-Aurora A small molecule inhibitor disrupts Alisertib-induced polyploidy in aggressive diffuse large B cell lymphoma

Chemotherapy induced polyploidy is a mechanism of inherited drug resistance resulting in an aggressive disease course in cancer patients. Alisertib, an Aurora Kinase A (AK-A) ATP site inhibitor, induces cell cycle disruption resulting in polyaneuploidy in Diffuse Large B Cell Lymphoma (DLBCL). Propidium iodide flow cytometry was utilized to quantify alisertib induced polyploidy in U2932 and VAL cell lines. In U2932 cells, 1µM alisertib generated 8n+ polyploidy in 48% of the total cell population after 5 days of treatment. Combination of Aurkin A an AK-A/TPX2 site inhibitor, plus alisertib disrupted alisertib induced polyploidy in a dose-dependent manner with associated increased apoptosis. We generated a stable FUCCI U2932 cell line expressing Geminin-clover (S/G2/M) and cdt1-mKO (G1), to monitor cell cycle progression. Using this system, we identified alisertib induces polyploidy through endomitosis, which was eliminated with Aurkin A treatment. In a VAL mouse xenograft model, we show polyploidy generation in alisertib treated mice versus vehicle control or Aurkin A. Aurkin A plus alisertib significantly reduced polyploidy to vehicle control levels. Our in vitro and in vivo studies show that Aurkin A synergizes with alisertib and significantly decreases the alisertib dose needed to disrupt polyploidy while increasing apoptosis in DLBCL cells.

Targeting AURKA to induce synthetic lethality in CREBBP-deficient B-cell malignancies via attenuation of MYC expression

Loss-of-function mutations in CREBBP, which encodes for a histone acetyltransferase, occur frequently in B-cell malignancies, highlighting CREBBP deficiency as an attractive therapeutic target. Using established isogenic cell models, we demonstrated that CREBBP-deficient cells are selectively vulnerable to AURKA inhibition. Mechanistically, we found that co-targeting CREBBP and AURKA suppressed MYC transcriptionally and post-translationally to induce replication stress and apoptosis. Inhibition of AURKA dramatically decreased MYC protein level in CREBBP-deficient cells, implying a dependency on AURKA to sustain MYC stability. Furthermore, in vivo studies showed that pharmacological inhibition of AURKA was efficacious in delaying tumor progression in CREBBP-deficient cells and was synergistic with CREBBP inhibitors in CREBBP-proficient cells. Our study sheds light on a novel synthetic lethal interaction between CREBBP and AURKA, indicating that targeting AURKA represents a potential therapeutic strategy for high-risk B-cell malignancies harboring CREBBP inactivating mutations.

Integrated Bioinformatic Analysis Reveals the Oncogenic, Survival, and Prognostic Characteristics of TPX2 in Hepatocellular Carcinoma

Targeting protein for Xenopus kinesin-like protein 2 (TPX2), a well-known mitotic protein, has been linked to carcinogenesis in several cancers. This study investigated the role of TPX2 in hepatocellular carcinoma (HCC) from various aspects using bioinformatic analyses. TPX2 expression and its prognostic value in pan-cancers were analyzed using SangerBox. TPX2 expression and its association with prognosis, immune infiltration, tumor mutations, and signaling pathways in HCC were analyzed using UALCAN, BoxKaplan-Meier Plotter, GEPIA, Human Protein Atlas, TIMER 2.0, and SangerBox. Genes co-expressed with TPX2 in HCC were analyzed using the HCCDB database, followed by functional enrichment using SangerBox. Clinical predictive models were established based on TPX2 and its co-expressed genes using the ACLBI database. TPX2 expression significantly increased in pan-cancers and was associated with survival in nearly half of the cancer types. High TPX2 expression has been linked to poor survival outcomes in patients with HCC. TPX2 expression was positively correlated with abundant infiltration of immune cells (including B cells, CD4 + /CD8 + T cells, macrophages, neutrophils, and dendritic cells), TP53 mutation, and carcinogenesis-related pathways, such as the PI3K/AKT/mTOR pathway, cellular response to hypoxia, and tumor proliferation signature. Nineteen genes were found to be co-expressed with TPX2 in HCC, and these genes showed close positive correlations and were mainly implicated in cell cycle-related functions. A prognostic model established using TPX2 and its expressed genes could stratify HCC patients into high- and low-risk groups, with a significantly shorter survival time in high-risk groups. The prognostic model performed well in predicting 1-, 3-, and 5-year survival of patients with HCC, with areas under the curve of 0.801, 0.725, and 0.711, respectively. TPX2 functions as an oncogene in HCC, and its high expression is detrimental to the survival of patients with HCC. Thus, TPX2 may be a prognostic biomarker and potential therapeutic target for HCC.

TPX2 overexpression promotes sensitivity to dasatinib in breast cancer by activating YAP transcriptional signaling

Chromosomal instability (CIN) is a hallmark of cancer aggressiveness, providing genetic plasticity and tumor heterogeneity that allows the tumor to evolve and adapt to stress conditions. CIN is considered a cancer therapeutic biomarker because healthy cells do not exhibit CIN. Despite recent efforts to identify therapeutic strategies related to CIN, the results obtained have been very limited. CIN is characterized by a genetic signature where a collection of genes, mostly mitotic regulators, are overexpressed in CIN-positive tumors, providing aggressiveness and poor prognosis. We attempted to identify new therapeutic strategies related to CIN genes by performing a drug screen, using cells that individually express CIN-associated genes in an inducible manner. We find that the overexpression of targeting protein for Xklp2 (TPX2) enhances sensitivity to the proto-oncogene c-Src (SRC) inhibitor dasatinib due to activation of the Yes-associated protein 1 (YAP) pathway. Furthermore, using breast cancer data from The Cancer Genome Atlas (TCGA) and a cohort of cancer-derived patient samples, we find that both TPX2 overexpression and YAP activation are present in a significant percentage of cancer tumor samples and are associated with poor prognosis; therefore, they are putative biomarkers for selection for dasatinib therapy.

Exploring the anticancer potential of fluoro flavone analogues: insights from molecular docking and dynamics studies with Aurora Kinase B

Aurora Kinase B belongs to the serine kinase family. It plays an essential role in cell division and participates in mitosis and chromatid segregation. Overexpression, polymorphism, and splicing variants in the protein lead to tumorigenesis, leading to cancer. Flavones belong to the class of flavonoids and are derived from plants and show anti-cancer activities. Fluoro flavones and their analogs are taken from the PubChem database, resulting in 3882 compounds which is 90% similar to the fluoro flavones. Lipinski's rule of five, REOS and PAINS drug-like filters were applied which resulted 2448 compounds. These compounds are docked with Aurora Kinase B using SP and XP modules of Glide software. The best binding scores for SP docking were - 9.153 kcal/mol for the compound with CID: 44298667, and XP docking was - 10.287 kcal/mol with CID: 101664315. Enrichment calculations were done using Aurora Kinase B's decoys to validate the docking result. The resulting R2 = 0.96 from enrichment calculations suggests that the docking protocol is valid. The SP and XP docking lead compounds and the Fluoro flavone were subjected to 100 ns MD simulation to probe the protein-ligand complex stability. Also, the binding free energies between the Aurora kinase B and lead compounds were computed by Prime MM/GBSA module. The result suggests that the lead compounds bind more strongly with Aurora Kinase B than the Fluoro flavone. These lead compounds can be further evaluated in vitro and in vivo and can be used as future novel drugs for the curation of cancer.

Graphical abstract

A comprehensive review on role of Aurora kinase inhibitors (AKIs) in cancer therapeutics

Aurora kinases (AURKs) are a family of serine /threonine protein kinases that have a crucial role in cell cycle process mainly in the event of chromosomal segregation, centrosome maturation and cytokinesis. The family consists of three members including Aurora kinase A (AURK-A), Aurora kinase B (AURK-B) and Aurora kinase C (AURK-C). All AURKs contain a conserved kinase domain for their activity but differ in their cellular localization and functions. AURK-A and AURK-B are expressed mainly in somatic cells while the expression of AURK-C is limited to germ cells. AURK-A promotes G2 to M transition of cell cycle by controlling centrosome maturation and mitotic spindle assembly. AURK-B and AURK-C form the chromosome passenger complex (CPC) that ensures proper chromosomal alignments and segregation. Aberrant expression of AURK-A and AURK-B has been detected in several solid tumours and malignancies. Hence, they have become an attractive therapeutic target against cancer. The first part of this review focuses on AURKs structure, functions, subcellular localization, and their role in tumorigenesis. The review also highlights the functional and clinical impact of selective as well as pan kinase inhibitors. Currently, >60 compounds that target AURKs are in preclinical and clinical studies. The drawbacks of existing inhibitors like selectivity, drug resistance and toxicity have also been addressed. Since, majority of inhibitors are Aurora kinase inhibitor (AKI) type-1 that bind to the active (DFGin and Cin) conformation of the kinase, this information may be utilized to design highly selective kinase inhibitors that can be combined with other therapeutic agents for better clinical outcomes.

AurkA/TPX2 co-overexpression in nontransformed cells promotes genome instability through induction of chromosome mis-segregation and attenuation of the p53 signalling pathway

The Aurora-A kinase (AurkA) and its major regulator TPX2 (Targeting Protein for Xklp2) are key mitotic players frequently co-overexpressed in human cancers, and the link between deregulation of the AurkA/TPX2 complex and tumourigenesis is actively investigated. Chromosomal instability, one of the hallmarks of cancer related to the development of intra-tumour heterogeneity, metastasis and chemo-resistance, has been frequently associated with TPX2-overexpressing tumours. In this study we aimed to investigate the actual contribution to chromosomal instability of deregulating the AurkA/TPX2 complex, by overexpressing it in nontransformed hTERT RPE-1 cells. Our results show that overexpression of both AurkA and TPX2 results in increased AurkA activation and severe mitotic defects, compared to AurkA overexpression alone. We also show that AurkA/TPX2 co-overexpression yields increased aneuploidy in daughter cells and the generation of micronucleated cells. Interestingly, the p53/p21 axis response is impaired in AurkA/TPX2 overexpressing cells subjected to different stimuli; consistently, cells acquire increased ability to proliferate after independent induction of mitotic errors, i.e. following nocodazole treatment. Based on our observation that increased levels of the AurkA/TPX2 complex affect chromosome segregation fidelity and interfere with the activation of a pivotal surveillance mechanism in response to altered cell division, we propose that co-overexpression of AurkA and TPX2 per se represents a condition promoting the generation of a genetically unstable context in nontransformed human cells.

Explore Key Genes and Mechanisms Involved in Colon Cancer Progression Based on Bioinformatics Analysis

To explore underlying mechanisms related to the progression of colon cancer and identify hub genes associated with the prognosis of patients with colon cancer. GSE10950 and GSE62932 were downloaded from the Gene Expression Omnibus (GEO) database. GEO2R was utilized to screen out the differentially expressed genes (DEGs). Gene ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were conducted on DEGs. Moreover, STRING and Cytoscape software were utilized for establishing the network of protein-protein interaction (PPI) and identifying hub genes. Afterward, data from The Cancer Genome Atlas (TCGA) was utilized for identifying prognosis-related hub genes by Kaplan-Meier survival analysis. Colon cancer cell line LOVO and human normal intestinal epithelial cell line NCM-460 were exploited to demonstrate the differential expression of selected hub genes through RT-qPCR and western blot. The LOVO cells were transfected to regulate expressions of prognosis-associated genes, followed by exploring the effects of those genes on prognosis by Cell Counting Kit-8 assay and colony-forming assay for cancer cell proliferation, cell scratch test and transwell migration assay for cancer cell migration and Annexin V-PE/7-AAD double staining as well as flow cytometry for cancer cell apoptosis. In this study, 266 common DEGs were obtained from the intersection of two datasets. The GO analysis suggested the common DEGs mainly participated in the one-carbon metabolic process, cell cycle G2/M phase transition, organelle fission, cell cycle phase transition regulation, and regulation of mitotic cell cycle phase transition. The KEGG analysis demonstrated the common DEGs were related to the p53 signaling pathway, nitrogen metabolism, mineral absorption, and cell cycle. 10 hub genes including CCNB1, KIF4A, TPX2, MT1F, PRC1, PLK4, CALD1, MMP9, CLCA1, and MMP1 were identified and CCNB1, CLCA1, and PLK4 were prognosis-related. Increased expression of CCNB1, CLCA1, and PLK4 restrained proliferation as well as migration of cancer cells and induced apoptosis of cancer cells. CCNB1, KIF4A, TPX2, MT1F, PRC1, PLK4, CALD1, MMP9, CLCA1, and MMP1 were identified as hub genes and CCNB1, CLCA1, and PLK4 could inhibit the progression of colon cancer through inhibiting proliferation as well as migration of the cancer cell and promoting apoptosis of cancer cell.

Identification of hub genes and microRNAs with prognostic values in esophageal cancer by integrated analysis

Esophageal cancer (EC) is the eighth most common cancer in the world, and the sixth most common cause of cancer-related mortality. The aim of the present study was to identify cell and molecular mechanisms involved in EC, and to provide the potential targets for diagnosis and treatment. Here, a microarray dataset (GSE20347) was screened to find differentially expressed genes (DEGs). Different bioinformatic methods were used to analyze the identified DEGs. The up-regulated DEGs were significantly involved in different biological processes and pathways including extracellular matrix organization and ECM-receptor interaction. FN1, CDK1, AURKA, TOP2A, FOXM1, BIRC5, CDC6, UBE2C, TTK, and TPX2 were identified as the most important genes among the up-regulated DEGs. Our analysis showed that has-miR-29a-3p, has-miR-29b-3p, has-miR-29c-3p, and has-miR-767-5p had the largest number of common targets among the up-regulated DEGs. These findings strengthen the understanding of EC development and progression, as well as representing potential markers for EC diagnosis and treatment.

AurkA nuclear localization is promoted by TPX2 and counteracted by protein degradation

The AurkA kinase is a well-known mitotic regulator, frequently overexpressed in tumors. The microtubule-binding protein TPX2 controls AurkA activity, localization, and stability in mitosis. Non-mitotic roles of AurkA are emerging, and increased nuclear localization in interphase has been correlated with AurkA oncogenic potential. Still, the mechanisms leading to AurkA nuclear accumulation are poorly explored. Here, we investigated these mechanisms under physiological or overexpression conditions. We observed that AurkA nuclear localization is influenced by the cell cycle phase and nuclear export, but not by its kinase activity. Importantly, AURKA overexpression is not sufficient to determine its accumulation in interphase nuclei, which is instead obtained when AURKA and TPX2 are co-overexpressed or, to a higher extent, when proteasome activity is impaired. Expression analyses show that AURKA, TPX2, and the import regulator CSE1L are co-overexpressed in tumors. Finally, using MCF10A mammospheres we show that TPX2 co-overexpression drives protumorigenic processes downstream of nuclear AurkA. We propose that AURKA/TPX2 co-overexpression in cancer represents a key determinant of AurkA nuclear oncogenic functions.

TPX2 Amplification-Driven Aberrant Mitosis in Culture Adapted Human Embryonic Stem Cells with gain of 20q11.21

BACKGROUND

Despite highly effective machinery for the maintenance of genome integrity in human embryonic stem cells (hESCs), the frequency of genetic aberrations during in-vitro culture has been a serious issue for future clinical applications.

METHOD

By passaging hESCs over a broad range of timepoints (up to 6 years), the isogenic hESC lines with different passage numbers with distinct cellular characteristics, were established.

RESULT

We found that mitotic aberrations, such as the delay of mitosis, multipolar centrosomes, and chromosome mis-segregation, were increased in parallel with polyploidy compared to early-passaged hESCs (EP-hESCs) with normal copy number. Through high-resolution genome-wide approaches and transcriptome analysis, we found that culture adapted-hESCs with a minimal amplicon in chromosome 20q11.21 highly expressed TPX2, a key protein for governing spindle assembly and cancer malignancy. Consistent with these findings, the inducible expression of TPX2 in EP-hESCs reproduced aberrant mitotic events, such as the delay of mitotic progression, spindle stabilization, misaligned chromosomes, and polyploidy.

CONCLUSION

These studies suggest that the increased transcription of TPX2 in culture adapted hESCs could contribute to an increase in aberrant mitosis due to altered spindle dynamics.

Emerging roles of Aurora-A kinase in cancer therapy resistance

Aurora kinase A (Aurora-A), a serine/threonine kinase, plays a pivotal role in various cellular processes, including mitotic entry, centrosome maturation and spindle formation. Overexpression or gene-amplification/mutation of Aurora-A kinase occurs in different types of cancer, including lung cancer, colorectal cancer, and breast cancer. Alteration of Aurora-A impacts multiple cancer hallmarks, especially, immortalization, energy metabolism, immune escape and cell death resistance which are involved in cancer progression and resistance. This review highlights the most recent advances in the oncogenic roles and related multiple cancer hallmarks of Aurora-A kinase-driving cancer therapy resistance, including chemoresistance (taxanes, cisplatin, cyclophosphamide), targeted therapy resistance (osimertinib, imatinib, sorafenib, etc.), endocrine therapy resistance (tamoxifen, fulvestrant) and radioresistance. Specifically, the mechanisms of Aurora-A kinase promote acquired resistance through modulating DNA damage repair, feedback activation bypass pathways, resistance to apoptosis, necroptosis and autophagy, metastasis, and stemness. Noticeably, our review also summarizes the promising synthetic lethality strategy for Aurora-A inhibitors in RB1, ARID1A and MYC gene mutation tumors, and potential synergistic strategy for mTOR, PAK1, MDM2, MEK inhibitors or PD-L1 antibodies combined with targeting Aurora-A kinase. In addition, we discuss the design and development of the novel class of Aurora-A inhibitors in precision medicine for cancer treatment.

Identification and Screening of Novel Anti-Cancer Compounds for Aurora Kinase-A from Chemical Database

Aurora kinase is a group of enzymes that belongs to a serine-threonine family and plays a critical role in cellular division. Aurora Kinase A is overexpressed and distributed beyond the nucleus and is involved in tumorigenesis. Flavones are a class of flavonoids that are present in plants that show anticancer activity. Similar compounds of 2'Fluoroflavones are retrieved from the PubChem database. Then drug-like filters viz. REOS and PAINS were applied to remove toxic compounds using Canvas software, resulting in 3882 compounds being subjected to Glide docking with Aurora kinase A. The lead compounds were selected on the merit of hydrogen bonding, salt bridge, as well as pi-pi interactions, 4-(6-Fluoro-4-oxychromen-2yl) benzoic acid, has been found one of the best molecules from docking studies. The binding mode of the lead compound with AURKA reveals that the amino acid residues viz, Lys162, Ala213, and His280 are more important for binding with the binding affinity of -11.760 kcal/mol. The molecular dynamics simulations of 100 ns were done, which shows the mean RMSD value of 1.77 Å for all 3 complexes of the protein and Fluoroflavone and its analogs. This shows that Fluoroflavone and its 2 best analogs are tightly attached to the active sites and thus have conformational stability. Our finding suggests that 4-(6-fluoro-4-oxochromen-2-yl)benzoic acid and 4-(4-Oxochromen-2-yl)benzoate can be further used in vitro and in vivo experiments and can probably serve as a novel drug for cancer treatment.

PARP1 proximity proteomics reveals interaction partners at stressed replication forks

PARP1 mediates poly-ADP-ribosylation of proteins on chromatin in response to different types of DNA lesions. PARP inhibitors are used for the treatment of BRCA1/2-deficient breast, ovarian, and prostate cancer. Loss of DNA replication fork protection is proposed as one mechanism that contributes to the vulnerability of BRCA1/2-deficient cells to PARP inhibitors. However, the mechanisms that regulate PARP1 activity at stressed replication forks remain poorly understood. Here, we performed proximity proteomics of PARP1 and isolation of proteins on stressed replication forks to map putative PARP1 regulators. We identified TPX2 as a direct PARP1-binding protein that regulates the auto-ADP-ribosylation activity of PARP1. TPX2 interacts with DNA damage response proteins and promotes homology-directed repair of DNA double-strand breaks. Moreover, TPX2 mRNA levels are increased in BRCA1/2-mutated breast and prostate cancers, and high TPX2 expression levels correlate with the sensitivity of cancer cells to PARP-trapping inhibitors. We propose that TPX2 confers a mitosis-independent function in the cellular response to replication stress by interacting with PARP1.

Radiosensitivity in patients affected by ARPC1B deficiency: a new disease trait?

Actin-related protein 2/3 complex subunit 1B (ARPC1B) deficiency is a recently described inborn error of immunity (IEI) presenting with combined immunodeficiency and characterized by recurrent infections and thrombocytopenia. Manifestations of immune dysregulation, including colitis, vasculitis, and severe dermatitis, associated with eosinophilia, hyper-IgA, and hyper-IgE are also described in ARPC1B-deficient patients. To date, hematopoietic stem cell transplantation seems to be the only curative option for patients. ARPC1B is part of the actin-related protein 2/3 complex (Arp2/3) and cooperates with the Wiskott–Aldrich syndrome protein (WASp) in the regulation of the actin cytoskeleton remodeling and in driving double-strand break clustering for homology-directed repair. In this study, we aimed to investigate radiosensitivity (RS) in ARPC1B-deficient patients to assess whether it can be considered an additional disease trait. First, we performed trio-based next-generation-sequencing studies to obtain the ARPC1B molecular diagnosis in our index case characterized by increased RS, and then we confirmed, using three different methods, an increment of radiosensitivity in all enrolled ARPC1B-deficient patients. In particular, higher levels of chromatid-type aberrations and γH2AX foci, with an increased number of cells arrested in the G2/M-phase of the cell cycle, were found in patients’ cells after ionizing radiation exposition and radiomimetic bleomycin treatment. Overall, our data suggest increased radiosensitivity as an additional trait in ARPC1B deficiency and support the necessity to investigate this feature in ARPC1B patients as well as in other IEI with cytoskeleton defects to address specific clinical follow-up and optimize therapeutic interventions.

Identifying General Tumor and Specific Lung Cancer Biomarkers by Transcriptomic Analysis

The bioinformatic pipeline previously developed in our research laboratory is used to identify potential general and specific deregulated tumor genes and transcription factors related to the establishment and progression of tumoral diseases, now comparing lung cancer with other two types of cancer. Twenty microarray datasets were selected and analyzed separately to identify hub differentiated expressed genes and compared to identify all the deregulated genes and transcription factors in common between the three types of cancer and those unique to lung cancer. The winning DEGs analysis allowed to identify an important number of TFs deregulated in the majority of microarray datasets, which can become key biomarkers of general tumors and specific to lung cancer. A coexpression network was constructed for every dataset with all deregulated genes associated with lung cancer, according to DAVID’s tool enrichment analysis, and transcription factors capable of regulating them, according to oPOSSUM´s tool. Several genes and transcription factors are coexpressed in the networks, suggesting that they could be related to the establishment or progression of the tumoral pathology in any tissue and specifically in the lung. The comparison of the coexpression networks of lung cancer and other types of cancer allowed the identification of common connectivity patterns with deregulated genes and transcription factors correlated to important tumoral processes and signaling pathways that have not been studied yet to experimentally validate their role in lung cancer. The Kaplan–Meier estimator determined the association of thirteen deregulated top winning transcription factors with the survival of lung cancer patients. The coregulatory analysis identified two top winning transcription factors networks related to the regulatory control of gene expression in lung and breast cancer. Our transcriptomic analysis suggests that cancer has an important coregulatory network of transcription factors related to the acquisition of the hallmarks of cancer. Moreover, lung cancer has a group of genes and transcription factors unique to pulmonary tissue that are coexpressed during tumorigenesis and must be studied experimentally to fully understand their role in the pathogenesis within its very complex transcriptomic scenario. Therefore, the downstream bioinformatic analysis developed was able to identify a coregulatory metafirm of cancer in general and specific to lung cancer taking into account the great heterogeneity of the tumoral process at cellular and population levels.

TPX2 Serves as a Cancer Susceptibility Gene and Is Closely Associated with the Poor Prognosis of Endometrial Cancer

Background. Endometrial cancer (EC) is a common tumor of the genital tract that affects the female reproductive system but with only limited treatment options. We aimed to discover new prognostic biomarkers for EC. Methods. We used mRNA-seq data to detect differentially expressed genes (DEGs) between EC and control tissues. Detailed clinicopathological information was collected, and changes in the mRNA and protein levels of hub DEGs were analyzed in EC. Copy number variation (CNV) was also evaluated for its association with the pathogenesis of EC. Gene set enrichment analysis (GSEA) was conducted to enrich significant pathways driven by the hub genes. Cox regression analysis was used to select variables to create a nomogram. The nomogram was calibrated by applying the concordance index (C-index), and net benefits of the nomogram at different threshold probabilities were quantified using decision curve analysis (DCA). Results. Differential expression analysis identified 24 DEGs as potential risk factors for EC. Survival analysis revealed that TPX2 expression was related to worsening overall survival in patients with advanced EC. A high CNV was associated with the overexpression of TPX2; this suggested that modifications in the cell-cycle pathway might be crucial in the advancement of EC. Moreover, an individualized nomogram was developed for TPX2 incorporating clinical factors; this was also evaluated for its ability to predict EC. Calibration and DCA analyses confirmed the robustness and clinical usefulness of the nomogram. Conclusion. We offer novel insights into the pathogenesis and molecular mechanisms of EC. The overexpression of TPX2 was related to a poorer prognosis and could serve as a biomarker for predicting prognostic outcomes in EC patients.

The cell cycle-related genes RHAMM, AURKA, TPX2, PLK1, and PLK4 are associated with the poor prognosis of breast cancer patients

Breast cancer is the third most common type of cancer diagnosed. Cell cycle is a complex but highly organized and controlled process, in which normal cells sense mitogenic growth signals that instruct them to enter and progress through their cell cycle. This process culminates in cell division generating two daughter cells with identical amounts of genetic material. Uncontrolled proliferation is one of the hallmarks of cancer. In this study, we analyzed the expression of the cell cycle-related genes receptor for hyaluronan (HA)-mediated motility (RHAMM), AURKA, TPX2, PLK1, and PLK4 and correlated them with the prognosis in a collective of 3952 breast cancer patients. A high messenger RNA expression of all studied genes correlated with a poor prognosis. Stratifying the patients according to the expression of hormonal receptors, we found that in patients with estrogen and progesterone receptor-positive and human epithelial growth factor receptor 2-negative tumors, and Luminal A and Luminal B tumors, the expression of the five analyzed genes correlates with worse survival. qPCR analysis of a panel of breast cancer cell lines representative of major molecular subtypes indicated a predominant expression in the luminal subtype. In vitro experiments showed that radiation influences the expression of the five analyzed genes both in luminal and triple-negative model cell lines. Functional analysis of MDA-MB-231 cells showed that small interfering RNA knockdown of PLK4 and TPX2 and pharmacological inhibition of PLK1 had an impact on the cell cycle and colony formation. Looking for a potential upstream regulation by microRNAs, we observed a differential expression of RHAMM, AURKA, TPX2, PLK1, and PLK4 after transfecting the MDA-MB-231 cells with three different microRNAs. Survival analysis of miR-34c-5p, miR-375, and miR-142-3p showed a different impact on the prognosis of breast cancer patients. Our study suggests that RHAMM, AURKA, TPX2, PLK1, and PLK4 can be used as potential targets for treatment or as a prognostic value in breast cancer patients.

Identification of biomarkers related to glycolysis with weighted gene co-expression network analysis in oral squamous cell carcinoma

BACKGROUND

Oral squamous cell carcinoma (OSCC) is the most common tumor in the oral cavity and maxillofacial region. Increasing evidence suggests that aerobic glycolysis plays an important role in the occurrence, development, and prognosis of OSCC. Therefore, the identification of biomarkers related to glycolysis in OSCC represents considerable potential for improving its treatment.

METHODS

In the present study, a single-sample gene-set enrichment analysis (ssGSEA) algorithm with weighted gene co-expression network analysis (WGCNA) were used to quantify the degree of glycolysis and identify key modules with the greatest correlation with glycolysis.

RESULTS

Glycolytic scores significantly correlated with prognosis. In the key module 5 HUB genes were finally selected, which displayed a robust predictive effect. The expressions of key genes were associated with glycolysis.

CONCLUSIONS

The research comprehensively analyzed the glycolysis of OSCC and identified several biomarkers related to glycolysis. These biomarkers may represent potential therapeutic targets for future OSCC therapy.

SUMO proteases SENP3 and SENP5 spatiotemporally regulate the kinase activity of Aurora A

Precise chromosome segregation is mediated by a well-assembled mitotic spindle, which requires balance of the kinase activity of Aurora A (AurA, also known as AURKA). However, how this kinase activity is regulated remains largely unclear. Here, using in vivo and in vitro assays, we report that conjugation of SUMO2 with AurA at K258 in early mitosis promotes the kinase activity of AurA and facilitates the binding with its activator Bora. Knockdown of the SUMO proteases SENP3 and SENP5 disrupts the deSUMOylation of AurA, leading to increased kinase activity and abnormalities in spindle assembly and chromosome segregation, which could be rescued by suppressing the kinase activity of AurA. Collectively, these results demonstrate that SENP3 and SENP5 deSUMOylate AurA to render spatiotemporal control on its kinase activity in mitosis. This article has an associated First Person interview with the first author of the paper.

Elevating CDCA3 levels in non-small cell lung cancer enhances sensitivity to platinum-based chemotherapy

Platinum-based chemotherapy remains the cornerstone of treatment for most non-small cell lung cancer (NSCLC) cases either as maintenance therapy or in combination with immunotherapy. However, resistance remains a primary issue. Our findings point to the possibility of exploiting levels of cell division cycle associated protein-3 (CDCA3) to improve response of NSCLC tumours to therapy. We demonstrate that in patients and in vitro analyses, CDCA3 levels correlate with measures of genome instability and platinum sensitivity, whereby CDCA3^high tumours are sensitive to cisplatin and carboplatin. In NSCLC, CDCA3 protein levels are regulated by the ubiquitin ligase APC/C and cofactor Cdh1. Here, we identified that the degradation of CDCA3 is modulated by activity of casein kinase 2 (CK2) which promotes an interaction between CDCA3 and Cdh1. Supporting this, pharmacological inhibition of CK2 with CX-4945 disrupts CDCA3 degradation, elevating CDCA3 levels and increasing sensitivity to platinum agents. We propose that combining CK2 inhibitors with platinum-based chemotherapy could enhance platinum efficacy in CDCA3^low NSCLC tumours and benefit patients.

Kildey et al find that high levels of mitotic regulator CDCA3 correlates with sensitivity to platinum agents in non-small cell lung cancer patients and cell lines. They show that interfering with CDCA3 degradation through CK2 inhibition enhances CDCA3 levels and increases sensitivity to platinum agents suggesting a therapeutic route.

Chromosomal copy number heterogeneity predicts survival rates across cancers

Survival rates of cancer patients vary widely within and between malignancies. While genetic aberrations are at the root of all cancers, individual genomic features cannot explain these distinct disease outcomes. In contrast, intra-tumour heterogeneity (ITH) has the potential to elucidate pan-cancer survival rates and the biology that drives cancer prognosis. Unfortunately, a comprehensive and effective framework to measure ITH across cancers is missing. Here, we introduce a scalable measure of chromosomal copy number heterogeneity (CNH) that predicts patient survival across cancers. We show that the level of ITH can be derived from a single-sample copy number profile. Using gene-expression data and live cell imaging we demonstrate that ongoing chromosomal instability underlies the observed heterogeneity. Analysing 11,534 primary cancer samples from 37 different malignancies, we find that copy number heterogeneity can be accurately deduced and predicts cancer survival across tissues of origin and stages of disease. Our results provide a unifying molecular explanation for the different survival rates observed between cancer types.

Integrated genomics and comprehensive validation reveal drivers of genomic evolution in esophageal adenocarcinoma

Esophageal adenocarcinoma (EAC) is associated with a marked genomic instability, which underlies disease progression and development of resistance to treatment. In this study, we used an integrated genomics approach to identify a genomic instability signature. Here we show that elevated expression of this signature correlates with poor survival in EAC as well as three other cancers. Knockout and overexpression screens establish the relevance of these genes to genomic instability. Indepth evaluation of three genes (TTK, TPX2 and RAD54B) confirms their role in genomic instability and tumor growth. Mutational signatures identified by whole genome sequencing and functional studies demonstrate that DNA damage and homologous recombination are common mechanisms of genomic instability induced by these genes. Our data suggest that the inhibitors of TTK and possibly other genes identified in this study have potential to inhibit/reduce growth and spontaneous as well as chemotherapy-induced genomic instability in EAC and possibly other cancers.

Subodh Kumar et al. identify a gene signature correlated with genomic instability and poor survival in esophageal adenocarcinoma (EAC), using a combination of integrative genomic analysis of patient data and laboratory validation in cell line models and mice. They find that inhibitors of some of the identified proteins, including TTK, could be used to reduce genomic evolution as well as inhibit growth of EAC cells.

High nuclear TPX2 expression correlates with TP53 mutation and poor clinical behavior in a large breast cancer cohort, but is not an independent predictor of chromosomal instability

BACKGROUND

Targeting Protein for Xenopus Kinesin Like Protein 2 (TPX2) is a microtubule associated protein that functions in mitotic spindle assembly. TPX2 also localizes to the nucleus where it functions in DNA damage repair during S-phase. We and others have previously shown that TPX2 RNA levels are strongly associated with chromosomal instability (CIN) in breast and other cancers, and TPX2 RNA levels have been demonstrated to correlate with aggressive behavior and poor clinical outcome across a range of solid malignancies, including breast cancer.

METHODS

We perform TPX2 IHC on a cohort of 253 primary breast cancers and adopt a clinically amenable scoring system to separate tumors into low, intermediate, or high TPX2 expression. We then correlate TPX2 expression against diverse pathologic parameters and important measures of clinical outcome, including disease-specific and overall survival. We link TPX2 expression to TP53 mutation and evaluate whether TPX2 is an independent predictor of chromosomal instability (CIN).

RESULTS

We find that TPX2 nuclear expression strongly correlates with high grade morphology, elevated clinical stage, negative ER and PR status, and both disease-specific and overall survival. We also show that increased TPX2 nuclear expression correlates with elevated ploidy, supernumerary centrosomes, and TP53 mutation. TPX2 nuclear expression correlates with CIN via univariate analyses but is not independently predictive when compared to ploidy, Ki67, TP53 mutational status, centrosome number, and patient age.

CONCLUSIONS

Our findings demonstrate a strong correlation between TPX2 nuclear expression and aggressive tumor behavior, and show that TPX2 overexpression frequently occurs in the setting of TP53 mutation and elevated ploidy. However, TPX2 expression is not an independent predictor of CIN where it fails to outperform existing clinical and pathologic metrics.

Use of signals of positive and negative selection to distinguish cancer genes and passenger genes

A major goal of cancer genomics is to identify all genes that play critical roles in carcinogenesis. Most approaches focused on genes positively selected for mutations that drive carcinogenesis and neglected the role of negative selection. Some studies have actually concluded that negative selection has no role in cancer evolution. We have re-examined the role of negative selection in tumor evolution through the analysis of the patterns of somatic mutations affecting the coding sequences of human genes. Our analyses have confirmed that tumor suppressor genes are positively selected for inactivating mutations, oncogenes, however, were found to display signals of both negative selection for inactivating mutations and positive selection for activating mutations. Significantly, we have identified numerous human genes that show signs of strong negative selection during tumor evolution, suggesting that their functional integrity is essential for the growth and survival of tumor cells.

The Aurora-A/TPX2 Axis Directs Spindle Orientation in Adherent Human Cells by Regulating NuMA and Microtubule Stability

Mitotic spindle orientation is a crucial process that defines the axis of cell division, contributing to daughter cell positioning and fate, and hence to tissue morphogenesis and homeostasis.^1,2 The trimeric NuMA/LGN/Gαi complex, the major determinant of spindle orientation, exerts pulling forces on the spindle poles by anchoring astral microtubules (MTs) and dynein motors to the cell cortex.^3,4 Mitotic kinases contribute to correct spindle orientation by regulating nuclear mitotic apparatus protein (NuMA) localization,^5-7 among which the Aurora-A centrosomal kinase regulates NuMA targeting to the cell cortex in metaphase.^8,9 Aurora-A and its activator targeting protein for Xklp2 (TPX2) are frequently overexpressed in cancer,^10-12 raising the question as to whether spindle orientation is among the processes downstream the Aurora-A/TPX2 signaling axis altered under pathological conditions. Here, we investigated the role of TPX2 in the Aurora-A- and NuMA-dependent spindle orientation. We show that, in cultured adherent human cells, the interaction with TPX2 is required for Aurora-A to exert this function. We also show that Aurora-A, TPX2, and NuMA are part of a complex at spindle MTs, where TPX2 acts as a platform for Aurora-A regulation of NuMA. Interestingly, excess TPX2 does not influence NuMA localization but induces a "super-alignment" of the spindle axis with respect to the substrate, although an excess of Aurora-A induces spindle misorientation. These opposite effects are both linked to altered MT stability. Overall, our results highlight the importance of TPX2 for spindle orientation and suggest that spindle orientation is differentially sensitive to unbalanced levels of Aurora-A, TPX2, or the Aurora-A/TPX2 complex.

Rebelled epigenome: histone H3S10 phosphorylation and H3S10 kinases in cancer biology and therapy

Background

With the discovery that more than half of human cancers harbor mutations in chromatin proteins, deregulation of epigenetic mechanisms has been recognized a hallmark of malignant transformation. Post-translational modifications (PTMs) of histone proteins, as main components of epigenetic regulatory machinery, are also broadly accepted as therapeutic target. Current “epigenetic” therapies target predominantly writers, erasers and readers of histone acetylation and (to a lesser extent) methylation, leaving other types of PTMs largely unexplored. One of them is the phosphorylation of serine 10 on histone H3 (H3S10ph).

Main body

H3S10ph is emerging as an important player in the initiation and propagation of cancer, as it facilitates cellular malignant transformation and participates in fundamental cellular functions. In normal cells this histone mark dictates the hierarchy of additional histone modifications involved in the formation of protein binding scaffolds, transcriptional regulation, blocking repressive epigenetic information and shielding gene regions from heterochromatin spreading. During cell division, this mark is essential for chromosome condensation and segregation. It is also involved in the function of specific DNA–RNA hybrids, called R-loops, which modulate transcription and facilitate chromosomal instability. Increase in H3S10ph is observed in numerous cancer types and its abundance has been associated with inferior prognosis. Many H3S10-kinases, including MSK1/2, PIM1, CDK8 and AURORA kinases, have been long considered targets in cancer therapy. However, since these proteins also participate in other critical processes, including signal transduction, apoptotic signaling, metabolic fitness and transcription, their chromatin functions are often neglected.

Conclusions

H3S10ph and enzymes responsible for deposition of this histone modification are important for chromatin activity and oncogenesis. Epigenetic-drugs targeting this axis of modifications, potentially in combination with conventional or targeted therapy, provide a promising angle in search for knowledge-driven therapeutic strategies in oncology.

The role of Aurora-A in human cancers and future therapeutics

Aurora-A is a mitotic serine/threonine-protein kinase and an oncogene. In normal cells, Aurora-A appears from G2 phase and localizes at the centrosome, where it participates in centrosome replication, isolation and maturation. Aurora-A also maintains Golgi apparatus structure and spindle assembly. Aurora-A undergoes ubiquitination-mediated degradation after the cell division phase. Aurora-A is abnormally expressed in tumor cells and promotes cell proliferation by regulating mitotic substrates, such as PP1, PLK1, TPX2, and LAST2, and affects other molecules through a non-mitotic pathway to promote cell invasion and metastasis. Some molecules in tumor cells also indirectly act on Aurora-A to regulate tumor cells. Aurora-A also mediates resistance to chemotherapy and radiotherapy and is involved in tumor immunotherapy. Clinical trials of Aurora-A molecular inhibitors are currently underway, and clinical transformation is just around the corner.

Co-expression Network Analysis Identifies Fourteen Hub Genes Associated with Prognosis in Clear Cell Renal Cell Carcinoma

Renal cancer is a common genitourinary malignance, of which clear cell renal cell carcinoma (ccRCC) has high aggressiveness and leads to most cancer-related deaths. Identification of sensitive and reliable biomarkers for predicting tumorigenesis and progression has great significance in guiding the diagnosis and treatment of ccRCC. Here, we identified 2397 common differentially expressed genes (DEGs) using paired normal and tumor ccRCC tissues from GSE53757 and The Cancer Genome Atlas (TCGA). Then, we performed weighted gene co-expression network analysis and protein-protein interaction network analysis, 17 candidate hub genes were identified. These candidate hub genes were further validated in GSE36895 and Oncomine database and 14 real hub genes were identified. All the hub genes were up-regulated and significantly positively correlated with pathological stage and histologic grade of ccRCC. Survival analysis showed that the higher expression level of each hub gene tended to predict a worse clinical outcome. ROC analysis showed that all the hub genes can accurately distinguish between tumor and normal samples, and between early stage and advanced stage ccRCC. Moreover, all the hub genes were positively associated with distant metastasis, lymph node infiltration, tumor recurrence and the expression of MKi67, suggesting these genes might promote tumor proliferation, invasion and metastasis. Furthermore, the functional annotation demonstrated that most genes were enriched in cell-cycle related biological function. In summary, our study identified 14 potential biomarkers for predicting tumorigenesis and progression, which might contribute to early diagnosis, prognosis prediction and therapeutic intervention.

Aurora A inhibition disrupts chromosome condensation and spindle assembly during the first embryonic division in pigs

As common overexpression of Aurora A in various tumours, much attention has focused on its function in inducing cancer, and its value in cancer therapeutics, considerably less is known regarding its role in the first cleavage division of mammalian embryos. Here, we highlight an indispensable role of Aurora A during the first mitotic division progression of pig embryos just after meiosis. The expression and spatiotemporal localization of Aurora A were initially assessed in pig embryos during the first mitotic division by Western blot analysis and indirect immunofluorescent staining. Then, the potential role of Aurora A was further evaluated using a highly selective Aurora A inhibitor, MLN8054, during this mitotic progression in pig embryos. Aurora A was found to express and exhibit a specific dynamic intracellular localization pattern during the first mitotic division in pig embryos. Aurora A was diffused in the cytoplasm at the prophase stage, and then exhibited a dynamic intracellular localization which was tightly associated with the chromosome and spindle dynamics throughout subsequent mitotic phases. Inhibition of Aurora A by MLN8054 treatment led to the failure of the first cleavage, with the majority of embryos being arrested in prophase of the mitotic division. Further subcellular structure examination showed that Aurora A inhibition not only led to the failure of spindle microtubule assembly, but also resulted in severe defects in chromosome condensation, accompanied by an obvious decrease in p-TACC3(S558) expression during the prophase of the first mitosis. Together, these results illustrated that Aurora A is crucial for both spindle assembly and chromosome condensation during the first mitotic division in pig embryos, and that the regulation of Aurora A may be associated with its effects on p-TACC3(S558) expression.

Screening and identification of genes associated with cell proliferation in cholangiocarcinoma

Cholangiocarcinoma (CCA), an aggressive tumor with poor prognosis, is a malignant cancer with increasing incidence and mortality rates. It is important to survey crucial genes in CCA to find and design potential drug targets, especially for those genes associated with cell proliferation that is a key biological process in tumorgenesis. Herein, we surveyed genes associated with cell proliferation via a comprehensive pan-cancer analysis. Candidate genes were further analyzed using multiple approaches, including cross-analysis from diverse molecular levels, examination of potential function and interactions, and additional experimental validation. We primarily screened 15 potential genes based on 11 validated genes, and these 26 genes were further examined to delineate their biological functions and potential roles in cancer treatment. Several of them were involved synthetically lethal genetic interactions, especially for RECQL4, TOP2A, MKI67 and ASPM, indicating their potential roles in drug design and cancer treatment. Further experimental validation indicated that some genes were significantly upregulated in several cancer cell lines, implying their important roles in tumorigenesis. Our study identifies some genes associated with cell proliferation, which may be potential future targets in molecular targeted therapy.

Excess TPX2 Interferes with Microtubule Disassembly and Nuclei Reformation at Mitotic Exit

The microtubule-associated protein TPX2 is a key mitotic regulator that contributes through distinct pathways to spindle assembly. A well-characterised function of TPX2 is the activation, stabilisation and spindle localisation of the Aurora-A kinase. High levels of TPX2 are reported in tumours and the effects of its overexpression have been investigated in cancer cell lines, while little is known in non-transformed cells. Here we studied TPX2 overexpression in hTERT RPE-1 cells, using either the full length TPX2 or a truncated form unable to bind Aurora-A, to identify effects that are dependent-or independent-on its interaction with the kinase. We observe significant defects in mitotic spindle assembly and progression through mitosis that are more severe when overexpressed TPX2 is able to interact with Aurora-A. Furthermore, we describe a peculiar, and Aurora-A-interaction-independent, phenotype in telophase cells, with aberrantly stable microtubules interfering with nuclear reconstitution and the assembly of a continuous lamin B1 network, resulting in daughter cells displaying doughnut-shaped nuclei. Our results using non-transformed cells thus reveal a previously uncharacterised consequence of abnormally high TPX2 levels on the correct microtubule cytoskeleton remodelling and G1 nuclei reformation, at the mitosis-to-interphase transition.

Identification of candidate cancer drivers by integrative Epi-DNA and Gene Expression (iEDGE) data analysis

The emergence of large-scale multi-omics data warrants method development for data integration. Genomic studies from cancer patients have identified epigenetic and genetic regulators - such as methylation marks, somatic mutations, and somatic copy number alterations (SCNAs), among others - as predictive features of cancer outcome. However, identification of "driver genes" associated with a given alteration remains a challenge. To this end, we developed a computational tool, iEDGE, to model cis and trans effects of (epi-)DNA alterations and identify potential cis driver genes, where cis and trans genes denote those genes falling within and outside the genomic boundaries of a given (epi-)genetic alteration, respectively. iEDGE first identifies the cis and trans gene expression signatures associated with the presence/absence of a particular epi-DNA alteration across samples. It then applies tests of statistical mediation to determine the cis genes predictive of the trans gene expression. Finally, cis and trans effects are annotated by pathway enrichment analysis to gain insights into the underlying regulatory networks. We used iEDGE to perform integrative analysis of SCNAs and gene expression data from breast cancer and 18 additional cancer types included in The Cancer Genome Atlas (TCGA). Notably, cis gene drivers identified by iEDGE were found to be significantly enriched for known driver genes from multiple compendia of validated oncogenes and tumor suppressors, suggesting that the remainder are of equal importance. Furthermore, predicted drivers were enriched for functionally relevant cancer genes with amplification-driven dependencies, which are of potential prognostic and therapeutic value. All the analyses results are accessible at https://montilab.bu.edu/iEDGE. In summary, integrative analysis of SCNAs and gene expression using iEDGE successfully identified known cancer driver genes and putative cancer therapeutic targets across 19 cancer types in the TCGA. The proposed method can easily be applied to the integration of gene expression profiles with other epi-DNA assays in a variety of disease contexts.

A Comprehensive Review on MAPK: A Promising Therapeutic Target in Cancer

The mitogen-activated protein kinase (MAPK) pathway is an important bridge in the switch from extracellular signals to intracellular responses. Alterations of signaling cascades are found in various diseases, including cancer, as a result of genetic and epigenetic changes. Numerous studies focused on both the homeostatic and the pathologic conduct of MAPK signaling; however, there is still much to be deciphered in terms of regulation and action models in both preclinical and clinical research. MAPK has implications in the response to cancer therapy, particularly the activation of the compensatory pathways in response to experimental MAPK inhibition. The present paper discusses new insights into MAPK as a complex cell signaling pathway with roles in the sustenance of cellular normal conduit, response to cancer therapy, and activation of compensatory pathways. Unfortunately, most MAPK inhibitors trigger resistance due to the activation of compensatory feed-back loops in tumor cells and tumor microenvironment components. Therefore, novel combinatorial therapies have to be implemented for cancer management in order to restrict the possibility of alternative pathway activation, as a perspective for developing novel therapies based on integration in translational studies.

Defining and Modulating 'BRCAness'

The concept of 'BRCAness' defines the pathogenesis and vulnerability of multiple cancers. The canonical definition of BRCAness is a defect in homologous recombination repair, mimicking BRCA1 or BRCA2 loss. In turn, BRCA-deficient cells utilize error-prone DNA-repair pathways, causing increased genomic instability, which may be responsible for their sensitivity to DNA damaging agents and poly-(ADP)-ribose polymerase inhibitors (PARPis). However, recent work has expanded the mechanistic basis of BRCAness, to include defects in replication fork protection (RFP). Here, we broaden the definition of BRCAness to include RFP and regulatory mechanisms that cause synthetic lethality with PARPis. We highlight these recent discoveries, which include mechanisms of RFP regulation, DNA damage checkpoint proteins, as well as kinases that regulate BRCA1/2 function. Importantly, many of these emerging mechanisms may be targeted for inhibition with small molecule inhibitors, thus inducing BRCAness in a much larger subset of BRCA-proficient tumors, with significant translational potential.

GECKO is a genetic algorithm to classify and explore high throughput sequencing data

Comparative analysis of high throughput sequencing data between multiple conditions often involves mapping of sequencing reads to a reference and downstream bioinformatics analyses. Both of these steps may introduce heavy bias and potential data loss. This is especially true in studies where patient transcriptomes or genomes may vary from their references, such as in cancer. Here we describe a novel approach and associated software that makes use of advances in genetic algorithms and feature selection to comprehensively explore massive volumes of sequencing data to classify and discover new sequences of interest without a mapping step and without intensive use of specialized bioinformatics pipelines. We demonstrate that our approach called GECKO for GEnetic Classification using k-mer Optimization is effective at classifying and extracting meaningful sequences from multiple types of sequencing approaches including mRNA, microRNA, and DNA methylome data.