Massively parallel in vivo Perturb-seq reveals cell-type-specific transcriptional networks in cortical development.

Leveraging AAVs' versatile tropism and labeling capacity, we expanded the scale of in vivo CRISPR screening with single-cell transcriptomic phenotyping across embryonic to adult brains and peripheral nervous systems. Through extensive tests of 86 vectors across AAV serotypes combined with a transposon system, we substantially amplified labeling efficacy and accelerated in vivo gene delivery from weeks to days. Our proof-of-principle in utero screen identified the pleiotropic effects of Foxg1, highlighting its tight regulation of distinct networks essential for cell fate specification of Layer 6 corticothalamic neurons. Notably, our platform can label >6% of cerebral cells, surpassing the current state-of-the-art efficacy at <0.1% by lentivirus, to achieve analysis of over 30,000 cells in one experiment and enable massively parallel in vivo Perturb-seq. Compatible with various phenotypic measurements (single-cell or spatial multi-omics), it presents a flexible approach to interrogate gene function across cell types in vivo, translating gene variants to their causal function.

Live birth of chimeric monkey with high contribution from embryonic stem cells.

A formal demonstration that mammalian pluripotent stem cells possess preimplantation embryonic cell-like (naive) pluripotency is the generation of chimeric animals through early embryo complementation with homologous cells. Whereas such naive pluripotency has been well demonstrated in rodents, poor chimerism has been achieved in other species including non-human primates due to the inability of the donor cells to match the developmental state of the host embryos. Here, we have systematically tested various culture conditions for establishing monkey naive embryonic stem cells and optimized the procedures for chimeric embryo culture. This approach generated an aborted fetus and a live chimeric monkey with high donor cell contribution. A stringent characterization pipeline demonstrated that donor cells efficiently (up to 90%) incorporated into various tissues (including the gonads and placenta) of the chimeric monkeys. Our results have major implications for the study of primate naive pluripotency and genetic engineering of non-human primates.

Proper acquisition of cell class identity in organoids allows definition of fate specification programs of the human cerebral cortex.

Realizing the full utility of brain organoids to study human development requires understanding whether organoids precisely replicate endogenous cellular and molecular events, particularly since acquisition of cell identity in organoids can be impaired by abnormal metabolic states. We present a comprehensive single-cell transcriptomic, epigenetic, and spatial atlas of human cortical organoid development, comprising over 610,000 cells, from generation of neural progenitors through production of differentiated neuronal and glial subtypes. We show that processes of cellular diversification correlate closely to endogenous ones, irrespective of metabolic state, empowering the use of this atlas to study human fate specification. We define longitudinal molecular trajectories of cortical cell types during organoid development, identify genes with predicted human-specific roles in lineage establishment, and uncover early transcriptional diversity of human callosal neurons. The findings validate this comprehensive atlas of human corticogenesis in vitro as a resource to prime investigation into the mechanisms of human cortical development.

An IGHG1 variant exhibits polarized prevalence and confers enhanced IgG1 antibody responses against life-threatening organisms.

Evolutionary pressures sculpt population genetics, whereas immune adaptation fortifies humans against life-threatening organisms. How the evolution of selective genetic variation in adaptive immune receptors orchestrates the adaptation of human populations to contextual perturbations remains elusive. Here, we show that the G396R coding variant within the human immunoglobulin G1 (IgG1) heavy chain presents a concentrated prevalence in Southeast Asian populations. We uncovered a 190-kb genomic linkage disequilibrium block peaked in close proximity to this variant, suggestive of potential Darwinian selection. This variant confers heightened immune resilience against various pathogens and viper toxins in mice. Mechanistic studies involving severe acute respiratory syndrome coronavirus 2 infection and vaccinated individuals reveal that this variant enhances pathogen-specific IgG1+ memory B cell activation and antibody production. This G396R variant may have arisen on a Neanderthal haplotype background. These findings underscore the importance of an IGHG1 variant in reinforcing IgG1 antibody responses against life-threatening organisms, unraveling the intricate interplay between human evolution and immune adaptation.

A pan-cancer single-cell transcriptional atlas of tumor infiltrating myeloid cells.

Tumor-infiltrating myeloid cells (TIMs) are key regulators in tumor progression, but the similarity and distinction of their fundamental properties across different tumors remain elusive. Here, by performing a pan-cancer analysis of single myeloid cells from 210 patients across 15 human cancer types, we identified distinct features of TIMs across cancer types. Mast cells in nasopharyngeal cancer were found to be associated with better prognosis and exhibited an anti-tumor phenotype with a high ratio of TNF+/VEGFA+ cells. Systematic comparison between cDC1- and cDC2-derived LAMP3+ cDCs revealed their differences in transcription factors and external stimulus. Additionally, pro-angiogenic tumor-associated macrophages (TAMs) were characterized with diverse markers across different cancer types, and the composition of TIMs appeared to be associated with certain features of somatic mutations and gene expressions. Our results provide a systematic view of the highly heterogeneous TIMs and suggest future avenues for rational, targeted immunotherapies.

Stress-Induced Metabolic Disorder in Peripheral CD4+ T Cells Leads to Anxiety-like Behavior.

Physical or mental stress leads to neuroplasticity in the brain and increases the risk of depression and anxiety. Stress exposure causes the dysfunction of peripheral T lymphocytes. However, the pathological role and underlying regulatory mechanism of peripheral T lymphocytes in mood disorders have not been well established. Here, we show that the lack of CD4+ T cells protects mice from stress-induced anxiety-like behavior. Physical stress-induced leukotriene B4 triggers severe mitochondrial fission in CD4+ T cells, which further leads to a variety of behavioral abnormalities including anxiety, depression, and social disorders. Metabolomic profiles and single-cell transcriptome reveal that CD4+ T cell-derived xanthine acts on oligodendrocytes in the left amygdala via adenosine receptor A1. Mitochondrial fission promotes the de novo synthesis of purine via interferon regulatory factor 1 accumulation in CD4+ T cells. Our study implicates a critical link between a purine metabolic disorder in CD4+ T cells and stress-driven anxiety-like behavior.

Optogenetic Editing Reveals the Hierarchical Organization of Learned Action Sequences.

The organization of action into sequences underlies complex behaviors that are essential for organismal survival and reproduction. Despite extensive studies of innate sequences in relation to central pattern generators, how learned action sequences are controlled and whether they are organized as a chain or a hierarchy remain largely unknown. By training mice to perform heterogeneous action sequences, we demonstrate that striatal direct and indirect pathways preferentially encode different behavioral levels of sequence structure. State-dependent closed-loop optogenetic stimulation of the striatal direct pathway can selectively insert a single action element into the sequence without disrupting the overall sequence length. Optogenetic manipulation of the striatal indirect pathway completely removes the ongoing subsequence while leaving the following subsequence to be executed with the appropriate timing and length. These results suggest that learned action sequences are not organized in a serial but rather a hierarchical structure that is distinctly controlled by basal ganglia pathways.

Genomic Analyses from Non-invasive Prenatal Testing Reveal Genetic Associations, Patterns of Viral Infections, and Chinese Population History.

We analyze whole-genome sequencing data from 141,431 Chinese women generated for non-invasive prenatal testing (NIPT). We use these data to characterize the population genetic structure and to investigate genetic associations with maternal and infectious traits. We show that the present day distribution of alleles is a function of both ancient migration and very recent population movements. We reveal novel phenotype-genotype associations, including several replicated associations with height and BMI, an association between maternal age and EMB, and between twin pregnancy and NRG1. Finally, we identify a unique pattern of circulating viral DNA in plasma with high prevalence of hepatitis B and other clinically relevant maternal infections. A GWAS for viral infections identifies an exceptionally strong association between integrated herpesvirus 6 and MOV10L1, which affects piwi-interacting RNA (piRNA) processing and PIWI protein function. These findings demonstrate the great value and potential of accumulating NIPT data for worldwide medical and genetic analyses.

Distinct Circuits for the Formation and Retrieval of an Imprinted Olfactory Memory.

Memories formed early in life are particularly stable and influential, representing privileged experiences that shape enduring behaviors. We show that exposing newly hatched C. elegans to pathogenic bacteria results in persistent aversion to those bacterial odors, whereas adult exposure generates only transient aversive memory. Long-lasting imprinted aversion has a critical period in the first larval stage and is specific to the experienced pathogen. Distinct groups of neurons are required during formation (AIB, RIM) and retrieval (AIY, RIA) of the imprinted memory. RIM synthesizes the neuromodulator tyramine, which is required in the L1 stage for learning. AIY memory retrieval neurons sense tyramine via the SER-2 receptor, which is essential for imprinted, but not for adult-learned, aversion. Odor responses in several neurons, most notably RIA, are altered in imprinted animals. These findings provide insight into neuronal substrates of different forms of memory, and lay a foundation for further understanding of early learning.

Complementary Contributions of Striatal Projection Pathways to Action Initiation and Execution.

The performance of an action relies on the initiation and execution of appropriate movement sequences. Two basal ganglia pathways have been classically hypothesized to regulate this process via opposing roles in movement facilitation and suppression. By using a series of state-dependent optogenetic manipulations, we dissected the contributions of each pathway and found that both the direct striatonigral pathway and the indirect striatopallidal pathway are necessary for smooth initiation and the execution of learned action sequences. Optogenetic inhibition or stimulation of each pathway before sequence initiation increased the latency for initiation: manipulations of the striatonigral pathway activity slowed action initiation, and those of the striatopallidal pathway aborted action initiation. The inhibition of each pathway after initiation also impaired ongoing execution. Furthermore, the subtle activation of striatonigral neurons sustained the performance of learned sequences, while striatopallidal manipulations aborted ongoing performance. These results suggest a supportive versus permissive model, where patterns of coordinated activity, rather than the relative amount of activity in these pathways, regulate movement initiation and execution.

Metastatic Latency and Immune Evasion through Autocrine Inhibition of WNT.

Metastasis frequently develops years after the removal of a primary tumor, from a minority of disseminated cancer cells that survived as latent entities through unknown mechanisms. We isolated latency competent cancer (LCC) cells from early stage human lung and breast carcinoma cell lines and defined the mechanisms that suppress outgrowth, support long-term survival, and maintain tumor-initiating potential in these cells during the latent metastasis stage. LCC cells show stem-cell-like characteristics and express SOX2 and SOX9 transcription factors, which are essential for their survival in host organs under immune surveillance and for metastatic outgrowth under permissive conditions. Through expression of the WNT inhibitor DKK1, LCC cells self-impose a slow-cycling state with broad downregulation of ULBP ligands for NK cells and evasion of NK-cell-mediated clearance. By expressing a Sox-dependent stem-like state and actively silencing WNT signaling, LCC cells can enter quiescence and evade innate immunity to remain latent for extended periods.

Serpins promote cancer cell survival and vascular co-option in brain metastasis.

Brain metastasis is an ominous complication of cancer, yet most cancer cells that infiltrate the brain die of unknown causes. Here, we identify plasmin from the reactive brain stroma as a defense against metastatic invasion, and plasminogen activator (PA) inhibitory serpins in cancer cells as a shield against this defense. Plasmin suppresses brain metastasis in two ways: by converting membrane-bound astrocytic FasL into a paracrine death signal for cancer cells, and by inactivating the axon pathfinding molecule L1CAM, which metastatic cells express for spreading along brain capillaries and for metastatic outgrowth. Brain metastatic cells from lung cancer and breast cancer express high levels of anti-PA serpins, including neuroserpin and serpin B2, to prevent plasmin generation and its metastasis-suppressive effects. By protecting cancer cells from death signals and fostering vascular co-option, anti-PA serpins provide a unifying mechanism for the initiation of brain metastasis in lung and breast cancers.

Brain-wide correspondence of neuronal epigenomics and distant projections.

Single-cell analyses parse the brain's billions of neurons into thousands of 'cell-type' clusters residing in different brain structures1. Many cell types mediate their functions through targeted long-distance projections allowing interactions between specific cell types. Here we used epi-retro-seq2 to link single-cell epigenomes and cell types to long-distance projections for 33,034 neurons dissected from 32 different regions projecting to 24 different targets (225 source-to-target combinations) across the whole mouse brain. We highlight uses of these data for interrogating principles relating projection types to transcriptomics and epigenomics, and for addressing hypotheses about cell types and connections related to genetics. We provide an overall synthesis with 926 statistical comparisons of discriminability of neurons projecting to each target for every source. We integrate this dataset into the larger BRAIN Initiative Cell Census Network atlas, composed of millions of neurons, to link projection cell types to consensus clusters. Integration with spatial transcriptomics further assigns projection-enriched clusters to smaller source regions than the original dissections. We exemplify this by presenting in-depth analyses of projection neurons from the hypothalamus, thalamus, hindbrain, amygdala and midbrain to provide insights into properties of those cell types, including differentially expressed genes, their associated cis-regulatory elements and transcription-factor-binding motifs, and neurotransmitter use.

PKM2 functions as a histidine kinase to phosphorylate PGAM1 and increase glycolysis shunts in cancer.

Phosphoglycerate mutase 1 (PGAM1) is a key node enzyme that diverts the metabolic reactions from glycolysis into its shunts to support macromolecule biosynthesis for rapid and sustainable cell proliferation. It is prevalent that PGAM1 activity is upregulated in various tumors; however, the underlying mechanism remains unclear. Here, we unveil that pyruvate kinase M2 (PKM2) moonlights as a histidine kinase in a phosphoenolpyruvate (PEP)-dependent manner to catalyze PGAM1 H11 phosphorylation, that is essential for PGAM1 activity. Moreover, monomeric and dimeric but not tetrameric PKM2 are efficient to phosphorylate and activate PGAM1. In response to epidermal growth factor signaling, Src-catalyzed PGAM1 Y119 phosphorylation is a prerequisite for PKM2 binding and the subsequent PGAM1 H11 phosphorylation, which constitutes a discrepancy between tumor and normal cells. A PGAM1-derived pY119-containing cell-permeable peptide or Y119 mutation disrupts the interaction of PGAM1 with PKM2 and PGAM1 H11 phosphorylation, dampening the glycolysis shunts and tumor growth. Together, these results identify a function of PKM2 as a histidine kinase, and illustrate the importance of enzyme crosstalk as a regulatory mode during metabolic reprogramming and tumorigenesis.

Selection of bone metastasis seeds by mesenchymal signals in the primary tumor stroma.

How organ-specific metastatic traits arise in primary tumors remains unknown. Here, we show a role of the breast tumor stroma in selecting cancer cells that are primed for metastasis in bone. Cancer-associated fibroblasts (CAFs) in triple-negative (TN) breast tumors skew heterogeneous cancer cell populations toward a predominance of clones that thrive on the CAF-derived factors CXCL12 and IGF1. Limiting concentrations of these factors select for cancer cells with high Src activity, a known clinical predictor of bone relapse and an enhancer of PI3K-Akt pathway activation by CXCL12 and IGF1. Carcinoma clones selected in this manner are primed for metastasis in the CXCL12-rich microenvironment of the bone marrow. The evidence suggests that stromal signals resembling those of a distant organ select for cancer cells that are primed for metastasis in that organ, thus illuminating the evolution of metastatic traits in a primary tumor and its distant metastases.

Neurotensin orchestrates valence assignment in the amygdala.

The ability to associate temporally segregated information and assign positive or negative valence to environmental cues is paramount for survival. Studies have shown that different projections from the basolateral amygdala (BLA) are potentiated following reward or punishment learning1-7. However, we do not yet understand how valence-specific information is routed to the BLA neurons with the appropriate downstream projections, nor do we understand how to reconcile the sub-second timescales of synaptic plasticity8-11 with the longer timescales separating the predictive cues from their outcomes. Here we demonstrate that neurotensin (NT)-expressing neurons in the paraventricular nucleus of the thalamus (PVT) projecting to the BLA (PVT-BLA:NT) mediate valence assignment by exerting NT concentration-dependent modulation in BLA during associative learning. We found that optogenetic activation of the PVT-BLA:NT projection promotes reward learning, whereas PVT-BLA projection-specific knockout of the NT gene (Nts) augments punishment learning. Using genetically encoded calcium and NT sensors, we further revealed that both calcium dynamics within the PVT-BLA:NT projection and NT concentrations in the BLA are enhanced after reward learning and reduced after punishment learning. Finally, we showed that CRISPR-mediated knockout of the Nts gene in the PVT-BLA pathway blunts BLA neural dynamics and attenuates the preference for active behavioural strategies to reward and punishment predictive cues. In sum, we have identified NT as a neuropeptide that signals valence in the BLA, and showed that NT is a critical neuromodulator that orchestrates positive and negative valence assignment in amygdala neurons by extending valence-specific plasticity to behaviourally relevant timescales.

MST2 methylation by PRMT5 inhibits Hippo signaling and promotes pancreatic cancer progression.

The Hippo signaling axis is a tumor suppressor pathway that is activated by various extra-pathway factors to regulate cell differentiation and organ development. Recent studies have reported that autophosphorylation of the core kinase cassette stimulates activation of the Hippo signaling cascade. Here, we demonstrate that protein arginine methyltransferase 5 (PRMT5) contributes to inactivation of the Hippo signaling pathway in pancreatic cancer. We show that the Hippo pathway initiator serine/threonine kinase 3 (STK3, also known as MST2) of Hippo signaling pathway can be symmetrically di-methylated by PRMT5 at arginine-461 (R461) and arginine-467 (R467) in its SARAH domain. Methylation suppresses MST2 autophosphorylation and kinase activity by blocking its homodimerization, thereby inactivating Hippo signaling pathway in pancreatic cancer. Moreover, we also show that the specific PRMT5 inhibitor GSK3326595 re-activates the dysregulated Hippo signaling pathway and inhibits the growth of human pancreatic cancer xenografts in immunodeficient mice, thus suggesting potential clinical application of PRMT5 inhibitors in pancreatic cancer.

Whole-genome sequencing in autism identifies hot spots for de novo germline mutation.

De novo mutation plays an important role in autism spectrum disorders (ASDs). Notably, pathogenic copy number variants (CNVs) are characterized by high mutation rates. We hypothesize that hypermutability is a property of ASD genes and may also include nucleotide-substitution hot spots. We investigated global patterns of germline mutation by whole-genome sequencing of monozygotic twins concordant for ASD and their parents. Mutation rates varied widely throughout the genome (by 100-fold) and could be explained by intrinsic characteristics of DNA sequence and chromatin structure. Dense clusters of mutations within individual genomes were attributable to compound mutation or gene conversion. Hypermutability was a characteristic of genes involved in ASD and other diseases. In addition, genes impacted by mutations in this study were associated with ASD in independent exome-sequencing data sets. Our findings suggest that regional hypermutation is a significant factor shaping patterns of genetic variation and disease risk in humans.

Corticosterone inhibits GAS6 to govern hair follicle stem-cell quiescence.

Chronic, sustained exposure to stressors can profoundly affect tissue homeostasis, although the mechanisms by which these changes occur are largely unknown. Here we report that the stress hormone corticosterone-which is derived from the adrenal gland and is the rodent equivalent of cortisol in humans-regulates hair follicle stem cell (HFSC) quiescence and hair growth in mice. In the absence of systemic corticosterone, HFSCs enter substantially more rounds of the regeneration cycle throughout life. Conversely, under chronic stress, increased levels of corticosterone prolong HFSC quiescence and maintain hair follicles in an extended resting phase. Mechanistically, corticosterone acts on the dermal papillae to suppress the expression of Gas6, a gene that encodes the secreted factor growth arrest specific 6. Restoring Gas6 expression overcomes the stress-induced inhibition of HFSC activation and hair growth. Our work identifies corticosterone as a systemic inhibitor of HFSC activity through its effect on the niche, and demonstrates that the removal of such inhibition drives HFSCs into frequent regeneration cycles, with no observable defects in the long-term.

Epigenomic diversity of cortical projection neurons in the mouse brain.

Neuronal cell types are classically defined by their molecular properties, anatomy and functions. Although recent advances in single-cell genomics have led to high-resolution molecular characterization of cell type diversity in the brain1, neuronal cell types are often studied out of the context of their anatomical properties. To improve our understanding of the relationship between molecular and anatomical features that define cortical neurons, here we combined retrograde labelling with single-nucleus DNA methylation sequencing to link neural epigenomic properties to projections. We examined 11,827 single neocortical neurons from 63 cortico-cortical and cortico-subcortical long-distance projections. Our results showed unique epigenetic signatures of projection neurons that correspond to their laminar and regional location and projection patterns. On the basis of their epigenomes, intra-telencephalic cells that project to different cortical targets could be further distinguished, and some layer 5 neurons that project to extra-telencephalic targets (L5 ET) formed separate clusters that aligned with their axonal projections. Such separation varied between cortical areas, which suggests that there are area-specific differences in L5 ET subtypes, which were further validated by anatomical studies. Notably, a population of cortico-cortical projection neurons clustered with L5 ET rather than intra-telencephalic neurons, which suggests that a population of L5 ET cortical neurons projects to both targets. We verified the existence of these neurons by dual retrograde labelling and anterograde tracing of cortico-cortical projection neurons, which revealed axon terminals in extra-telencephalic targets including the thalamus, superior colliculus and pons. These findings highlight the power of single-cell epigenomic approaches to connect the molecular properties of neurons with their anatomical and projection properties.