Multimodal single-cell and whole-genome sequencing of small, frozen clinical specimens

Single-cell genomics enables dissection of tumor heterogeneity and molecular underpinnings of drug response at an unprecedented resolution1-11. However, broad clinical application of these methods remains challenging, due to several practical and preanalytical challenges that are incompatible with typical clinical care workflows, namely the need for relatively large, fresh tissue inputs. In the present study, we show that multimodal, single-nucleus (sn)RNA/T cell receptor (TCR) sequencing, spatial transcriptomics and whole-genome sequencing (WGS) are feasible from small, frozen tissues that approximate routinely collected clinical specimens (for example, core needle biopsies). Compared with data from sample-matched fresh tissue, we find a similar quality in the biological outputs of snRNA/TCR-seq data, while reducing artifactual signals and compositional biases introduced by fresh tissue processing. Profiling sequentially collected melanoma samples from a patient treated in the KEYNOTE-001 trial12, we resolved cellular, genomic, spatial and clonotype dynamics that represent molecular patterns of heterogeneous intralesional evolution during anti-programmed cell death protein 1 therapy. To demonstrate applicability to banked biospecimens of rare diseases13, we generated a single-cell atlas of uveal melanoma liver metastasis with matched WGS data. These results show that single-cell genomics from archival, clinical specimens is feasible and provides a framework for translating these methods more broadly to the clinical arena

UV irradiation induces a postreplication DNA damage checkpoint

Eukaryotic cells irradiated with high doses of UV exhibit cell-cycle responses referred to as G(1)/S, intraS, and G(2)/M checkpoints. After a moderate UV dose that approximates sunlight exposure and is lethal to fission yeast checkpoint mutants, we found unexpectedly that these cell-cycle responses do not occur. Instead, cells at all stages of the cell cycle carry lesions into S phase and delay cell-cycle progression for hours after the completion of bulk DNA synthesis. Both DNA replication and the checkpoint kinase, Chk1, are required to generate this cell-cycle response. UV-irradiation of Δchk1 cells causes chromosome damage and loss of viability only after cells have replicated irradiated DNA and entered mitosis. These data suggest that an important physiological role of the cell-cycle response to UV is to provide time for postreplication repair

Dissecting the treatment-naive ecosystem of human melanoma brain metastasis

Melanoma brain metastasis (MBM) frequently occurs in patients with advanced melanoma; yet, our understanding of the underlying salient biology is rudimentary. Here, we performed single-cell/nucleus RNA-seq in 22 treatment-naive MBMs and 10 extracranial melanoma metastases (ECMs) and matched spatial single-cell transcriptomics and T cell receptor (TCR)-seq. Cancer cells from MBM were more chromosomally unstable, adopted a neuronal-like cell state, and enriched for spatially variably expressed metabolic pathways. Key observations were validated in independent patient cohorts, patient-derived MBM/ECM xenograft models, RNA/ATAC-seq, proteomics, and multiplexed imaging. Integrated spatial analyses revealed distinct geography of putative cancer immune evasion and evidence for more abundant intra-tumoral B to plasma cell differentiation in lymphoid aggregates in MBM. MBM harbored larger fractions of monocyte-derived macrophages and dysfunctional TOX+CD8+ T cells with distinct expression of immune checkpoints. This work provides comprehensive insights into MBM biology and serves as a foundational resource for further discovery and therapeutic exploration

A continental phenology model for monitoring vegetation responses to interannual climatic variability

Regional phenology is important in ecosystem simulation models and coupled biosphere/atmosphere models. In the continental United States, the timing of the onset of greenness in the spring (leaf expansion, grass green-up) and offset of greenness in the fall (leaf abscission, cessation of height growth, grass brown-oft) are strongly influenced by meteorological and climatological conditions. We developed predictive phenology models based on traditional phenology research using commonly available meteorological and climatological data. Predictions were compared with satellite phenology observations at numerous 20 km x 20 km contiguous landcover sites. Onset mean absolute error was 7.2 days in the deciduous broadleaf forest (DBF) biome and 6.1 days in the grassland biome. Offset mean absolute error was 5.3 days in the DBF biome and 6.3 days in the grassland biome. Maximum expected errors at a 95 % probability level ranged from 10 to 14 days. Onset was strongly associated with temperature summations in both grassland and DBF biomes; DBF offset was best predicted with a photoperiod function, while grassland offset required a combination of precipitation and temperature controls. A long-term regional test of the DBF onset model captured field-measured interannual variability trends in lila

Demographic history and rare allele sharing among human populations

High-throughput sequencing technology enables population-level surveys of human genomic variation. Here, we examine the joint allele frequency distributions across continental human populations and present an approach for combining complementary aspects of whole-genome, low-coverage data and targeted high-coverage data. We apply this approach to data generated by the pilot phase of the Thousand Genomes Project, including whole-genome 2–4× coverage data for 179 samples from HapMap European, Asian, and African panels as well as high-coverage target sequencing of the exons of 800 genes from 697 individuals in seven populations. We use the site frequency spectra obtained from these data to infer demographic parameters for an Out-of-Africa model for populations of African, European, and Asian descent and to predict, by a jackknife-based approach, the amount of genetic diversity that will be discovered as sample sizes are increased. We predict that the number of discovered nonsynonymous coding variants will reach 100,000 in each population after ∼1,000 sequenced chromosomes per population, whereas ∼2,500 chromosomes will be needed for the same number of synonymous variants. Beyond this point, the number of segregating sites in the European and Asian panel populations is expected to overcome that of the African panel because of faster recent population growth. Overall, we find that the majority of human genomic variable sites are rare and exhibit little sharing among diverged populations. Our results emphasize that replication of disease association for specific rare genetic variants across diverged populations must overcome both reduced statistical power because of rarity and higher population divergence