Linking chromatin acylation mark-defined proteome and genome in living cells.

A generalizable strategy with programmable site specificity for in situ profiling of histone modifications on unperturbed chromatin remains highly desirable but challenging. We herein developed a single-site-resolved multi-omics (SiTomics) strategy for systematic mapping of dynamic modifications and subsequent profiling of chromatinized proteome and genome defined by specific chromatin acylations in living cells. By leveraging the genetic code expansion strategy, our SiTomics toolkit revealed distinct crotonylation (e.g., H3K56cr) and ß-hydroxybutyrylation (e.g., H3K56bhb) upon short chain fatty acids stimulation and established linkages for chromatin acylation mark-defined proteome, genome, and functions. This led to the identification of GLYR1 as a distinct interacting protein in modulating H3K56cr's gene body localization as well as the discovery of an elevated super-enhancer repertoire underlying bhb-mediated chromatin modulations. SiTomics offers a platform technology for elucidating the "metabolites-modification-regulation" axis, which is widely applicable for multi-omics profiling and functional dissection of modifications beyond acylations and proteins beyond histones.

Multi-omics analysis of mucosal and systemic immunity to SARS-CoV-2 after birth.

The dynamics of immunity to infection in infants remain obscure. Here, we used a multi-omics approach to perform a longitudinal analysis of immunity to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in infants and young children by analyzing blood samples and weekly nasal swabs collected before, during, and after infection with Omicron and non-Omicron variants. Infection stimulated robust antibody titers that, unlike in adults, showed no sign of decay for up to 300 days. Infants mounted a robust mucosal immune response characterized by inflammatory cytokines, interferon (IFN) α, and T helper (Th) 17 and neutrophil markers (interleukin [IL]-17, IL-8, and CXCL1). The immune response in blood was characterized by upregulation of activation markers on innate cells, no inflammatory cytokines, but several chemokines and IFNα. The latter correlated with viral load and expression of interferon-stimulated genes (ISGs) in myeloid cells measured by single-cell multi-omics. Together, these data provide a snapshot of immunity to infection during the initial weeks and months of life.

Neurulation of the cynomolgus monkey embryo achieved from 3D blastocyst culture.

Neural tube (NT) defects arise from abnormal neurulation and result in the most common birth defects worldwide. Yet, mechanisms of primate neurulation remain largely unknown due to prohibitions on human embryo research and limitations of available model systems. Here, we establish a three-dimensional (3D) prolonged in vitro culture (pIVC) system supporting cynomolgus monkey embryo development from 7 to 25 days post-fertilization. Through single-cell multi-omics analyses, we demonstrate that pIVC embryos form three germ layers, including primordial germ cells, and establish proper DNA methylation and chromatin accessibility through advanced gastrulation stages. In addition, pIVC embryo immunofluorescence confirms neural crest formation, NT closure, and neural progenitor regionalization. Finally, we demonstrate that the transcriptional profiles and morphogenetics of pIVC embryos resemble key features of similarly staged in vivo cynomolgus and human embryos. This work therefore describes a system to study non-human primate embryogenesis through advanced gastrulation and early neurulation.

A blood atlas of COVID-19 defines hallmarks of disease severity and specificity.

Treatment of severe COVID-19 is currently limited by clinical heterogeneity and incomplete description of specific immune biomarkers. We present here a comprehensive multi-omic blood atlas for patients with varying COVID-19 severity in an integrated comparison with influenza and sepsis patients versus healthy volunteers. We identify immune signatures and correlates of host response. Hallmarks of disease severity involved cells, their inflammatory mediators and networks, including progenitor cells and specific myeloid and lymphocyte subsets, features of the immune repertoire, acute phase response, metabolism, and coagulation. Persisting immune activation involving AP-1/p38MAPK was a specific feature of COVID-19. The plasma proteome enabled sub-phenotyping into patient clusters, predictive of severity and outcome. Systems-based integrative analyses including tensor and matrix decomposition of all modalities revealed feature groupings linked with severity and specificity compared to influenza and sepsis. Our approach and blood atlas will support future drug development, clinical trial design, and personalized medicine approaches for COVID-19.

Multiple early factors anticipate post-acute COVID-19 sequelae.

Post-acute sequelae of COVID-19 (PASC) represent an emerging global crisis. However, quantifiable risk factors for PASC and their biological associations are poorly resolved. We executed a deep multi-omic, longitudinal investigation of 309 COVID-19 patients from initial diagnosis to convalescence (2-3 months later), integrated with clinical data and patient-reported symptoms. We resolved four PASC-anticipating risk factors at the time of initial COVID-19 diagnosis: type 2 diabetes, SARS-CoV-2 RNAemia, Epstein-Barr virus viremia, and specific auto-antibodies. In patients with gastrointestinal PASC, SARS-CoV-2-specific and CMV-specific CD8+ T cells exhibited unique dynamics during recovery from COVID-19. Analysis of symptom-associated immunological signatures revealed coordinated immunity polarization into four endotypes, exhibiting divergent acute severity and PASC. We find that immunological associations between PASC factors diminish over time, leading to distinct convalescent immune states. Detectability of most PASC factors at COVID-19 diagnosis emphasizes the importance of early disease measurements for understanding emergent chronic conditions and suggests PASC treatment strategies.

Time-resolved systems immunology reveals a late juncture linked to fatal COVID-19.

COVID-19 exhibits extensive patient-to-patient heterogeneity. To link immune response variation to disease severity and outcome over time, we longitudinally assessed circulating proteins as well as 188 surface protein markers, transcriptome, and T cell receptor sequence simultaneously in single peripheral immune cells from COVID-19 patients. Conditional-independence network analysis revealed primary correlates of disease severity, including gene expression signatures of apoptosis in plasmacytoid dendritic cells and attenuated inflammation but increased fatty acid metabolism in CD56dimCD16hi NK cells linked positively to circulating interleukin (IL)-15. CD8+ T cell activation was apparent without signs of exhaustion. Although cellular inflammation was depressed in severe patients early after hospitalization, it became elevated by days 17-23 post symptom onset, suggestive of a late wave of inflammatory responses. Furthermore, circulating protein trajectories at this time were divergent between and predictive of recovery versus fatal outcomes. Our findings stress the importance of timing in the analysis, clinical monitoring, and therapeutic intervention of COVID-19.

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.

Fundamental Biological Features of Spaceflight: Advancing the Field to Enable Deep-Space Exploration.

Research on astronaut health and model organisms have revealed six features of spaceflight biology that guide our current understanding of fundamental molecular changes that occur during space travel. The features include oxidative stress, DNA damage, mitochondrial dysregulation, epigenetic changes (including gene regulation), telomere length alterations, and microbiome shifts. Here we review the known hazards of human spaceflight, how spaceflight affects living systems through these six fundamental features, and the associated health risks of space exploration. We also discuss the essential issues related to the health and safety of astronauts involved in future missions, especially planned long-duration and Martian missions.

Multi-Omics Resolves a Sharp Disease-State Shift between Mild and Moderate COVID-19.

We present an integrated analysis of the clinical measurements, immune cells, and plasma multi-omics of 139 COVID-19 patients representing all levels of disease severity, from serial blood draws collected during the first week of infection following diagnosis. We identify a major shift between mild and moderate disease, at which point elevated inflammatory signaling is accompanied by the loss of specific classes of metabolites and metabolic processes. Within this stressed plasma environment at moderate disease, multiple unusual immune cell phenotypes emerge and amplify with increasing disease severity. We condensed over 120,000 immune features into a single axis to capture how different immune cell classes coordinate in response to SARS-CoV-2. This immune-response axis independently aligns with the major plasma composition changes, with clinical metrics of blood clotting, and with the sharp transition between mild and moderate disease. This study suggests that moderate disease may provide the most effective setting for therapeutic intervention.

An Organoid Biobank of Neuroendocrine Neoplasms Enables Genotype-Phenotype Mapping.

Gastroenteropancreatic (GEP) neuroendocrine neoplasm (NEN) that consists of neuroendocrine tumor and neuroendocrine carcinoma (NEC) is a lethal but under-investigated disease owing to its rarity. To fill the scarcity of clinically relevant models of GEP-NEN, we here established 25 lines of NEN organoids and performed their comprehensive molecular characterization. GEP-NEN organoids recapitulated pathohistological and functional phenotypes of the original tumors. Whole-genome sequencing revealed frequent genetic alterations in TP53 and RB1 in GEP-NECs, and characteristic chromosome-wide loss of heterozygosity in GEP-NENs. Transcriptome analysis identified molecular subtypes that are distinguished by the expression of distinct transcription factors. GEP-NEN organoids gained independence from the stem cell niche irrespective of genetic mutations. Compound knockout of TP53 and RB1, together with overexpression of key transcription factors, conferred on the normal colonic epithelium phenotypes that are compatible with GEP-NEN biology. Altogether, our study not only provides genetic understanding of GEP-NEN, but also connects its genetics and biological phenotypes.

Multimodal Analysis of Composition and Spatial Architecture in Human Squamous Cell Carcinoma.

To define the cellular composition and architecture of cutaneous squamous cell carcinoma (cSCC), we combined single-cell RNA sequencing with spatial transcriptomics and multiplexed ion beam imaging from a series of human cSCCs and matched normal skin. cSCC exhibited four tumor subpopulations, three recapitulating normal epidermal states, and a tumor-specific keratinocyte (TSK) population unique to cancer, which localized to a fibrovascular niche. Integration of single-cell and spatial data mapped ligand-receptor networks to specific cell types, revealing TSK cells as a hub for intercellular communication. Multiple features of potential immunosuppression were observed, including T regulatory cell (Treg) co-localization with CD8 T cells in compartmentalized tumor stroma. Finally, single-cell characterization of human tumor xenografts and in vivo CRISPR screens identified essential roles for specific tumor subpopulation-enriched gene networks in tumorigenesis. These data define cSCC tumor and stromal cell subpopulations, the spatial niches where they interact, and the communicating gene networks that they engage in cancer.

Molecular Choreography of Acute Exercise.

Acute physical activity leads to several changes in metabolic, cardiovascular, and immune pathways. Although studies have examined selected changes in these pathways, the system-wide molecular response to an acute bout of exercise has not been fully characterized. We performed longitudinal multi-omic profiling of plasma and peripheral blood mononuclear cells including metabolome, lipidome, immunome, proteome, and transcriptome from 36 well-characterized volunteers, before and after a controlled bout of symptom-limited exercise. Time-series analysis revealed thousands of molecular changes and an orchestrated choreography of biological processes involving energy metabolism, oxidative stress, inflammation, tissue repair, and growth factor response, as well as regulatory pathways. Most of these processes were dampened and some were reversed in insulin-resistant participants. Finally, we discovered biological pathways involved in cardiopulmonary exercise response and developed prediction models revealing potential resting blood-based biomarkers of peak oxygen consumption.

Rewiring of the Fruit Metabolome in Tomato Breeding.

Humans heavily rely on dozens of domesticated plant species that have been further improved through intensive breeding. To evaluate how breeding changed the tomato fruit metabolome, we have generated and analyzed a dataset encompassing genomes, transcriptomes, and metabolomes from hundreds of tomato genotypes. The combined results illustrate how breeding globally altered fruit metabolite content. Selection for alleles of genes associated with larger fruits altered metabolite profiles as a consequence of linkage with nearby genes. Selection of five major loci reduced the accumulation of anti-nutritional steroidal glycoalkaloids in ripened fruits, rendering the fruit more edible. Breeding for pink tomatoes modified the content of over 100 metabolites. The introgression of resistance genes from wild relatives in cultivars also resulted in major and unexpected metabolic changes. The study reveals a multi-omics view of the metabolic breeding history of tomato, as well as provides insights into metabolome-assisted breeding and plant biology.

New horizons in the stormy sea of multimodal single-cell data integration.

While measurements of RNA expression have dominated the world of single-cell analyses, new single-cell techniques increasingly allow collection of different data modalities, measuring different molecules, structural connections, and intermolecular interactions. Integrating the resulting multimodal single-cell datasets is a new bioinformatics challenge. Equally important, it is a new experimental design challenge for the bench scientist, who is not only choosing from a myriad of techniques for each data modality but also faces new challenges in experimental design. The ultimate goal is to design, execute, and analyze multimodal single-cell experiments that are more than just descriptive but enable the learning of new causal and mechanistic biology. This objective requires strict consideration of the goals behind the analysis, which might range from mapping the heterogeneity of a cellular population to assembling system-wide causal networks that can further our understanding of cellular functions and eventually lead to models of tissues and organs. We review steps and challenges toward this goal. Single-cell transcriptomics is now a mature technology, and methods to measure proteins, lipids, small-molecule metabolites, and other molecular phenotypes at the single-cell level are rapidly developing. Integrating these single-cell readouts so that each cell has measurements of multiple types of data, e.g., transcriptomes, proteomes, and metabolomes, is expected to allow identification of highly specific cellular subpopulations and to provide the basis for inferring causal biological mechanisms.

Single-Cell Analyses Reveal Megakaryocyte-Biased Hematopoiesis in Myelofibrosis and Identify Mutant Clone-Specific Targets.

Myelofibrosis is a severe myeloproliferative neoplasm characterized by increased numbers of abnormal bone marrow megakaryocytes that induce fibrosis, destroying the hematopoietic microenvironment. To determine the cellular and molecular basis for aberrant megakaryopoiesis in myelofibrosis, we performed single-cell transcriptome profiling of 135,929 CD34+ lineage- hematopoietic stem and progenitor cells (HSPCs), single-cell proteomics, genomics, and functional assays. We identified a bias toward megakaryocyte differentiation apparent from early multipotent stem cells in myelofibrosis and associated aberrant molecular signatures. A sub-fraction of myelofibrosis megakaryocyte progenitors (MkPs) are transcriptionally similar to healthy-donor MkPs, but the majority are disease specific, with distinct populations expressing fibrosis- and proliferation-associated genes. Mutant-clone HSPCs have increased expression of megakaryocyte-associated genes compared to wild-type HSPCs, and we provide early validation of G6B as a potential immunotherapy target. Our study paves the way for selective targeting of the myelofibrosis clone and illustrates the power of single-cell multi-omics to discover tumor-specific therapeutic targets and mediators of tissue fibrosis.

A Highly Scalable Method for Joint Whole-Genome Sequencing and Gene-Expression Profiling of Single Cells.

To address how genetic variation alters gene expression in complex cell mixtures, we developed direct nuclear tagmentation and RNA sequencing (DNTR-seq), which enables whole-genome and mRNA sequencing jointly in single cells. DNTR-seq readily identified minor subclones within leukemia patients. In a large-scale DNA damage screen, DNTR-seq was used to detect regions under purifying selection and identified genes where mRNA abundance was resistant to copy-number alteration, suggesting strong genetic compensation. mRNA sequencing (mRNA-seq) quality equals RNA-only methods, and the low positional bias of genomic libraries allowed detection of sub-megabase aberrations at ultra-low coverage. Each cell library is individually addressable and can be re-sequenced at increased depth, allowing multi-tiered study designs. Additionally, the direct tagmentation protocol enables coverage-independent estimation of ploidy, which can be used to identify cell singlets. Thus, DNTR-seq directly links each cell's state to its corresponding genome at scale, enabling routine analysis of heterogeneous tumors and other complex tissues.

Deciphering dynamic interactions between spermatozoa and the ovarian microenvironment through integrated multi-omics approaches in viviparous Sebastes schlegelii.

The ovarian microenvironment plays a crucial role in ensuring the reproductive success of viviparous teleosts. However, the molecular mechanism underlying the interaction between spermatozoa and the ovarian microenvironment has remained elusive. This study aimed to contribute to a better understanding of this process in black rockfish (Sebastes schlegelii) using integrated multi-omics approaches. The results demonstrated significant upregulation of ovarian complement-related proteins and pattern recognition receptors, along with remodeling of glycans on the surface of spermatozoa at the early spermatozoa-storage stage (1 month after mating). As spermatozoa were stored over time, ovarian complement proteins were progressively repressed by tryptophan and hippurate, indicating a remarkable adaptation of spermatozoa to the ovarian microenvironment. Before fertilization, a notable upregulation of cellular junction proteins was observed. The study revealed that spermatozoa bind to ZPB2a protein through GSTM3 and that ZPB2a promotes spermatozoa survival and movement in a GSTM3-dependent manner. These findings shed light on a key mechanism that influences the dynamics of spermatozoa in the female reproductive tract, providing valuable insights into the molecular networks regulating spermatozoa adaptation and survival in species with internal fertilization.

Open-source, high-throughput targeted in situ transcriptomics for developmental and tissue biology.

Multiplexed spatial profiling of mRNAs has recently gained traction as a tool to explore the cellular diversity and the architecture of tissues. We propose a sensitive, open-source, simple and flexible method for the generation of in situ expression maps of hundreds of genes. We use direct ligation of padlock probes on mRNAs, coupled with rolling circle amplification and hybridization-based in situ combinatorial barcoding, to achieve high detection efficiency, high-throughput and large multiplexing. We validate the method across a number of species and show its use in combination with orthogonal methods such as antibody staining, highlighting its potential value for developmental and tissue biology studies. Finally, we provide an end-to-end computational workflow that covers the steps of probe design, image processing, data extraction, cell segmentation, clustering and annotation of cell types. By enabling easier access to high-throughput spatially resolved transcriptomics, we hope to encourage a diversity of applications and the exploration of a wide range of biological questions.

Melanocyte lineage dynamics in development, growth and disease.

Melanocytes evolved to produce the melanin that gives colour to our hair, eyes and skin. The melanocyte lineage also gives rise to melanoma, the most lethal form of skin cancer. The melanocyte lineage differentiates from neural crest cells during development, and most melanocytes reside in the skin and hair, where they are replenished by melanocyte stem cells. Because the molecular mechanisms necessary for melanocyte specification, migration, proliferation and differentiation are co-opted during melanoma initiation and progression, studying melanocyte development is directly relevant to human disease. Here, through the lens of advances in cellular omic and genomic technologies, we review the latest findings in melanocyte development and differentiation, and how these developmental pathways become dysregulated in disease.

Systems analysis of human innate immunity in COVID-19.

Recent developments in sequencing technologies, the computer and data sciences, as well as increasingly high-throughput immunological measurements have made it possible to derive holistic views on pathophysiological processes of disease and treatment effects directly in humans. We and others have illustrated that incredibly predictive data for immune cell function can be generated by single cell multi-omics (SCMO) technologies and that these technologies are perfectly suited to dissect pathophysiological processes in a new disease such as COVID-19, triggered by SARS-CoV-2 infection. Systems level interrogation not only revealed the different disease endotypes, highlighted the differential dynamics in context of disease severity, and pointed towards global immune deviation across the different arms of the immune system, but was already instrumental to better define long COVID phenotypes, suggest promising biomarkers for disease and therapy outcome predictions and explains treatment responses for the widely used corticosteroids. As we identified SCMO to be the most informative technologies in the vest to better understand COVID-19, we propose to routinely include such single cell level analysis in all future clinical trials and cohorts addressing diseases with an immunological component.