Chromosomal instability is a major driver of intratumoral heterogeneity (ITH), promoting tumor progression. In the present study, we combined structural variant discovery and nucleosome occupancy profiling with transcriptomic and immunophenotypic changes in single cells to study ITH in complex karyotype acute myeloid leukemia (CK-AML). We observed complex structural variant landscapes within individual cells of patients with CK-AML characterized by linear and circular breakage–fusion–bridge cycles and chromothripsis. We identified three clonal evolution patterns in diagnosis or salvage CK-AML (monoclonal, linear and branched polyclonal), with 75% harboring multiple subclones that frequently displayed ongoing karyotype remodeling. Using patient-derived xenografts, we demonstrated varied clonal evolution of leukemic stem cells (LSCs) and further dissected subclone-specific drug–response profiles to identify LSC-targeting therapies, including BCL-xL inhibition. In paired longitudinal patient samples, we further revealed genetic evolution and cell-type plasticity as mechanisms of disease progression. By dissecting dynamic genomic, phenotypic and functional complexity of CK-AML, our findings offer clinically relevant avenues for characterizing and targeting disease-driving LSCs.

Somatic structural variants (SVs) are widespread in cancer, but their impact on disease evolution is understudied due to a lack of methods to directly characterize their functional consequences. We present a computational method, scNOVA, which uses Strand-seq to perform haplotype-aware integration of SV discovery and molecular phenotyping in single cells by using nucleosome occupancy to infer gene expression as a readout. Application to leukemias and cell lines identifies local effects of copy-balanced rearrangements on gene deregulation, and consequences of SVs on aberrant signaling pathways in subclones. We discovered distinct SV subclones with dysregulated Wnt signaling in a chronic lymphocytic leukemia patient. We further uncovered the consequences of subclonal chromothripsis in T cell acute lymphoblastic leukemia, which revealed c-Myb activation, enrichment of a primitive cell state and informed successful targeting of the subclone in cell culture, using a Notch inhibitor. By directly linking SVs to their functional effects, scNOVA enables systematic single-cell multiomic studies of structural variation in heterogeneous cell populations.

The functional impact and cellular context of mosaic structural variants (mSVs) in normal tissues is understudied. Utilizing Strand-seq, we sequenced 1,133 single-cell genomes from 19 human donors of increasing age, and discovered the heterogeneous mSV landscapes of hematopoietic stem and progenitor cells. While mSVs are continuously acquired throughout life, expanded subclones in our cohort are confined to individuals >60. Cells already harboring mSVs are more likely to acquire additional somatic structural variants, including megabase-scale segmental aneuploidies. Capitalizing on comprehensive single-cell micrococcal nuclease digestion with sequencing reference data, we conducted high-resolution cell-typing for eight hematopoietic stem and progenitor cells. Clonally expanded mSVs disrupt normal cellular function by dysregulating diverse cellular pathways, and enriching for myeloid progenitors. Our findings underscore the contribution of mSVs to the cellular and molecular phenotypes associated with the aging hematopoietic system, and establish a foundation for deciphering the molecular links between mSVs, aging and disease susceptibility in normal tissues.

Leaf senescence, the last stage of leaf development, is essential for whole-plant fitness as it marks the relocation of nutrients from senescing leaves to reproductive or other developing organs. Temporally coordinated physiological and functional changes along leaf aging are fine-tuned by a highly regulated genetic program involving multi-layered regulatory mechanisms. Long noncoding RNAs (lncRNAs) are newly emerging as hidden players in many biological processes; however, their contribution to leaf senescence has been largely unknown. Here, we performed comprehensive analyses of RNA-seq data representing all developmental stages of leaves to determine the genome-wide lncRNA landscape along leaf aging. A total of 771 lncRNAs, including 232 unannotated lncRNAs, were identified. Time-course analysis revealed 446 among 771 developmental age-related lncRNAs (AR-lncRNAs). Intriguingly, the expression of AR-lncRNAs was regulated more dynamically in senescing leaves than in growing leaves, revealing the relevant contribution of these lncRNAs to leaf senescence. Further analyses enabled us to infer the function of lncRNAs, based on their interacting miRNA or mRNA partners. We considered functionally diverse lncRNAs including antisense lncRNAs (which regulate overlapping protein-coding genes), competitive endogenous RNAs (ceRNAs; which regulate paired mRNAs using miRNAs as anchors), and mRNA-interacting lncRNAs (which affect the stability of mRNAs). Furthermore, we experimentally validated the senescence regulatory function of three novel AR-lncRNAs including one antisense lncRNA and two mRNA-interacting lncRNAs through molecular and phenotypic analyses. Our study provides a valuable resource of AR-lncRNAs and potential regulatory networks that link the function of coding mRNA and AR-lncRNAs. Together, our results reveal AR-lncRNAs as important elements in the leaf senescence process.

Naïve CD4 T cells are traditionally viewed as a quiescent, homogeneous, resting population, but emerging evidence reveals their heterogeneity, which can be crucial for understanding disease contexts and therapeutic outcomes. In this study, we identify distinct subpopulations within both murine and human naïve CD4 T cells by single cell-RNA-sequencing (scRNA-seq), particularly focusing on a subpopulation that expresses super-high levels of interleukin-7 receptor (IL-7R^sup-hi), along with CD97, IL-18R, and Ly6C. This subpopulation, absent in the thymus and peripherally induced, exhibits type 1 helper T cell (Th1)-poised characteristics and contributes to the inhibition of cancer progression in B16F10 tumor-bearing mice. In humans, this IL-7R^sup-hi subpopulation expressing CD97 correlates with the responsiveness to anti-PD-1 therapy in cancer patients and the disease state of multiple sclerosis. By elucidating the heterogeneity of naive CD4 T cells and identifying a Th1-poised subpopulation capable of robust type 1 responses, we highlight the importance of this heterogeneity in inflammatory conditions for defining the disease states and predicting drug responsiveness.

Time-of-flight secondary ion mass spectrometry (TOF-SIMS) emerges as a promising tool to identify the ions (small molecules) indicative of disease states from the surface of patient tissues. In TOF-SIMS analysis, an enhanced ionization of surface molecules is critical to increase the number of detected ions. Several methods have been developed to enhance ionization capability. However, how these methods improve identification of disease-related ions has not been systematically explored. Here, we present a multi-dimensional SIMS (MD-SIMS) that combines conventional TOF-SIMS and metal-assisted SIMS (MetA-SIMS). Using this approach, we analyzed cancer and adjacent normal tissues first by TOF-SIMS and subsequently by MetA-SIMS. In total, TOF- and MetA-SIMS detected 632 and 959 ions, respectively. Among them, 426 were commonly detected by both methods, while 206 and 533 were detected uniquely by TOF- and MetA-SIMS, respectively. Of the 426 commonly detected ions, 250 increased in their intensities by MetA-SIMS, whereas 176 decreased. The integrated analysis of the ions detected by the two methods resulted in an increased number of discriminatory ions leading to an enhanced separation between cancer and normal tissues. Therefore, the results show that MD-SIMS can be a useful approach to provide a comprehensive list of discriminatory ions indicative of disease states.