Spatially exploring RNA biology in archival formalin-fixed paraffin-embedded tissues.

The capability to spatially explore RNA biology in formalin-fixed paraffin-embedded (FFPE) tissues holds transformative potential for histopathology research. Here, we present pathology-compatible deterministic barcoding in tissue (Patho-DBiT) by combining in situ polyadenylation and computational innovation for spatial whole transcriptome sequencing, tailored to probe the diverse RNA species in clinically archived FFPE samples. It permits spatial co-profiling of gene expression and RNA processing, unveiling region-specific splicing isoforms, and high-sensitivity transcriptomic mapping of clinical tumor FFPE tissues stored for 5 years. Furthermore, genome-wide single-nucleotide RNA variants can be captured to distinguish malignant subclones from non-malignant cells in human lymphomas. Patho-DBiT also maps microRNA regulatory networks and RNA splicing dynamics, decoding their roles in spatial tumorigenesis. Single-cell level Patho-DBiT dissects the spatiotemporal cellular dynamics driving tumor clonal architecture and progression. Patho-DBiT stands poised as a valuable platform to unravel rich RNA biology in FFPE tissues to aid in clinical pathology evaluation.

High-Spatial-Resolution Multi-Omics Sequencing via Deterministic Barcoding in Tissue.

We present deterministic barcoding in tissue for spatial omics sequencing (DBiT-seq) for co-mapping of mRNAs and proteins in a formaldehyde-fixed tissue slide via next-generation sequencing (NGS). Parallel microfluidic channels were used to deliver DNA barcodes to the surface of a tissue slide, and crossflow of two sets of barcodes, A1-50 and B1-50, followed by ligation in situ, yielded a 2D mosaic of tissue pixels, each containing a unique full barcode AB. Application to mouse embryos revealed major tissue types in early organogenesis as well as fine features like microvasculature in a brain and pigmented epithelium in an eye field. Gene expression profiles in 10-µm pixels conformed into the clusters of single-cell transcriptomes, allowing for rapid identification of cell types and spatial distributions. DBiT-seq can be adopted by researchers with no experience in microfluidics and may find applications in a range of fields including developmental biology, cancer biology, neuroscience, and clinical pathology.

Generation and Application of Mouse-Rat Allodiploid Embryonic Stem Cells.

Mammalian interspecific hybrids provide unique advantages for mechanistic studies of speciation, gene expression regulation, and X chromosome inactivation (XCI) but are constrained by their limited natural resources. Previous artificially generated mammalian interspecific hybrid cells are usually tetraploids with unstable genomes and limited developmental abilities. Here, we report the generation of mouse-rat allodiploid embryonic stem cells (AdESCs) by fusing haploid ESCs of the two species. The AdESCs have a stable allodiploid genome and are capable of differentiating into all three germ layers and early-stage germ cells. Both the mouse and rat alleles have comparable contributions to the expression of most genes. We have proven AdESCs as a powerful tool to study the mechanisms regulating X chromosome inactivation and to identify X inactivation-escaping genes, as well as to efficiently identify genes regulating phenotypic differences between species. A similar method could be used to create hybrid AdESCs of other distantly related species.

An intrinsically disordered region controlling condensation of a circadian clock component and rhythmic transcription in the liver.

Circadian gene transcription is fundamental to metabolic physiology. Here we report that the nuclear receptor REV-ERBα, a repressive component of the molecular clock, forms circadian condensates in the nuclei of mouse liver. These condensates are dictated by an intrinsically disordered region (IDR) located in the protein's hinge region which specifically concentrates nuclear receptor corepressor 1 (NCOR1) at the genome. IDR deletion diminishes the recruitment of NCOR1 and disrupts rhythmic gene transcription in vivo. REV-ERBα condensates are located at high-order transcriptional repressive hubs in the liver genome that are highly correlated with circadian gene repression. Deletion of the IDR disrupts transcriptional repressive hubs and diminishes silencing of target genes by REV-ERBα. This work demonstrates physiological circadian protein condensates containing REV-ERBα whose IDR is required for hub formation and the control of rhythmic gene expression.

The chemokine CCL1 triggers an AMFR-SPRY1 pathway that promotes differentiation of lung fibroblasts into myofibroblasts and drives pulmonary fibrosis.

Recruitment of immune cells to the site of inflammation by the chemokine CCL1 is important in the pathology of inflammatory diseases. Here, we examined the role of CCL1 in pulmonary fibrosis (PF). Bronchoalveolar lavage fluid from PF mouse models contained high amounts of CCL1, as did lung biopsies from PF patients. Immunofluorescence analyses revealed that alveolar macrophages and CD4+ T cells were major producers of CCL1 and targeted deletion of Ccl1 in these cells blunted pathology. Deletion of the CCL1 receptor Ccr8 in fibroblasts limited migration, but not activation, in response to CCL1. Mass spectrometry analyses of CCL1 complexes identified AMFR as a CCL1 receptor, and deletion of Amfr impaired fibroblast activation. Mechanistically, CCL1 binding triggered ubiquitination of the ERK inhibitor Spry1 by AMFR, thus activating Ras-mediated profibrotic protein synthesis. Antibody blockade of CCL1 ameliorated PF pathology, supporting the therapeutic potential of targeting this pathway for treating fibroproliferative lung diseases.

Isoform-specific functions of PPARγ in gene regulation and metabolism.

Peroxisome proliferator-activated receptor γ (PPARγ) is a nuclear receptor that is a vital regulator of adipogenesis, insulin sensitivity, and lipid metabolism. Activation of PPARγ by antidiabetic thiazolidinediones (TZD) reverses insulin resistance but also leads to weight gain that limits the use of these drugs. There are two main PPARγ isoforms, but the specific functions of each are not established. Here we generated mouse lines in which endogenous PPARγ1 and PPARγ2 were epitope-tagged to interrogate isoform-specific genomic binding, and mice deficient in either PPARγ1 or PPARγ2 to assess isoform-specific gene regulation. Strikingly, although PPARγ1 and PPARγ2 contain identical DNA binding domains, we uncovered isoform-specific genomic binding sites in addition to shared sites. Moreover, PPARγ1 and PPARγ2 regulated a different set of genes in adipose tissue depots, suggesting distinct roles in adipocyte biology. Indeed, mice with selective deficiency of PPARγ1 maintained body temperature better than wild-type or PPARγ2-deficient mice. Most remarkably, although TZD treatment improved glucose tolerance in mice lacking either PPARγ1 or PPARγ2, the PPARγ1-deficient mice were protected from TZD-induced body weight gain compared with PPARγ2-deficient mice. Thus, PPARγ isoforms have specific and separable metabolic functions that may be targeted to improve therapy for insulin resistance and diabetes.

Spatial epigenome-transcriptome co-profiling of mammalian tissues.

Emerging spatial technologies, including spatial transcriptomics and spatial epigenomics, are becoming powerful tools for profiling of cellular states in the tissue context1-5. However, current methods capture only one layer of omics information at a time, precluding the possibility of examining the mechanistic relationship across the central dogma of molecular biology. Here, we present two technologies for spatially resolved, genome-wide, joint profiling of the epigenome and transcriptome by cosequencing chromatin accessibility and gene expression, or histone modifications (H3K27me3, H3K27ac or H3K4me3) and gene expression on the same tissue section at near-single-cell resolution. These were applied to embryonic and juvenile mouse brain, as well as adult human brain, to map how epigenetic mechanisms control transcriptional phenotype and cell dynamics in tissue. Although highly concordant tissue features were identified by either spatial epigenome or spatial transcriptome we also observed distinct patterns, suggesting their differential roles in defining cell states. Linking epigenome to transcriptome pixel by pixel allows the uncovering of new insights in spatial epigenetic priming, differentiation and gene regulation within the tissue architecture. These technologies are of great interest in life science and biomedical research.

Enantioconvergent Cu-catalysed N-alkylation of aliphatic amines.

Chiral amines are commonly used in the pharmaceutical and agrochemical industries1. The strong demand for unnatural chiral amines has driven the development of catalytic asymmetric methods1,2. Although the N-alkylation of aliphatic amines with alkyl halides has been widely adopted for over 100 years, catalyst poisoning and unfettered reactivity have been preventing the development of a catalyst-controlled enantioselective version3-5. Here we report the use of chiral tridentate anionic ligands to enable the copper-catalysed chemoselective and enantioconvergent N-alkylation of aliphatic amines with α-carbonyl alkyl chlorides. This method can directly convert feedstock chemicals, including ammonia and pharmaceutically relevant amines, into unnatural chiral α-amino amides under mild and robust conditions. Excellent enantioselectivity and functional-group tolerance were observed. The power of the method is demonstrated in a number of complex settings, including late-stage functionalization and in the expedited synthesis of diverse amine drug molecules. The current method indicates that multidentate anionic ligands are a general solution for overcoming transition-metal-catalyst poisoning.

An electric molecular motor.

Macroscopic electric motors continue to have a large impact on almost every aspect of modern society. Consequently, the effort towards developing molecular motors1-3 that can be driven by electricity could not be more timely. Here we describe an electric molecular motor based on a [3]catenane4,5, in which two cyclobis(paraquat-p-phenylene)6 (CBPQT4+) rings are powered by electricity in solution to circumrotate unidirectionally around a 50-membered loop. The constitution of the loop ensures that both rings undergo highly (85%) unidirectional movement under the guidance of a flashing energy ratchet7,8, whereas the interactions between the two rings give rise to a two-dimensional potential energy surface (PES) similar to that shown by FOF1 ATP synthase9. The unidirectionality is powered by an oscillating10 voltage11,12 or external modulation of the redox potential13. Initially, we focused our attention on the homologous [2]catenane, only to find that the kinetic asymmetry was insufficient to support unidirectional movement of the sole ring. Accordingly, we incorporated a second CBPQT4+ ring to provide further symmetry breaking by interactions between the two mobile rings. This demonstration of electrically driven continual circumrotatory motion of two rings around a loop in a [3]catenane is free from the production of waste products and represents an important step towards surface-bound14 electric molecular motors.

Spatial profiling of chromatin accessibility in mouse and human tissues.

Cellular function in tissue is dependent on the local environment, requiring new methods for spatial mapping of biomolecules and cells in the tissue context1. The emergence of spatial transcriptomics has enabled genome-scale gene expression mapping2-5, but the ability to capture spatial epigenetic information of tissue at the cellular level and genome scale is lacking. Here we describe a method for spatially resolved chromatin accessibility profiling of tissue sections using next-generation sequencing (spatial-ATAC-seq) by combining in situ Tn5 transposition chemistry6 and microfluidic deterministic barcoding5. Profiling mouse embryos using spatial-ATAC-seq delineated tissue-region-specific epigenetic landscapes and identified gene regulators involved in the development of the central nervous system. Mapping the accessible genome in the mouse and human brain revealed the intricate arealization of brain regions. Applying spatial-ATAC-seq to tonsil tissue resolved the spatially distinct organization of immune cell types and states in lymphoid follicles and extrafollicular zones. This technology progresses spatial biology by enabling spatially resolved chromatin accessibility profiling to improve our understanding of cell identity, cell state and cell fate decision in relation to epigenetic underpinnings in development and disease.

Electron-catalysed molecular recognition.

Molecular recognition1-4 and supramolecular assembly5-8 cover a broad spectrum9-11 of non-covalently orchestrated phenomena between molecules. Catalysis12 of such processes, however, unlike that for the formation of covalent bonds, is limited to approaches13-16 that rely on sophisticated catalyst design. Here we establish a simple and versatile strategy to facilitate molecular recognition by extending electron catalysis17, which is widely applied18-21 in synthetic covalent chemistry, into the realm of supramolecular non-covalent chemistry. As a proof of principle, we show that the formation of a trisradical complex22 between a macrocyclic host and a dumbbell-shaped guest-a molecular recognition process that is kinetically forbidden under ambient conditions-can be accelerated substantially on the addition of catalytic amounts of a chemical electron source. It is, therefore, electrochemically possible to control23 the molecular recognition temporally and produce a nearly arbitrary molar ratio between the substrates and complexes ranging between zero and the equilibrium value. Such kinetically stable supramolecular systems24 are difficult to obtain precisely by other means. The use of the electron as a catalyst in molecular recognition will inspire chemists and biologists to explore strategies that can be used to fine-tune non-covalent events, control assembly at different length scales25-27 and ultimately create new forms of complex matter28-30.

Hippo pathway deletion in adult resting cardiac fibroblasts initiates a cell state transition with spontaneous and self-sustaining fibrosis.

Cardiac fibroblasts (CFs) respond to injury by transitioning through multiple cell states, including resting CFs, activated CFs, and myofibroblasts. We report here that Hippo signaling cell-autonomously regulates CF fate transitions and proliferation, and non-cell-autonomously regulates both myeloid and CF activation in the heart. Conditional deletion of Hippo pathway kinases, Lats1 and Lats2, in uninjured CFs initiated a self-perpetuating fibrotic response in the adult heart that was exacerbated by myocardial infarction (MI). Single cell transcriptomics showed that uninjured Lats1/2 mutant CFs spontaneously transitioned to a myofibroblast cell state. Through gene regulatory network reconstruction, we found that Hippo-deficient myofibroblasts deployed a network of transcriptional regulators of endoplasmic reticulum (ER) stress, and the unfolded protein response (UPR) consistent with elevated secretory activity. We observed an expansion of myeloid cell heterogeneity in uninjured Lats1/2 CKO hearts with similarity to cells recovered from control hearts post-MI. Integrated genome-wide analysis of Yap chromatin occupancy revealed that Yap directly activates myofibroblast cell identity genes, the proto-oncogene Myc, and an array of genes encoding pro-inflammatory factors through enhancer-promoter looping. Our data indicate that Lats1/2 maintain the resting CF cell state through restricting the Yap-induced injury response.

A global phylogenomic analysis of the shiitake genus Lentinula.

Lentinula is a broadly distributed group of fungi that contains the cultivated shiitake mushroom, L. edodes. We sequenced 24 genomes representing eight described species and several unnamed lineages of Lentinula from 15 countries on four continents. Lentinula comprises four major clades that arose in the Oligocene, three in the Americas and one in Asia-Australasia. To expand sampling of shiitake mushrooms, we assembled 60 genomes of L. edodes from China that were previously published as raw Illumina reads and added them to our dataset. Lentinula edodes sensu lato (s. lat.) contains three lineages that may warrant recognition as species, one including a single isolate from Nepal that is the sister group to the rest of L. edodes s. lat., a second with 20 cultivars and 12 wild isolates from China, Japan, Korea, and the Russian Far East, and a third with 28 wild isolates from China, Thailand, and Vietnam. Two additional lineages in China have arisen by hybridization among the second and third groups. Genes encoding cysteine sulfoxide lyase (lecsl) and γ-glutamyl transpeptidase (leggt), which are implicated in biosynthesis of the organosulfur flavor compound lenthionine, have diversified in Lentinula. Paralogs of both genes that are unique to Lentinula (lecsl 3 and leggt 5b) are coordinately up-regulated in fruiting bodies of L. edodes. The pangenome of L. edodes s. lat. contains 20,308 groups of orthologous genes, but only 6,438 orthogroups (32%) are shared among all strains, whereas 3,444 orthogroups (17%) are found only in wild populations, which should be targeted for conservation.

Deep ancestral introgressions between ovine species shape sheep genomes via argali-mediated gene flow.

Previous studies reveal extensive genetic introgression between Ovis species, which affects genetic adaptation and morphological traits. However, the exact evolutionary scenarios underlying the hybridization between sheep and allopatric wild relatives remain unknown. To address this problem, we here integrate the reference genomes of several ovine and caprine species: domestic sheep, argali, bighorn sheep, snow sheep, and domestic goats. Additionally, we use 856 whole genomes representing 169 domestic sheep populations and their 6 wild relatives: Asiatic mouflon, urial, argali, snow sheep, thinhorn sheep and bighorn sheep. We implement a comprehensive set of analyses to test introgression among these species. We infer that the argali lineage originated ca. 3.08-3.35 Mya and hybridized with the ancestor of Pachyceriforms (e.g., bighorn sheep and snow sheep) at ∼1.56 Mya. Previous studies show apparent introgression from North American Pachyceriforms into the Bashibai sheep, a Chinese native sheep breed, despite their wide geographic separation. We show here that, in fact, the apparent introgression from the Pachyceriforms into Bashibai can be explained by the old introgression from Pachyceriforms into argali, and subsequent recent introgression from argali into Bashibai. Our results illustrate the challenges of estimating complex introgression histories and provide an example of how indirect and direct introgression can be distinguished.

A novel recessive mutation in OXR1 is identified in patient with hearing loss recapitulated by the knockdown zebrafish.

Hereditary hearing loss is a highly genetically heterogeneous disorder. More than 150 genes have been identified to link to human non-syndromic hearing impairment. To identify genetic mutations and underlying molecular mechanisms in affected individuals and families with congenital hearing loss, we recruited a cohort of 389 affected individuals in 354 families for whole-exome sequencing analysis. In this study, we report a novel homozygous missense variant (c.233A > G, p.Lys78Arg) in the OXR1 gene, which was identified in a 4-year-old girl with sensorineural hearing loss. OXR1 encodes Oxidation Resistance 1 and is evolutionarily conserved from zebrafish to human. We found that the ortholog oxr1b gene is expressed in the statoacoustic ganglion (SAG, a sensory ganglion of ear) and posterior lateral line ganglion (pLL) in zebrafish. Knockdown of oxr1b in zebrafish resulted in a significant developmental defect of SAG and pLL. This phenotype can be rescued by co-injection of wild-type human OXR1 mRNAs, but not mutant OXR1 (c.233A > G) mRNAs. OXR1-associated pathway analysis revealed that mutations of TBC1D24, a TLDc-domain-containing homolog gene of OXR1, have previously been identified in patients with hearing loss. Interestingly, mutations or knockout of OXR1 interacting molecules such as ATP6V1B1 and ESR1 are also associated with hearing loss in patients or animal models, hinting an important role of OXR1 and associated partners in cochlear development and hearing function.

CXCL8/CXCR2 signaling mediates bone marrow fibrosis and is a therapeutic target in myelofibrosis.

Proinflammatory signaling is a hallmark feature of human cancer, including in myeloproliferative neoplasms (MPNs), most notably myelofibrosis (MF). Dysregulated inflammatory signaling contributes to fibrotic progression in MF; however, the individual cytokine mediators elicited by malignant MPN cells to promote collagen-producing fibrosis and disease evolution are yet to be fully elucidated. Previously, we identified a critical role for combined constitutive JAK/STAT and aberrant NF-κB proinflammatory signaling in MF development. Using single-cell transcriptional and cytokine-secretion studies of primary cells from patients with MF and the human MPLW515L (hMPLW515L) murine model of MF, we extend our previous work and delineate the role of CXCL8/CXCR2 signaling in MF pathogenesis and bone marrow fibrosis progression. Hematopoietic stem/progenitor cells from patients with MF are enriched for a CXCL8/CXCR2 gene signature and display enhanced proliferation and fitness in response to an exogenous CXCL8 ligand in vitro. Genetic deletion of Cxcr2 in the hMPLW515L-adoptive transfer model abrogates fibrosis and extends overall survival, and pharmacologic inhibition of the CXCR1/2 pathway improves hematologic parameters, attenuates bone marrow fibrosis, and synergizes with JAK inhibitor therapy. Our mechanistic insights provide a rationale for therapeutic targeting of the CXCL8/CXCR2 pathway among patients with MF.

Upregulation of proBDNF/p75NTR signaling in immune cells and its correlation with inflammatory markers in patients with major depression.

ProBDNF is the precursor protein of brain-derived neurotrophic factor (BDNF) expressed in the central nervous system and peripheral tissues. Previous studies showed that the blood levels of both proBDNF and p75 neurotrophic receptors (p75NTR) in major depressive disorder (MDD) were increased, but which blood cell types express proBDNF and its receptors is not known. Furthermore, the relationship between proBDNF/p75NTR and inflammatory cytokines in peripheral blood of MDD is unclear. Peripheral blood mononuclear cells (PBMCs) and serum were obtained from depressive patients (n = 32) and normal donors (n = 20). We examined the expression of proBDNF and inflammatory markers and their correlative relationship in patients with major depression. Using flow cytometry analysis, we examined which blood cells express proBDNF and its receptors. Finally, the role of proBDNF/p75NTR signal in inflammatory immune activity of PBMCs was verified in vitro experiments. Inflammatory cytokines in PBMC from MDD patients were increased and correlated with the major depression scores. The levels of IL-1ß and IL-10 were also positively correlated with the major depression scores, while the levels of TNF-α and IL-6 were negatively correlated with the major depression scores. Intriguingly, the levels of sortilin were positively correlated with IL-1ß. Q-PCR and Western blots showed proBDNF, p75NTR, and sortilin levels were significantly increased in PBMCs from MDD patients compared with that from the normal donors. Flow cytometry studies showed that proBDNF and p75NTR were present mainly in CD4+ and CD8+ T cells. The number of proBDNF and p75NTR positive CD4+ and CD8+ T cells from MDD patients was increased and subsequently reversed after therapeutic management. Exogenous proBDNF protein or p75ECD-Fc treatment of cultured PBMC affected the release of inflammatory cytokines in vitro. ProBDNF promoted the expression of inflammatory cytokines, while p75ECD-Fc inhibited the expression of inflammatory cytokines. Given there was an inflammatory response of lymphocytes to proBDNF, it is suggested that proBDNF/p75NTR signaling may upstream inflammatory cytokines in MDD. Our data suggest that proBDNF/p75NTR signaling may not only serve as biomarkers but also may be a potential therapeutic target for MDD.

14-3-3τ drives estrogen receptor loss via ERα36 induction and GATA3 inhibition in breast cancer.

About one-fourth of recurrent estrogen receptor-positive (ER+) breast cancers lose ER expression, leading to endocrine therapy failure. However, the mechanisms underlying ER loss remain to be fully explored. We now show that 14-3-3τ, up-regulated in ∼60% of breast cancer, drives the conversion of ER+ to ER- and epithelial-to-mesenchymal transition (EMT). We identify ERα36, an isoform of ERα66, as a downstream effector of 14-3-3τ. Overexpression of 14-3-3τ induces ERα36 in xenografts and tumor spheroids. The regulation is further supported by a positive correlation between ERα36 and 14-3-3τ expression in human breast cancers. ERα36 can antagonize ERα66 and inhibit ERα66 expression. Isoform-specific depletion of ERα36 blocks the ER conversion and EMT induced by 14-3-3τ overexpression in tumor spheroids, thus establishing ERα36 as a key mediator in 14-3-3τ-driven ER loss and EMT. ERα36 promoter is repressed by GATA3, which can be phosphorylated by AKT at consensus binding sites for 14-3-3. Upon AKT activation, 14-3-3τ binds phosphorylated GATA3 and facilitates the degradation of GATA3 causing GATA3 to lose transcriptional control over its target genes ERα66 and ERα36. We also demonstrate a role for the collaboration between 14-3-3τ and AKT in ERα36 induction and endocrine therapy resistance by three-dimensional spheroid and tamoxifen treatment models in MCF7 and T47D ER+ breast cancer cells. Thus, the 14-3-3τ-ERα36 regulation provides a previously unrecognized mechanism for ER loss and endocrine therapy failure.

Syntheses of three-dimensional catenanes under kinetic control.

Although catenanes comprising two ring-shaped components can be made in large quantities by templation, the preparation of three-dimensional (3D) catenanes with cage-shaped components is still in its infancy. Here, we report the design and syntheses of two 3D catenanes by a sequence of SN2 reactions in one pot. The resulting triply mechanically interlocked molecules were fully characterized in both the solution and solid states. Mechanistic studies have revealed that a suit[3]ane, which contains a threefold symmetric cage component as the suit and a tribromide component as the body, is formed at elevated temperatures. This suit[3]ane was identified as the key reactive intermediate for the selective formation of the two 3D catenanes which do not represent thermodynamic minima. We foresee a future in which this particular synthetic strategy guides the rational design and production of mechanically interlocked molecules under kinetic control.