Whole-genome doubling in tissues and tumors.

The overwhelming majority of proliferating somatic human cells are diploid, and this genomic state is typically maintained across successive cell divisions. However, failures in cell division can induce a whole-genome doubling (WGD) event, in which diploid cells transition to a tetraploid state. While some WGDs are developmentally programmed to produce nonproliferative tetraploid cells with specific cellular functions, unscheduled WGDs can be catastrophic: erroneously arising tetraploid cells are ill-equipped to cope with their doubled cellular and chromosomal content and quickly become genomically unstable and tumorigenic. Deciphering the genetics that underlie the genesis, physiology, and evolution of whole-genome doubled (WGD+) cells may therefore reveal therapeutic avenues to selectively eliminate pathological WGD+ cells.

Cancer recurrence and lethality are enabled by enhanced survival and reversible cell cycle arrest of polyaneuploid cells.

We present a unifying theory to explain cancer recurrence, therapeutic resistance, and lethality. The basis of this theory is the formation of simultaneously polyploid and aneuploid cancer cells, polyaneuploid cancer cells (PACCs), that avoid the toxic effects of systemic therapy by entering a state of cell cycle arrest. The theory is independent of which of the classically associated oncogenic mutations have already occurred. PACCs have been generally disregarded as senescent or dying cells. Our theory states that therapeutic resistance is driven by PACC formation that is enabled by accessing a polyploid program that allows an aneuploid cancer cell to double its genomic content, followed by entry into a nondividing cell state to protect DNA integrity and ensure cell survival. Upon removal of stress, e.g., chemotherapy, PACCs undergo depolyploidization and generate resistant progeny that make up the bulk of cancer cells within a tumor.

Cancer cells employ an evolutionarily conserved polyploidization program to resist therapy.

Unusually large cancer cells with abnormal nuclei have been documented in the cancer literature since 1858. For more than 100 years, they have been generally disregarded as irreversibly senescent or dying cells, too morphologically misshapen and chromatin too disorganized to be functional. Cell enlargement, accompanied by whole genome doubling or more, is observed across organisms, often associated with mitigation strategies against environmental change, severe stress, or the lack of nutrients. Our comparison of the mechanisms for polyploidization in other organisms and non-transformed tissues suggest that cancer cells draw from a conserved program for their survival, utilizing whole genome doubling and pausing proliferation to survive stress. These polyaneuploid cancer cells (PACCs) are the source of therapeutic resistance, responsible for cancer recurrence and, ultimately, cancer lethality.

Characteristics of DNA macro-alterations in breast cancer with liver metastasis before treatment.

BACKGROUND: Whole-genome doubling (WGD) has been observed in 30% of cancers, followed by a highly complex rearranged karyotype unfavourable to breast cancer's outcome. However, the macro-alterations that characterise liver metastasis in breast cancer(BC) are poorly understood. Here, we conducted a whole-genome sequencing analysis of liver metastases to explore the status and the time frame model of these macro-alterations in pre-treatment patients with metastatic breast cancer. RESULTS: Whole-genome sequencing was conducted in 11 paired primary tumours, lymph node metastasis, and liver metastasis fresh samples from four patients with late-stage breast cancer. We also chose five postoperative frozen specimens from patients with early-stage breast cancer before any treatment as control. Surprisingly, all four liver metastasis samples were classified as WGD + . However, the previous study reported that WGD happened in 30% of cancers and 2/5 in our early-stage samples. WGD was not observed in the two separate primary tumours and one lymph node metastasis of one patient with metastatic BC, but her liver metastasis showed an early burst of bi-allelic copy number gain. The phylogenetic tree proves her 4 tumour samples were the polyclonal origin and only one WGD + clone metastasis to the liver. Another 3 metastatic BC patients' primary tumour and lymph node metastasis experienced WGD as well as liver metastasis, and they all showed similar molecular time-frame of copy number(CN) gain across locations within the same patient. These patients' tumours were of monoclonal origin, and WGD happened in a founding clone before metastasis, explaining that all samples share the CN-gain time frame. After WGD, the genomes usually face instability to evolve other macro-alterations. For example, a greater quantity and variety of complex structural variations (SVs) were detected in WGD + samples. The breakpoints were enriched in the chr17: 39 Mb-40 Mb tile, which contained the HER2 gene, resulting in the formation of tyfonas, breakage-fusion-bridge cycles, and double minutes. These complex SVs may be involved in the evolutionary mechanisms of the dramatic increase of HER2 copy number. CONCLUSION: Our work revealed that the WGD + clone might be a critical evolution step for liver metastasis and favoured following complex SV of breast cancer.

The evolution and expansion of RWP-RK gene family improve the heat adaptability of elephant grass (Pennisetum purpureum Schum.).

BACKGROUND: Along with global warming, resulting in crop production, exacerbating the global food crisis. Therefore, it is urgent to study the mechanism of plant heat resistance. However, crop resistance genes were lost due to long-term artificial domestication. By analyzing the potential heat tolerance genes and molecular mechanisms in other wild materials, more genetic resources can be provided for improving the heat tolerance of crops. Elephant grass (Pennisetum purpureum Schum.) has strong adaptability to heat stress and contains abundant heat-resistant gene resources. RESULTS: Through sequence structure analysis, a total of 36 RWP-RK members were identified in elephant grass. Functional analysis revealed their close association with heat stress. Four randomly selected RKDs (RKD1.1, RKD4.3, RKD6.6, and RKD8.1) were analyzed for expression, and the results showed upregulation under high temperature conditions, suggesting their active role in response to heat stress. The members of RWP-RK gene family (36 genes) in elephant grass were 2.4 times higher than that of related tropical crops, rice (15 genes) and sorghum (15 genes). The 36 RWPs of elephant grass contain 15 NLPs and 21 RKDs, and 73% of RWPs are related to WGD. Among them, combined with the DAP-seq results, it was found that RWP-RK gene family expansion could improve the heat adaptability of elephant grass by enhancing nitrogen use efficiency and peroxidase gene expression. CONCLUSIONS: RWP-RK gene family expansion in elephant grass is closely related to thermal adaptation evolution and speciation. The RKD subgroup showed a higher responsiveness than the NLP subgroup when exposed to high temperature stress. The promoter region of the RKD subgroup contains a significant number of MeJA and ABA responsive elements, which may contribute to their positive response to heat stress. These results provided a scientific basis for analyzing the heat adaptation mechanism of elephant grass and improving the heat tolerance of other crops.

Thermodynamic energetics underlying genomic instability and whole-genome doubling in cancer.

Genomic instability contributes to tumorigenesis through the amplification and deletion of cancer driver genes. DNA copy number (CN) profiling of ensembles of tumors allows a thermodynamic analysis of the profile for each tumor. The free energy of the distribution of CNs is found to be a monotonically increasing function of the average chromosomal ploidy. The dependence is universal across several cancer types. Surprisal analysis distinguishes two main known subgroups: tumors with cells that have or have not undergone whole-genome duplication (WGD). The analysis uncovers that CN states having a narrower distribution are energetically more favorable toward the WGD transition. Surprisal analysis also determines the deviations from a fully stable-state distribution. These deviations reflect constraints imposed by tumor fitness selection pressures. The results point to CN changes that are more common in high-ploidy tumors and thus support altered selection pressures upon WGD.

Autophagy modulating therapeutics inhibit ovarian cancer colony generation by polyploid giant cancer cells (PGCCs).

BACKGROUND: Genomic instability and chemoresistance can arise in cancer due to a unique form of plasticity: that of polyploid giant cancer cells (PGCCs). These cells form under the stress of chemotherapy and have higher than diploid chromosome content. PGCCs are able to then repopulate tumors through an asymmetric daughter cell budding process. PGCCs have been observed in ovarian cancer histology, including the deadly and common form high-grade serous ovarian carcinoma (HGSC). We previously discovered that drugs which disrupt the cellular recycling process of autophagy are uniquely efficacious in pre-clinical HGSC models. While autophagy induction has been associated with PGCCs, it has never been previously investigated if autophagy modulation interacts with the PGCC life cycle and this form of tumor cell plasticity. METHODS: CAOV3 and OVCAR3 ovarian cancer cell lines were treated with carboplatin or docetaxel to induce PGCC formation. Microscopy was used to characterize and quantify PGCCs formed by chemotherapy. Two clinically available drugs that inhibit autophagy, hydroxychloroquine and nelfinavir, and a clinically available activator of autophagy, rapamycin, were employed to test the effect of these autophagy modulators on PGCC induction and subsequent colony formation from PGCCs. Crystal violet-stained colony formation assays were used to quantify the tumor-repopulating stage of the PGCC life cycle. RESULTS: Autophagy inhibitors did not prevent PGCC formation in OVCAR3 or CAOV3 cells. Rapamycin did not induce PGCC formation on its own nor did it exacerbate PGCC formation by chemotherapy. However, hydroxychloroquine prevented efficient colony formation in CAOV3 PGCCs induced by carboplatin (27% inhibition) or docetaxel (41% inhibition), as well as in OVCAR3 cells (95% and 77%, respectively). Nelfinavir similarly prevented colony formation in CAOV3 PGCCs induced by carboplatin (64% inhibition) or docetaxel (94% inhibition) as well as in OVCAR3 cells (89% and 80%, respectively). Rapamycin surprisingly also prevented PGCC colony outgrowth (52-84% inhibition). CONCLUSIONS: While the autophagy previously observed to correlate with PGCC formation is unlikely necessary for PGCCs to form, autophagy modulating drugs severely impair the ability of HGSC PGCCs to form colonies. Clinical trials which utilize hydroxychloroquine, nelfinavir, and/or rapamycin after chemotherapy may be of future interest.

Poorly differentiated chordoma with whole-genome doubling evolving from a SMARCB1-deficient conventional chordoma: A case report.

Evolution of poorly differentiated chordoma from conventional chordoma has not been previously reported. We encountered a case of a poorly differentiated chordoma with evidence of whole-genome doubling arising from a SMARCB1-deficient conventional chordoma. The tumor presented as a destructive sacral mass in a 43-year-old man and was comprised of a highly cellular poorly differentiated chordoma with small, morphologically distinct nodules of conventional chordoma accounting for <5% of the total tumor volume. Immunohistochemistry (IHC) revealed both components were strongly reactive for brachyury and lacked normal staining for INI1. Single nucleotide polymorphism (SNP) array analysis identified multiple genomic imbalances in the conventional component, including deletions of 1p, 3p, and 22q (involving SMARCB1) and loss of chromosomes 5 and 15, while the poorly differentiated component exhibited the same aberrations at a more profound level with additional loss of chromosome 4, low level focal deletion of 17p (involving TP53), and tetraploidy. Homozygous deletion of SMARCB1 was present in both components. Fluorescence in situ hybridization (FISH) analysis confirmed the relevant deletions in both components as well as genome doubling in the poorly differentiated tumor. This case suggests that SMARCB1 loss is an early event in rare conventional chordomas that could potentially evolve into poorly differentiated chordoma through additional genomic aberrations such as genome doubling. Further studies with additional patients will be needed to determine if genome doubling is a consistent pathway for evolution of poorly differentiated chordoma.

A primary thymic adenocarcinoma with two components that traced distinct evolutionary trajectories.

Even though it is a rare subtype, identifying the genetic features of thymic adenocarcinoma is valuable for a multifaceted understanding of thymic epithelial tumors. We experienced a female patient with thymic adenocarcinoma associated with thymic cysts. The tumor consisted of a solid whitish lesion (lesion-1) and a large cystic lesion with small papillary nodules (lesion-2). Microscopically, lesion-1 exhibited poorly differentiated adenocarcinoma accompanying numerous inflammatory cell infiltrates, and lesion-2 (the nodules within the cystic lesion) exhibited enteric-type adenocarcinoma. Consistent with the histological difference, whole-exome sequencing revealed that these two components exhibited distinct genetic features, except for only a few shared mutations, including CDKN2A truncation. Lesion-1 exhibited microsatellite instability-high signature with high mutation burden, for which immune checkpoint inhibitors might apply; and lesion-2 exhibited whole-genome doubling with KRAS hotspot mutation. Our case presents novel genetic features of thymic adenocarcinoma and demonstrates that distinct mutational processes can be operative within a single tumor.

Distinguishing successive ancient polyploidy levels based on genome-internal syntenic alignment.

BACKGROUND: A basic tool for studying the polyploidization history of a genome, especially in plants, is the distribution of duplicate gene similarities in syntenically aligned regions of a genome. This distribution can usually be decomposed into two or more components identifiable by peaks, or local maxima, each representing a different polyploidization event. The distributions may be generated by means of a discrete time branching process, followed by a sequence divergence model. The branching process, as well as the inference of fractionation rates based on it, requires knowledge of the ploidy level of each event, which cannot be directly inferred from the pair similarity distribution. RESULTS: For a sequence of two events of unknown ploidy, either tetraploid, giving rise to whole genome doubling (WGD), or hexaploid, giving rise to whole genome tripling (WGT), we base our analysis on triples of similar genes. We calculate the probability of the four triplet types with origins in one or the other event, or both, and impose a mutational model so that the distribution resembles the original data. Using a ML transition point in the similarities between the two events as a discriminator for the hypothesized origin of each similarity, we calculate the predicted number of triplets of each type for each model combining WGT and/or WGD. This yields a predicted profile of triplet types for each model. We compare the observed and predicted triplet profiles for each model to confirm the polyploidization history of durian, poplar and cabbage. CONCLUSIONS: We have developed a way of inferring the ploidy of up to three successive WGD and/or WGT events by estimating the time of origin of each of the similarities in triples of genes. This may be generalized to a larger number of events and to higher ploidies.

Differential whole-genome doubling based signatures for improvement on clinical outcomes and drug response in patients with breast cancer.

Whole genome doublings (WGD), a hallmark of human cancer, is pervasive in breast cancer patients. However, the molecular mechanism of the complete impact of WGD on survival and treatment response in breast cancer remains unclear. To address this, we performed a comprehensive and systematic analysis of WGD, aiming to identify distinct genetic alterations linked to WGD and highlight its improvement on clinical outcomes and treatment response for breast cancer. A linear regression model along with weighted gene co-expression network analysis (WGCNA) was applied on The Cancer Genome Atlas (TCGA) dataset to identify critical genes related to WGD. Further Cox regression models with random selection were used to optimize the most useful prognostic markers in the TCGA dataset. The clinical implication of the risk model was further assessed through prognostic impact evaluation, tumor stratification, functional analysis, genomic feature difference analysis, drug response analysis, and multiple independent datasets for validation. Our findings revealed a high aneuploidy burden, chromosomal instability (CIN), copy number variation (CNV), and mutation burden in breast tumors exhibiting WGD events. Moreover, 247 key genes associated with WGD were identified from the distinct genomic patterns in the TCGA dataset. A risk model consisting of 22 genes was optimized from the key genes. High-risk breast cancer patients were more prone to WGD and exhibited greater genomic diversity compared to low-risk patients. Some oncogenic signaling pathways were enriched in the high-risk group, while primary immune deficiency pathways were enriched in the low-risk group. We also identified a risk gene, ANLN (anillin), which displayed a strong positive correlation with two crucial WGD genes, KIF18A and CCNE2. Tumors with high expression of ANLN were more prone to WGD events and displayed worse clinical survival outcomes. Furthermore, the expression levels of these risk genes were significantly associated with the sensitivities of BRCA cell lines to multiple drugs, providing valuable insights for targeted therapies. These findings will be helpful for further improvement on clinical outcomes and contribution to drug development in breast cancer.

Exploration of inhibitors targeting KIF18A with ploidy-specific lethality.

Currently, various antimitotic inhibitors applied in tumor therapy. However, these inhibitors exhibit targeted toxicity to some extent. As a motor protein, kinesin family member 18A (KIF18A) is crucial to spindle formation and is associated with tumors exhibiting ploidy-specific characteristics such as chromosomal aneuploidy, whole-genome doubling (WGD), and chromosomal instability (CIN). Differing from traditional antimitotic targets, KIF18A exhibits tumor-specific selectivity. The functional loss or attenuation of KIF18A results in vulnerability of tumor cells with ploidy-specific characteristics, with lesser effects on diploid cells. Research on inhibitors targeting KIF18A with ploidy-specific lethality holds significant importance. This review provides a brief overview of the regulatory mechanisms of the ploidy-specific lethality target KIF18A and the research advancements in its inhibitors, aiming to facilitate the development of KIF18A inhibitors.

scAbsolute: measuring single-cell ploidy and replication status.

Cancer cells often exhibit DNA copy number aberrations and can vary widely in their ploidy. Correct estimation of the ploidy of single-cell genomes is paramount for downstream analysis. Based only on single-cell DNA sequencing information, scAbsolute achieves accurate and unbiased measurement of single-cell ploidy and replication status, including whole-genome duplications. We demonstrate scAbsolute's capabilities using experimental cell multiplets, a FUCCI cell cycle expression system, and a benchmark against state-of-the-art methods. scAbsolute provides a robust foundation for single-cell DNA sequencing analysis across different technologies and has the potential to enable improvements in a number of downstream analyses.

Profiling Numerical and Structural Chromosomal Instability in Different Cancer Types.

Many cancers display whole chromosome instability (W-CIN) and structural chromosomal instability (S-CIN), referring to increased rates of acquiring numerically and structurally abnormal chromosome changes. This protocol provides detailed steps to analyze the W-CIN and S-CIN across cancer types, intending to leverage large-scale bulk sequencing and SNP array data complemented with the computational models to gain a better understanding of W-CIN and S-CIN.

Whole-Genome Doubling as a source of cancer: how, when, where, and why?

Chromosome instability is a well-known hallmark of cancer, leading to increased genetic plasticity of tumoral cells, which favors cancer aggressiveness, and poor prognosis. One of the main sources of chromosomal instability are events that lead to a Whole-Genome Duplication (WGD) and the subsequently generated cell polyploidy. In recent years, several studies showed that WGD occurs at the early stages of cell transformation, which allows cells to later become aneuploid, thus leading to cancer progression. On the other hand, other studies convey that polyploidy plays a tumor suppressor role, by inducing cell cycle arrest, cell senescence, apoptosis, and even prompting cell differentiation, depending on the tissue cell type. There is still a gap in understanding how cells that underwent WGD can overcome the deleterious effect on cell fitness and evolve to become tumoral. Some laboratories in the chromosomal instability field recently explored this paradox, finding biomarkers that modulate polyploid cells to become oncogenic. This review brings a historical view of how WGD and polyploidy impact cell fitness and cancer progression, and bring together the last studies that describe the genes helping cells to adapt to polyploidy.

Experimental Approaches to Generate and Isolate Human Tetraploid Cells.

Cancer cells are frequently affected by large-scale chromosome copy number changes, such as polyploidy or whole chromosome aneuploidy, and thus understanding the consequences of these changes is important for cancer research. In the past, it has been difficult to study the consequences of large-scale genomic changes, especially in pure isogenic populations. Here, we describe two methods to generate tetraploid cells induced either by cytokinesis failure or mitotic slippage. These treatments result in mixed population of diploids and tetraploids that can be analyzed directly. Alternatively, tetraploid populations can be established by single cell clone selection or by fluorescence activated cell sorting. These methods enable to analyze and compare the consequences of whole-genome doubling between the parental cell line, freshly arising tetraploid cells, and post-tetraploid aneuploid clones.

A case of vasculogenic mesenchymal tumor in the mediastinum: whole-exome sequencing reveals origin from pre-existing germ cell tumor.

Vasculogenic mesenchymal tumor (VMT), a primitive mesenchymal neoplasm enriched by various-sized atypical vessels, is a new entity that develops in mediastinal germ cell tumors (GCTs) with yolk sac tumor (YST) components after chemotherapy. Notably, patients with VMT in the residual GCT have increased risk of developing sarcomas or hematopoietic malignancies. Here, we report a late-teenage male patient with residual teratoma and high-grade VMT after chemotherapy for a mediastinal mixed GCT, including YST. Whole-exome sequencing revealed biallelic inactivation of TP53 and extensive copy number alterations that suggested whole-genome doubling. The biopsy tissue of the mixed GCT before chemotherapy exhibited overlapping genetic alterations to those in the VMT. Immunohistochemical analyses of the VMT showed that the abnormal vessels were positive for cytokeratin, glypican 3, EZH2, and IMP3. The findings that VMT inherits the genetic alterations of pre-existing mixed GCT and exhibits a partly YST-like immunophenotype might contribute to its clinical aggressiveness.

Whole genome doubling-induced the enrichment of H3K27me3 in genes carrying specific TEs in Aegilops tauschii.

Polyploidization plays important roles in the evolution and breeding of the common wheat. Aegilops tauschii, the D-genome progenitor of the common wheat, provides a valuable pool of resistance genes to multiple diseases. Extensive studies focus on the exploration of these genes for wheat improvement. However, few studies have unveiled alternations on genome-wide expression pattern and histone modifications induced by whole-genome doubling (WGD) process. In this study, we conducted transcriptome analysis for the diploid and tetraploid Ae. taushcii lines using the leaf and root tissues. Both lines tend to display similar tissue-specific pattern. Interestingly, we found that TEs located in genic regions were depleted of the repressive histone mark H3K27me3, whereas their adjacent chromatin was enriched with H3K27me3. The tetraploid line exhibited higher levels of H3K27me3 in those regions than the diploid line, particularly for genic regions associated with TEs of the long interspersed nuclear elements (LINEs), CACTA, PIF/Harbinger, Tc1/Mariner and unclassed DNA transposon. Surprisingly, the expression levels of these TEs cognate genes were negatively associated with the levels of H3K27me3 between the tetraploid and diploid lines, suggesting the five types of TEs located within genic regions might be involved in the regulation of the ploidy-related gene expression, possibly through differential enrichment of H3K27me3 in the genic regions. These findings will help to understand the potential role of specific types of TEs on transcription in response to WGD.