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 among media used for ALI-cultured human airway epithelial cells, the aim of our study was to evaluate the impact of several media (BEGM, PneumaCult, 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 PneumaCult-ALI and Half & Half, stronger EGF signaling from basal cells in BEGM-ALI, differential expression of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) entry factor ACE2, and distinct secretome transcripts depending on the media used. We also established that proliferation in PneumaCult-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 aid in choosing 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.

Novel dynamics of human mucociliary differentiation revealed by single-cell RNA sequencing of nasal epithelial cultures.

The upper airway epithelium, which is mainly composed of multiciliated, goblet, club and basal cells, ensures proper mucociliary function and can regenerate in response to assaults. In chronic airway diseases, defective repair leads to tissue remodeling. Delineating key drivers of differentiation dynamics can help understand how normal or pathological regeneration occurs. Using single-cell transcriptomics and lineage inference, we have unraveled trajectories from basal to luminal cells, providing novel markers for specific populations. We report that: (1) a precursor subgroup of multiciliated cells, which we have entitled deuterosomal cells, is defined by specific markers, such as DEUP1, FOXN4, YPEL1, HES6 and CDC20B; (2) goblet cells can be precursors of multiciliated cells, thus explaining the presence of hybrid cells that co-express markers of goblet and multiciliated cells; and (3) a repertoire of molecules involved in the regeneration process, such as keratins or components of the Notch, Wnt or BMP/TGFß pathways, can be identified. Confirmation of our results on fresh human and pig airway samples, and on mouse tracheal cells, extend and confirm our conclusions regarding the molecular and cellular choreography at work during mucociliary epithelial differentiation.

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.

CD4+ T Cells Affect the Thyroid Hormone Transport at the Choroid Plexus in Mice Raised in Enriched Environment.

BACKGROUND: Others and we have shown that T cells have an important role in hippocampal synaptic plasticity, including neurogenesis in the dentate gyrus, spinogenesis, and glutamatergic synaptic function in the CA of the hippocampus. Hippocampus plasticity is particularly involved in the brain effects of the enriched environment (EE), and interestingly CD4+ and CD8+ T cells play essential and differential roles in these effects. However, the precise mechanisms by which they act on the brain remain elusive. OBJECTIVES: We searched for a putative mechanism of action by which CD4+ T cells could influence brain plasticity and hypothesized that they could regulate protein transport at the level of the blood-CSF barrier in the choroid plexus. METHOD: We compared mice housed in EE and deprived of CD4+ T cells using a depleting antibody with a control group injected with the control isotype. We analyzed in the hippocampus the gene expression profiles using the Agilent system, and the expression of target proteins in plasma, CSF, and the choroid plexus using ELISA. RESULTS: We show that CD4+ T cells may influence EE-induced hippocampus plasticity via thyroid hormone signaling by regulating in the choroid plexus the expression of transthyretin, the major transporter of thyroxine (T4) to the brain parenchyma. CONCLUSIONS: Our study highlights the contribution of close interactions between the immune and neuroendocrine systems in brain plasticity and function.

A cost effective 5΄ selective single cell transcriptome profiling approach with improved UMI design.

Single cell RNA sequencing approaches are instrumental in studies of cell-to-cell variability. 5΄ selective transcriptome profiling approaches allow simultaneous definition of the transcription start size and have advantages over 3΄ selective approaches which just provide internal sequences close to the 3΄ end. The only currently existing 5΄ selective approach requires costly and labor intensive fragmentation and cell barcoding after cDNA amplification. We developed an optimized 5΄ selective workflow where all the cell indexing is done prior to fragmentation. With our protocol, cell indexing can be performed in the Fluidigm C1 microfluidic device, resulting in a significant reduction of cost and labor. We also designed optimized unique molecular identifiers that show less sequence bias and vulnerability towards sequencing errors resulting in an improved accuracy of molecule counting. We provide comprehensive experimental workflows for Illumina and Ion Proton sequencers that allow single cell sequencing in a cost range comparable to qPCR assays.

Comparative Genomic Analysis of Drechmeria coniospora Reveals Core and Specific Genetic Requirements for Fungal Endoparasitism of Nematodes.

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.

CD8+ T cells are essential for the effects of enriched environment on hippocampus-dependent behavior, hippocampal neurogenesis and synaptic plasticity.

Enriched environment (EE) induces plasticity changes in the brain. Recently, CD4+ T cells have been shown to be involved in brain plasticity processes. Here, we show that CD8+ T cells are required for EE-induced brain plasticity in mice, as revealed by measurements of hippocampal volume, neurogenesis in the DG of the hippocampus, spinogenesis and glutamatergic synaptic function in the CA of the hippocampus. As a consequence, EE-induced behavioral benefits depend, at least in part, on CD8+ T cells. In addition, we show that spleen CD8+ T cells from mice housed in standard environment (SE) and EE have different properties in terms of 1) TNFα release after in vitro CD3/CD28 or PMA/Iono stimulation 2) in vitro proliferation properties 3) CD8+ CD44+ CD62Llow and CD62Lhi T cells repartition 4) transcriptomic signature as revealed by RNA sequencing. CD8+ T cells purified from the choroid plexus of SE and EE mice also exhibit different transcriptomic profiles as highlighted by single-cell mRNA sequencing. We show that CD8+ T cells are essential mediators of beneficial EE effects on brain plasticity and cognition. Additionally, we propose that EE differentially primes CD8+ T cells leading to behavioral improvement.

Identification of novel target genes for safer and more specific control of root-knot nematodes from a pan-genome mining.

Root-knot nematodes are globally the most aggressive and damaging plant-parasitic nematodes. Chemical nematicides have so far constituted the most efficient control measures against these agricultural pests. Because of their toxicity for the environment and danger for human health, these nematicides have now been banned from use. Consequently, new and more specific control means, safe for the environment and human health, are urgently needed to avoid worldwide proliferation of these devastating plant-parasites. Mining the genomes of root-knot nematodes through an evolutionary and comparative genomics approach, we identified and analyzed 15,952 nematode genes conserved in genomes of plant-damaging species but absent from non target genomes of chordates, plants, annelids, insect pollinators and mollusks. Functional annotation of the corresponding proteins revealed a relative abundance of putative transcription factors in this parasite-specific set compared to whole proteomes of root-knot nematodes. This may point to important and specific regulators of genes involved in parasitism. Because these nematodes are known to secrete effector proteins in planta, essential for parasitism, we searched and identified 993 such effector-like proteins absent from non-target species. Aiming at identifying novel targets for the development of future control methods, we biologically tested the effect of inactivation of the corresponding genes through RNA interference. A total of 15 novel effector-like proteins and one putative transcription factor compatible with the design of siRNAs were present as non-redundant genes and had transcriptional support in the model root-knot nematode Meloidogyne incognita. Infestation assays with siRNA-treated M. incognita on tomato plants showed significant and reproducible reduction of the infestation for 12 of the 16 tested genes compared to control nematodes. These 12 novel genes, showing efficient reduction of parasitism when silenced, constitute promising targets for the development of more specific and safer control means.

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.

An integrated cell atlas of the lung in health and disease.

Single-cell technologies have transformed our understanding of human tissues. Yet, studies typically capture only a limited number of donors and disagree on cell type definitions. Integrating many single-cell datasets can address these limitations of individual studies and capture the variability present in the population. Here we present the integrated Human Lung Cell Atlas (HLCA), combining 49 datasets of the human respiratory system into a single atlas spanning over 2.4 million cells from 486 individuals. The HLCA presents a consensus cell type re-annotation with matching marker genes, including annotations of rare and previously undescribed cell types. Leveraging the number and diversity of individuals in the HLCA, we identify gene modules that are associated with demographic covariates such as age, sex and body mass index, as well as gene modules changing expression along the proximal-to-distal axis of the bronchial tree. Mapping new data to the HLCA enables rapid data annotation and interpretation. Using the HLCA as a reference for the study of disease, we identify shared cell states across multiple lung diseases, including SPP1+ profibrotic monocyte-derived macrophages in COVID-19, pulmonary fibrosis and lung carcinoma. Overall, the HLCA serves as an example for the development and use of large-scale, cross-dataset organ atlases within the Human Cell Atlas.

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

Genes regulated in MPTP-treated macaques and human Parkinson's disease suggest a common signature in prefrontal cortex.

The presymptomatic phase of Parkinson's disease (PD) is now recognized as a prodromal phase, with compensatory mechanism masking its progression and non-motor early manifestations, such as depression, cognitive disturbances and apathy. Those mechanisms were thought to be strictly dopamine-mediated until recent advances have shed light upon involvement of putative outside-basal ganglia, i.e. cortical, structures. We took advantage of our progressive 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP)-treated macaque model to monitor whole genome transcriptional changes in several brain areas. Our data reveals that transcriptomic activity changes take place from early stages, suggesting very early compensatory mechanisms or pathological activity outside the basal ganglia, including the PFC. Specific transcriptomic changes occurring in the PFC of fully parkinsonian MPTP-treated macaques have been identified. Interestingly, a large part of these transcriptomic changes were also observed in human post-mortem samples of patients with neurodegenerative diseases analysed by quantitative PCR. These results suggest that the PFC is able to detect the progression of dopamine denervation even at very early time points. There are therefore mechanisms, within the PFC, leading to compensatory alterations and/or participating to pathophysiology of prodromal PD manifestations.

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.

Population genomics supports clonal reproduction and multiple independent gains and losses of parasitic abilities in the most devastating nematode pest.

The root-knot nematodes are the most devastating worms to worldwide agriculture with Meloidogyne incognita being the most widely distributed and damaging species. This parasitic and ecological success seems surprising given its supposed obligatory clonal reproduction. Clonal reproduction has been suspected based on cytological observations but, so far, never confirmed by population genomics data. As a species, M. incognita is highly polyphagous with thousands of host plants. However, different M. incognita isolates present distinct and overlapping patterns of host compatibilities. Historically, four "host races" had been defined as a function of ranges of compatible and incompatible plants. In this study, we used population genomics to assess whether (a) reproduction is actually clonal in this species, (b) the host races follow an underlying phylogenetic signal or, rather represent multiple independent transitions, and (c) how genome variations associate with other important biological traits such as the affected crops and geographical distribution. We sequenced the genomes of 11 M. incognita isolates across Brazil that covered the four host races in replicates. By aligning the genomic reads of these isolates to the M. incognita reference genome assembly, we identified point variations. Analysis of linkage disequilibrium and 4-gametes test showed no evidence for recombination, corroborating the clonal reproduction of M. incognita. The few point variations between the isolates showed no significant association with the host races, the geographical origin of the samples, or the crop on which they have been collected. Addition of isolates from other locations around the world confirmed this lack of underlying phylogenetic signal. This suggests multiple gains and losses of parasitic abilities and adaptations to different environments account for the broad host spectrum and wide geographical distribution of M. incognita and thus to its high economic impact. This surprising adaptability without sex poses both evolutionary and agro-economic challenges.

Primate-specific spliced PMCHL RNAs are non-protein coding in human and macaque tissues.

BACKGROUND: Brain-expressed genes that were created in primate lineage represent obvious candidates to investigate molecular mechanisms that contributed to neural reorganization and emergence of new behavioural functions in Homo sapiens. PMCHL1 arose from retroposition of a pro-melanin-concentrating hormone (PMCH) antisense mRNA on the ancestral human chromosome 5p14 when platyrrhines and catarrhines diverged. Mutations before divergence of hylobatidae led to creation of new exons and finally PMCHL1 duplicated in an ancestor of hominids to generate PMCHL2 at the human chromosome 5q13. A complex pattern of spliced and unspliced PMCHL RNAs were found in human brain and testis. RESULTS: Several novel spliced PMCHL transcripts have been characterized in human testis and fetal brain, identifying an additional exon and novel splice sites. Sequencing of PMCHL genes in several non-human primates allowed to carry out phylogenetic analyses revealing that the initial retroposition event took place within an intron of the brain cadherin (CDH12) gene, soon after platyrrhine/catarrhine divergence, i.e. 30-35 Mya, and was concomitant with the insertion of an AluSg element. Sequence analysis of the spliced PMCHL transcripts identified only short ORFs of less than 300 bp, with low (VMCH-p8 and protein variants) or no evolutionary conservation. Western blot analyses of human and macaque tissues expressing PMCHL RNA failed to reveal any protein corresponding to VMCH-p8 and protein variants encoded by spliced transcripts. CONCLUSION: Our present results improve our knowledge of the gene structure and the evolutionary history of the primate-specific chimeric PMCHL genes. These genes produce multiple spliced transcripts, bearing short, non-conserved and apparently non-translated ORFs that may function as mRNA-like non-coding RNAs.

CD4+ T Cells Have a Permissive Effect on Enriched Environment-Induced Hippocampus Synaptic Plasticity.

Living in an enriched environment (EE) benefits health by acting synergistically on various biological systems including the immune and the central nervous systems. The dialog between the brain and the immune cells has recently gained interest and is thought to play a pivotal role in beneficial effects of EE. Recent studies show that T lymphocytes have an important role in hippocampal plasticity, learning, and memory, although the precise mechanisms by which they act on the brain remain elusive. Using a mouse model of EE, we show here that CD4+ T cells are essential for spinogenesis and glutamatergic synaptic function in the CA of the hippocampus. However, CD4+ lymphocytes do not influence EE-induced neurogenesis in the DG of the hippocampus, by contrast to what we previously demonstrated for CD8+ T cells. Importantly, CD4+ T cells located in the choroid plexus have a specific transcriptomic signature as a function of the living environment. Our study highlights the contribution of CD4+ T cells in the brain plasticity and function.

Elucidation of the unexplored biodiversity of ant venom peptidomes via MALDI-TOF mass spectrometry and its application for chemotaxonomy.

The rise of integrative taxonomy, a multi-criteria approach used in characterizing species, fosters the development of new tools facilitating species delimitation. Mass spectrometric (MS) analysis of venom peptides from venomous animals has previously been demonstrated to be a valid method for identifying species. Here we aimed to develop a rapid chemotaxonomic tool for identifying ants based on venom peptide mass fingerprinting. The study focused on the biodiversity of ponerine ants (Hymenoptera: Formicidae: Ponerinae) in French Guiana. Initial experiments optimized the use of automated matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS) to determine variations in the mass profiles of ant venoms using several MALDI matrices and additives. Data were then analyzed via a hierarchical cluster analysis to classify the venoms of 17 ant species. In addition, phylogenetic relationships were assessed and were highly correlated with methods using DNA sequencing of the mitochondrial gene cytochrome c oxidase subunit 1. By combining a molecular genetics approach with this chemotaxonomic approach, we were able to improve the accuracy of the taxonomic findings to reveal cryptic ant species within species complexes. This chemotaxonomic tool can therefore contribute to more rapid species identification and more accurate taxonomies. BIOLOGICAL SIGNIFICANCE: This is the first extensive study concerning the peptide analysis of the venom of both Pachycondyla and Odontomachus ants. We studied the venoms of 17 ant species from French Guiana that permitted us to fine-tune the venom analysis of ponerine ants via MALDI-TOF mass spectrometry. We explored the peptidomes of crude ant venom and demonstrated that venom peptides can be used in the identification of ant species. In addition, the application of this novel chemotaxonomic method combined with a parallel genetic approach using COI sequencing permitted us to reveal the presence of cryptic ants within both the Pachycondyla apicalis and Pachycondyla stigma species complexes. This adds a new dimension to the search for means of exploiting the enormous biodiversity of venomous ants as a source for novel therapeutic drugs or biopesticides. This article is part of a Special Issue entitled: Proteomics of non-model organisms.

siRNAs Trigger Efficient Silencing of a Parasitism Gene in Plant Parasitic Root-Knot Nematodes.

Expanding genomic data on plant pathogens open new perspectives for the development of specific and environment friendly pest management strategies based on the inhibition of parasitism genes that are essential for the success of infection. Identifying such genes relies on accurate reverse genetics tools and the screening of pathogen knock-down phenotypes. Root-knot nematodes are major cosmopolitan crop pests that feed on a wide range of host plants. Small interfering RNAs (siRNAs) would provide a powerful tool for reverse genetics of nematode parasitism genes provided that they could (1) target genes expressed in inner tissues of infective nematodes and (2) target genes expressed during parasitism. In this study, we show that siRNAs can access inner tissues of the infective juveniles during soaking and accumulate in the esophagus, amphidial pouches and related neurons of the nematode. We provide evidence that siRNAs can trigger knock-down of the parasitism gene Mi-CRT, a calreticulin gene expressed in the esophageal glands of Meloidogyne incognita. Mi-CRT knock-down in infective juveniles affected nematode virulence. However, Mi-CRT knock-down was not persistent after plant infection, indicating that siRNA-mediated RNAi is best suited for functional analysis of genes involved in pre-parasitic stages or in the early steps of infection.