Unravelling epigenetic mechanisms in Cerastoderma edule genome: a comparison of healthy and neoplastic cockles.

Cancer is a multifaceted genetic disease characterized by the acquisition of several essential hallmarks. Notably, certain cancers exhibit horizontal transmissibility, observed across mammalian species and diverse bivalves, the latter referred to as hemic neoplasia. Within this complex landscape, epigenetic mechanisms such as histone modifications and cytosine methylation emerge as fundamental contributors to the pathogenesis of these transmissible cancers. Our study delves into the epigenetic landscape of Cerastoderma edule, focusing on whole-genome methylation and hydroxymethylation profiles in heathy specimens and transmissible neoplasias by means of Nanopore long-read sequencing. Our results unveiled a global hypomethylation in the neoplastic specimens compared to their healthy counterparts, emphasizing the role of DNA methylation in these tumorigenic processes. Furthermore, we verified that intragenic CpG methylation positively correlated with gene expression, emphasizing its role in modulating transcription in healthy and neoplastic cockles, as also highlighted by some up-methylated oncogenic genes. Hydroxymethylation levels were significantly more elevated in the neoplastic samples, particularly within satellites and complex repeats, likely related to structural functions. Additionally, our analysis also revealed distinct methylation and activity patterns in retrotransposons, providing additional insights into bivalve neoplastic processes. Altogether, these findings contribute to understanding the epigenetic dynamics of bivalve neoplasias and shed light on the roles of DNA methylation and hydroxymethylation in tumorigenesis. Understanding these epigenetic alterations holds promise for advancing our broader understanding of cancer epigenetics.

A three-gene expression score for predicting clinical benefit to anti-PD-1 blockade in advanced renal cell carcinoma.

In the advanced renal cell carcinoma (RCC) scenario, there are no consistent biomarkers to predict the clinical benefit patients derived from immune checkpoint blockade (ICB). Taking this into consideration, herein, we conducted a retrospective study in order to develop and validate a gene expression score for predicting clinical benefit to the anti-PD-1 antibody nivolumab in the context of patients diagnosed with advanced clear cell RCC enrolled in the CheckMate-009, CheckMate-010, and CheckMate-025 clinical trials. First, a three-gene expression score (3GES) with prognostic value for overall survival integrating HMGA1, NUP62, and ARHGAP42 transcripts was developed in a cohort of patients treated with nivolumab. Its prognostic value was then validated in the TCGA-KIRC cohort. Second, the predictive value for nivolumab was confirmed in a set of patients from the CheckMate-025 phase 3 clinical trial. Lastly, we explored the correlation of our 3GES with different clinical, molecular, and immune tumor characteristics. If the results of this study are definitively validated in other retrospective and large-scale, prospective studies, the 3GES will represent a valuable tool for guiding the design of ICB-based clinical trials in the aRCC scenario in the near future.

Dispersal history of SARS-CoV-2 in Galicia, Spain.

The dynamics of SARS-CoV-2 transmission are influenced by a variety of factors, including social restrictions and the emergence of distinct variants. In this study, we delve into the origins and dissemination of the Alpha, Delta, and Omicron variants of concern in Galicia, northwest Spain. For this, we leveraged genomic data collected by the EPICOVIGAL Consortium and from the GISAID database, along with mobility information from other Spanish regions and foreign countries. Our analysis indicates that initial introductions during the Alpha phase were predominantly from other Spanish regions and France. However, as the pandemic progressed, introductions from Portugal and the USA became increasingly significant. Notably, Galicia's major coastal cities emerged as critical hubs for viral transmission, highlighting their role in sustaining and spreading the virus. This research emphasizes the critical role of regional connectivity in the spread of SARS-CoV-2 and offers essential insights for enhancing public health strategies and surveillance measures.

Somatic evolution of marine transmissible leukemias in the common cockle, Cerastoderma edule.

Transmissible cancers are malignant cell lineages that spread clonally between individuals. Several such cancers, termed bivalve transmissible neoplasia (BTN), induce leukemia-like disease in marine bivalves. This is the case of BTN lineages affecting the common cockle, Cerastoderma edule, which inhabits the Atlantic coasts of Europe and northwest Africa. To investigate the evolution of cockle BTN, we collected 6,854 cockles, diagnosed 390 BTN tumors, generated a reference genome and assessed genomic variation across 61 tumors. Our analyses confirmed the existence of two BTN lineages with hemocytic origins. Mitochondrial variation revealed mitochondrial capture and host co-infection events. Mutational analyses identified lineage-specific signatures, one of which likely reflects DNA alkylation. Cytogenetic and copy number analyses uncovered pervasive genomic instability, with whole-genome duplication, oncogene amplification and alkylation-repair suppression as likely drivers. Satellite DNA distributions suggested ancient clonal origins. Our study illuminates long-term cancer evolution under the sea and reveals tolerance of extreme instability in neoplastic genomes.

MSIMEP: Predicting microsatellite instability from microarray DNA methylation tumor profiles.

Deficiency in DNA MMR activity results in tumors with a hypermutator phenotype, termed microsatellite instability (MSI). Beyond its utility in Lynch syndrome screening algorithms, today MSI has gained importance as predictive biomarker for various anti-PD-1 therapies across many different tumor types. Over the past years, many computational methods have emerged to infer MSI using either DNA- or RNA-based approaches. Considering this together with the fact that MSI-high tumors frequently exhibit a hypermethylated phenotype, herein we developed and validated MSIMEP, a computational tool for predicting MSI status from microarray DNA methylation tumor profiles of colorectal cancer samples. We demonstrated that MSIMEP optimized and reduced models have high performance in predicting MSI in different colorectal cancer cohorts. Moreover, we tested its consistency in other tumor types with high prevalence of MSI such as gastric and endometrial cancers. Finally, we demonstrated better performance of both MSIMEP models vis-à-vis a MLH1 promoter methylation-based one in colorectal cancer.

Clinical significance of baseline Pan-Immune-Inflammation Value and its dynamics in metastatic colorectal cancer patients under first-line chemotherapy.

Pan-Immune-Inflammation Value (PIV) has been recently proposed as a new blood-based prognostic biomarker in metastatic colorectal cancer (mCRC). Herein we aimed to validate its prognostic significance and to evaluate its utility for disease monitoring in patients with mCRC receiving first-line chemotherapy. We conducted a single-centre retrospective study involving 130 previously untreated mCRC patients under first-line standard chemotherapy in a real-world scenario. PIV was calculated as (neutrophil count × platelet count × monocyte count)/lymphocyte count at three different time-points: baseline, week 4 after therapy initiation, and at disease progression. We analyzed the influence of baseline PIV on overall survival (OS), progression-free survival (PFS), disease control rate (DCR), and overall response rate (ORR). We also explored the utility of PIV dynamics for disease monitoring. Baseline PIV high was significantly associated with worse OS in univariate [hazard ratio (HR) = 2.10, 95% CI, 1.41-3.15; p = 0.000299] and multivariate (HR = 1.82, 95% CI, 1.15-2.90; p = 0.011) analyses. Baseline PIV was also associated with worse PFS in univariate (HR = 2.04, 95% CI, 1.40-2.97; p = 0.000187) and multivariate (HR = 1.56, 95% CI, 1.05-2.31; p = 0.026) analyses. Baseline PIV was not correlated either with DCR or ORR. Regarding PIV dynamics, there was a statistically significant increase from week 4 to disease progression (p = 0.0003), which was at the expense of cases with disease control as best response (p < 0.0001). In conclusion, this study validates the prognostic significance of baseline PIV in patients with mCRC receiving first-line standard chemotherapy in a real-world scenario. Moreover, it suggests the potential utility of PIV monitoring to anticipate the disease progression among those patients who achieve initial disease control.

PanCancer analysis of somatic mutations in repetitive regions reveals recurrent mutations in snRNA U2.

Current somatic mutation callers are biased against repetitive regions, preventing the identification of potential driver alterations in these loci. We developed a mutation caller for repetitive regions, and applied it to study repetitive non protein-coding genes in more than 2200 whole-genome cases. We identified a recurrent mutation at position c.28 in the gene encoding the snRNA U2. This mutation is present in B-cell derived tumors, as well as in prostate and pancreatic cancer, suggesting U2 c.28 constitutes a driver candidate associated with worse prognosis. We showed that the GRCh37 reference genome is incomplete, lacking the U2 cluster in chromosome 17, preventing the identification of mutations in this gene. Furthermore, the 5'-flanking region of WDR74, previously described as frequently mutated in cancer, constitutes a functional copy of U2. These data reinforce the relevance of non-coding mutations in cancer, and highlight current challenges of cancer genomic research in characterizing mutations affecting repetitive genes.

Mitochondrial genome sequencing of marine leukaemias reveals cancer contagion between clam species in the Seas of Southern Europe.

Clonally transmissible cancers are tumour lineages that are transmitted between individuals via the transfer of living cancer cells. In marine bivalves, leukaemia-like transmissible cancers, called hemic neoplasia (HN), have demonstrated the ability to infect individuals from different species. We performed whole-genome sequencing in eight warty venus clams that were diagnosed with HN, from two sampling points located more than 1000 nautical miles away in the Atlantic Ocean and the Mediterranean Sea Coasts of Spain. Mitochondrial genome sequencing analysis from neoplastic animals revealed the coexistence of haplotypes from two different clam species. Phylogenies estimated from mitochondrial and nuclear markers confirmed this leukaemia originated in striped venus clams and later transmitted to clams of the species warty venus, in which it survives as a contagious cancer. The analysis of mitochondrial and nuclear gene sequences supports all studied tumours belong to a single neoplastic lineage that spreads in the Seas of Southern Europe.

In humans and other animals, cancer cells divide excessively, forming tumours or flooding the blood, but they rarely spread to other individuals. However, some animals, including dogs, Tasmanian devils and bivalve molluscs like clams, cockles and mussels, can develop cancers that are transmitted from one individual to another. Despite these cancers being contagious, each one originates in a single animal, meaning that even when the cancer has spread to many individuals, its origins can be traced through its DNA. Cancer contagion is rare, but transmissible cancers seem to be particularly common in the oceans. In fact, 7 types of contagious cancer have been described in bivalve species so far. These cancers are known as 'hemic neoplasias', and are characterized by the uncontrolled division of blood-like cells, which can be released by the host they developed in, and survive in ocean water. When these cells encounter individuals from the same species, they can infect them, causing them to develop hemic neoplasia too There are still many unanswered questions about contagious cancers in bivalves. For example, how many species do the cancers affect, and which species do the cancers originate in? To address these questions, Garcia-Souto, Bruzos, Díaz et al. gathered over 400 specimens of a species of clam called the warty venus clam from the coastlines of Europe and examined them for signs of cancer. Clams collected in two regions of Spain showed signs of hemic neoplasia: one of the populations was from the Balearic Islands in the Mediterranean Sea, while the other came from the Atlantic coast of northwestern Spain. Analyzing the genomes of the tumours from each population showed that the cancer cells from both regions had likely originated in the same animal, indicating that the cancer is contagious and had spread through different populations. The analysis also revealed that the cancer did not originally develop in warty venus clams: the cancer cells contained DNA from both warty venus clams and another species called striped venus clams. These two species live close together in the Mediterranean Sea, suggesting that the cancer started in a striped venus clam and then spread to a warty venus clam. To determine whether the cancer still affected both species, Garcia-Souto, Bruzos, Díaz et al. screened 200 striped venus clams from the same areas, but no signs of cancer were found in these clams. This suggests that currently the cancer only affects the warty venus clam. These findings confirm that contagious cancers can jump between clam species, which could be threat to the marine environment. The fact that the cancer was so similar in clams from the Atlantic coast and from the Mediterranean Sea, however, suggests that it may have emerged very recently, or that human activity helped it to spread from one place to another. If the latter is the case, it may be possible to prevent further spread of these sea-borne cancers through human intervention.

Aberrant integration of Hepatitis B virus DNA promotes major restructuring of human hepatocellular carcinoma genome architecture.

Most cancers are characterized by the somatic acquisition of genomic rearrangements during tumour evolution that eventually drive the oncogenesis. Here, using multiplatform sequencing technologies, we identify and characterize a remarkable mutational mechanism in human hepatocellular carcinoma caused by Hepatitis B virus, by which DNA molecules from the virus are inserted into the tumour genome causing dramatic changes in its configuration, including non-homologous chromosomal fusions, dicentric chromosomes and megabase-size telomeric deletions. This aberrant mutational mechanism, present in at least 8% of all HCC tumours, can provide the driver rearrangements that a cancer clone requires to survive and grow, including loss of relevant tumour suppressor genes. Most of these events are clonal and occur early during liver cancer evolution. Real-time timing estimation reveals some HBV-mediated rearrangements occur as early as two decades before cancer diagnosis. Overall, these data underscore the importance of characterising liver cancer genomes for patterns of HBV integration.

Reliability and validity of the reduced Spanish version of the Prenatal Breastfeeding Self-Efficacy Scale.

Breastfeeding self-efficacy is a relevant predictor of exclusive breastfeeding during the first six months of life. The Prenatal Breastfeeding Self-Efficacy Scale (PBSES) is a 20-item tool developed to determine breastfeeding self-efficacy during pregnancy. Our study aimed to assess the structural validity and psychometric characteristics of the PBSES and to explore item reduction according to the statistical criteria for parsimony and incremental validity. In this study, conducted in six hospitals in eastern Spain, we recruited 1183 women with healthy, full-term, single-birth newborns. Data on sociodemographic, breastfeeding-related variables, and the PBSES, were obtained from self-administered questionnaires during the third trimester of pregnancy, at postpartum discharge, and 5 months postpartum, and from a phone survey 12 months postpartum. Item reduction was conducted after revising the PBSES item floor and ceiling effects, interitem correlations, and item-score relationships with breastfeeding-status variables during follow-up. The factorial structure of the short form of the PBSES (PBSES-SF) was tested using both exploratory and confirmatory approaches. After item reduction, the confirmatory factor analysis of the 12 remaining items of the PBSES-SF revealed adequate fit statistics for a three-factor structure and a second-order factor. Internal consistency was measured using the Cronbach's α coefficient of the PBSES-SF (0.86). We provided evidence on the discriminant validity of the PBSES-SF by comparing its scores between known groups, convergent validity by examining its correlations with other variables, and predictive validity by assessing the association of PBSES-SF scores with breastfeeding behavior at critical points in time during the first postpartum year.

Subclone Eradication Analysis Identifies Targets for Enhanced Cancer Therapy and Reveals L1 Retrotransposition as a Dynamic Source of Cancer Heterogeneity.

Treatment-eradicated cancer subclones have been reported in leukemia and have recently been detected in solid tumors. Here we introduce Differential Subclone Eradication and Resistance (DSER) analysis, a method developed to identify molecular targets for improved therapy by direct comparison of genomic features of eradicated and resistant subclones in pre- and posttreatment samples from a patient with BRCA2-deficient metastatic prostate cancer. FANCI and EYA4 were identified as candidate DNA repair-related targets for converting subclones from resistant to eradicable, and RNAi-mediated depletion of FANCI confirmed it as a potential target. The EYA4 alteration was associated with adjacent L1 transposon insertion during cancer evolution upon treatment, raising questions surrounding the role of therapy in L1 activation. Both carboplatin and enzalutamide turned on L1 transposon machinery in LNCaP and VCaP but not in PC3 and 22Rv1 prostate cancer cell lines. L1 activation in LNCaP and VCaP was inhibited by the antiretroviral drug azidothymidine. L1 activation was also detected postcastration in LuCaP 77 and LuCaP 105 xenograft models and postchemotherapy in previously published time-series transcriptomic data from SCC25 head and neck cancer cells. In conclusion, DSER provides an informative intermediate step toward effective precision cancer medicine and should be tested in future studies, especially those including dramatic but temporary metastatic tumor regression. L1 transposon activation may be a modifiable source of cancer genomic heterogeneity, suggesting the potential of leveraging newly discovered triggers and blockers of L1 activity to overcome therapy resistance. SIGNIFICANCE: Differential analysis of eradicated and resistant subclones following cancer treatment identifies that L1 activity associated with resistance is induced by current therapies and blocked by the antiretroviral drug azidothymidine.

Haplotype-resolved diverse human genomes and integrated analysis of structural variation.

Long-read and strand-specific sequencing technologies together facilitate the de novo assembly of high-quality haplotype-resolved human genomes without parent-child trio data. We present 64 assembled haplotypes from 32 diverse human genomes. These highly contiguous haplotype assemblies (average minimum contig length needed to cover 50% of the genome: 26 million base pairs) integrate all forms of genetic variation, even across complex loci. We identified 107,590 structural variants (SVs), of which 68% were not discovered with short-read sequencing, and 278 SV hotspots (spanning megabases of gene-rich sequence). We characterized 130 of the most active mobile element source elements and found that 63% of all SVs arise through homology-mediated mechanisms. This resource enables reliable graph-based genotyping from short reads of up to 50,340 SVs, resulting in the identification of 1526 expression quantitative trait loci as well as SV candidates for adaptive selection within the human population.

CDKN2A deletion is a frequent event associated with poor outcome in patients with peripheral T-cell lymphoma not otherwise specified (PTCL-NOS).

Nodal peripheral T-cell lymphoma not otherwise specified (PTCL-NOS) remains a diagnosis encompassing a heterogenous group of PTCL cases not fitting criteria for more homogeneous subtypes. They are characterized by a poor clinical outcome when treated with anthracycline-containing regimens. A better understanding of their biology could improve prognostic stratification and foster the development of novel therapeutic approaches. Recent targeted and whole exome sequencing studies have shown recurrent copy number abnormalities (CNAs) with prognostic significance. Here, investigating 5 formalin-fixed, paraffin embedded cases of PTCL-NOS by whole genome sequencing (WGS), we found a high prevalence of structural variants and complex events, such as chromothripsis likely responsible for the observed CNAs. Among them, CDKN2A and PTEN deletions emerged as the most frequent aberration, as confirmed in a final cohort of 143 patients with nodal PTCL. The incidence of CDKN2A and PTEN deletions among PTCL-NOS was 46% and 26%, respectively. Furthermore, we found that co-occurrence of CDKN2A and PTEN deletions is an event associated with PTCL-NOS with absolute specificity. In contrast, these deletions were rare and never co-occurred in angioimmunoblastic and anaplastic lymphomas. CDKN2A deletion was associated with shorter overall survival in multivariate analysis corrected by age, IPI, transplant eligibility and GATA3 expression (adjusted HR =2.53; 95% CI 1.006-6.3; p=0.048). These data suggest that CDKN2A deletions may be relevant for refining the prognosis of PTCL-NOS and their significance should be evaluated in prospective trials.

Extensive heterogeneity in somatic mutation and selection in the human bladder.

The extent of somatic mutation and clonal selection in the human bladder remains unknown. We sequenced 2097 bladder microbiopsies from 20 individuals using targeted (n = 1914 microbiopsies), whole-exome (n = 655), and whole-genome (n = 88) sequencing. We found widespread positive selection in 17 genes. Chromatin remodeling genes were frequently mutated, whereas mutations were absent in several major bladder cancer genes. There was extensive interindividual variation in selection, with different driver genes dominating the clonal landscape across individuals. Mutational signatures were heterogeneous across clones and individuals, which suggests differential exposure to mutagens in the urine. Evidence of APOBEC mutagenesis was found in 22% of the microbiopsies. Sequencing multiple microbiopsies from five patients with bladder cancer enabled comparisons with cancer-free individuals and across histological features. This study reveals a rich landscape of mutational processes and selection in normal urothelium with large heterogeneity across clones and individuals.

VAV2 signaling promotes regenerative proliferation in both cutaneous and head and neck squamous cell carcinoma.

Regenerative proliferation capacity and poor differentiation are histological features usually linked to poor prognosis in head and neck squamous cell carcinoma (hnSCC). However, the pathways that regulate them remain ill-characterized. Here, we show that those traits can be triggered by the RHO GTPase activator VAV2 in keratinocytes present in the skin and oral mucosa. VAV2 is also required to maintain those traits in hnSCC patient-derived cells. This function, which is both catalysis- and RHO GTPase-dependent, is mediated by c-Myc- and YAP/TAZ-dependent transcriptomal programs associated with regenerative proliferation and cell undifferentiation, respectively. High levels of VAV2 transcripts and VAV2-regulated gene signatures are both associated with poor hnSCC patient prognosis. These results unveil a druggable pathway linked to the malignancy of specific SCC subtypes.

Author Correction: The evolutionary history of lethal metastatic prostate cancer.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Genomic landscape and chronological reconstruction of driver events in multiple myeloma.

The multiple myeloma (MM) genome is heterogeneous and evolves through preclinical and post-diagnosis phases. Here we report a catalog and hierarchy of driver lesions using sequences from 67 MM genomes serially collected from 30 patients together with public exome datasets. Bayesian clustering defines at least 7 genomic subgroups with distinct sets of co-operating events. Focusing on whole genome sequencing data, complex structural events emerge as major drivers, including chromothripsis and a novel replication-based mechanism of templated insertions, which typically occur early. Hyperdiploidy also occurs early, with individual trisomies often acquired in different chronological windows during evolution, and with a preferred order of acquisition. Conversely, positively selected point mutations, whole genome duplication and chromoplexy events occur in later disease phases. Thus, initiating driver events, drawn from a limited repertoire of structural and numerical chromosomal changes, shape preferred trajectories of evolution that are biologically relevant but heterogeneous across patients.

Characterizing Mutational Signatures in Human Cancer Cell Lines Reveals Episodic APOBEC Mutagenesis.

Multiple signatures of somatic mutations have been identified in cancer genomes. Exome sequences of 1,001 human cancer cell lines and 577 xenografts revealed most common mutational signatures, indicating past activity of the underlying processes, usually in appropriate cancer types. To investigate ongoing patterns of mutational-signature generation, cell lines were cultured for extended periods and subsequently DNA sequenced. Signatures of discontinued exposures, including tobacco smoke and ultraviolet light, were not generated in vitro. Signatures of normal and defective DNA repair and replication continued to be generated at roughly stable mutation rates. Signatures of APOBEC cytidine deaminase DNA-editing exhibited substantial fluctuations in mutation rate over time with episodic bursts of mutations. The initiating factors for the bursts are unclear, although retrotransposon mobilization may contribute. The examined cell lines constitute a resource of live experimental models of mutational processes, which potentially retain patterns of activity and regulation operative in primary human cancers.

An integrated genomic analysis of anaplastic meningioma identifies prognostic molecular signatures.

Anaplastic meningioma is a rare and aggressive brain tumor characterised by intractable recurrences and dismal outcomes. Here, we present an integrated analysis of the whole genome, transcriptome and methylation profiles of primary and recurrent anaplastic meningioma. A key finding was the delineation of distinct molecular subgroups that were associated with diametrically opposed survival outcomes. Relative to lower grade meningiomas, anaplastic tumors harbored frequent driver mutations in SWI/SNF complex genes, which were confined to the poor prognosis subgroup. Aggressive disease was further characterised by transcriptional evidence of increased PRC2 activity, stemness and epithelial-to-mesenchymal transition. Our analyses discern biologically distinct variants of anaplastic meningioma with prognostic and therapeutic significance.

The Origins and Vulnerabilities of Two Transmissible Cancers in Tasmanian Devils.

Transmissible cancers are clonal lineages that spread through populations via contagious cancer cells. Although rare in nature, two facial tumor clones affect Tasmanian devils. Here we perform comparative genetic and functional characterization of these lineages. The two cancers have similar patterns of mutation and show no evidence of exposure to exogenous mutagens or viruses. Genes encoding PDGF receptors have copy number gains and are present on extrachromosomal double minutes. Drug screening indicates causative roles for receptor tyrosine kinases and sensitivity to inhibitors of DNA repair. Y chromosome loss from a male clone infecting a female host suggests immunoediting. These results imply that Tasmanian devils may have inherent susceptibility to transmissible cancers and present a suite of therapeutic compounds for use in conservation.