Highlighting fibroblast plasticity in lung fibrosis: the WI-38 cell line as a model for investigating the myofibroblast and lipofibroblast switch.

Background: Myofibroblasts (MYFs) are generally considered the principal culprits in excessive extracellular matrix deposition and scar formation in the pathogenesis of lung fibrosis. Lipofibroblasts (LIFs), on the other hand, are defined by their lipid-storing capacity and are predominantly found in the alveolar regions of the lung. They have been proposed to play a protective role in lung fibrosis. We previously reported that a LIF to MYF reversible differentiation switch occurred during fibrosis formation and resolution. In this study, we tested whether WI-38 cells, a human embryonic lung fibroblast cell line, could be used to study fibroblast differentiation towards the LIF or MYF phenotype and whether this could be relevant for idiopathic pulmonary fibrosis (IPF). Methods: Using WI-38 cells, Fibroblast (FIB) to MYF differentiation was triggered using TGF-ß1 treatment and FIB to LIF differentiation using Metformin treatment. We also analyzed the MYF to LIF and LIF to MYF differentiation by pre-treating the WI-38 cells with TGF-ß1 or Metformin respectively. We used IF, qPCR and bulk RNA-Seq to analyze the phenotypic and transcriptomic changes in the cells. We correlated our in vitro transcriptome data from WI-38 cells (obtained via bulk RNA sequencing) with the transcriptomic signature of LIFs and MYFs derived from the IPF cell atlas as well as with our own single-cell transcriptomic data from IPF patients-derived lung fibroblasts (LF-IPF) cultured in vitro. We also carried out alveolosphere assays to evaluate the ability of the proposed LIF and MYF cells to support the growth of alveolar epithelial type 2 cells. Results: WI-38 cells and LF-IPF display similar phenotypical and gene expression responses to TGF-ß1 and Metformin treatment. Bulk RNA-Seq analysis of WI-38 cells and LF-IPF treated with TGF-ß1, or Metformin indicate similar transcriptomic changes. We also show the partial conservation of the LIF and MYF signature extracted from the Habermann et al. scRNA-seq dataset in WI-38 cells treated with Metformin or TGF-ß1, respectively. Alveolosphere assays indicate that LIFs enhance organoid growth, while MYFs inhibit organoid growth. Finally, we provide evidence supporting the MYF to LIF and LIF to MYF reversible switch using WI-38 cells. Conclusions: WI-38 cells represent a versatile and reliable model to study the intricate dynamics of fibroblast differentiation towards the MYF or LIF phenotype associated with lung fibrosis formation and resolution, providing valuable insights to drive future research.

Early transcriptional similarities between two distinct neural lineages during ascidian embryogenesis.

In chordates, the central nervous system arises from precursors that have distinct developmental and transcriptional trajectories. Anterior nervous systems are ontogenically associated with ectodermal lineages while posterior nervous systems are associated with mesoderm. Taking advantage of the well-documented cell lineage of ascidian embryos, we asked to what extent the transcriptional states of the different neural lineages become similar during the course of progressive lineage restriction. We performed single-cell RNA sequencing (scRNA-seq) analyses on hand-dissected neural precursor cells of the two distinct lineages, together with those of their sister cell lineages, with a high temporal resolution covering five successive cell cycles from the 16-cell to neural plate stages. A transcription factor binding site enrichment analysis of neural specific genes at the neural plate stage revealed limited evidence for shared transcriptional control between the two neural lineages, consistent with their different ontogenies. Nevertheless, PCA analysis and hierarchical clustering showed that, by neural plate stages, the two neural lineages cluster together. Consistent with this, we identified a set of genes enriched in both neural lineages at the neural plate stage, including miR-124, Celf3.a, Zic.r-b, and Ets1/2. Altogether, the current study has revealed genome-wide transcriptional dynamics of neural progenitor cells of two distinct developmental origins. Our scRNA-seq dataset is unique and provides a valuable resource for future analyses, enabling a precise temporal resolution of cell types not previously described from dissociated embryos.

Cell Culture Differentiation and Proliferation Conditions Influence the In Vitro Regeneration of the Human Airway Epithelium.

The human airway mucociliary epithelium can be recapitulated in vitro using primary cells cultured in an Air-Liquid Interface (ALI), a reliable surrogate to perform pathophysiological studies. As tremendous variations exist between media used for ALI-cultured human airway epithelial cells, our study aimed to evaluate the impact of several media (BEGMTM, PneumaCultTM, "Half&Half" and "Clancy") on cell type distribution using single-cell RNA sequencing and imaging. Our work revealed the impact of these media on cell composition, gene expression profile, cell signaling and epithelial morphology. We found higher proportions of multiciliated cells in PneumaCultTM-ALI and Half&Half, stronger EGF signaling from basal cells in BEGMTM-ALI, differential expression of the SARS-CoV-2 entry factor ACE2, and distinct secretome transcripts depending on media used. We also established that proliferation in PneumaCultTM-Ex Plus favored secretory cell fate, showing the key influence of proliferation media on late differentiation epithelial characteristics. Altogether, our data offer a comprehensive repertoire for evaluating the effects of culture conditions on airway epithelial differentiation and will help to choose the most relevant medium according to the processes to be investigated such as cilia, mucus biology or viral infection. We detail useful parameters that should be explored to document airway epithelial cell fate and morphology. This article is open access and distributed under the terms of the Creative Commons Attribution Non-Commercial No Derivatives License 4.0 (http://creativecommons.org/licenses/by-nc-nd/4.0/).

Identifying Candidate Gene Drivers Associated with Relapse in Pediatric T-Cell Acute Lymphoblastic Leukemia Using a Gene Co-Expression Network Approach.

Pediatric T-cell Acute Lymphoblastic Leukemia (T-ALL) relapses are still associated with a dismal outcome, justifying the search for new therapeutic targets and relapse biomarkers. Using single-cell RNA sequencing (scRNAseq) data from three paired samples of pediatric T-ALL at diagnosis and relapse, we first conducted a high-dimensional weighted gene co-expression network analysis (hdWGCNA). This analysis highlighted several gene co-expression networks (GCNs) and identified relapse-associated hub genes, which are considered potential driver genes. Shared relapse-expressed genes were found to be related to antigen presentation (HLA, B2M), cytoskeleton remodeling (TUBB, TUBA1B), translation (ribosomal proteins, EIF1, EEF1B2), immune responses (MIF, EMP3), stress responses (UBC, HSP90AB1/AA1), metabolism (FTH1, NME1/2, ARCL4C), and transcriptional remodeling (NF-κB family genes, FOS-JUN, KLF2, or KLF6). We then utilized sparse partial least squares discriminant analysis to select from a pool of 481 unique leukemic hub genes, which are the genes most discriminant between diagnosis and relapse states (comprising 44, 35, and 31 genes, respectively, for each patient). Applying a Cox regression method to these patient-specific genes, along with transcriptomic and clinical data from the TARGET-ALL AALL0434 cohort, we generated three model gene signatures that efficiently identified relapsed patients within the cohort. Overall, our approach identified new potential relapse-associated genes and proposed three model gene signatures associated with lower survival rates for high-score patients.

Detecting subtle transcriptomic perturbations induced by lncRNAs knock-down in single-cell CRISPRi screening using a new sparse supervised autoencoder neural network.

Single-cell CRISPR-based transcriptome screens are potent genetic tools for concomitantly assessing the expression profiles of cells targeted by a set of guides RNA (gRNA), and inferring target gene functions from the observed perturbations. However, due to various limitations, this approach lacks sensitivity in detecting weak perturbations and is essentially reliable when studying master regulators such as transcription factors. To overcome the challenge of detecting subtle gRNA induced transcriptomic perturbations and classifying the most responsive cells, we developed a new supervised autoencoder neural network method. Our Sparse supervised autoencoder (SSAE) neural network provides selection of both relevant features (genes) and actual perturbed cells. We applied this method on an in-house single-cell CRISPR-interference-based (CRISPRi) transcriptome screening (CROP-Seq) focusing on a subset of long non-coding RNAs (lncRNAs) regulated by hypoxia, a condition that promote tumor aggressiveness and drug resistance, in the context of lung adenocarcinoma (LUAD). The CROP-seq library of validated gRNA against a subset of lncRNAs and, as positive controls, HIF1A and HIF2A, the 2 main transcription factors of the hypoxic response, was transduced in A549 LUAD cells cultured in normoxia or exposed to hypoxic conditions during 3, 6 or 24 h. We first validated the SSAE approach on HIF1A and HIF2 by confirming the specific effect of their knock-down during the temporal switch of the hypoxic response. Next, the SSAE method was able to detect stable short hypoxia-dependent transcriptomic signatures induced by the knock-down of some lncRNAs candidates, outperforming previously published machine learning approaches. This proof of concept demonstrates the relevance of the SSAE approach for deciphering weak perturbations in single-cell transcriptomic data readout as part of CRISPR-based screening.

Cholesterol efflux pathways hinder KRAS-driven lung tumor progenitor cell expansion.

Cholesterol efflux pathways could be exploited in tumor biology to unravel cancer vulnerabilities. A mouse model of lung-tumor-bearing KRASG12D mutation with specific disruption of cholesterol efflux pathways in epithelial progenitor cells promoted tumor growth. Defective cholesterol efflux in epithelial progenitor cells governed their transcriptional landscape to support their expansion and create a pro-tolerogenic tumor microenvironment (TME). Overexpression of the apolipoprotein A-I, to raise HDL levels, protected these mice from tumor development and dire pathologic consequences. Mechanistically, HDL blunted a positive feedback loop between growth factor signaling pathways and cholesterol efflux pathways that cancer cells hijack to expand. Cholesterol removal therapy with cyclodextrin reduced tumor burden in progressing tumor by suppressing the proliferation and expansion of epithelial progenitor cells of tumor origin. Local and systemic perturbations of cholesterol efflux pathways were confirmed in human lung adenocarcinoma (LUAD). Our results position cholesterol removal therapy as a putative metabolic target in lung cancer progenitor cells.

The MIR34B/C genomic region contains multiple potential regulators of multiciliogenesis.

The MIR449 genomic locus encompasses several regulators of multiciliated cell (MCC) formation (multiciliogenesis). The miR-449 homologs miR-34b/c represent additional regulators of multiciliogenesis that are transcribed from another locus. Here, we characterized the expression of BTG4, LAYN, and HOATZ, located in the MIR34B/C locus using single-cell RNA-seq and super-resolution microscopy from human, mouse, or pig multiciliogenesis models. BTG4, LAYN, and HOATZ transcripts were expressed in both precursors and mature MCCs. The Layilin/LAYN protein was absent from primary cilia, but it was expressed in apical membrane regions or throughout motile cilia. LAYN silencing altered apical actin cap formation and multiciliogenesis. HOATZ protein was detected in primary cilia or throughout motile cilia. Altogether, our data suggest that the MIR34B/C locus may gather potential actors of multiciliogenesis.

The spatial landscape of gene expression isoforms in tissue sections.

In situ capturing technologies add tissue context to gene expression data, with the potential of providing a greater understanding of complex biological systems. However, splicing variants and full-length sequence heterogeneity cannot be characterized at spatial resolution with current transcriptome profiling methods. To that end, we introduce spatial isoform transcriptomics (SiT), an explorative method for characterizing spatial isoform variation and sequence heterogeneity using long-read sequencing. We show in mouse brain how SiT can be used to profile isoform expression and sequence heterogeneity in different areas of the tissue. SiT reveals regional isoform switching of Plp1 gene between different layers of the olfactory bulb, and the use of external single-cell data allows the nomination of cell types expressing each isoform. Furthermore, SiT identifies differential isoform usage for several major genes implicated in brain function (Snap25, Bin1, Gnas) that are independently validated by in situ sequencing. SiT also provides for the first time an in-depth A-to-I RNA editing map of the adult mouse brain. Data exploration can be performed through an online resource (https://www.isomics.eu), where isoform expression and RNA editing can be visualized in a spatial context.

GLI1+ cells are a source of repair-supportive mesenchymal cells (RSMCs) during airway epithelial regeneration.

Repair-supportive mesenchymal cells (RSMCs) have been recently reported in the context of naphthalene (NA)-induced airway injury and regeneration. These cells transiently express smooth muscle actin (Acta2) and are enriched with platelet-derived growth factor receptor alpha (Pdgfra) and fibroblast growth factor 10 (Fgf10) expression. Genetic deletion of Ctnnb1 (gene coding for beta catenin) or Fgf10 in these cells using the Acta2-Cre-ERT2 driver line after injury (defined as NA-Tam condition; Tam refers to tamoxifen) led to impaired repair of the airway epithelium. In this study, we demonstrate that RSMCs are mostly captured using the Acta2-Cre-ERT2 driver when labeling occurs after (NA-Tam condition) rather than before injury (Tam-NA condition), and that their expansion occurs mostly between days 3 and 7 following NA treatment. Previous studies have shown that lineage-traced peribronchial GLI1+ cells are transiently amplified after NA injury. Here, we report that Gli1 expression is enriched in RSMCs. Using lineage tracing with Gli1Cre-ERT2 mice combined with genetic inactivation of Fgf10, we show that GLI1+ cells with Fgf10 deletion fail to amplify around the injured airways, thus resulting in impaired airway epithelial repair. Interestingly, Fgf10 expression is not upregulated in GLI1+ cells following NA treatment, suggesting that epithelial repair is mostly due to the increased number of Fgf10-expressing GLI1+ cells. Co-culture of SCGB1A1+ cells with GLI1+ cells isolated from non-injured or injured lungs showed that GLI1+ cells from these two conditions are similarly capable of supporting bronchiolar organoid (or bronchiolosphere) formation. Single-cell RNA sequencing on sorted lineage-labeled cells showed that the RSMC signature resembles that of alveolar fibroblasts. Altogether, our study provides strong evidence for the involvement of mesenchymal progenitors in airway epithelial regeneration and highlights the critical role played by Fgf10-expressing GLI1+ cells in this context.

Coupling live-cell imaging and in situ isolation of the same single cell to profile the transient states of predicted drug-tolerant cells.

Cell response variability is a starting point in cancer drug resistance that has been difficult to analyze because the tolerant cell states are short lived. Here, we present fate-seq, an approach to isolate single cells in their transient states of drug sensitivity or tolerance before profiling. The drug response is predicted in live cells, which are laser-captured by microdissection before any drug-induced change can alter their states. This framework enables the identification of the cell-state signatures causing differential cell decisions upon treatment. For complete details on the use and execution of this protocol, please refer to Meyer et al. (2020).

Persistent Properties of a Subpopulation of Cancer Cells Overexpressing the Hedgehog Receptor Patched.

Despite the development of new therapeutic strategies, cancer remains one of the leading causes of mortality worldwide. One of the current major challenges is the resistance of cancers to chemotherapy treatments inducing metastases and relapse of the tumor. The Hedgehog receptor Patched (Ptch1) is overexpressed in many types of cancers. We showed that Ptch1 contributes to the efflux of doxorubicin and plays an important role in the resistance to chemotherapy in adrenocortical carcinoma (ACC), a rare cancer which presents strong resistance to the standard of care chemotherapy treatment. In the present study, we isolated and characterized a subpopulation of the ACC cell line H295R in which Ptch1 is overexpressed and more present at the cell surface. This cell subpopulation is more resistant to doxorubicin, grows as spheroids, and has a greater capability of clonogenicity, migration, and invasion than the parental cells. Xenograft experiments performed in mice and in ovo showed that this cell subpopulation is more tumorigenic and metastatic than the parental cells. These results suggest that this cell subpopulation has cancer stem-like or persistent cell properties which were strengthened by RNA-seq. If present in tumors from ACC patients, these cells could be responsible for therapy resistance, relapse, and metastases.

The FibromiR miR-214-3p Is Upregulated in Duchenne Muscular Dystrophy and Promotes Differentiation of Human Fibro-Adipogenic Muscle Progenitors.

Fibrosis is a deleterious invasion of tissues associated with many pathological conditions, such as Duchenne muscular dystrophy (DMD) for which no cure is at present available for its prevention or its treatment. Fibro-adipogenic progenitors (FAPs) are resident cells in the human skeletal muscle and can differentiate into myofibroblasts, which represent the key cell population responsible for fibrosis. In this study, we delineated the pool of microRNAs (miRNAs) that are specifically modulated by TGFß1 in FAPs versus myogenic progenitors (MPs) by a global miRNome analysis. A subset of candidates, including several "FibromiRs", was found differentially expressed between FAPs and MPs and was also deregulated in DMD versus healthy biopsies. Among them, the expression of the TGFß1-induced miR-199a~214 cluster was strongly correlated with the fibrotic score in DMD biopsies. Loss-of-function experiments in FAPs indicated that a miR-214-3p inhibitor efficiently blocked expression of fibrogenic markers in both basal conditions and following TGFß1 stimulation. We found that FGFR1 is a functional target of miR-214-3p, preventing the signaling of the anti-fibrotic FGF2 pathway during FAP fibrogenesis. Overall, our work demonstrates that the « FibromiR ¼ miR-214-3p is a key activator of FAP fibrogenesis by modulating the FGF2/FGFR1/TGFß axis, opening new avenues for the treatment of DMD.

Monitoring SARS-CoV-2 variants alterations in Nice neighborhoods by wastewater nanopore sequencing.

BACKGROUND: Wastewater surveillance was proposed as an epidemiological tool to define the prevalence and evolution of the SARS-CoV-2 epidemics. However, most implemented SARS-CoV-2 wastewater surveillance projects were based on qPCR measurement of virus titers and did not address the mutational spectrum of SARS-CoV-2 circulating in the population. METHODS: We have implemented a nanopore RNA sequencing monitoring system in the city of Nice (France, 550,000 inhabitants). Between October 2020 and March 2021, we monthly analyzed the SARS-CoV-2 variants in 113 wastewater samples collected in the main wastewater treatment plant and 20 neighborhoods. FINDINGS: We initially detected the lineages predominant in Europe at the end of 2020 (B.1.160, B.1.177, B.1.367, B.1.474, and B.1.221). In January, a localized emergence of a variant (Spike:A522S) of the B.1.1.7 lineage occurred in one neighborhood. It rapidly spread and became dominant all over the city. Other variants of concern (B.1.351, P.1) were also detected in some neighborhoods, but at low frequency. Comparison with individual clinical samples collected during the same week showed that wastewater sequencing correctly identified the same lineages as those found in COVID-19 patients. INTERPRETATION: Wastewater sequencing allowed to document the diversity of SARS-CoV-2 sequences within the different neighborhoods of the city of Nice. Our results illustrate how sequencing of sewage samples can be used to track pathogen sequence diversity in the current pandemics and in future infectious disease outbreaks. TRANSLATION: For the French translation of the abstract see Supplementary Materials section.

Correction: The nuclear hypoxia-regulated NLUCAT1 long non-coding RNA contributes to an aggressive phenotype in lung adenocarcinoma through regulation of oxidative stress.

Lung cancer is the leading cause of cancer death worldwide, with poor prognosis and a high rate of recurrence despite early surgical removal. Hypoxic regions within tumors represent sources of aggressiveness and resistance to therapy. Although long non-coding RNAs (lncRNAs) are increasingly recognized as major gene expression regulators, their regulation and function following hypoxic stress are still largely unexplored. Combining profiling studies on early-stage lung adenocarcinoma (LUAD) biopsies and on A549 LUAD cell lines cultured in normoxic or hypoxic conditions, we identified a subset of lncRNAs that are both correlated with the hypoxic status of tumors and regulated by hypoxia in vitro. We focused on a new transcript, Nuclear LUCAT1 (NLUCAT1), which is strongly upregulated by hypoxia in vitro and correlated with hypoxic markers and poor prognosis in LUADs. Full molecular characterization showed that NLUCAT1 is a large nuclear transcript composed of six exons and mainly regulated by NF-κB and NRF2 transcription factors. CRISPR-Cas9-mediated invalidation of NLUCAT1 revealed a decrease in proliferative and invasive properties, an increase in oxidative stress and a higher sensitivity to cisplatin-induced apoptosis. Transcriptome analysis of NLUCAT1-deficient cells showed repressed genes within the antioxidant and/or cisplatin-response networks. We demonstrated that the concomitant knockdown of four of these genes products, GPX2, GLRX, ALDH3A1, and PDK4, significantly increased ROS-dependent caspase activation, thus partially mimicking the consequences of NLUCAT1 inactivation in LUAD cells. Overall, we demonstrate that NLUCAT1 contributes to an aggressive phenotype in early-stage hypoxic tumors, suggesting it may represent a new potential therapeutic target in LUADs.

Single-cell RNA sequencing reveals intratumoral heterogeneity in primary uveal melanomas and identifies HES6 as a driver of the metastatic disease.

Intratumor heterogeneity has been recognized in numerous cancers as a major source of metastatic dissemination. In uveal melanomas, the existence and identity of specific subpopulations, their biological function and their contribution to metastasis remain unknown. Here, in multiscale analyses using single-cell RNA sequencing of six different primary uveal melanomas, we uncover an intratumoral heterogeneity at the genomic and transcriptomic level. We identify distinct transcriptional cell states and diverse tumor-associated populations in a subset of the samples. We also decipher a gene regulatory network underlying an invasive and poor prognosis state driven in part by the transcription factor HES6. HES6 heterogenous expression has been validated by RNAscope assays within primary human uveal melanomas, which further unveils the existence of these cells conveying a dismal prognosis in tumors diagnosed with a favorable outcome using bulk analyses. Depletion of HES6 impairs proliferation, migration and metastatic dissemination in vitro and in vivo using the chick chorioallantoic membrane assay, demonstrating the essential role of HES6 in uveal melanomas. Thus, single-cell analysis offers an unprecedented view of primary uveal melanoma heterogeneity, identifies bona fide biomarkers for metastatic cells in the primary tumor, and reveals targetable modules driving growth and metastasis formation. Significantly, our findings demonstrate that HES6 is a valid target to stop uveal melanoma progression.

Evidence of early increased sialylation of airway mucins and defective mucociliary clearance in CFTR-deficient piglets.

BACKGROUND: Bacterial colonization in cystic fibrosis (CF) lungs has been directly associated to the loss of CFTR function, and/or secondarily linked to repetitive cycles of chronic inflammation/infection. We hypothesized that altered molecular properties of mucins could contribute to this process. METHODS: Newborn CFTR+/+ and CFTR-/- were sacrificed before and 6 h after inoculation with luminescent Pseudomonas aeruginosa into the tracheal carina. Tracheal mucosa and the bronchoalveolar lavage (BAL) fluid were collected to determine the level of mucin O-glycosylation, bacteria binding to mucins and the airways transcriptome. Disturbances in mucociliary transport were determined by ex-vivo imaging of luminescent Pseudomonas aeruginosa. RESULTS: We provide evidence of an increased sialylation of CF airway mucins and impaired mucociliary transport that occur before the onset of inflammation. Hypersialylation of mucins was reproduced on tracheal explants from non CF animals treated with GlyH101, an inhibitor of CFTR channel activity, indicating a causal relationship between the absence of CFTR expression and the sialylation of mucins. This increased sialylation was correlated to an increased adherence of P. aeruginosa to mucins. In vivo infection of newborn CF piglets by live luminescent P. aeruginosa demonstrated an impairment of mucociliary transport of this bacterium, with no evidence of pre-existing inflammation. CONCLUSIONS: Our results document for the first time in a well-defined CF animal model modifications that affect the O-glycan chains of mucins. These alterations precede infection and inflammation of airway tissues, and provide a favorable context for microbial development in CF lung that hallmarks this disease.

Identification of a Repair-Supportive Mesenchymal Cell Population during Airway Epithelial Regeneration.

Tissue regeneration requires coordinated and dynamic remodeling of stem and progenitor cells and the surrounding niche. Although the plasticity of epithelial cells has been well explored in many tissues, the dynamic changes occurring in niche cells remain elusive. Here, we show that, during lung repair after naphthalene injury, a population of PDGFRα+ cells emerges in the non-cartilaginous conducting airway niche, which is normally populated by airway smooth muscle cells (ASMCs). This cell population, which we term "repair-supportive mesenchymal cells" (RSMCs), is distinct from conventional ASMCs, which have previously been shown to contribute to epithelial repair. Gene expression analysis on sorted lineage-labeled cells shows that RSMCs express low levels of ASMC markers, but high levels of the pro-regenerative marker Fgf10. Organoid co-cultures demonstrate an enhanced ability for RSMCs in supporting club-cell growth. Our study highlights the dynamics of mesenchymal cells in the airway niche and has implications for chronic airway-injury-associated diseases.

Longevity strategies in response to light in the reef coral Stylophora pistillata.

Aging is a multifactorial process that results in progressive loss of regenerative capacity and tissue function while simultaneously favoring the development of a large array of age-related diseases. Evidence suggests that the accumulation of senescent cells in tissue promotes both normal and pathological aging. Oxic stress is a key driver of cellular senescence. Because symbiotic long-lived reef corals experience daily hyperoxic and hypoxic transitions, we hypothesized that these long-lived animals have developed specific longevity strategies in response to light. We analyzed transcriptome variation in the reef coral Stylophora pistillata during the day-night cycle and revealed a signature of the FoxO longevity pathway. We confirmed this pathway by immunofluorescence using antibodies against coral FoxO to demonstrate its nuclear translocation. Through qPCR analysis of nycthemeral variations of candidate genes under different light regimens, we found that, among genes that were specifically up- or downregulated upon exposure to light, human orthologs of two "light-up" genes (HEY1 and LONF3) exhibited anti-senescence properties in primary human fibroblasts. Therefore, these genes are interesting candidates for counteracting skin aging. We propose a large screen for other light-up genes and an investigation of the biological response of reef corals to light (e.g., metabolic switching) to elucidate these processes and identify effective interventions for promoting healthy aging in humans.

Identification of oncolytic vaccinia restriction factors in canine high-grade mammary tumor cells using single-cell transcriptomics.

Mammary carcinoma, including triple-negative breast carcinomas (TNBC) are tumor-types for which human and canine pathologies are closely related at the molecular level. The efficacy of an oncolytic vaccinia virus (VV) was compared in low-passage primary carcinoma cells from TNBC versus non-TNBC. Non-TNBC cells were 28 fold more sensitive to VV than TNBC cells in which VV replication is impaired. Single-cell RNA-seq performed on two different TNBC cell samples, infected or not with VV, highlighted three distinct populations: naïve cells, bystander cells, defined as cells exposed to the virus but not infected and infected cells. The transcriptomes of these three populations showed striking variations in the modulation of pathways regulated by cytokines and growth factors. We hypothesized that the pool of genes expressed in the bystander populations was enriched in antiviral genes. Bioinformatic analysis suggested that the reduced activity of the virus was associated with a higher mesenchymal status of the cells. In addition, we demonstrated experimentally that high expression of one gene, DDIT4, is detrimental to VV production. Considering that DDIT4 is associated with a poor prognosis in various cancers including TNBC, our data highlight DDIT4 as a candidate resistance marker for oncolytic poxvirus therapy. This information could be used to design new generations of oncolytic poxviruses. Beyond the field of gene therapy, this study demonstrates that single-cell transcriptomics can be used to identify cellular factors influencing viral replication.

Profiling the Non-genetic Origins of Cancer Drug Resistance with a Single-Cell Functional Genomics Approach Using Predictive Cell Dynamics.

Non-genetic heterogeneity observed in clonal cell populations is an immediate cause of drug resistance that remains challenging to profile because of its transient nature. Here, we coupled three single-cell technologies to link the predicted drug response of a cell to its own genome-wide transcriptomic profile. As a proof of principle, we analyzed the response to tumor-necrosis-factor-related apoptosis-inducing ligand (TRAIL) in HeLa cells to demonstrate that cell dynamics can discriminate the transient transcriptional states at the origin of cell decisions such as sensitivity and resistance. Our same-cell approach, named fate-seq, can reveal the molecular factors regulating the efficacy of a drug in clonal cells, providing therapeutic targets of non-genetic drug resistance otherwise confounded in gene expression noise. A record of this paper's transparent peer review process is included in the Supplemental Information.