Integration of a Cross-Ancestry Polygenic Model With Clinical Risk Factors Improves Breast Cancer Risk Stratification.

PURPOSE: To develop and validate a cross-ancestry integrated risk score (caIRS) that combines a cross-ancestry polygenic risk score (caPRS) with a clinical estimator for breast cancer (BC) risk. We hypothesized that the caIRS is a better predictor of BC risk than clinical risk factors across diverse ancestry groups. METHODS: We used diverse retrospective cohort data with longitudinal follow-up to develop a caPRS and integrate it with the Tyrer-Cuzick (T-C) clinical model. We tested the association between the caIRS and BC risk in two validation cohorts including > 130,000 women. We compared model discrimination for 5-year and remaining lifetime BC risk between the caIRS and T-C and assessed how the caIRS would affect screening in the clinic. RESULTS: The caIRS outperformed T-C alone for all populations tested in both validation cohorts and contributed significantly to risk prediction beyond T-C. The area under the receiver operating characteristic curve improved from 0.57 to 0.65, and the odds ratio per standard deviation increased from 1.35 (95% CI, 1.27 to 1.43) to 1.79 (95% CI, 1.70 to 1.88) in validation cohort 1 with similar improvements observed in validation cohort 2. We observed the largest gain in positive predictive value using the caIRS in Black/African American women across both validation cohorts, with an approximately two-fold increase and an equivalent negative predictive value as the T-C. In a multivariate, age-adjusted logistic regression model including both caIRS and T-C, caIRS remained significant, indicating that caIRS provides information over T-C alone. CONCLUSION: Adding a caPRS to the T-C model improves BC risk stratification for women of multiple ancestries, which could have implications for screening recommendations and prevention.

Addressing health disparities in Hispanic breast cancer: accurate and inexpensive sequencing of BRCA1 and BRCA2.

BACKGROUND: Germline mutations in the BRCA1 and BRCA2 genes account for 20-25 % of inherited breast cancers and about 10 % of all breast cancer cases. Detection of BRCA mutation carriers can lead to therapeutic interventions such as mastectomy, oophorectomy, hormonal prevention therapy, improved screening, and targeted therapies such as PARP-inhibition. We estimate that African Americans and Hispanics are 4-5 times less likely to receive BRCA screening, despite having similar mutation frequencies as non-Jewish Caucasians, who have higher breast cancer mortality. To begin addressing this health disparity, we initiated a nationwide trial of BRCA testing of Latin American women with breast cancer. Patients were recruited through community organizations, clinics, public events, and by mail and Internet. Subjects completed the consent process and questionnaire, and provided a saliva sample by mail or in person. DNA from 120 subjects was used to sequence the entirety of BRCA1 and BRCA2 coding regions and splice sites, and validate pathogenic mutations, with a total material cost of $85/subject. Subjects ranged in age from 23 to 81 years (mean age, 51 years), 6 % had bilateral disease, 57 % were ER/PR+, 23 % HER2+, and 17 % had triple-negative disease. RESULTS: A total of seven different predicted deleterious mutations were identified, one newly described and the rest rare. In addition, four variants of unknown effect were found. CONCLUSIONS: Application of this strategy on a larger scale could lead to improved cancer care of minority and underserved populations.

Genome Analysis of Latin American Cervical Cancer: Frequent Activation of the PIK3CA Pathway.

PURPOSE: Cervical cancer is one of the most common causes of cancer mortality for women living in poverty, causing more than 28,000 deaths annually in Latin America and 266,000 worldwide. To better understand the molecular basis of the disease, we ascertained blood and tumor samples from Guatemala and Venezuela and performed genomic characterization. EXPERIMENTAL DESIGN: We performed human papillomavirus (HPV) typing and identified somatically mutated genes using exome and ultra-deep targeted sequencing with confirmation in samples from Mexico. Copy number changes were also assessed in the exome sequence. RESULTS: Cervical cancer cases in Guatemala and Venezuela have an average age of diagnosis of 50 years and 5.6 children. Analysis of 675 tumors revealed activation of PIK3CA and other PI3K/AKT pathway genes in 31% of squamous carcinomas and 24% of adeno- and adenosquamous tumors, predominantly at two sites (E542K, E545K) in the helical domain of the PIK3CA gene. This distribution of PIK3CA mutations is distinct from most other cancer types and does not result in the in vitro phosphorylation of AKT. Somatic mutations were more frequent in squamous carcinomas diagnosed after the age of 50 years. Frequent gain of chromosome 3q was found, and low PIK3CA mutation fractions in many tumors suggest that PI3K mutation can be a late event in tumor progression. CONCLUSIONS: PI3K pathway mutation is important to cervical carcinogenesis in Latin America. Therapeutic agents that directly target PI3K could play a role in the therapy of this common malignancy.

Concurrent alterations in TERT, KDM6A, and the BRCA pathway in bladder cancer.

PURPOSE: Genetic analysis of bladder cancer has revealed a number of frequently altered genes, including frequent alterations of the telomerase (TERT) gene promoter, although few altered genes have been functionally evaluated. Our objective is to characterize alterations observed by exome sequencing and sequencing of the TERT promoter, and to examine the functional relevance of histone lysine (K)-specific demethylase 6A (KDM6A/UTX), a frequently mutated histone demethylase, in bladder cancer. EXPERIMENTAL DESIGN: We analyzed bladder cancer samples from 54 U.S. patients by exome and targeted sequencing and confirmed somatic variants using normal tissue from the same patient. We examined the biologic function of KDM6A using in vivo and in vitro assays. RESULTS: We observed frequent somatic alterations in BRCA1 associated protein-1 (BAP1) in 15% of tumors, including deleterious alterations to the deubiquitinase active site and the nuclear localization signal. BAP1 mutations contribute to a high frequency of tumors with breast cancer (BRCA) DNA repair pathway alterations and were significantly associated with papillary histologic features in tumors. BAP1 and KDM6A mutations significantly co-occurred in tumors. Somatic variants altering the TERT promoter were found in 69% of tumors but were not correlated with alterations in other bladder cancer genes. We examined the function of KDM6A, altered in 24% of tumors, and show depletion in human bladder cancer cells, enhanced in vitro proliferation, in vivo tumor growth, and cell migration. CONCLUSIONS: This study is the first to identify frequent BAP1 and BRCA pathway alterations in bladder cancer, show TERT promoter alterations are independent of other bladder cancer gene alterations, and show KDM6A loss is a driver of the bladder cancer phenotype.

The new sequencer on the block: comparison of Life Technology's Proton sequencer to an Illumina HiSeq for whole-exome sequencing.

We assessed the performance of the new Life Technologies Proton sequencer by comparing whole-exome sequence data in a Centre d'Etude du Polymorphisme Humain trio (family 1463) to the Illumina HiSeq instrument. To simulate a typical user's results, we utilized the standard capture, alignment and variant calling methods specific to each platform. We restricted data analysis to include the capture region common to both methods. The Proton produced high quality data at a comparable average depth and read length, and the Ion Reporter variant caller identified 96 % of single nucleotide polymorphisms (SNPs) detected by the HiSeq and GATK pipeline. However, only 40 % of small insertion and deletion variants (indels) were identified by both methods. Usage of the trio structure and segregation of platform-specific alleles supported this result. Further comparison of the trio data with Complete Genomics sequence data and Illumina SNP microarray genotypes documented high concordance and accurate SNP genotyping of both Proton and Illumina platforms. However, our study underscored the problem of accurate detection of indels for both the Proton and HiSeq platforms.

Somatic alterations contributing to metastasis of a castration-resistant prostate cancer.

Metastatic castration-resistant prostate cancer (mCRPC) is a lethal disease, and molecular markers that differentiate indolent from aggressive subtypes are needed. We sequenced the exomes of five metastatic tumors and healthy kidney tissue from an index case with mCRPC to identify lesions associated with disease progression and metastasis. An Ashkenazi Jewish (AJ) germline founder mutation, del185AG in BRCA1, was observed and AJ ancestry was confirmed. Sixty-two somatic variants altered proteins in tumors, including cancer-associated genes, TMPRSS2-ERG, PBRM1, and TET2. The majority (n = 53) of somatic variants were present in all metastases and only a subset (n = 31) was observed in the primary tumor. Integrating tumor next-generation sequencing and DNA copy number showed somatic loss of BRCA1 and TMPRSS2-ERG. We sequenced 19 genes with deleterious mutations in the index case in additional mCRPC samples and detected a frameshift, two somatic missense alterations, tumor loss of heterozygosity, and combinations of germline missense SNPs in TET2. In summary, genetic analysis of metastases from an index case permitted us to infer a chronology for the clonal spread of disease based on sequential accrual of somatic lesions. The role of TET2 in mCRPC deserves additional analysis and may define a subset of metastatic disease.

Analysis of the ABCA4 gene by next-generation sequencing.

PURPOSE: To find all possible disease-associated variants in coding sequences of the ABCA4 gene in a large cohort of patients diagnosed with ABCA4-associated diseases. METHODS: One hundred sixty-eight patients who had been clinically diagnosed with Stargardt disease, cone-rod dystrophy, and other ABCA4-associated phenotypes were prescreened for mutations in ABCA4 with the ABCA4 microarray, resulting in finding 1 of 2 expected mutations in 111 patients and 0 of 2 mutations in 57 patients. The next-generation sequencing (NGS) strategy was applied to these patients to sequence the entire coding region and the splice sites of the ABCA4 gene. Identified new variants were confirmed or rejected by Sanger sequencing and analyzed for possible pathogenicity by in silico programs and, where possible, by segregation analyses. RESULTS: Sequencing was successful in 159 of 168 patients and identified the second disease-associated allele in 49 of 103 (~48%) of patients with one previously identified mutation. Among those with no mutations, both disease-associated alleles were detected in 4 of 56 patients, and one mutation was detected in 10 of 56 patients. The authors detected a total of 57 previously unknown, possibly pathogenic, variants: 29 missense, 4 nonsense, 9 small deletions and 15 splice-site-altering variants. Of these, 55 variants were deemed pathogenic by a combination of predictive methods and segregation analyses. CONCLUSIONS: Many mutations in the coding sequences of the ABCA4 gene are still unknown, and many possibly reside in noncoding regions of the ABCA4 locus. Although the ABCA4 array remains a good first-pass screening option, the NGS platform is a time- and cost-efficient tool for screening large cohorts.

Association assessment of copy number polymorphism and risk of age-related macular degeneration.

PURPOSE: We previously identified a genetic copy number polymorphism (CNP147) that was statistically associated with age-related macular degeneration (AMD) and that resides downstream of the complement factor H (CFH) gene. Factor H protein is polymorphic at amino acid 402, in which the resulting histidine containing moiety has been established to impart significant risk of AMD. We present a method to precisely determine the exact copy number of CNP147 and examine in more detail the association with AMD. DESIGN: Case-control study. PARTICIPANTS: A total of 421 Age-Related Eye Disease Study (AREDS) subjects, of whom approximately 35% were diagnosed with neovascular disease, 19% were diagnosed with geographic atrophy, 16% were diagnosed with both, 30% were diagnosed with large drusen, and 215 were controls. METHODS: By using copy number assays available from Applied Biosystems Inc. (Carlsbad, CA), we examined 4 loci spanning CNP147 and neighboring CNP148 in an AREDS matched case-control sample set. We analyzed these data by copy number while controlling for 2 high-risk CFH variants, rs1061170 (Y402H) and rs1410996. We phased the high-risk CFH variants with CNP147 and analyzed haplotype frequencies in cases and controls. To further validate copy numbers, 6 Utah Centre D'etude du Polymorphism Humaine (CEPH) families were typed for CNP147, and the segregation was assessed. MAIN OUTCOME MEASURES: Increased or decreased risk of AMD from genetic loci. RESULTS: Having fewer than 2 copies of CNP147 was associated with an estimated 43% reduction in odds of having AMD in this sample set (adjusted odds ratio [OR] = 0.57, P=0.006). CNP148 variation is rare in Caucasians and was not statistically significant. Common haplotypes reveal that the risk alleles for rs1061170 and rs1410996 most frequently segregate with higher copy numbers for CNP147, but not exclusively, and that 1 haplotype that carried a deletion of CNP147 was highly protective (OR = 0.25 P=1.3×10(-13)) when compared with the reference. CONCLUSIONS: In this matched subset of AREDS subjects, after adjusting for 2 known risk variants in CFH, CNP147 deletion statistically associates with diminished risk for AMD. FINANCIAL DISCLOSURE(S): The author(s) have no proprietary or commercial interest in any materials discussed in this article.

Haplotype structure in Ashkenazi Jewish BRCA1 and BRCA2 mutation carriers.

Three founder mutations in BRCA1 and BRCA2 contribute to the risk of hereditary breast and ovarian cancer in Ashkenazi Jews (AJ). They are observed at increased frequency in the AJ compared to other BRCA mutations in Caucasian non-Jews (CNJ). Several authors have proposed that elevated allele frequencies in the surrounding genomic regions reflect adaptive or balancing selection. Such proposals predict long-range linkage disequilibrium (LD) resulting from a selective sweep, although genetic drift in a founder population may also act to create long-distance LD. To date, few studies have used the tools of statistical genomics to examine the likelihood of long-range LD at a deleterious locus in a population that faced a genetic bottleneck. We studied the genotypes of hundreds of women from a large international consortium of BRCA1 and BRCA2 mutation carriers and found that AJ women exhibited long-range haplotypes compared to CNJ women. More than 50% of the AJ chromosomes with the BRCA1 185delAG mutation share an identical 2.1 Mb haplotype and nearly 16% of AJ chromosomes carrying the BRCA2 6174delT mutation share a 1.4 Mb haplotype. Simulations based on the best inference of Ashkenazi population demography indicate that long-range haplotypes are expected in the context of a genome-wide survey. Our results are consistent with the hypothesis that a local bottleneck effect from population size constriction events could by chance have resulted in the large haplotype blocks observed at high frequency in the BRCA1 and BRCA2 regions of Ashkenazi Jews.