Whole-genome Mutational Analysis for Tumor-informed Detection of Circulating Tumor DNA in Patients with Urothelial Carcinoma.

BACKGROUND AND OBJECTIVE: Circulating tumor DNA (ctDNA) can be used for sensitive detection of minimal residual disease (MRD). However, the probability of detecting ctDNA in settings of low tumor burden is limited by the number of mutations analyzed and the plasma volume available. We used a whole-genome sequencing (WGS) approach for ctDNA detection in patients with urothelial carcinoma. METHODS: We used a tumor-informed WGS approach for ctDNA-based detection of MRD and evaluation of treatment responses. We analyzed 916 longitudinally collected plasma samples from 112 patients with localized muscle-invasive bladder cancer who received neoadjuvant chemotherapy (NAC) before radical cystectomy. Recurrence-free survival (primary endpoint), overall survival, and ctDNA dynamics during NAC were assessed. KEY FINDINGS AND LIMITATIONS: We found that WGS-based ctDNA detection is prognostic for patient outcomes with a median lead time of 131 d over radiographic imaging. WGS-based ctDNA assessment after radical cystectomy identified recurrence with sensitivity of 91% and specificity of 92%. In addition, genomic characterization of post-treatment plasma samples with a high ctDNA level revealed acquisition of platinum therapy-associated mutational signatures and copy number variations not present in the primary tumors. The sequencing depth is a limitation for studying tumor evolution. CONCLUSIONS AND CLINICAL IMPLICATIONS: Our results support the use of WGS for ultrasensitive ctDNA detection and highlight the possibility of plasma-based tracking of tumor evolution. WGS-based ctDNA detection represents a promising option for clinical use owing to the low volume of plasma needed and the ease of performing WGS, eliminating the need for personalized assay design.

Host Nuclear Genome Copy Number Variations Identify High-Risk Anal Precancers in People Living With HIV.

BACKGROUND: People living with HIV (PLWH) have substantially increased incidence of anal precancer and cancer. There are very little data regarding genomic disturbances in anal precancers among PLWH. In this study, specific chromosomal variants were identified in anal squamous intraepithelial lesions. METHODS: Overall, 63 anal biopsy specimens (27 low-grade intraepithelial lesions [LSIL] and 36 high-grade intraepithelial lesions [HSIL]) were collected from PLWH obtained as part of anal cancer screening in our NYC-based health system. Data on patient demographics, anal cytological, and high-risk human papillomavirus (HR-HPV) diagnoses were collected. Specimens were tested for a panel of chromosomal alterations associated with HPV-induced oncogenesis using fluorescence in situ hybridization, and analyses compared the associations of these alterations with clinical characteristics. RESULTS: Gains of 3q26, 5p15, 20q13, and cen7 were detected in 42%, 31%, 31%, and 19% of HSIL compared with 7%, 0%, 4%, and 0% of LSIL, respectively. If at least 1 abnormality was observed, 89% had a 3q26 gain. In lesions with 5p15 gains, 20q13 gains co-occurred in 91% of cases, while cen7 gain only co-occurred with the other 3 alterations. The sensitivity and specificity of any alteration to predict HSIL were 47% (95% CI: 30%-65%) and 93% (95% CI: 76%-99%), respectively. CONCLUSIONS: Genomic alterations seen in HPV-associated cancers may help distinguish anal LSIL from HSIL. 3q26 amplification may be an early component of anal carcinogenesis, preceding 5p16, 20q13, and/or chr7. IMPACT: Insights into potential genomic biomarkers for discriminating high-risk anal precancers are shared.

Whole-exome sequencing of Nigerian benign prostatic hyperplasia reveals increased alterations in apoptotic pathways.

BACKGROUND: Through whole-exome sequencing of 60 formalin-fixed paraffin-embedded Nigerian (NGRn) benign prostatic hyperplasia (BPH) samples, we identified germline and somatic alterations in apoptotic pathways impacting BPH development and progression. Prostate enlargement is a common occurrence in male aging; however, this enlargement can lead to lower urinary tract symptoms that negatively impact quality of life. This impact is disproportionately present in men of African ancestry. BPH pathophysiology is poorly understood and studies examining non-European populations are lacking. METHODS: In this study, NGRn BPH, normal prostate, and prostate cancer (PCa) tumor samples were sequenced and compared to characterize genetic alterations in NGRn BPH. RESULTS: Two hundred and two nonbenign, ClinVar-annotated germline variants were present in NGRn BPH samples. Six genes [BRCA1 (92%), HSD3B1 (85%), TP53 (37%), PMS2 (23%), BARD1 (20%), and BRCA2 (17%)] were altered in at least 10% of samples; however, compared to NGRn normal and tumor, the frequency of alterations in BPH samples showed no significant differences at the gene or variant level. BRCA2_rs11571831 and TP53_rs1042522 germline alterations had a statistically significant co-occurrence interaction in BPH samples. In at least two BPH samples, 173 genes harbored somatic variants known to be clinically actionable. Three genes (COL18A1, KIF16B, and LRP1) showed a statistically significant (p < 0.05) higher frequency in BPH. NGRn BPH also had five gene pairs (PKD1/KIAA0100, PKHD1/PKD1, DNAH9/LRP1B, NWD1/DCHS2, and TCERG1/LMTK2) with statistically significant co-occurring interactions. Two hundred and seventy-nine genes contained novel somatic variants in NGRn BPH. Three genes (CABP1, FKBP1C, and RP11-595B24.2) had a statistically significant (p < 0.05) higher alteration frequency in NGRn BPH and three were significantly higher in NGRn tumor (CACNA1A, DMKN, and CACNA2D2). Pairwise Fisher's exact tests showed 14 gene pairs with statistically significant (p < 0.05) interactions and four interactions approaching significance (p < 0.10). Mutational patterns in NGRn BPH were similar to COSMIC (Catalog of Somatic Mutations in Cancer) signatures associated with aging and dysfunctional DNA damage repair. CONCLUSIONS: NGRn BPH contained significant germline alteration interactions (BRCA2_rs11571831 and TP53_rs1042522) and increased somatic alteration frequencies (LMTK2, LRP1, COL18A1, CABP1, and FKBP1C) that impact apoptosis. Normal prostate development is maintained by balancing apoptotic and proliferative activity. Dysfunction in either mechanism can lead to abnormal prostate growth. This work is the first to examine genomic sequencing in NGRn BPH and provides data that fill known gaps in the understanding BPH and how it impacts men of African ancestry.

Exercise-Induced Cell-Free DNA Correlates with Energy Expenditure in Multiple Exercise Protocols.

PURPOSE: Exercise-induced cell-free DNA (ei-cfDNA) has been studied in response to various types of exercise. Its correlation with exercise intensity and duration has been observed consistently. However, comprehensive measurements and exploration of the tissue of origin are lacking. The aim of this study is to establish precise connections between exercise variables and the distribution of tissue of origin, aiming to provide further evidence supporting its use as a biomarker for exercise. METHODS: Twelve self-identified active adults (six men and six women) performed a crossover study starting with either endurance testing or resistance testing under different intensities and protocols. We obtained blood before and after each exercise session and measured the levels of cfDNA and determined its tissue of origin utilizing cell type-specific DNA methylation patterns in plasma. RESULTS: We found that when duration and intensity are fixed, ei-cfDNA fold change correlates with energy expenditure ( P = 0.001) in endurance testing and years trained ( P = 0.001) in resistance testing. Most of the ei-cfDNA comes from increases in white blood cells (~95%) where neutrophils make up the majority (~74%) and the distribution is different between exercise modalities and protocols. CONCLUSIONS: This study highlights the potential of exercise-induced cfDNA as a biomarker for exercise, showing correlations with energy expenditure and a consistent pattern of tissue origin. Additional research is needed to investigate potential sex differences in the response of cfDNA to exercise, further exploring its clinical implications.

Elevated cfDNA after exercise is derived primarily from mature polymorphonuclear neutrophils, with a minor contribution of cardiomyocytes.

Strenuous physical exercise causes a massive elevation in the concentration of circulating cell-free DNA (cfDNA), which correlates with effort intensity and duration. The cellular sources and physiological drivers of this phenomenon are unknown. Using methylation patterns of cfDNA and associated histones, we show that cfDNA in exercise originates mostly in extramedullary polymorphonuclear neutrophils. Strikingly, cardiomyocyte cfDNA concentration increases after a marathon, consistent with elevated troponin levels and indicating low-level, delayed cardiac cell death. Physical impact, low oxygen levels, and elevated core body temperature contribute to neutrophil cfDNA release, while muscle contraction, increased heart rate, ß-adrenergic signaling, or steroid treatment fail to cause elevation of cfDNA. Physical training reduces neutrophil cfDNA release after a standard exercise, revealing an inverse relationship between exercise-induced cfDNA release and training level. We speculate that the release of cfDNA from neutrophils in exercise relates to the activation of neutrophils in the context of exercise-induced muscle damage.

A framework for clinical cancer subtyping from nucleosome profiling of cell-free DNA.

Cell-free DNA (cfDNA) has the potential to inform tumor subtype classification and help guide clinical precision oncology. Here we develop Griffin, a framework for profiling nucleosome protection and accessibility from cfDNA to study the phenotype of tumors using as low as 0.1x coverage whole genome sequencing data. Griffin employs a GC correction procedure tailored to variable cfDNA fragment sizes, which generates a better representation of chromatin accessibility and improves the accuracy of cancer detection and tumor subtype classification. We demonstrate estrogen receptor subtyping from cfDNA in metastatic breast cancer. We predict estrogen receptor subtype in 139 patients with at least 5% detectable circulating tumor DNA with an area under the receive operator characteristic curve (AUC) of 0.89 and validate performance in independent cohorts (AUC = 0.96). In summary, Griffin is a framework for accurate tumor subtyping and can be generalizable to other cancer types for precision oncology applications.

Characterization of genetic predisposition to molecular subtypes of breast cancer in Brazilian patients.

Introduction: BRCA1 and BRCA2 germline pathogenic variants (GPVs) account for most of the 5-10% of breast cancer (BC) that is attributable to inherited genetic variants. BRCA1 GPVs are associated with the triple negative subtype, whereas BRCA2 GPVs are likely to result in higher grade, estrogen-receptor positive BCs. The contribution of other genes of high and moderate risk for BC has not been well defined and risk estimates to specific BC subtypes is lacking, especially for an admixed population like Brazilian. Objective: The aim of this study is to evaluate the value of a multigene panel in detecting germline mutations in cancer-predisposing genes for Brazilian BC patients and its relation with molecular subtypes and the predominant molecular ancestry. Patients and methods: A total of 321 unrelated BC patients who fulfilled NCCN criteria for BRCA1/2 testing between 2016-2018 were investigated with a 94-genes panel. Molecular subtypes were retrieved from medical records and ancestry-specific variants were obtained from off-target reads obtained from the sequencing data. Results: We detected 83 GPVs in 81 patients (positivity rate of 25.2%). Among GPVs, 47% (39/83) were identified in high-risk BC genes (BRCA1/2, PALB2 and TP53) and 18% (15/83) in moderate-penetrance genes (ATM, CHEK2 and RAD51C). The remainder of the GPVs (35% - 29/83), were identified in lower-risk genes. As for the molecular subtypes, triple negative BC had a mutation frequency of 31.6% (25/79), with predominance in BRCA1 (12.6%; 10/79). Among the luminal subtypes, except Luminal B HER2-positive, 18.7% (29/155) had GPV with BRCA1/2 genes contributing 7.1% (11/155) and non-BRCA1/2 genes, 12.9% (20/155). For Luminal B HER2-positive subtype, 40% (16/40) had GPVs, with a predominance of ATM gene (15% - 6/40) and BRCA2 with only 2.5% (1/40). Finally, HER2-enriched subtype presented a mutation rate of 30.8% (4/13) with contribution of BRCA2 of 7.5% (1/13) and non-BRCA1/2 of 23% (3/13). Variants of uncertain significance (VUS) were identified in 77.6% (249/321) of the patients and the number of VUS was increased in patients with Asian and Native American ancestry. Conclusion: The multigene panel contributed to identify GPVs in genes other than BRCA1/2, increasing the positivity of the genetic test from 9.6% (BRCA1/2) to 25.2% and, considering only the most clinically relevant BC predisposing genes, to 16.2%. These results indicate that women with clinical criteria for hereditary BC may benefit from a multigene panel testing, as it allows identifying GPVs in genes that directly impact the clinical management of these patients and family members.

Diverse mutational landscapes in human lymphocytes.

The lymphocyte genome is prone to many threats, including programmed mutation during differentiation1, antigen-driven proliferation and residency in diverse microenvironments. Here, after developing protocols for expansion of single-cell lymphocyte cultures, we sequenced whole genomes from 717 normal naive and memory B and T cells and haematopoietic stem cells. All lymphocyte subsets carried more point mutations and structural variants than haematopoietic stem cells, with higher burdens in memory cells than in naive cells, and with T cells accumulating mutations at a higher rate throughout life. Off-target effects of immunological diversification accounted for approximately half of the additional differentiation-associated mutations in lymphocytes. Memory B cells acquired, on average, 18 off-target mutations genome-wide for every on-target IGHV mutation during the germinal centre reaction. Structural variation was 16-fold higher in lymphocytes than in stem cells, with around 15% of deletions being attributable to off-target recombinase-activating gene activity. DNA damage from ultraviolet light exposure and other sporadic mutational processes generated hundreds to thousands of mutations in some memory cells. The mutation burden and signatures of normal B cells were broadly similar to those seen in many B-cell cancers, suggesting that malignant transformation of lymphocytes arises from the same mutational processes that are active across normal ontogeny. The mutational landscape of normal lymphocytes chronicles the off-target effects of programmed genome engineering during immunological diversification and the consequences of differentiation, proliferation and residency in diverse microenvironments.

The Great Majority of Homologous Recombination Repair-Deficient Tumors Are Accounted for by Established Causes.

Background: Gene-agnostic genomic biomarkers were recently developed to identify homologous recombination deficiency (HRD) tumors that are likely to respond to treatment with PARP inhibitors. Two machine-learning algorithms that predict HRD status, CHORD, and HRDetect, utilize various HRD-associated features extracted from whole-genome sequencing (WGS) data and show high sensitivity in detecting patients with BRCA1/2 bi-allelic inactivation in all cancer types. When using only DNA mutation data for the detection of potential causes of HRD, both HRDetect and CHORD find that 30-40% of cases that have been classified as HRD are due to unknown causes. Here, we examined the impact of tumor-specific thresholds and measurement of promoter methylation of BRCA1 and RAD51C on unexplained proportions of HRD cases across various tumor types. Methods: We gathered published CHORD and HRDetect probability scores for 828 samples from breast, ovarian, and pancreatic cancer from previous studies, as well as evidence of their biallelic inactivation (by either DNA alterations or promoter methylation) in HR-related genes. ROC curve analysis evaluated the performance of each classifier in specific cancer. Tenfold nested cross-validation was used to find the optimal threshold values of HRDetect and CHORD for classifying HR-deficient samples within each cancer type. Results: With the universal threshold, HRDetect has higher sensitivity in the detection of biallelic inactivation in BRCA1/2 than CHORD and resulted in a higher proportion of unexplained cases. When promoter methylation was excluded, in ovarian carcinoma, the proportion of unexplained cases increased from 26.8 to 48.8% for HRDetect and from 14.7 to 41.2% for CHORD. A similar increase was observed in breast cancer. Applying cancer-type-specific thresholds led to similar sensitivity and specificity for both methods. The cancer-type-specific thresholds for HRDetect reduced the number of unexplained cases from 21 to 12.3% without reducing the 96% sensitivity to known events. For CHORD, unexplained cases were reduced from 10 to 9% while sensitivity increased from 85.3 to 93.9%. Conclusion: These results suggest that WGS-based HRD classifiers should be adjusted for tumor types. When applied, only ∼10% of breast, ovarian, and pancreas cancer cases are not explained by known events in our dataset.

Machine-learning of complex evolutionary signals improves classification of SNVs.

Conservation is a strong predictor for the pathogenicity of single-nucleotide variants (SNVs). However, some positions that present complex conservation patterns across vertebrates stray from this paradigm. Here, we analyzed the association between complex conservation patterns and the pathogenicity of SNVs in the 115 disease-genes that had sufficient variant data. We show that conservation is not a one-rule-fits-all solution since its accuracy highly depends on the analyzed set of species and genes. For example, pairwise comparisons between the human and 99 vertebrate species showed that species differ in their ability to predict the clinical outcomes of variants among different genes using conservation. Furthermore, certain genes were less amenable for conservation-based variant prediction, while others demonstrated species that optimize prediction. These insights led to developing EvoDiagnostics, which uses the conservation against each species as a feature within a random-forest machine-learning classification algorithm. EvoDiagnostics outperformed traditional conservation algorithms, deep-learning based methods and most ensemble tools in every prediction-task, highlighting the strength of optimizing conservation analysis per-species and per-gene. Overall, we suggest a new and a more biologically relevant approach for analyzing conservation, which improves prediction of variant pathogenicity.

Cell-free DNA profiling informs all major complications of hematopoietic cell transplantation.

Allogeneic hematopoietic cell transplantation (HCT) provides effective treatment for hematologic malignancies and immune disorders. Monitoring of posttransplant complications is critical, yet current diagnostic options are limited. Here, we show that cell-free DNA (cfDNA) in blood is a versatile analyte for monitoring of the most important complications that occur after HCT: graft-versus-host disease (GVHD), a frequent immune complication of HCT, infection, relapse of underlying disease, and graft failure. We demonstrate that these therapeutic complications are informed from a single assay, low-coverage bisulfite sequencing of cfDNA, followed by disease-specific bioinformatic analyses. To inform GVHD, we profile cfDNA methylation marks to trace the cfDNA tissues-of-origin and to quantify tissue-specific injury. To inform infection, we implement metagenomic cfDNA profiling. To inform cancer relapse, we implement analyses of tumor-specific genomic aberrations. Finally, to detect graft failure, we quantify the proportion of donor- and recipient-specific cfDNA. We applied this assay to 170 plasma samples collected from 27 HCT recipients at predetermined timepoints before and after allogeneic HCT. We found that the abundance of solid-organ-derived cfDNA in the blood at 1 mo after HCT is predictive of acute GVHD (area under the curve, 0.88). Metagenomic profiling of cfDNA revealed the frequent occurrence of viral reactivation in this patient population. The fraction of donor-specific cfDNA was indicative of relapse and remission, and the fraction of tumor-specific cfDNA was informative of cancer relapse. This proof-of-principle study shows that cfDNA has the potential to improve the care of allogeneic HCT recipients by enabling earlier detection and better prediction of the complex array of complications that occur after HCT.

COOBoostR: An Extreme Gradient Boosting-Based Tool for Robust Tissue or Cell-of-Origin Prediction of Tumors.

We present here COOBoostR, a computational method designed for the putative prediction of the tissue- or cell-of-origin of various cancer types. COOBoostR leverages regional somatic mutation density information and chromatin mark features to be applied to an extreme gradient boosting-based machine-learning algorithm. COOBoostR ranks chromatin marks from various tissue and cell types, which best explain the somatic mutation density landscape of any sample of interest. A specific tissue or cell type matching the chromatin mark feature with highest explanatory power is designated as a potential tissue- or cell-of-origin. Through integrating either ChIP-seq based chromatin data, along with regional somatic mutation density data derived from normal cells/tissue, precancerous lesions, and cancer types, we show that COOBoostR outperforms existing random forest-based methods in prediction speed, with comparable or better tissue or cell-of-origin prediction performance (prediction accuracy-normal cells/tissue: 76.99%, precancerous lesions: 95.65%, cancer cells: 89.39%). In addition, our results suggest a dynamic somatic mutation accumulation at the normal tissue or cell stage which could be intertwined with the changes in open chromatin marks and enhancer sites. These results further represent chromatin marks shaping the somatic mutation landscape at the early stage of mutation accumulation, possibly even before the initiation of precancerous lesions or neoplasia.

Whole-exome Sequencing of Nigerian Prostate Tumors from the Prostate Cancer Transatlantic Consortium (CaPTC) Reveals DNA Repair Genes Associated with African Ancestry.

In this study, we used whole-exome sequencing of a cohort of 45 advanced-stage, treatment-naïve Nigerian (NG) primary prostate cancer tumors and 11 unmatched nontumor tissues to compare genomic mutations with African American (AA) and European American (EA) The Cancer Genome Atlas (TCGA) prostate cancer. NG samples were collected from six sites in central and southwest Nigeria. After whole-exome sequencing, samples were processed using GATK best practices. BRCA1 (100%), BARD1 (45%), BRCA2 (27%), and PMS2(18%) had germline alterations in at least two NG nontumor samples. Across 111 germline variants, the AA cohort reflected a pattern [BRCA1 (68%), BARD1 (34%), BRCA2 (28%), and PMS2 (16%)] similar to NG samples. Of the most frequently mutated genes, BRCA1 showed a statistically (P ≤ 0.05) higher germline mutation frequency in men of African ancestry (MAA) and increasing variant frequency with increased African ancestry. Disaggregating gene-level mutation frequencies by variants revealed both ancestry-linked and NG-specific germline variant patterns. Driven by rs799917 (T>C), BRCA1 showed an increasing mutation frequency as African ancestry increased. BRCA2_rs11571831 was present only in MAA, and BRCA2_rs766173 was elevated in NG men. A total of 133 somatic variants were present in 26 prostate cancer-associated genes within the NG tumor cohort. BRCA2 (27%), APC (20%), ATM (20%), BRCA1 (13%), DNAJC6 (13%), EGFR (13%), MAD1L1 (13%), MLH1 (11%), and PMS2 (11%) showed mutation frequencies >10%. Compared with TCGA cohorts, NG tumors showed statistically significant elevated frequencies of BRCA2, APC, and BRCA1. The NG cohort variant pattern shared similarities (cosign similarities ≥0.734) with Catalogue of Somatic Mutations in Cancer signatures 5 and 6, and mutated genes showed significant (q < 0.001) gene ontology (GO) and functional enrichment in mismatch repair and non-homologous repair deficiency pathways. Here, we showed that mutations in DNA damage response genes were higher in NG prostate cancer samples and that a portion of those mutations correlate with African ancestry. Moreover, we identified variants of unknown significance that may contribute to population-specific routes of tumorigenesis and treatment. These results present the most comprehensive characterization of the NG prostate cancer exome to date and highlight the need to increase diversity of study populations. Significance: MAA have higher rates of prostate cancer incidence and mortality, however, are severely underrepresented in genomic studies. This is the first study utilizing whole-exome sequencing in NG men to identify West African ancestry-linked variant patterns that impact DNA damage repair pathways.