The gene dose makes the disease.

A long-standing challenge in genomics has been to identify causal genes within rare copy-number variant regions that are intolerant to altered dosage. In this issue, Collins et al. perform a meta-analysis of almost a million individuals to identify dosage-sensitive segments and genes conferring risk for a range of disease phenotypes.

A cross-disorder dosage sensitivity map of the human genome.

Rare copy-number variants (rCNVs) include deletions and duplications that occur infrequently in the global human population and can confer substantial risk for disease. In this study, we aimed to quantify the properties of haploinsufficiency (i.e., deletion intolerance) and triplosensitivity (i.e., duplication intolerance) throughout the human genome. We harmonized and meta-analyzed rCNVs from nearly one million individuals to construct a genome-wide catalog of dosage sensitivity across 54 disorders, which defined 163 dosage sensitive segments associated with at least one disorder. These segments were typically gene dense and often harbored dominant dosage sensitive driver genes, which we were able to prioritize using statistical fine-mapping. Finally, we designed an ensemble machine-learning model to predict probabilities of dosage sensitivity (pHaplo & pTriplo) for all autosomal genes, which identified 2,987 haploinsufficient and 1,559 triplosensitive genes, including 648 that were uniquely triplosensitive. This dosage sensitivity resource will provide broad utility for human disease research and clinical genetics.

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.

Regulation of ribosomal RNA gene copy number, transcription and nucleolus organization in eukaryotes.

One of the first biological machineries to be created seems to have been the ribosome. Since then, organisms have dedicated great efforts to optimize this apparatus. The ribosomal RNA (rRNA) contained within ribosomes is crucial for protein synthesis and maintenance of cellular function in all known organisms. In eukaryotic cells, rRNA is produced from ribosomal DNA clusters of tandem rRNA genes, whose organization in the nucleolus, maintenance and transcription are strictly regulated to satisfy the substantial demand for rRNA required for ribosome biogenesis. Recent studies have elucidated mechanisms underlying the integrity of ribosomal DNA and regulation of its transcription, including epigenetic mechanisms and a unique recombination and copy-number control system to stably maintain high rRNA gene copy number. In this Review, we disucss how the crucial maintenance of rRNA gene copy number through control of gene amplification and of rRNA production by RNA polymerase I are orchestrated. We also discuss how liquid-liquid phase separation controls the architecture and function of the nucleolus and the relationship between rRNA production, cell senescence and disease.

Population Structure, Stratification, and Introgression of Human Structural Variation.

Structural variants contribute substantially to genetic diversity and are important evolutionarily and medically, but they are still understudied. Here we present a comprehensive analysis of structural variation in the Human Genome Diversity panel, a high-coverage dataset of 911 samples from 54 diverse worldwide populations. We identify, in total, 126,018 variants, 78% of which were not identified in previous global sequencing projects. Some reach high frequency and are private to continental groups or even individual populations, including regionally restricted runaway duplications and putatively introgressed variants from archaic hominins. By de novo assembly of 25 genomes using linked-read sequencing, we discover 1,643 breakpoint-resolved unique insertions, in aggregate accounting for 1.9 Mb of sequence absent from the GRCh38 reference. Our results illustrate the limitation of a single human reference and the need for high-quality genomes from diverse populations to fully discover and understand human genetic variation.

Smc3 dosage regulates B cell transit through germinal centers and restricts their malignant transformation.

During the germinal center (GC) reaction, B cells undergo extensive redistribution of cohesin complex and three-dimensional reorganization of their genomes. Yet, the significance of cohesin and architectural programming in the humoral immune response is unknown. Herein we report that homozygous deletion of Smc3, encoding the cohesin ATPase subunit, abrogated GC formation, while, in marked contrast, Smc3 haploinsufficiency resulted in GC hyperplasia, skewing of GC polarity and impaired plasma cell (PC) differentiation. Genome-wide chromosomal conformation and transcriptional profiling revealed defects in GC B cell terminal differentiation programs controlled by the lymphoma epigenetic tumor suppressors Tet2 and Kmt2d and failure of Smc3-haploinsufficient GC B cells to switch from B cell- to PC-defining transcription factors. Smc3 haploinsufficiency preferentially impaired the connectivity of enhancer elements controlling various lymphoma tumor suppressor genes, and, accordingly, Smc3 haploinsufficiency accelerated lymphomagenesis in mice with constitutive Bcl6 expression. Collectively, our data indicate a dose-dependent function for cohesin in humoral immunity to facilitate the B cell to PC phenotypic switch while restricting malignant transformation.

Timing the Landmark Events in the Evolution of Clear Cell Renal Cell Cancer: TRACERx Renal.

Clear cell renal cell carcinoma (ccRCC) is characterized by near-universal loss of the short arm of chromosome 3, deleting several tumor suppressor genes. We analyzed whole genomes from 95 biopsies across 33 patients with clear cell renal cell carcinoma. We find hotspots of point mutations in the 5' UTR of TERT, targeting a MYC-MAX-MAD1 repressor associated with telomere lengthening. The most common structural abnormality generates simultaneous 3p loss and 5q gain (36% patients), typically through chromothripsis. This event occurs in childhood or adolescence, generally as the initiating event that precedes emergence of the tumor's most recent common ancestor by years to decades. Similar genomic changes drive inherited ccRCC. Modeling differences in age incidence between inherited and sporadic cancers suggests that the number of cells with 3p loss capable of initiating sporadic tumors is no more than a few hundred. Early development of ccRCC follows well-defined evolutionary trajectories, offering opportunity for early intervention.

Quantitative analysis of T cell proteomes and environmental sensors during T cell differentiation.

Quantitative mass spectrometry reveals how CD4+ and CD8+ T cells restructure proteomes in response to antigen and mammalian target of rapamycin complex 1 (mTORC1). Analysis of copy numbers per cell of >9,000 proteins provides new understanding of T cell phenotypes, exposing the metabolic and protein synthesis machinery and environmental sensors that shape T cell fate. We reveal that lymphocyte environment sensing is controlled by immune activation, and that CD4+ and CD8+ T cells differ in their intrinsic nutrient transport and biosynthetic capacity. Our data also reveal shared and divergent outcomes of mTORC1 inhibition in naïve versus effector T cells: mTORC1 inhibition impaired cell cycle progression in activated naïve cells, but not effector cells, whereas metabolism was consistently impacted in both populations. This study provides a comprehensive map of naïve and effector T cell proteomes, and a resource for exploring and understanding T cell phenotypes and cell context effects of mTORC1.

Integrated Proteogenomic Characterization of Human High-Grade Serous Ovarian Cancer.

To provide a detailed analysis of the molecular components and underlying mechanisms associated with ovarian cancer, we performed a comprehensive mass-spectrometry-based proteomic characterization of 174 ovarian tumors previously analyzed by The Cancer Genome Atlas (TCGA), of which 169 were high-grade serous carcinomas (HGSCs). Integrating our proteomic measurements with the genomic data yielded a number of insights into disease, such as how different copy-number alternations influence the proteome, the proteins associated with chromosomal instability, the sets of signaling pathways that diverse genome rearrangements converge on, and the ones most associated with short overall survival. Specific protein acetylations associated with homologous recombination deficiency suggest a potential means for stratifying patients for therapy. In addition to providing a valuable resource, these findings provide a view of how the somatic genome drives the cancer proteome and associations between protein and post-translational modification levels and clinical outcomes in HGSC. VIDEO ABSTRACT.

A Landscape of Pharmacogenomic Interactions in Cancer.

Systematic studies of cancer genomes have provided unprecedented insights into the molecular nature of cancer. Using this information to guide the development and application of therapies in the clinic is challenging. Here, we report how cancer-driven alterations identified in 11,289 tumors from 29 tissues (integrating somatic mutations, copy number alterations, DNA methylation, and gene expression) can be mapped onto 1,001 molecularly annotated human cancer cell lines and correlated with sensitivity to 265 drugs. We find that cell lines faithfully recapitulate oncogenic alterations identified in tumors, find that many of these associate with drug sensitivity/resistance, and highlight the importance of tissue lineage in mediating drug response. Logic-based modeling uncovers combinations of alterations that sensitize to drugs, while machine learning demonstrates the relative importance of different data types in predicting drug response. Our analysis and datasets are rich resources to link genotypes with cellular phenotypes and to identify therapeutic options for selected cancer sub-populations.

Mutations in Human Accelerated Regions Disrupt Cognition and Social Behavior.

Comparative analyses have identified genomic regions potentially involved in human evolution but do not directly assess function. Human accelerated regions (HARs) represent conserved genomic loci with elevated divergence in humans. If some HARs regulate human-specific social and behavioral traits, then mutations would likely impact cognitive and social disorders. Strikingly, rare biallelic point mutations-identified by whole-genome and targeted "HAR-ome" sequencing-showed a significant excess in individuals with ASD whose parents share common ancestry compared to familial controls, suggesting a contribution in 5% of consanguineous ASD cases. Using chromatin interaction sequencing, massively parallel reporter assays (MPRA), and transgenic mice, we identified disease-linked, biallelic HAR mutations in active enhancers for CUX1, PTBP2, GPC4, CDKL5, and other genes implicated in neural function, ASD, or both. Our data provide genetic evidence that specific HARs are essential for normal development, consistent with suggestions that their evolutionary changes may have altered social and/or cognitive behavior. PAPERCLIP.

Functional Genomic Landscape of Human Breast Cancer Drivers, Vulnerabilities, and Resistance.

Large-scale genomic studies have identified multiple somatic aberrations in breast cancer, including copy number alterations and point mutations. Still, identifying causal variants and emergent vulnerabilities that arise as a consequence of genetic alterations remain major challenges. We performed whole-genome small hairpin RNA (shRNA) "dropout screens" on 77 breast cancer cell lines. Using a hierarchical linear regression algorithm to score our screen results and integrate them with accompanying detailed genetic and proteomic information, we identify vulnerabilities in breast cancer, including candidate "drivers," and reveal general functional genomic properties of cancer cells. Comparisons of gene essentiality with drug sensitivity data suggest potential resistance mechanisms, effects of existing anti-cancer drugs, and opportunities for combination therapy. Finally, we demonstrate the utility of this large dataset by identifying BRD4 as a potential target in luminal breast cancer and PIK3CA mutations as a resistance determinant for BET-inhibitors.

Chromosomal Arrangement of Phosphorelay Genes Couples Sporulation and DNA Replication.

Genes encoding proteins in a common regulatory network are frequently located close to one another on the chromosome to facilitate co-regulation or couple gene expression to growth rate. Contrasting with these observations, here, we demonstrate a functional role for the arrangement of Bacillus subtilis sporulation network genes on opposite sides of the chromosome. We show that the arrangement of two sporulation network genes, one located close to the origin and the other close to the terminus, leads to a transient gene dosage imbalance during chromosome replication. This imbalance is detected by the sporulation network to produce cell-cycle coordinated pulses of the sporulation master regulator Spo0A∼P. This pulsed response allows cells to decide between sporulation and continued vegetative growth during each cell cycle spent in starvation. The simplicity of this coordination mechanism suggests that it may be widely applicable in a variety of gene regulatory and stress-response settings. VIDEO ABSTRACT.

Extensive strain-level copy-number variation across human gut microbiome species.

Within each bacterial species, different strains may vary in the set of genes they encode or in the copy number of these genes. Yet, taxonomic characterization of the human microbiota is often limited to the species level or to previously sequenced strains, and accordingly, the prevalence of intra-species variation, its functional role, and its relation to host health remain unclear. Here, we present a comprehensive large-scale analysis of intra-species copy-number variation in the gut microbiome, introducing a rigorous computational pipeline for detecting such variation directly from shotgun metagenomic data. We uncover a large set of variable genes in numerous species and demonstrate that this variation has significant functional and clinically relevant implications. We additionally infer intra-species compositional profiles, identifying population structure shifts and the presence of yet uncharacterized variants. Our results highlight the complex relationship between microbiome composition and functional capacity, linking metagenome-level compositional shifts to strain-level variation.

Differential Requirements for eIF4E Dose in Normal Development and Cancer.

eIF4E, the major cap-binding protein, has long been considered limiting for translating the mammalian genome. However, the eIF4E dose requirement at an organismal level remains unexplored. By generating an Eif4e haploinsufficient mouse, we found that a 50% reduction in eIF4E expression, while compatible with normal development and global protein synthesis, significantly impeded cellular transformation. Genome-wide translational profiling uncovered a translational program induced by oncogenic transformation and revealed a critical role for the dose of eIF4E, specifically in translating a network of mRNAs enriched for a unique 5' UTR signature. In particular, we demonstrate that the dose of eIF4E is essential for translating mRNAs that regulate reactive oxygen species, fueling transformation and cancer cell survival in vivo. Our findings indicate eIF4E is maintained at levels in excess for normal development that are hijacked by cancer cells to drive a translational program supporting tumorigenesis.

Transcription-replication interactions reveal bacterial genome regulation.

Organisms determine the transcription rates of thousands of genes through a few modes of regulation that recur across the genome1. In bacteria, the relationship between the regulatory architecture of a gene and its expression is well understood for individual model gene circuits2,3. However, a broader perspective of these dynamics at the genome scale is lacking, in part because bacterial transcriptomics has hitherto captured only a static snapshot of expression averaged across millions of cells4. As a result, the full diversity of gene expression dynamics and their relation to regulatory architecture remains unknown. Here we present a novel genome-wide classification of regulatory modes based on the transcriptional response of each gene to its own replication, which we term the transcription-replication interaction profile (TRIP). Analysing single-bacterium RNA-sequencing data, we found that the response to the universal perturbation of chromosomal replication integrates biological regulatory factors with biophysical molecular events on the chromosome to reveal the local regulatory context of a gene. Whereas the TRIPs of many genes conform to a gene dosage-dependent pattern, others diverge in distinct ways, and this is shaped by factors such as intra-operon position and repression state. By revealing the underlying mechanistic drivers of gene expression heterogeneity, this work provides a quantitative, biophysical framework for modelling replication-dependent expression dynamics.

Regulatory T Cell-Derived TGF-ß1 Controls Multiple Checkpoints Governing Allergy and Autoimmunity.

The mechanisms by which regulatory T (Treg) cells differentially control allergic and autoimmune responses remain unclear. We show that Treg cells in food allergy (FA) had decreased expression of transforming growth factor beta 1 (TGF-ß1) because of interleukin-4 (IL-4)- and signal transducer and activator of transciription-6 (STAT6)-dependent inhibition of Tgfb1 transcription. These changes were modeled by Treg cell-specific Tgfb1 monoallelic inactivation, which induced allergic dysregulation by impairing microbiota-dependent retinoic acid receptor-related orphan receptor gamma t (ROR-γt)+ Treg cell differentiation. This dysregulation was rescued by treatment with Clostridiales species, which upregulated Tgfb1 expression in Treg cells. Biallelic deficiency precipitated fatal autoimmunity with intense autoantibody production and dysregulated T follicular helper and B cell responses. These results identify a privileged role of Treg cell-derived TGF-ß1 in regulating allergy and autoimmunity at distinct checkpoints in a Tgfb1 gene dose- and microbiota-dependent manner.

The transcription factor titration effect dictates level of gene expression.

Models of transcription are often built around a picture of RNA polymerase and transcription factors (TFs) acting on a single copy of a promoter. However, most TFs are shared between multiple genes with varying binding affinities. Beyond that, genes often exist at high copy number-in multiple identical copies on the chromosome or on plasmids or viral vectors with copy numbers in the hundreds. Using a thermodynamic model, we characterize the interplay between TF copy number and the demand for that TF. We demonstrate the parameter-free predictive power of this model as a function of the copy number of the TF and the number and affinities of the available specific binding sites; such predictive control is important for the understanding of transcription and the desire to quantitatively design the output of genetic circuits. Finally, we use these experiments to dynamically measure plasmid copy number through the cell cycle.

Antibiotic-induced replication stress triggers bacterial competence by increasing gene dosage near the origin.

Streptococcus pneumoniae (pneumococcus) kills nearly 1 million children annually, and the emergence of antibiotic-resistant strains poses a serious threat to human health. Because pneumococci can take up DNA from their environment by a process called competence, genes associated with antibiotic resistance can rapidly spread. Remarkably, competence is activated in response to several antibiotics. Here, we demonstrate that antibiotics targeting DNA replication cause an increase in the copy number of genes proximal to the origin of replication (oriC). As the genes required for competence initiation are located near oriC, competence is thereby activated. Transcriptome analyses show that antibiotics targeting DNA replication also upregulate origin-proximal gene expression in other bacteria. This mechanism is a direct, intrinsic consequence of replication fork stalling. Our data suggest that evolution has conserved the oriC-proximal location of important genes in bacteria to allow for a robust response to replication stress without the need for complex gene-regulatory pathways. PAPERCLIP:

Mono- and biallelic variant effects on disease at biobank scale.

Identifying causal factors for Mendelian and common diseases is an ongoing challenge in medical genetics1. Population bottleneck events, such as those that occurred in the history of the Finnish population, enrich some homozygous variants to higher frequencies, which facilitates the identification of variants that cause diseases with recessive inheritance2,3. Here we examine the homozygous and heterozygous effects of 44,370 coding variants on 2,444 disease phenotypes using data from the nationwide electronic health records of 176,899 Finnish individuals. We find associations for homozygous genotypes across a broad spectrum of phenotypes, including known associations with retinal dystrophy and novel associations with adult-onset cataract and female infertility. Of the recessive disease associations that we identify, 13 out of 20 would have been missed by the additive model that is typically used in genome-wide association studies. We use these results to find many known Mendelian variants whose inheritance cannot be adequately described by a conventional definition of dominant or recessive. In particular, we find variants that are known to cause diseases with recessive inheritance with significant heterozygous phenotypic effects. Similarly, we find presumed benign variants with disease effects. Our results show how biobanks, particularly in founder populations, can broaden our understanding of complex dosage effects of Mendelian variants on disease.