Cell-type-specific consequences of mosaic structural variants in hematopoietic stem and progenitor cells.

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.

Functional analysis of structural variants in single cells using Strand-seq.

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.

Improved assembly and variant detection of a haploid human genome using single-molecule, high-fidelity long reads.

The sequence and assembly of human genomes using long-read sequencing technologies has revolutionized our understanding of structural variation and genome organization. We compared the accuracy, continuity, and gene annotation of genome assemblies generated from either high-fidelity (HiFi) or continuous long-read (CLR) datasets from the same complete hydatidiform mole human genome. We find that the HiFi sequence data assemble an additional 10% of duplicated regions and more accurately represent the structure of tandem repeats, as validated with orthogonal analyses. As a result, an additional 5 Mbp of pericentromeric sequences are recovered in the HiFi assembly, resulting in a 2.5-fold increase in the NG50 within 1 Mbp of the centromere (HiFi 480.6 kbp, CLR 191.5 kbp). Additionally, the HiFi genome assembly was generated in significantly less time with fewer computational resources than the CLR assembly. Although the HiFi assembly has significantly improved continuity and accuracy in many complex regions of the genome, it still falls short of the assembly of centromeric DNA and the largest regions of segmental duplication using existing assemblers. Despite these shortcomings, our results suggest that HiFi may be the most effective standalone technology for de novo assembly of human genomes.

Single-cell analysis of structural variations and complex rearrangements with tri-channel processing.

Structural variation (SV), involving deletions, duplications, inversions and translocations of DNA segments, is a major source of genetic variability in somatic cells and can dysregulate cancer-related pathways. However, discovering somatic SVs in single cells has been challenging, with copy-number-neutral and complex variants typically escaping detection. Here we describe single-cell tri-channel processing (scTRIP), a computational framework that integrates read depth, template strand and haplotype phase to comprehensively discover SVs in individual cells. We surveyed SV landscapes of 565 single cells, including transformed epithelial cells and patient-derived leukemic samples, to discover abundant SV classes, including inversions, translocations and complex DNA rearrangements. Analysis of the leukemic samples revealed four times more somatic SVs than cytogenetic karyotyping, submicroscopic copy-number alterations, oncogenic copy-neutral rearrangements and a subclonal chromothripsis event. Advancing current methods, single-cell tri-channel processing can directly measure SV mutational processes in individual cells, such as breakage-fusion-bridge cycles, facilitating studies of clonal evolution, genetic mosaicism and SV formation mechanisms, which could improve disease classification for precision medicine.

Direct chromosome-length haplotyping by single-cell sequencing.

Haplotypes are fundamental to fully characterize the diploid genome of an individual, yet methods to directly chart the unique genetic makeup of each parental chromosome are lacking. Here we introduce single-cell DNA template strand sequencing (Strand-seq) as a novel approach to phasing diploid genomes along the entire length of all chromosomes. We demonstrate this by building a complete haplotype for a HapMap individual (NA12878) at high accuracy (concordance 99.3%), without using generational information or statistical inference. By use of this approach, we mapped all meiotic recombination events in a family trio with high resolution (median range ∼14 kb) and phased larger structural variants like deletions, indels, and balanced rearrangements like inversions. Lastly, the single-cell resolution of Strand-seq allowed us to observe loss of heterozygosity regions in a small number of cells, a significant advantage for studies of heterogeneous cell populations, such as cancer cells. We conclude that Strand-seq is a unique and powerful approach to completely phase individual genomes and map inheritance patterns in families, while preserving haplotype differences between single cells.

Characterizing polymorphic inversions in human genomes by single-cell sequencing.

Identifying genomic features that differ between individuals and cells can help uncover the functional variants that drive phenotypes and disease susceptibilities. For this, single-cell studies are paramount, as it becomes increasingly clear that the contribution of rare but functional cellular subpopulations is important for disease prognosis, management, and progression. Until now, studying these associations has been challenged by our inability to map structural rearrangements accurately and comprehensively. To overcome this, we coupled single-cell sequencing of DNA template strands (Strand-seq) with custom analysis software to rapidly discover, map, and genotype genomic rearrangements at high resolution. This allowed us to explore the distribution and frequency of inversions in a heterogeneous cell population, identify several polymorphic domains in complex regions of the genome, and locate rare alleles in the reference assembly. We then mapped the entire genomic complement of inversions within two unrelated individuals to characterize their distinct inversion profiles and built a nonredundant global reference of structural rearrangements in the human genome. The work described here provides a powerful new framework to study structural variation and genomic heterogeneity in single-cell samples, whether from individuals for population studies or tissue types for biomarker discovery.

The Prognostic Impact of CD163-Positive Macrophages in Follicular Lymphoma: A Study from the BC Cancer Agency and the Lymphoma Study Association.

PURPOSE: We aimed to assess the prognostic significance of follicular lymphoma-associated macrophages in the era of rituximab treatment and maintenance. EXPERIMENTAL DESIGN: We applied immunohistochemistry for CD68 and CD163 to two large tissue microarrays (TMA). The first TMA included samples from 186 patients from the BC Cancer Agency (BCCA) who had been treated with first-line systemic treatment including rituximab, cyclophosphamide, vincristine, and prednisone. The second contained 395 samples from PRIMA trial patients treated with rituximab, cyclophosphamide, doxorubicin, vincristine, and prednisone, and randomized to rituximab maintenance or observation. Macrophage infiltration was assessed using Aperio image analysis. Each of the two cohorts was randomly split into training/validation sets. RESULTS: An increased CD163-positive pixel count was predictive of adverse outcome in the BCCA dataset [5-year progression-free survival (PFS) 38% vs. 72%, respectively, P = 0.004 in the training cohort and 5-year PFS 29% vs. 61%, respectively, P = 0.004 in the validation cohort]. In the PRIMA trial, an increased CD163 pixel count was associated with favorable outcome (5-year PFS 60% vs. 44%, respectively, P = 0.011 in the training cohort and 5-year PFS 55% vs. 37%, respectively, P = 0.030 in the validation cohort). CONCLUSIONS: CD163-positive macrophages predict outcome in follicular lymphoma, but their prognostic impact is highly dependent on treatment received.

Essential role for Ptpn11 in survival of hematopoietic stem and progenitor cells.

Src homology 2 domain-containing phosphatase 2 (Shp2), encoded by Ptpn11, is a member of the nonreceptor protein-tyrosine phosphatase family, and functions in cell survival, proliferation, migration, and differentiation in many tissues. Here we report that loss of Ptpn11 in murine hematopoietic cells leads to bone marrow aplasia and lethality. Mutant mice show rapid loss of hematopoietic stem cells (HSCs) and immature progenitors of all hematopoietic lineages in a gene dosage-dependent and cell-autonomous manner. Ptpn11-deficient HSCs and progenitors undergo apoptosis concomitant with increased Noxa expression. Mutant HSCs/progenitors also show defective Erk and Akt activation in response to stem cell factor and diminished thrombopoietin-evoked Erk activation. Activated Kras alleviates the Ptpn11 requirement for colony formation by progenitors and cytokine/growth factor responsiveness of HSCs, indicating that Ras is functionally downstream of Shp2 in these cells. Thus, Shp2 plays a critical role in controlling the survival and maintenance of HSCs and immature progenitors in vivo.