Chemotherapy induces myeloid-driven spatial T-cell exhaustion in ovarian cancer.

To uncover the intricate, chemotherapy-induced spatiotemporal remodeling of the tumor microenvironment, we conducted integrative spatial and molecular characterization of 97 high-grade serous ovarian cancer (HGSC) samples collected before and after chemotherapy. Using single-cell and spatial analyses, we identify increasingly versatile immune cell states, which form spatiotemporally dynamic microcommunities at the tumor-stroma interface. We demonstrate that chemotherapy triggers spatial redistribution and exhaustion of CD8+ T cells due to prolonged antigen presentation by macrophages, both within interconnected myeloid networks termed "Myelonets" and at the tumor stroma interface. Single-cell and spatial transcriptomics identifies prominent TIGIT-NECTIN2 ligand-receptor interactions induced by chemotherapy. Using a functional patient-derived immuno-oncology platform, we show that CD8+T-cell activity can be boosted by combining immune checkpoint blockade with chemotherapy. Our discovery of chemotherapy-induced myeloid-driven spatial T-cell exhaustion paves the way for novel immunotherapeutic strategies to unleash CD8+ T-cell-mediated anti-tumor immunity in HGSC.

Tracing back primed resistance in cancer via sister cells.

Exploring non-genetic evolution of cell states during cancer treatments has become attainable by recent advances in lineage-tracing methods. However, transcriptional changes that drive cells into resistant fates may be subtle, necessitating high resolution analysis. Here, we present ReSisTrace that uses shared transcriptomic features of sister cells to predict the states priming treatment resistance. Applying ReSisTrace in ovarian cancer cells perturbed with olaparib, carboplatin or natural killer (NK) cells reveals pre-resistant phenotypes defined by proteostatic and mRNA surveillance features, reflecting traits enriched in the upcoming subclonal selection. Furthermore, we show that DNA repair deficiency renders cells susceptible to both DNA damaging agents and NK killing in a context-dependent manner. Finally, we leverage the obtained pre-resistance profiles to predict and validate small molecules driving cells to sensitive states prior to treatment. In summary, ReSisTrace resolves pre-existing transcriptional features of treatment vulnerability, facilitating both molecular patient stratification and discovery of synergistic pre-sensitizing therapies.

Functional Homologous Recombination Assay on FFPE Specimens of Advanced High-Grade Serous Ovarian Cancer Predicts Clinical Outcomes.

PURPOSE: Deficiency in homologous recombination (HR) repair of DNA damage is characteristic of many high-grade serous ovarian cancers (HGSC). It is imperative to identify patients with homologous recombination-deficient (HRD) tumors as they are most likely to benefit from platinum-based chemotherapy and PARP inhibitors (PARPi). Existing methods measure historical, not necessarily current HRD and/or require high tumor cell content, which is not achievable for many patients. We set out to develop a clinically feasible assay for identifying functionally HRD tumors that can predict clinical outcomes. EXPERIMENTAL DESIGN: We quantified RAD51, a key HR protein, in immunostained formalin-fixed, paraffin-embedded (FFPE) tumor samples obtained from chemotherapy-naïve and neoadjuvant chemotherapy (NACT)-treated HGSC patients. We defined cutoffs for functional HRD separately for these sample types, classified the patients accordingly as HRD or HR-proficient, and analyzed correlations with clinical outcomes. From the same specimens, genomics-based HRD estimates (HR gene mutations, genomic signatures, and genomic scars) were also determined, and compared with functional HR (fHR) status. RESULTS: fHR status significantly predicted several clinical outcomes, including progression-free survival (PFS) and overall survival (OS), when determined from chemo-naïve (PFS, P < 0.0001; OS, P < 0.0001) as well as NACT-treated (PFS, P < 0.0001; OS, P = 0.0033) tumor specimens. The fHR test also identified as HRD those PARPi-at-recurrence-treated patients with longer OS (P = 0.0188). CONCLUSIONS: We developed an fHR assay performed on routine FFPE specimens, obtained from either chemo-naïve or NACT-treated HGSC patients, that can significantly predict real-world platinum-based chemotherapy and PARPi response. See related commentary by Garg and Oza, p. 2957.

Optimized detection of homologous recombination deficiency improves the prediction of clinical outcomes in cancer.

Homologous recombination DNA-repair deficiency (HRD) is a common driver of genomic instability and confers a therapeutic vulnerability in cancer. The accurate detection of somatic allelic imbalances (AIs) has been limited by methods focused on BRCA1/2 mutations and using mixtures of cancer types. Using pan-cancer data, we revealed distinct patterns of AIs in high-grade serous ovarian cancer (HGSC). We used machine learning and statistics to generate improved criteria to identify HRD in HGSC (ovaHRDscar). ovaHRDscar significantly predicted clinical outcomes in three independent patient cohorts with higher precision than previous methods. Characterization of 98 spatiotemporally distinct metastatic samples revealed low intra-patient variation and indicated the primary tumor as the preferred site for clinical sampling in HGSC. Further, our approach improved the prediction of clinical outcomes in triple-negative breast cancer (tnbcHRDscar), validated in two independent patient cohorts. In conclusion, our tumor-specific, systematic approach has the potential to improve patient selection for HR-targeted therapies.

Correction: Whole genome variation in 27 Mexican indigenous populations, demographic and biomedical insights.

[This corrects the article DOI: 10.1371/journal.pone.0249773.].

Fanconi Anemia Patients from an Indigenous Community in Mexico Carry a New Founder Pathogenic Variant in FANCG.

Fanconi anemia (FA) is a rare genetic disorder caused by pathogenic variants (PV) in at least 22 genes, which cooperate in the Fanconi anemia/Breast Cancer (FA/BRCA) pathway to maintain genome stability. PV in FANCA, FANCC, and FANCG account for most cases (~90%). This study evaluated the chromosomal, molecular, and physical phenotypic findings of a novel founder FANCG PV, identified in three patients with FA from the Mixe community of Oaxaca, Mexico. All patients presented chromosomal instability and a homozygous PV, FANCG: c.511-3_511-2delCA, identified by next-generation sequencing analysis. Bioinformatic predictions suggest that this deletion disrupts a splice acceptor site promoting the exon 5 skipping. Analysis of Cytoscan 750 K arrays for haplotyping and global ancestry supported the Mexican origin and founder effect of the variant, reaffirming the high frequency of founder PV in FANCG. The degree of bone marrow failure and physical findings (described through the acronyms VACTERL-H and PHENOS) were used to depict the phenotype of the patients. Despite having a similar frequency of chromosomal aberrations and genetic constitution, the phenotype showed a wide spectrum of severity. The identification of a founder PV could help for a systematic and accurate genetic screening of patients with FA suspicion in this population.

Whole genome variation in 27 Mexican indigenous populations, demographic and biomedical insights.

There has been limited study of Native American whole genome diversity to date, which impairs effective implementation of personalized medicine and a detailed description of its demographic history. Here we report high coverage whole genome sequencing of 76 unrelated individuals, from 27 indigenous groups across Mexico, with more than 97% average Native American ancestry. On average, each individual has 3.26 million Single Nucleotide Variants and short indels, that together comprise a catalog of 9,737,152 variants, 44,118 of which are novel. We report 497 common Single Nucleotide Variants (with allele frequency > 5%) mapped to drug responses and 316,577 in enhancer or promoter elements; interestingly we found some of these enhancer variants in PPARG, a nuclear receptor involved in highly prevalent health problems in Mexican population, such as obesity, diabetes, and insulin resistance. By detecting signals of positive selection we report 24 enriched key pathways under selection, most of them related to immune mechanisms. No missense variants in ACE2, the receptor responsible for the entry of the SARS CoV-2 virus, were found in any individual. Population genomics and phylogenetic analyses demonstrated stratification in a Northern-Central-Southern axis, with major substructure in the Central region. The Seri, a northern group with the most genetic divergence in our study, showed a distinctive genomic context with the most novel variants, and the most population specific genotypes. Genome-wide analysis showed that the average haplotype blocks are longer in Native Mexicans than in other world populations. With this dataset we describe previously undetected population level variation in Native Mexicans, helping to reduce the gap in genomic data representation of such groups.

Heterogeneity and Clonal Evolution of Acquired PARP Inhibitor Resistance in TP53- and BRCA1-Deficient Cells.

Homologous recombination (HR)-deficient cancers are sensitive to poly-ADP ribose polymerase inhibitors (PARPi), which have shown clinical efficacy in the treatment of high-grade serous cancers (HGSC). However, the majority of patients will relapse, and acquired PARPi resistance is emerging as a pressing clinical problem. Here we generated seven single-cell clones with acquired PARPi resistance derived from a PARPi-sensitive TP53 -/- and BRCA1 -/- epithelial cell line generated using CRISPR/Cas9. These clones showed diverse resistance mechanisms, and some clones presented with multiple mechanisms of resistance at the same time. Genomic analysis of the clones revealed unique transcriptional and mutational profiles and increased genomic instability in comparison with a PARPi-sensitive cell line. Clonal evolutionary analyses suggested that acquired PARPi resistance arose via clonal selection from an intrinsically unstable and heterogenous cell population in the sensitive cell line, which contained preexisting drug-tolerant cells. Similarly, clonal and spatial heterogeneity in tumor biopsies from a clinical patient with BRCA1-mutant HGSC with acquired PARPi resistance was observed. In an imaging-based drug screening, the clones showed heterogenous responses to targeted therapeutic agents, indicating that not all PARPi-resistant clones can be targeted with just one therapy. Furthermore, PARPi-resistant clones showed mechanism-dependent vulnerabilities to the selected agents, demonstrating that a deeper understanding on the mechanisms of resistance could lead to improved targeting and biomarkers for HGSC with acquired PARPi resistance. SIGNIFICANCE: This study shows that BRCA1-deficient cells can give rise to multiple genomically and functionally heterogenous PARPi-resistant clones, which are associated with various vulnerabilities that can be targeted in a mechanism-specific manner.

Genetically Defined Syngeneic Mouse Models of Ovarian Cancer as Tools for the Discovery of Combination Immunotherapy.

Despite advances in immuno-oncology, the relationship between tumor genotypes and response to immunotherapy remains poorly understood, particularly in high-grade serous tubo-ovarian carcinomas (HGSC). We developed a series of mouse models that carry genotypes of human HGSCs and grow in syngeneic immunocompetent hosts to address this gap. We transformed murine-fallopian tube epithelial cells to phenocopy homologous recombination-deficient tumors through a combined loss of Trp53, Brca1, Pten, and Nf1 and overexpression of Myc and Trp53 R172H, which was contrasted with an identical model carrying wild-type Brca1. For homologous recombination-proficient tumors, we constructed genotypes combining loss of Trp53 and overexpression of Ccne1, Akt2, and Trp53 R172H, and driven by KRAS G12V or Brd4 or Smarca4 overexpression. These lines form tumors recapitulating human disease, including genotype-driven responses to treatment, and enabled us to identify follistatin as a driver of resistance to checkpoint inhibitors. These data provide proof of concept that our models can identify new immunotherapy targets in HGSC. SIGNIFICANCE: We engineered a panel of murine fallopian tube epithelial cells bearing mutations typical of HGSC and capable of forming tumors in syngeneic immunocompetent hosts. These models recapitulate tumor microenvironments and drug responses characteristic of human disease. In a Ccne1-overexpressing model, immune-checkpoint resistance was driven by follistatin.This article is highlighted in the In This Issue feature, p. 211.

Delivery mode-associated gut microbiota in the first 3 months of life in a country with high obesity rates: A descriptive study.

Delivery methods during childbirth and their related gut microbiota profiles have important impacts on health later in life, they can contribute to the development of diseases such as obesity, whose highest prevalence rate is found among the Mexican child population. Coincidentally, Mexico has one of the highest global average annual rate increase in cesarean births (C-section). Since Mexico leads the world in childhood obesity, studying the relationship between childbirth delivery methods and gut microbiota profiles in this vulnerable population may be used to identify early risk factors for obesity in other developed and developing countries. The objective of this study is to determine the association between child delivery method and gut microbiota profiles in healthy Mexican newborns.Fecal samples of 57 term infants who participated in a randomized clinical trial in 2013 to study the safety of Agave fructans in newborns, were used in this study. DNA samples were extracted and used to characterize the microbiota composition using high-throughput 16S rRNA gene sequencing. The samples were further divided based on childbirth delivery method, as well as early diet. Gut microbiota profiles were determined and analyzed using cluster analysis followed by multiple correspondence analysis.An unusual high abundance of Proteobacteria was found in the gut microbiota of all Mexican infants studied, regardless of delivery method. Feces from infants born by C-section had low levels of Bacteroidetes, high levels of Firmicutes, especially Clostridium and Enterococcus, and a strikingly high ratio of Firmicutes/Bacteroidetes (F:B). Profiles enriched in Bacteroidetes and low F:B ratios, were strongly associated with vaginal delivery.The profile of gut microbiota associated with feces from Mexican infants born by C-section, may be added to the list of boosting factors for the worrying obesity epidemic in Mexico.

Genomic study of dilated cardiomyopathy in a group of Mexican patients using site-directed next generation sequencing.

BACKGROUND: Dilated cardiomyopathy (DCM) is a major cause of nonischemic heart failure and death in young adults. Next generation sequencing (NGS) has become part of the diagnostic workup in idiopathic and familial DCM. More than 50 DCM genes have been identified, revealing great molecular heterogeneity and variable diagnostic yield. Interpretation of variant pathogenicity is challenging particularly in underrepresented populations, as pathogenic variant databases include studies mainly from European/Caucasian populations. To date, no studies on genomic diagnosis of DCM have been conducted in Mexico. METHODS: We recruited 55 unrelated DCM patients, 22 familial (F-DCM), and 33 idiopathic (I-DCM), and performed site-directed NGS seeking causal mutations. Diagnostic yield was defined as the proportion of individuals with at least one pathogenic (P) or likely pathogenic (LP) variant in DCM genes. RESULTS: Overall diagnostic yield was 47.3%, and higher in F-DCM (63.6%) than in I-DCM (36.4%, p = 0.047). Overall, NGS disclosed 41 variants of clinical interest (61.0% novel), 27 were classified as P/LP and 14 of unknown clinical significance. Of P/LP variants, 10 were A-band region TTN truncating variants, five were found in DSP (18.5%), five in MYH7 (18.5%), two in LMNA (7.4%), and one in RBM20, ABCC9, FKTN, ACTA1, and TNNT2. NGS findings suggested autosomal recessive inheritance in three families, two with DSP loss of function mutations in affected individuals. The increasing number of mutation reports in DCM, increasing knowledge on the functional consequences of mutations, mutational hotspots and functional domains of DCM-related proteins, the recent refinement ACMG/ClinGen Guidelines, and co-segregation analysis in DCM families helped increase the diagnostic yield. CONCLUSION: This is the first NGS study performed in a group of Mexican DCM patients, contributing to understand the mutational spectrum and complexity of DCM molecular diagnosis.

Marine Archaeon Methanosarcina acetivorans Enhances Polyphosphate Metabolism Under Persistent Cadmium Stress.

Phosphate metabolism was studied to determine whether polyphosphate (polyP) pools play a role in the enhanced resistance against Cd2+ and metal-removal capacity of Cd2+-preadapted (CdPA) Methanosarcina acetivorans. Polyphosphate kinase (PPK), exopolyphosphatase (PPX) and phosphate transporter transcript levels and their activities increased in CdPA cells compared to control (Cnt) cells. K+ inhibited recombinant Ma-PPK and activated Ma-PPX, whereas divalent cations activated both enzymes. Metal-binding polyP and thiol-containing molecule contents, Cd2+-removal, and biofilm synthesis were significantly higher in CdPA cells >Cnt cells plus a single addition of Cd2+>Cnt cells. Also, CdPA cells showed a higher number of cadmium, sulfur, and phosphorus enriched-acidocalcisomes than control cells. Biochemical and physiological phenotype exhibited by CdPA cells returned to that of Cnt cells when cultured without Cd2+. Furthermore, no differences in the sequenced genomes upstream and downstream of the genes involved in Cd2+ resistance were found between CdPA and Cnt cells, suggesting phenotype loss rather than genome mutations induced by chronic Cd2+-exposure. Instead, a metabolic adaptation induced by Cd2+ stress was apparent. The dynamic ability of M. acetivorans to change its metabolism, depending on the environmental conditions, may be advantageous to remove cadmium in nature and biodigesters.

Demographic history and biologically relevant genetic variation of Native Mexicans inferred from whole-genome sequencing.

Understanding the genetic structure of Native American populations is important to clarify their diversity, demographic history, and to identify genetic factors relevant for biomedical traits. Here, we show a demographic history reconstruction from 12 Native American whole genomes belonging to six distinct ethnic groups representing the three main described genetic clusters of Mexico (Northern, Southern, and Maya). Effective population size estimates of all Native American groups remained below 2,000 individuals for up to 10,000 years ago. The proportion of missense variants predicted as damaging is higher for undescribed (~ 30%) than for previously reported variants (~ 15%). Several variants previously associated with biological traits are highly frequent in the Native American genomes. These findings suggest that the demographic and adaptive processes that occurred in these groups shaped their genetic architecture and could have implications in biological processes of the Native Americans and Mestizos of today.