Genomic characterization of metastatic patterns from prospective clinical sequencing of 25,000 patients.

Metastatic progression is the main cause of death in cancer patients, whereas the underlying genomic mechanisms driving metastasis remain largely unknown. Here, we assembled MSK-MET, a pan-cancer cohort of over 25,000 patients with metastatic diseases. By analyzing genomic and clinical data from this cohort, we identified associations between genomic alterations and patterns of metastatic dissemination across 50 tumor types. We found that chromosomal instability is strongly correlated with metastatic burden in some tumor types, including prostate adenocarcinoma, lung adenocarcinoma, and HR+/HER2+ breast ductal carcinoma, but not in others, including colorectal cancer and high-grade serous ovarian cancer, where copy-number alteration patterns may be established early in tumor development. We also identified somatic alterations associated with metastatic burden and specific target organs. Our data offer a valuable resource for the investigation of the biological basis for metastatic spread and highlight the complex role of chromosomal instability in cancer progression.

Transition to invasive breast cancer is associated with progressive changes in the structure and composition of tumor stroma.

Ductal carcinoma in situ (DCIS) is a pre-invasive lesion that is thought to be a precursor to invasive breast cancer (IBC). To understand the changes in the tumor microenvironment (TME) accompanying transition to IBC, we used multiplexed ion beam imaging by time of flight (MIBI-TOF) and a 37-plex antibody staining panel to interrogate 79 clinically annotated surgical resections using machine learning tools for cell segmentation, pixel-based clustering, and object morphometrics. Comparison of normal breast with patient-matched DCIS and IBC revealed coordinated transitions between four TME states that were delineated based on the location and function of myoepithelium, fibroblasts, and immune cells. Surprisingly, myoepithelial disruption was more advanced in DCIS patients that did not develop IBC, suggesting this process could be protective against recurrence. Taken together, this HTAN Breast PreCancer Atlas study offers insight into drivers of IBC relapse and emphasizes the importance of the TME in regulating these processes.

Discovering dominant tumor immune archetypes in a pan-cancer census.

Cancers display significant heterogeneity with respect to tissue of origin, driver mutations, and other features of the surrounding tissue. It is likely that individual tumors engage common patterns of the immune system-here "archetypes"-creating prototypical non-destructive tumor immune microenvironments (TMEs) and modulating tumor-targeting. To discover the dominant immune system archetypes, the University of California, San Francisco (UCSF) Immunoprofiler Initiative (IPI) processed 364 individual tumors across 12 cancer types using standardized protocols. Computational clustering of flow cytometry and transcriptomic data obtained from cell sub-compartments uncovered dominant patterns of immune composition across cancers. These archetypes were profound insofar as they also differentiated tumors based upon unique immune and tumor gene-expression patterns. They also partitioned well-established classifications of tumor biology. The IPI resource provides a template for understanding cancer immunity as a collection of dominant patterns of immune organization and provides a rational path forward to learn how to modulate these to improve therapy.

Limb development genes underlie variation in human fingerprint patterns.

Fingerprints are of long-standing practical and cultural interest, but little is known about the mechanisms that underlie their variation. Using genome-wide scans in Han Chinese cohorts, we identified 18 loci associated with fingerprint type across the digits, including a genetic basis for the long-recognized "pattern-block" correlations among the middle three digits. In particular, we identified a variant near EVI1 that alters regulatory activity and established a role for EVI1 in dermatoglyph patterning in mice. Dynamic EVI1 expression during human development supports its role in shaping the limbs and digits, rather than influencing skin patterning directly. Trans-ethnic meta-analysis identified 43 fingerprint-associated loci, with nearby genes being strongly enriched for general limb development pathways. We also found that fingerprint patterns were genetically correlated with hand proportions. Taken together, these findings support the key role of limb development genes in influencing the outcome of fingerprint patterning.

SARS-CoV-2 mRNA vaccination elicits a robust and persistent T follicular helper cell response in humans.

SARS-CoV-2 mRNA vaccines induce robust anti-spike (S) antibody and CD4+ T cell responses. It is not yet clear whether vaccine-induced follicular helper CD4+ T (TFH) cell responses contribute to this outstanding immunogenicity. Using fine-needle aspiration of draining axillary lymph nodes from individuals who received the BNT162b2 mRNA vaccine, we evaluated the T cell receptor sequences and phenotype of lymph node TFH. Mining of the responding TFH T cell receptor repertoire revealed a strikingly immunodominant HLA-DPB1∗04-restricted response to S167-180 in individuals with this allele, which is among the most common HLA alleles in humans. Paired blood and lymph node specimens show that while circulating S-specific TFH cells peak one week after the second immunization, S-specific TFH persist at nearly constant frequencies for at least six months. Collectively, our results underscore the key role that robust TFH cell responses play in establishing long-term immunity by this efficacious human vaccine.

Developmental and temporal characteristics of clonal sperm mosaicism.

Throughout development and aging, human cells accumulate mutations resulting in genomic mosaicism and genetic diversity at the cellular level. Mosaic mutations present in the gonads can affect both the individual and the offspring and subsequent generations. Here, we explore patterns and temporal stability of clonal mosaic mutations in male gonads by sequencing ejaculated sperm. Through 300× whole-genome sequencing of blood and sperm from healthy men, we find each ejaculate carries on average 33.3 ± 12.1 (mean ± SD) clonal mosaic variants, nearly all of which are detected in serial sampling, with the majority absent from sampled somal tissues. Their temporal stability and mutational signature suggest origins during embryonic development from a largely immutable stem cell niche. Clonal mosaicism likely contributes a transmissible, predicted pathogenic exonic variant for 1 in 15 men, representing a life-long threat of transmission for these individuals and a significant burden on human population health.

COVID-19 immune features revealed by a large-scale single-cell transcriptome atlas.

A dysfunctional immune response in coronavirus disease 2019 (COVID-19) patients is a recurrent theme impacting symptoms and mortality, yet a detailed understanding of pertinent immune cells is not complete. We applied single-cell RNA sequencing to 284 samples from 196 COVID-19 patients and controls and created a comprehensive immune landscape with 1.46 million cells. The large dataset enabled us to identify that different peripheral immune subtype changes are associated with distinct clinical features, including age, sex, severity, and disease stages of COVID-19. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) RNA was found in diverse epithelial and immune cell types, accompanied by dramatic transcriptomic changes within virus-positive cells. Systemic upregulation of S100A8/A9, mainly by megakaryocytes and monocytes in the peripheral blood, may contribute to the cytokine storms frequently observed in severe patients. Our data provide a rich resource for understanding the pathogenesis of and developing effective therapeutic strategies for COVID-19.

A role of PIEZO1 in iron metabolism in mice and humans.

Iron overload causes progressive organ damage and is associated with arthritis, liver damage, and heart failure. Elevated iron levels are present in 1%-5% of individuals; however, iron overload is undermonitored and underdiagnosed. Genetic factors affecting iron homeostasis are emerging. Individuals with hereditary xerocytosis, a rare disorder with gain-of-function (GOF) mutations in mechanosensitive PIEZO1 ion channel, develop age-onset iron overload. We show that constitutive or macrophage expression of a GOF Piezo1 allele in mice disrupts levels of the iron regulator hepcidin and causes iron overload. We further show that PIEZO1 is a key regulator of macrophage phagocytic activity and subsequent erythrocyte turnover. Strikingly, we find that E756del, a mild GOF PIEZO1 allele present in one-third of individuals of African descent, is strongly associated with increased plasma iron. Our study links macrophage mechanotransduction to iron metabolism and identifies a genetic risk factor for increased iron levels in African Americans.

Meta-analysis of tumor- and T cell-intrinsic mechanisms of sensitization to checkpoint inhibition.

Checkpoint inhibitors (CPIs) augment adaptive immunity. Systematic pan-tumor analyses may reveal the relative importance of tumor-cell-intrinsic and microenvironmental features underpinning CPI sensitization. Here, we collated whole-exome and transcriptomic data for >1,000 CPI-treated patients across seven tumor types, utilizing standardized bioinformatics workflows and clinical outcome criteria to validate multivariable predictors of CPI sensitization. Clonal tumor mutation burden (TMB) was the strongest predictor of CPI response, followed by total TMB and CXCL9 expression. Subclonal TMB, somatic copy alteration burden, and histocompatibility leukocyte antigen (HLA) evolutionary divergence failed to attain pan-cancer significance. Dinucleotide variants were identified as a source of immunogenic epitopes associated with radical amino acid substitutions and enhanced peptide hydrophobicity/immunogenicity. Copy-number analysis revealed two additional determinants of CPI outcome supported by prior functional evidence: 9q34 (TRAF2) loss associated with response and CCND1 amplification associated with resistance. Finally, single-cell RNA sequencing (RNA-seq) of clonal neoantigen-reactive CD8 tumor-infiltrating lymphocytes (TILs), combined with bulk RNA-seq analysis of CPI-responding tumors, identified CCR5 and CXCL13 as T-cell-intrinsic markers of CPI sensitivity.

Translocated microbiome composition determines immunological outcome in treated HIV infection.

The impact of the microbiome on HIV disease is widely acknowledged although the mechanisms downstream of fluctuations in microbial composition remain speculative. We detected rapid, dynamic changes in translocated microbial constituents during two years after cART initiation. An unbiased systems biology approach revealed two distinct pathways driven by changes in the abundance ratio of Serratia to other bacterial genera. Increased CD4 T cell numbers over the first year were associated with high Serratia abundance, pro-inflammatory innate cytokines, and metabolites that drive Th17 gene expression signatures and restoration of mucosal integrity. Subsequently, decreased Serratia abundance and downregulation of innate cytokines allowed re-establishment of systemic T cell homeostasis promoting restoration of Th1 and Th2 gene expression signatures. Analyses of three other geographically distinct cohorts of treated HIV infection established a more generalized principle that changes in diversity and composition of translocated microbial species influence systemic inflammation and consequently CD4 T cell recovery.

Proteogenomic characterization of pancreatic ductal adenocarcinoma.

Pancreatic ductal adenocarcinoma (PDAC) is a highly aggressive cancer with poor patient survival. Toward understanding the underlying molecular alterations that drive PDAC oncogenesis, we conducted comprehensive proteogenomic analysis of 140 pancreatic cancers, 67 normal adjacent tissues, and 9 normal pancreatic ductal tissues. Proteomic, phosphoproteomic, and glycoproteomic analyses were used to characterize proteins and their modifications. In addition, whole-genome sequencing, whole-exome sequencing, methylation, RNA sequencing (RNA-seq), and microRNA sequencing (miRNA-seq) were performed on the same tissues to facilitate an integrated proteogenomic analysis and determine the impact of genomic alterations on protein expression, signaling pathways, and post-translational modifications. To ensure robust downstream analyses, tumor neoplastic cellularity was assessed via multiple orthogonal strategies using molecular features and verified via pathological estimation of tumor cellularity based on histological review. This integrated proteogenomic characterization of PDAC will serve as a valuable resource for the community, paving the way for early detection and identification of novel therapeutic targets.

Impaired local intrinsic immunity to SARS-CoV-2 infection in severe COVID-19.

SARS-CoV-2 infection can cause severe respiratory COVID-19. However, many individuals present with isolated upper respiratory symptoms, suggesting potential to constrain viral pathology to the nasopharynx. Which cells SARS-CoV-2 primarily targets and how infection influences the respiratory epithelium remains incompletely understood. We performed scRNA-seq on nasopharyngeal swabs from 58 healthy and COVID-19 participants. During COVID-19, we observe expansion of secretory, loss of ciliated, and epithelial cell repopulation via deuterosomal cell expansion. In mild and moderate COVID-19, epithelial cells express anti-viral/interferon-responsive genes, while cells in severe COVID-19 have muted anti-viral responses despite equivalent viral loads. SARS-CoV-2 RNA+ host-target cells are highly heterogenous, including developing ciliated, interferon-responsive ciliated, AZGP1high goblet, and KRT13+ "hillock"-like cells, and we identify genes associated with susceptibility, resistance, or infection response. Our study defines protective and detrimental responses to SARS-CoV-2, the direct viral targets of infection, and suggests that failed nasal epithelial anti-viral immunity may underlie and precede severe COVID-19.

Spatially organized multicellular immune hubs in human colorectal cancer.

Immune responses to cancer are highly variable, with mismatch repair-deficient (MMRd) tumors exhibiting more anti-tumor immunity than mismatch repair-proficient (MMRp) tumors. To understand the rules governing these varied responses, we transcriptionally profiled 371,223 cells from colorectal tumors and adjacent normal tissues of 28 MMRp and 34 MMRd individuals. Analysis of 88 cell subsets and their 204 associated gene expression programs revealed extensive transcriptional and spatial remodeling across tumors. To discover hubs of interacting malignant and immune cells, we identified expression programs in different cell types that co-varied across tumors from affected individuals and used spatial profiling to localize coordinated programs. We discovered a myeloid cell-attracting hub at the tumor-luminal interface associated with tissue damage and an MMRd-enriched immune hub within the tumor, with activated T cells together with malignant and myeloid cells expressing T cell-attracting chemokines. By identifying interacting cellular programs, we reveal the logic underlying spatially organized immune-malignant cell networks.

Single-cell protein activity analysis identifies recurrence-associated renal tumor macrophages.

Clear cell renal carcinoma (ccRCC) is a heterogeneous disease with a variable post-surgical course. To assemble a comprehensive ccRCC tumor microenvironment (TME) atlas, we performed single-cell RNA sequencing (scRNA-seq) of hematopoietic and non-hematopoietic subpopulations from tumor and tumor-adjacent tissue of treatment-naive ccRCC resections. We leveraged the VIPER algorithm to quantitate single-cell protein activity and validated this approach by comparison to flow cytometry. The analysis identified key TME subpopulations, as well as their master regulators and candidate cell-cell interactions, revealing clinically relevant populations, undetectable by gene-expression analysis. Specifically, we uncovered a tumor-specific macrophage subpopulation characterized by upregulation of TREM2/APOE/C1Q, validated by spatially resolved, quantitative multispectral immunofluorescence. In a large clinical validation cohort, these markers were significantly enriched in tumors from patients who recurred following surgery. The study thus identifies TREM2/APOE/C1Q-positive macrophage infiltration as a potential prognostic biomarker for ccRCC recurrence, as well as a candidate therapeutic target.

NNT mediates redox-dependent pigmentation via a UVB- and MITF-independent mechanism.

Ultraviolet (UV) light and incompletely understood genetic and epigenetic variations determine skin color. Here we describe an UV- and microphthalmia-associated transcription factor (MITF)-independent mechanism of skin pigmentation. Targeting the mitochondrial redox-regulating enzyme nicotinamide nucleotide transhydrogenase (NNT) resulted in cellular redox changes that affect tyrosinase degradation. These changes regulate melanosome maturation and, consequently, eumelanin levels and pigmentation. Topical application of small-molecule inhibitors yielded skin darkening in human skin, and mice with decreased NNT function displayed increased pigmentation. Additionally, genetic modification of NNT in zebrafish alters melanocytic pigmentation. Analysis of four diverse human cohorts revealed significant associations of skin color, tanning, and sun protection use with various single-nucleotide polymorphisms within NNT. NNT levels were independent of UVB irradiation and redox modulation. Individuals with postinflammatory hyperpigmentation or lentigines displayed decreased skin NNT levels, suggesting an NNT-driven, redox-dependent pigmentation mechanism that can be targeted with NNT-modifying topical drugs for medical and cosmetic purposes.

SARS-CoV-2 infection triggers profibrotic macrophage responses and lung fibrosis.

COVID-19-induced "acute respiratory distress syndrome" (ARDS) is associated with prolonged respiratory failure and high mortality, but the mechanistic basis of lung injury remains incompletely understood. Here, we analyze pulmonary immune responses and lung pathology in two cohorts of patients with COVID-19 ARDS using functional single-cell genomics, immunohistology, and electron microscopy. We describe an accumulation of CD163-expressing monocyte-derived macrophages that acquired a profibrotic transcriptional phenotype during COVID-19 ARDS. Gene set enrichment and computational data integration revealed a significant similarity between COVID-19-associated macrophages and profibrotic macrophage populations identified in idiopathic pulmonary fibrosis. COVID-19 ARDS was associated with clinical, radiographic, histopathological, and ultrastructural hallmarks of pulmonary fibrosis. Exposure of human monocytes to SARS-CoV-2, but not influenza A virus or viral RNA analogs, was sufficient to induce a similar profibrotic phenotype in vitro. In conclusion, we demonstrate that SARS-CoV-2 triggers profibrotic macrophage responses and pronounced fibroproliferative ARDS.

The Magnifying GLASS: Longitudinal Analysis of Adult Diffuse Gliomas.

Diffuse gliomas inevitably progress, but our understanding of the molecular events associated with recurrence is limited. Recent work from the Glioma Longitudinal Analysis (GLASS) consortium (Barthel et al., 2019) reports temporal DNA sequencing on a large cohort of primary and recurrent glioma pairs, establishing the evolutionary molecular characteristics of adult diffuse gliomas.

Key Parameters of Tumor Epitope Immunogenicity Revealed Through a Consortium Approach Improve Neoantigen Prediction.

Many approaches to identify therapeutically relevant neoantigens couple tumor sequencing with bioinformatic algorithms and inferred rules of tumor epitope immunogenicity. However, there are no reference data to compare these approaches, and the parameters governing tumor epitope immunogenicity remain unclear. Here, we assembled a global consortium wherein each participant predicted immunogenic epitopes from shared tumor sequencing data. 608 epitopes were subsequently assessed for T cell binding in patient-matched samples. By integrating peptide features associated with presentation and recognition, we developed a model of tumor epitope immunogenicity that filtered out 98% of non-immunogenic peptides with a precision above 0.70. Pipelines prioritizing model features had superior performance, and pipeline alterations leveraging them improved prediction performance. These findings were validated in an independent cohort of 310 epitopes prioritized from tumor sequencing data and assessed for T cell binding. This data resource enables identification of parameters underlying effective anti-tumor immunity and is available to the research community.

Identification of ALK in Thinness.

There is considerable inter-individual variability in susceptibility to weight gain despite an equally obesogenic environment in large parts of the world. Whereas many studies have focused on identifying the genetic susceptibility to obesity, we performed a GWAS on metabolically healthy thin individuals (lowest 6th percentile of the population-wide BMI spectrum) in a uniquely phenotyped Estonian cohort. We discovered anaplastic lymphoma kinase (ALK) as a candidate thinness gene. In Drosophila, RNAi mediated knockdown of Alk led to decreased triglyceride levels. In mice, genetic deletion of Alk resulted in thin animals with marked resistance to diet- and leptin-mutation-induced obesity. Mechanistically, we found that ALK expression in hypothalamic neurons controls energy expenditure via sympathetic control of adipose tissue lipolysis. Our genetic and mechanistic experiments identify ALK as a thinness gene, which is involved in the resistance to weight gain.

Tau PTM Profiles Identify Patient Heterogeneity and Stages of Alzheimer's Disease.

To elucidate the role of Tau isoforms and post-translational modification (PTM) stoichiometry in Alzheimer's disease (AD), we generated a high-resolution quantitative proteomics map of 95 PTMs on multiple isoforms of Tau isolated from postmortem human tissue from 49 AD and 42 control subjects. Although Tau PTM maps reveal heterogeneity across subjects, a subset of PTMs display high occupancy and frequency for AD, suggesting importance in disease. Unsupervised analyses indicate that PTMs occur in an ordered manner, leading to Tau aggregation. The processive addition and minimal set of PTMs associated with seeding activity was further defined by analysis of size-fractionated Tau. To summarize, features in the Tau protein critical for disease intervention at different stages of disease are identified, including enrichment of 0N and 4R isoforms, underrepresentation of the C terminus, an increase in negative charge in the proline-rich region (PRR), and a decrease in positive charge in the microtubule binding domain (MBD).