Germline Sequencing Analysis to Inform Clinical Gene Panel Testing for Aggressive Prostate Cancer

Importance

Germline gene panel testing is recommended for men with advanced prostate cancer (PCa) or a family history of cancer. While evidence is limited for some genes currently included in panel testing, gene panels are also likely to be incomplete and missing genes that influence PCa risk and aggressive disease.

Objective

To identify genes associated with aggressive PCa.

Design, Setting, and Participants

A 2-stage exome sequencing case-only genetic association study was conducted including men of European ancestry from 18 international studies. Data analysis was performed from January 2021 to March 2023. Participants were 9185 men with aggressive PCa (including 6033 who died of PCa and 2397 with confirmed metastasis) and 8361 men with nonaggressive PCa.

Exposure

Sequencing data were evaluated exome-wide and in a focused investigation of 29 DNA repair pathway and cancer susceptibility genes, many of which are included on gene panels.

Main Outcomes and Measures

The primary study outcomes were aggressive (category T4 or both T3 and Gleason score ≥8 tumors, metastatic PCa, or PCa death) vs nonaggressive PCa (category T1 or T2 and Gleason score ≤6 tumors without known recurrence), and metastatic vs nonaggressive PCa.

Results

A total of 17 546 men of European ancestry were included in the analyses; mean (SD) age at diagnosis was 65.1 (9.2) years in patients with aggressive PCa and 63.7 (8.0) years in those with nonaggressive disease. The strongest evidence of association with aggressive or metastatic PCa was noted for rare deleterious variants in known PCa risk genes BRCA2 and ATM (P ≤ 1.9 × 10-6), followed by NBN (P = 1.7 × 10-4). This study found nominal evidence (P < .05) of association with rare deleterious variants in MSH2, XRCC2, and MRE11A. Five other genes had evidence of greater risk (OR≥2) but carrier frequency differences between aggressive and nonaggressive PCa were not statistically significant: TP53, RAD51D, BARD1, GEN1, and SLX4. Deleterious variants in these 11 candidate genes were carried by 2.3% of patients with nonaggressive, 5.6% with aggressive, and 7.0% with metastatic PCa.

Conclusions and Relevance

The findings of this study provide further support for DNA repair and cancer susceptibility genes to better inform disease management in men with PCa and for extending testing to men with nonaggressive disease, as men carrying deleterious alleles in these genes are likely to develop more advanced disease.

Genetic Modifiers of Cystic Fibrosis Lung Disease Severity: Whole-Genome Analysis of 7,840 Patients

Rationale: Lung disease is the major cause of morbidity and mortality in persons with cystic fibrosis (pwCF). Variability in CF lung disease has substantial non-CFTR (CF transmembrane conductance regulator) genetic influence. Identification of genetic modifiers has prognostic and therapeutic importance. Objectives: Identify genetic modifier loci and genes/pathways associated with pulmonary disease severity. Methods: Whole-genome sequencing data on 4,248 unique pwCF with pancreatic insufficiency and lung function measures were combined with imputed genotypes from an additional 3,592 patients with pancreatic insufficiency from the United States, Canada, and France. This report describes association of approximately 15.9 million SNPs using the quantitative Kulich normal residual mortality-adjusted (KNoRMA) lung disease phenotype in 7,840 pwCF using premodulator lung function data. Measurements and Main Results: Testing included common and rare SNPs, transcriptome-wide association, gene-level, and pathway analyses. Pathway analyses identified novel associations with genes that have key roles in organ development, and we hypothesize that these genes may relate to dysanapsis and/or variability in lung repair. Results confirmed and extended previous genome-wide association study findings. These whole-genome sequencing data provide finely mapped genetic information to support mechanistic studies. No novel primary associations with common single variants or rare variants were found. Multilocus effects at chr5p13 (SLC9A3/CEP72) and chr11p13 (EHF/APIP) were identified. Variant effect size estimates at associated loci were consistently ordered across the cohorts, indicating possible age or birth cohort effects. Conclusions: This premodulator genomic, transcriptomic, and pathway association study of 7,840 pwCF will facilitate mechanistic and postmodulator genetic studies and the development of novel therapeutics for CF lung disease.

Clonal hematopoiesis and risk of prostate cancer in large samples of European ancestry men

Little is known regarding the potential relationship between clonal hematopoiesis (CH) of indeterminate potential (CHIP), which is the expansion of hematopoietic stem cells with somatic mutations, and risk of prostate cancer, the fifth leading cause of cancer death of men worldwide. We evaluated the association of age-related CHIP with overall and aggressive prostate cancer risk in two large whole-exome sequencing studies of 75 047 European ancestry men, including 7663 prostate cancer cases, 2770 of which had aggressive disease, and 3266 men carrying CHIP variants. We found that CHIP, defined by over 50 CHIP genes individually and in aggregate, was not significantly associated with overall (aggregate HR = 0.93, 95% CI = 0.76-1.13, P = 0.46) or aggressive (aggregate OR = 1.14, 95% CI = 0.92-1.41, P = 0.22) prostate cancer risk. CHIP was weakly associated with genetic risk of overall prostate cancer, measured using a polygenic risk score (OR = 1.05 per unit increase, 95% CI = 1.01-1.10, P = 0.01). CHIP was not significantly associated with carrying pathogenic/likely pathogenic/deleterious variants in DNA repair genes, which have previously been found to be associated with aggressive prostate cancer. While findings from this study suggest that CHIP is likely not a risk factor for prostate cancer, it will be important to investigate other types of CH in association with prostate cancer risk.

Abstract 688: Multi-stage exome sequencing study of 17,546 aggressive and non-aggressive prostate cancer cases

Rare pathogenic variants in DNA repair genes have been found to influence risk of aggressive prostate cancer. We conducted a large case-only exome sequencing study to further understand the role of rare coding variation in aggressive prostate cancer in a study of 9,185 aggressive (prostate cancer death, metastatic disease, T4, or both T3 and Gleason≥8) and 8,361 non-aggressive cases (T1/T2 and Gleason≤6) of European ancestry from 19 international studies. Stage 1 samples (N=5,545) had whole exome-sequencing, and stage 2 samples (N=12,001) had targeted exome sequencing for 1,459 genes selected based on stage 1 results and previous evidence. Logistic regression models were used to evaluate gene-based tests and the aggregate effect of multiple genes to investigate whether carrying pathogenic/likely pathogenic/deleterious (P/LP/D) variants (18,759 identified) was associated with risk of aggressive prostate cancer, prostate cancer death (N=6,033), or metastatic disease (N=1,730) compared to non-aggressive disease. Gene-based tests were meta-analyzed across stages 1 and 2. BRCA2, ATM, and NBN had the most statistically significant gene-based results: BRCA2 P/LP/D variant carriers had 4.3-fold higher odds of aggressive disease (95% CI=3.15-5.86, P=4x10-20), 4.7-fold higher odds of prostate cancer death (95% CI=3.41-6.59, P=2x10-20), and 5.7-fold higher odds of metastatic disease (95% CI=3.71-8.76, P=2x10-15); ATM P/LP/D variant carriers had 2.2-fold higher odds of aggressive disease (95% CI=1.58-2.99, P=2x10-6), 2.2-fold higher odds of prostate cancer death (95% CI=1.52-3.05, P=2x10-5), and 3.0-fold higher odds of metastatic disease (95% CI=1.93-4.61, P=9x10-7); and NBN P/LP/D variant carriers had 5.9-fold higher odds of metastatic disease (95% CI=2.56-13.84, P=3x10-5). Among potentially novel genes with strong but not exome-wide significant statistical evidence were MMP19, involved in reproduction and metastasis, with carriers having 2.8-fold higher odds of prostate cancer death (95% CI=1.53-5.05, P=8x10-4); PKD2L2, involved in fertility, with carriers having 3.5-fold higher odds of prostate cancer death (95% CI=1.76-7.04, P=5x10-4); and SMPD1, involved in converting sphingomyelin to ceramide, with carriers having 5.3-fold higher odds of metastatic disease (95% CI=1.85-14.98, P=0.002). At least one P/LP/D variant within a subset of 24 previously curated candidate prostate cancer DNA repair genes was carried by 12.8% of aggressive cases (OR=1.48, 95% CI=1.34-1.64, P=3x10-14), 12.6% of cases who died due to prostate cancer (OR=1.47, 95% CI=1.31-1.65, P=3x10-11), and 15.1% of metastatic cases (OR=2.16, 95% CI=1.57-2.16, P=5x10-14) compared to 9.4% of non-aggressive cases. These findings support the importance of rare genetic variation in aggressive prostate cancer risk and may have important implications for prostate cancer risk stratification and screening.

Citation Format: Burcu F. Darst, Ed Saunders, Tokhir Dadaev, Xin Sheng, Peggy Wan, Loreall Pooler, Lucy Y. Xia, Stephen Chanock, Sonja I. Berndt, Susan M. Gapstur, Victoria Stevens, Demetrius Albanes, Stephanie J. Weinstein, Vincent Gnanapragasam, Graham G. Giles, Tu Nguyen-Dumont, Roger L. Milne, Mark M. Pomerantz, Julie A. Schmidt, Ruth C. Travis, Timothy J. Key, Konrad H. Stopsack, Lorelei A. Mucci, William J. Catalona, Beth Marosy, Kurt N. Hetrick, Kimberly F. Doheny, Robert J. MacInnis, Melissa C. Southey, Rosalind A. Eeles, Fredrik Wiklund, Zsofia Kote-Jarai, David V. Conti, Christopher A. Haiman. Multi-stage exome sequencing study of 17,546 aggressive and non-aggressive prostate cancer cases [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 688.

Germline Sequencing DNA Repair Genes in 5545 Men With Aggressive and Nonaggressive Prostate Cancer

BACKGROUND

There is an urgent need to identify factors specifically associated with aggressive prostate cancer (PCa) risk. We investigated whether rare pathogenic, likely pathogenic, or deleterious (P/LP/D) germline variants in DNA repair genes are associated with aggressive PCa risk in a case-case study of aggressive vs nonaggressive disease.

METHODS

Participants were 5545 European-ancestry men, including 2775 nonaggressive and 2770 aggressive PCa cases, which included 467 metastatic cases (16.9%). Samples were assembled from 12 international studies and germline sequenced together. Rare (minor allele frequency < 0.01) P/LP/D variants were analyzed for 155 DNA repair genes. We compared single variant, gene-based, and DNA repair pathway-based burdens by disease aggressiveness. All statistical tests are 2-sided.

RESULTS

BRCA2 and PALB2 had the most statistically significant gene-based associations, with 2.5% of aggressive and 0.8% of nonaggressive cases carrying P/LP/D BRCA2 alleles (odds ratio [OR] = 3.19, 95% confidence interval [CI] = 1.94 to 5.25, P = 8.58 × 10-7) and 0.65% of aggressive and 0.11% of nonaggressive cases carrying P/LP/D PALB2 alleles (OR = 6.31, 95% CI = 1.83 to 21.68, P = 4.79 × 10-4). ATM had a nominal association, with 1.6% of aggressive and 0.8% of nonaggressive cases carrying P/LP/D ATM alleles (OR = 1.88, 95% CI = 1.10 to 3.22, P = .02). In aggregate, P/LP/D alleles within 24 literature-curated candidate PCa DNA repair genes were more common in aggressive than nonaggressive cases (carrier frequencies = 14.2% vs 10.6%, respectively; P = 5.56 × 10-5). However, this difference was non-statistically significant (P = .18) on excluding BRCA2, PALB2, and ATM. Among these 24 genes, P/LP/D carriers had a 1.06-year younger diagnosis age (95% CI = -1.65 to 0.48, P = 3.71 × 10-4).

CONCLUSIONS

Risk conveyed by DNA repair genes is largely driven by rare P/LP/D alleles within BRCA2, PALB2, and ATM. These findings support the importance of these genes in both screening and disease management considerations.

Multi-Omics Analysis Reveals a HIF Network and Hub Gene EPAS1 Associated with Lung Adenocarcinoma

Recent technological advancements have permitted high-throughput measurement of the human genome, epigenome, metabolome, transcriptome, and proteome at the population level. We hypothesized that subsets of genes identified from omic studies might have closely related biological functions and thus might interact directly at the network level. Therefore, we conducted an integrative analysis of multi-omic datasets of non-small cell lung cancer (NSCLC) to search for association patterns beyond the genome and transcriptome. A large, complex, and robust gene network containing well-known lung cancer-related genes, including EGFR and TERT, was identified from combined gene lists for lung adenocarcinoma. Members of the hypoxia-inducible factor (HIF) gene family were at the center of this network. Subsequent sequencing of network hub genes within a subset of samples from the Transdisciplinary Research in Cancer of the Lung-International Lung Cancer Consortium (TRICL-ILCCO) consortium revealed a SNP (rs12614710) in EPAS1 associated with NSCLC that reached genome-wide significance (OR = 1.50; 95% CI: 1.31-1.72; p = 7.75 × 10^-9). Using imputed data, we found that this SNP remained significant in the entire TRICL-ILCCO consortium (p = .03). Additional functional studies are warranted to better understand interrelationships among genetic polymorphisms, DNA methylation status, and EPAS1 expression.

Generating Exome Enriched Sequencing Libraries from Formalin-Fixed, Paraffin-Embedded Tissue DNA for Next-Generation Sequencing

This unit describes a technique for generating exome-enriched sequencing libraries using DNA extracted from formalin-fixed paraffin-embedded (FFPE) samples. Utilizing commercially available kits, we present a low-input FFPE workflow starting with 50 ng of DNA. This procedure includes a repair step to address damage caused by FFPE preservation that improves sequence quality. Subsequently, libraries undergo an in-solution-targeted selection for exons, followed by sequencing using the Illumina next-generation short-read sequencing platform. © 2017 by John Wiley & Sons, Inc.

Novel Deletion of SERPINF1 Causes Autosomal Recessive Osteogenesis Imperfecta Type VI in Two Brazilian Families

Autosomal recessive osteogenesis imperfecta (OI) accounts for 10% of all OI cases, and, currently, mutations in 10 genes (CRTAP, LEPRE1, PPIB, SERPINH1, FKBP10, SERPINF1, SP7, BMP1, TMEM38B, and WNT1) are known to be responsible for this form of the disease. PEDF is a secreted glycoprotein of the serpin superfamily that maintains bone homeostasis and regulates osteoid mineralization, and it is encoded by SERPINF1, currently associated with OI type VI (MIM 172860). Here, we report a consanguineous Brazilian family in which multiple individuals from at least 4 generations are affected with a severe form of OI, and we also report an unrelated individual from the same small city in Brazil with a similar but more severe phenotype. In both families the same homozygous SERPINF1 19-bp deletion was identified which is not known in the literature yet. We described intra- and interfamilial clinical and radiological phenotypic variability of OI type VI caused by the same homozygous SERPINF1 19-bp deletion and suggest a founder effect. Furthermore, the SERPINF1 genotypes/phenotypes reported so far in the literature are reviewed.

Mutations in Alström protein impair terminal differentiation of cardiomyocytes

Cardiomyocyte cell division and replication in mammals proceed through embryonic development and abruptly decline soon after birth. The process governing cardiomyocyte cell cycle arrest is poorly understood. Here we carry out whole exome sequencing in an infant with evidence of persistent postnatal cardiomyocyte replication to determine the genetic risk factors. We identify compound heterozygous ALMS1 mutations in the proband, and confirm their presence in her affected sibling, one copy inherited from each heterozygous parent. Next, we recognise homozygous or compound heterozygous truncating mutations in ALMS1 in four other children with high levels of postnatal cardiomyocyte proliferation. Alms1 mRNA knockdown increases multiple markers of proliferation in cardiomyocytes, the percentage of cardiomyocytes in G2/M phases, and the number of cardiomyocytes by 10% in cultured cells. Homozygous Alms1-mutant mice have increased cardiomyocyte proliferation at two weeks postnatal compared to wild-type littermates. We conclude that deficiency of Alström protein impairs postnatal cardiomyocyte cell cycle arrest.

Detecting and estimating contamination of human DNA samples in sequencing and array-based genotype data

DNA sample contamination is a serious problem in DNA sequencing studies and may result in systematic genotype misclassification and false positive associations. Although methods exist to detect and filter out cross-species contamination, few methods to detect within-species sample contamination are available. In this paper, we describe methods to identify within-species DNA sample contamination based on (1) a combination of sequencing reads and array-based genotype data, (2) sequence reads alone, and (3) array-based genotype data alone. Analysis of sequencing reads allows contamination detection after sequence data is generated but prior to variant calling; analysis of array-based genotype data allows contamination detection prior to generation of costly sequence data. Through a combination of analysis of in silico and experimentally contaminated samples, we show that our methods can reliably detect and estimate levels of contamination as low as 1%. We evaluate the impact of DNA contamination on genotype accuracy and propose effective strategies to screen for and prevent DNA contamination in sequencing studies.

Copy number variation in familial Parkinson disease

Copy number variants (CNVs) are known to cause Mendelian forms of Parkinson disease (PD), most notably in SNCA and PARK2. PARK2 has a recessive mode of inheritance; however, recent evidence demonstrates that a single CNV in PARK2 (but not a single missense mutation) may increase risk for PD. We recently performed a genome-wide association study for PD that excluded individuals known to have either a LRRK2 mutation or two PARK2 mutations. Data from the Illumina370Duo arrays were re-clustered using only white individuals with high quality intensity data, and CNV calls were made using two algorithms, PennCNV and QuantiSNP. After quality assessment, the final sample included 816 cases and 856 controls. Results varied between the two CNV calling algorithms for many regions, including the PARK2 locus (genome-wide p = 0.04 for PennCNV and p = 0.13 for QuantiSNP). However, there was consistent evidence with both algorithms for two novel genes, USP32 and DOCK5 (empirical, genome-wide p-values<0.001). PARK2 CNVs tended to be larger, and all instances that were molecularly tested were validated. In contrast, the CNVs in both novel loci were smaller and failed to replicate using real-time PCR, MLPA, and gel electrophoresis. The DOCK5 variation is more akin to a VNTR than a typical CNV and the association is likely caused by artifact due to DNA source. DNA for all the cases was derived from whole blood, while the DNA for all controls was derived from lymphoblast cell lines. The USP32 locus contains many SNPs with low minor allele frequency leading to a loss of heterozygosity that may have been spuriously interpreted by the CNV calling algorithms as support for a deletion. Thus, only the CNVs within the PARK2 locus could be molecularly validated and associated with PD susceptibility.

R/Bioconductor software for Illumina's Infinium whole-genome genotyping BeadChips

Summary: Illumina produces a number of microarray-based technologies for human genotyping. An Infinium BeadChip is a two-color platform that types between 10⁵ and 10⁶ single nucleotide polymorphisms (SNPs) per sample. Despite being widely used, there is a shortage of open source software to process the raw intensities from this platform into genotype calls. To this end, we have developed the R/Bioconductor package crlmm for analyzing BeadChip data. After careful preprocessing, our software applies the CRLMM algorithm to produce genotype calls, confidence scores and other quality metrics at both the SNP and sample levels. We provide access to the raw summary-level intensity data, allowing users to develop their own methods for genotype calling or copy number analysis if they wish.

Availability and Implementation: The crlmm Bioconductor package is available from http://www.bioconductor.org. Data packages and documentation are available from http://rafalab.jhsph.edu/software.html.

Contact:mritchie@wehi.edu.au; rafa@jhu.edu

A genome-wide association study of type 2 diabetes in Finns detects multiple susceptibility variants

Identifying the genetic variants that increase the risk of type 2 diabetes (T2D) in humans has been a formidable challenge. Adopting a genome-wide association strategy, we genotyped 1161 Finnish T2D cases and 1174 Finnish normal glucose-tolerant (NGT) controls with >315,000 single-nucleotide polymorphisms (SNPs) and imputed genotypes for an additional >2 million autosomal SNPs. We carried out association analysis with these SNPs to identify genetic variants that predispose to T2D, compared our T2D association results with the results of two similar studies, and genotyped 80 SNPs in an additional 1215 Finnish T2D cases and 1258 Finnish NGT controls. We identify T2D-associated variants in an intergenic region of chromosome 11p12, contribute to the identification of T2D-associated variants near the genes IGF2BP2 and CDKAL1 and the region of CDKN2A and CDKN2B, and confirm that variants near TCF7L2, SLC30A8, HHEX, FTO, PPARG, and KCNJ11 are associated with T2D risk. This brings the number of T2D loci now confidently identified to at least 10.