Dynamics of the most common pathogenic mtDNA variant m.3243A > G demonstrate frequency-dependency in blood and positive selection in the germline.

The A-to-G point mutation at position 3243 in the human mitochondrial genome (m.3243A > G) is the most common pathogenic mtDNA variant responsible for disease in humans. It is widely accepted that m.3243A > G levels decrease in blood with age, and an age correction representing ~ 2% annual decline is often applied to account for this change in mutation level. Here we report that recent data indicate that the dynamics of m.3243A > G are more complex and depend on the mutation level in blood in a bi-phasic way. Consequently, the traditional 2% correction, which is adequate 'on average', creates opposite predictive biases at high and low mutation levels. Unbiased age correction is needed to circumvent these drawbacks of the standard model. We propose to eliminate both biases by using an approach where age correction depends on mutation level in a biphasic way to account for the dynamics of m.3243A > G in blood. The utility of this approach was further tested in estimating germline selection of m.3243A > G. The biphasic approach permitted us to uncover patterns consistent with the possibility of positive selection for m.3243A > G. Germline selection of m.3243A > G shows an 'arching' profile by which selection is positive at intermediate mutant fractions and declines at high and low mutant fractions. We conclude that use of this biphasic approach will greatly improve the accuracy of modelling changes in mtDNA mutation frequencies in the germline and in somatic cells during aging.

Increased KRAS G12C Prevalence, High Tumor Mutational Burden, and Specific Mutational Signatures Are Associated With MUTYH Mutations: A Pan-Cancer Analysis.

The aim of this study was to determine the pan-cancer landscape of MUTYH alterations and the relationship between MUTYH mutations and potentially actionable biomarkers such as specific genomic alterations, tumor mutational burden, and mutational signatures. We used a large pan-cancer comprehensive genomic dataset from patients profiled (tissue next generation sequencing) during routine clinical care. Overall, 2.8% of 229 120 solid tumors had MUTYH alterations, of which 55% were predicted germline. Thirty tumor types had a 2% or greater MUTYH mutation rate. MUTYH-altered versus -WT cancers had significantly higher tumor mutational burden and more frequent alterations in KRAS G12C, but not in KRAS in general; these observations were statistically significant, especially in colorectal cancers. Across cancers, PD-L1 expression levels (immunohistochemistry) were not associated with MUTYH alteration status. In silico computation demonstrated that MUTYH mutational signatures are associated with higher levels of hydrophobicity (which may reflect higher immunogenicity of neoantigens) relative to several other signature types such as microsatellite instability. Survival of patients with MUTYH-altered versus -WT tumors was similar. In conclusion, comprehensive genomic profiling suggests that several features of MUTYH-altered cancers may be pharmacologically targetable. Drugs such as sotorasib (targeting KRAS G12C) and immune checkpoint inhibitors, targeting the increased mutational load and higher neo-antigen hydrophobicity/immunogenicity merit investigation in MUTYH-mutated malignancies.

A mitochondria-specific mutational signature of aging: increased rate of A > G substitutions on the heavy strand.

The mutational spectrum of the mitochondrial DNA (mtDNA) does not resemble any of the known mutational signatures of the nuclear genome and variation in mtDNA mutational spectra between different organisms is still incomprehensible. Since mitochondria are responsible for aerobic respiration, it is expected that mtDNA mutational spectrum is affected by oxidative damage. Assuming that oxidative damage increases with age, we analyse mtDNA mutagenesis of different species in regards to their generation length. Analysing, (i) dozens of thousands of somatic mtDNA mutations in samples of different ages (ii) 70053 polymorphic synonymous mtDNA substitutions reconstructed in 424 mammalian species with different generation lengths and (iii) synonymous nucleotide content of 650 complete mitochondrial genomes of mammalian species we observed that the frequency of AH > GH substitutions (H: heavy strand notation) is twice bigger in species with high versus low generation length making their mtDNA more AH poor and GH rich. Considering that AH > GH substitutions are also sensitive to the time spent single-stranded (TSSS) during asynchronous mtDNA replication we demonstrated that AH > GH substitution rate is a function of both species-specific generation length and position-specific TSSS. We propose that AH > GH is a mitochondria-specific signature of oxidative damage associated with both aging and TSSS.

A Novel HRD Signature Is Predictive of FOLFIRINOX Benefit in Metastatic Pancreatic Cancer.

BACKGROUND: Pancreatic cancer (PC) represents an aggressive disease with median overall survival (OS) of less than 1 year in the front-line setting. FOLFIRINOX and gemcitabine and paclitaxel (GP) are standard of care options for these patients; however, optimal selection of therapy is challenging. METHODS: Comprehensive genomic profiling was performed on 8358 PC patients. Outcomes were available for 1149 metastatic PC patients treated with 1L FOLFIRINOX or GP. A scar-based measure of HRD was called using a machine learning-based algorithm incorporating copy number and indel features. RESULTS: A scar-based HRD signature (HRDsig) was identified in 9% of patients. HRDsig significantly co-occurred with biallelic alterations in BRCA1/2, PALB2, BARD1, and RAD51C/D, but encompassed a larger population than that defined by BRCA1/BRCA2/PALB2 (9% vs. 6%). HRDsig was predictive of 1L FOLFIRNOX chemotherapy benefit with doubled OS relative to gemcitabine and paclitaxel (GP) (rwOS aHR 0.37 [0.22-0.62]), including 25% of the population with long-term (2 year+) survival in a real-world cohort of patients. Less benefit from FOLFIRINOX was observed in the HRDsig(-) population. Predictive value was seen in both the BRCA1/2/PALB2 mutant and wildtype populations, suggesting additional value to mutational profiling. CONCLUSION: A scar-based HRD biomarker was identified in a significant fraction of PC patients and is predictive of FOLFIRINOX benefit. Incorporating a biomarker like HRDsig could identify the right patients for platinum chemotherapy and potentially reduce FOLFIRINOX use by over 40%, minimizing toxicities with similar survival outcomes. Confirmatory studies should be performed.

Workflow Optimization for Identification of Female Germline or Oogonial Stem Cells in Human Ovarian Cortex Using Single-Cell RNA Sequence Analysis.

In 2004, the identification of female germline or oogonial stem cells (OSCs) that can support post-natal oogenesis in ovaries of adult mice sparked a major paradigm shift in reproductive biology. Although these findings have been independently verified, and further extended to include identification of OSCs in adult ovaries of many species ranging from pigs and cows to non-human primates and humans, a recent study rooted in single-cell RNA sequence analysis (scRNA-seq) of adult human ovarian cortical tissue claimed that OSCs do not exist, and that other groups working with OSCs following isolation by magnetic-assisted or fluorescence-activated cell sorting have mistaken perivascular cells (PVCs) for germ cells. Here we report that rare germ lineage cells with a gene expression profile matched to OSCs but distinct from that of other cells, including oocytes and PVCs, can be identified in adult human ovarian cortical tissue by scRNA-seq after optimization of analytical workflow parameters. Deeper cell-by-cell expression profiling also uncovered evidence of germ cells undergoing meiosis-I in adult human ovaries. Lastly, we show that, if not properly controlled for, PVCs can be inadvertently isolated during flow cytometry protocols designed to sort OSCs because of inherently high cellular autofluorescence. However, human PVCs and human germ cells segregate into distinct clusters following scRNA-seq due to non-overlapping gene expression profiles, which would preclude the mistaken identification and use of PVCs as OSCs during functional characterization studies.

3GOLD: optimized Levenshtein distance for clustering third-generation sequencing data.

BACKGROUND: Third-generation sequencing offers some advantages over next-generation sequencing predecessors, but with the caveat of harboring a much higher error rate. Clustering-related sequences is an essential task in modern biology. To accurately cluster sequences rich in errors, error type and frequency need to be accounted for. Levenshtein distance is a well-established mathematical algorithm for measuring the edit distance between words and can specifically weight insertions, deletions and substitutions. However, there are drawbacks to using Levenshtein distance in a biological context and hence has rarely been used for this purpose. We present novel modifications to the Levenshtein distance algorithm to optimize it for clustering error-rich biological sequencing data. RESULTS: We successfully introduced a bidirectional frameshift allowance with end-user determined accommodation caps combined with weighted error discrimination. Furthermore, our modifications dramatically improved the computational speed of Levenstein distance. For simulated ONT MinION and PacBio Sequel datasets, the average clustering sensitivity for 3GOLD was 41.45% (S.D. 10.39) higher than Sequence-Levenstein distance, 52.14% (S.D. 9.43) higher than Levenshtein distance, 55.93% (S.D. 8.67) higher than Starcode, 42.68% (S.D. 8.09) higher than CD-HIT-EST and 61.49% (S.D. 7.81) higher than DNACLUST. For biological ONT MinION data, 3GOLD clustering sensitivity was 27.99% higher than Sequence-Levenstein distance, 52.76% higher than Levenshtein distance, 56.39% higher than Starcode, 48% higher than CD-HIT-EST and 70.4% higher than DNACLUST. CONCLUSION: Our modifications to Levenshtein distance have improved its speed and accuracy compared to the classic Levenshtein distance, Sequence-Levenshtein distance and other commonly used clustering approaches on simulated and biological third-generation sequenced datasets. Our clustering approach is appropriate for datasets of unknown cluster centroids, such as those generated with unique molecular identifiers as well as known centroids such as barcoded datasets. A strength of our approach is high accuracy in resolving small clusters and mitigating the number of singletons.

Genomic landscape of pleural and peritoneal mesothelioma tumours.

BACKGROUND: Malignant pleural and peritoneal mesotheliomas are rare malignancies with unacceptable poor prognoses and limited treatment options. The genomic landscape is mainly characterised by the loss of tumour suppressor genes and mutations in DNA repair genes. Currently, data from next-generation sequencing (NGS) of mesothelioma tumours is restricted to a limited number of cases; moreover, data comparing molecular features of mesothelioma from the pleural and peritoneal origin with NGS are lacking. METHODS: We analysed 1113 pleural mesothelioma and 355 peritoneal mesothelioma samples. All tumours were sequenced with the FoundationOne® or FoundationOne®CDx assay for detection of substitutions, insertion-deletions, copy-number alterations and selected rearrangements in at least 324 cancer genes. RESULTS: This analysis revealed alterations in 19 genes with an overall prevalence of at least 2%. Alterations in BAP1, CDKN2A, CDKN2B, NF2, MTAP, TP53 and SETD2 occurred with a prevalence of at least 10%. Peritoneal, compared to pleural mesothelioma, was characterised by a lower prevalence of alterations in CDKN2A, CDKN2B and MTAP. Moreover, we could define four distinct subgroups according to alterations in BAP1 and CDKN2A/B. Alterations in Hedgehog pathway-related genes (PTCH1/2 and SUFU) and Hippo pathway-related gene (NF2) as well as KRAS, EGFR, PDGFRA/B, ERBB2 and FGFR3 were detected in both cohorts. CONCLUSION: Here, we report the molecular aberrations from the largest cohort of patients with mesothelioma. This analysis identified a proportion of patients with targetable alterations and suggests that molecular profiling can identify new treatment options for patients with mesothelioma.

Reanalysis of mtDNA mutations of human primordial germ cells (PGCs) reveals NUMT contamination and suggests that selection in PGCs may be positive.

The resilience of the mitochondrial genome (mtDNA) to a high mutational pressure depends, in part, on negative purifying selection in the germline. A paradigm in the field has been that such selection, at least in part, takes place in primordial germ cells (PGCs). Specifically, Floros et al. (Nature Cell Biology 20: 144-51) reported an increase in the synonymity of mtDNA mutations (a sign of purifying selection) between early-stage and late-stage PGCs. We re-analyzed Floros' et al. data and determined that their mutational dataset was significantly contaminated with single nucleotide variants (SNVs) derived from a nuclear sequence of mtDNA origin (NUMT) located on chromosome 5. Contamination was caused by co-amplification of the NUMT sequence by cross-specific PCR primers. Importantly, when we removed NUMT-derived SNVs, the evidence of purifying selection was abolished. In addition to bulk PGCs, Floros et al. reported the analysis of single-cell late-stage PGCs, which were amplified with different sets of PCR primers that cannot amplify the NUMT sequence. Accordingly, there were no NUMT-derived SNVs among single PGC mutations. Interestingly, single PGC mutations show adecreaseof synonymity with increased intracellular mutant fraction. More specifically, nonsynonymous mutations show faster intracellular genetic drift towards higher mutant fraction than synonymous ones. This pattern is incompatible with predominantly negative selection. This suggests that germline selection of mtDNA mutations is a complex phenomenon and that the part of this process that takes place in PGCs may be predominantly positive. However counterintuitive, positive germline selection of detrimental mtDNA mutations has been reported previously andpotentially may be evolutionarily advantageous.

Pan-Cancer Analysis of Copy-Number Features Identifies Recurrent Signatures and a Homologous Recombination Deficiency Biomarker to Predict Poly (ADP-Ribose) Polymerase Inhibitor Response.

PURPOSE: Copy-number (CN) features reveal the molecular state of cancers and may have predictive and prognostic value in the treatment of cancer. We sought to apply published CN analysis methods to a large pan-cancer data set and characterize ubiquitous CN signatures across tumor types, including potential utility for treatment selection. METHODS: We analyzed the landscape of CN features in 260,333 pan-cancer samples. We examined the association of 10 signatures with genomic alterations and clinical characteristics and trained a machine learning classifier using CN and insertion and deletion features to detect homologous recombination deficiency signature (HRDsig) positivity. Clinical outcomes were assessed using a real-world clinicogenomic database (CGDB) of comprehensive genomic profiling linked to deidentified, electronic health record-derived clinical data. RESULTS: CN signatures were prevalent across cancer types and associated with diverse processes including focal tandem duplications, seismic amplifications, genome-wide loss of heterozygosity (gLOH), and HRD. Our novel HRDsig outperformed gLOH in predicting BRCAness and effectively distinguished biallelic BRCA and homologous recombination-repair wild-type (HRRwt) samples pan-tumor, demonstrating high sensitivity to detect biallelic BRCA in ovarian (93%) and other HRD-associated cancers (80%-87%). Pan-tumor prevalence of HRDsig was 6.4%. HRRwt cases represented a significant fraction of the HRDsig-positive cohort, likely reflecting a population with nongenomic mechanisms of HRD. In ovarian and prostate CGDBs, HRDsig identified more patients than gLOH and had predictive value for poly (ADP-ribose) polymerase inhibitor (PARPi) benefit. CONCLUSION: Tumor CN profiles are informative, revealing diverse processes active in cancer. We describe the landscape of 10 CN signatures in a large pan-cancer cohort, including two associated with HRD. We trained a machine learning-based HRDsig that robustly identified BRCAness and associated with biallelic BRCA pan-tumor, and was predictive of PARPi benefit in real-world ovarian and prostate data sets.

Co-occurring BRCA2/SPOP Mutations Predict Exceptional Poly (ADP-ribose) Polymerase Inhibitor Sensitivity in Metastatic Castration-Resistant Prostate Cancer.

BACKGROUND AND OBJECTIVE: BRCA2 mutations in metastatic castration-resistant prostate cancer (mCRPC) confer sensitivity to poly (ADP-ribose) polymerase (PARP) inhibitors. However, additional factors predicting PARP inhibitor efficacy in mCRPC are needed. Preclinical studies support a relationship between speckle-type POZ protein (SPOP) inactivation and PARP inhibitor sensitivity. We hypothesized that SPOP mutations may predict enhanced PARP inhibitor response in BRCA2-altered mCRPC. METHODS: We conducted a multicenter retrospective study involving 13 sites. We identified 131 patients with BRCA2-altered mCRPC treated with PARP inhibitors, 14 of which also carried concurrent SPOP mutations. The primary efficacy endpoint was prostate-specific antigen (PSA) response rate (≥50% PSA decline). The secondary endpoints were biochemical progression-free survival (PSA-PFS), clinical/radiographic progression-free survival (PFS), and overall survival (OS). These were compared by multivariable Cox proportional hazard models adjusting for age, tumor stage, baseline PSA level, Gleason sum, prior therapies, BRCA2 alteration types, and co-occurring mutations. KEY FINDINGS AND LIMITATIONS: Baseline characteristics were similar between groups. PSA responses were observed in 60% (70/117) of patients with BRCA2mut/SPOPwt disease and in 86% (12/14) of patients with BRCA2mut/SPOPmut disease (p = 0.06). The median time on PARP inhibitor treatment was 24.0 mo (95% confidence interval [CI] 19.2 mo to not reached) in this group versus 8.0 mo (95% CI 6.1-10.9 mo) in patients with BRCA2 mutation alone (p = 0.05). In an unadjusted analysis, patients with BRCA2mut/SPOPmut disease experienced longer PSA-PFS (hazard ratio [HR] 0.33 [95% CI 0.15-0.72], p = 0.005) and clinical/radiographic PFS (HR 0.4 [95% CI 0.18-0.86], p = 0.02), and numerically longer OS (HR 0.4 [95% CI 0.15-1.12], p = 0.08). In a multivariable analysis including histology, Gleason sum, prior taxane, prior androgen receptor pathway inhibitor, stage, PSA, BRCA2 alteration characteristics, and other co-mutations, patients with BRCA2mut/SPOPmut disease experienced longer PSA-PFS (HR 0.16 [95% CI 0.05-0.47], adjusted p = 0.001), clinical/radiographic PFS (HR 0.28 [95% CI 0.1-0.81], adjusted p = 0.019), and OS (HR 0.19 [95% CI 0.05-0.69], adjusted p = 0.012). In a separate cohort of patients not treated with a PARP inhibitor, there was no difference in OS between patients with BRCA2mut/SPOPmut versus BRCA2mut/SPOPwt disease (HR 0.97 [95% CI 0.40-2.4], p = 0.94). In a genomic signature analysis, Catalog of Somatic Mutations in Cancer (COSMIC) SBS3 scores predictive of homologous recombination repair (HRR) defects were higher for BRCA2mut/SPOPmut than for BRCA2mut/SPOPwt disease (p = 0.04). This was a retrospective study, and additional prospective validation cohorts are needed. CONCLUSIONS AND CLINICAL IMPLICATIONS: In this retrospective analysis, PARP inhibitors appeared more effective in patients with BRCA2mut/SPOPmut than in patients with BRCA2mut/SPOPwt mCRPC. This may be related to an increase in HRR defects in coaltered disease. PATIENT SUMMARY: In this study, we demonstrate that co-alteration of both BRCA2 and SPOP predicts superior clinical outcomes to treatment with poly (ADP-ribose) polymerase (PARP) inhibitors than BRCA2 alteration without SPOP mutation.

Mitochondrial Pseudogenes Suggest Repeated Inter-Species Hybridization among Direct Human Ancestors.

The hypothesis that the evolution of humans involves hybridization between diverged species has been actively debated in recent years. We present the following novel evidence in support of this hypothesis: the analysis of nuclear pseudogenes of mtDNA ("NUMTs"). NUMTs are considered "mtDNA fossils" as they preserve sequences of ancient mtDNA and thus carry unique information about ancestral populations. Our comparison of a NUMT sequence shared by humans, chimpanzees, and gorillas with their mtDNAs implies that, around the time of divergence between humans and chimpanzees, our evolutionary history involved the interbreeding of individuals whose mtDNA had diverged as much as ~4.5 Myr prior. This large divergence suggests a distant interspecies hybridization. Additionally, analysis of two other NUMTs suggests that such events occur repeatedly. Our findings suggest a complex pattern of speciation in primate/human ancestors and provide one potential explanation for the mosaic nature of fossil morphology found at the emergence of the hominin lineage. A preliminary version of this manuscript was uploaded to the preprint server BioRxiv in 2017 (10.1101/134502).

LUCS: a high-resolution nucleic acid sequencing tool for accurate long-read analysis of individual DNA molecules.

Nucleic acid sequence analyses are fundamental to all aspects of biological research, spanning aging, mitochondrial DNA (mtDNA) and cancer, as well as microbial and viral evolution. Over the past several years, significant improvements in DNA sequencing, including consensus sequence analysis, have proven invaluable for high-throughput studies. However, all current DNA sequencing platforms have limited utility for studies of complex mixtures or of individual long molecules, the latter of which is crucial to understanding evolution and consequences of single nucleotide variants and their combinations. Here we report a new technology termed LUCS (Long-molecule UMI-driven Consensus Sequencing), in which reads from third-generation sequencing are aggregated by unique molecular identifiers (UMIs) specific for each individual DNA molecule. This enables in-silico reconstruction of highly accurate consensus reads of each DNA molecule independent of other molecules in the sample. Additionally, use of two UMIs enables detection of artificial recombinants (chimeras). As proof of concept, we show that application of LUCS to assessment of mitochondrial genomes in complex mixtures from single cells was associated with an error rate of 1X10-4 errors/nucleotide. Thus, LUCS represents a major step forward in DNA sequencing that offers high-throughput capacity and high-accuracy reads in studies of long DNA templates and nucleotide variants in heterogenous samples.