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.

Intrinsically disordered domain of transcription factor TCF-1 is required for T cell developmental fidelity.

In development, pioneer transcription factors access silent chromatin to reveal lineage-specific gene programs. The structured DNA-binding domains of pioneer factors have been well characterized, but whether and how intrinsically disordered regions affect chromatin and control cell fate is unclear. Here, we report that deletion of an intrinsically disordered region of the pioneer factor TCF-1 (termed L1) leads to an early developmental block in T cells. The few T cells that develop from progenitors expressing TCF-1 lacking L1 exhibit lineage infidelity distinct from the lineage diversion of TCF-1-deficient cells. Mechanistically, L1 is required for activation of T cell genes and repression of GATA2-driven genes, normally reserved to the mast cell and dendritic cell lineages. Underlying this lineage diversion, L1 mediates binding of TCF-1 to its earliest target genes, which are subject to repression as T cells develop. These data suggest that the intrinsically disordered N terminus of TCF-1 maintains T cell lineage fidelity.

The Human Tumor Atlas Network: Charting Tumor Transitions across Space and Time at Single-Cell Resolution.

Crucial transitions in cancer-including tumor initiation, local expansion, metastasis, and therapeutic resistance-involve complex interactions between cells within the dynamic tumor ecosystem. Transformative single-cell genomics technologies and spatial multiplex in situ methods now provide an opportunity to interrogate this complexity at unprecedented resolution. The Human Tumor Atlas Network (HTAN), part of the National Cancer Institute (NCI) Cancer Moonshot Initiative, will establish a clinical, experimental, computational, and organizational framework to generate informative and accessible three-dimensional atlases of cancer transitions for a diverse set of tumor types. This effort complements both ongoing efforts to map healthy organs and previous large-scale cancer genomics approaches focused on bulk sequencing at a single point in time. Generating single-cell, multiparametric, longitudinal atlases and integrating them with clinical outcomes should help identify novel predictive biomarkers and features as well as therapeutically relevant cell types, cell states, and cellular interactions across transitions. The resulting tumor atlases should have a profound impact on our understanding of cancer biology and have the potential to improve cancer detection, prevention, and therapeutic discovery for better precision-medicine treatments of cancer patients and those at risk for cancer.

Distinct Classes of Complex Structural Variation Uncovered across Thousands of Cancer Genome Graphs.

Cancer genomes often harbor hundreds of somatic DNA rearrangement junctions, many of which cannot be easily classified into simple (e.g., deletion) or complex (e.g., chromothripsis) structural variant classes. Applying a novel genome graph computational paradigm to analyze the topology of junction copy number (JCN) across 2,778 tumor whole-genome sequences, we uncovered three novel complex rearrangement phenomena: pyrgo, rigma, and tyfonas. Pyrgo are "towers" of low-JCN duplications associated with early-replicating regions, superenhancers, and breast or ovarian cancers. Rigma comprise "chasms" of low-JCN deletions enriched in late-replicating fragile sites and gastrointestinal carcinomas. Tyfonas are "typhoons" of high-JCN junctions and fold-back inversions associated with expressed protein-coding fusions, breakend hypermutation, and acral, but not cutaneous, melanomas. Clustering of tumors according to genome graph-derived features identified subgroups associated with DNA repair defects and poor prognosis.

Illuminating G-Protein-Coupling Selectivity of GPCRs.

Heterotrimetic G proteins consist of four subfamilies (Gs, Gi/o, Gq/11, and G12/13) that mediate signaling via G-protein-coupled receptors (GPCRs), principally by receptors binding Gα C termini. G-protein-coupling profiles govern GPCR-induced cellular responses, yet receptor sequence selectivity determinants remain elusive. Here, we systematically quantified ligand-induced interactions between 148 GPCRs and all 11 unique Gα subunit C termini. For each receptor, we probed chimeric Gα subunit activation via a transforming growth factor-α (TGF-α) shedding response in HEK293 cells lacking endogenous Gq/11 and G12/13 proteins, and complemented G-protein-coupling profiles through a NanoBiT-G-protein dissociation assay. Interrogation of the dataset identified sequence-based coupling specificity features, inside and outside the transmembrane domain, which we used to develop a coupling predictor that outperforms previous methods. We used the predictor to engineer designer GPCRs selectively coupled to G12. This dataset of fine-tuned signaling mechanisms for diverse GPCRs is a valuable resource for research in GPCR signaling.

Excess Translation of Epigenetic Regulators Contributes to Fragile X Syndrome and Is Alleviated by Brd4 Inhibition.

Fragile X syndrome (FXS) is a leading genetic cause of intellectual disability and autism. FXS results from the loss of function of fragile X mental retardation protein (FMRP), which represses translation of target transcripts. Most of the well-characterized target transcripts of FMRP are synaptic proteins, yet targeting these proteins has not provided effective treatments. We examined a group of FMRP targets that encode transcriptional regulators, particularly chromatin-associated proteins. Loss of FMRP in mice results in widespread changes in chromatin regulation and aberrant gene expression. To determine if targeting epigenetic factors could reverse phenotypes associated with the disorder, we focused on Brd4, a BET protein and chromatin reader targeted by FMRP. Inhibition of Brd4 function alleviated many of the phenotypes associated with FXS. We conclude that loss of FMRP results in significant epigenetic misregulation and that targeting transcription via epigenetic regulators like Brd4 may provide new treatments for FXS.

Genomic Patterns of De Novo Mutation in Simplex Autism.

To further our understanding of the genetic etiology of autism, we generated and analyzed genome sequence data from 516 idiopathic autism families (2,064 individuals). This resource includes >59 million single-nucleotide variants (SNVs) and 9,212 private copy number variants (CNVs), of which 133,992 and 88 are de novo mutations (DNMs), respectively. We estimate a mutation rate of ∼1.5 × 10-8 SNVs per site per generation with a significantly higher mutation rate in repetitive DNA. Comparing probands and unaffected siblings, we observe several DNM trends. Probands carry more gene-disruptive CNVs and SNVs, resulting in severe missense mutations and mapping to predicted fetal brain promoters and embryonic stem cell enhancers. These differences become more pronounced for autism genes (p = 1.8 × 10-3, OR = 2.2). Patients are more likely to carry multiple coding and noncoding DNMs in different genes, which are enriched for expression in striatal neurons (p = 3 × 10-3), suggesting a path forward for genetically characterizing more complex cases of autism.

Core Clinical Data Elements for Cancer Genomic Repositories: A Multi-stakeholder Consensus.

The Center for Medical Technology Policy and the Molecular Evidence Development Consortium gathered a diverse group of more than 50 stakeholders to develop consensus on a core set of data elements and values essential to understanding the clinical utility of molecularly targeted therapies in oncology.

Crystal Structure of the Human Cannabinoid Receptor CB1.

Cannabinoid receptor 1 (CB1) is the principal target of Δ9-tetrahydrocannabinol (THC), a psychoactive chemical from Cannabis sativa with a wide range of therapeutic applications and a long history of recreational use. CB1 is activated by endocannabinoids and is a promising therapeutic target for pain management, inflammation, obesity, and substance abuse disorders. Here, we present the 2.8 Å crystal structure of human CB1 in complex with AM6538, a stabilizing antagonist, synthesized and characterized for this structural study. The structure of the CB1-AM6538 complex reveals key features of the receptor and critical interactions for antagonist binding. In combination with functional studies and molecular modeling, the structure provides insight into the binding mode of naturally occurring CB1 ligands, such as THC, and synthetic cannabinoids. This enhances our understanding of the molecular basis for the physiological functions of CB1 and provides new opportunities for the design of next-generation CB1-targeting pharmaceuticals.

Fibrin-targeting immunotherapy protects against neuroinflammation and neurodegeneration.

Activation of innate immunity and deposition of blood-derived fibrin in the central nervous system (CNS) occur in autoimmune and neurodegenerative diseases, including multiple sclerosis (MS) and Alzheimer's disease (AD). However, the mechanisms that link disruption of the blood-brain barrier (BBB) to neurodegeneration are poorly understood, and exploration of fibrin as a therapeutic target has been limited by its beneficial clotting functions. Here we report the generation of monoclonal antibody 5B8, targeted against the cryptic fibrin epitope γ377-395, to selectively inhibit fibrin-induced inflammation and oxidative stress without interfering with clotting. 5B8 suppressed fibrin-induced nicotinamide adenine dinucleotide phosphate (NADPH) oxidase activation and the expression of proinflammatory genes. In animal models of MS and AD, 5B8 entered the CNS and bound to parenchymal fibrin, and its therapeutic administration reduced the activation of innate immunity and neurodegeneration. Thus, fibrin-targeting immunotherapy inhibited autoimmunity- and amyloid-driven neurotoxicity and might have clinical benefit without globally suppressing innate immunity or interfering with coagulation in diverse neurological diseases.

Hepatitis C virus RNA functionally sequesters miR-122.

Hepatitis C virus (HCV) uniquely requires the liver-specific microRNA-122 for replication, yet global effects on endogenous miRNA targets during infection are unexplored. Here, high-throughput sequencing and crosslinking immunoprecipitation (HITS-CLIP) experiments of human Argonaute (AGO) during HCV infection showed robust AGO binding on the HCV 5'UTR at known and predicted miR-122 sites. On the human transcriptome, we observed reduced AGO binding and functional mRNA de-repression of miR-122 targets during virus infection. This miR-122 "sponge" effect was relieved and redirected to miR-15 targets by swapping the miRNA tropism of the virus. Single-cell expression data from reporters containing miR-122 sites showed significant de-repression during HCV infection depending on expression level and site number. We describe a quantitative mathematical model of HCV-induced miR-122 sequestration and propose that such miR-122 inhibition by HCV RNA may result in global de-repression of host miR-122 targets, providing an environment fertile for the long-term oncogenic potential of HCV.

3D genomic mapping reveals multifocality of human pancreatic precancers.

Pancreatic intraepithelial neoplasias (PanINs) are the most common precursors of pancreatic cancer, but their small size and inaccessibility in humans make them challenging to study1. Critically, the number, dimensions and connectivity of human PanINs remain largely unknown, precluding important insights into early cancer development. Here, we provide a microanatomical survey of human PanINs by analysing 46 large samples of grossly normal human pancreas with a machine-learning pipeline for quantitative 3D histological reconstruction at single-cell resolution. To elucidate genetic relationships between and within PanINs, we developed a workflow in which 3D modelling guides multi-region microdissection and targeted and whole-exome sequencing. From these samples, we calculated a mean burden of 13 PanINs per cm3 and extrapolated that the normal intact adult pancreas harbours hundreds of PanINs, almost all with oncogenic KRAS hotspot mutations. We found that most PanINs originate as independent clones with distinct somatic mutation profiles. Some spatially continuous PanINs were found to contain multiple KRAS mutations; computational and in situ analyses demonstrated that different KRAS mutations localize to distinct cell subpopulations within these neoplasms, indicating their polyclonal origins. The extensive multifocality and genetic heterogeneity of PanINs raises important questions about mechanisms that drive precancer initiation and confer differential progression risk in the human pancreas. This detailed 3D genomic mapping of molecular alterations in human PanINs provides an empirical foundation for early detection and rational interception of pancreatic cancer.

EBF1 nuclear repositioning instructs chromatin refolding to promote therapy resistance in T leukemic cells.

Chromatin misfolding has been implicated in cancer pathogenesis; yet, its role in therapy resistance remains unclear. Here, we systematically integrated sequencing and imaging data to examine the spatial and linear chromatin structures in targeted therapy-sensitive and -resistant human T cell acute lymphoblastic leukemia (T-ALL). We found widespread alterations in successive layers of chromatin organization including spatial compartments, contact domain boundaries, and enhancer positioning upon the emergence of targeted therapy resistance. The reorganization of genome folding structures closely coincides with the restructuring of chromatin activity and redistribution of architectural proteins. Mechanistically, the derepression and repositioning of the B-lineage-determining transcription factor EBF1 from the heterochromatic nuclear envelope to the euchromatic interior instructs widespread genome refolding and promotes therapy resistance in leukemic T cells. Together, our findings suggest that lineage-determining transcription factors can instruct changes in genome topology as a driving force for epigenetic adaptations in targeted therapy resistance.

YTHDC2 control of gametogenesis requires helicase activity but not m6A binding.

Mechanisms regulating meiotic progression in mammals are poorly understood. The N6-methyladenosine (m6A) reader and 3' → 5' RNA helicase YTHDC2 switches cells from mitotic to meiotic gene expression programs and is essential for meiotic entry, but how this critical cell fate change is accomplished is unknown. Here, we provide insight into its mechanism and implicate YTHDC2 in having a broad role in gene regulation during multiple meiotic stages. Unexpectedly, mutation of the m6A-binding pocket of YTHDC2 had no detectable effect on gametogenesis and mouse fertility, suggesting that YTHDC2 function is m6A-independent. Supporting this conclusion, CLIP data defined YTHDC2-binding sites on mRNA as U-rich and UG-rich motif-containing regions within 3' UTRs and coding sequences, distinct from the sites that contain m6A during spermatogenesis. Complete loss of YTHDC2 during meiotic entry did not substantially alter translation of its mRNA binding targets in whole-testis ribosome profiling assays but did modestly affect their steady-state levels. Mutation of the ATPase motif in the helicase domain of YTHDC2 did not affect meiotic entry, but it blocked meiotic prophase I progression, causing sterility. Our findings inform a model in which YTHDC2 binds transcripts independent of m6A status and regulates gene expression during multiple stages of meiosis by distinct mechanisms.

Genetic Variation in Type 1 Diabetes Reconfigures the 3D Chromatin Organization of T Cells and Alters Gene Expression.

Genetics is a major determinant of susceptibility to autoimmune disorders. Here, we examined whether genome organization provides resilience or susceptibility to sequence variations, and how this would contribute to the molecular etiology of an autoimmune disease. We generated high-resolution maps of linear and 3D genome organization in thymocytes of NOD mice, a model of type 1 diabetes (T1D), and the diabetes-resistant C57BL/6 mice. Multi-enhancer interactions formed at genomic regions harboring genes with prominent roles in T cell development in both strains. However, diabetes risk-conferring loci coalesced enhancers and promoters in NOD, but not C57BL/6 thymocytes. 3D genome mapping of NODxC57BL/6 F1 thymocytes revealed that genomic misfolding in NOD mice is mediated in cis. Moreover, immune cells infiltrating the pancreas of humans with T1D exhibited increased expression of genes located on misfolded loci in mice. Thus, genetic variation leads to altered 3D chromatin architecture and associated changes in gene expression that may underlie autoimmune pathology.

Extreme dynamics in a biomolecular condensate.

Proteins and nucleic acids can phase-separate in the cell to form concentrated biomolecular condensates1-4. The functions of condensates span many length scales: they modulate interactions and chemical reactions at the molecular scale5, organize biochemical processes at the mesoscale6 and compartmentalize cells4. Understanding the underlying mechanisms of these processes will require detailed knowledge of the rich dynamics across these scales7. The mesoscopic dynamics of biomolecular condensates have been extensively characterized8, but their behaviour at the molecular scale has remained more elusive. Here, as an example of biomolecular phase separation, we study complex coacervates of two highly and oppositely charged disordered human proteins9. Their dense phase is 1,000 times more concentrated than the dilute phase, and the resulting percolated interaction network10 leads to a bulk viscosity 300 times greater than that of water. However, single-molecule spectroscopy optimized for measurements within individual droplets reveals that at the molecular scale, the disordered proteins remain exceedingly dynamic, with their chain configurations interconverting on submicrosecond timescales. Massive all-atom molecular dynamics simulations reproduce the experimental observations and explain this apparent discrepancy: the underlying interactions between individual charged side chains are short-lived and exchange on a pico- to nanosecond timescale. Our results indicate that, despite the high macroscopic viscosity of phase-separated systems, local biomolecular rearrangements required for efficient reactions at the molecular scale can remain rapid.

Extensive global wetland loss over the past three centuries.

Wetlands have long been drained for human use, thereby strongly affecting greenhouse gas fluxes, flood control, nutrient cycling and biodiversity1,2. Nevertheless, the global extent of natural wetland loss remains remarkably uncertain3. Here, we reconstruct the spatial distribution and timing of wetland loss through conversion to seven human land uses between 1700 and 2020, by combining national and subnational records of drainage and conversion with land-use maps and simulated wetland extents. We estimate that 3.4 million km2 (confidence interval 2.9-3.8) of inland wetlands have been lost since 1700, primarily for conversion to croplands. This net loss of 21% (confidence interval 16-23%) of global wetland area is lower than that suggested previously by extrapolations of data disproportionately from high-loss regions. Wetland loss has been concentrated in Europe, the United States and China, and rapidly expanded during the mid-twentieth century. Our reconstruction elucidates the timing and land-use drivers of global wetland losses, providing an improved historical baseline to guide assessment of wetland loss impact on Earth system processes, conservation planning to protect remaining wetlands and prioritization of sites for wetland restoration4.

Normative spatiotemporal fetal brain maturation with satisfactory development at 2 years.

Maturation of the human fetal brain should follow precisely scheduled structural growth and folding of the cerebral cortex for optimal postnatal function1. We present a normative digital atlas of fetal brain maturation based on a prospective international cohort of healthy pregnant women2, selected using World Health Organization recommendations for growth standards3. Their fetuses were accurately dated in the first trimester, with satisfactory growth and neurodevelopment from early pregnancy to 2 years of age4,5. The atlas was produced using 1,059 optimal quality, three-dimensional ultrasound brain volumes from 899 of the fetuses and an automated analysis pipeline6-8. The atlas corresponds structurally to published magnetic resonance images9, but with finer anatomical details in deep grey matter. The between-study site variability represented less than 8.0% of the total variance of all brain measures, supporting pooling data from the eight study sites to produce patterns of normative maturation. We have thereby generated an average representation of each cerebral hemisphere between 14 and 31 weeks' gestation with quantification of intracranial volume variability and growth patterns. Emergent asymmetries were detectable from as early as 14 weeks, with peak asymmetries in regions associated with language development and functional lateralization between 20 and 26 weeks' gestation. These patterns were validated in 1,487 three-dimensional brain volumes from 1,295 different fetuses in the same cohort. We provide a unique spatiotemporal benchmark of fetal brain maturation from a large cohort with normative postnatal growth and neurodevelopment.

Design and testing of a humanized porcine donor for xenotransplantation.

Recent human decedent model studies1,2 and compassionate xenograft use3 have explored the promise of porcine organs for human transplantation. To proceed to human studies, a clinically ready porcine donor must be engineered and its xenograft successfully tested in nonhuman primates. Here we describe the design, creation and long-term life-supporting function of kidney grafts from a genetically engineered porcine donor transplanted into a cynomolgus monkey model. The porcine donor was engineered to carry 69 genomic edits, eliminating glycan antigens, overexpressing human transgenes and inactivating porcine endogenous retroviruses. In vitro functional analyses showed that the edited kidney endothelial cells modulated inflammation to an extent that was indistinguishable from that of human endothelial cells, suggesting that these edited cells acquired a high level of human immune compatibility. When transplanted into cynomolgus monkeys, the kidneys with three glycan antigen knockouts alone experienced poor graft survival, whereas those with glycan antigen knockouts and human transgene expression demonstrated significantly longer survival time, suggesting the benefit of human transgene expression in vivo. These results show that preclinical studies of renal xenotransplantation could be successfully conducted in nonhuman primates and bring us closer to clinical trials of genetically engineered porcine renal grafts.

ZMYND8-regulated IRF8 transcription axis is an acute myeloid leukemia dependency.

The transformed state in acute leukemia requires gene regulatory programs involving transcription factors and chromatin modulators. Here, we uncover an IRF8-MEF2D transcriptional circuit as an acute myeloid leukemia (AML)-biased dependency. We discover and characterize the mechanism by which the chromatin "reader" ZMYND8 directly activates IRF8 in parallel with the MYC proto-oncogene through their lineage-specific enhancers. ZMYND8 is essential for AML proliferation in vitro and in vivo and associates with MYC and IRF8 enhancer elements that we define in cell lines and in patient samples. ZMYND8 occupancy at IRF8 and MYC enhancers requires BRD4, a transcription coactivator also necessary for AML proliferation. We show that ZMYND8 binds to the ET domain of BRD4 via its chromatin reader cassette, which in turn is required for proper chromatin occupancy and maintenance of leukemic growth in vivo. Our results rationalize ZMYND8 as a potential therapeutic target for modulating essential transcriptional programs in AML.