Characterizing genetic intra-tumor heterogeneity across 2,658 human cancer genomes.

Intra-tumor heterogeneity (ITH) is a mechanism of therapeutic resistance and therefore an important clinical challenge. However, the extent, origin, and drivers of ITH across cancer types are poorly understood. To address this, we extensively characterize ITH across whole-genome sequences of 2,658 cancer samples spanning 38 cancer types. Nearly all informative samples (95.1%) contain evidence of distinct subclonal expansions with frequent branching relationships between subclones. We observe positive selection of subclonal driver mutations across most cancer types and identify cancer type-specific subclonal patterns of driver gene mutations, fusions, structural variants, and copy number alterations as well as dynamic changes in mutational processes between subclonal expansions. Our results underline the importance of ITH and its drivers in tumor evolution and provide a pan-cancer resource of comprehensively annotated subclonal events from whole-genome sequencing data.

Proteogenomic Characterization of Endometrial Carcinoma.

We undertook a comprehensive proteogenomic characterization of 95 prospectively collected endometrial carcinomas, comprising 83 endometrioid and 12 serous tumors. This analysis revealed possible new consequences of perturbations to the p53 and Wnt/ß-catenin pathways, identified a potential role for circRNAs in the epithelial-mesenchymal transition, and provided new information about proteomic markers of clinical and genomic tumor subgroups, including relationships to known druggable pathways. An extensive genome-wide acetylation survey yielded insights into regulatory mechanisms linking Wnt signaling and histone acetylation. We also characterized aspects of the tumor immune landscape, including immunogenic alterations, neoantigens, common cancer/testis antigens, and the immune microenvironment, all of which can inform immunotherapy decisions. Collectively, our multi-omic analyses provide a valuable resource for researchers and clinicians, identify new molecular associations of potential mechanistic significance in the development of endometrial cancers, and suggest novel approaches for identifying potential therapeutic targets.

Identification, Characterization, and Heritability of Murine Metastable Epialleles: Implications for Non-genetic Inheritance.

Generally repressed by epigenetic mechanisms, retrotransposons represent around 40% of the murine genome. At the Agouti viable yellow (Avy) locus, an endogenous retrovirus (ERV) of the intracisternal A particle (IAP) class retrotransposed upstream of the agouti coat-color locus, providing an alternative promoter that is variably DNA methylated in genetically identical individuals. This results in variable expressivity of coat color that is inherited transgenerationally. Here, a systematic genome-wide screen identifies multiple C57BL/6J murine IAPs with Avy epigenetic properties. Each exhibits a stable methylation state within an individual but varies between individuals. Only in rare instances do they act as promoters controlling adjacent gene expression. Their methylation state is locus-specific within an individual, and their flanking regions are enriched for CTCF. Variably methylated IAPs are reprogrammed after fertilization and re-established as variable loci in the next generation, indicating reconstruction of metastable epigenetic states and challenging the generalizability of non-genetic inheritance at these regions.

The BRAF pseudogene functions as a competitive endogenous RNA and induces lymphoma in vivo.

Research over the past decade has suggested important roles for pseudogenes in physiology and disease. In vitro experiments demonstrated that pseudogenes contribute to cell transformation through several mechanisms. However, in vivo evidence for a causal role of pseudogenes in cancer development is lacking. Here, we report that mice engineered to overexpress either the full-length murine B-Raf pseudogene Braf-rs1 or its pseudo "CDS" or "3' UTR" develop an aggressive malignancy resembling human diffuse large B cell lymphoma. We show that Braf-rs1 and its human ortholog, BRAFP1, elicit their oncogenic activity, at least in part, as competitive endogenous RNAs (ceRNAs) that elevate BRAF expression and MAPK activation in vitro and in vivo. Notably, we find that transcriptional or genomic aberrations of BRAFP1 occur frequently in multiple human cancers, including B cell lymphomas. Our engineered mouse models demonstrate the oncogenic potential of pseudogenes and indicate that ceRNA-mediated microRNA sequestration may contribute to the development of cancer.

De novo variants in the RNU4-2 snRNA cause a frequent neurodevelopmental syndrome.

Around 60% of individuals with neurodevelopmental disorders (NDD) remain undiagnosed after comprehensive genetic testing, primarily of protein-coding genes1. Large genome-sequenced cohorts are improving our ability to discover new diagnoses in the non-coding genome. Here, we identify the non-coding RNA RNU4-2 as a syndromic NDD gene. RNU4-2 encodes the U4 small nuclear RNA (snRNA), which is a critical component of the U4/U6.U5 tri-snRNP complex of the major spliceosome2. We identify an 18 bp region of RNU4-2 mapping to two structural elements in the U4/U6 snRNA duplex (the T-loop and Stem III) that is severely depleted of variation in the general population, but in which we identify heterozygous variants in 115 individuals with NDD. Most individuals (77.4%) have the same highly recurrent single base insertion (n.64_65insT). In 54 individuals where it could be determined, the de novo variants were all on the maternal allele. We demonstrate that RNU4-2 is highly expressed in the developing human brain, in contrast to RNU4-1 and other U4 homologs. Using RNA-sequencing, we show how 5' splice site usage is systematically disrupted in individuals with RNU4-2 variants, consistent with the known role of this region during spliceosome activation. Finally, we estimate that variants in this 18 bp region explain 0.4% of individuals with NDD. This work underscores the importance of non-coding genes in rare disorders and will provide a diagnosis to thousands of individuals with NDD worldwide.

Autonomous healing of fatigue cracks via cold welding.

Fatigue in metals involves gradual failure through incremental propagation of cracks under repetitive mechanical load. In structural applications, fatigue accounts for up to 90% of in-service failure1,2. Prevention of fatigue relies on implementation of large safety factors and inefficient overdesign3. In traditional metallurgical design for fatigue resistance, microstructures are developed to either arrest or slow the progression of cracks. Crack growth is assumed to be irreversible. By contrast, in other material classes, there is a compelling alternative based on latent healing mechanisms and damage reversal4-9. Here, we report that fatigue cracks in pure metals can undergo intrinsic self-healing. We directly observe the early progression of nanoscale fatigue cracks, and as expected, the cracks advance, deflect and arrest at local microstructural barriers. However, unexpectedly, cracks were also observed to heal by a process that can be described as crack flank cold welding induced by a combination of local stress state and grain boundary migration. The premise that fatigue cracks can autonomously heal in metals through local interaction with microstructural features challenges the most fundamental theories on how engineers design and evaluate fatigue life in structural materials. We discuss the implications for fatigue in a variety of service environments.

Cooperativity and rapid evolution of cobound transcription factors in closely related mammals.

To mechanistically characterize the microevolutionary processes active in altering transcription factor (TF) binding among closely related mammals, we compared the genome-wide binding of three tissue-specific TFs that control liver gene expression in six rodents. Despite an overall fast turnover of TF binding locations between species, we identified thousands of TF regions of highly constrained TF binding intensity. Although individual mutations in bound sequence motifs can influence TF binding, most binding differences occur in the absence of nearby sequence variations. Instead, combinatorial binding was found to be significant for genetic and evolutionary stability; cobound TFs tend to disappear in concert and were sensitive to genetic knockout of partner TFs. The large, qualitative differences in genomic regions bound between closely related mammals, when contrasted with the smaller, quantitative TF binding differences among Drosophila species, illustrate how genome structure and population genetics together shape regulatory evolution.

Genome sequencing reveals loci under artificial selection that underlie disease phenotypes in the laboratory rat.

Large numbers of inbred laboratory rat strains have been developed for a range of complex disease phenotypes. To gain insights into the evolutionary pressures underlying selection for these phenotypes, we sequenced the genomes of 27 rat strains, including 11 models of hypertension, diabetes, and insulin resistance, along with their respective control strains. Altogether, we identified more than 13 million single-nucleotide variants, indels, and structural variants across these rat strains. Analysis of strain-specific selective sweeps and gene clusters implicated genes and pathways involved in cation transport, angiotensin production, and regulators of oxidative stress in the development of cardiovascular disease phenotypes in rats. Many of the rat loci that we identified overlap with previously mapped loci for related traits in humans, indicating the presence of shared pathways underlying these phenotypes in rats and humans. These data represent a step change in resources available for evolutionary analysis of complex traits in disease models.

Genome-wide generation and systematic phenotyping of knockout mice reveals new roles for many genes.

Mutations in whole organisms are powerful ways of interrogating gene function in a realistic context. We describe a program, the Sanger Institute Mouse Genetics Project, that provides a step toward the aim of knocking out all genes and screening each line for a broad range of traits. We found that hitherto unpublished genes were as likely to reveal phenotypes as known genes, suggesting that novel genes represent a rich resource for investigating the molecular basis of disease. We found many unexpected phenotypes detected only because we screened for them, emphasizing the value of screening all mutants for a wide range of traits. Haploinsufficiency and pleiotropy were both surprisingly common. Forty-two percent of genes were essential for viability, and these were less likely to have a paralog and more likely to contribute to a protein complex than other genes. Phenotypic data and more than 900 mutants are openly available for further analysis. PAPERCLIP:

A polymorphism in IRF4 affects human pigmentation through a tyrosinase-dependent MITF/TFAP2A pathway.

Sequence polymorphisms linked to human diseases and phenotypes in genome-wide association studies often affect noncoding regions. A SNP within an intron of the gene encoding Interferon Regulatory Factor 4 (IRF4), a transcription factor with no known role in melanocyte biology, is strongly associated with sensitivity of skin to sun exposure, freckles, blue eyes, and brown hair color. Here, we demonstrate that this SNP lies within an enhancer of IRF4 transcription in melanocytes. The allele associated with this pigmentation phenotype impairs binding of the TFAP2A transcription factor that, together with the melanocyte master regulator MITF, regulates activity of the enhancer. Assays in zebrafish and mice reveal that IRF4 cooperates with MITF to activate expression of Tyrosinase (TYR), an essential enzyme in melanin synthesis. Our findings provide a clear example of a noncoding polymorphism that affects a phenotype by modulating a developmental gene regulatory network.

Two defence systems eliminate plasmids from seventh pandemic Vibrio cholerae.

Horizontal gene transfer can trigger rapid shifts in bacterial evolution. Driven by a variety of mobile genetic elements-in particular bacteriophages and plasmids-the ability to share genes within and across species underpins the exceptional adaptability of bacteria. Nevertheless, invasive mobile genetic elements can also present grave risks to the host; bacteria have therefore evolved a vast array of defences against these elements1. Here we identify two plasmid defence systems conserved in the Vibrio cholerae El Tor strains responsible for the ongoing seventh cholera pandemic2-4. These systems, termed DdmABC and DdmDE, are encoded on two major pathogenicity islands that are a hallmark of current pandemic strains. We show that the modules cooperate to rapidly eliminate small multicopy plasmids by degradation. Moreover, the DdmABC system is widespread and can defend against bacteriophage infection by triggering cell suicide (abortive infection, or Abi). Notably, we go on to show that, through an Abi-like mechanism, DdmABC increases the burden of large low-copy-number conjugative plasmids, including a broad-host IncC multidrug resistance plasmid, which creates a fitness disadvantage that counterselects against plasmid-carrying cells. Our results answer the long-standing question of why plasmids, although abundant in environmental strains, are rare in pandemic strains; have implications for understanding the dissemination of antibiotic resistance plasmids; and provide insights into how the interplay between two defence systems has shaped the evolution of the most successful lineage of pandemic V. cholerae.

Identification of bacteria-derived HLA-bound peptides in melanoma.

A variety of species of bacteria are known to colonize human tumours1-11, proliferate within them and modulate immune function, which ultimately affects the survival of patients with cancer and their responses to treatment12-14. However, it is not known whether antigens derived from intracellular bacteria are presented by the human leukocyte antigen class I and II (HLA-I and HLA-II, respectively) molecules of tumour cells, or whether such antigens elicit a tumour-infiltrating T cell immune response. Here we used 16S rRNA gene sequencing and HLA peptidomics to identify a peptide repertoire derived from intracellular bacteria that was presented on HLA-I and HLA-II molecules in melanoma tumours. Our analysis of 17 melanoma metastases (derived from 9 patients) revealed 248 and 35 unique HLA-I and HLA-II peptides, respectively, that were derived from 41 species of bacteria. We identified recurrent bacterial peptides in tumours from different patients, as well as in different tumours from the same patient. Our study reveals that peptides derived from intracellular bacteria can be presented by tumour cells and elicit immune reactivity, and thus provides insight into a mechanism by which bacteria influence activation of the immune system and responses to therapy.

Loss of the endoplasmic reticulum protein Tmem208 affects cell polarity, development, and viability.

Nascent proteins destined for the cell membrane and the secretory pathway are targeted to the endoplasmic reticulum (ER) either posttranslationally or cotranslationally. The signal-independent pathway, containing the protein TMEM208, is one of three pathways that facilitates the translocation of nascent proteins into the ER. The in vivo function of this protein is ill characterized in multicellular organisms. Here, we generated a CRISPR-induced null allele of the fruit fly ortholog CG8320/Tmem208 by replacing the gene with the Kozak-GAL4 sequence. We show that Tmem208 is broadly expressed in flies and that its loss causes lethality, although a few short-lived flies eclose. These animals exhibit wing and eye developmental defects consistent with impaired cell polarity and display mild ER stress. Tmem208 physically interacts with Frizzled (Fz), a planar cell polarity (PCP) receptor, and is required to maintain proper levels of Fz. Moreover, we identified a child with compound heterozygous variants in TMEM208 who presents with developmental delay, skeletal abnormalities, multiple hair whorls, cardiac, and neurological issues, symptoms that are associated with PCP defects in mice and humans. Additionally, fibroblasts of the proband display mild ER stress. Expression of the reference human TMEM208 in flies fully rescues the loss of Tmem208, and the two proband-specific variants fail to rescue, suggesting that they are loss-of-function alleles. In summary, our study uncovers a role of TMEM208 in development, shedding light on its significance in ER homeostasis and cell polarity.

Haplotype-based analysis resolves missing heritability in oculocutaneous albinism type 1B.

Oculocutaneous albinism (OCA) is a rare disorder of pigment production. Affected individuals have variably decreased global pigmentation and visual-developmental changes that lead to low vision. OCA is notable for significant missing heritability, particularly among individuals with residual pigmentation. Tyrosinase (TYR) is the rate-limiting enzyme in melanin pigment biosynthesis and mutations that decrease enzyme function are one of the most common causes of OCA. We present the analysis of high-depth short-read TYR sequencing data for a cohort of 352 OCA probands, ∼50% of whom were previously sequenced without yielding a definitive diagnostic result. Our analysis identified 66 TYR single-nucleotide variants (SNVs) and small insertion/deletions (indels), 3 structural variants, and a rare haplotype comprised of two common frequency variants (p.Ser192Tyr and p.Arg402Gln) in cis-orientation, present in 149/352 OCA probands. We further describe a detailed analysis of the disease-causing haplotype, p.[Ser192Tyr; Arg402Gln] ("cis-YQ"). Haplotype analysis suggests that the cis-YQ allele arose by recombination and that multiple cis-YQ haplotypes are segregating in OCA-affected individuals and control populations. The cis-YQ allele is the most common disease-causing allele in our cohort, representing 19.1% (57/298) of TYR pathogenic alleles in individuals with type 1 (TYR-associated) OCA. Finally, among the 66 TYR variants, we found several additional alleles defined by a cis-oriented combination of minor, potentially hypomorph-producing alleles at common variant sites plus a second, rare pathogenic variant. Together, these results suggest that identification of phased variants for the full TYR locus are required for an exhaustive assessment for potentially disease-causing alleles.

Bi-allelic SNAPC4 variants dysregulate global alternative splicing and lead to neuroregression and progressive spastic paraparesis.

The vast majority of human genes encode multiple isoforms through alternative splicing, and the temporal and spatial regulation of those isoforms is critical for organismal development and function. The spliceosome, which regulates and executes splicing reactions, is primarily composed of small nuclear ribonucleoproteins (snRNPs) that consist of small nuclear RNAs (snRNAs) and protein subunits. snRNA gene transcription is initiated by the snRNA-activating protein complex (SNAPc). Here, we report ten individuals, from eight families, with bi-allelic, deleterious SNAPC4 variants. SNAPC4 encoded one of the five SNAPc subunits that is critical for DNA binding. Most affected individuals presented with delayed motor development and developmental regression after the first year of life, followed by progressive spasticity that led to gait alterations, paraparesis, and oromotor dysfunction. Most individuals had cerebral, cerebellar, or basal ganglia volume loss by brain MRI. In the available cells from affected individuals, SNAPC4 abundance was decreased compared to unaffected controls, suggesting that the bi-allelic variants affect SNAPC4 accumulation. The depletion of SNAPC4 levels in HeLa cell lines via genomic editing led to decreased snRNA expression and global dysregulation of alternative splicing. Analysis of available fibroblasts from affected individuals showed decreased snRNA expression and global dysregulation of alternative splicing compared to unaffected cells. Altogether, these data suggest that these bi-allelic SNAPC4 variants result in loss of function and underlie the neuroregression and progressive spasticity in these affected individuals.

Transcatheter Repair for Patients with Tricuspid Regurgitation.

BACKGROUND: Severe tricuspid regurgitation is a debilitating condition that is associated with substantial morbidity and often with poor quality of life. Decreasing tricuspid regurgitation may reduce symptoms and improve clinical outcomes in patients with this disease. METHODS: We conducted a prospective randomized trial of percutaneous tricuspid transcatheter edge-to-edge repair (TEER) for severe tricuspid regurgitation. Patients with symptomatic severe tricuspid regurgitation were enrolled at 65 centers in the United States, Canada, and Europe and were randomly assigned in a 1:1 ratio to receive either TEER or medical therapy (control). The primary end point was a hierarchical composite that included death from any cause or tricuspid-valve surgery; hospitalization for heart failure; and an improvement in quality of life as measured with the Kansas City Cardiomyopathy Questionnaire (KCCQ), with an improvement defined as an increase of at least 15 points in the KCCQ score (range, 0 to 100, with higher scores indicating better quality of life) at the 1-year follow-up. The severity of tricuspid regurgitation and safety were also assessed. RESULTS: A total of 350 patients were enrolled; 175 were assigned to each group. The mean age of the patients was 78 years, and 54.9% were women. The results for the primary end point favored the TEER group (win ratio, 1.48; 95% confidence interval, 1.06 to 2.13; P = 0.02). The incidence of death or tricuspid-valve surgery and the rate of hospitalization for heart failure did not appear to differ between the groups. The KCCQ quality-of-life score changed by a mean (±SD) of 12.3±1.8 points in the TEER group, as compared with 0.6±1.8 points in the control group (P<0.001). At 30 days, 87.0% of the patients in the TEER group and 4.8% of those in the control group had tricuspid regurgitation of no greater than moderate severity (P<0.001). TEER was found to be safe; 98.3% of the patients who underwent the procedure were free from major adverse events at 30 days. CONCLUSIONS: Tricuspid TEER was safe for patients with severe tricuspid regurgitation, reduced the severity of tricuspid regurgitation, and was associated with an improvement in quality of life. (Funded by Abbott; TRILUMINATE Pivotal ClinicalTrials.gov number, NCT03904147.).

A previously unrecognized superfamily of macro-conotoxins includes an inhibitor of the sensory neuron calcium channel Cav2.3.

Animal venom peptides represent valuable compounds for biomedical exploration. The venoms of marine cone snails constitute a particularly rich source of peptide toxins, known as conotoxins. Here, we identify the sequence of an unusually large conotoxin, Mu8.1, which defines a new class of conotoxins evolutionarily related to the well-known con-ikot-ikots and 2 additional conotoxin classes not previously described. The crystal structure of recombinant Mu8.1 displays a saposin-like fold and shows structural similarity with con-ikot-ikot. Functional studies demonstrate that Mu8.1 curtails calcium influx in defined classes of murine somatosensory dorsal root ganglion (DRG) neurons. When tested on a variety of recombinantly expressed voltage-gated ion channels, Mu8.1 displayed the highest potency against the R-type (Cav2.3) calcium channel. Ca2+ signals from Mu8.1-sensitive DRG neurons were also inhibited by SNX-482, a known spider peptide modulator of Cav2.3 and voltage-gated K+ (Kv4) channels. Our findings highlight the potential of Mu8.1 as a molecular tool to identify and study neuronal subclasses expressing Cav2.3. Importantly, this multidisciplinary study showcases the potential of uncovering novel structures and bioactivities within the largely unexplored group of macro-conotoxins.

In vivo identification of tumor- suppressive PTEN ceRNAs in an oncogenic BRAF-induced mouse model of melanoma.

We recently proposed that competitive endogenous RNAs (ceRNAs) sequester microRNAs to regulate mRNA transcripts containing common microRNA recognition elements (MREs). However, the functional role of ceRNAs in cancer remains unknown. Loss of PTEN, a tumor suppressor regulated by ceRNA activity, frequently occurs in melanoma. Here, we report the discovery of significant enrichment of putative PTEN ceRNAs among genes whose loss accelerates tumorigenesis following Sleeping Beauty insertional mutagenesis in a mouse model of melanoma. We validated several putative PTEN ceRNAs and further characterized one, the ZEB2 transcript. We show that ZEB2 modulates PTEN protein levels in a microRNA-dependent, protein coding-independent manner. Attenuation of ZEB2 expression activates the PI3K/AKT pathway, enhances cell transformation, and commonly occurs in human melanomas and other cancers expressing low PTEN levels. Our study genetically identifies multiple putative microRNA decoys for PTEN, validates ZEB2 mRNA as a bona fide PTEN ceRNA, and demonstrates that abrogated ZEB2 expression cooperates with BRAF(V600E) to promote melanomagenesis.

Repeat expansions nested within tandem CNVs: a unique structural change in GLS exemplifies the diagnostic challenges of non-coding pathogenic variation.

Glutaminase deficiency has recently been associated with ataxia and developmental delay due to repeat expansions in the 5'UTR of the glutaminase (GLS) gene. Patients with the described GLS repeat expansion may indeed remain undiagnosed due to the rarity of this variant, the challenge of its detection and the recency of its discovery. In this study, we combined advanced bioinformatics screening of ~3000 genomes and ~1500 exomes with optical genome mapping and long-read sequencing for confirmation studies. We identified two GLS families, previously intensely and unsuccessfully analyzed. One family carries an unusual and complex structural change involving a homozygous repeat expansion nested within a quadruplication event in the 5'UTR of GLS. Glutaminase deficiency and its metabolic consequences were validated by in-depth biochemical analysis. The identified GLS patients showed progressive early-onset ataxia, cognitive deficits, pyramidal tract damage and optic atrophy, thus demonstrating susceptibility of several specific neuron populations to glutaminase deficiency. This large-scale screening study demonstrates the ability of bioinformatics analysis-validated by latest state-of-the-art technologies (optical genome mapping and long-read sequencing)-to effectively flag complex repeat expansions using short-read datasets and thus facilitate diagnosis of ultra-rare disorders.