Whole-exome sequencing of ovarian cancer families uncovers putative predisposition genes

A prognostic model for ovarian neoplasms established by an integrated analysis of 1580 transcriptomic profiles

Even after debulking surgery combined with chemotherapy or new adjuvant chemotherapy paired with internal surgery, the average year of disease free survival in advanced ovarian cancer was approximately 1.7 years1. The development of a molecular predictor of early recurrence would allow for the identification of ovarian cancer (OC) patients with high risk of relapse. The Ovarian Cancer Disease Free Survival Predictor (ODFSP), a predictive model constructed from a special set of 1580 OC tumors in which gene expression was assessed using both microarray and sequencing platforms, was created by our team. To construct gene expression barcodes that were resistant to biases caused by disparate profiling platforms and batch effects, we employed a meta-analysis methodology that was based on the binary gene pair technique. We demonstrate that ODFSP is a reliable single-sample predictor of early recurrence (1 year or less) using the largest pool of OC transcriptome data sets available to date. The ODFSP model showed significantly high prognostic value for binary recurrence prediction unaffected by clinicopathologic factors, with a meta-estimate of the area under the receiver operating curve of 0.64 (P  =  4.6E-05) and a D-index (robust hazard ratio) of 1.67 (P  =  9.2E-06), respectively. GO analysis of ODFSP's 2040 gene pairs (collapsed to 886 distinct genes) revealed the involvement in small molecular catabolic process, sulfur compound metabolic process, organic acid catabolic process, sulfur compound biosynthetic process, glycosaminoglycan metabolic process and aminometabolic process. Kyoto encyclopedia of genes and genomes pathway analysis of ODFSP's signature genes identified prominent pathways that included cAMP signaling pathway and FoxO signaling pathway. By identifying individuals who might benefit from a more aggressive treatment plan or enrolment in a clinical trial but who will not benefit from standard surgery or chemotherapy, ODFSP could help with treatment decisions.

Hereditary cancer syndromes

Hereditary cancer syndromes (HCSs) are arguably the most frequent category of Mendelian genetic diseases, as at least 2% of presumably healthy subjects carry highly-penetrant tumor-predisposing pathogenic variants (PVs). Hereditary breast-ovarian cancer and Lynch syndrome make the highest contribution to cancer morbidity; in addition, there are several dozen less frequent types of familial tumors. The development of the majority albeit not all hereditary malignancies involves two-hit mechanism, i.e. the somatic inactivation of the remaining copy of the affected gene. Earlier studies on cancer families suggested nearly fatal penetrance for the majority of HCS genes; however, population-based investigations and especially large-scale next-generation sequencing data sets demonstrate that the presence of some highly-penetrant PVs is often compatible with healthy status. Hereditary cancer research initially focused mainly on cancer detection and prevention. Recent studies identified multiple HCS-specific drug vulnerabilities, which translated into the development of highly efficient therapeutic options.

Genetic analyses of DNA repair pathway associated genes implicate new candidate cancer predisposing genes in ancestrally defined ovarian cancer cases

Not all familial ovarian cancer (OC) cases are explained by pathogenic germline variants in known risk genes. A candidate gene approach involving DNA repair pathway genes was applied to identify rare recurring pathogenic variants in familial OC cases not associated with known OC risk genes from a population exhibiting genetic drift. Whole exome sequencing (WES) data of 15 OC cases from 13 families tested negative for pathogenic variants in known OC risk genes were investigated for candidate variants in 468 DNA repair pathway genes. Filtering and prioritization criteria were applied to WES data to select top candidates for further analyses. Candidates were genotyped in ancestry defined study groups of 214 familial and 998 sporadic OC or breast cancer (BC) cases and 1025 population-matched controls and screened for additional carriers in 605 population-matched OC cases. The candidate genes were also analyzed in WES data from 937 familial or sporadic OC cases of diverse ancestries. Top candidate variants in ERCC5, EXO1, FANCC, NEIL1 and NTHL1 were identified in 5/13 (39%) OC families. Collectively, candidate variants were identified in 7/435 (1.6%) sporadic OC cases and 1/566 (0.2%) sporadic BC cases versus 1/1025 (0.1%) controls. Additional carriers were identified in 6/605 (0.9%) OC cases. Tumour DNA from ERCC5, NEIL1 and NTHL1 variant carriers exhibited loss of the wild-type allele. Carriers of various candidate variants in these genes were identified in 31/937 (3.3%) OC cases of diverse ancestries versus 0-0.004% in cancer-free controls. The strategy of applying a candidate gene approach in a population exhibiting genetic drift identified new candidate OC predisposition variants in DNA repair pathway genes.

Mutational Signature and Integrative Genomic Analysis of Human Papillomavirus-Associated Penile Squamous Cell Carcinomas from Latin American Patients

Simple Summary

DNA sequencing has been crucial to comprehending cancer mutational patterns, leading to the identification of driver genes and altered signaling pathways. Thus, identifying new pathogenic variants and their impact on tumor onset, progression, and treatment response has fueled tumor biology research. Here, we present novel findings addressing the first whole-exome sequencing (WES) of human papillomavirus (HPV)-associated penile squamous cell carcinoma (PSCC) from Latin Americans and its association with pathogenesis. We also compared the molecular profile of the tumors to that of three previous studies from populations with different genetic and socioeconomic backgrounds, the majority of which was HPV-negative. We describe the most altered genes and the main pathogenic variants found in the Latin Americans, ten of which are exclusive to our study sample. The data allowed us to identify molecular pathways and druggable targets with potential treatment value for this still-neglected HPV-associated carcinoma.

Abstract

High-throughput DNA sequencing has allowed for the identification of genomic alterations and their impact on tumor development, progression, and therapeutic responses. In PSCC, for which the incidence has progressively increased worldwide, there are still limited data on the molecular mechanisms involved in the disease pathogenesis. In this study, we characterized the mutational signature of 30 human papillomavirus (HPV)-associated PSCC cases from Latin Americans, using whole-exome sequencing. Copy number variations (CNVs) were also identified and compared to previous array-generated data. Enrichment analyses were performed to reveal disrupted pathways and to identify alterations mapped to HPV integration sites (HPVis) and miRNA–mRNA hybridization regions. Among the most frequently mutated genes were NOTCH1, TERT, TTN, FAT1, TP53, CDKN2A, RYR2, CASP8, FBXW7, HMCN2, and ITGA8. Of note, 92% of these altered genes were localized at HPVis. We also found mutations in ten novel genes (KMT2C, SMARCA4, PTPRB, AJUBA, CR1, KMT2D, NBEA, FAM135B, GTF2I, and CIC), thus increasing our understanding of the potential HPV-disrupted pathways. Therefore, our study reveals innovative targets with potential therapeutic benefits for HPV-associated PSCCs. The CNV analysis by sequencing (CNV-seq) revealed five cancer-associated genes as the most frequent with gains (NOTCH1, MYC, NUMA1, PLAG1, and RAD21), while 30% of the tumors showed SMARCA4 with loss. Additionally, four cancer-associated genes (CARD11, CSMD3, KDR, and TLX3) carried untranslated regions (UTRs) variants, which may impact gene regulation by affecting the miRNAs hybridization regions. Altogether, these data contribute to the characterization of the mutational spectrum and its impact on cellular signaling pathways in PSCC, thus reinforcing the pivotal role of HPV infection in the molecular pathogenesis of these tumors.

Whole-Genome Sequencing Identifies PPARGC1A as a Putative Modifier of Cancer Risk in BRCA1/2 Mutation Carriers

While BRCA1 and BRCA2 mutations are known to confer the largest risk of breast cancer and ovarian cancer, the incomplete penetrance of the mutations and the substantial variability in age at cancer onset among carriers suggest additional factors modifying the risk of cancer in BRCA1/2 mutation carriers. To identify genetic modifiers of BRCA1/2, we carried out a whole-genome sequencing study of 66 ovarian cancer patients that were enriched with BRCA carriers, followed by validation using data from the Pan-Cancer Analysis of Whole Genomes Consortium. We found PPARGC1A, a master regulator of mitochondrial biogenesis and function, to be highly mutated in BRCA carriers, and patients with both PPARGC1A and BRCA1/2 mutations were diagnosed with breast or ovarian cancer at significantly younger ages, while the mutation status of each gene alone did not significantly associate with age of onset. Our study suggests PPARGC1A as a possible BRCA modifier gene. Upon further validation, this finding can help improve cancer risk prediction and provide personalized preventive care for BRCA carriers.

Identification of Tumor Microenvironment-Related Prognostic Biomarkers for Ovarian Serous Cancer 3-Year Mortality Using Targeted Maximum Likelihood Estimation: A TCGA Data Mining Study

Ovarian serous cancer (OSC) is one of the leading causes of death across the world. The role of the tumor microenvironment (TME) in OSC has received increasing attention. Targeted maximum likelihood estimation (TMLE) is developed under a counterfactual framework to produce effect estimation for both the population level and individual level. In this study, we aim to identify TME-related genes and using the TMLE method to estimate their effects on the 3-year mortality of OSC. In total, 285 OSC patients from the TCGA database constituted the studying population. ESTIMATE algorithm was implemented to evaluate immune and stromal components in TME. Differential analysis between high-score and low-score groups regarding ImmuneScore and StromalScore was performed to select shared differential expressed genes (DEGs). Univariate logistic regression analysis was followed to evaluate associations between DEGs and clinical pathologic factors with 3-year mortality. TMLE analysis was conducted to estimate the average effect (AE), individual effect (IE), and marginal odds ratio (MOR). The validation was performed using three datasets from Gene Expression Omnibus (GEO) database. Additionally, 355 DEGs were selected after differential analysis, and 12 genes from DEGs were significant after univariate logistic regression. Four genes remained significant after TMLE analysis. In specific, ARID3C and FREM2 were negatively correlated with OSC 3-year mortality. CROCC2 and PTF1A were positively correlated with OSC 3-year mortality. Combining of ESTIMATE algorithm and TMLE algorithm, we identified four TME-related genes in OSC. AEs were estimated to provide averaged effects based on the population level, while IEs were estimated to provide individualized effects and may be helpful for precision medicine.

Population-based targeted sequencing of 54 candidate genes identifies PALB2 as a susceptibility gene for high-grade serous ovarian cancer

PURPOSE

The known epithelial ovarian cancer (EOC) susceptibility genes account for less than 50% of the heritable risk of ovarian cancer suggesting that other susceptibility genes exist. The aim of this study was to evaluate the contribution to ovarian cancer susceptibility of rare deleterious germline variants in a set of candidate genes.

METHODS

We sequenced the coding region of 54 candidate genes in 6385 invasive EOC cases and 6115 controls of broad European ancestry. Genes with an increased frequency of putative deleterious variants in cases versus controls were further examined in an independent set of 14 135 EOC cases and 28 655 controls from the Ovarian Cancer Association Consortium and the UK Biobank. For each gene, we estimated the EOC risks and evaluated associations between germline variant status and clinical characteristics.

RESULTS

The ORs associated for high-grade serous ovarian cancer were 3.01 for PALB2 (95% CI 1.59 to 5.68; p=0.00068), 1.99 for POLK (95% CI 1.15 to 3.43; p=0.014) and 4.07 for SLX4 (95% CI 1.34 to 12.4; p=0.013). Deleterious mutations in FBXO10 were associated with a reduced risk of disease (OR 0.27, 95% CI 0.07 to 1.00, p=0.049). However, based on the Bayes false discovery probability, only the association for PALB2 in high-grade serous ovarian cancer is likely to represent a true positive.

CONCLUSIONS

We have found strong evidence that carriers of PALB2 deleterious mutations are at increased risk of high-grade serous ovarian cancer. Whether the magnitude of risk is sufficiently high to warrant the inclusion of PALB2 in cancer gene panels for ovarian cancer risk testing is unclear; much larger sample sizes will be needed to provide sufficiently precise estimates for clinical counselling.

Whole-Exome Sequencing Reveals Novel Variations in Patients with Familial Von Hippel-Lindau Syndrome

OBJECTIVE

Von Hippel-Lindau (VHL) syndrome is a rare disease that occurs in an autosomal-dominant genetic pattern. Due to the high genetic variability of VHL diseases, current studies have limited clinical value. Moreover, casual genetic variations in patients with VHL syndrome are still unclear.

METHODS

Here, we performed whole-exome sequencing of 25 individuals to identify reliable disease-related variations. Systemic computational analysis was performed for variant detection, and Sanger sequencing was used to validate detected mutations.

RESULTS

Most of the known mutations in the VHL gene were observed in the studied population. In addition, a large fragment deletion in VHL exon 2 in the immediate family members of the last family was detected. This had not been reported earlier. Moreover, we identified 3 novel mutation sites in the MAP2K3 gene that may be involved in the occurrence and development of the VHL disease.

CONCLUSIONS

These results demonstrated that the heterogeneous nature of VHL syndrome and novel mutational signatures may help to improve the diagnostic ability of VHL syndrome.

Peripheral blood BRCA1 methylation profiling to predict familial ovarian cancer

Objective

Familial cancer appears at a young age and its incidence is increasing. About 12% of familial ovarian cancer cases are associated with BRCA1/2 mutations (BRCAm). In this study, we investigated BRCA1 methylation may predict ovarian cancer in those with a family history of cancer (FHC) but without BRCA1/2 mutations (BRCAwt).

Methods

Using peripheral blood DNA from 55 subjects without a history of cancer [cancer(−)] and 52 ovarian cancer patients, we examined BRCA1 promoter methylation through bisulfite sequencing of the promoter and expressed the results as the cumulative methylation index. Then, we evaluated the BRCA1 promoter methylation according to BRCA1/2 germline mutations.

Results

BRCA1 methylation was more prevalent in the BRCAm cancer(−) group than in the BRCAwt cancer(−) group and ovarian cancer patients (p=0.031 and p=0.019, respectively). In the BRCAwt cancer(−) group, BRCA1 methylation was more prevalent in those with an FHC than in those without one and in the BRCAm cancer(−) group with an FHC (p=0.001 and p<0.001, respectively).

Conclusion

Our data suggest a predictive role of BRCA1 methylation profile for ovarian cancer in those without a history of cancer but with an FHC. BRCA1 methylation has important implications for diagnostic and predictive testing of those with BRCAwt cancer(−) status with FHC.

An Overview of Candidate Therapeutic Target Genes in Ovarian Cancer

Ovarian cancer (OC) shows the highest mortality rate among gynecological malignancies and, because of the absence of specific symptoms, it is frequently diagnosed at an advanced stage, mainly due to the lack of specific and early biomarkers, such as those based on cancer molecular signature identification. Indeed, although significant progress has been made toward improving the clinical outcome of other cancers, rates of mortality for OC are essentially unchanged since 1980, suggesting the need of new approaches to identify and characterize the molecular mechanisms underlying pathogenesis and progression of these malignancies. In addition, due to the low response rate and the high frequency of resistance to current treatments, emerging therapeutic strategies against OC focus on targeting single factors and pathways specifically involved in tumor growth and metastasis. To date, loss-of-function screenings are extensively applied to identify key drug targets in cancer, seeking for more effective, disease-tailored treatments to overcome lack of response or resistance to current therapies. We review here the information relative to essential genes and functional pathways recently discovered in OC, often strictly interconnected with each other and representing promising biomarkers and molecular targets to treat these malignancies.

Exome sequencing of familial high-grade serous ovarian carcinoma reveals heterogeneity for rare candidate susceptibility genes

High-grade serous ovarian carcinoma (HGSOC) has a significant hereditary component, approximately half of which cannot be explained by known genes. To discover genes, we analyse germline exome sequencing data from 516 BRCA1/2-negative women with HGSOC, focusing on genes enriched with rare, protein-coding loss-of-function (LoF) variants. Overall, there is a significant enrichment of rare protein-coding LoF variants in the cases (p < 0.0001, chi-squared test). Only thirty-four (6.6%) have a pathogenic variant in a known or proposed predisposition gene. Few genes have LoF mutations in more than four individuals and the majority are detected in one individual only. Forty-three highly-ranked genes are identified with three or more LoF variants that are enriched by three-fold or more compared to GnomAD. These genes represent diverse functional pathways with relatively few involved in DNA repair, suggesting that much of the remaining heritability is explained by previously under-explored genes and pathways.