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.

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.

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:

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.

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.

Germline MBD4 deficiency causes a multi-tumor predisposition syndrome.

We report an autosomal recessive, multi-organ tumor predisposition syndrome, caused by bi-allelic loss-of-function germline variants in the base excision repair (BER) gene MBD4. We identified five individuals with bi-allelic MBD4 variants within four families and these individuals had a personal and/or family history of adenomatous colorectal polyposis, acute myeloid leukemia, and uveal melanoma. MBD4 encodes a glycosylase involved in repair of G:T mismatches resulting from deamination of 5'-methylcytosine. The colorectal adenomas from MBD4-deficient individuals showed a mutator phenotype attributable to mutational signature SBS1, consistent with the function of MBD4. MBD4-deficient polyps harbored somatic mutations in similar driver genes to sporadic colorectal tumors, although AMER1 mutations were more common and KRAS mutations less frequent. Our findings expand the role of BER deficiencies in tumor predisposition. Inclusion of MBD4 in genetic testing for polyposis and multi-tumor phenotypes is warranted to improve disease management.

Prolidase Deficiency Causes Spontaneous T Cell Activation and Lupus-like Autoimmunity.

Prolidase deficiency (PD) is a multisystem disorder caused by mutations in the PEPD gene, which encodes a ubiquitously expressed metallopeptidase essential for the hydrolysis of dipeptides containing C-terminal proline or hydroxyproline. PD typically presents in childhood with developmental delay, skin ulcers, recurrent infections, and, in some patients, autoimmune features that can mimic systemic lupus erythematosus. The basis for the autoimmune association is uncertain, but might be due to self-antigen exposure with tissue damage, or indirectly driven by chronic infection and microbial burden. In this study, we address the question of causation and show that Pepd-null mice have increased antinuclear autoantibodies and raised serum IgA, accompanied by kidney immune complex deposition, consistent with a systemic lupus erythematosus-like disease. These features are associated with an accumulation of CD4 and CD8 effector T cells in the spleen and liver. Pepd deficiency leads to spontaneous T cell activation and proliferation into the effector subset, which is cell intrinsic and independent of Ag receptor specificity or antigenic stimulation. However, an increase in KLRG1+ effector CD8 cells is not observed in mixed chimeras, in which the autoimmune phenotype is also absent. Our findings link autoimmune susceptibility in PD to spontaneous T cell dysfunction, likely to be acting in combination with immune activators that lie outside the hemopoietic system but result from the abnormal metabolism or loss of nonenzymatic prolidase function. This knowledge provides insight into the role of prolidase in the maintenance of self-tolerance and highlights the importance of treatment to control T cell activation.

Multitarget nociceptor sensitization by a promiscuous peptide from the venom of the King Baboon spider.

The King Baboon spider, Pelinobius muticus, is a burrowing African tarantula. Its impressive size and appealing coloration are tempered by reports describing severe localized pain, swelling, itchiness, and muscle cramping after accidental envenomation. Hyperalgesia is the most prominent symptom after bites from P. muticus, but the molecular basis by which the venom induces pain is unknown. Proteotranscriptomic analysis of P. muticus venom uncovered a cysteine-rich peptide, δ/κ-theraphotoxin-Pm1a (δ/κ-TRTX-Pm1a), that elicited nocifensive behavior when injected into mice. In small dorsal root ganglion neurons, synthetic δ/κ-TRTX-Pm1a (sPm1a) induced hyperexcitability by enhancing tetrodotoxin-resistant sodium currents, impairing repolarization and lowering the threshold of action potential firing, consistent with the severe pain associated with envenomation. The molecular mechanism of nociceptor sensitization by sPm1a involves multimodal actions over several ion channel targets, including NaV1.8, KV2.1, and tetrodotoxin-sensitive NaV channels. The promiscuous targeting of peptides like δ/κ-TRTX-Pm1a may be an evolutionary adaptation in pain-inducing defensive venoms.

Defining novel causal SNPs and linked phenotypes at melanoma-associated loci.

A number of genomic regions have been associated with melanoma risk through genome-wide association studies; however, the causal variants underlying the majority of these associations remain unknown. Here, we sequenced either the full locus or the functional regions including exons of 19 melanoma-associated loci in 1959 British melanoma cases and 737 controls. Variant filtering followed by Fisher's exact test analyses identified 66 variants associated with melanoma risk. Sequential conditional logistic regression identified the distinct haplotypes on which variants reside, and massively parallel reporter assays provided biological insights into how these variants influence gene function. We performed further analyses to link variants to melanoma risk phenotypes and assessed their association with melanoma-specific survival. Our analyses replicate previously known associations in the melanocortin 1 receptor (MC1R) and tyrosinase (TYR) loci, while identifying novel potentially causal variants at the MTAP/CDKN2A and CASP8 loci. These results improve our understanding of the architecture of melanoma risk and outcome.

In situ CRISPR-Cas9 base editing for the development of genetically engineered mouse models of breast cancer.

Genetically engineered mouse models (GEMMs) of cancer have proven to be of great value for basic and translational research. Although CRISPR-based gene disruption offers a fast-track approach for perturbing gene function and circumvents certain limitations of standard GEMM development, it does not provide a flexible platform for recapitulating clinically relevant missense mutations in vivo. To this end, we generated knock-in mice with Cre-conditional expression of a cytidine base editor and tested their utility for precise somatic engineering of missense mutations in key cancer drivers. Upon intraductal delivery of sgRNA-encoding vectors, we could install point mutations with high efficiency in one or multiple endogenous genes in situ and assess the effect of defined allelic variants on mammary tumorigenesis. While the system also produces bystander insertions and deletions that can stochastically be selected for when targeting a tumor suppressor gene, we could effectively recapitulate oncogenic nonsense mutations. We successfully applied this system in a model of triple-negative breast cancer, providing the proof of concept for extending this flexible somatic base editing platform to other tissues and tumor types.

CRLF3 plays a key role in the final stage of platelet genesis and is a potential therapeutic target for thrombocythemia.

The process of platelet production has so far been understood to be a 2-stage process: megakaryocyte maturation from hematopoietic stem cells followed by proplatelet formation, with each phase regulating the peripheral blood platelet count. Proplatelet formation releases into the bloodstream beads-on-a-string preplatelets, which undergo fission into mature platelets. For the first time, we show that preplatelet maturation is a third, tightly regulated, critical process akin to cytokinesis that regulates platelet count. We show that deficiency in cytokine receptor-like factor 3 (CRLF3) in mice leads to an isolated and sustained 25% to 48% reduction in the platelet count without any effect on other blood cell lineages. We show that Crlf3-/- preplatelets have increased microtubule stability, possibly because of increased microtubule glutamylation via the interaction of CRLF3 with key members of the Hippo pathway. Using a mouse model of JAK2 V617F essential thrombocythemia, we show that a lack of CRLF3 leads to long-term lineage-specific normalization of the platelet count. We thereby postulate that targeting CRLF3 has therapeutic potential for treatment of thrombocythemia.

Breathing in danger: how particulate matter pollution is putting the public at risk of lung cancer†.

We are constantly exposed to chemicals and other agents in our environment that can influence our risk of tumorigenesis, but exactly how these factors contribute to cancer development is largely unknown. Fine particulate matter measuring ≤2.5 µm (PM2.5 ) from air pollution can accumulate in alveoli, contributing to inflammation and tissue damage. Despite prior correlative studies highlighting the mortality risk, there has been a historical reluctance to lower national standards for safe PM2.5 exposure. A recent publication further highlights the attributable risk of PM2.5 exposure with lung cancer - particularly in 'never-smokers' with EGFR-driven non-small cell lung cancer. Importantly, it also elucidates a mechanistic link between PM2.5 exposure and tumorigenesis using in vivo models of EGFR non-small cell lung cancer. © 2023 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Mutually exclusive genetic interactions and gene essentiality shape the genomic landscape of primary melanoma.

Melanoma is a heterogenous malignancy with an unpredictable clinical course. Most patients who present in the clinic are diagnosed with primary melanoma, yet large-scale sequencing efforts have focused primarily on metastatic disease. In this study we sequence-profiled 524 American Joint Committee on Cancer Stage I-III primary tumours. Our analysis of these data reveals recurrent driver mutations, mutually exclusive genetic interactions, where two genes were never or rarely co-mutated, and an absence of co-occurring genetic events. Further, we intersected copy number calls from our primary melanoma data with whole-genome CRISPR screening data to identify the transcription factor interferon regulatory factor 4 (IRF4) as a melanoma-associated dependency. © 2022 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

The IFITM proteins mediate cellular resistance to influenza A H1N1 virus, West Nile virus, and dengue virus.

Influenza viruses exploit host cell machinery to replicate, resulting in epidemics of respiratory illness. In turn, the host expresses antiviral restriction factors to defend against infection. To find host cell modifiers of influenza A H1N1 viral infection, we used a functional genomic screen and identified over 120 influenza A virus-dependency factors with roles in endosomal acidification, vesicular trafficking, mitochondrial metabolism, and RNA splicing. We discovered that the interferon-inducible transmembrane proteins IFITM1, 2, and 3 restrict an early step in influenza A viral replication. The IFITM proteins confer basal resistance to influenza A virus but are also inducible by interferons type I and II and are critical for interferon's virustatic actions. Further characterization revealed that the IFITM proteins inhibit the early replication of flaviviruses, including dengue virus and West Nile virus. Collectively this work identifies a family of antiviral restriction factors that mediate cellular innate immunity to at least three major human pathogens.

Placentation defects are highly prevalent in embryonic lethal mouse mutants.

Large-scale phenotyping efforts have demonstrated that approximately 25-30% of mouse gene knockouts cause intrauterine lethality. Analysis of these mutants has largely focused on the embryo and not the placenta, despite the crucial role of this extraembryonic organ for developmental progression. Here we screened 103 embryonic lethal and sub-viable mouse knockout lines from the Deciphering the Mechanisms of Developmental Disorders program for placental phenotypes. We found that 68% of knockout lines that are lethal at or after mid-gestation exhibited placental dysmorphologies. Early lethality (embryonic days 9.5-14.5) is almost always associated with severe placental malformations. Placental defects correlate strongly with abnormal brain, heart and vascular development. Analysis of mutant trophoblast stem cells and conditional knockouts suggests that a considerable number of factors that cause embryonic lethality when ablated have primary gene function in trophoblast cells. Our data highlight the hugely under-appreciated importance of placental defects in contributing to abnormal embryo development and suggest key molecular nodes that govern placenta formation.

Population-based analysis of POT1 variants in a cutaneous melanoma case-control cohort.

Pathogenic germline variants in the protection of telomeres 1 gene (POT1) have been associated with predisposition to a range of tumour types, including melanoma, glioma, leukaemia and cardiac angiosarcoma. We sequenced all coding exons of the POT1 gene in 2928 European-descent melanoma cases and 3298 controls, identifying 43 protein-changing genetic variants. We performed POT1-telomere binding assays for all missense and stop-gained variants, finding nine variants that impair or disrupt protein-telomere complex formation, and we further define the role of variants in the regulation of telomere length and complex formation through molecular dynamics simulations. We determine that POT1 coding variants are a minor contributor to melanoma burden in the general population, with only about 0.5% of melanoma cases carrying germline pathogenic variants in this gene, but should be screened in individuals with a strong family history of melanoma and/or multiple malignancies.

Rational Design of Potent α-Conotoxin PeIA Analogues with Non-Natural Amino Acids for the Inhibition of Human α9α10 Nicotinic Acetylcholine Receptors.

α-Conotoxins (α-CTxs) are structurally related peptides that antagonize nicotinic acetylcholine receptors (nAChRs), which may serve as new alternatives to opioid-based treatment for pain-related conditions. The non-natural amino acid analogues of α-CTxs have been demonstrated with improved potency compared to the native peptide. In this study, we chemically synthesized Dab/Dap-substituted analogues of α-CTx PeIA and evaluated their activity at heterologously expressed human α9α10 nAChRs. PeIA[S4Dap, S9Dap] had the most potent half-maximal inhibitory concentration (IC50) of 0.93 nM. Molecular dynamic simulations suggested that the side chain amino group of Dap4 formed additional hydrogen bonds with S168 and D169 of the receptor and Dap9 formed an extra hydrogen bond interaction with Q34, which is distinctive to PeIA. Overall, our findings provide new insights into further development of more potent analogues of α-CTxs, and PeIA[S4Dap, S9Dap] has potential as a drug candidate for the treatment of chronic neuropathic pain.