Tissue-resident FOLR2+ macrophages associate with CD8+ T cell infiltration in human breast cancer.

Macrophage infiltration is a hallmark of solid cancers, and overall macrophage infiltration correlates with lower patient survival and resistance to therapy. Tumor-associated macrophages, however, are phenotypically and functionally heterogeneous. Specific subsets of tumor-associated macrophage might be endowed with distinct roles on cancer progression and antitumor immunity. Here, we identify a discrete population of FOLR2+ tissue-resident macrophages in healthy mammary gland and breast cancer primary tumors. FOLR2+ macrophages localize in perivascular areas in the tumor stroma, where they interact with CD8+ T cells. FOLR2+ macrophages efficiently prime effector CD8+ T cells ex vivo. The density of FOLR2+ macrophages in tumors positively correlates with better patient survival. This study highlights specific roles for tumor-associated macrophage subsets and paves the way for subset-targeted therapeutic interventions in macrophages-based cancer therapies.

Macrophages Maintain Epithelium Integrity by Limiting Fungal Product Absorption.

The colon is primarily responsible for absorbing fluids. It contains a large number of microorganisms including fungi, which are enriched in its distal segment. The colonic mucosa must therefore tightly regulate fluid influx to control absorption of fungal metabolites, which can be toxic to epithelial cells and lead to barrier dysfunction. How this is achieved remains unknown. Here, we describe a mechanism by which the innate immune system allows rapid quality check of absorbed fluids to avoid intoxication of colonocytes. This mechanism relies on a population of distal colon macrophages that are equipped with "balloon-like" protrusions (BLPs) inserted in the epithelium, which sample absorbed fluids. In the absence of macrophages or BLPs, epithelial cells keep absorbing fluids containing fungal products, leading to their death and subsequent loss of epithelial barrier integrity. These results reveal an unexpected and essential role of macrophages in the maintenance of colon-microbiota interactions in homeostasis. VIDEO ABSTRACT.

Epithelial colonization by gut dendritic cells promotes their functional diversification.

Dendritic cells (DCs) patrol tissues and transport antigens to lymph nodes to initiate adaptive immune responses. Within tissues, DCs constitute a complex cell population composed of distinct subsets that can exhibit different activation states and functions. How tissue-specific cues orchestrate DC diversification remains elusive. Here, we show that the small intestine included two pools of cDC2s originating from common pre-DC precursors: (1) lamina propria (LP) CD103+CD11b+ cDC2s that were mature-like proinflammatory cells and (2) intraepithelial cDC2s that exhibited an immature-like phenotype as well as tolerogenic properties. These phenotypes resulted from the action of food-derived retinoic acid (ATRA), which enhanced actomyosin contractility and promoted LP cDC2 transmigration into the epithelium. There, cDC2s were imprinted by environmental cues, including ATRA itself and the mucus component Muc2. Hence, by reaching distinct subtissular niches, DCs can exist as immature and mature cells within the same tissue, revealing an additional mechanism of DC functional diversification.

A Conserved Noncoding Locus Regulates Random Monoallelic Xist Expression across a Topological Boundary.

cis-Regulatory communication is crucial in mammalian development and is thought to be restricted by the spatial partitioning of the genome in topologically associating domains (TADs). Here, we discovered that the Xist locus is regulated by sequences in the neighboring TAD. In particular, the promoter of the noncoding RNA Linx (LinxP) acts as a long-range silencer and influences the choice of X chromosome to be inactivated. This is independent of Linx transcription and independent of any effect on Tsix, the antisense regulator of Xist that shares the same TAD as Linx. Unlike Tsix, LinxP is well conserved across mammals, suggesting an ancestral mechanism for random monoallelic Xist regulation. When introduced in the same TAD as Xist, LinxP switches from a silencer to an enhancer. Our study uncovers an unsuspected regulatory axis for X chromosome inactivation and a class of cis-regulatory effects that may exploit TAD partitioning to modulate developmental decisions.

m6A RNA methylation regulates the fate of endogenous retroviruses.

Endogenous retroviruses (ERVs) are abundant and heterogenous groups of integrated retroviral sequences that affect genome regulation and cell physiology throughout their RNA-centred life cycle1. Failure to repress ERVs is associated with cancer, infertility, senescence and neurodegenerative diseases2,3. Here, using an unbiased genome-scale CRISPR knockout screen in mouse embryonic stem cells, we identify m6A RNA methylation as a way to restrict ERVs. Methylation of ERV mRNAs is catalysed by the complex of methyltransferase-like METTL3-METTL144 proteins, and we found that depletion of METTL3-METTL14, along with their accessory subunits WTAP and ZC3H13, led to increased mRNA abundance of intracisternal A-particles (IAPs) and related ERVK elements specifically, by targeting their 5' untranslated region. Using controlled auxin-dependent degradation of the METTL3-METTL14 enzymatic complex, we showed that IAP mRNA and protein abundance is dynamically and inversely correlated with m6A catalysis. By monitoring chromatin states and mRNA stability upon METTL3-METTL14 double depletion, we found that m6A methylation mainly acts by reducing the half-life of IAP mRNA, and this occurs by the recruitment of the YTHDF family of m6A reader proteins5. Together, our results indicate that RNA methylation provides a protective effect in maintaining cellular integrity by clearing reactive ERV-derived RNA species, which may be especially important when transcriptional silencing is less stringent.

Unraveling the features of somatic transposition in the Drosophila intestine.

Transposable elements (TEs) play a significant role in evolution, contributing to genetic variation. However, TE mobilization in somatic cells is not well understood. Here, we address the prevalence of transposition in a somatic tissue, exploiting the Drosophila midgut as a model. Using whole-genome sequencing of in vivo clonally expanded gut tissue, we have mapped hundreds of high-confidence somatic TE integration sites genome-wide. We show that somatic retrotransposon insertions are associated with inactivation of the tumor suppressor Notch, likely contributing to neoplasia formation. Moreover, applying Oxford Nanopore long-read sequencing technology we provide evidence for tissue-specific differences in retrotransposition. Comparing somatic TE insertional activity with transcriptomic and small RNA sequencing data, we demonstrate that transposon mobility cannot be simply predicted by whole tissue TE expression levels or by small RNA pathway activity. Finally, we reveal that somatic TE insertions in the adult fly intestine are enriched in genic regions and in transcriptionally active chromatin. Together, our findings provide clear evidence of ongoing somatic transposition in Drosophila and delineate previously unknown features underlying somatic TE mobility in vivo.

Artificial intelligence-based pathology as a biomarker of sensitivity to atezolizumab-bevacizumab in patients with hepatocellular carcinoma: a multicentre retrospective study.

BACKGROUND: Clinical benefits of atezolizumab plus bevacizumab (atezolizumab-bevacizumab) are observed only in a subset of patients with hepatocellular carcinoma and the development of biomarkers is needed to improve therapeutic strategies. The atezolizumab-bevacizumab response signature (ABRS), assessed by molecular biology profiling techniques, has been shown to be associated with progression-free survival after treatment initiation. The primary objective of our study was to develop an artificial intelligence (AI) model able to estimate ABRS expression directly from histological slides, and to evaluate if model predictions were associated with progression-free survival. METHODS: In this multicentre retrospective study, we developed a model (ABRS-prediction; ABRS-P), which was derived from the previously published clustering-constrained attention multiple instance learning (or CLAM) pipeline. We trained the model fit for regression analysis using a multicentre dataset from The Cancer Genome Atlas (patients treated by surgical resection, n=336). The ABRS-P model was externally validated on two independent series of samples from patients with hepatocellular carcinoma (a surgical resection series, n=225; and a biopsy series, n=157). The predictive value of the model was further tested in a series of biopsy samples from a multicentre cohort of patients with hepatocellular carcinoma treated with atezolizumab-bevacizumab (n=122). All samples in the study were from adults (aged ≥18 years). The validation sets were sampled between Jan 1, 2008, to Jan 1, 2023. For the multicentre validation set, the primary objective was to assess the association of high versus low ABRS-P values, defined relative to cross-validation median split thresholds in the first biopsy series, with progression-free survival after treatment initiation. Finally, we performed spatial transcriptomics and matched prediction heatmaps with in situ expression profiles. FINDINGS: Of the 840 patients sampled, 641 (76%) were male and 199 (24%) were female. Across the development and validation datasets, hepatocellular carcinoma risk factors included alcohol intake, hepatitis B and C virus infections, and non-alcoholic steatohepatitis. Using cross-validation in the development series, the mean Pearson's correlation between ABRS-P values and ABRS score (mean expression of ABRS genes) was r=0·62 (SD 0·09; mean p<0·0001, SD<0·0001). The ABRS-P generalised well on the external validation series (surgical resection series, r=0·60 [95% CI 0·51-0·68], p<0·0001; biopsy series, r=0·53 [0·40-0·63], p<0·0001). In the 122 patients treated with atezolizumab-bevacizumab, those with ABRS-P-high tumours (n=74) showed significantly longer median progression-free survival than those with ABRS-P-low tumours (n=48) after treatment initiation (12 months [95% CI 7-not reached] vs 7 months [4-9]; p=0·014). Spatial transcriptomics showed significantly higher ABRS score, along with upregulation of various other immune effectors, in tumour areas with high ABRS-P values versus areas with low ABRS-P values. INTERPRETATION: Our study indicates that AI applied on hepatocellular carcinoma digital slides is able to serve as a biomarker for progression-free survival in patients treated with atezolizumab-bevacizumab. This approach could be used in the development of inexpensive and fast biomarkers for targeted therapies. The combination of AI heatmaps with spatial transcriptomics provides insight on the molecular features associated with predictions. This methodology could be applied to other cancers or diseases and improve understanding of the biological mechanisms that drive responses to treatments. FUNDING: Institut National du Cancer, Fondation ARC, China Scholarship Council, Ligue Contre le Cancer du Val de Marne, Fondation de l'Avenir, Ipsen, and Fondation Bristol Myers Squibb Pour la Recherche en Immuno-Oncologie.

Human papilloma virus (HPV) integration signature in Cervical Cancer: identification of MACROD2 gene as HPV hot spot integration site.

BACKGROUND: Cervical cancer (CC) remains a leading cause of gynaecological cancer-related mortality with infection by human papilloma virus (HPV) being the most important risk factor. We analysed the association between different viral integration signatures, clinical parameters and outcome in pre-treated CCs. METHODS: Different integration signatures were identified using HPV double capture followed by next-generation sequencing (NGS) in 272 CC patients from the BioRAIDs study [NCT02428842]. Correlations between HPV integration signatures and clinical, biological and molecular features were assessed. RESULTS: Episomal HPV was much less frequent in CC as compared to anal carcinoma (p < 0.0001). We identified >300 different HPV-chromosomal junctions (inter- or intra-genic). The most frequent integration site in CC was in MACROD2 gene followed by MIPOL1/TTC6 and TP63. HPV integration signatures were not associated with histological subtype, FIGO staging, treatment or PFS. HPVs were more frequently episomal in PIK3CA mutated tumours (p = 0.023). Viral integration type was dependent on HPV genotype (p < 0.0001); HPV18 and HPV45 being always integrated. High HPV copy number was associated with longer PFS (p = 0.011). CONCLUSIONS: This is to our knowledge the first study assessing the prognostic value of HPV integration in a prospectively annotated CC cohort, which detects a hotspot of HPV integration at MACROD2; involved in impaired PARP1 activity and chromosome instability.

Functional variants of human papillomavirus type 16 demonstrate host genome integration and transcriptional alterations corresponding to their unique cancer epidemiology.

BACKGROUND: Human papillomaviruses (HPVs) are a worldwide burden as they are a widespread group of tumour viruses in humans. Having a tropism for mucosal tissues, high-risk HPVs are detected in nearly all cervical cancers. HPV16 is the most common high-risk type but not all women infected with high-risk HPV develop a malignant tumour. Likely relevant, HPV genomes are polymorphic and some HPV16 single nucleotide polymorphisms (SNPs) are under evolutionary constraint instigating variable oncogenicity and immunogenicity in the infected host. RESULTS: To investigate the tumourigenicity of two common HPV16 variants, we used our recently developed, three-dimensional organotypic model reminiscent of the natural HPV infectious cycle and conducted various "omics" and bioinformatics approaches. Based on epidemiological studies we chose to examine the HPV16 Asian-American (AA) and HPV16 European Prototype (EP) variants. They differ by three non-synonymous SNPs in the transforming and virus-encoded E6 oncogene where AAE6 is classified as a high- and EPE6 as a low-risk variant. Remarkably, the high-risk AAE6 variant genome integrated into the host DNA, while the low-risk EPE6 variant genome remained episomal as evidenced by highly sensitive Capt-HPV sequencing. RNA-seq experiments showed that the truncated form of AAE6, integrated in chromosome 5q32, produced a local gene over-expression and a large variety of viral-human fusion transcripts, including long distance spliced transcripts. In addition, differential enrichment of host cell pathways was observed between both HPV16 E6 variant-containing epithelia. Finally, in the high-risk variant, we detected a molecular signature of host chromosomal instability, a common property of cancer cells. CONCLUSIONS: We show how naturally occurring SNPs in the HPV16 E6 oncogene cause significant changes in the outcome of HPV infections and subsequent viral and host transcriptome alterations prone to drive carcinogenesis. Host genome instability is closely linked to viral integration into the host genome of HPV-infected cells, which is a key phenomenon for malignant cellular transformation and the reason for uncontrolled E6 oncogene expression. In particular, the finding of variant-specific integration potential represents a new paradigm in HPV variant biology.

The pelvic organs receive no parasympathetic innervation.

The pelvic organs (bladder, rectum, and sex organs) have been represented for a century as receiving autonomic innervation from two pathways - lumbar sympathetic and sacral parasympathetic - by way of a shared relay, the pelvic ganglion, conceived as an assemblage of sympathetic and parasympathetic neurons. Using single-cell RNA sequencing, we find that the mouse pelvic ganglion is made of four classes of neurons, distinct from both sympathetic and parasympathetic ones, albeit with a kinship to the former, but not the latter, through a complex genetic signature. We also show that spinal lumbar preganglionic neurons synapse in the pelvic ganglion onto equal numbers of noradrenergic and cholinergic cells, both of which therefore serve as sympathetic relays. Thus, the pelvic viscera receive no innervation from parasympathetic or typical sympathetic neurons, but instead from a divergent tail end of the sympathetic chains, in charge of its idiosyncratic functions.

An interactive murine single-cell atlas of the lung responses to radiation injury.

Radiation Induced Lung Injury (RILI) is one of the main limiting factors of thorax irradiation, which can induce acute pneumonitis as well as pulmonary fibrosis, the latter being a life-threatening condition. The order of cellular and molecular events in the progression towards fibrosis is key to the physiopathogenesis of the disease, yet their coordination in space and time remains largely unexplored. Here, we present an interactive murine single cell atlas of the lung response to irradiation, generated from C57BL6/J female mice. This tool opens the door for exploration of the spatio-temporal dynamics of the mechanisms that lead to radiation-induced pulmonary fibrosis. It depicts with unprecedented detail cell type-specific radiation-induced responses associated with either lung regeneration or the failure thereof. A better understanding of the mechanisms leading to lung fibrosis will help finding new therapeutic options that could improve patients' quality of life.

Breast cancer remotely imposes a myeloid bias on haematopoietic stem cells by reprogramming the bone marrow niche.

Myeloid cell infiltration of solid tumours generally associates with poor patient prognosis and disease severity1-13. Therefore, understanding the regulation of myeloid cell differentiation during cancer is crucial to counteract their pro-tumourigenic role. Bone marrow (BM) haematopoiesis is a tightly regulated process for the production of all immune cells in accordance to tissue needs14. Myeloid cells differentiate during haematopoiesis from multipotent haematopoietic stem and progenitor cells (HSPCs)15-17. HSPCs can sense inflammatory signals from the periphery during infections18-21 or inflammatory disorders22-27. In these settings, HSPC expansion is associated with increased myeloid differentiation28,29. During carcinogenesis, the elevation of haematopoietic growth factors supports the expansion and differentiation of committed myeloid progenitors5,30. However, it is unclear whether cancer-related inflammation also triggers demand-adapted haematopoiesis at the level of multipotent HSPCs. In the BM, HSPCs reside within the haematopoietic niche which delivers HSC maintenance and differentiation cues31-35. Mesenchymal stem cells (MSCs) are a major cellular component of the BM niche and contribute to HSC homeostasis36-41. Modifications of MSCs in systemic disorders have been associated with HSC differentiation towards myeloid cells22,42. It is unknown if MSCs are regulated in the context of solid tumours and if their myeloid supportive activity is impacted by cancer-induced systemic changes. Here, using unbiased transcriptomic analysis and in situ imaging of HSCs and the BM niche during breast cancer, we show that both HSCs and MSCs are transcriptionally and spatially modified. We demonstrate that breast tumour can distantly remodel the cellular cross-talks in the BM niche leading to increased myelopoiesis.

viroCapt: A Bioinformatics Pipeline for Identifying Viral Insertion in Human Host Genome.

INTRODUCTION: The implication of viruses in human cancers, as well as the emergence of next generation sequencing has permitted to investigate further their role and pathophysiology in the development of this disease. One such mechanism is the integration of portions of viral genomes in the human genome, as well as the specific action of viral oncogenes.inding integration sites and preserved oncogenes is still relying on heavy manual intervention. METHODS: We developed an analysis and interpretation pipeline to determine viral insertions. Using data from directed viral capture, the pipeline conducts a crude genotyping phase to select reference viral genomes, identifies chimeric reads, extracts the putative human sequences to locate in the human reference genome, scores and ranks candidate junctions, and exports tabular and visual results. RESULTS: We leverage common bioinformatics tools (bowtie2, samtools, blat), and a dedicated filtering and ranking algorithm, implemented in R, to infer candidate junctions and insertions. Static results (tables, figures) are produced, as well as an interactive interpretation tool developed as a shiny web app. DISCUSSION: We validated this pipeline against published results of HPV, HBV, and AAV2 insertions and show good information retrieval.

DNMT3A-dependent DNA methylation is required for spermatogonial stem cells to commit to spermatogenesis.

DNA methylation plays a critical role in spermatogenesis, as evidenced by the male sterility of DNA methyltransferase (DNMT) mutant mice. Here, we report a division of labor in the establishment of the methylation landscape of male germ cells and its functions in spermatogenesis. Although DNMT3C is essential for preventing retrotransposons from interfering with meiosis, DNMT3A broadly methylates the genome (with the exception of DNMT3C-dependent retrotransposons) and controls spermatogonial stem cell (SSC) plasticity. By reconstructing developmental trajectories through single-cell RNA sequencing and profiling chromatin states, we found that Dnmt3A mutant SSCs can only self-renew and no longer differentiate in association with spurious enhancer activation that enforces an irreversible stem cell gene program. Our findings therefore highlight a key function of DNA methylation in male fertility: the epigenetic programming of SSC commitment to differentiation and lifelong spermatogenesis supply.

Human papilloma virus integration sites and genomic signatures in head and neck squamous cell carcinoma.

A prevalence of around 26% of human papillomavirus (HPV) in head and neck squamous cell carcinoma (HNSCC) has been previously reported. HPV induced oncogenesis mainly involving E6 and E7 viral oncoproteins. In some cases, HPV viral DNA has been detected to integrate with the host genome and possibly contributes to carcinogenesis by affecting the gene expression. We retrospectively assessed HPV integration sites and signatures in 80 HPV positive patients with HNSCC, by using a double capture-HPV method followed by next-generation Sequencing. We detected HPV16 in 90% of the analyzed cohort and confirmed five previously described mechanistic signatures of HPV integration [episomal (EPI), integrated in a truncated form revealing two HPV-chromosomal junctions colinear (2J-COL) or nonlinear (2J-NL), multiple hybrid junctions clustering in a single chromosomal region (MJ-CL) or scattered over different chromosomal regions (MJ-SC) of the human genome]. Our results suggested that HPV remained episomal in 38.8% of the cases or was integrated/mixed in the remaining 61.2% of patients with HNSCC. We showed a lack of association of HPV genomic signatures to tumour and patient characteristics, as well as patient survival. Similar to other HPV associated cancers, low HPV copy number was associated with worse prognosis. We identified 267 HPV-human junctions scattered on most chromosomes. Remarkably, we observed four recurrent integration regions: PDL1/PDL2/PLGRKT (8.2%), MYC/PVT1 (6.1%), MACROD2 (4.1%) and KLF5/KLF12 regions (4.1%). We detected the overexpression of PDL1 and MYC upon integration by gene expression analysis. In conclusion, we identified recurrent targeting of several cancer genes such as PDL1 and MYC upon HPV integration, suggesting a role of altered gene expression by HPV integration during HNSCC carcinogenesis.

Contribution of resident and circulating precursors to tumor-infiltrating CD8+ T cell populations in lung cancer.

Tumor-infiltrating lymphocytes (TILs), in general, and especially CD8+ TILs, represent a favorable prognostic factor in non-small cell lung cancer (NSCLC). The tissue origin, regenerative capacities, and differentiation pathways of TIL subpopulations remain poorly understood. Using a combination of single-cell RNA and T cell receptor (TCR) sequencing, we investigate the functional organization of TIL populations in primary NSCLC. We identify two CD8+ TIL subpopulations expressing memory-like gene modules: one is also present in blood (circulating precursors) and the other one in juxtatumor tissue (tissue-resident precursors). In tumors, these two precursor populations converge through a unique transitional state into terminally differentiated cells, often referred to as dysfunctional or exhausted. Differentiation is associated with TCR expansion, and transition from precursor to late-differentiated states correlates with intratumor T cell cycling. These results provide a coherent working model for TIL origin, ontogeny, and functional organization in primary NSCLC.

Episomal HPV16 responsible for aggressive and deadly metastatic anal squamous cell carcinoma evidenced in peripheral blood.

Archival tissue samples collected longitudinally from a patient who died from HPV16-induced high-grade anal intraepithelial squamous cell carcinoma with vertebral HPV16-positive metastasis were retrospectively analyzed by the Capture-HPV method (Capt-HPV) followed by Next-Generation Sequencing (NGS). Full length nucleotide sequences of the same HPV16 were identified from the initial and second anal biopsy samples, from plasma sample and from vertebral metastasis biopsy. Remarkably, HPV was episomal in each sample. The HPV genome sequence was closest to the HPV16 Qv18158E variant subtype (A1 lineage) exhibiting base substitutions and deletions in 7 and 2 HPV loci, respectively. In conclusion, the powerful Capt-HPV followed by NGS allows evidencing the detailed cartography of tumoral and circulating HPV DNA, giving rise to a unique and unexpected episomal virus molecular status in a context of aggressive carcinoma, underlying the importance of HPV status and its association with clinical features for further prospective studies.

Classification of the Zoonotic Hepatitis E Virus Genotype 3 Into Distinct Subgenotypes.

Hepatitis E virus (HEV) genotype 3 is the most common genotype linked to HEV infections in Europe and America. Three major clades (HEV-3.1, HEV-3.2, and HEV-3.3) have been identified but the overlaps between intra-subtype and inter-subtype p-distances make subtype classification inconsistent. Reference sequences have been proposed to facilitate communication between researchers and new putative subtypes have been identified recently. We have used the full or near full-length HEV-3 genome sequences available in the Genbank database (April 2020; n = 503) and distance analyses of clades HEV-3.1 and HEV-3.2 to determine a p-distance cut-off (0.093 nt substitutions/site) in order to define subtypes. This could help to harmonize HEV-3 genotyping, facilitate molecular epidemiology studies and investigations of the biological and clinical differences between HEV-3 subtypes.