Positive selection of somatically mutated clones identifies adaptive pathways in metabolic liver disease.

Somatic mutations in nonmalignant tissues accumulate with age and injury, but whether these mutations are adaptive on the cellular or organismal levels is unclear. To interrogate genes in human metabolic disease, we performed lineage tracing in mice harboring somatic mosaicism subjected to nonalcoholic steatohepatitis (NASH). Proof-of-concept studies with mosaic loss of Mboat7, a membrane lipid acyltransferase, showed that increased steatosis accelerated clonal disappearance. Next, we induced pooled mosaicism in 63 known NASH genes, allowing us to trace mutant clones side by side. This in vivo tracing platform, which we coined MOSAICS, selected for mutations that ameliorate lipotoxicity, including mutant genes identified in human NASH. To prioritize new genes, additional screening of 472 candidates identified 23 somatic perturbations that promoted clonal expansion. In validation studies, liver-wide deletion of Tbx3, Bcl6, or Smyd2 resulted in protection against hepatic steatosis. Selection for clonal fitness in mouse and human livers identifies pathways that regulate metabolic disease.

FOXP1 and KLF2 reciprocally regulate checkpoints of stem-like to effector transition in CAR T cells.

In cancer and infections, self-renewing stem-like CD8+ T cells mediate the response of immunotherapies and replenish terminally exhausted T cells and effector-like T cells. However, the programs governing the lineage choice in chimeric antigen receptor (CAR) T cells are unclear. Here, by simultaneously profiling single-cell chromatin accessibility and transcriptome in the same CAR T cells, we identified heterogeneous chromatin states within CD8+ T cell subsets that foreshadowed transcriptional changes and were primed for regulation by distinct transcription factors. Transcription factors that controlled each CD8+ T cell subset were regulated by high numbers of enhancers and positioned as hubs of gene networks. FOXP1, a hub in the stem-like network, promoted expansion and stemness of CAR T cells and limited excessive effector differentiation. In the effector network, KLF2 enhanced effector CD8+ T cell differentiation and prevented terminal exhaustion. Thus, we identified gene networks and hub transcription factors that controlled the differentiation of stem-like CD8+ CAR T cells into effector or exhausted CD8+ CAR T cells.

Coupling of N7-methyltransferase and 3'-5' exoribonuclease with SARS-CoV-2 polymerase reveals mechanisms for capping and proofreading.

The capping of mRNA and the proofreading play essential roles in SARS-CoV-2 replication and transcription. Here, we present the cryo-EM structure of the SARS-CoV-2 replication-transcription complex (RTC) in a form identified as Cap(0)-RTC, which couples a co-transcriptional capping complex (CCC) composed of nsp12 NiRAN, nsp9, the bifunctional nsp14 possessing an N-terminal exoribonuclease (ExoN) and a C-terminal N7-methyltransferase (N7-MTase), and nsp10 as a cofactor of nsp14. Nsp9 and nsp12 NiRAN recruit nsp10/nsp14 into the Cap(0)-RTC, forming the N7-CCC to yield cap(0) (7MeGpppA) at 5' end of pre-mRNA. A dimeric form of Cap(0)-RTC observed by cryo-EM suggests an in trans backtracking mechanism for nsp14 ExoN to facilitate proofreading of the RNA in concert with polymerase nsp12. These results not only provide a structural basis for understanding co-transcriptional modification of SARS-CoV-2 mRNA but also shed light on how replication fidelity in SARS-CoV-2 is maintained.

Cryo-EM Structure of an Extended SARS-CoV-2 Replication and Transcription Complex Reveals an Intermediate State in Cap Synthesis.

Transcription of SARS-CoV-2 mRNA requires sequential reactions facilitated by the replication and transcription complex (RTC). Here, we present a structural snapshot of SARS-CoV-2 RTC as it transitions toward cap structure synthesis. We determine the atomic cryo-EM structure of an extended RTC assembled by nsp7-nsp82-nsp12-nsp132-RNA and a single RNA-binding protein, nsp9. Nsp9 binds tightly to nsp12 (RdRp) NiRAN, allowing nsp9 N terminus inserting into the catalytic center of nsp12 NiRAN, which then inhibits activity. We also show that nsp12 NiRAN possesses guanylyltransferase activity, catalyzing the formation of cap core structure (GpppA). The orientation of nsp13 that anchors the 5' extension of template RNA shows a remarkable conformational shift, resulting in zinc finger 3 of its ZBD inserting into a minor groove of paired template-primer RNA. These results reason an intermediate state of RTC toward mRNA synthesis, pave a way to understand the RTC architecture, and provide a target for antiviral development.

Influenza-trained mucosal-resident alveolar macrophages confer long-term antitumor immunity in the lungs.

Respiratory viral infections reprogram pulmonary macrophages with altered anti-infectious functions. However, the potential function of virus-trained macrophages in antitumor immunity in the lung, a preferential target of both primary and metastatic malignancies, is not well understood. Using mouse models of influenza and lung metastatic tumors, we show here that influenza trains respiratory mucosal-resident alveolar macrophages (AMs) to exert long-lasting and tissue-specific antitumor immunity. Trained AMs infiltrate tumor lesions and have enhanced phagocytic and tumor cell cytotoxic functions, which are associated with epigenetic, transcriptional and metabolic resistance to tumor-induced immune suppression. Generation of antitumor trained immunity in AMs is dependent on interferon-γ and natural killer cells. Notably, human AMs with trained immunity traits in non-small cell lung cancer tissue are associated with a favorable immune microenvironment. These data reveal a function for trained resident macrophages in pulmonary mucosal antitumor immune surveillance. Induction of trained immunity in tissue-resident macrophages might thereby be a potential antitumor strategy.

Single-Cell Analyses Inform Mechanisms of Myeloid-Targeted Therapies in Colon Cancer.

Single-cell RNA sequencing (scRNA-seq) is a powerful tool for defining cellular diversity in tumors, but its application toward dissecting mechanisms underlying immune-modulating therapies is scarce. We performed scRNA-seq analyses on immune and stromal populations from colorectal cancer patients, identifying specific macrophage and conventional dendritic cell (cDC) subsets as key mediators of cellular cross-talk in the tumor microenvironment. Defining comparable myeloid populations in mouse tumors enabled characterization of their response to myeloid-targeted immunotherapy. Treatment with anti-CSF1R preferentially depleted macrophages with an inflammatory signature but spared macrophage populations that in mouse and human expresses pro-angiogenic/tumorigenic genes. Treatment with a CD40 agonist antibody preferentially activated a cDC population and increased Bhlhe40+ Th1-like cells and CD8+ memory T cells. Our comprehensive analysis of key myeloid subsets in human and mouse identifies critical cellular interactions regulating tumor immunity and defines mechanisms underlying myeloid-targeted immunotherapies currently undergoing clinical testing.

Somatic Mutations Increase Hepatic Clonal Fitness and Regeneration in Chronic Liver Disease.

Normal tissues accumulate genetic changes with age, but it is unknown if somatic mutations promote clonal expansion of non-malignant cells in the setting of chronic degenerative diseases. Exome sequencing of diseased liver samples from 82 patients revealed a complex mutational landscape in cirrhosis. Additional ultra-deep sequencing identified recurrent mutations in PKD1, PPARGC1B, KMT2D, and ARID1A. The number and size of mutant clones increased as a function of fibrosis stage and tissue damage. To interrogate the functional impact of mutated genes, a pooled in vivo CRISPR screening approach was established. In agreement with sequencing results, examination of 147 genes again revealed that loss of Pkd1, Kmt2d, and Arid1a promoted clonal expansion. Conditional heterozygous deletion of these genes in mice was also hepatoprotective in injury assays. Pre-malignant somatic alterations are often viewed through the lens of cancer, but we show that mutations can promote regeneration, likely independent of carcinogenesis.

Colony stimulating factor-1 producing endothelial cells and mesenchymal stromal cells maintain monocytes within a perivascular bone marrow niche.

Macrophage colony stimulating factor-1 (CSF-1) plays a critical role in maintaining myeloid lineage cells. However, congenital global deficiency of CSF-1 (Csf1op/op) causes severe musculoskeletal defects that may indirectly affect hematopoiesis. Indeed, we show here that osteolineage-derived Csf1 prevented developmental abnormalities but had no effect on monopoiesis in adulthood. However, ubiquitous deletion of Csf1 conditionally in adulthood decreased monocyte survival, differentiation, and migration, independent of its effects on bone development. Bone histology revealed that monocytes reside near sinusoidal endothelial cells (ECs) and leptin receptor (Lepr)-expressing perivascular mesenchymal stromal cells (MSCs). Targeted deletion of Csf1 from sinusoidal ECs selectively reduced Ly6C- monocytes, whereas combined depletion of Csf1 from ECs and MSCs further decreased Ly6Chi cells. Moreover, EC-derived CSF-1 facilitated recovery of Ly6C- monocytes and protected mice from weight loss following induction of polymicrobial sepsis. Thus, monocytes are supported by distinct cellular sources of CSF-1 within a perivascular BM niche.

SARS-CoV-2 exacerbates proinflammatory responses in myeloid cells through C-type lectin receptors and Tweety family member 2.

Despite mounting evidence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) engagement with immune cells, most express little, if any, of the canonical receptor of SARS-CoV-2, angiotensin-converting enzyme 2 (ACE2). Here, using a myeloid cell receptor-focused ectopic expression screen, we identified several C-type lectins (DC-SIGN, L-SIGN, LSECtin, ASGR1, and CLEC10A) and Tweety family member 2 (TTYH2) as glycan-dependent binding partners of the SARS-CoV-2 spike. Except for TTYH2, these molecules primarily interacted with spike via regions outside of the receptor-binding domain. Single-cell RNA sequencing analysis of pulmonary cells from individuals with coronavirus disease 2019 (COVID-19) indicated predominant expression of these molecules on myeloid cells. Although these receptors do not support active replication of SARS-CoV-2, their engagement with the virus induced robust proinflammatory responses in myeloid cells that correlated with COVID-19 severity. We also generated a bispecific anti-spike nanobody that not only blocked ACE2-mediated infection but also the myeloid receptor-mediated proinflammatory responses. Our findings suggest that SARS-CoV-2-myeloid receptor interactions promote immune hyperactivation, which represents potential targets for COVID-19 therapy.

Fossils document evolutionary changes of jaw joint to mammalian middle ear.

The dual jaw joint of Morganucodon1,2 consists of the dentary-squamosal joint laterally and the articular-quadrate one medially. The articular-quadrate joint and its associated post-dentary bones constitute the precursor of the mammalian middle ear. Fossils documenting the transition from such a precursor to the mammalian middle ear are poor, resulting in inconsistent interpretations of this hallmark apparatus in the earliest stage of mammaliaform evolution1-5. Here we report mandibular middle ears from two Jurassic mammaliaforms: a new morganucodontan-like species and a pseudotribosphenic shuotheriid species6. The morganucodontan-like species shows many previously unknown post-dentary bone morphologies1,2 and exhibits features that suggest a loss of load-bearing function in its articular-quadrate joint. The middle ear of the shuotheriid approaches the mammalian condition in that it has features that are suitable for an exclusively auditory function, although the post-dentary bones are still attached to the dentary. With size reduction of the jaw-joint bones, the quadrate shifts medially at different degrees in relation to the articular in the two mammaliaforms. These changes provide evidence of a gradual loss of load-bearing function in the articular-quadrate jaw joint-a prerequisite for the detachment of the post-dentary bones from the dentary7-12 and the eventual breakdown of the Meckel's cartilage13-15 during the evolution of mammaliaforms.

Black holes regulate cool gas accretion in massive galaxies.

The nucleus of almost all massive galaxies contains a supermassive black hole (BH)1. The feedback from the accretion of these BHs is often considered to have crucial roles in establishing the quiescence of massive galaxies2-14, although some recent studies show that even galaxies hosting the most active BHs do not exhibit a reduction in their molecular gas reservoirs or star formation rates15-17. Therefore, the influence of BHs on galaxy star formation remains highly debated and lacks direct evidence. Here, based on a large sample of nearby galaxies with measurements of masses of both BHs and atomic hydrogen (HI), the main component of the interstellar medium18, we show that the HI gas mass to stellar masses ratio (µHI = MHI/M⋆) is more strongly correlated with BH masses (MBH) than with any other galaxy parameters, including stellar mass, stellar mass surface density and bulge masses. Moreover, once the µHI-MBH correlation is considered, µHI loses dependence on other galactic parameters, demonstrating that MBH serves as the primary driver of µHI. These findings provide important evidence for how the accumulated energy from BH accretion regulates the cool gas content in galaxies, by ejecting interstellar medium gas and/or suppressing gas cooling from the circumgalactic medium.

Hybrid 2D-CMOS microchips for memristive applications.

Exploiting the excellent electronic properties of two-dimensional (2D) materials to fabricate advanced electronic circuits is a major goal for the semiconductor industry1,2. However, most studies in this field have been limited to the fabrication and characterization of isolated large (more than 1 µm2) devices on unfunctional SiO2-Si substrates. Some studies have integrated monolayer graphene on silicon microchips as a large-area (more than 500 µm2) interconnection3 and as a channel of large transistors (roughly 16.5 µm2) (refs. 4,5), but in all cases the integration density was low, no computation was demonstrated and manipulating monolayer 2D materials was challenging because native pinholes and cracks during transfer increase variability and reduce yield. Here, we present the fabrication of high-integration-density 2D-CMOS hybrid microchips for memristive applications-CMOS stands for complementary metal-oxide-semiconductor. We transfer a sheet of multilayer hexagonal boron nitride onto the back-end-of-line interconnections of silicon microchips containing CMOS transistors of the 180 nm node, and finalize the circuits by patterning the top electrodes and interconnections. The CMOS transistors provide outstanding control over the currents across the hexagonal boron nitride memristors, which allows us to achieve endurances of roughly 5 million cycles in memristors as small as 0.053 µm2. We demonstrate in-memory computation by constructing logic gates, and measure spike-timing dependent plasticity signals that are suitable for the implementation of spiking neural networks. The high performance and the relatively-high technology readiness level achieved represent a notable advance towards the integration of 2D materials in microelectronic products and memristive applications.

Histamine synthesis and transport are coupled in axon terminals via a dual quality control system.

Monoamine neurotransmitters generated by de novo synthesis are rapidly transported and stored into synaptic vesicles at axon terminals. This transport is essential both for sustaining synaptic transmission and for limiting the toxic effects of monoamines. Here, synthesis of the monoamine histamine by histidine decarboxylase (HDC) and subsequent loading of histamine into synaptic vesicles are shown to be physically and functionally coupled within Drosophila photoreceptor terminals. This process requires HDC anchoring to synaptic vesicles via interactions with N-ethylmaleimide-sensitive fusion protein 1 (NSF1). Disassociating HDC from synaptic vesicles disrupts visual synaptic transmission and causes somatic accumulation of histamine, which leads to retinal degeneration. We further identified a proteasome degradation system mediated by the E3 ubiquitin ligase, purity of essence (POE), which clears mislocalized HDC from the soma, thus eliminating the cytotoxic effects of histamine. Taken together, our results reveal a dual mechanism for translocation and degradation of HDC that ensures restriction of histamine synthesis to axonal terminals and at the same time rapid loading into synaptic vesicles. This is crucial for sustaining neurotransmission and protecting against cytotoxic monoamines.

Epigenomic analysis of multilineage differentiation of human embryonic stem cells.

Epigenetic mechanisms have been proposed to play crucial roles in mammalian development, but their precise functions are only partially understood. To investigate epigenetic regulation of embryonic development, we differentiated human embryonic stem cells into mesendoderm, neural progenitor cells, trophoblast-like cells, and mesenchymal stem cells and systematically characterized DNA methylation, chromatin modifications, and the transcriptome in each lineage. We found that promoters that are active in early developmental stages tend to be CG rich and mainly engage H3K27me3 upon silencing in nonexpressing lineages. By contrast, promoters for genes expressed preferentially at later stages are often CG poor and primarily employ DNA methylation upon repression. Interestingly, the early developmental regulatory genes are often located in large genomic domains that are generally devoid of DNA methylation in most lineages, which we termed DNA methylation valleys (DMVs). Our results suggest that distinct epigenetic mechanisms regulate early and late stages of ES cell differentiation.

Haplotype-resolved 3D chromatin architecture of the hybrid pig.

In diploid mammals, allele-specific three-dimensional (3D) genome architecture may lead to imbalanced gene expression. Through ultradeep in situ Hi-C sequencing of three representative somatic tissues (liver, skeletal muscle, and brain) from hybrid pigs generated by reciprocal crosses of phenotypically and physiologically divergent Berkshire and Tibetan pigs, we uncover extensive chromatin reorganization between homologous chromosomes across multiple scales. Haplotype-based interrogation of multi-omic data revealed the tissue dependence of 3D chromatin conformation, suggesting that parent-of-origin-specific conformation may drive gene imprinting. We quantify the effects of genetic variations and histone modifications on allelic differences of long-range promoter-enhancer contacts, which likely contribute to the phenotypic differences between the parental pig breeds. We also observe the fine structure of somatically paired homologous chromosomes in the pig genome, which has a functional implication genome-wide. This work illustrates how allele-specific chromatin architecture facilitates concomitant shifts in allele-biased gene expression, as well as the possible consequential phenotypic changes in mammals.

Mapping cellular interactions from spatially resolved transcriptomics data.

Cell-cell communication (CCC) is essential to how life forms and functions. However, accurate, high-throughput mapping of how expression of all genes in one cell affects expression of all genes in another cell is made possible only recently through the introduction of spatially resolved transcriptomics (SRT) technologies, especially those that achieve single-cell resolution. Nevertheless, substantial challenges remain to analyze such highly complex data properly. Here, we introduce a multiple-instance learning framework, Spacia, to detect CCCs from data generated by SRTs, by uniquely exploiting their spatial modality. We highlight Spacia's power to overcome fundamental limitations of popular analytical tools for inference of CCCs, including losing single-cell resolution, limited to ligand-receptor relationships and prior interaction databases, high false positive rates and, most importantly, the lack of consideration of the multiple-sender-to-one-receiver paradigm. We evaluated the fitness of Spacia for three commercialized single-cell resolution SRT technologies: MERSCOPE/Vizgen, CosMx/NanoString and Xenium/10x. Overall, Spacia represents a notable step in advancing quantitative theories of cellular communications.

A cargo sorting receptor mediates chloroplast protein trafficking through the secretory pathway.

Nucleus-encoded chloroplast proteins can be transported via the secretory pathway. The molecular mechanisms underlying the trafficking of chloroplast proteins between the intracellular compartments are largely unclear, and a cargo sorting receptor has not previously been identified in the secretory pathway. Here, we report a cargo sorting receptor that is specifically present in Viridiplantae and mediates the transport of cargo proteins to the chloroplast. Using a forward genetic analysis, we identified a gene encoding a transmembrane protein (MtTP930) in barrel medic (Medicago truncatula). Mutation of MtTP930 resulted in impaired chloroplast function and a dwarf phenotype. MtTP930 is highly expressed in the aerial parts of the plant and is localized to the endoplasmic reticulum (ER) exit sites and Golgi. MtTP930 contains typical cargo sorting receptor motifs, interacts with Sar1, Sec12, and Sec24, and participates in coat protein complex II vesicular transport. Importantly, MtTP930 can recognize the cargo proteins plastidial N-glycosylated nucleotide pyrophosphatase/phosphodiesterase (MtNPP) and α-carbonic anhydrase (MtCAH) in the ER and then transport them to the chloroplast via the secretory pathway. Mutation of a homolog of MtTP930 in Arabidopsis (Arabidopsis thaliana) resulted in a similar dwarf phenotype. Furthermore, MtNPP-GFP failed to localize to chloroplasts when transgenically expressed in Attp930 protoplasts, implying that these cargo sorting receptors are conserved in plants. These findings fill a gap in our understanding of the mechanism by which chloroplast proteins are sorted and transported via the secretory pathway.

The S-acylation cycle of transcription factor MtNAC80 influences cold stress responses in Medicago truncatula.

S-acylation is a reversible post-translational modification catalyzed by protein S-acyltransferases (PATs), and acyl protein thioesterases (APTs) mediate de-S-acylation. Although many proteins are S-acylated, how the S-acylation cycle modulates specific biological functions in plants is poorly understood. In this study, we report that the S-acylation cycle of transcription factor MtNAC80 is involved in the Medicago truncatula cold stress response. Under normal conditions, MtNAC80 localized to membranes through MtPAT9-induced S-acylation. In contrast, under cold stress conditions, MtNAC80 translocated to the nucleus through de-S-acylation mediated by thioesterases such as MtAPT1. MtNAC80 functions in the nucleus by directly binding the promoter of the glutathione S-transferase gene MtGSTU1 and promoting its expression, which enables plants to survive under cold stress by removing excess malondialdehyde and H2O2. Our findings reveal an important function of the S-acylation cycle in plants and provide insight into stress response and tolerance mechanisms.

Morphological diversity of single neurons in molecularly defined cell types.

Dendritic and axonal morphology reflects the input and output of neurons and is a defining feature of neuronal types1,2, yet our knowledge of its diversity remains limited. Here, to systematically examine complete single-neuron morphologies on a brain-wide scale, we established a pipeline encompassing sparse labelling, whole-brain imaging, reconstruction, registration and analysis. We fully reconstructed 1,741 neurons from cortex, claustrum, thalamus, striatum and other brain regions in mice. We identified 11 major projection neuron types with distinct morphological features and corresponding transcriptomic identities. Extensive projectional diversity was found within each of these major types, on the basis of which some types were clustered into more refined subtypes. This diversity follows a set of generalizable principles that govern long-range axonal projections at different levels, including molecular correspondence, divergent or convergent projection, axon termination pattern, regional specificity, topography, and individual cell variability. Although clear concordance with transcriptomic profiles is evident at the level of major projection type, fine-grained morphological diversity often does not readily correlate with transcriptomic subtypes derived from unsupervised clustering, highlighting the need for single-cell cross-modality studies. Overall, our study demonstrates the crucial need for quantitative description of complete single-cell anatomy in cell-type classification, as single-cell morphological diversity reveals a plethora of ways in which different cell types and their individual members may contribute to the configuration and function of their respective circuits.