Single-cell genome and transcriptome sequencing without upfront whole-genome amplification reveals cell state plasticity of melanoma subclones.

Single-cell multi-omics methods enable the study of cell state diversity, which is largely determined by the interplay of the genome, epigenome, and transcriptome. Here, we describe Gtag&T-seq, a genome-and-transcriptome sequencing (G&T-seq) protocol of the same single cells that omits whole-genome amplification (WGA) by using direct genomic tagmentation (Gtag). Gtag drastically decreases the cost and improves coverage uniformity at single-cell and pseudo-bulk levels compared to WGA-based G&T-seq.

Crowd-sourced benchmarking of single-sample tumor subclonal reconstruction.

Subclonal reconstruction algorithms use bulk DNA sequencing data to quantify parameters of tumor evolution, allowing an assessment of how cancers initiate, progress and respond to selective pressures. We launched the ICGC-TCGA (International Cancer Genome Consortium-The Cancer Genome Atlas) DREAM Somatic Mutation Calling Tumor Heterogeneity and Evolution Challenge to benchmark existing subclonal reconstruction algorithms. This 7-year community effort used cloud computing to benchmark 31 subclonal reconstruction algorithms on 51 simulated tumors.

Single-cell spatial multi-omics and deep learning dissect enhancer-driven gene regulatory networks in liver zonation.

In the mammalian liver, hepatocytes exhibit diverse metabolic and functional profiles based on their location within the liver lobule. However, it is unclear whether this spatial variation, called zonation, is governed by a well-defined gene regulatory code. Here, using a combination of single-cell multiomics, spatial omics, massively parallel reporter assays and deep learning, we mapped enhancer-gene regulatory networks across mouse liver cell types.

System-wide mapping of peptide-GPCR interactions in C. elegans.

Neuropeptides and peptide hormones are ancient, widespread signaling molecules that underpin almost all brain functions. They constitute a broad ligand-receptor network, mainly by binding to G protein-coupled receptors (GPCRs). However, the organization of the peptidergic network and roles of many peptides remain elusive, as our insight into peptide-receptor interactions is limited and many peptide GPCRs are still orphan receptors.

Single-cell somatic copy number variants in brain using different amplification methods and reference genomes.

The presence of somatic mutations, including copy number variants (CNVs), in the brain is well recognized. Comprehensive study requires single-cell whole genome amplification, with several methods available, prior to sequencing. We compared PicoPLEX with two recent adaptations of multiple displacement amplification (MDA): primary template-directed amplification (PTA) and droplet MDA, across 93 human brain cortical nuclei.

The Clinicogenomic Landscape of Induction Failure in Childhood and Young Adult T-Cell Acute Lymphoblastic Leukemia.

Purpose: Failure to respond to induction chemotherapy portends a poor outcome in childhood acute lymphoblastic leukemia (ALL) and is more frequent in T-cell ALL (T-ALL) than B-cell ALL. We aimed to address the limited understanding of clinical and genetic factors that influence outcome in a cohort of patients with T-ALL induction failure (IF).

Methods: We studied all cases of T-ALL IF on two consecutive multinational randomized trials, UKALL2003 and UKALL2011, to define risk factors, treatment, and outcomes.

Evolutionary characterization of lung adenocarcinoma morphology in TRACERx.

Lung adenocarcinomas (LUADs) display a broad histological spectrum from low-grade lepidic tumors through to mid-grade acinar and papillary and high-grade solid, cribriform and micropapillary tumors. How morphology reflects tumor evolution and disease progression is poorly understood. Whole-exome sequencing data generated from 805 primary tumor regions and 121 paired metastatic samples across 248 LUADs from the TRACERx 421 cohort, together with RNA-sequencing data from 463 primary tumor regions, were integrated with detailed whole-tumor and regional histopathological analysis.

Estimation of tumor cell total mRNA expression in 15 cancer types predicts disease progression.

Single-cell RNA sequencing studies have suggested that total mRNA content correlates with tumor phenotypes. Technical and analytical challenges, however, have so far impeded at-scale pan-cancer examination of total mRNA content. Here we present a method to quantify tumor-specific total mRNA expression (TmS) from bulk sequencing data, taking into account tumor transcript proportion, purity and ploidy, which are estimated through transcriptomic/genomic deconvolution.

Biallelic mutations in cancer genomes reveal local mutational determinants.

The infinite sites model of molecular evolution posits that every position in the genome is mutated at most once. By restricting the number of possible mutation histories, haplotypes and alleles, it forms a cornerstone of tumor phylogenetic analysis and is often implied when calling, phasing and interpreting variants or studying the mutational landscape as a whole. Here we identify 18,295 biallelic mutations, where the same base is mutated independently on both parental copies, in 559 (21%) bulk sequencing samples from the Pan-Cancer Analysis of Whole Genomes study.

Analysis of long and short enhancers in melanoma cell states.

Understanding how enhancers drive cell-type specificity and efficiently identifying them is essential for the development of innovative therapeutic strategies. In melanoma, the melanocytic (MEL) and the mesenchymal-like (MES) states present themselves with different responses to therapy, making the identification of specific enhancers highly relevant. Using massively parallel reporter assays (MPRAs) in a panel of patient-derived melanoma lines (MM lines), we set to identify and decipher melanoma enhancers by first focusing on regions with state-specific H3K27 acetylation close to differentially expressed genes.

A ubiquitous disordered protein interaction module orchestrates transcription elongation.

During eukaryotic transcription elongation, RNA polymerase II (RNAP2) is regulated by a chorus of factors. Here, we identified a common binary interaction module consisting of TFIIS N-terminal domains (TNDs) and natively unstructured TND-interacting motifs (TIMs). This module was conserved among the elongation machinery and linked complexes including transcription factor TFIIS, Mediator, super elongation complex, elongin, IWS1, SPT6, PP1-PNUTS phosphatase, H3K36me3 readers, and other factors.

A pan-cancer landscape of somatic mutations in non-unique regions of the human genome.

A substantial fraction of the human genome displays high sequence similarity with at least one other genomic sequence, posing a challenge for the identification of somatic mutations from short-read sequencing data. Here we annotate genomic variants in 2,658 cancers from the Pan-Cancer Analysis of Whole Genomes (PCAWG) cohort with links to similar sites across the human genome. We train a machine learning model to use signals distributed over multiple genomic sites to call somatic events in non-unique regions and validate the data against linked-read sequencing in an independent dataset.

Interpretation of allele-specific chromatin accessibility using cell state-aware deep learning.

Genomic sequence variation within enhancers and promoters can have a significant impact on the cellular state and phenotype. However, sifting through the millions of candidate variants in a personal genome or a cancer genome, to identify those that impact -regulatory function, remains a major challenge. Interpretation of noncoding genome variation benefits from explainable artificial intelligence to predict and interpret the impact of a mutation on gene regulation.

Transcriptional Profiling of STAT1 Gain-of-Function Reveals Common and Mutation-Specific Fingerprints.

STAT1 gain-of-function (GOF) is a primary immunodeficiency typically characterized by chronic mucocutaneous candidiasis (CMC), recurrent respiratory infections, and autoimmunity. Less commonly, also immunodysregulation polyendocrinopathy enteropathy X-linked (IPEX)-like syndromes with CMC, and combined immunodeficiency without CMC have been described. Recently, our group and others have shown that different mutation-specific mechanisms underlie STAT1 GOF , including faster nuclear accumulation (R274W), and reduced mobility (R321, N574I) to near immobility in the nucleus (T419R) upon IFNγ stimulation.