Search

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.

Abstract The variant 20I/501Y.V1, associated to a higher risk of transmissibility, emerged in Nice city (Southeast of France, French Riviera) during January 2021. The pandemic has resumed late December 2020 in this area. A high incidence rate together with a fast turn-over of the main circulating variants, provided us the opportunity to analyze modifications in clinical profile and outcome traits. We performed an observational study in the University hospital of Nice from December 2020 to February 2021. We analyzed data of sequencing of SARS-CoV-2 from the sewage collector and PCR screening from all positive samples at the hospital. Then, we described the characteristics of all COVID-19 patients admitted in the emergency department (ED) (n = 1247) and those hospitalized in the infectious diseases ward or ICU (n = 232). The UK-variant was absent in this area in December, then increasingly spread in January representing 59% of the PCR screening performed mid-February. The rate of patients over 65 years admitted to the ED decreased from 63 to 50% (p = 0.001). The mean age of hospitalized patients in the infectious diseases ward decreased from 70.7 to 59.2 (p

Adipose tissue is composed of a number of adipocytes and a number of other cells including immune cells. Here the authors use single-cell sequencing of murine brown adipose tissue immune cells and describe multiple macrophage and monocyte subsets and show that monocytes contribute to brown adipose tissue expansion.

Summary: 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). : Publisher’s note: Undertaking any experimental protocol requires adherence to local institutional guidelines for laboratory safety and ethics.

Summary: 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.

Droplet-based high throughput single cell sequencing techniques can often lose information on transcript splicing and heterogenity. Here the authors introduce ScNaUmi-seq, which uses Oxford Nanopore sequencing and barcoding to generate high accuracy full length sequences.

Abstract Current demand for SARS‐CoV‐2 testing is straining material resource and labor capacity around the globe. As a result, the public health and clinical community are hindered in their ability to monitor and contain the spread of COVID‐19. Despite broad consensus that more testing is needed, pragmatic guidance toward realizing this objective has been limited. This paper addresses this limitation by proposing a novel and geographically agnostic framework (the 4Ps framework) to guide multidisciplinary, scalable, resource‐efficient, and achievable efforts toward enhanced testing capacity. The 4Ps (Prioritize, Propagate, Partition, and Provide) are described in terms of specific opportunities to enhance the volume, diversity, characterization, and implementation of SARS‐CoV‐2 testing to benefit public health. Coordinated deployment of the strategic and tactical recommendations described in this framework has the potential to rapidly expand available testing capacity, improve public health decision‐making in response to the COVID‐19 pandemic, and/or to be applied in future emergent disease outbreaks.

The emergence and quick spread of SARS-CoV-2 has pointed at a low capacity response for testing large populations in many countries, in line of material, technical and staff limitations. The traditional RT-qPCR diagnostic test remains the reference method and is by far the most widely used test. These assays are limited to a few probe sets, require large sample PCR reaction volumes, along with an expensive and time-consuming RNA extraction step. Here we describe a quantitative nanofluidic assay that overcomes some of these shortcomings, based on the BiomarkTM instrument from Fluidigm. This system offers the possibility of performing 4608 qPCR end-points in a single run, equivalent to 192 clinical samples combined with 12 pairs of primers/probe sets in duplicate, thus allowing the monitoring of SARS-CoV-2 including the detection of specific SARS-CoV-2 variants, as well as the detection other pathogens and/or host cellular responses (virus receptors, response markers, microRNAs). The 10 nL-range volume of BiomarkTM reactions is compatible with sensitive and reproducible reactions that can be easily and cost-effectively adapted to various RT-qPCR configurations and sets of primers/probe. Finally, we also evaluated the use of inactivating lysis buffers composed of various detergents in the presence or absence of proteinase K to assess the compatibility of these buffers with a direct reverse transcription enzymatic step and we propose several protocols, bypassing the need for RNA purification. We advocate that the combined utilization of an optimized processing buffer and a high-throughput real-time PCR device would contribute to improve the turn-around-time to deliver the test results to patients and increase the SARS-CoV-2 testing capacities.

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.

In vertebrate multiciliated cells, most centrioles are synthesized and released by poorly characterized structures called deuterosomes. Here, the authors report that the previously-uncharacterized Cdc20b in mouse and Xenopus associates with deuterosomes and contributes to centriole release.

MicroRNAs have been best characterized for their functions in the cytoplasm; however, there is growing evidence of a nuclear localized role. Here, the authors identify Sfpq as an Ago2-interacting protein that modulates miRNA activity in both the nucleus and cytoplasm.

Summary: Breast cancer stem cells (bCSCs) have been implicated in tumor progression and therapeutic resistance; however, the molecular mechanisms that define this state are unclear. We have performed two microRNA (miRNA) gain- and loss-of-function screens to identify miRNAs that regulate the choice between bCSC self-renewal and differentiation. We find that micro-RNA (miR)-600 silencing results in bCSC expansion, while its overexpression reduces bCSC self-renewal, leading to decreased in vivo tumorigenicity. miR-600 targets stearoyl desaturase 1 (SCD1), an enzyme required to produce active, lipid-modified WNT proteins. In the absence of miR-600, WNT signaling is active and promotes self-renewal, whereas overexpression of miR-600 inhibits the production of active WNT and promotes bCSC differentiation. In a series of 120 breast tumors, we found that a low level of miR-600 is correlated with active WNT signaling and a poor prognosis. These findings highlight a miR-600-centered signaling network that governs bCSC-fate decisions and influences tumor progression. : El Helou et al. identify miRNAs that are able to balance bCSC fate. They find that miR-600 silencing results in bCSC expansion, while its overexpression reduces bCSC self-renewal. miR-600 was further found to regulate WNT signaling through SCD1, and miR-600 expression correlates with clinical outcome.

The ribosome profiling technique (Ribo-seq) allows the selective sequencing of translated RNA regions. Recently, the analysis of genomic sequences associated to Ribo-seq reads has been widely employed to assess their coding potential. These analyses led to the identification of differentially translated transcripts under different experimental conditions, and/or ribosome pausing on codon motifs. In the context of the ever-growing need for tools analyzing Ribo-seq reads, we have developed ‘RiboProfiling’, a new Bioconductor open-source package. ‘RiboProfiling’ provides a full pipeline to cover all key steps for the analysis of ribosome footprints. This pipeline has been implemented in a single R workflow. The package takes an alignment (BAM) file as input and performs ribosome footprint quantification at a transcript level. It also identifies footprint accumulation on particular amino acids or multi amino-acids motifs. Report summary graphs and data quantification are generated automatically. The package facilitates quality assessment and quantification of Ribo-seq experiments. Its implementation in Bioconductor enables the modeling and statistical analysis of its output through the vast choice of packages available in R. This article illustrates how to identify codon-motifs accumulating ribosome footprints, based on data from Escherichia coli.

Drechmeria coniospora is an obligate fungal pathogen that infects nematodes via the adhesion of specialized spores to the host cuticle. D. coniospora is frequently found associated with Caenorhabditis elegans in environmental samples. It is used in the study of the nematode's response to fungal infection. Full understanding of this bi-partite interaction requires knowledge of the pathogen's genome, analysis of its gene expression program and a capacity for genetic engineering. The acquisition of all three is reported here. A phylogenetic analysis placed D. coniospora close to the truffle parasite Tolypocladium ophioglossoides, and Hirsutella minnesotensis, another nematophagous fungus. Ascomycete nematopathogenicity is polyphyletic; D. coniospora represents a branch that has not been molecularly characterized. A detailed in silico functional analysis, comparing D. coniospora to 11 fungal species, revealed genes and gene families potentially involved in virulence and showed it to be a highly specialized pathogen. A targeted comparison with nematophagous fungi highlighted D. coniospora-specific genes and a core set of genes associated with nematode parasitism. A comparative gene expression analysis of samples from fungal spores and mycelia, and infected C. elegans, gave a molecular view of the different stages of the D. coniospora lifecycle. Transformation of D. coniospora allowed targeted gene knock-out and the production of fungus that expresses fluorescent reporter genes. It also permitted the initial characterisation of a potential fungal counter-defensive strategy, involving interference with a host antimicrobial mechanism. This high-quality annotated genome for D. coniospora gives insights into the evolution and virulence of nematode-destroying fungi. Coupled with genetic transformation, it opens the way for molecular dissection of D. coniospora physiology, and will allow both sides of the interaction between D. coniospora and C. elegans, as well as the evolutionary arms race that exists between pathogen and host, to be studied.

The prevalence of obesity has steadily increased over the last few decades. During this time, populations of industrialized countries have been exposed to diets rich in fat with a high content of linoleic acid and a low content of α-linolenic acid compared with recommended intake. To assess the contribution of dietary fatty acids, male and female mice fed a high-fat diet (35% energy as fat, linoleic acid:α-linolenic acid ratio of 28) were mated randomly and maintained after breeding on the same diet for successive generations. Offspring showed, over four generations, a gradual enhancement in fat mass due to combined hyperplasia and hypertrophy with no change in food intake. Transgenerational alterations in adipokine levels were accompanied by hyperinsulinemia. Gene expression analyses of the stromal vascular fraction of adipose tissue, over generations, revealed discrete and steady changes in certain important players, such as CSF3 and Nocturnin. Thus, under conditions of genome stability and with no change in the regimen over four generations, we show that a Western-like fat diet induces a gradual fat mass enhancement, in accordance with the increasing prevalence of obesity observed in humans.

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.

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.

Specificity of interaction between a microRNA (miRNA) and its targets crucially depends on the seed region located in its 5'-end. It is often implicitly considered that two miRNAs sharing the same biological activity should display similarity beyond the strict six nucleotide region that forms the seed, in order to form specific complexes with the same mRNA targets. We have found that expression of hsa-miR-147b and hsa-miR-210, though triggered by different stimuli (i.e. lipopolysaccharides and hypoxia, respectively), induce very similar cellular effects in term of proliferation, migration and apoptosis. Hsa-miR-147b only shares a "minimal" 6-nucleotides seed sequence with hsa-miR-210, but is identical with hsa-miR-147a over 20 nucleotides, except for one base located in the seed region. Phenotypic changes induced after heterologous expression of miR-147a strikingly differ from those induced by miR-147b or miR-210. In particular, miR-147a behaves as a potent inhibitor of cell proliferation and migration. These data fit well with the gene expression profiles observed for miR-147b and miR-210, which are very similar, and the gene expression profile of miR-147a, which is distinct from the two others. Bioinformatics analysis of all human miRNA sequences indicates multiple cases of miRNAs from distinct families exhibiting the same kind of similarity that would need to be further characterized in terms of putative functional redundancy. Besides, it implies that functional impact of some miRNAs can be masked by robust expression of miRNAs belonging to distinct families.

C. burnetii is a Gram-negative intracellular Y-proteobacteria that causes the zoonotic disease Q fever. Q fever can manifest as an acute or chronic illness. Different typing methods have been previously developed to classify C. burnetii isolates to explore its pathogenicity. Here, we report a comprehensive genomotyping method based on the presence or absence of genes using microarrays. The genomotyping method was then tested in 52 isolates obtained from different geographic areas, different hosts and patients with different clinical manifestations. The analysis revealed the presence of 10 genomotypes organized into 3 groups, with a topology congruent with that obtained through multi-spacer typing. We also found that only 4 genomotypes were specifically associated with acute Q fever, whereas all of the genomotypes could be associated to chronic human infection. Serendipitously, the genomotyping results revealed that all hard tick isolates, including the Nine Mile strain, belong to the same genomotype.

Coxiella burnetii, the causative agent of the zoonotic disease Q fever, is mainly transmitted to humans through an aerosol route. A spore-like form allows C. burnetii to resist different environmental conditions. Because of this, analysis of the survival strategies used by this bacterium to adapt to new environmental conditions is critical for our understanding of C. burnetii pathogenicity. Here, we report the early transcriptional response of C. burnetii under temperature stresses. Our data show that C. burnetii exhibited minor changes in gene regulation under short exposure to heat or cold shock. While small differences were observed, C. burnetii seemed to respond similarly to cold and heat shock. The expression profiles obtained using microarrays produced in-house were confirmed by quantitative RT-PCR. Under temperature stresses, 190 genes were differentially expressed in at least one condition, with a fold change of up to 4. Globally, the differentially expressed genes in C. burnetii were associated with bacterial division, (p)ppGpp synthesis, wall and membrane biogenesis and, especially, lipopolysaccharide and peptidoglycan synthesis. These findings could be associated with growth arrest and witnessed transformation of the bacteria to a spore-like form. Unexpectedly, clusters of neighboring genes were differentially expressed. These clusters do not belong to operons or genetic networks; they have no evident associated functions and are not under the control of the same promoters. We also found undescribed but comparable clusters of regulation in previously reported transcriptomic analyses of intracellular bacteria, including Rickettsia sp. and Listeria monocytogenes. The transcriptomic patterns of C. burnetii observed under temperature stresses permits the recognition of unpredicted clusters of regulation for which the trigger mechanism remains unidentified but which may be the result of a new mechanism of epigenetic regulation.

BackgroundEpithelial-mesenchymal interactions are critical in regulating many aspects of vertebrate embryo development, and for the maintenance of homeostatic equilibrium in adult tissues. The interactions between epithelium and mesenchyme are believed to be mediated by paracrine signals such as cytokines and extracellular matrix components secreted from fibroblasts that affect adjacent epithelia. In this study, we sought to identify the repertoire of microRNAs (miRNAs) in normal lung human fibroblasts and their potential regulation by the cytokines TNF-alpha, IL-1beta and TGF-beta.Methodology/principal findingsMiR-155 was significantly induced by inflammatory cytokines TNF-alpha and IL-1beta while it was down-regulated by TGF-beta. Ectopic expression of miR-155 in human fibroblasts induced modulation of a large set of genes related to "cell to cell signalling", "cell morphology" and "cellular movement". This was consistent with an induction of caspase-3 activity and with an increase in cell migration in fibroblasts tranfected with miR-155. Using different miRNA bioinformatic target prediction tools, we found a specific enrichment for miR-155 predicted targets among the population of down-regulated transcripts. Among fibroblast-selective targets, one interesting hit was keratinocyte growth factor (KGF, FGF-7), a member of the fibroblast growth factor (FGF) family, which owns two potential binding sites for miR-155 in its 3'-UTR. Luciferase assays experimentally validated that miR-155 can efficiently target KGF 3'-UTR. Site-directed mutagenesis revealed that only one out of the 2 potential sites was truly functional. Functional in vitro assays experimentally validated that miR-155 can efficiently target KGF 3'-UTR. Furthermore, in vivo experiments using a mouse model of lung fibrosis showed that miR-155 expression level was correlated with the degree of lung fibrosis.Conclusions/significanceOur results strongly suggest a physiological function of miR-155 in lung fibroblasts. Altogether, this study implicates this miRNA in the regulation by mesenchymal cells of surrounding lung epithelium, making it a potential key player during tissue injury.

The ribosome profiling technique (Ribo-seq) allows the selective sequencing of translated RNA regions. Recently, the analysis of genomic sequences associated to Ribo-seq reads has been widely employed to assess their coding potential. These analyses led to the identification of differentially translated transcripts under different experimental conditions, and/or ribosome pausing on codon motifs. In the context of the ever-growing need for tools analyzing Ribo-seq reads, we have developed ‘RiboProfiling’, a new Bioconductor open-source package. ‘RiboProfiling’ provides a full pipeline to cover all key steps for the analysis of ribosome footprints. This pipeline has been implemented in a single R workflow. The package takes an alignment (BAM) file as input and performs ribosome footprint quantification at a transcript level. It also identifies footprint accumulation on particular amino acids or multi amino-acids motifs. Report summary graphs and data quantification are generated automatically. The package facilitates quality assessment and quantification of Ribo-seq experiments. Its implementation in Bioconductor enables the modeling and statistical analysis of its output through the vast choice of packages available in R. This article illustrates how to identify codon-motifs accumulating ribosome footprints, based on data from Escherichia coli.