Daniel J. Serie

Loss of BAP1 Protein Expression Is an Independent Marker of Poor Prognosis in patients with Low Risk Clear Cell Renal Cell Carcinoma

The majority of patients diagnosed with clear cell renal cell carcinoma (ccRCC) have low‐risk disease with a < 10% chance of ccRCC‐specific death. DNA sequencing revealed that mutations in BAP1 (BRCA1 associated protein‐1) occur in 5% to 15% of ccRCC cases and are associated with poor outcomes. The vast majority of BAP1 mutations abolish protein expression. In this study, we used a highly sensitive and specific immunohistochemistry (IHC) assay to test whether BAP1 expression is an independent marker of ccRCC‐specific survival, particularly in patients with low‐risk disease.

Genome-wide association study of corticobasal degeneration identifies risk variants shared with progressive supranuclear palsy

Corticobasal degeneration (CBD) is a neurodegenerative disorder affecting movement and cognition, definitively diagnosed only at autopsy. Here, we conduct a genome-wide association study (GWAS) in CBD cases (n=152) and 3,311 controls, and 67 CBD cases and 439 controls in a replication stage. Associations with meta-analysis were 17q21 at MAPT (P=1.42 × 10−12), 8p12 at lnc-KIF13B-1, a long non-coding RNA (rs643472; P=3.41 × 10−8), and 2p22 at SOS1 (rs963731; P=1.76 × 10−7). Testing for association of CBD with top progressive supranuclear palsy (PSP) GWAS single-nucleotide polymorphisms (SNPs) identified associations at MOBP (3p22; rs1768208; P=2.07 × 10−7) and MAPT H1c (17q21; rs242557; P=7.91 × 10−6). We previously reported SNP/transcript level associations with rs8070723/MAPT, rs242557/MAPT, and rs1768208/MOBP and herein identified association with rs963731/SOS1. We identify new CBD susceptibility loci and show that CBD and PSP share a genetic risk factor other than MAPT at 3p22 MOBP (myelin-associated oligodendrocyte basic protein).

A Novel Automated Mammographic Density Measure and Breast Cancer Risk

BACKGROUND Mammographic breast density is a strong breast cancer risk factor but is not used in the clinical setting, partly because of a lack of standardization and automation. We developed an automated and objective measurement of the grayscale value variation within a mammogram, evaluated its association with breast cancer, and compared its performance with that of percent density (PD). METHODS Three clinic-based studies were included: a case-cohort study of 217 breast cancer case subjects and 2094 non-case subjects and two case-control studies comprising 928 case subjects and 1039 control subjects and 246 case subjects and 516 control subjects, respectively. Percent density was estimated from digitized mammograms using the computer-assisted Cumulus thresholding program, and variation was estimated from an automated algorithm. We estimated hazards ratios (HRs), odds ratios (ORs), the area under the receiver operating characteristic curve (AUC), and 95% confidence intervals (CIs) using Cox proportional hazards models for the cohort and logistic regression for case-control studies, with adjustment for age and body mass index. We performed a meta-analysis using random study effects to obtain pooled estimates of the associations between the two mammographic measures and breast cancer. All statistical tests were two-sided. RESULTS The variation measure was statistically significantly associated with the risk of breast cancer in all three studies (highest vs lowest quartile: HR = 2.0 [95% CI = 1.3 to 3.1]; OR = 2.7 [95% CI = 2.1 to 3.6]; OR = 2.4 [95% CI = 1.4 to 3.9]; [corrected] all P (trend) < .001). [corrected]. The risk estimates and AUCs for the variation measure were similar to [corrected] those for percent density (AUCs for variation = 0.60-0.62 and [corrected] AUCs for percent density = 0.61-0.65). [corrected]. A meta-analysis of the three studies demonstrated similar associations [corrected] between variation and breast cancer (highest vs lowest quartile: RR = 1.8, 95% CI = 1.4 to 2.3) and [corrected] percent density and breast cancer (highest vs lowest quartile: RR = 2.3, 95% CI = 1.9 to 2.9). CONCLUSION The association between the automated variation measure and the risk of breast cancer is at least as strong as that for percent density. Efforts to further evaluate and translate the variation measure to the clinical setting are warranted.

Clear Cell Renal Cell Carcinoma Subtypes Identified by BAP1 and PBRM1 Expression.

PURPOSE In clear cell renal cell carcinoma BAP1 and PBRM1 are 2 of the most commonly mutated genes (10% to 15% and 40% to 50%, respectively). We sought to determine the prognostic significance of PBRM1 and BAP1 expression in clear cell renal cell carcinoma. MATERIALS AND METHODS We used immunohistochemistry to assess PBRM1 protein expression in 1,479 primary clear cell renal cell carcinoma tumors that were previously stained for BAP1. A centralized pathologist reviewed all cases and categorized tumors as positive or deficient for PBRM1 and BAP1. Kaplan-Meier and Cox regression models were used to evaluate association of PBRM1 and BAP1 expression with the risk of death from renal cell carcinoma and the risk of metastasis after adjustment for age and the Mayo Clinic SSIGN (stage, size, grade and necrosis) score. RESULTS PBRM1 and BAP1 expression was PBRM1+ BAP1+ in 40.1% of tumors, PBRM1- BAP1+ in 48.6%, PBRM1+ BAP1- in 8.7% and PBRM1- BAP1- in 1.8%. The incidence of PBRM1 and BAP1 loss in the same tumor was significantly lower than expected (actual 1.8% vs expected 5.3%, p <0.0001). Compared to patients with PBRM1+ BAP1+ tumors those with PBRM1- BAP1+ lesions were more likely to die of renal cell carcinoma (HR 1.39, p = 0.035), followed by those with PBRM1+ BAP1- and PBRM1- BAP1- tumors (HR 3.25 and 5.2, respectively, each p <0.001). PBRM1 and BAP1 expression did not add independent prognostic information to the SSIGN score. CONCLUSIONS PBRM1 and BAP1 expression identified 4 clinical subgroups of patients with clear cell renal cell carcinoma who had divergent clinical outcomes. The clinical value of these biomarkers will be fully realized when therapies targeting pathways downstream of PBRM1 and BAP1 are developed.

Common variation at 6p21.31 (BAK1) influences the risk of chronic lymphocytic leukemia

We performed a meta-analysis of 3 genome-wide association studies to identify additional common variants influencing chronic lymphocytic leukemia (CLL) risk. The discovery phase was composed of genome-wide association study data from 1121 cases and 3745 controls. Replication analysis was performed in 861 cases and 2033 controls. We identified a novel CLL risk locus at 6p21.33 (rs210142; intronic to the BAK1 gene, BCL2 antagonist killer 1; P = 9.47 × 10(-16)). A strong relationship between risk genotype and reduced BAK1 expression was shown in lymphoblastoid cell lines. This finding provides additional support for polygenic inheritance to CLL and provides further insight into the biologic basis of disease development.

The influence of mammogram acquisition on the mammographic density and breast cancer association in the mayo mammography health study cohort

IntroductionMammographic density is a strong risk factor for breast cancer. Image acquisition technique varies across mammograms to limit radiation and produce a clinically useful image. We examined whether acquisition technique parameters at the time of mammography were associated with mammographic density and whether the acquisition parameters confounded the density and breast cancer association.MethodsWe examined this question within the Mayo Mammography Health Study (MMHS) cohort, comprised of 19,924 women (51.2% of eligible) seen in the Mayo Clinic mammography screening practice from 2003 to 2006. A case-cohort design, comprising 318 incident breast cancers diagnosed through December 2009 and a random subcohort of 2,259, was used to examine potential confounding of mammogram acquisition technique parameters (x-ray tube voltage peak (kVp), milliampere-seconds (mAs), thickness and compression force) on the density and breast cancer association. The Breast Imaging Reporting and Data System four-category tissue composition measure (BI-RADS) and percent density (PD) (Cumulus program) were estimated from screen-film mammograms at time of enrollment. Spearman correlation coefficients (r) and means (standard deviations) were used to examine the relationship of density measures with acquisition parameters. Hazard ratios (HR) and C-statistics were estimated using Cox proportional hazards regression, adjusting for age, menopausal status, body mass index and postmenopausal hormones. A change in the HR of at least 15% indicated confounding.ResultsAdjusted PD and BI-RADS density were associated with breast cancer (p-trends < 0.001), with a 3 to 4-fold increased risk in the extremely dense vs. fatty BI-RADS categories (HR: 3.0, 95% CI, 1.7 - 5.1) and the ≥ 25% vs. ≤ 5% PD categories (HR: 3.8, 95% CI, 2.5 - 5.9). Of the acquisition parameters, kVp was not correlated with PD (r = 0.04, p = 0.07). Although thickness (r = -0.27, p < 0.001), compression force (r = -0.16, p < 0.001), and mAs (r = -0.06, p = 0.008) were inversely correlated with PD, they did not confound the PD or BI-RADS associations with breast cancer and their inclusion did not improve discriminatory accuracy. Results were similar for associations of dense and non-dense area with breast cancer.ConclusionsWe confirmed a strong association between mammographic density and breast cancer risk that was not confounded by mammogram acquisition technique.

Gene expression, single nucleotide variant and fusion transcript discovery in archival material from breast tumors.

Advantages of RNA-Seq over array based platforms are quantitative gene expression and discovery of expressed single nucleotide variants (eSNVs) and fusion transcripts from a single platform, but the sensitivity for each of these characteristics is unknown. We measured gene expression in a set of manually degraded RNAs, nine pairs of matched fresh-frozen, and FFPE RNA isolated from breast tumor with the hybridization based, NanoString nCounter (226 gene panel) and with whole transcriptome RNA-Seq using RiboZeroGold ScriptSeq V2 library preparation kits. We performed correlation analyses of gene expression between samples and across platforms. We then specifically assessed whole transcriptome expression of lincRNA and discovery of eSNVs and fusion transcripts in the FFPE RNA-Seq data. For gene expression in the manually degraded samples, we observed Pearson correlations of >0.94 and >0.80 with NanoString and ScriptSeq protocols, respectively. Gene expression data for matched fresh-frozen and FFPE samples yielded mean Pearson correlations of 0.874 and 0.783 for NanoString (226 genes) and ScriptSeq whole transcriptome protocols respectively, p<2x10-16. Specifically for lincRNAs, we observed superb Pearson correlation (0.988) between matched fresh-frozen and FFPE pairs. FFPE samples across NanoString and RNA-Seq platforms gave a mean Pearson correlation of 0.838. In FFPE libraries, we detected 53.4% of high confidence SNVs and 24% of high confidence fusion transcripts. Sensitivity of fusion transcript detection was not overcome by an increase in depth of sequencing up to 3-fold (increase from ~56 to ~159 million reads). Both NanoString and ScriptSeq RNA-Seq technologies yield reliable gene expression data for degraded and FFPE material. The high degree of correlation between NanoString and RNA-Seq platforms suggests discovery based whole transcriptome studies from FFPE material will produce reliable expression data. The RiboZeroGold ScriptSeq protocol performed particularly well for lincRNA expression from FFPE libraries, but detection of eSNV and fusion transcripts was less sensitive.

Identification of Novel Variants in Colorectal Cancer Families by High-Throughput Exome Sequencing

Background: Colorectal cancer (CRC) in densely affected families without Lynch Syndrome may be due to mutations in undiscovered genetic loci. Familial linkage analyses have yielded disparate results; the use of exome sequencing in coding regions may identify novel segregating variants. Methods: We completed exome sequencing on 40 affected cases from 16 multicase pedigrees to identify novel loci. Variants shared among all sequenced cases within each family were identified and filtered to exclude common variants and single-nucleotide variants (SNV) predicted to be benign. Results: We identified 32 nonsense or splice-site SNVs, 375 missense SNVs, 1,394 synonymous or noncoding SNVs, and 50 indels in the 16 families. Of particular interest are two validated and replicated missense variants in CENPE and KIF23, which are both located within previously reported CRC linkage regions, on chromosomes 1 and 15, respectively. Conclusions: Whole-exome sequencing identified DNA variants in multiple genes. Additional sequencing of these genes in additional samples will further elucidate the role of variants in these regions in CRC susceptibility. Impact: Exome sequencing of familial CRC cases can identify novel rare variants that may influence disease risk. Cancer Epidemiol Biomarkers Prev; 22(7); 1239–51. ©2013 AACR.

RNA Sequencing Identifies Multiple Fusion Transcripts, Differentially Expressed Genes, and Reduced Expression of Immune Function Genes in BRAF (V600E) Mutant vs BRAF Wild-Type Papillary Thyroid Carcinoma

CONTEXT The BRAF V600E mutation (BRAF-MUT) confers an aggressive phenotype in papillary thyroid carcinoma, but unidentified additional genomic abnormalities may be required for full phenotypic expression. OBJECTIVE RNA sequencing (RNA-Seq) was performed to identify genes differentially expressed between BRAF-MUT and BRAF wild-type (BRAF-WT) tumors and to correlate changes to patient clinical status. DESIGN BRAF-MUT and BRAF-WT tumors were identified in patients with T1N0 and T2-3N1 tumors evaluated in a referral medical center. Gene expression levels were determined (RNA-Seq) and fusion transcripts were detected. Multiplexed capture/detection and digital counting of mRNA transcripts (nCounter, NanoString Technologies) validated RNA-Seq data for immune system-related genes. PATIENTS BRAF-MUT patients included nine women, three men; nine were TNM stage I and three were stage III. Three (25%) had tumor infiltrating lymphocytes. BRAF-WT included five women, three men; all were stage I, and five (62.5%) had tumor infiltrating lymphocytes. RESULTS RNA-Seq identified 560 of 13 085 genes differentially expressed between BRAF-MUT and BRAF-WT tumors. Approximately 10% of these genes were related to MetaCore immune function pathways; 51 were underexpressed in BRAF-MUT tumors, whereas 4 (HLAG, CXCL14, TIMP1, IL1RAP) were overexpressed. The four most differentially overexpressed immune genes in BRAF-WT tumors (IL1B; CCL19; CCL21; CXCR4) correlated with lymphocyte infiltration. nCounter confirmed the RNA-Seq expression level data. Eleven different high-confidence fusion transcripts were detected (four interchromosomal; seven intrachromosomal) in 13 of 20 tumors. All in-frame fusions were validated by RT-PCR. CONCLUSION BRAF-MUT papillary thyroid cancers have reduced expression of immune/inflammatory response genes compared with BRAF-WT tumors and correlate with lymphocyte infiltration. In contrast, HLA-G and CXCL14 are overexpressed in BRAF-MUT tumors. Sixty-five percent of tumors had between one and three fusion transcripts. Functional studies will be required to determine the potential role of these newly identified genomic abnormalities in contributing to the aggressiveness of BRAF-MUT and BRAF-WT tumors.

Evaluation of memory endophenotypes for association with CLU, CR1, and PICALM variants in black and white subjects

Genetic variants at the CLU, CR1, and PICALM loci associate with risk for late‐onset Alzheimers disease (LOAD) in genomewide association studies. In this study, our aim was to determine whether the LOAD risk variants at these three loci influence memory endophenotypes in black and white subjects.