CXCL13 Positive Cells Localization Predict Response to Anti-PD-1/PD-L1 in Pulmonary Non-Small Cell Carcinoma

Background: Immune checkpoint inhibitors (ICIs) have revolutionized non-small cell lung cancers (NSCLCs) treatment, but only 20-30% of patients benefit from these treatments. Currently, PD-L1 expression in tumor cells is the only clinically approved predictor of ICI response in lung cancer, but concerns arise due to its low negative and positive predictive value. Recent studies suggest that CXCL13+ T cells in the tumor microenvironment (TME) may be a good predictor of response. We aimed to assess if CXCL13+ cell localization within the TME can predict ICI response in advanced NSCLC patients. Methods: This retrospective study included 65 advanced NSCLC patients treated with Nivolumab/Pembrolizumab at IUCPQ or CHUM and for whom a pretreatment surgical specimen was available. Good responders were defined as having a complete radiologic response at 1 year, and bad responders were defined as showing cancer progression at 1 year. IHC staining for CXCL13 was carried out on a representative slide from a resection specimen, and CXCL13+ cell density was evaluated in tumor (T), invasive margin (IM), non-tumor (NT), and tertiary lymphoid structure (TLS) compartments. Cox models were used to analyze progression-free survival (PFS) and overall survival (OS) probability, while the Mann-Whitney test was used to compare CXCL13+ cell density between responders and non-responders. Results: We showed that CXCL13+ cell density localization within the TME is associated with ICI efficacy. An increased density of CXCL13+ cells across all compartments was associated with a poorer prognostic (OS; HR = 1.22; 95%CI = 1.04-1.42; p = 0.01, PFS; HR = 1.16; p = 0.02), or a better prognostic when colocalized within TLSs (PFS; HR = 0.84, p = 0.03). Conclusion: Our results support the role of CXCL13+ cells in advanced NSCLC patients, with favorable prognosis when localized within TLSs and unfavorable prognosis when present elsewhere. The concomitant proximity of CXCL13+ and CD20+ cells within TLSs may favor antigen presentation to T cells, thus enhancing the effect of PD-1/PD-L1 axis inhibition. Further validation is warranted to confirm the potential relevance of this biomarker in a clinical setting.

The Cross-Scale Association between Pathomics and Radiomics Features in Immunotherapy-Treated NSCLC Patients: A Preliminary Study

BACKGROUND

Recent advances in cancer biomarker development have led to a surge of distinct data modalities, such as medical imaging and histopathology. To develop predictive immunotherapy biomarkers, these modalities are leveraged independently, despite their orthogonality. This study aims to explore the cross-scale association between radiological scans and digitalized pathology images for immunotherapy-treated non-small cell lung cancer (NSCLC) patients.

METHODS

This study involves 36 NSCLC patients who were treated with immunotherapy and for whom both radiology and pathology images were available. A total of 851 and 260 features were extracted from CT scans and cell density maps of histology images at different resolutions. We investigated the radiopathomics relationship and their association with clinical and biological endpoints. We used the Kolmogorov-Smirnov (KS) method to test the differences between the distributions of correlation coefficients with the two imaging modality features. Unsupervised clustering was done to identify which imaging modality captures poor and good survival patients.

RESULTS

Our results demonstrated a significant correlation between cell density pathomics and radiomics features. Furthermore, we also found a varying distribution of correlation values between imaging-derived features and clinical endpoints. The KS test revealed that the two imaging feature distributions were different for PFS and CD8 counts, while similar for OS. In addition, clustering analysis resulted in significant differences in the two clusters generated from the radiomics and pathomics features with respect to patient survival and CD8 counts.

CONCLUSION

The results of this study suggest a cross-scale association between CT scans and pathology H&E slides among ICI-treated patients. These relationships can be further explored to develop multimodal immunotherapy biomarkers to advance personalized lung cancer care.

A Radiomics-Clinical Model Predicts Overall Survival of Non-Small Cell Lung Cancer Patients Treated with Immunotherapy: A Multicenter Study

BACKGROUND

Immune checkpoint inhibitors (ICIs) are a great breakthrough in cancer treatments and provide improved long-term survival in a subset of non-small cell lung cancer (NSCLC) patients. However, prognostic and predictive biomarkers of immunotherapy still remain an unmet clinical need. In this work, we aim to leverage imaging data and clinical variables to develop survival risk models among advanced NSCLC patients treated with immunotherapy.

METHODS

This retrospective study includes a total of 385 patients from two institutions who were treated with ICIs. Radiomics features extracted from pretreatment CT scans were used to build predictive models. The objectives were to predict overall survival (OS) along with building a classifier for short- and long-term survival groups. We employed the XGBoost learning method to build radiomics and integrated clinical-radiomics predictive models. Feature selection and model building were developed and validated on a multicenter cohort.

RESULTS

We developed parsimonious models that were associated with OS and a classifier for short- and long-term survivor groups. The concordance indices (C-index) of the radiomics model were 0.61 and 0.57 to predict OS in the discovery and validation cohorts, respectively. While the area under the curve (AUC) values of the radiomic models for short- and long-term groups were found to be 0.65 and 0.58 in the discovery and validation cohorts. The accuracy of the combined radiomics-clinical model resulted in 0.63 and 0.62 to predict OS and in 0.77 and 0.62 to classify the survival groups in the discovery and validation cohorts, respectively.

CONCLUSIONS

We developed and validated novel radiomics and integrated radiomics-clinical survival models among NSCLC patients treated with ICIs. This model has important translational implications, which can be used to identify a subset of patients who are not likely to benefit from immunotherapy. The developed imaging biomarkers may allow early prediction of low-group survivors, though additional validation of these radiomics models is warranted.

Recombination affects allele-specific expression of deleterious variants in human populations

How the genetic composition of a population changes through stochastic processes, such as genetic drift, in combination with deterministic processes, such as selection, is critical to understanding how phenotypes vary in space and time. Here, we show how evolutionary forces affecting selection, including recombination and effective population size, drive genomic patterns of allele-specific expression (ASE). Integrating tissue-specific genotypic and transcriptomic data from 1500 individuals from two different cohorts, we demonstrate that ASE is less often observed in regions of low recombination, and loci in high or normal recombination regions are more efficient at using ASE to underexpress harmful mutations. By tracking genetic ancestry, we discriminate between ASE variability due to past demographic effects, including subsequent bottlenecks, versus local environment. We observe that ASE is not randomly distributed along the genome and that population parameters influencing the efficacy of natural selection alter ASE levels genome wide.

Pan-cancer analysis of whole genomes

Cancer is driven by genetic change, and the advent of massively parallel sequencing has enabled systematic documentation of this variation at the whole-genome scale1-3. Here we report the integrative analysis of 2,658 whole-cancer genomes and their matching normal tissues across 38 tumour types from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). We describe the generation of the PCAWG resource, facilitated by international data sharing using compute clouds. On average, cancer genomes contained 4-5 driver mutations when combining coding and non-coding genomic elements; however, in around 5% of cases no drivers were identified, suggesting that cancer driver discovery is not yet complete. Chromothripsis, in which many clustered structural variants arise in a single catastrophic event, is frequently an early event in tumour evolution; in acral melanoma, for example, these events precede most somatic point mutations and affect several cancer-associated genes simultaneously. Cancers with abnormal telomere maintenance often originate from tissues with low replicative activity and show several mechanisms of preventing telomere attrition to critical levels. Common and rare germline variants affect patterns of somatic mutation, including point mutations, structural variants and somatic retrotransposition. A collection of papers from the PCAWG Consortium describes non-coding mutations that drive cancer beyond those in the TERT promoter4; identifies new signatures of mutational processes that cause base substitutions, small insertions and deletions and structural variation5,6; analyses timings and patterns of tumour evolution7; describes the diverse transcriptional consequences of somatic mutation on splicing, expression levels, fusion genes and promoter activity8,9; and evaluates a range of more-specialized features of cancer genomes8,10-18.

Aberrant PRDM9 expression impacts the pan-cancer genomic landscape

The binding of PRDM9 to chromatin is a key step in the induction of DNA double-strand breaks associated with meiotic recombination hotspots; it is normally expressed solely in germ cells. We interrogated 1879 cancer samples in 39 different cancer types and found that PRDM9 is unexpectedly expressed in 20% of these tumors even after stringent gene homology correction. The expression levels of PRDM9 in tumors are significantly higher than those found in healthy neighboring tissues and in healthy nongerm tissue databases. Recurrently mutated regions located within 5 Mb of the PRDM9 loci, as well as differentially expressed genes in meiotic pathways, correlate with PRDM9 expression. In samples with aberrant PRDM9 expression, structural variant breakpoints frequently neighbor the DNA motif recognized by PRDM9, and there is an enrichment of structural variants at sites of known meiotic PRDM9 activity. This study is the first to provide evidence of an association between aberrant expression of the meiosis-specific gene PRDM9 with genomic instability in cancer.

FOXA and master transcription factors recruit Mediator and Cohesin to the core transcriptional regulatory circuitry of cancer cells

Controlling the transcriptional program is essential to maintain the identity and the biological functions of a cell. The Mediator and Cohesin complexes have been established as central cofactors controlling the transcriptional program in normal cells. However, the distribution, recruitment and importance of these complexes in cancer cells have not been fully investigated. Here we show that FOXA and master transcription factors are part of the core transcriptional regulatory circuitry of cancer cells and are essential to recruit M ediator and Cohesin. Indeed, Mediator and Cohesin occupied the enhancer and promoter regions of actively transcribed genes and maintained the proliferation and colony forming potential. Through integration of publically available ChIP-Seq datasets, we predicted the core transcriptional regulatory circuitry of each cancer cell. Unexpectedly, for all cells investigated, the pioneer transcription factors FOXA1 and/or FOXA2 were identified in addition to cell-specific master transcription factors. Loss of both types of transcription factors phenocopied the loss of Mediator and Cohesin. Lastly, the master and pioneer transcription factors were essential to recruit Mediator and Cohesin to regulatory regions of actively transcribed genes. Our study proposes that maintenance of the cancer cell state is dependent on recruitment of Mediator and Cohesin through FOXA and master transcription factors.

Standing chromosomal variation in Lake Whitefish species pairs: the role of historical contingency and relevance for speciation

The role of chromosome changes in speciation remains a debated topic, although demographic conditions associated with divergence should promote their appearance. We tested a potential relationship between chromosome changes and speciation by studying two Lake Whitefish (Coregonus clupeaformis) lineages that recently colonized postglacial lakes following allopatry. A dwarf limnetic species evolved repeatedly from the normal benthic species, becoming reproductively isolated. Lake Whitefish hybrids experience mitotic and meiotic instability, which may result from structurally divergent chromosomes. Motivated by this observation, we test the hypothesis that chromosome organization differs between Lake Whitefish species pairs using cytogenetics. While chromosome and fundamental numbers are conserved between the species (2n = 80, NF = 98), we observe extensive polymorphism of subtle karyotype traits. We describe intrachromosomal differences associated with heterochromatin and repetitive DNA, and test for parallelism among three sympatric species pairs. Multivariate analyses support the hypothesis that differentiation at the level of subchromosomal markers mostly appeared during allopatry. Yet we find no evidence for parallelism between species pairs among lakes, consistent with colonization effect or postcolonization differentiation. The reported intrachromosomal polymorphisms do not appear to play a central role in driving adaptive divergence between normal and dwarf Lake Whitefish. We discuss how chromosomal differentiation in the Lake Whitefish system may contribute to the destabilization of mitotic and meiotic chromosome segregation in hybrids, as documented previously. The chromosome structures detected here are still difficult to sequence and assemble, demonstrating the value of cytogenetics as a complementary approach to understand the genomic bases of speciation.

metagene Profiles Analyses Reveal Regulatory Element's Factor-Specific Recruitment Patterns

ChIP-Sequencing (ChIP-Seq) provides a vast amount of information regarding the localization of proteins across the genome. The aggregation of ChIP-Seq enrichment signal in a metagene plot is an approach commonly used to summarize data complexity and to obtain a high level visual representation of the general occupancy pattern of a protein. Here we present the R package metagene, the graphical interface Imetagene and the companion package similaRpeak. Together, they provide a framework to integrate, summarize and compare the ChIP-Seq enrichment signal from complex experimental designs. Those packages identify and quantify similarities or dissimilarities in patterns between large numbers of ChIP-Seq profiles. We used metagene to investigate the differential occupancy of regulatory factors at noncoding regulatory regions (promoters and enhancers) in relation to transcriptional activity in GM12878 B-lymphocytes. The relationships between occupancy patterns and transcriptional activity suggest two different mechanisms of action for transcriptional control: i) a “gradient effect” where the regulatory factor occupancy levels follow transcription and ii) a “threshold effect” where the regulatory factor occupancy levels max out prior to reaching maximal transcription. metagene, Imetagene and similaRpeak are implemented in R under the Artistic license 2.0 and are available on Bioconductor.

Mutant cohesin affects RNA polymerase II regulation in Cornelia de Lange syndrome

In addition to its role in sister chromatid cohesion, genome stability and integrity, the cohesin complex is involved in gene transcription. Mutations in core cohesin subunits SMC1A, SMC3 and RAD21, or their regulators NIPBL and HDAC8, cause Cornelia de Lange syndrome (CdLS). Recent evidence reveals that gene expression dysregulation could be the underlying mechanism for CdLS. These findings raise intriguing questions regarding the potential role of cohesin-mediated transcriptional control and pathogenesis. Here, we identified numerous dysregulated genes occupied by cohesin by combining the transcriptome of CdLS cell lines carrying mutations in SMC1A gene and ChIP-Seq data. Genome-wide analyses show that genes changing in expression are enriched for cohesin-binding. In addition, our results indicate that mutant cohesin impairs both RNA polymerase II (Pol II) transcription initiation at promoters and elongation in the gene body. These findings highlight the pivotal role of cohesin in transcriptional regulation and provide an explanation for the typical gene dysregulation observed in CdLS patients.

Neutral and selective processes shape MHC gene diversity and expression in stocked brook charr populations (Salvelinus fontinalis)

The capacity of an individual to battle infection is an important fitness determinant in wild vertebrate populations. The major histocompatibility complex (MHC) genes are crucial for a host's adaptive immune system to detect pathogens. However, anthropogenic activities may disrupt natural cycles of co-evolution between hosts and pathogens. In this study, we investigated the dynamic sequence and expression variation of host parasite interactions in brook charr (Salvelinus fontinalis) in a context of past human disturbance via population supplementation from domestic individuals. To do so, we developed a new method to examine selection shaping MHC diversity within and between populations and found a complex interplay between neutral and selective processes that varied between lakes that were investigated. We provided evidence for a lower introgression rate of domestic alleles and found that parasite infection increased with domestic genomic background of individuals. We also documented an association between individual MHC alleles and parasite taxa. Finally, longer cis-regulatory minisatellites were positively correlated with MHC II down-regulation and domestic admixture, suggesting that inadvertent selection during domestication resulted in a lower immune response capacity, through a trade-off between growth and immunity, which explained the negative selection of domestic alleles at least under certain circumstances.

Stocking impacts the expression of candidate genes and physiological condition in introgressed brook charr (Salvelinus fontinalis) populations

Translocation of plants and animal populations between environments is one of the major forms of anthropogenic perturbation experienced by pristine populations, and consequently, human-mediated hybridization by stocking practices between wild and exogenous conspecifics is of increasing concern. In this study, we compared the expression of seven candidate genes involved in multifactorial traits and regulatory pathways for growth as a function of level of introgressive hybridization between wild and domestic brook charr to test the null hypothesis of no effect of introgression on wild fish. Our analyses revealed that the expression of two of the genes tested, cytochrome c oxidase VIIa and the growth hormone receptor isoform I, was positively correlated with the level of introgression. We also observed a positive relationship between the extent of introgression and physiological status quantified by the Fulton's condition index. The expression of other genes was influenced by other variables, including year of sampling (reflecting different thermal conditions), sampling method and lake of origin. This is the first demonstration in nature that introgression from stocked populations has an impact on the expression of genes playing a role in important biological functions that may be related with fitness in wild introgressed populations.