New technologies for improved relevance in miRNA research.

After a number of years of research in the field of miRNA, the robustness and biological relevance of many published articles is increasingly being questioned. We propose the use of new RNA-seq approaches, genome editing technologies, and updated public databases to improve the quality, reliability, and relevance of published data.

CD206+ macrophages are relevant non-invasive imaging biomarkers and therapeutic targets in experimental lung fibrosis.

BACKGROUND: Interstitial lung diseases (ILDs) include a large number of diseases associated with progressive pulmonary fibrosis (PPF), including idiopathic pulmonary fibrosis (IPF). Despite the rarity of each of the fibrotic ILDs individually, they cumulatively affect a considerable number of patients. PPF is characterised by an excessive collagen deposition leading to functional decline. OBJECTIVES: Therapeutic options are limited to nintedanib and pirfenidone which are only able to reduce fibrosis progression. CD206-expressing M2 macrophages are involved in fibrosis progression, and whether they may be relevant therapeutic targets or biomarkers remains an open question. RESULTS: In our study, CD206+ lung macrophages were monitored in bleomycin-induced lung fibrosis in mice by combining flow cytometry, scRNAseq and in vivo molecular imaging using a single photon emission computed tomography (SPECT) radiopharmaceutical, 99mTc-tilmanocept. The antifibrotic effect of the inhibition of M2 macrophage polarisation with a JAK inhibitor, tofacitinib, was assessed in vivo. We demonstrate that CD206-targeted in vivo SPECT imaging with 99mTc-tilmanocept was able to accurately detect and quantify the increase in CD206+ macrophages from early to advanced stages of experimental fibrosis and ex vivo in lung biopsies from patients with IPF. CD206-targeted imaging also specifically detected a decrease in CD206+ lung macrophages on nintedanib and tofacitinib treatment. Importantly, early in vivo imaging of CD206+ macrophages allowed the prediction of experimental lung fibrosis progression as well as nintedanib and tofacitinib efficacy. CONCLUSIONS: These findings indicate that M2 macrophages may be relevant theranostic targets for personalised medicine for patients with PPF.

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.

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 critical role for miR-142 in alveolar epithelial lineage formation in mouse lung development.

The respiratory epithelium arises from alveolar epithelial progenitors which differentiate into alveolar epithelial type 1 (AT1) and type 2 (AT2) cells. AT2 cells are stem cells in the lung critical for the repair process after injury. Mechanisms regulating AT1 and AT2 cell maturation are poorly defined. We report that the activation of the glucocorticoid pathway in an in vitro alveolar epithelial lineage differentiation assay led to increased AT2 marker Sftpc and decreased miR-142 expression. Using miR-142 KO mice, we demonstrate an increase in the AT2/AT1 cell number ratio. Overexpression of miR-142 in alveolar progenitor cells in vivo led to the opposite effect. Examination of the KO lungs at E18.5 revealed enhanced expression of miR-142 targets Apc, Ep300 and Kras associated with increased ß-catenin and p-Erk signaling. Silencing of miR-142 expression in lung explants grown in vitro triggers enhanced Sftpc expression as well as increased AT2/AT1 cell number ratio. Pharmacological inhibition of Ep300-ß-catenin but not Erk in vitro prevented the increase in Sftpc expression triggered by loss of miR-142. These results suggest that the glucocorticoid-miR-142-Ep300-ß-catenin signaling axis controls pneumocyte maturation.

The Long Noncoding RNA DNM3OS Is a Reservoir of FibromiRs with Major Functions in Lung Fibroblast Response to TGF-ß and Pulmonary Fibrosis.

Rationale: Given the paucity of effective treatments for idiopathic pulmonary fibrosis (IPF), new insights into the deleterious mechanisms controlling lung fibroblast activation, the key cell type driving the fibrogenic process, are essential to develop new therapeutic strategies. TGF-ß (transforming growth factor-ß) is the main profibrotic factor, but its inhibition is associated with severe side effects because of its pleiotropic role. Objectives: To determine if downstream noncoding effectors of TGF-ß in fibroblasts may represent new effective therapeutic targets whose modulation may be well tolerated. Methods: We investigated the whole noncoding fraction of TGF-ß-stimulated lung fibroblast transcriptome to identify new genomic determinants of lung fibroblast differentiation into myofibroblasts. Differential expression of the long noncoding RNA (lncRNA) DNM3OS (dynamin 3 opposite strand) and its associated microRNAs (miRNAs) was validated in a murine model of pulmonary fibrosis and in IPF tissue samples. Distinct and complementary antisense oligonucleotide-based strategies aiming at interfering with DNM3OS were used to elucidate the role of DNM3OS and its associated miRNAs in IPF pathogenesis. Measurements and Main Results: We identified DNM3OS as a fibroblast-specific critical downstream effector of TGF-ß-induced lung myofibroblast activation. Mechanistically, DNM3OS regulates this process in trans by giving rise to three distinct profibrotic mature miRNAs (i.e., miR-199a-5p/3p and miR-214-3p), which influence SMAD and non-SMAD components of TGF-ß signaling in a multifaceted way. In vivo, we showed that interfering with DNM3OS function not only prevents lung fibrosis but also improves established pulmonary fibrosis. Conclusions: Pharmacological approaches aiming at interfering with the lncRNA DNM3OS may represent new effective therapeutic strategies in IPF.

Rapid decay of engulfed extracellular miRNA by XRN1 exonuclease promotes transient epithelial-mesenchymal transition.

Extracellular vesicles (EVs) have been shown to play an important role in intercellular communication as carriers of DNA, RNA and proteins. While the intercellular transfer of miRNA through EVs has been extensively studied, the stability of extracellular miRNA (ex-miRNA) once engulfed by a recipient cell remains to be determined. Here, we identify the ex-miRNA-directed phenotype to be transient due to the rapid decay of ex-miRNA. We demonstrate that the ex-miR-223-3p transferred from polymorphonuclear leukocytes to cancer cells were functional, as demonstrated by the decreased expression of its target FOXO1 and the occurrence of epithelial-mesenchymal transition reprogramming. We showed that the engulfed ex-miRNA, unlike endogenous miRNA, was unstable, enabling dynamic regulation and a return to a non-invasive phenotype within 8 h. This transient phenotype could be modulated by targeting XRN1/PACMAN exonuclease. Indeed, its silencing was associated with slower decay of ex-miR-223-3p and subsequently prolonged the invasive properties. In conclusion, we showed that the 'steady step' level of engulfed miRNA and its subsequent activity was dependent on the presence of a donor cell in the surroundings to constantly fuel the recipient cell with ex-miRNAs and of XRN1 exonuclease, which is involved in the decay of these imported miRNA.

A new long noncoding RNA (lncRNA) is induced in cutaneous squamous cell carcinoma and down-regulates several anticancer and cell differentiation genes in mouse.

Keratinocyte-derived cutaneous squamous cell carcinoma (cSCC) is the most common metastatic skin cancer. Although some of the early events involved in this pathology have been identified, the subsequent steps leading to tumor development are poorly defined. We demonstrate here that the development of mouse tumors induced by the concomitant application of a carcinogen and a tumor promoter (7,12-dimethylbenz[a]anthracene (DMBA) and 12-O-tetradecanoylphorbol-13-acetate (TPA), respectively) is associated with the up-regulation of a previously uncharacterized long noncoding RNA (lncRNA), termed AK144841. We found that AK144841 expression was absent from normal skin and was specifically stimulated in tumors and highly tumorigenic cells. We also found that AK144841 exists in two variants, one consisting of a large 2-kb transcript composed of four exons and one consisting of a 1.8-kb transcript lacking the second exon. Gain- and loss-of-function studies indicated that AK144841 mainly inhibited gene expression, specifically down-regulating the expression of genes of the late cornified envelope-1 (Lce1) family involved in epidermal terminal differentiation and of anticancer genes such as Cgref1, Brsk1, Basp1, Dusp5, Btg2, Anpep, Dhrs9, Stfa2, Tpm1, SerpinB2, Cpa4, Crct1, Cryab, Il24, Csf2, and Rgs16 Interestingly, the lack of the second exon significantly decreased AK144841's inhibitory effect on gene expression. We also noted that high AK144841 expression correlated with a low expression of the aforementioned genes and with the tumorigenic potential of cell lines. These findings suggest that AK144841 could contribute to the dedifferentiation program of tumor-forming keratinocytes and to molecular cascades leading to tumor development.

Tacrolimus-induced nephrotoxicity in mice is associated with microRNA deregulation.

Although Tacrolimus is an immunosuppressive drug widely used in renal transplantation, its chronic use paradoxically induces nephrotoxic effects, in particular renal fibrosis, which is responsible for chronic allograft dysfunction and represents a major prognostic factor of allograft survival. As molecular pathways and mechanisms involved in Tacrolimus-induced fibrogenic response are poorly elucidated, we assessed whether miRNAs are involved in the nephrotoxic effects mediated by Tacrolimus. Treatment of CD-1 mice with Tacrolimus (1 mg/kg/d for 28 days) resulted in kidney injury and was associated with alteration of a gene expression signature associated with cellular stress, fibrosis and inflammation. Tacrolimus also affected renal miRNA expression, including miRNAs previously involved in fibrotic and inflammatory processes as "fibromirs" such as miR-21-5p, miR-199a-5p and miR-214-3p. In agreement with in vivo data, Renal Proximal Tubular Epithelial cells exposed to Tacrolimus (25 and 50 µM) showed upregulation of miR-21-5p and the concomitant induction of epithelial phenotypic changes, inflammation and oxidative stress. In conclusion, this study suggests for the first time that miRNAs, especially fibromiRs, participate to Tacrolimus-induced nephrotoxic effects. Therefore, targeting miRNAs may be a new therapeutic option to counteract Tacrolimus deleterious effects on kidney.

MicroRNA-142 is a multifaceted regulator in organogenesis, homeostasis, and disease.

Over the past decade, microRNA-142 (miR-142) is emerging as a major regulator of cell fate decision in the hematopoietic system. However, miR-142 is expressed in many other tissues, and recent evidence suggests that it may play a more pleiotropic role during embryonic development. In addition, miR-142 has been shown to play important functions in disease. miR-142 displays a functional role in cancer, virus infection, inflammation, and immune tolerance. Both a guide strand (miR-142-3p) and passenger strand (miR-142-5p) are generated from the miR-142 hairpin. miR-142-3p and -5p display overlapping but also independent target genes. Loss of function mouse models (genetrap, global knock out [KO], and conditional KO) have been reported and support the important role of miR-142 in different biological processes. This review will summarize the abundant literature already available for miR-142 and will lay the foundation for future works on this important microRNA. Developmental Dynamics 246:285-290, 2017. © 2016 Wiley Periodicals, Inc.

A novel role for the RNA-binding protein FXR1P in myoblasts cell-cycle progression by modulating p21/Cdkn1a/Cip1/Waf1 mRNA stability.

The Fragile X-Related 1 gene (FXR1) is a paralog of the Fragile X Mental Retardation 1 gene (FMR1), whose absence causes the Fragile X syndrome, the most common form of inherited intellectual disability. FXR1P plays an important role in normal muscle development, and its absence causes muscular abnormalities in mice, frog, and zebrafish. Seven alternatively spliced FXR1 transcripts have been identified and two of them are skeletal muscle-specific. A reduction of these isoforms is found in myoblasts from Facio-Scapulo Humeral Dystrophy (FSHD) patients. FXR1P is an RNA-binding protein involved in translational control; however, so far, no mRNA target of FXR1P has been linked to the drastic muscular phenotypes caused by its absence. In this study, gene expression profiling of C2C12 myoblasts reveals that transcripts involved in cell cycle and muscular development pathways are modulated by Fxr1-depletion. We observed an increase of p21--a regulator of cell-cycle progression--in Fxr1-knocked-down mouse C2C12 and FSHD human myoblasts. Rescue of this molecular phenotype is possible by re-expressing human FXR1P in Fxr1-depleted C2C12 cells. FXR1P muscle-specific isoforms bind p21 mRNA via direct interaction with a conserved G-quadruplex located in its 3' untranslated region. The FXR1P/G-quadruplex complex reduces the half-life of p21 mRNA. In the absence of FXR1P, the upregulation of p21 mRNA determines the elevated level of its protein product that affects cell-cycle progression inducing a premature cell-cycle exit and generating a pool of cells blocked at G0. Our study describes a novel role of FXR1P that has crucial implications for the understanding of its role during myogenesis and muscle development, since we show here that in its absence a reduced number of myoblasts will be available for muscle formation/regeneration, shedding new light into the pathophysiology of FSHD.

miR-199a-5p Is upregulated during fibrogenic response to tissue injury and mediates TGFbeta-induced lung fibroblast activation by targeting caveolin-1.

As miRNAs are associated with normal cellular processes, deregulation of miRNAs is thought to play a causative role in many complex diseases. Nevertheless, the precise contribution of miRNAs in fibrotic lung diseases, especially the idiopathic form (IPF), remains poorly understood. Given the poor response rate of IPF patients to current therapy, new insights into the pathogenic mechanisms controlling lung fibroblasts activation, the key cell type driving the fibrogenic process, are essential to develop new therapeutic strategies for this devastating disease. To identify miRNAs with potential roles in lung fibrogenesis, we performed a genome-wide assessment of miRNA expression in lungs from two different mouse strains known for their distinct susceptibility to develop lung fibrosis after bleomycin exposure. This led to the identification of miR-199a-5p as the best miRNA candidate associated with bleomycin response. Importantly, miR-199a-5p pulmonary expression was also significantly increased in IPF patients (94 IPF versus 83 controls). In particular, levels of miR-199a-5p were selectively increased in myofibroblasts from injured mouse lungs and fibroblastic foci, a histologic feature associated with IPF. Therefore, miR-199a-5p profibrotic effects were further investigated in cultured lung fibroblasts: miR-199a-5p expression was induced upon TGFß exposure, and ectopic expression of miR-199a-5p was sufficient to promote the pathogenic activation of pulmonary fibroblasts including proliferation, migration, invasion, and differentiation into myofibroblasts. In addition, we demonstrated that miR-199a-5p is a key effector of TGFß signaling in lung fibroblasts by regulating CAV1, a critical mediator of pulmonary fibrosis. Remarkably, aberrant expression of miR-199a-5p was also found in unilateral ureteral obstruction mouse model of kidney fibrosis, as well as in both bile duct ligation and CCl4-induced mouse models of liver fibrosis, suggesting that dysregulation of miR-199a-5p represents a general mechanism contributing to the fibrotic process. MiR-199a-5p thus behaves as a major regulator of tissue fibrosis with therapeutic potency to treat fibroproliferative diseases.

Tetraspanin CD63 acts as a pro-metastatic factor via ß-catenin stabilization.

The tetraspanin CD63 is implicated in pro-metastatic signaling pathways but, so far, it is unclear, how CD63 levels affect the tumor cell phenotype. Here, we investigated the effect of CD63 modulation in different metastatic tumor cell lines. In vitro, knock down of CD63 induced a more epithelial-like phenotype concomitant with increased E-cadherin expression, downregulation of its repressors Slug and Zeb1, and decreased N-cadherin. In addition, ß-catenin protein was markedly reduced, negatively affecting expression of the target genes MMP-2 and PAI-1. ß-catenin inhibitors mimicked the epithelial phenotype induced by CD63 knock down. Inhibition of ß-catenin upstream regulators PI3K/AKT or GSK3ß could rescue the mesenchymal phenotype underlining the importance of the ß-catenin pathway in CD63-regulated cell plasticity. CD63 knock down-induced phenotypical changes correlated with a decrease of experimental metastasis whereas CD63 overexpression enhanced the tumor cell-intrinsic metastatic potential. Taken together, our data show that CD63 is a crucial player in the regulation of the tumor cell-intrinsic metastatic potential by affecting cell plasticity.

The 3' UTR of FMR1 mRNA is a target of miR-101, miR-129-5p and miR-221: implications for the molecular pathology of FXTAS at the synapse.

While FMR1 is silenced in Fragile X syndrome (FXS) patients carrying the full mutation, its expression is elevated (2-8 fold) in premutated individuals. These people may develop the Fragile X-associated Tremor/Ataxia syndrome (FXTAS), a late onset neurodegenerative disorder characterized by ataxia and parkinsonism. In addition, people carrying the premutation can be affected by a set of neurological and behavioral disorders during young age. Problems of memory have been detected in these patients as well as in the mouse models for FXTAS. To date little is known concerning the metabolism of FMR1 mRNA, notwithstanding the importance of the finely tuned regulation of the expression of this gene. In the present study, we identified three microRNAs that specifically target the 3' UTR of FMR1 and can modulate its expression throughout the brain particularly at the synapse where their expression is very high. The expression level of miR-221 is reduced in synaptosomal preparations of young FXTAS mice suggesting a general deregulation of transcripts located at the synapse of these mice. By transcriptome analysis, we show here a robust deregulation of the expression levels of genes involved in learning, memory and autistic behavior, Parkinson disease and neurodegeneration. These findings suggest the presence of a synaptopathy in these animals. Interestingly, many of those deregulated mRNAs are target of the same microRNAs that modulate the expression of FMR1 at the synapse.

miR-193b/365a cluster controls progression of epidermal squamous cell carcinoma.

Incidence of cutaneous squamous cell carcinomas (cSCCs) constantly increases in the Caucasian population. Developing preferentially on precancerous lesions such as actinic keratoses due to chronic sunlight exposure, cSCCs result from the malignant transformation of keratinocytes. Although a resection of the primary tumor is usually curative, a subset of aggressive cSCCs shows a high risk of recurrence and metastases. The characterization of the molecular dysfunctions involved in cSCC development should help to identify new relevant targets against these aggressive cSCCs. In that context, we have used small RNA sequencing to identify 100 microRNAs (miRNAs) whose expression was altered during chemically induced mouse skin tumorigenesis. The decreased expression of the miR-193b/365a cluster during tumor progression suggests a tumor suppressor role. Ectopic expression of these miRNAs in tumor cells indeed inhibited their proliferation, clonogenic potential and migration, which were stimulated in normal keratinocytes when these miRNAs were blocked with antisense oligonucleotides. A combination of in silico predictions and transcriptome analyses identified several target genes of interest. We validated KRAS and MAX as direct targets of miR-193b and miR-365a. Repression of these targets using siRNAs mimicked the effects of miR-193b and miR-365a, suggesting that these genes might mediate, at least in part, the tumor-suppressive action of these miRNAs.

Forkhead Box F1 represses cell growth and inhibits COL1 and ARPC2 expression in lung fibroblasts in vitro.

Aberrant expression of master phenotype regulators or alterations in their downstream pathways in lung fibroblasts may play a central role in idiopathic pulmonary fibrosis (IPF). Interrogating IPF fibroblast transcriptome datasets, we identified Forkhead Box F1 (FOXF1), a DNA-binding protein required for lung development, as a candidate actor in IPF. Thus we determined FOXF1 expression levels in fibroblasts cultured from normal or IPF lungs in vitro, and explored FOXF1 functions in these cells using transient and stable loss-of-function and gain-of-function models. FOXF1 mRNA and protein were expressed at higher levels in IPF fibroblasts compared with normal fibroblasts (mRNA: +44%, protein: +77%). Immunohistochemistry showed FOXF1 expression in nuclei of bronchial smooth muscle cells, endothelial cells, and lung fibroblasts including fibroblastic foci of IPF lungs. In normal lung fibroblasts, FOXF1 repressed cell growth and expression of collagen-1 (COL1) and actin-related protein 2/3 complex, subunit 2 (ARPC2). ARPC2 knockdown inhibited cell growth and COL1 expression, consistent with FOXF1 acting in part through ARPC2 repression. In IPF fibroblasts, COL1 and ARPC2 repression by FOXF1 was blunted, and FOXF1 did not repress growth. FOXF1 expression was induced by the antifibrotic mediator prostaglandin E2 and repressed by the profibrotic cytokine transforming growth factor-ß1 in both normal and IPF lung fibroblasts. Ex vivo, FOXF1 knockdown conferred CCL-210 lung fibroblasts the ability to implant in uninjured mouse lungs. In conclusion, FOXF1 functions and regulation were consistent with participation in antifibrotic pathways. Alterations of pathways downstream of FOXF1 may participate to fibrogenesis in IPF fibroblasts.

Dkk3 is a component of the genetic circuitry regulating aldosterone biosynthesis in the adrenal cortex.

Primary aldosteronism (PA, autonomous aldosterone production from the adrenal cortex) causes the most common form of secondary arterial hypertension (HT), which is also the most common curable form of HT. Recent studies have highlighted an important role of mutations in genes encoding potassium channels in the pathogenesis of PA, both in human disease and in animal models. Here, we have exploited the unique features of the hyperaldosteronemic phenotype of Kcnk3 null mice, which is dependent on sexual hormones, to identify genes whose expression is modulated in the adrenal gland according to the dynamic hyperaldosteronemic phenotype of those animals. Genetic inactivation of one of the genes identified by our strategy, dickkopf-3 (Dkk3), whose expression is increased by calcium influx into adrenocortical cells, in the Kcnk3 null background results in the extension of the low-renin, potassium-rich diet insensitive hyperaldosteronemic phenotype to the male sex. Compound Kcnk3/Dkk3 animals display an increased expression of Cyp11b2, the rate-limiting enzyme for aldosterone biosynthesis in the adrenal zona glomerulosa (ZG). Our data show that Dkk3 can act as a modifier gene in a mouse model for altered potassium channel function and suggest its potential involvement in human PA syndromes.

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.