Deep sequencing of Phox2a nuclei reveals five classes of anterolateral system neurons.

The anterolateral system (ALS) is a major ascending pathway from the spinal cord that projects to multiple brain areas and underlies the perception of pain, itch, and skin temperature. Despite its importance, our understanding of this system has been hampered by the considerable functional and molecular diversity of its constituent cells. Here, we use fluorescence-activated cell sorting to isolate ALS neurons belonging to the Phox2a-lineage for single-nucleus RNA sequencing. We reveal five distinct clusters of ALS neurons (ALS1-5) and document their laminar distribution in the spinal cord using in situ hybridization. We identify three clusters of neurons located predominantly in laminae I-III of the dorsal horn (ALS1-3) and two clusters with cell bodies located in deeper laminae (ALS4 and ALS5). Our findings reveal the transcriptional logic that underlies ALS neuronal diversity in the adult mouse and uncover the molecular identity of two previously identified classes of projection neurons. We also show that these molecular signatures can be used to target groups of ALS neurons using retrograde viral tracing. Overall, our findings provide a valuable resource for studying somatosensory biology and targeting subclasses of ALS neurons.

Single-cell nanobiopsy enables multigenerational longitudinal transcriptomics of cancer cells.

Single-cell RNA sequencing has revolutionized our understanding of cellular heterogeneity, but routine methods require cell lysis and fail to probe the dynamic trajectories responsible for cellular state transitions, which can only be inferred. Here, we present a nanobiopsy platform that enables the injection of exogenous molecules and multigenerational longitudinal cytoplasmic sampling from a single cell and its progeny. The technique is based on scanning ion conductance microscopy (SICM) and, as a proof of concept, was applied to longitudinally profile the transcriptome of single glioblastoma (GBM) brain tumor cells in vitro over 72 hours. The GBM cells were biopsied before and after exposure to chemotherapy and radiotherapy, and our results suggest that treatment either induces or selects for more transcriptionally stable cells. We envision the nanobiopsy will contribute to transforming standard single-cell transcriptomics from a static analysis into a dynamic assay.

Automated Purification of DNA Origami with SPRI Beads.

DNA origami synthesis is a well-established technique with wide-ranging applications. In most cases, the synthesized origami must be purified to remove excess materials such as DNA oligos and other functional molecules. While several purification techniques are routinely used, all have limitations, and cannot be integrated with robotic systems. Here the use of solid-phase reversible immobilization (SPRI) beads as a scalable, high-throughput, and automatable method to purify DNA origami is demonstrated. Not only can this method remove unreacted oligos and biomolecules with yields comparable to existing methods while maintaining the high structural integrity of the origami, but it can also be integrated into an automated workflow to purify simultaneously large numbers and quantities of samples. It is envisioned that the SPRI beads purification method will improve the scalability of DNA nanostructures synthesis both for research and commercial applications.

Absence of the primary cilia formation gene Talpid3 impairs muscle stem cell function.

Skeletal muscle stem cells (MuSC) are crucial for tissue homoeostasis and repair after injury. Following activation, they proliferate to generate differentiating myoblasts. A proportion of cells self-renew, re-enter the MuSC niche under the basal lamina outside the myofiber and become quiescent. Quiescent MuSC have a primary cilium, which is disassembled upon cell cycle entry. Ex vivo experiments suggest cilia are important for MuSC self-renewal, however, their requirement for muscle regeneration in vivo remains poorly understood. Talpid3 (TA3) is essential for primary cilia formation and Hedgehog (Hh) signalling. Here we use tamoxifen-inducible conditional deletion of TA3 in MuSC (iSC-KO) and show that regeneration is impaired in response to cytotoxic injury. Depletion of MuSC after regeneration suggests impaired self-renewal, also consistent with an exacerbated phenotype in TA3iSC-KO mice after repeat injury. Single cell transcriptomics of MuSC progeny isolated from myofibers identifies components of several signalling pathways, which are deregulated in absence of TA3, including Hh and Wnt. Pharmacological activation of Wnt restores muscle regeneration, while purmorphamine, an activator of the Smoothened (Smo) co-receptor in the Hh pathway, has no effect. Together, our data show that TA3 and primary cilia are important for MuSC self-renewal and pharmacological treatment can efficiently restore muscle regeneration.

Deep sequencing of Phox2a nuclei reveals five classes of anterolateral system neurons.

The anterolateral system (ALS) is a major ascending pathway from the spinal cord that projects to multiple brain areas and underlies the perception of pain, itch and skin temperature. Despite its importance, our understanding of this system has been hampered by the considerable functional and molecular diversity of its constituent cells. Here we use fluorescence-activated cell sorting to isolate ALS neurons belonging to the Phox2a-lineage for single-nucleus RNA sequencing. We reveal five distinct clusters of ALS neurons (ALS1-5) and document their laminar distribution in the spinal cord using in situ hybridization. We identify 3 clusters of neurons located predominantly in laminae I-III of the dorsal horn (ALS1-3) and two clusters with cell bodies located in deeper laminae (ALS4 & ALS5). Our findings reveal the transcriptional logic that underlies ALS neuronal diversity in the adult mouse and uncover the molecular identity of two previously identified classes of projection neurons. We also show that these molecular signatures can be used to target groups of ALS neurons using retrograde viral tracing. Overall, our findings provide a valuable resource for studying somatosensory biology and targeting subclasses of ALS neurons.

Identification of a non-canonical ciliate nuclear genetic code where UAA and UAG code for different amino acids.

The genetic code is one of the most highly conserved features across life. Only a few lineages have deviated from the "universal" genetic code. Amongst the few variants of the genetic code reported to date, the codons UAA and UAG virtually always have the same translation, suggesting that their evolution is coupled. Here, we report the genome and transcriptome sequencing of a novel uncultured ciliate, belonging to the Oligohymenophorea class, where the translation of the UAA and UAG stop codons have changed to specify different amino acids. Genomic and transcriptomic analyses revealed that UAA has been reassigned to encode lysine, while UAG has been reassigned to encode glutamic acid. We identified multiple suppressor tRNA genes with anticodons complementary to the reassigned codons. We show that the retained UGA stop codon is enriched in the 3'UTR immediately downstream of the coding region of genes, suggesting that there is functional drive to maintain tandem stop codons. Using a phylogenomics approach, we reconstructed the ciliate phylogeny and mapped genetic code changes, highlighting the remarkable number of independent genetic code changes within the Ciliophora group of protists. According to our knowledge, this is the first report of a genetic code variant where UAA and UAG encode different amino acids.

Perivascular niche cells sense thrombocytopenia and activate hematopoietic stem cells in an IL-1 dependent manner.

Hematopoietic stem cells (HSCs) residing in specialized niches in the bone marrow are responsible for the balanced output of multiple short-lived blood cell lineages in steady-state and in response to different challenges. However, feedback mechanisms by which HSCs, through their niches, sense acute losses of specific blood cell lineages remain to be established. While all HSCs replenish platelets, previous studies have shown that a large fraction of HSCs are molecularly primed for the megakaryocyte-platelet lineage and are rapidly recruited into proliferation upon platelet depletion. Platelets normally turnover in an activation-dependent manner, herein mimicked by antibodies inducing platelet activation and depletion. Antibody-mediated platelet activation upregulates expression of Interleukin-1 (IL-1) in platelets, and in bone marrow extracellular fluid in vivo. Genetic experiments demonstrate that rather than IL-1 directly activating HSCs, activation of bone marrow Lepr+ perivascular niche cells expressing IL-1 receptor is critical for the optimal activation of quiescent HSCs upon platelet activation and depletion. These findings identify a feedback mechanism by which activation-induced depletion of a mature blood cell lineage leads to a niche-dependent activation of HSCs to reinstate its homeostasis.

Panhematopoietic RNA barcoding enables kinetic measurements of nucleate and anucleate lineages and the activation of myeloid clones following acute platelet depletion.

BACKGROUND: Platelets and erythrocytes constitute over 95% of all hematopoietic stem cell output. However, the clonal dynamics of HSC contribution to these lineages remains largely unexplored. RESULTS: We use lentiviral genetic labeling of mouse hematopoietic stem cells to quantify output from all lineages, nucleate, and anucleate, simultaneously linking these with stem and progenitor cell transcriptomic phenotypes using single-cell RNA-sequencing. We observe dynamic shifts of clonal behaviors through time in same-animal peripheral blood and demonstrate that acute platelet depletion shifts the output of multipotent hematopoietic stem cells to the exclusive production of platelets. Additionally, we observe the emergence of new myeloid-biased clones, which support short- and long-term production of blood cells. CONCLUSIONS: Our approach enables kinetic studies of multi-lineage output in the peripheral blood and transcriptional heterogeneity of individual hematopoietic stem cells. Our results give a unique insight into hematopoietic stem cell reactivation upon platelet depletion and of clonal dynamics in both steady state and under stress.

Single-cell gene and isoform expression analysis reveals signatures of ageing in haematopoietic stem and progenitor cells.

Single-cell approaches have revealed that the haematopoietic hierarchy is a continuum of differentiation, from stem cell to committed progenitor, marked by changes in gene expression. However, many of these approaches neglect isoform-level information and thus do not capture the extent of alternative splicing within the system. Here, we present an integrated short- and long-read single-cell RNA-seq analysis of haematopoietic stem and progenitor cells. We demonstrate that over half of genes detected in standard short-read single-cell analyses are expressed as multiple, often functionally distinct, isoforms, including many transcription factors and key cytokine receptors. We observe global and HSC-specific changes in gene expression with ageing but limited impact of ageing on isoform usage. Integrating single-cell and cell-type-specific isoform landscape in haematopoiesis thus provides a new reference for comprehensive molecular profiling of heterogeneous tissues, as well as novel insights into transcriptional complexity, cell-type-specific splicing events and consequences of ageing.

Computational prediction and experimental validation of Salmonella Typhimurium SopE-mediated fine-tuning of autophagy in intestinal epithelial cells.

Macroautophagy is a ubiquitous homeostasis and health-promoting recycling process of eukaryotic cells, targeting misfolded proteins, damaged organelles and intracellular infectious agents. Some intracellular pathogens such as Salmonella enterica serovar Typhimurium hijack this process during pathogenesis. Here we investigate potential protein-protein interactions between host transcription factors and secreted effector proteins of Salmonella and their effect on host gene transcription. A systems-level analysis identified Salmonella effector proteins that had the potential to affect core autophagy gene regulation. The effect of a SPI-1 effector protein, SopE, that was predicted to interact with regulatory proteins of the autophagy process, was investigated to validate our approach. We then confirmed experimentally that SopE can directly bind to SP1, a host transcription factor, which modulates the expression of the autophagy gene MAP1LC3B. We also revealed that SopE might have a double role in the modulation of autophagy: Following initial increase of MAP1LC3B transcription triggered by Salmonella infection, subsequent decrease in MAP1LC3B transcription at 6h post-infection was SopE-dependent. SopE also played a role in modulation of the autophagy flux machinery, in particular MAP1LC3B and p62 autophagy proteins, depending on the level of autophagy already taking place. Upon typical infection of epithelial cells, the autophagic flux is increased. However, when autophagy was chemically induced prior to infection, SopE dampened the autophagic flux. The same was also observed when most of the intracellular Salmonella cells were not associated with the SCV (strain lacking sifA) regardless of the autophagy induction status before infection. We demonstrated how regulatory network analysis can be used to better characterise the impact of pathogenic effector proteins, in this case, Salmonella. This study complements previous work in which we had demonstrated that specific pathogen effectors can affect the autophagy process through direct interaction with autophagy proteins. Here we show that effector proteins can also influence the upstream regulation of the process. Such interdisciplinary studies can increase our understanding of the infection process and point out targets important in intestinal epithelial cell defense.

Single-cell genomics reveals population structures from in vitro evolutionary studies of Salmonella.

Single-cell DNA sequencing has the potential to reveal detailed hierarchical structures in evolving populations of cells. Single cell approaches are increasingly used to study clonal evolution in human ageing and cancer but have not yet been deployed to study evolving clonal microbial populations. Here, we present an approach for single bacterial genomic analysis for in vitro evolution experiments using FACS isolation of individual bacteria followed by whole-genome amplification and sequencing. We apply this to the experimental evolution of a hypermutator strain of Salmonella in response to antibiotic stress (ciprofloxacin). By analysing sequence polymorphisms in individual cells from populations we identified the presence and prevalence of sub-populations which have acquired polymorphisms in genes previously demonstrated to be associated with ciprofloxacin susceptibility. We were also able to identify that the population exposed to antibiotic stress was able to develop resistance whilst maintaining diversity. This population structure could not be resolved from bulk sequence data, and our results show how high-throughput single-cell sequencing can enhance experimental studies of bacterial evolution.

Into the multiverse: advances in single-cell multiomic profiling.

Single-cell transcriptomic approaches have revolutionised the study of complex biological systems, with the routine measurement of gene expression in thousands of cells enabling construction of whole-organism cell atlases. However, the transcriptome is just one layer amongst many that coordinate to define cell type and state and, ultimately, function. In parallel with the widespread uptake of single-cell RNA-seq (scRNA-seq), there has been a rapid emergence of methods that enable multiomic profiling of individual cells, enabling parallel measurement of intercellular heterogeneity in the genome, epigenome, transcriptome, and proteomes. Linking measurements from each of these layers has the potential to reveal regulatory and functional mechanisms underlying cell behaviour in healthy development and disease.

High expression of the MADS-box gene VRT2 increases the number of rudimentary basal spikelets in wheat.

Spikelets are the fundamental building blocks of Poaceae inflorescences, and their development and branching patterns determine the various inflorescence architectures and grain yield of grasses. In wheat (Triticum aestivum), the central spikelets produce the most and largest grains, while spikelet size gradually decreases acropetally and basipetally, giving rise to the characteristic lanceolate shape of wheat spikes. The acropetal gradient corresponds with the developmental age of spikelets; however, the basal spikelets are developed first, and the cause of their small size and rudimentary development is unclear. Here, we adapted G&T-seq, a low-input transcriptomics approach, to characterize gene expression profiles within spatial sections of individual spikes before and after the establishment of the lanceolate shape. We observed larger differences in gene expression profiles between the apical, central, and basal sections of a single spike than between any section belonging to consecutive developmental time points. We found that SHORT VEGETATIVE PHASE MADS-box transcription factors, including VEGETATIVE TO REPRODUCTIVE TRANSITION 2 (VRT-A2), are expressed highest in the basal section of the wheat spike and display the opposite expression gradient to flowering E-class SEPALLATA1 genes. Based on multi-year field trials and transgenic lines, we show that higher expression of VRT-A2 in the basal sections of the spike is associated with increased numbers of rudimentary basal spikelets. Our results, supported by computational modeling, suggest that the delayed transition of basal spikelets from vegetative to floral developmental programs results in the lanceolate shape of wheat spikes. This study highlights the value of spatially resolved transcriptomics to gain insights into developmental genetics pathways of grass inflorescences.

Free fatty-acid transport via CD36 drives ß-oxidation-mediated hematopoietic stem cell response to infection.

Acute infection is known to induce rapid expansion of hematopoietic stem cells (HSCs), but the mechanisms supporting this expansion remain incomplete. Using mouse models, we show that inducible CD36 is required for free fatty acid uptake by HSCs during acute infection, allowing the metabolic transition from glycolysis towards ß-oxidation. Mechanistically, high CD36 levels promote FFA uptake, which enables CPT1A to transport fatty acyl chains from the cytosol into the mitochondria. Without CD36-mediated FFA uptake, the HSCs are unable to enter the cell cycle, subsequently enhancing mortality in response to bacterial infection. These findings enhance our understanding of HSC metabolism in the bone marrow microenvironment, which supports the expansion of HSCs during pathogenic challenge.

Chemotherapy induces canalization of cell state in childhood B-cell precursor acute lymphoblastic leukemia.

Comparison of intratumor genetic heterogeneity in cancer at diagnosis and relapse suggests that chemotherapy induces bottleneck selection of subclonal genotypes. However, evolutionary events subsequent to chemotherapy could also explain changes in clonal dominance seen at relapse. We, therefore, investigated the mechanisms of selection in childhood B-cell precursor acute lymphoblastic leukemia (BCP-ALL) during induction chemotherapy where maximal cytoreduction occurs. To distinguish stochastic versus deterministic events, individual leukemias were transplanted into multiple xenografts and chemotherapy administered. Analyses of the immediate post-treatment leukemic residuum at single-cell resolution revealed that chemotherapy has little impact on genetic heterogeneity. Rather, it acts on extensive, previously unappreciated, transcriptional and epigenetic heterogeneity in BCP-ALL, dramatically reducing the spectrum of cell states represented, leaving a genetically polyclonal but phenotypically uniform population with hallmark signatures relating to developmental stage, cell cycle and metabolism. Hence, canalization of cell state accounts for a significant component of bottleneck selection during induction chemotherapy.

Meiosis and beyond - understanding the mechanistic and evolutionary processes shaping the germline genome.

The separation of germ cell populations from the soma is part of the evolutionary transition to multicellularity. Only genetic information present in the germ cells will be inherited by future generations, and any molecular processes affecting the germline genome are therefore likely to be passed on. Despite its prevalence across taxonomic kingdoms, we are only starting to understand details of the underlying micro-evolutionary processes occurring at the germline genome level. These include segregation, recombination, mutation and selection and can occur at any stage during germline differentiation and mitotic germline proliferation to meiosis and post-meiotic gamete maturation. Selection acting on germ cells at any stage from the diploid germ cell to the haploid gametes may cause significant deviations from Mendelian inheritance and may be more widespread than previously assumed. The mechanisms that affect and potentially alter the genomic sequence and allele frequencies in the germline are pivotal to our understanding of heritability. With the rise of new sequencing technologies, we are now able to address some of these unanswered questions. In this review, we comment on the most recent developments in this field and identify current gaps in our knowledge.

Common clonal origin of conventional T cells and induced regulatory T cells in breast cancer patients.

Regulatory CD4+ T cells (Treg) prevent tumor clearance by conventional T cells (Tconv) comprising a major obstacle of cancer immune-surveillance. Hitherto, the mechanisms of Treg repertoire formation in human cancers remain largely unclear. Here, we analyze Treg clonal origin in breast cancer patients using T-Cell Receptor and single-cell transcriptome sequencing. While Treg in peripheral blood and breast tumors are clonally distinct, Tconv clones, including tumor-antigen reactive effectors (Teff), are detected in both compartments. Tumor-infiltrating CD4+ cells accumulate into distinct transcriptome clusters, including early activated Tconv, uncommitted Teff, Th1 Teff, suppressive Treg and pro-tumorigenic Treg. Trajectory analysis suggests early activated Tconv differentiation either into Th1 Teff or into suppressive and pro-tumorigenic Treg. Importantly, Tconv, activated Tconv and Treg share highly-expanded clones contributing up to 65% of intratumoral Treg. Here we show that Treg in human breast cancer may considerably stem from antigen-experienced Tconv converting into secondary induced Treg through intratumoral activation.

Transcriptional characterization of human megakaryocyte polyploidization and lineage commitment.

BACKGROUND: Megakaryocytes (MKs) originate from cells immuno-phenotypically indistinguishable from hematopoietic stem cells (HSCs), bypassing intermediate progenitors. They mature within the adult bone marrow and release platelets into the circulation. Until now, there have been no transcriptional studies of primary human bone marrow MKs. OBJECTIVES: To characterize MKs and HSCs from human bone marrow using single-cell RNA sequencing, to investigate MK lineage commitment, maturation steps, and thrombopoiesis. RESULTS: We show that MKs at different levels of polyploidization exhibit distinct transcriptional states. Although high levels of platelet-specific gene expression occur in the lower ploidy classes, as polyploidization increases, gene expression is redirected toward translation and posttranslational processing transcriptional programs, in preparation for thrombopoiesis. Our findings are in keeping with studies of MK ultrastructure and supersede evidence generated using in vitro cultured MKs. Additionally, by analyzing transcriptional signatures of a single HSC, we identify two MK-biased HSC subpopulations exhibiting unique differentiation kinetics. We show that human bone marrow MKs originate from these HSC subpopulations, supporting the notion that they display priming for MK differentiation. Finally, to investigate transcriptional changes in MKs associated with stress thrombopoiesis, we analyzed bone marrow MKs from individuals with recent myocardial infarction and found a specific gene expression signature. Our data support the modulation of MK differentiation in this thrombotic state. CONCLUSIONS: Here, we use single-cell sequencing for the first time to characterize the human bone marrow MK transcriptome at different levels of polyploidization and investigate their differentiation from the HSC.

A low-cost pipeline for soil microbiome profiling.

Common bottlenecks in environmental and crop microbiome studies are the consumable and personnel costs necessary for genomic DNA extraction and sequencing library construction. This is harder for challenging environmental samples such as soil, which is rich in Polymerase Chain Reaction (PCR) inhibitors. To address this, we have established a low-cost genomic DNA extraction method for soil samples. We also present an Illumina-compatible 16S and ITS rRNA gene amplicon library preparation workflow that uses common laboratory equipment. We evaluated the performance of our genomic DNA extraction method against two leading commercial soil genomic DNA kits (MoBio PowerSoil® and MP Biomedicals™ FastDNA™ SPIN) and a recently published non-commercial extraction method by Zou et al. (PLoS Biology, 15, e2003916, 2017). Our benchmarking experiment used four different soil types (coniferous, broad-leafed, and mixed forest plus a standardized cereal crop compost mix) assessing the quality and quantity of the extracted genomic DNA by analyzing sequence variants of 16S V4 and ITS rRNA amplicons. We found that our genomic DNA extraction method compares well to both commercially available genomic DNA extraction kits in DNA quality and quantity. The MoBio PowerSoil® kit, which relies on silica column-based DNA extraction with extensive washing, delivered the cleanest genomic DNA, for example, best A260:A280 and A260:A230 absorbance ratios. The MP Biomedicals™ FastDNA™ SPIN kit, which uses a large amount of binding material, yielded the most genomic DNA. Our method fits between the two commercial kits, producing both good yields and clean genomic DNA with fragment sizes of approximately 10 kb. Comparative analysis of detected amplicon sequence variants shows that our method correlates well with the two commercial kits. Here, we present a low-cost genomic DNA extraction method for soil samples that can be coupled to an Illumina-compatible simple two-step amplicon library construction workflow for 16S V4 and ITS marker genes. Our method delivers high-quality genomic DNA at a fraction of the cost of commercial kits and enables cost-effective, large-scale amplicon sequencing projects. Notably, our extracted gDNA molecules are long enough to be suitable for downstream techniques such as full gene sequencing or even metagenomics shotgun approaches using long reads (PacBio or Nanopore), 10x Genomics linked reads, and Dovetail genomics.

Non-neuronal expression of SARS-CoV-2 entry genes in the olfactory system suggests mechanisms underlying COVID-19-associated anosmia.

Altered olfactory function is a common symptom of COVID-19, but its etiology is unknown. A key question is whether SARS-CoV-2 (CoV-2) - the causal agent in COVID-19 - affects olfaction directly, by infecting olfactory sensory neurons or their targets in the olfactory bulb, or indirectly, through perturbation of supporting cells. Here we identify cell types in the olfactory epithelium and olfactory bulb that express SARS-CoV-2 cell entry molecules. Bulk sequencing demonstrated that mouse, non-human primate and human olfactory mucosa expresses two key genes involved in CoV-2 entry, ACE2 and TMPRSS2. However, single cell sequencing revealed that ACE2 is expressed in support cells, stem cells, and perivascular cells, rather than in neurons. Immunostaining confirmed these results and revealed pervasive expression of ACE2 protein in dorsally-located olfactory epithelial sustentacular cells and olfactory bulb pericytes in the mouse. These findings suggest that CoV-2 infection of non-neuronal cell types leads to anosmia and related disturbances in odor perception in COVID-19 patients.