A curated census of pathogenic and likely pathogenic UTR variants and evaluation of deep learning models for variant effect prediction.

Introduction: Variants in 5' and 3' untranslated regions (UTR) contribute to rare disease. While predictive algorithms to assist in classifying pathogenicity can potentially be highly valuable, the utility of these tools is often unclear, as it depends on carefully selected training and validation conditions. To address this, we developed a high confidence set of pathogenic (P) and likely pathogenic (LP) variants and assessed deep learning (DL) models for predicting their molecular effects. Methods: 3' and 5' UTR variants documented as P or LP (P/LP) were obtained from ClinVar and refined by reviewing the annotated variant effect and reassessing evidence of pathogenicity following published guidelines. Prediction scores from sequence-based DL models were compared between three groups: P/LP variants acting though the mechanism for which the model was designed (model-matched), those operating through other mechanisms (model-mismatched), and putative benign variants. PhyloP was used to compare conservation scores between P/LP and putative benign variants. Results: 295 3' and 188 5' UTR variants were obtained from ClinVar, of which 26 3' and 68 5' UTR variants were classified as P/LP. Predictions by DL models achieved statistically significant differences when comparing modelmatched P/LP variants to both putative benign variants and modelmismatched P/LP variants, as well as when comparing all P/LP variants to putative benign variants. PhyloP conservation scores were significantly higher among P/LP compared to putative benign variants for both the 3' and 5' UTR. Discussion: In conclusion, we present a high-confidence set of P/LP 3' and 5' UTR variants spanning a range of mechanisms and supported by detailed pathogenicity and molecular mechanism evidence curation. Predictions from DL models further substantiate these classifications. These datasets will support further development and validation of DL algorithms designed to predict the functional impact of variants that may be implicated in rare disease.

Single-cell genomic variation induced by mutational processes in cancer.

How cell-to-cell copy number alterations that underpin genomic instability1 in human cancers drive genomic and phenotypic variation, and consequently the evolution of cancer2, remains understudied. Here, by applying scaled single-cell whole-genome sequencing3 to wild-type, TP53-deficient and TP53-deficient;BRCA1-deficient or TP53-deficient;BRCA2-deficient mammary epithelial cells (13,818 genomes), and to primary triple-negative breast cancer (TNBC) and high-grade serous ovarian cancer (HGSC) cells (22,057 genomes), we identify three distinct 'foreground' mutational patterns that are defined by cell-to-cell structural variation. Cell- and clone-specific high-level amplifications, parallel haplotype-specific copy number alterations and copy number segment length variation (serrate structural variations) had measurable phenotypic and evolutionary consequences. In TNBC and HGSC, clone-specific high-level amplifications in known oncogenes were highly prevalent in tumours bearing fold-back inversions, relative to tumours with homologous recombination deficiency, and were associated with increased clone-to-clone phenotypic variation. Parallel haplotype-specific alterations were also commonly observed, leading to phylogenetic evolutionary diversity and clone-specific mono-allelic expression. Serrate variants were increased in tumours with fold-back inversions and were highly correlated with increased genomic diversity of cellular populations. Together, our findings show that cell-to-cell structural variation contributes to the origins of phenotypic and evolutionary diversity in TNBC and HGSC, and provide insight into the genomic and mutational states of individual cancer cells.

Clonal fitness inferred from time-series modelling of single-cell cancer genomes.

Progress in defining genomic fitness landscapes in cancer, especially those defined by copy number alterations (CNAs), has been impeded by lack of time-series single-cell sampling of polyclonal populations and temporal statistical models1-7. Here we generated 42,000 genomes from multi-year time-series single-cell whole-genome sequencing of breast epithelium and primary triple-negative breast cancer (TNBC) patient-derived xenografts (PDXs), revealing the nature of CNA-defined clonal fitness dynamics induced by TP53 mutation and cisplatin chemotherapy. Using a new Wright-Fisher population genetics model8,9 to infer clonal fitness, we found that TP53 mutation alters the fitness landscape, reproducibly distributing fitness over a larger number of clones associated with distinct CNAs. Furthermore, in TNBC PDX models with mutated TP53, inferred fitness coefficients from CNA-based genotypes accurately forecast experimentally enforced clonal competition dynamics. Drug treatment in three long-term serially passaged TNBC PDXs resulted in cisplatin-resistant clones emerging from low-fitness phylogenetic lineages in the untreated setting. Conversely, high-fitness clones from treatment-naive controls were eradicated, signalling an inversion of the fitness landscape. Finally, upon release of drug, selection pressure dynamics were reversed, indicating a fitness cost of treatment resistance. Together, our findings define clonal fitness linked to both CNA and therapeutic resistance in polyclonal tumours.

Ubiquitin-mediated DNA damage response is synthetic lethal with G-quadruplex stabilizer CX-5461.

CX-5461 is a G-quadruplex (G4) ligand currently in trials with initial indications of clinical activity in cancers with defects in homologous recombination repair. To identify more genetic defects that could sensitize tumors to CX-5461, we tested synthetic lethality for 480 DNA repair and genome maintenance genes to CX-5461, pyridostatin (PDS), a structurally unrelated G4-specific stabilizer, and BMH-21, which binds GC-rich DNA but not G4 structures. We identified multiple members of HRD, Fanconi Anemia pathways, and POLQ, a polymerase with a helicase domain important for G4 structure resolution. Significant synthetic lethality was observed with UBE2N and RNF168, key members of the DNA damage response associated ubiquitin signaling pathway. Loss-of-function of RNF168 and UBE2N resulted in significantly lower cell survival in the presence of CX-5461 and PDS but not BMH-21. RNF168 recruitment and histone ubiquitination increased with CX-5461 treatment, and nuclear ubiquitination response frequently co-localized with G4 structures. Pharmacological inhibition of UBE2N acted synergistically with CX-5461. In conclusion, we have uncovered novel genetic vulnerabilities to CX-5461 with potential significance for patient selection in future clinical trials.

Epiclomal: Probabilistic clustering of sparse single-cell DNA methylation data.

We present Epiclomal, a probabilistic clustering method arising from a hierarchical mixture model to simultaneously cluster sparse single-cell DNA methylation data and impute missing values. Using synthetic and published single-cell CpG datasets, we show that Epiclomal outperforms non-probabilistic methods and can handle the inherent missing data characteristic that dominates single-cell CpG genome sequences. Using newly generated single-cell 5mCpG sequencing data, we show that Epiclomal discovers sub-clonal methylation patterns in aneuploid tumour genomes, thus defining epiclones that can match or transcend copy number-determined clonal lineages and opening up an important form of clonal analysis in cancer. Epiclomal is written in R and Python and is available at https://github.com/shahcompbio/Epiclomal.

Clonal Decomposition and DNA Replication States Defined by Scaled Single-Cell Genome Sequencing.

Accurate measurement of clonal genotypes, mutational processes, and replication states from individual tumor-cell genomes will facilitate improved understanding of tumor evolution. We have developed DLP+, a scalable single-cell whole-genome sequencing platform implemented using commodity instruments, image-based object recognition, and open source computational methods. Using DLP+, we have generated a resource of 51,926 single-cell genomes and matched cell images from diverse cell types including cell lines, xenografts, and diagnostic samples with limited material. From this resource we have defined variation in mitotic mis-segregation rates across tissue types and genotypes. Analysis of matched genomic and image measurements revealed correlations between cellular morphology and genome ploidy states. Aggregation of cells sharing copy number profiles allowed for calculation of single-nucleotide resolution clonal genotypes and inference of clonal phylogenies and avoided the limitations of bulk deconvolution. Finally, joint analysis over the above features defined clone-specific chromosomal aneuploidy in polyclonal populations.

Divergent modes of clonal spread and intraperitoneal mixing in high-grade serous ovarian cancer.

We performed phylogenetic analysis of high-grade serous ovarian cancers (68 samples from seven patients), identifying constituent clones and quantifying their relative abundances at multiple intraperitoneal sites. Through whole-genome and single-nucleus sequencing, we identified evolutionary features including mutation loss, convergence of the structural genome and temporal activation of mutational processes that patterned clonal progression. We then determined the precise clonal mixtures comprising each tumor sample. The majority of sites were clonally pure or composed of clones from a single phylogenetic clade. However, each patient contained at least one site composed of polyphyletic clones. Five patients exhibited monoclonal and unidirectional seeding from the ovary to intraperitoneal sites, and two patients demonstrated polyclonal spread and reseeding. Our findings indicate that at least two distinct modes of intraperitoneal spread operate in clonal dissemination and highlight the distribution of migratory potential over clonal populations comprising high-grade serous ovarian cancers.

ATOH7 mutations cause autosomal recessive persistent hyperplasia of the primary vitreous.

The vertebrate basic helix-loop-helix (bHLH) transcription factor ATOH7 (Math5) is specifically expressed in the embryonic neural retina and is required for the genesis of retinal ganglion cells (RGCs) and optic nerves. In Atoh7 mutant mice, the absence of trophic factors secreted by RGCs prevents the development of the intrinsic retinal vasculature and the regression of fetal blood vessels, causing persistent hyperplasia of the primary vitreous (PHPV). We therefore screened patients with hereditary PHPV, as well as bilateral optic nerve aplasia (ONA) or hypoplasia (ONH), for mutations in ATOH7. We identified a homozygous ATOH7 mutation (N46H) in a large family with an autosomal recessive PHPV disease trait linked to 10q21, and a heterozygous variant (R65G, p.Arg65Gly) in one of five sporadic ONA patients. High-density single-nucleotide polymorphism analysis also revealed a CNTN4 duplication and an OTX2 deletion in the ONA cohort. Functional analysis of ATOH7 bHLH domain substitutions, by electrophoretic mobility shift and luciferase cotransfection assays, revealed that the N46H variant cannot bind DNA or activate transcription, consistent with structural modeling. The N46H variant also failed to rescue RGC development in mouse Atoh7-/- retinal explants. The R65G variant retains all of these activities, similar to wild-type human ATOH7. Our results strongly suggest that autosomal recessive persistent hyperplastic primary vitreous is caused by N46H and is etiologically related to nonsyndromic congenital retinal nonattachment. The R65G allele, however, cannot explain the ONA phenotype. Our study firmly establishes ATOH7 as a retinal disease gene and provides a functional basis to analyze new coding variants.

A hypomorphic Artemis human disease allele causes aberrant chromosomal rearrangements and tumorigenesis.

The Artemis gene encodes a DNA nuclease that plays important roles in non-homologous end-joining (NHEJ), a major double-strand break (DSB) repair pathway in mammalian cells. NHEJ factors repair general DSBs as well as programmed breaks generated during the lymphoid-specific DNA rearrangement, V(D)J recombination, which is required for lymphocyte development. Mutations that inactivate Artemis cause a human severe combined immunodeficiency syndrome associated with cellular radiosensitivity. In contrast, hypomorphic Artemis mutations result in combined immunodeficiency syndromes of varying severity, but, in addition, are hypothesized to predispose to lymphoid malignancy. To elucidate the distinct molecular defects caused by hypomorphic compared with inactivating Artemis mutations, we examined tumor predisposition in a mouse model harboring a targeted partial loss-of-function disease allele. We find that, in contrast to Artemis nullizygosity, the hypomorphic mutation leads to increased aberrant intra- and interchromosomal V(D)J joining events. We also observe that dysfunctional Artemis activity combined with p53 inactivation predominantly predisposes to thymic lymphomas harboring clonal translocations distinct from those observed in Artemis nullizygosity. Thus, the Artemis hypomorphic allele results in unique molecular defects, tumor spectrum and oncogenic chromosomal rearrangements. Our findings have significant implications for disease outcomes and treatment of patients with different Artemis mutations.

Leaky severe combined immunodeficiency and aberrant DNA rearrangements due to a hypomorphic RAG1 mutation.

The RAG1/2 endonuclease initiates programmed DNA rearrangements in progenitor lymphocytes by generating double-strand breaks at specific recombination signal sequences. This process, known as V(D)J recombination, assembles the vastly diverse antigen receptor genes from numerous V, D, and J coding segments. In vitro biochemical and cellular transfection studies suggest that RAG1/2 may also play postcleavage roles by forming complexes with the recombining ends to facilitate DNA end processing and ligation. In the current study, we examine the in vivo consequences of a mutant form of RAG1, RAG1-S723C, that is proficient for DNA cleavage, yet exhibits defects in postcleavage complex formation and end joining in vitro. We generated a knockin mouse model harboring the RAG1-S723C hypomorphic mutation and examined the immune system in this fully in vivo setting. RAG1-S723C homozygous mice exhibit impaired lymphocyte development and decreased V(D)J rearrangements. Distinct from RAG nullizygosity, the RAG1-S723C hypomorph results in aberrant DNA double-strand breaks within rearranging loci. RAG1-S723C also predisposes to thymic lymphomas associated with chromosomal translocations in a p53 mutant background, and heterozygosity for the mutant allele accelerates age-associated immune system dysfunction. Thus, our study provides in vivo evidence that implicates aberrant RAG1/2 activity in lymphoid tumor development and premature immunosenescence.