BreastPRS is a gene expression assay that stratifies intermediate-risk Oncotype DX patients into high- or low-risk for disease recurrence.

Molecular prognostic assays, such as Oncotype DX, are increasingly incorporated into the management of patients with invasive breast carcinoma. BreastPRS is a new molecular assay developed and validated from a meta-analysis of publically available genomic datasets. We applied the assay to matched fresh-frozen (FF) and formalin-fixed paraffin-embedded (FFPE) tumor samples to translate the assay to FFPE. A linear relationship of the BreastPRS prognostic score was observed between tissue preservation formats. BreastPRS recurrence scores were compared with Oncotype DX recurrence scores from 246 patients with invasive breast carcinoma and known Oncotype DX results. Using this series, a 120-gene Oncotype DX approximation algorithm was trained to predict Oncotype DX risk groups and then applied to series of untreated, node-negative, estrogen receptor (ER)-positive patients from previously published studies with known clinical outcomes. Correlation of recurrence score and risk group between Oncotype DX and BreastPRS was statistically significant (P < 0.0001). 59 of 260 (23 %) patients from four previously published studies were classified as intermediate-risk when the 120-gene Oncotype DX approximation algorithm was applied. BreastPRS reclassified the 59 patients into binary risk groups (high- vs. low-risk). 23 (39 %) patients were classified as low-risk and 36 (61 %) as high-risk (P = 0.029, HR: 3.64, 95 % CI: 1.40-9.50). At 10 years from diagnosis, the low-risk group had a 90 % recurrence-free survival (RFS) rate compared to 60 % for the high-risk group. BreastPRS recurrence score is comparable with Oncotype DX and can reclassify Oncotype DX intermediate-risk patients into two groups with significant differences in RFS. Further studies are needed to validate these findings.

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.

Exome sequencing identifies recurrent SPOP, FOXA1 and MED12 mutations in prostate cancer.

Prostate cancer is the second most common cancer in men worldwide and causes over 250,000 deaths each year. Overtreatment of indolent disease also results in significant morbidity. Common genetic alterations in prostate cancer include losses of NKX3.1 (8p21) and PTEN (10q23), gains of AR (the androgen receptor gene) and fusion of ETS family transcription factor genes with androgen-responsive promoters. Recurrent somatic base-pair substitutions are believed to be less contributory in prostate tumorigenesis but have not been systematically analyzed in large cohorts. Here, we sequenced the exomes of 112 prostate tumor and normal tissue pairs. New recurrent mutations were identified in multiple genes, including MED12 and FOXA1. SPOP was the most frequently mutated gene, with mutations involving the SPOP substrate-binding cleft in 6-15% of tumors across multiple independent cohorts. Prostate cancers with mutant SPOP lacked ETS family gene rearrangements and showed a distinct pattern of genomic alterations. Thus, SPOP mutations may define a new molecular subtype of prostate cancer.

A computational framework discovers new copy number variants with functional importance.

Structural variants which cause changes in copy numbers constitute an important component of genomic variability. They account for 0.7% of genomic differences in two individual genomes, of which copy number variants (CNVs) are the largest component. A recent population-based CNV study revealed the need of better characterization of CNVs, especially the small ones (<500 bp).We propose a three step computational framework (Identification of germline Changes in Copy Number or IgC2N) to discover and genotype germline CNVs. First, we detect candidate CNV loci by combining information across multiple samples without imposing restrictions to the number of coverage markers or to the variant size. Secondly, we fine tune the detection of rare variants and infer the putative copy number classes for each locus. Last, for each variant we combine the relative distance between consecutive copy number classes with genetic information in a novel attempt to estimate the reference model bias. This computational approach is applied to genome-wide data from 1250 HapMap individuals. Novel variants were discovered and characterized in terms of size, minor allele frequency, type of polymorphism (gains, losses or both), and mechanism of formation. Using data generated for a subset of individuals by a 42 million marker platform, we validated the majority of the variants with the highest validation rate (66.7%) was for variants of size larger than 1 kb. Finally, we queried transcriptomic data from 129 individuals determined by RNA-sequencing as further validation and to assess the functional role of the new variants. We investigated the possible enrichment for variant's regulatory effect and found that smaller variants (<1 Kb) are more likely to regulate gene transcript than larger variants (p-value = 2.04e-08). Our results support the validity of the computational framework to detect novel variants relevant to disease susceptibility studies and provide evidence of the importance of genetic variants in regulatory network studies.