Genetic architecture of prostate cancer in the Ashkenazi Jewish population

SNPs at SMG7 Associated with Time from Biochemical Recurrence to Prostate Cancer Death

BACKGROUND

A previous genome-wide association study identified several loci with genetic variants associated with prostate cancer survival time in two cohorts from Sweden. Whether these variants have an effect in other populations or if their effect is homogenous across the course of disease is unknown.

METHODS

These variants were genotyped in a cohort of 1,298 patients. Samples were linked with age, PSA level, Gleason score, cancer stage at surgery, and times from surgery to biochemical recurrence to death from prostate cancer. SNPs rs2702185 and rs73055188 were tested for association with prostate cancer-specific survival time using a multivariate Cox proportional hazard model. SNP rs2702185 was further tested for association with time to biochemical recurrence and time from biochemical recurrence to death with a multi-state model.

RESULTS

SNP rs2702185 at SMG7 was associated with prostate cancer-specific survival time, specifically the time from biochemical recurrence to prostate cancer death (HR, 2.5; 95% confidence interval, 1.4-4.5; P = 0.0014). Nine variants were in linkage disequilibrium (LD) with rs2702185; one, rs10737246, was found to be most likely to be functional based on LD patterns and overlap with open chromatin. Patterns of open chromatin and correlation with gene expression suggest that this SNP may affect expression of SMG7 in T cells.

CONCLUSIONS

The SNP rs2702185 at the SMG7 locus is associated with time from biochemical recurrence to prostate cancer death, and its LD partner rs10737246 is predicted to be functional.

IMPACT

These results suggest that future association studies of prostate cancer survival should consider various intervals over the course of disease.

Clinical correlations and prognostic value of Nudix hydroxylase 10 in patients with gastric cancer

Gastric cancer (GC) is one of the most common and lethal cancers worldwide. The Nudix hydroxylase (NUDT) genes have been reported to play notable roles in tumor progression. However, the role of NUDT10 in GC has not been reported. In this study, we investigated the expression of NUDT10 in GC and its association with clinicopathological characteristics. Quantitative real-time polymerase chain reaction and analyses of The Cancer Genome Atlas and Human Protein Atlas databases were performed to determine NUDT10 mRNA and protein expression. Receiver operating characteristic curve analysis was used to assess the diagnostic value of NUDT10 in patients with GC. We used Cox regression and the Kaplan-Meier method to assess the correlations between clinicopathological factors and survival outcomes of patients with GC. Gene set enrichment analysis (GSEA) was performed to identify the underlying signaling pathways. NUDT10 mRNA and protein expression was significantly lower in GC tissues compared to normal tissues. Interestingly, higher NUDT10 expression was correlated with advanced tumor stage, deeper local invasion, and worse survival outcomes. Patients with higher NUDT10 expression had a significantly worse prognosis than those with lower NUDT10 expression. Multivariate analysis showed that high NUDT10 expression was an independent predictor of survival outcome. Several pathways, including mismatch repair, nucleotide excision repair, extracellular matrix receptor interaction, and cancer signaling, were identified as enriched pathways in GC through GSEA. To our knowledge, this study is the first to characterize NUDT10 expression in GC. Our study demonstrates that NUDT10 is a promising independent biomarker for GC prognosis.

PRNCR1: a long non-coding RNA with a pivotal oncogenic role in cancer

Long non-coding RNAs (lncRNAs) have been gaining importance in the field of cancer research in recent years. PRNCR1 (prostate cancer-associated non-coding RNA1) is a 12.7 kb, intron-less lncRNA found to play an oncogenic role in malignancy of diverse organs including prostate, breast, lung, oral cavity, colon and rectum. Single-nucleotide polymorphisms (SNPs) of PRNCR1 locus have been found to be associated with cancer susceptibility in different populations. In this review, an attempt has been made for the first time to summarize all sorts of available data on PRNCR1 to date from relevant databases (GeneCard, LncExpDB, Ensembl genome browser, and PubMed). As functional roles of PRNCR1, miRNA (microRNA) sponging was mostly highlighted in the pathogenesis of different cancer; in addition, an association of the lncRNA with chromatin-modifying complex to enhance androgen receptor-mediated gene transcription was reported in prostate cancer. Diagnostic and prognostic importance of PRNCR1 was found in some malignancies suggesting potency of the lncRNA to serve as a clinical biomarker. For PRNCR1 SNPs, although cancer susceptibility of the risk alleles/genotypes was reported in different populations, majorities of the findings were not replicated and underlying molecular mechanisms remained unexplored. Therapeutic implication of PRNCR1 was not studied well and future research may come up in this direction for intervening novel strategies to fight against cancer.

Validation of prostate cancer risk variants rs10993994 and rs7098889 by CRISPR/Cas9 mediated genome editing

GWAS have identified numerous SNPs associated with prostate cancer risk. One such SNP is rs10993994. It is located in the β-microseminoprotein (MSMB) promoter region, mediates MSMB prostate secretion levels, and is linked to mRNA expression changes in both MSMB and the adjacent gene NCOA4. In addition, our previous work showed a second SNP, rs7098889, is in positive linkage disequilibrium with rs10993994 and associated with MSMB expression independent of rs10993994. Here, we generate a series of clones with single alleles removed by double guide RNA (gRNA) mediated CRISPR/Cas9 deletions, through which we demonstrate that each of these SNPs independently and greatly alters MSMB expression in an allele-specific manner. We further show that these SNPs have no substantial effect on the expression of NCOA4. These data demonstrate that a single SNP can have a large effect on gene expression and illustrate the importance of functional validation studies to deconvolute observed correlations. The method we have developed is generally applicable to test any SNP for which a relevant heterozygous cell line is available. AUTHOR SUMMARY: In pursuing the underlying biological mechanism of prostate cancer pathogenesis, scientists utilized the existence of common single nucleotide polymorphisms (SNPs) in the human genome as genetic markers to perform large scale genome wide association studies (GWAS) and have so far identified more than a hundred prostate cancer risk variants. Such variants provide an unbiased and systematic new venue to study the disease mechanism, and the next big challenge is to translate these genetic associations to the causal role of altered gene function in oncogenesis. The majority of these variants are waiting to be studied and lots of them may act in oncogenesis through gene expression regulation. To prove the concept, we took rs10993994 and its linked rs7098889 as an example and engineered single cell clones by allelic-specific CRISPR/Cas9 deletion to separate the effect of each allele. We observed that a single nucleotide difference would lead to surprisingly high level of MSMB gene expression change in a gene specific and cell-type specific manner. Our study strongly supports the notion that differential level of gene expression caused by risk variants and their associated genetic locus play a major role in oncogenesis and also highlights the importance of studying the function of MSMB encoded β-MSP in prostate cancer pathogenesis.

Fumarate hydratase FH c.1431_1433dupAAA (p.Lys477dup) variant is not associated with cancer including renal cell carcinoma

Fumarate hydratase (FH) mutations underpin the autosomal recessive syndrome. FH deficiency and the autosomal dominant syndrome hereditary leiomyomatosis and renal cell carcinoma (HLRCC). The FH c.1431_1433dupAAA (p.Lys477dup) genomic alteration has been conclusively shown to contribute to FH deficiency when occurring with another FH germline alteration. However, a sufficiently large dataset has been lacking to conclusively determine its clinical significance to cancer predisposition in the heterozygous state. We reviewed a series of 7,571 patients with cancer who received germline results through MSK-IMPACT testing at the Memorial Sloan Kettering Cancer Center. The FH c.1431_1433dupAAA (p.Lys477dup) variant was detected in 24 individuals, none of whom was affected with renal cancer. Eleven of the 372 patients with renal cancer were identified to carried pathogenic FH variants associated with HLRCC. None of these 372 patients with renal cancer carried the FH c.1431_1433dupAAA variant. Our data indicate the FH c.1431_1433dupAAA is not associated with cancer including renal cell carcinoma.

NUDT11 rs5945572 polymorphism and prostate cancer risk: a meta-analysis

The association between the NUDT10 rs5945572 polymorphism and prostate cancer (PCa) was not clear. We thus conducted a meta-analysis to assess the association between NUDT10 rs5945572 polymorphism and PCa risk. A literature search was carried out using PUBMED, EMBASE, and Cochrane Library Central database before Dec 2014. The strength of the associations between the NUDT10 rs5945572 polymorphism and PCa risk was measured by odds ratios (OR) with 95% confidence intervals (CI). The random-effects model was used. NUDT10 rs5945572 polymorphism was significantly associated with PCa risk (OR = 1.22, 95% CI 1.19-1.26, P < 0.001, I(2) = 0%, Figure 2). In the subgroup analysis by ethnicity, a significant association was found among Caucasians (OR = 1.25, 95% CI 1.00-1.57, P = 0.05, I(2) = 0%), and Asians (OR = 1.23, 95% CI 1.19-1.28, P < 0.001, I(2) = 0%), and Africans (OR = 1.22, 95% CI 1.03-1.45, P = 0.02, I(2) = 48%). In conclusion, this meta-analysis found a significant association between NUDT10 rs5945572 polymorphism and prostate cancer.

A genome-wide association study of prostate cancer in West African men

Age-adjusted mortality rates for prostate cancer are higher for African-American men compared with those of European ancestry. Recent data suggest that West African men also have elevated risk for prostate cancer relative to European men. Genetic susceptibility to prostate cancer could account for part of this difference. We conducted a genome-wide association study (GWAS) of prostate cancer in West African men in the Ghana Prostate Study. Association testing was performed using multivariable logistic regression adjusted for age and genetic ancestry for 474 prostate cancer cases and 458 population-based controls on the Illumina HumanOmni-5 Quad BeadChip. The most promising association was at 10p14 within an intron of a long non-coding RNA (lncRNA RP11-543F8.2) 360 kb centromeric of GATA3 (p = 1.29E-7). In sub-analyses, SNPs at 5q31.3 were associated with high Gleason score (≥7) cancers, the strongest of which was a missense SNP in PCDHA1 (rs34575154, p = 3.66E-8), and SNPs at Xq28 (rs985081, p = 8.66E-9) and 6q21 (rs2185710, p = 5.95E-8) were associated with low Gleason score (<7) cancers. We sought to validate our findings in silico in the African Ancestry Prostate Cancer GWAS Consortium, but only one SNP, at 10p14, replicated at p < 0.05. Of the 90 prostate cancer loci reported from studies of men of European, Asian or African-American ancestry, we were able to test 81 in the Ghana Prostate Study, and 10 of these replicated at p < 0.05. Further genetic studies of prostate cancer in West African men are needed to confirm our promising susceptibility loci.

Variants at IRX4 as prostate cancer expression quantitative trait loci

Genome-wide association studies (GWAS) have identified numerous prostate cancer-associated risk loci. Some variants at these loci may be regulatory and influence expression of nearby genes. Such loci are known as cis-expression quantitative trait loci (cis-eQTL). As cis-eQTLs are highly tissue-specific, we asked if GWAS-identified prostate cancer risk loci are cis-eQTLs in human prostate tumor tissues. We investigated 50 prostate cancer samples for their genotype at 59 prostate cancer risk-associated single-nucleotide polymorphisms (SNPs) and performed cis-eQTL analysis of transcripts from paired primary tumors within two megabase windows. We tested 586 transcript-genotype associations, of which 27 were significant (false discovery rate ≤10%). An equivalent eQTL analysis of the same prostate cancer risk loci in lymphoblastoid cell lines did not result in any significant associations. The top-ranked cis-eQTL involved the IRX4 (Iroquois homeobox protein 4) transcript and rs12653946, tagged by rs10866528 in our study (P=4.91 × 10(-5)). Replication studies, linkage disequilibrium, and imputation analyses highlight population specificity at this locus. We independently validated IRX4 as a potential prostate cancer risk gene through cis-eQTL analysis of prostate cancer risk variants. Cis-eQTL analysis in relevant tissues, even with a small sample size, can be a powerful method to expedite functional follow-up of GWAS.

Characterization of SNPs associated with prostate cancer in men of Ashkenazic descent from the set of GWAS identified SNPs: impact of cancer family history and cumulative SNP risk prediction

BACKGROUND

Genome-wide association studies (GWAS) have identified multiple SNPs associated with prostate cancer (PrCa). Population isolates may have different sets of risk alleles for PrCa constituting unique population and individual risk profiles.

METHODS

To test this hypothesis, associations between 31 GWAS SNPs of PrCa were examined among 979 PrCa cases and 1,251 controls of Ashkenazic descent using logistic regression. We also investigated risks by age at diagnosis, pathological features of PrCa, and family history of cancer. Moreover, we examined associations between cumulative number of risk alleles and PrCa and assessed the utility of risk alleles in PrCa risk prediction by comparing the area under the curve (AUC) for different logistic models.

RESULTS

Of the 31 genotyped SNPs, 8 were associated with PrCa at p ≤ 0.002 (corrected p-value threshold) with odds ratios (ORs) ranging from 1.22 to 1.42 per risk allele. Four SNPs were associated with aggressive PrCa, while three other SNPs showed potential interactions for PrCa by family history of PrCa (rs8102476; 19q13), lung cancer (rs17021918; 4q22), and breast cancer (rs10896449; 11q13). Men in the highest vs. lowest quartile of cumulative number of risk alleles had ORs of 3.70 (95% CI 2.76-4.97); 3.76 (95% CI 2.57-5.50), and 5.20 (95% CI 2.94-9.19) for overall PrCa, aggressive cancer and younger age at diagnosis, respectively. The addition of cumulative risk alleles to the model containing age at diagnosis and family history of PrCa yielded a slightly higher AUC (0.69 vs. 0.64).

CONCLUSION

These data define a set of risk alleles associated with PrCa in men of Ashkenazic descent and indicate possible genetic differences for PrCa between populations of European and Ashkenazic ancestry. Use of genetic markers might provide an opportunity to identify men at highest risk for younger age of onset PrCa; however, their clinical utility in identifying men at highest risk for aggressive cancer remains limited.