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.

Recessively Inherited Deficiency of Secreted WFDC2 (HE4) Causes Nasal Polyposis and Bronchiectasis.

RATIONALE: Bronchiectasis is a pathological dilatation of the bronchi in the respiratory airways associated with environmental or genetic causes (e.g., cystic fibrosis, primary ciliary dyskinesia and primary immunodeficiency disorders), but most cases remain idiopathic. OBJECTIVES: To identify novel genetic defects in unsolved cases of bronchiectasis presenting with severe rhinosinusitis, nasal polyposis, and pulmonary Pseudomonas aeruginosa infection. METHODS: DNA was analyzed by next-generation or targeted Sanger sequencing. RNA was analyzed by quantitative PCR and single-cell RNA sequencing. Patient-derived, cells, cell cultures and secretions (mucus, saliva, seminal fluid) were analyzed by Western blotting and immunofluorescence microscopy, and mucociliary activity was measured. Blood serum was analyzed by electrochemiluminescence immunoassay. Protein structure and proteomic analyses were used to assess the impact of a disease-causing founder variant. MEASUREMENTS AND MAIN RESULTS: We identified bi-allelic pathogenic variants in WFDC2 in 11 individuals from 10 unrelated families originating from the United States, Europe, Asia, and Africa. Expression of WFDC2 was detected predominantly in secretory cells of control airway epithelium and also in submucosal glands. We demonstrate that WFDC2 is below the limit of detection in blood serum and hardly detectable in samples of saliva, seminal fluid, and airway surface liquid from WFDC2-deficient individuals. Computer simulations and deglycosylation assays indicate that the disease-causing founder variant p.Cys49Arg structurally hampers glycosylation and thus secretion of mature WFDC2. CONCLUSIONS: WFDC2 dysfunction defines a novel molecular etiology of bronchiectasis characterized by the deficiency of a secreted component of the airways. A commercially available blood test combined with genetic testing allows its diagnosis. This article is open access and distributed under the terms of the Creative Commons Attribution 4.0 International License (https://creativecommons.org/licenses/by/4.0/).

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.

BABA-induced pathogen resistance: a multi-omics analysis of the tomato response reveals a hyper-receptive status involving ethylene.

Prior exposure to microbial-associated molecular patterns or specific chemical compounds can promote plants into a primed state with stronger defence responses. ß-aminobutyric acid (BABA) is an endogenous stress metabolite that induces resistance protecting various plants towards diverse stresses. In this study, by integrating BABA-induced changes in selected metabolites with transcriptome and proteome data, we generated a global map of the molecular processes operating in BABA-induced resistance (BABA-IR) in tomato. BABA significantly restricts the growth of the pathogens Oidium neolycopersici and Phytophthora parasitica but not Botrytis cinerea. A cluster analysis of the upregulated processes showed that BABA acts mainly as a stress factor in tomato. The main factor distinguishing BABA-IR from other stress conditions was the extensive induction of signaling and perception machinery playing a key role in effective resistance against pathogens. Interestingly, the signalling processes and immune response activated during BABA-IR in tomato differed from those in Arabidopsis with substantial enrichment of genes associated with jasmonic acid (JA) and ethylene (ET) signalling and no change in Asp levels. Our results revealed key differences between the effect of BABA on tomato and other model plants studied until now. Surprisingly, salicylic acid (SA) is not involved in BABA downstream signalization whereas ET and JA play a crucial role.

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.

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).

Unravelling the sex-specific diversity and functions of adrenal gland macrophages.

Despite the ubiquitous function of macrophages across the body, the diversity, origin, and function of adrenal gland macrophages remain largely unknown. We define the heterogeneity of adrenal gland immune cells using single-cell RNA sequencing and use genetic models to explore the developmental mechanisms yielding macrophage diversity. We define populations of monocyte-derived and embryonically seeded adrenal gland macrophages and identify a female-specific subset with low major histocompatibility complex (MHC) class II expression. In adulthood, monocyte recruitment dominates adrenal gland macrophage maintenance in female mice. Adrenal gland macrophage sub-tissular distribution follows a sex-dimorphic pattern, with MHC class IIlow macrophages located at the cortico-medullary junction. Macrophage sex dimorphism depends on the presence of the cortical X-zone. Adrenal gland macrophage depletion results in altered tissue homeostasis, modulated lipid metabolism, and decreased local aldosterone production during stress exposure. Overall, these data reveal the heterogeneity of adrenal gland macrophages and point toward sex-restricted distribution and functions of these cells.

TLR5 signalling is hyper-responsive in porcine cystic fibrosis airways epithelium.

Excessive lung inflammation and airway epithelium damage are hallmarks of cystic fibrosis (CF) disease. It is unclear whether lung inflammation is related to an intrinsic defect in the immune response or to chronic infection. We aimed to determine whether TLR5-mediated response is defective in the CF airway epithelium. We used a newborn CF pig model to study intrinsic alterations in CF airway epithelium innate immune response. Airway epithelial cells (AECs) were stimulated with flagellin or lipopolysaccharide to determine responses specific for TLR5 and TLR4, respectively. We observed a significant increase in cytokine secretion when CF AECs were stimulated with flagellin compared to wild type (WT) AECs. These results were recapitulated when AECs were treated with an inhibitor of CFTR channel activity. We show that TLR5-signalling is altered in CF lung epithelium at birth. Modulation of TLR5 signalling could contribute to better control the excessive inflammatory response observed in CF lungs.

Brown adipose tissue monocytes support tissue expansion.

Monocytes are part of the mononuclear phagocytic system. Monocytes play a central role during inflammatory conditions and a better understanding of their dynamics might open therapeutic opportunities. In the present study, we focused on the characterization and impact of monocytes on brown adipose tissue (BAT) functions during tissue remodeling. Single-cell RNA sequencing analysis of BAT immune cells uncovered a large diversity in monocyte and macrophage populations. Fate-mapping experiments demonstrated that the BAT macrophage pool requires constant replenishment from monocytes. Using a genetic model of BAT expansion, we found that brown fat monocyte numbers were selectively increased in this scenario. This observation was confirmed using a CCR2-binding radiotracer and positron emission tomography. Importantly, in line with their tissue recruitment, blood monocyte counts were decreased while bone marrow hematopoiesis was not affected. Monocyte depletion prevented brown adipose tissue expansion and altered its architecture. Podoplanin engagement is strictly required for BAT expansion. Together, these data redefine the diversity of immune cells in the BAT and emphasize the role of monocyte recruitment for tissue remodeling.

Airway Administration of Flagellin Regulates the Inflammatory Response to Pseudomonas aeruginosa.

Excessive lung inflammation and airway epithelial damage are hallmarks of human inflammatory lung diseases, such as cystic fibrosis (CF). Enhancement of innate immunity provides protection against pathogens while reducing lung-damaging inflammation. However, the mechanisms underlying innate immunity-mediated protection in the lung remain mysterious, in part because of the lack of appropriate animal models for these human diseases. TLR5 (Toll-like receptor 5) stimulation by its specific ligand, the bacterial protein flagellin, has been proposed to enhance protection against several respiratory infectious diseases, although other cellular events, such as calcium signaling, may also control the intensity of the innate immune response. Here, we investigated the molecular events prompted by stimulation with flagellin and its role in regulating innate immunity in the lung of the pig, which is anatomically and genetically more similar to humans than rodent models. We found that flagellin treatment modulated NF-κB signaling and intracellular calcium homeostasis in airway epithelial cells. Flagellin pretreatment reduced the NF-κB nuclear translocation and the expression of proinflammatory cytokines to a second flagellin stimulus as well as to Pseudomonas aeruginosa infection. Moreover, in vivo administration of flagellin decreased the severity of P. aeruginosa-induced pneumonia. Then we confirmed these beneficial effects of flagellin in a pathological model of CF by using ex vivo precision-cut lung slices from a CF pigz model. These results provide evidence that flagellin treatment contributes to a better regulation of the inflammatory response in inflammatory lung diseases such as CF.

Cutaneous Squamous Cell Carcinoma Development Is Associated with a Temporal Infiltration of ILC1 and NK Cells with Immune Dysfunctions.

NK cells and tissue-resident innate lymphoid cells (ILCs) are innate effectors found in the skin. To investigate their temporal dynamics and specific functions throughout the development of cutaneous squamous cell carcinoma (cSCC), we combined transcriptomic and immunophenotyping analyses in mouse and human cSCCs. We identified an infiltration of NK cells and ILC1s as well as the presence of a few ILC3s. Adoptive transfer of NK cells in NK cell‒ and ILC-deficient Nfil3-/- mice revealed a role for NK cells in early control of cSCC. During tumor progression, we identified a population skewing with the infiltration of atypical ILC1 secreting inflammatory cytokines but reduced levels of IFN-γ at the papilloma stage. NK cells and ILC1s were functionally impaired, with reduced cytotoxicity and IFN-γ secretion associated with the downregulation of activating receptors. They also showed a high degree of heterogeneity in mouse and human cSCCs with the expression of several markers of exhaustion, including TIGIT on NK cells and PD-1 and TIM-3 on ILC1s. Our data show an enrichment in inflammatory ILC1 at the precancerous stage together with impaired antitumor functions in NK cells and ILC1 that could contribute to the development of cSCC and thus suggest that future immunotherapies should take both ILC populations into account.

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.

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.

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.

A Single-Cell Atlas of the Human Healthy Airways.

Rationale: The respiratory tract constitutes an elaborate line of defense that is based on a unique cellular ecosystem.Objectives: We aimed to investigate cell population distributions and transcriptional changes along the airways by using single-cell RNA profiling.Methods: We have explored the cellular heterogeneity of the human airway epithelium in 10 healthy living volunteers by single-cell RNA profiling. A total of 77,969 cells were collected at 35 distinct locations, from the nose to the 12th division of the airway tree.Measurements and Main Results: The resulting atlas is composed of a high percentage of epithelial cells (89.1%) but also immune (6.2%) and stromal (4.7%) cells with distinct cellular proportions in different regions of the airways. It reveals differential gene expression between identical cell types (suprabasal, secretory, and multiciliated cells) from the nose (MUC4, PI3, SIX3) and tracheobronchial (SCGB1A1, TFF3) airways. By contrast, cell-type-specific gene expression is stable across all tracheobronchial samples. Our atlas improves the description of ionocytes, pulmonary neuroendocrine cells, and brush cells and identifies a related population of NREP-positive cells. We also report the association of KRT13 with dividing cells that are reminiscent of previously described mouse "hillock" cells and with squamous cells expressing SCEL and SPRR1A/B.Conclusions: Robust characterization of a single-cell cohort in healthy airways establishes a valuable resource for future investigations. The precise description of the continuum existing from the nasal epithelium to successive divisions of the airways and the stable gene expression profile of these regions better defines conditions under which relevant tracheobronchial proxies of human respiratory diseases can be developed.

Tumor-Associated Neutrophils Dampen Adaptive Immunity and Promote Cutaneous Squamous Cell Carcinoma Development.

Cutaneous squamous cell carcinoma (cSCC) development has been linked to immune dysfunctions but the mechanisms are still unclear. Here, we report a progressive infiltration of tumor-associated neutrophils (TANs) in precancerous and established cSCC lesions from chemically induced skin carcinogenesis. Comparative in-depth gene expression analyses identified a predominant protumor gene expression signature of TANs in lesions compared to their respective surrounding skin. In addition, in vivo depletion of neutrophils delayed tumor growth and significantly increased the frequency of proliferating IFN-γ (interferon-γ)-producing CD8+ T cells. Mechanisms that limited antitumor responses involved high arginase activity, production of reactive oxygen species (ROS) and nitrite (NO), and the expression of programmed death-ligand 1 (PD-L1) on TAN, concomitantly with an induction of PD-1 on CD8+ T cells, which correlated with tumor size. Our data highlight the relevance of targeting neutrophils and PD-L1-PD-1 (programmed death-1) interaction in the treatment of cSCC.

Regeneration of airway epithelial cells to study rare cell states in cystic fibrosis.

Pathological remodeling of the airway epithelium is commonly observed in cystic fibrosis (CF). Thus, tissue repair is critical to restore integrity and maintenance of the epithelial barrier function. Epithelial repair is a multi-step process initiated by progenitor cell migration into the injured area, proliferation, and re-differentiation into all of the cell types that contribute to the function of a normal airway epithelium. Recent technological advances applied to relevant animal and cell injury models have helped in understanding the complexity of progenitor cell differentiation. This short review will introduce the current knowledge of the mechanisms regulating airway epithelial cell (AEC) regeneration and repair, with a focus on the specification of two rare cell types/states: ionocytes and deuterosomal cells.