Multi-ancestry epigenome-wide analyses identify methylated sites associated with aortic augmentation index in TOPMed MESA

Despite the prognostic value of arterial stiffness (AS) and pulsatile hemodynamics (PH) for cardiovascular morbidity and mortality, epigenetic modifications that contribute to AS/PH remain unknown. To gain a better understanding of the link between epigenetics (DNA methylation) and AS/PH, we examined the relationship of eight measures of AS/PH with CpG sites and co-methylated regions using multi-ancestry participants from Trans-Omics for Precision Medicine (TOPMed) Multi-Ethnic Study of Atherosclerosis (MESA) with sample sizes ranging from 438 to 874. Epigenome-wide association analysis identified one genome-wide significant CpG (cg20711926-CYP1B1) associated with aortic augmentation index (AIx). Follow-up analyses, including gene set enrichment analysis, expression quantitative trait methylation analysis, and functional enrichment analysis on differentially methylated positions and regions, further prioritized three CpGs and their annotated genes (cg23800023-ETS1, cg08426368-TGFB3, and cg17350632-HLA-DPB1) for AIx. Among these, ETS1 and TGFB3 have been previously prioritized as candidate genes. Furthermore, both ETS1 and HLA-DPB1 have significant tissue correlations between Whole Blood and Aorta in GTEx, which suggests ETS1 and HLA-DPB1 could be potential biomarkers in understanding pathophysiology of AS/PH. Overall, our findings support the possible role of epigenetic regulation via DNA methylation of specific genes associated with AIx as well as identifying potential targets for regulation of AS/PH.

Distinct germline genetic susceptibility profiles identified for common non-Hodgkin lymphoma subtypes

Lymphoma risk is elevated for relatives with common non-Hodgkin lymphoma (NHL) subtypes, suggesting shared genetic susceptibility across subtypes. To evaluate the extent of mutual heritability among NHL subtypes and discover novel loci shared among subtypes, we analyzed data from eight genome-wide association studies within the InterLymph Consortium, including 10,629 cases and 9505 controls. We utilized Association analysis based on SubSETs (ASSET) to discover loci for subsets of NHL subtypes and evaluated shared heritability across the genome using Genome-wide Complex Trait Analysis (GCTA) and polygenic risk scores. We discovered 17 genome-wide significant loci (P < 5 × 10-8) for subsets of NHL subtypes, including a novel locus at 10q23.33 (HHEX) (P = 3.27 × 10-9). Most subset associations were driven primarily by only one subtype. Genome-wide genetic correlations between pairs of subtypes varied broadly from 0.20 to 0.86, suggesting substantial heterogeneity in the extent of shared heritability among subtypes. Polygenic risk score analyses of established loci for different lymphoid malignancies identified strong associations with some NHL subtypes (P < 5 × 10-8), but weak or null associations with others. Although our analyses suggest partially shared heritability and biological pathways, they reveal substantial heterogeneity among NHL subtypes with each having its own distinct germline genetic architecture.

Pulmonary Function and Blood DNA Methylation: A Multiancestry Epigenome-Wide Association Meta-analysis

Rationale: Methylation integrates factors present at birth and modifiable across the lifespan that can influence pulmonary function. Studies are limited in scope and replication. Objectives: To conduct large-scale epigenome-wide meta-analyses of blood DNA methylation and pulmonary function. Methods: Twelve cohorts analyzed associations of methylation at cytosine-phosphate-guanine probes (CpGs), using Illumina 450K or EPIC/850K arrays, with FEV₁, FVC, and FEV₁/FVC. We performed multiancestry epigenome-wide meta-analyses (total of 17,503 individuals; 14,761 European, 2,549 African, and 193 Hispanic/Latino ancestries) and interpreted results using integrative epigenomics. Measurements and Main Results: We identified 1,267 CpGs (1,042 genes) differentially methylated (false discovery rate, <0.025) in relation to FEV₁, FVC, or FEV₁/FVC, including 1,240 novel and 73 also related to chronic obstructive pulmonary disease (1,787 cases). We found 294 CpGs unique to European or African ancestry and 395 CpGs unique to never or ever smokers. The majority of significant CpGs correlated with nearby gene expression in blood. Findings were enriched in key regulatory elements for gene function, including accessible chromatin elements, in both blood and lung. Sixty-nine implicated genes are targets of investigational or approved drugs. One example novel gene highlighted by integrative epigenomic and druggable target analysis is TNFRSF4. Mendelian randomization and colocalization analyses suggest that epigenome-wide association study signals capture causal regulatory genomic loci. Conclusions: We identified numerous novel loci differentially methylated in relation to pulmonary function; few were detected in large genome-wide association studies. Integrative analyses highlight functional relevance and potential therapeutic targets. This comprehensive discovery of potentially modifiable, novel lung function loci expands knowledge gained from genetic studies, providing insights into lung pathogenesis.

Epigenome-wide association study of mitochondrial genome copy number

We conducted cohort- and race-specific epigenome-wide association analyses of mitochondrial deoxyribonucleic acid (mtDNA) copy number (mtDNA CN) measured in whole blood from participants of African and European origins in five cohorts (n = 6182, mean age = 57-67 years, 65% women). In the meta-analysis of all the participants, we discovered 21 mtDNA CN-associated DNA methylation sites (CpG) (P < 1 × 10-7), with a 0.7-3.0 standard deviation increase (3 CpGs) or decrease (18 CpGs) in mtDNA CN corresponding to a 1% increase in DNA methylation. Several significant CpGs have been reported to be associated with at least two risk factors (e.g. chronological age or smoking) for cardiovascular disease (CVD). Five genes [PR/SET domain 16, nuclear receptor subfamily 1 group H member 3 (NR1H3), DNA repair protein, DNA polymerase kappa and decaprenyl-diphosphate synthase subunit 2], which harbor nine significant CpGs, are known to be involved in mitochondrial biosynthesis and functions. For example, NR1H3 encodes a transcription factor that is differentially expressed during an adipose tissue transition. The methylation level of cg09548275 in NR1H3 was negatively associated with mtDNA CN (effect size = -1.71, P = 4 × 10-8) and was positively associated with the NR1H3 expression level (effect size = 0.43, P = 0.0003), which indicates that the methylation level in NR1H3 may underlie the relationship between mtDNA CN, the NR1H3 transcription factor and energy expenditure. In summary, the study results suggest that mtDNA CN variation in whole blood is associated with DNA methylation levels in genes that are involved in a wide range of mitochondrial activities. These findings will help reveal molecular mechanisms between mtDNA CN and CVD.

Genome-wide association studies identify 137 genetic loci for DNA methylation biomarkers of aging

BACKGROUND

Biological aging estimators derived from DNA methylation data are heritable and correlate with morbidity and mortality. Consequently, identification of genetic and environmental contributors to the variation in these measures in populations has become a major goal in the field.

RESULTS

Leveraging DNA methylation and SNP data from more than 40,000 individuals, we identify 137 genome-wide significant loci, of which 113 are novel, from genome-wide association study (GWAS) meta-analyses of four epigenetic clocks and epigenetic surrogate markers for granulocyte proportions and plasminogen activator inhibitor 1 levels, respectively. We find evidence for shared genetic loci associated with the Horvath clock and expression of transcripts encoding genes linked to lipid metabolism and immune function. Notably, these loci are independent of those reported to regulate DNA methylation levels at constituent clock CpGs. A polygenic score for GrimAge acceleration showed strong associations with adiposity-related traits, educational attainment, parental longevity, and C-reactive protein levels.

CONCLUSION

This study illuminates the genetic architecture underlying epigenetic aging and its shared genetic contributions with lifestyle factors and longevity.

Discovery and fine-mapping of height loci via high-density imputation of GWASs in individuals of African ancestry

Although many loci have been associated with height in European ancestry populations, very few have been identified in African ancestry individuals. Furthermore, many of the known loci have yet to be generalized to and fine-mapped within a large-scale African ancestry sample. We performed sex-combined and sex-stratified meta-analyses in up to 52,764 individuals with height and genome-wide genotyping data from the African Ancestry Anthropometry Genetics Consortium (AAAGC). We additionally combined our African ancestry meta-analysis results with published European genome-wide association study (GWAS) data. In the African ancestry analyses, we identified three novel loci (SLC4A3, NCOA2, ECD/FAM149B1) in sex-combined results and two loci (CRB1, KLF6) in women only. In the African plus European sex-combined GWAS, we identified an additional three novel loci (RCCD1, G6PC3, CEP95) which were equally driven by AAAGC and European results. Among 39 genome-wide significant signals at known loci, conditioning index SNPs from European studies identified 20 secondary signals. Two of the 20 new secondary signals and none of the 8 novel loci had minor allele frequencies (MAF) < 5%. Of 802 known European height signals, 643 displayed directionally consistent associations with height, of which 205 were nominally significant (p < 0.05) in the African ancestry sex-combined sample. Furthermore, 148 of 241 loci contained ≤20 variants in the credible sets that jointly account for 99% of the posterior probability of driving the associations. In summary, trans-ethnic meta-analyses revealed novel signals and further improved fine-mapping of putative causal variants in loci shared between African and European ancestry populations.

Sequencing of 53,831 diverse genomes from the NHLBI TOPMed Program

The Trans-Omics for Precision Medicine (TOPMed) programme seeks to elucidate the genetic architecture and biology of heart, lung, blood and sleep disorders, with the ultimate goal of improving diagnosis, treatment and prevention of these diseases. The initial phases of the programme focused on whole-genome sequencing of individuals with rich phenotypic data and diverse backgrounds. Here we describe the TOPMed goals and design as well as the available resources and early insights obtained from the sequence data. The resources include a variant browser, a genotype imputation server, and genomic and phenotypic data that are available through dbGaP (Database of Genotypes and Phenotypes)1. In the first 53,831 TOPMed samples, we detected more than 400 million single-nucleotide and insertion or deletion variants after alignment with the reference genome. Additional previously undescribed variants were detected through assembly of unmapped reads and customized analysis in highly variable loci. Among the more than 400 million detected variants, 97% have frequencies of less than 1% and 46% are singletons that are present in only one individual (53% among unrelated individuals). These rare variants provide insights into mutational processes and recent human evolutionary history. The extensive catalogue of genetic variation in TOPMed studies provides unique opportunities for exploring the contributions of rare and noncoding sequence variants to phenotypic variation. Furthermore, combining TOPMed haplotypes with modern imputation methods improves the power and reach of genome-wide association studies to include variants down to a frequency of approximately 0.01%.

Abstract PR05: A meta-analysis of genome-wide association study and eQTL analysis of multiple myeloma among African Americans

Background: Persons of African ancestry (AA) experience a 1.5-2-fold risk of multiple myeloma (MM) compared to persons of European ancestry (EA). We assembled a set of MM patients with self-reported AA in order to evaluate the contribution of genetics to etiology in this high-risk group.

Methods: Here we present the results of a meta-analysis of two GWAS in 1,813 cases and 8,871 controls of AA. We also conducted an admixture mapping scan to identify risk alleles associated with local ancestry, fine-mapped the 23 known susceptibility loci to find markers that could better capture MM risk in individuals of AA, and constructed a polygenic risk score (PRS) to assess the aggregated effect of known MM risk alleles. Finally, we conducted an eQTL analysis measuring gene expression in those genes harboring a risk variant in malignant plasma cells from 292 of the patients from a single site.

Results: In GWAS analysis, we identified two suggestive novel loci located at 9p24.3 and 9p13.1 at P&lt;1 × 10-6, but no genome-wide significant association was noted. In admixture mapping, we observed a genome-wide significant inverse association between local AA at 2p24.1-23.1 and MM risk in AA individuals. 20 of the 23 known EA risk variants showed directional consistency and 9 replicated at P&lt;0.05 in AA individuals. In eight regions, we identified markers that better capture MM risk in persons of AA. AA individuals with a PRS in the top 10% had a 1.82-fold (95%CI: 1.56, 2.11) increased MM risk compared to those with average risk (25-75%). The strongest functional association was between the risk allele for variant rs56219066 at 5q15 and lower ELL2 expression (P= 5.1 × 10–12).

Conclusion: Our study shows that common genetic variation contributes to MM risk individuals of AA.

This abstract is also being presented as Poster C040.

Citation Format: Zhaohui Du, Niels Weinhold, Gregory Chi Song, Kristen A. Rand, David J. Van Den Berg, Amie E. Hwang, Xin Sheng, Victor Hom, Sikander Ailawadhi, Ajay K. Nooka, Seema Singhal, Karen Pawlish, Edward Peters, Cathryn Bock, Ann Mohrbacher, Alexander Stram, Sonja I. Berndt, William J. Blot, Graham Casey, Victoria L. Stevens, Rick Kittles, Phyllis J. Goodman, W. Ryan Diver, Anselm Hennis, Barbara Nemesure, Eric A. Klein, Benjamin A. Rybicki, Janet L. Stanford, John S. Witte, Lisa Signorello, Esther M. John, Leslie Bernstein, Antoinette Stroup, Owen W. Stephens, Maurizio Zangari, Frits Van Rhee, Andrew Olshan, Wei Zheng, Jennifer J. Hu, Regina Ziegler, Sarah J. Nyante, Sue Ann Ingles, Michael Press, John David Carpten, Stephen Chanock, Jayesh Mehta, Graham A Colditz, Jeffrey Wolf, Thomas G. Martin, Michael Tomasson, Mark A. Fiala, Howard Terebelo, Nalini Janakiraman, Laurence Kolonel, Kenneth C. Anderson, Loic Le Marchand, Daniel Auclair, Brian C.-H. Chiu, Elad Ziv, Daniel Stram, Ravi Vij, Leon Bernal-Mizrachi, Gareth J. Morgan, Jeffrey A. Zonder, Carol Ann Huff, Sagar Lonial, Robert Z. Orlowski, David V. Conti, Christopher A. Haiman, Wendy Cozen. A meta-analysis of genome-wide association study and eQTL analysis of multiple myeloma among African Americans [abstract]. In: Proceedings of the Twelfth AACR Conference on the Science of Cancer Health Disparities in Racial/Ethnic Minorities and the Medically Underserved; 2019 Sep 20-23; San Francisco, CA. Philadelphia (PA): AACR; Cancer Epidemiol Biomarkers Prev 2020;29(6 Suppl_2):Abstract nr PR05.

Variability in Cytogenetic Testing for Multiple Myeloma: A Comprehensive Analysis From Across the United States

PURPOSE

Multiple myeloma (MM) treatment has changed tremendously, with significant improvement in patient out-comes. One group with a suboptimal benefit is patients with high-risk cytogenetics, as tested by conventional karyotyping or fluorescence in situ hybridization (FISH). Methodology for these tests has been published, but not necessarily standardized.

METHODS

We address variability in the testing and reporting methodology for MM cytogenetics in the United States using the ongoing African American Multiple Myeloma Study (AAMMS). We evaluated clinical and cytogenetic data from 1,221 patients (1,161 with conventional karyotyping and 976 with FISH) tested between 1998 and 2016 across 58 laboratories nationwide.

RESULTS

Interlab and intralab variability was noted for the number of cells analyzed for karyotyping, with a significantly higher number of cells analyzed in patients in whom cytogenetics were normal (P 5.0025). For FISH testing, CD138-positive cell enrichment was used in 29.7% of patients and no enrichment in 50% of patients, whereas the remainder had unknown status. A significantly smaller number of cells was analyzed for patients in which CD138 cell enrichment was used compared with those without such enrichment (median, 50 v 200; P, .0001). A median of 7 loci probes (range, 1-16) were used for FISH testing across all laboratories, with variability in the loci probed even within a given laboratory. Chromosome 13-related abnormalities were the most frequently tested abnormality (n5956; 97.9%), and t(14;16) was the least frequently tested abnormality (n 5 119; 12.2%).

CONCLUSIONS

We report significant variability in cytogenetic testing across the United States for MM, potentially leading to variability in risk stratification, with possible clinical implications and personalized treatment approaches.

Pan-cancer analysis of whole genomes

Cancer is driven by genetic change, and the advent of massively parallel sequencing has enabled systematic documentation of this variation at the whole-genome scale1-3. Here we report the integrative analysis of 2,658 whole-cancer genomes and their matching normal tissues across 38 tumour types from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). We describe the generation of the PCAWG resource, facilitated by international data sharing using compute clouds. On average, cancer genomes contained 4-5 driver mutations when combining coding and non-coding genomic elements; however, in around 5% of cases no drivers were identified, suggesting that cancer driver discovery is not yet complete. Chromothripsis, in which many clustered structural variants arise in a single catastrophic event, is frequently an early event in tumour evolution; in acral melanoma, for example, these events precede most somatic point mutations and affect several cancer-associated genes simultaneously. Cancers with abnormal telomere maintenance often originate from tissues with low replicative activity and show several mechanisms of preventing telomere attrition to critical levels. Common and rare germline variants affect patterns of somatic mutation, including point mutations, structural variants and somatic retrotransposition. A collection of papers from the PCAWG Consortium describes non-coding mutations that drive cancer beyond those in the TERT promoter4; identifies new signatures of mutational processes that cause base substitutions, small insertions and deletions and structural variation5,6; analyses timings and patterns of tumour evolution7; describes the diverse transcriptional consequences of somatic mutation on splicing, expression levels, fusion genes and promoter activity8,9; and evaluates a range of more-specialized features of cancer genomes8,10-18.

Data-adaptive multi-locus association testing in subjects with arbitrary genealogical relationships

Genome-wide sequencing enables evaluation of associations between traits and combinations of variants in genes and pathways. But such evaluation requires multi-locus association tests with good power, regardless of the variant and trait characteristics. And since analyzing families may yield more power than analyzing unrelated individuals, we need multi-locus tests applicable to both related and unrelated individuals. Here we describe such tests, and we introduce SKAT-X, a new test statistic that uses genome-wide data obtained from related or unrelated subjects to optimize power for the specific data at hand. Simulations show that: a) SKAT-X performs well regardless of variant and trait characteristics; and b) for binary traits, analyzing affected relatives brings more power than analyzing unrelated individuals, consistent with previous findings for single-locus tests. We illustrate the methods by application to rare unclassified missense variants in the tumor suppressor gene BRCA2, as applied to combined data from prostate cancer families and unrelated prostate cancer cases and controls in the Multi-ethnic Cohort (MEC). The methods can be implemented using open-source code for public use as the R-package GATARS (Genetic Association Tests for Arbitrarily Related Subjects) <https://gailg.github.io/gatars/>.

Novel Common Genetic Susceptibility Loci for Colorectal Cancer

BACKGROUND

Previous genome-wide association studies (GWAS) have identified 42 loci (P < 5 × 10-8) associated with risk of colorectal cancer (CRC). Expanded consortium efforts facilitating the discovery of additional susceptibility loci may capture unexplained familial risk.

METHODS

We conducted a GWAS in European descent CRC cases and control subjects using a discovery-replication design, followed by examination of novel findings in a multiethnic sample (cumulative n = 163 315). In the discovery stage (36 948 case subjects/30 864 control subjects), we identified genetic variants with a minor allele frequency of 1% or greater associated with risk of CRC using logistic regression followed by a fixed-effects inverse variance weighted meta-analysis. All novel independent variants reaching genome-wide statistical significance (two-sided P < 5 × 10-8) were tested for replication in separate European ancestry samples (12 952 case subjects/48 383 control subjects). Next, we examined the generalizability of discovered variants in East Asians, African Americans, and Hispanics (12 085 case subjects/22 083 control subjects). Finally, we examined the contributions of novel risk variants to familial relative risk and examined the prediction capabilities of a polygenic risk score. All statistical tests were two-sided.

RESULTS

The discovery GWAS identified 11 variants associated with CRC at P < 5 × 10-8, of which nine (at 4q22.2/5p15.33/5p13.1/6p21.31/6p12.1/10q11.23/12q24.21/16q24.1/20q13.13) independently replicated at a P value of less than .05. Multiethnic follow-up supported the generalizability of discovery findings. These results demonstrated a 14.7% increase in familial relative risk explained by common risk alleles from 10.3% (95% confidence interval [CI] = 7.9% to 13.7%; known variants) to 11.9% (95% CI = 9.2% to 15.5%; known and novel variants). A polygenic risk score identified 4.3% of the population at an odds ratio for developing CRC of at least 2.0.

CONCLUSIONS

This study provides insight into the architecture of common genetic variation contributing to CRC etiology and improves risk prediction for individualized screening.

Discovery of common and rare genetic risk variants for colorectal cancer

To further dissect the genetic architecture of colorectal cancer (CRC), we performed whole-genome sequencing of 1,439 cases and 720 controls, imputed discovered sequence variants and Haplotype Reference Consortium panel variants into genome-wide association study data, and tested for association in 34,869 cases and 29,051 controls. Findings were followed up in an additional 23,262 cases and 38,296 controls. We discovered a strongly protective 0.3% frequency variant signal at CHD1. In a combined meta-analysis of 125,478 individuals, we identified 40 new independent signals at P < 5 × 10-8, bringing the number of known independent signals for CRC to ~100. New signals implicate lower-frequency variants, Krüppel-like factors, Hedgehog signaling, Hippo-YAP signaling, long noncoding RNAs and somatic drivers, and support a role for immune function. Heritability analyses suggest that CRC risk is highly polygenic, and larger, more comprehensive studies enabling rare variant analysis will improve understanding of biology underlying this risk and influence personalized screening strategies and drug development.

The OncoArray Consortium: A Network for Understanding the Genetic Architecture of Common Cancers

BACKGROUND

Common cancers develop through a multistep process often including inherited susceptibility. Collaboration among multiple institutions, and funding from multiple sources, has allowed the development of an inexpensive genotyping microarray, the OncoArray. The array includes a genome-wide backbone, comprising 230,000 SNPs tagging most common genetic variants, together with dense mapping of known susceptibility regions, rare variants from sequencing experiments, pharmacogenetic markers, and cancer-related traits.

METHODS

The OncoArray can be genotyped using a novel technology developed by Illumina to facilitate efficient genotyping. The consortium developed standard approaches for selecting SNPs for study, for quality control of markers, and for ancestry analysis. The array was genotyped at selected sites and with prespecified replicate samples to permit evaluation of genotyping accuracy among centers and by ethnic background.

RESULTS

The OncoArray consortium genotyped 447,705 samples. A total of 494,763 SNPs passed quality control steps with a sample success rate of 97% of the samples. Participating sites performed ancestry analysis using a common set of markers and a scoring algorithm based on principal components analysis.

CONCLUSIONS

Results from these analyses will enable researchers to identify new susceptibility loci, perform fine-mapping of new or known loci associated with either single or multiple cancers, assess the degree of overlap in cancer causation and pleiotropic effects of loci that have been identified for disease-specific risk, and jointly model genetic, environmental, and lifestyle-related exposures.

IMPACT

Ongoing analyses will shed light on etiology and risk assessment for many types of cancer. Cancer Epidemiol Biomarkers Prev; 26(1); 126-35. ©2016 AACR.

A Meta-analysis of Multiple Myeloma Risk Regions in African and European Ancestry Populations Identifies Putatively Functional Loci

BACKGROUND

Genome-wide association studies (GWAS) in European populations have identified genetic risk variants associated with multiple myeloma.

METHODS

We performed association testing of common variation in eight regions in 1,318 patients with multiple myeloma and 1,480 controls of European ancestry and 1,305 patients with multiple myeloma and 7,078 controls of African ancestry and conducted a meta-analysis to localize the signals, with epigenetic annotation used to predict functionality.

RESULTS

We found that variants in 7p15.3, 17p11.2, 22q13.1 were statistically significantly (P < 0.05) associated with multiple myeloma risk in persons of African ancestry and persons of European ancestry, and the variant in 3p22.1 was associated in European ancestry only. In a combined African ancestry-European ancestry meta-analysis, variation in five regions (2p23.3, 3p22.1, 7p15.3, 17p11.2, 22q13.1) was statistically significantly associated with multiple myeloma risk. In 3p22.1, the correlated variants clustered within the gene body of ULK4 Correlated variants in 7p15.3 clustered around an enhancer at the 3' end of the CDCA7L transcription termination site. A missense variant at 17p11.2 (rs34562254, Pro251Leu, OR, 1.32; P = 2.93 × 10-7) in TNFRSF13B encodes a lymphocyte-specific protein in the TNF receptor family that interacts with the NF-κB pathway. SNPs correlated with the index signal in 22q13.1 cluster around the promoter and enhancer regions of CBX7 CONCLUSIONS: We found that reported multiple myeloma susceptibility regions contain risk variants important across populations, supporting the use of multiple racial/ethnic groups with different underlying genetic architecture to enhance the localization and identification of putatively functional alleles.

IMPACT

A subset of reported risk loci for multiple myeloma has consistent effects across populations and is likely to be functional. Cancer Epidemiol Biomarkers Prev; 25(12); 1609-18. ©2016 AACR.

Abstract 4629: Multiple myeloma susceptibility loci examined in African and European ancestry populations

Genome-wide association studies (GWAS) of multiple myeloma (MM) in Northern Europeans have identified seven novel risk loci (2p23.3, 3p22.1, 3q26.2, 6p21.33, 7p15.3, 17p11.2, 22q13.1). We performed a multiethnic meta-analysis of these regions in 1,274 MM patients and 1,486 controls of European ancestry (EA) and 1,049 MM patients and 7,080 controls of African ancestry (AA), leveraging the differential linkage-disequilibrium of these populations in order to better localize the putative functional variants. We observed directionally consistent effects for all seven index SNPs in both populations, with four significantly associated (p&lt;0.05) with risk in EAs (3p22.1, 7p15.3, 17p11.2, 22q13.1), and two significantly associated with risk in AAs (7p15.3 and 22q13.1). In a fixed effects meta-analysis of six regions (excluding the HLA region on chromosome 6), variation in five of the regions (2p33.3, 3p22.1, 7p15.3, 17p11.2, 22q13.1) had statistically significant associations with risk (Table 1). In one region, the index variant had the strongest association [rs4487645 at 7p15.3, (OR = 1.30, p = 8.7×10−8)]. Five of the six most significantly associated variants identified in the multiethnic analyses overlapped with biologically relevant features indicating regulatory activity based on CD20+ (B lymphocyte) cells, showing evidence of potential function; those included a missense variant in (17p11.2, rs34562254, Pro251Leu) in TNFRSF13B, which encodes a lymphocyte-specific protein in the tumor necrosis factor receptor family that interacts with the NF-kb pathway. Our study shows that these regions are important in MM risk across ethnicities and further supports the use of multiple ethnic groups in genetic studies to enhance identification of risk variants.

Table 1.Most significant associations for each region in the multiethnic meta-analysis.Individuals of European AncestryIndividuals of African AncestryMulitethnic Metar2 with IndexcCHRSNPRAaFreqbORP-valueFreqbORP-valueORP-valueP-het2rs732075G0.591.222.0×10−30.621.122.0×10−21.162.6×10−40.280.09/0.283rs73069394A0.191.243.0×10−30.621.181.5×10−21.201.3×10−50.550.77/0.963rs12637184G0.761.136.0×10−20.921.192.7×10−11.151.0×10−20.640.94/1.007rs4487645C0.671.237.0×10−40.891.485.5×10−51.308.7×10−80.07-d17rs34562254A0.121.452.4×10−50.131.212.2×10−31.312.5×10−60.120.33/0.9022rs139400T0.491.224.0×10−40.531.172.1×10−31.191.2×10−60.630.63/0.96aRisk allelebFrequency of the risk allele in European and African ancestry studiescr2 metrics based on 1000 Genomes Project AFR/EUR populationsdIndex SNP

Citation Format: Kristin A. Rand, Chi Song, Eric Dean, Daniel Serie, Karen Curtin, Dennis Hazelett, Amie E. Hwang, Xin Sheng, Alex Stram, David J. Van Den Berg, Carol Ann Huff, Leon Bernal-Mizrachi, Michael H. Tomasson, Sikander Ailawadhi, Anneclaire De Roos, Seema Singhal, Karen Pawlish, Edward Peters, Catherine Bock, David V. Conti, Graham Colditz, Todd Zimmerman, Scott Huntsman, John Graff, African Ancestry Prostate Cancer GWAS Consortium,African Ancestry Breast Cancer GWAS Consortium, Stephen J. Chanock, Michael Lieber, Jayesh Mehta, Eric A. Klein, Nalini Janakiraman, Richard K. Severson, Angela R. Brooks-Wilson, Vincent Rajkumar, Elizabeth E. Brown, Laurence Kolonel, Susan Slager, Brian E. Henderson, Graham G. Giles, John J. Spinelli, Brian Chiu, Kenneth C. Anderson, Jeffrey Zonder, Robert Z. Orlowski, Sagar Lonial, Nicola Camp, Celine Vachon, Elad Ziv, Dan O. Stram, Christopher A. Haiman, Wendy Cozen. Multiple myeloma susceptibility loci examined in African and European ancestry populations. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 4629. doi:10.1158/1538-7445.AM2015-4629

Comprehensive, Integrative Genomic Analysis of Diffuse Lower-Grade Gliomas

BACKGROUND

Diffuse low-grade and intermediate-grade gliomas (which together make up the lower-grade gliomas, World Health Organization grades II and III) have highly variable clinical behavior that is not adequately predicted on the basis of histologic class. Some are indolent; others quickly progress to glioblastoma. The uncertainty is compounded by interobserver variability in histologic diagnosis. Mutations in IDH, TP53, and ATRX and codeletion of chromosome arms 1p and 19q (1p/19q codeletion) have been implicated as clinically relevant markers of lower-grade gliomas.

METHODS

We performed genomewide analyses of 293 lower-grade gliomas from adults, incorporating exome sequence, DNA copy number, DNA methylation, messenger RNA expression, microRNA expression, and targeted protein expression. These data were integrated and tested for correlation with clinical outcomes.

RESULTS

Unsupervised clustering of mutations and data from RNA, DNA-copy-number, and DNA-methylation platforms uncovered concordant classification of three robust, nonoverlapping, prognostically significant subtypes of lower-grade glioma that were captured more accurately by IDH, 1p/19q, and TP53 status than by histologic class. Patients who had lower-grade gliomas with an IDH mutation and 1p/19q codeletion had the most favorable clinical outcomes. Their gliomas harbored mutations in CIC, FUBP1, NOTCH1, and the TERT promoter. Nearly all lower-grade gliomas with IDH mutations and no 1p/19q codeletion had mutations in TP53 (94%) and ATRX inactivation (86%). The large majority of lower-grade gliomas without an IDH mutation had genomic aberrations and clinical behavior strikingly similar to those found in primary glioblastoma.

CONCLUSIONS

The integration of genomewide data from multiple platforms delineated three molecular classes of lower-grade gliomas that were more concordant with IDH, 1p/19q, and TP53 status than with histologic class. Lower-grade gliomas with an IDH mutation either had 1p/19q codeletion or carried a TP53 mutation. Most lower-grade gliomas without an IDH mutation were molecularly and clinically similar to glioblastoma. (Funded by the National Institutes of Health.).

Abstract 3852: Obesity is associated with clinical characteristics in African American multiple myeloma patients

African-American ethnicity, male sex, older age and obesity are accepted risk factors for multiple myeloma (MM). Obesity early in life is a risk factor for many cancers, including MM; most studies have focused on populations of European origin. African-Americans have a higher prevalence of obesity than other populations, and may have a distinct genetic contribution to this condition. We established a multi-center study to investigate possible explanations for the excess risk of MM among African-Americans. The aim of the present case-case analysis was to determine whether body mass index (BMI) was associated with clinical characteristics at presentation in African-American MM patients.

1,065 patients diagnosed with active MM from January 1, 2009 through September 30, 2013 were recruited through nine outpatient centers at academic medical centers and 3 SEER population-based cancer registries. Information on weight and height at 20 years of age and at 5 years prior to diagnosis was obtained from questionnaires. Clinical information collected from the medical records and interviews included age at diagnosis, stage, percent plasmacytosis on bone marrow biopsy,

B2 microglobulin level, immunoglobulin class and presence of lytic bone lesions. Molecular profiling abstracted from pathology reports included chromosomal gain/deletions (1P-, 1Q+, +3, +5, +7, +9, +11, 12P-, Del(13), +15, Del(17P), Del(17P13)), translocations (4;14), (11;14), and (14;16 ). Analysis of Variance (ANOVA) was used to test the null hypothesis that age at diagnosis was similar across BMI categories. A T-test was used to test the null hypothesis that the mean age at diagnosis in underweight/normal patients was not different than that for overweight and obese patients. The Mantel-Haenszel Chi square test for trend was used to test the null hypothesis that the proportion of cytogenetic and molecular abnormalities was similar across BMI categories.

The study population consists of 306 African American male and 396 African American female MM patients. In both, obesity at 20 years of age was associated with younger age at diagnosis, compared to a normal/underweight (p=0.0015). Among males only, obesity at 5 years prior to diagnosis was also associated with younger age at diagnosis (p=0.03). A trend between increasing BMI and molecular characteristics associated with a poor prognosis was suggested.

We observed a strong association between obesity at age 20 and younger age at diagnosis among males and females, and in males only at 5 years prior to diagnosis. Obesity is one of the few known potentially modifiable risk factors for MM. Younger age at diagnosis reflects an earlier accumulation of either or both genetic and environmental risk factors. We also observed a link between obesity at each of the two time points and specific molecular markers of poor prognosis. Additional studies are needed to determine biological significance of the results.

Citation Format: Amie E. Hwang, Sikander Ailwadhi, Carol Ann Huff, Leon Bernal-Mizrachi, Christopher A. Haiman, Edward Peters, Seema Singhal, Karen Pawlish, Cathryn Bock, Todd Zimmerman, David J. Van Den Berg, David V. Conti, Brenda B. Birmann, Jayesh Mehta, John J. Graff, Daniel O. Stram, Niquelle Brown, Yang Yu, Moosa Azadian, Laurence Kolonel, Brian E. Henderson, Ann Mohrbacher, Graham Colditz, Brian-C Chiu, Michael Tomasson, Jeffrey Zonder, Robert Z. Orlowski, Sagar Lonial, Wendy Cozen. Obesity is associated with clinical characteristics in African American multiple myeloma patients. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 3852. doi:10.1158/1538-7445.AM2014-3852

Breast cancer susceptibility variants and mammographic density phenotypes in norwegian postmenopausal women

BACKGROUND

Mammographic density (MD) is one of the strongest known breast cancer risk factors. Twin studies have suggested that a large part of the variation in MD is genetically determined. We hypothesized that breast cancer susceptibility variants may affect MD, and that their effects may be modified by nongenetic factors.

METHODS

We assessed MD, using a computer-assisted method, on 2,348 postmenopausal Caucasian women (50-69 years) who participated in the Norwegian Breast Cancer Screening Program (NBCSP) in 2004 or 2006-07. We used linear regression (additive models) to determine the association between each SNP and MD, adjusting for age, body mass index (BMI), and study. We evaluated MD associations with 17 established breast cancer SNPs, overall, and by strata defined by non-genetic factors.

RESULTS

Two variants, 6q25.1-rs9383938 and TXNRD2-rs8141691, were statistically significantly associated with percent MD (P = 0.019 and 0.03, respectively), with the 6q25.1-rs9383938 association being consistent with the SNP effect on breast cancer risk. The effect of 6q25.1-rs3734805 on percent MD varied between parous and nulliparous women (Pinteraction = 0.02), whereas the effects of 9q31.2-rs865686 and MRPS30:FGF10-rs4415084 differed across strata of BMI (Pinteraction = 0.01 and 0.005, respectively). There was no evidence of effect modification by estrogen and progestin therapy use or alcohol consumption.

CONCLUSION

This study provides novel evidence of shared genetic risk factors between MD and breast cancer and of possible MD genetic-environmental interactions.

IMPACT

Although the results may be chance findings, they nevertheless highlight the need to investigate interactions with nongenetic factors in studies on the genetics of MD.

Reproducibility and reliability of SNP analysis using human cellular DNA at or near nanogram levels

BACKGROUND

Illumina SNP arrays have been routinely used for genome-wide association studies to identify potential biomarkers for various diseases. The recommended 200 ng of DNA for high-quality results is a roadblock to utilizing this assay when such quantities of DNA are not available. The goal of this study is to determine the reproducibility and reliability of the assay when reduced amounts of DNA are used for the SNP arrays.

FINDINGS

A serial 3-fold reduction of DNA from 200 ng to 0.8 ng was used for an Illumina SNP array in duplicates (200 ng, 66.7 ng, 22,2 ng, and 7.4 ng) or triplicates (2.47 ng and 0.8 ng). The reproducibility of the assay was determined by comparing allele calls (genotypes) at each locus within the duplicates or triplicates. The reliability of samples of reduced quantity was determined by comparing allele calls from samples of different quantities. As expected, the reproducibility and reliability both decrease with decreasing amounts of DNA used for the arrays. However, results of comparable quality to the 200 ng DNA recommended by Illumina can be obtained with much reduced amounts of DNA.

CONCLUSION

Reasonably reproducible and reliable results can be obtained with quantities of DNA, as low as 0.8 ng (equivalent to 133 human cells), well below the manufacturer's recommendation. Results of nearly equal quality to that of using 200 ng DNA can be obtained with 22.2 ng of DNA reliably, and clearly acceptable data can be obtained using 7.4 ng of DNA for Illumina SNP arrays.

Testicular germ cell tumor susceptibility associated with the UCK2 locus on chromosome 1q23

Genome-wide association studies (GWASs) have identified multiple common genetic variants associated with an increased risk of testicular germ cell tumors (TGCTs). A previous GWAS reported a possible TGCT susceptibility locus on chromosome 1q23 in the UCK2 gene, but failed to reach genome-wide significance following replication. We interrogated this region by conducting a meta-analysis of two independent GWASs including a total of 940 TGCT cases and 1559 controls for 122 single-nucleotide polymorphisms (SNPs) on chromosome 1q23 and followed up the most significant SNPs in an additional 2202 TGCT cases and 2386 controls from four case-control studies. We observed genome-wide significant associations for several UCK2 markers, the most significant of which was for rs3790665 (PCombined = 6.0 × 10(-9)). Additional support is provided from an independent familial study of TGCT where a significant over-transmission for rs3790665 with TGCT risk was observed (PFBAT = 2.3 × 10(-3)). Here, we provide substantial evidence for the association between UCK2 genetic variation and TGCT risk.

A genome-wide scan for breast cancer risk haplotypes among African American women

Genome-wide association studies (GWAS) simultaneously investigating hundreds of thousands of single nucleotide polymorphisms (SNP) have become a powerful tool in the investigation of new disease susceptibility loci. Haplotypes are sometimes thought to be superior to SNPs and are promising in genetic association analyses. The application of genome-wide haplotype analysis, however, is hindered by the complexity of haplotypes themselves and sophistication in computation. We systematically analyzed the haplotype effects for breast cancer risk among 5,761 African American women (3,016 cases and 2,745 controls) using a sliding window approach on the genome-wide scale. Three regions on chromosomes 1, 4 and 18 exhibited moderate haplotype effects. Furthermore, among 21 breast cancer susceptibility loci previously established in European populations, 10p15 and 14q24 are likely to harbor novel haplotype effects. We also proposed a heuristic of determining the significance level and the effective number of independent tests by the permutation analysis on chromosome 22 data. It suggests that the effective number was approximately half of the total (7,794 out of 15,645), thus the half number could serve as a quick reference to evaluating genome-wide significance if a similar sliding window approach of haplotype analysis is adopted in similar populations using similar genotype density.

Polymorphisms in hormone metabolism and growth factor genes and mammographic density in Norwegian postmenopausal hormone therapy users and non-users

INTRODUCTION

Mammographic density (MD) is one of the strongest known breast cancer risk factors. Estrogen and progestin therapy (EPT) has been associated with increases in MD. Dense breast tissue is characterized by increased stromal tissue and (to a lesser degree) increased numbers of breast epithelial cells. It is possible that genetic factors modify the association between EPT and MD, and that certain genetic variants are particularly important in determining MD in hormone users. We evaluated the association between MD and 340 tagging single nucleotide polymorphisms (SNPs) from about 30 candidate genes in hormone metabolism/growth factor pathways among women who participated in the Norwegian Breast Cancer Screening Program (NBCSP) in 2004.

METHODS

We assessed MD on 2,036 postmenopausal women aged 50 to 69 years using a computer-assisted method (Madena, University of Southern California) in a cross-sectional study. We used linear regression to determine the association between each SNP and MD, adjusting for potential confounders. The postmenopausal women were stratified into HT users (EPT and estrogen-only) and non-users (never HT).

RESULTS

For current EPT users, there was an association between a variant in the prolactin gene (PRL; rs10946545) and MD (dominant model, Bonferroni-adjusted P (Pb) = 0.0144). This association remained statistically significant among current users of norethisterone acetate (NETA)-based EPT, a regimen common in Nordic countries. Among current estrogen-only users (ET), there was an association between rs4670813 in the cytochrome P450 gene (CYP1B1) and MD (dominant model, Pb = 0.0396). In never HT users, rs769177 in the tumor necrosis factor (TNF) gene and rs1968752 in the region of the sulfotransferase gene (SULT1A1/SULT1A2), were significantly associated with MD (Pb = 0.0202; Pb = 0.0349).

CONCLUSIONS

We found some evidence that variants in the PRL gene were associated with MD in current EPT and NETA users. In never HT users, variants in the TNF and SULT1A1/SULT1A2 genes were significantly associated with MD. These findings may suggest that several genes in the hormone metabolism and growth factor pathways are implicated in determining MD.

Abstract 2634: Polymorphisms in DNA repair genes and risk of multiple myeloma

Background and Objective: DNA repair is an essential process for development and differentiation of all cells, including plasma cells. Associations between single nucleotide polymorphisms (SNPs) in DNA repair genes and multiple myeloma risk have been reported but generally not replicated. In a pilot study, we observed suggestive evidence of an association between multiple myeloma risk and five nonsynonymous SNPs in DNA repair genes (rs25489 in XRCC1, rs1801516 in ATM, rs2227999 in XPC, rs2228528 in CSB and rs2228615 in ICAM5). The objective of this study was to examine the effect of these five SNPs on multiple myeloma risk in a large multi-center analysis. Methods: We conducted a pooled analysis in a multi-ethnic sample of 489 multiple myeloma cases and 939 controls. Studies contributing samples included two population-based case-control studies conducted in Los Angeles and Seattle/Detroit that used the Surveillance, Epidemiology and End Results (SEER) cancer registries to identify cases, and three nested case-control studies using subjects from the Multiethnic Cohort (MEC), the Health Professionals Follow Up Study (HPFS), and the Nurses’ Health Study (NHS). The SNPs coded for an amino acid substitution, and were assumed to have the same effect across ethnic groups, and thus subjects from all racial/ethnic groups were included in the analysis. Assuming a log-additive inheritance model, odds ratios (OR) and 95% confidence intervals (95% CIs) were calculated for allele-specific risk comparing cases to controls using unconditional logistic regression, adjusting for age, race/ethnicity, sex, and study site. Results: We observed a 40% decreased risk per variant T allele in the XRCC1 SNP rs25489 (OR=0.60, 95% CI=0.40, 0.91, p-value=0.02). The minor allele frequency in cases and controls was 3.4% and 5.6%, respectively. Non-significant OR's &lt; 1.0 were observed when whites (OR=0.67, 95% CI=0.38, 1.17) and African Americans (OR=0.33, 95% CI= 0.09, 1.20) were considered separately. The SNPs in other DNA repair genes, ATM (rs1801516), XPC (rs2227999), CSB (rs2228528) and ICAM5 (rs2228615), showed no association with multiple myeloma risk. Conclusion: We observed an association between multiple myeloma and XRCC1, a gene that encodes a protein involved in the repair of single-stranded DNA breaks caused by ionizing radiation and alkylating agents. The variant allele of rs25489 results in a missense mutation substituting an amino acid (Arg &gt; His). Several studies have shown diminished DNA repair capacity associated with the variant allele of this SNP, an effect which may be relevant to multiple myeloma etiology. Replication in additional studies is needed to validate this finding.

Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 103rd Annual Meeting of the American Association for Cancer Research; 2012 Mar 31-Apr 4; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2012;72(8 Suppl):Abstract nr 2634. doi:1538-7445.AM2012-2634

Meta-analysis of genome-wide association studies of asthma in ethnically diverse North American populations

Asthma is a common disease with a complex risk architecture including both genetic and environmental factors. We performed a meta-analysis of North American genome-wide association studies of asthma in 5,416 individuals with asthma (cases) including individuals of European American, African American or African Caribbean, and Latino ancestry, with replication in an additional 12,649 individuals from the same ethnic groups. We identified five susceptibility loci. Four were at previously reported loci on 17q21, near IL1RL1, TSLP and IL33, but we report for the first time, to our knowledge, that these loci are associated with asthma risk in three ethnic groups. In addition, we identified a new asthma susceptibility locus at PYHIN1, with the association being specific to individuals of African descent (P = 3.9 × 10(-9)). These results suggest that some asthma susceptibility loci are robust to differences in ancestry when sufficiently large samples sizes are investigated, and that ancestry-specific associations also contribute to the complex genetic architecture of asthma.

Genetic variation in insulin-like growth factor 2 may play a role in ovarian cancer risk

The insulin-like growth factor (IGF) signaling axis plays an important role in cancer biology. We hypothesized that genetic variation in this pathway may influence risk of ovarian cancer. A three-center study of non-Hispanic whites including 1880 control women, 1135 women with invasive epithelial ovarian cancer and 321 women with borderline epithelial ovarian tumors was carried out to test the association between tag single-nucleotide polymorphisms (tSNPs) (n=58) in this pathway and risk of ovarian cancer. We found no association between variation in IGF1, IGFBP1 or IGFBP3 and risk of invasive disease, whereas five tSNPs in IGF2 were associated with risk of invasive epithelial ovarian cancer at P<0.05 and followed-up one of the associated SNPs. We conducted genotyping in 3216 additional non-Hispanic white cases and 5382 additional controls and were able to independently replicate our initial findings. In the combined set of studies, rs4320932 was associated with a 13% decreased risk of ovarian cancer per copy of the minor allele carried (95% confidence interval 0.81-0.93, P-trend=7.4 × 10(-5)). No heterogeneity of effect across study centers was observed (p(het)=0.25). IGF2 is emerging as an important gene for ovarian cancer; additional genotyping is warranted to further confirm these associations with IGF2 and to narrow down the region harboring the causal SNP.

Common alleles in candidate susceptibility genes associated with risk and development of epithelial ovarian cancer

Common germline genetic variation in the population is associated with susceptibility to epithelial ovarian cancer. Microcell-mediated chromosome transfer and expression microarray analysis identified nine genes associated with functional suppression of tumorogenicity in ovarian cancer cell lines; AIFM2, AKTIP, AXIN2, CASP5, FILIP1L, RBBP8, RGC32, RUVBL1 and STAG3. Sixty-three tagging single nucleotide polymorphisms (tSNPs) in these genes were genotyped in 1,799 invasive ovarian cancer cases and 3,045 controls to look for associations with disease risk. Two SNPs in RUVBL1, rs13063604 and rs7650365, were associated with increased risk of serous ovarian cancer [HetOR = 1.42 (1.15-1.74) and the HomOR = 1.63 (1.10-1.42), p-trend = 0.0002] and [HetOR = 0.97 (0.80-1.17), HomOR = 0.74 (0.58-0.93), p-trend = 0.009], respectively. We genotyped rs13063604 and rs7650365 in an additional 4,590 cases and 6,031 controls from ten sites from the United States, Europe and Australia; however, neither SNP was significant in Stage 2. We also evaluated the potential role of tSNPs in these nine genes in ovarian cancer development by testing for allele-specific loss of heterozygosity (LOH) in 286 primary ovarian tumours. We found frequent LOH for tSNPs in AXIN2, AKTIP and RGC32 (64, 46 and 34%, respectively) and one SNP, rs1637001, in STAG3 showed significant allele-specific LOH with loss of the common allele in 94% of informative tumours (p = 0.015). Array comparative genomic hybridisation indicated that this nonrandom allelic imbalance was due to amplification of the rare allele. In conclusion, we show evidence for the involvement of a common allele of STAG3 in the development of epithelial ovarian cancer.

Abstract LB-173: Genome-scale analysis of aberrant DNA methylation in colorectal cancer

Colorectal cancer (CRC) is a heterogeneous disease in which unique subtypes are characterized by distinct genetic and epigenetic alterations. Here we performed comprehensive genome-scale DNA methylation profiling of 125 colorectal tumors and 29 adjacent normal tissues. We identified four DNA methylation-based subgroups of CRC using model-based cluster analyses. Each subtype shows characteristic genetic and clinical features, indicating that they represent biologically distinct subgroups. A CIMP-high (CIMP-H) subgroup includes all CIMP-positive tumors characterized by the MethyLight five-marker panel, and is strongly associated with MLH1 DNA hypermethylation and the BRAFV600E mutation. A CIMP-low subgroup is enriched for KRAS mutations and characterized by DNA hypermethylation of a subset of CIMP-associated markers rather than a unique group of CpG islands. Non-CIMP tumors are separated into two distinct clusters. One non-CIMP subgroup is distinguished by a significantly higher frequency of TP53 mutations and frequent occurrence in the distal colon, while the tumors that belong to the fourth group exhibit a low frequency of both cancer-specific DNA hypermethylation and gene mutations and are significantly enriched for rectal tumors. Furthermore, we identified 114 genes that were downregulated more than 2-fold in CIMP-H tumors together with promoter DNA hypermethylation. These represent approximately 7% of genes that acquired promoter DNA methylation in CIMP-H tumors. Intriguingly, 48/114 genes were also transcriptionally silent in non-CIMP subgroups, but this was not attributable to promoter DNA hypermethylation. Together, we identified four distinct DNA methylation subgroups of CRC and provided novel insight regarding the role of CIMP-specific DNA hypermethylation in gene silencing.

Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 102nd Annual Meeting of the American Association for Cancer Research; 2011 Apr 2-6; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2011;71(8 Suppl):Abstract nr LB-173. doi:10.1158/1538-7445.AM2011-LB-173

Genetic variations on chromosomes 5p15 and 15q25 and bladder cancer risk: findings from the Los Angeles-Shanghai bladder case-control study

Genome-wide association studies have associated common variations at chromosomes 5p15 and 15q25 with lung cancer risk. The 5p15 locus has also been associated with increased bladder cancer risk in a recent report. The 15q25 locus has been associated with nicotine dependence and self-reported number of cigarettes smoked per day in some studies and it was proposed that its association with lung cancer may be mediated through differences in smoking behavior. Here, we investigated the roles of variations at 5p15 (rs401681, rs402710, rs2736098 and rs2736100) and 15q25 (rs1051730 and rs8034191) in bladder cancer etiology in two case-control studies conducted separately in Los Angeles County, CA, USA (498 cases and 588 controls) and in Shanghai, China (506 cases and 530 controls). We replicated the association between the 5p15 locus and bladder cancer among non-Hispanic whites (NHW) in Los Angeles [for rs2736100, per C allele odds ratio (OR) = 1.23; 95% confidence interval (CI), 1.02-1.48; P = 0.029] and among Chinese in Shanghai (OR = 1.22; 95% CI, 1.02-1.47; P = 0.033). Both rs1051730 and rs8034191 at 15q25 were rare among Chinese. Among NHW, a significant association was found between rs8034191 and bladder cancer which persisted after adjustment for cigarette smoking status, number of cigarettes smoked per day and number of years of smoking (per C allele OR = 1.26; 95% CI, 1.04-1.54; P = 0.017). Our results support 5p15 and 15q25 as susceptibility regions for bladder cancer risk.

Common variants at 19p13 are associated with susceptibility to ovarian cancer

Epithelial ovarian cancer (EOC) is the leading cause of death from gynecological malignancy in the developed world, accounting for 4% of the deaths from cancer in women. We performed a three-phase genome-wide association study of EOC survival in 8,951 individuals with EOC (cases) with available survival time data and a parallel association analysis of EOC susceptibility. Two SNPs at 19p13.11, rs8170 and rs2363956, showed evidence of association with survival (overall P = 5 × 10⁻⁴ and P = 6 × 10⁻⁴, respectively), but they did not replicate in phase 3. However, the same two SNPs demonstrated genome-wide significance for risk of serous EOC (P = 3 × 10⁻⁹ and P = 4 × 10⁻¹¹, respectively). Expression analysis of candidate genes at this locus in ovarian tumors supported a role for the BRCA1-interacting gene C19orf62, also known as MERIT40, which contains rs8170, in EOC development.

Progesterone receptor gene polymorphisms and risk of endometriosis: results from an international collaborative effort

OBJECTIVE

To investigate the association between self-reported endometriosis and the putative functional promoter +331C/T single nucleotide polymorphism and the PROGINS allele.

DESIGN

Control subjects from ovarian cancer case-control studies participating in the international Ovarian Cancer Association Consortium. The majority of controls are drawn from population-based studies.

SETTING

An international ovarian cancer consortium including studies from Australia, Europe, and the United States.

PATIENT(S)

Five thousand eight hundred twelve white female controls, of whom 348 had endometriosis, from eight ovarian cancer case-control studies.

INTERVENTION(S)

None.

MAIN OUTCOME MEASURE(S)

Genotypes for the +331C/T single nucleotide polymorphism and PROGINS allele and a history of endometriosis.

RESULT(S)

The occurrence of endometriosis was reduced in women carrying one or more copies of the +331 T allele (odds ratio=0.65; 95% confidence interval: 0.43-0.98), whereas there was no association between the PROGINS allele and endometriosis (odds ratio=0.94, 95% confidence interval 0.76-1.16).

CONCLUSION(S)

Additional studies of the +331C/T variant are warranted given the current finding and the equivocal results of previous studies. The +331 T allele has been shown to result in a reduced progesterone (P) receptor A to P receptor B ratio, and if the observed association with endometriosis is confirmed it would suggest that this ratio is important for this disease.

ESR1/SYNE1 polymorphism and invasive epithelial ovarian cancer risk: an Ovarian Cancer Association Consortium study

We genotyped 13 single nucleotide polymorphisms (SNPs) in the estrogen receptor alpha gene (ESR1) region in three population-based case-control studies of epithelial ovarian cancer conducted in the United States, comprising a total of 1,128 and 1,866 non-Hispanic white invasive cases and controls, respectively. A SNP 19 kb downstream of ESR1 (rs2295190, G-to-T change) was associated with invasive ovarian cancer risk, with a per-T-allele odds ratio (OR) of 1.24 [95% confidence interval (CI), 1.06-1.44, P = 0.006]. rs2295190 is a nonsynonymous coding SNP in a neighboring gene called spectrin repeat containing, nuclear envelope 1 (SYNE1), which is involved in nuclear organization and structural integrity, function of the Golgi apparatus, and cytokinesis. An isoform encoded by SYNE1 has been reported to be downregulated in ovarian and other cancers. rs2295190 was genotyped in an additional 12 studies through the Ovarian Cancer Association Consortium, with 5,279 invasive epithelial cases and 7,450 controls. The per-T-allele OR for this 12-study set was 1.09 (95% CI, 1.02-1.17; P = 0.017). Results for the serous subtype in the 15 combined studies were similar to those overall (n = 3,545; OR, 1.09; 95% CI, 1.01-1.18; P = 0.025), and our findings were strongest for the mucinous subtype (n = 447; OR, 1.32; 95% CI, 1.11-1.58; P = 0.002). No association was observed for the endometrioid subtype. In an additional analysis of 1,459 borderline ovarian cancer cases and 7,370 controls, rs2295190 was not associated with risk. These data provide suggestive evidence that the rs2295190 T allele, or another allele in linkage disequilibrium with it, may be associated with increased risk of invasive ovarian cancer.

Polymorphism in the GALNT1 gene and epithelial ovarian cancer in non-Hispanic white women: the Ovarian Cancer Association Consortium

Aberrant glycosylation is a well-described hallmark of cancer. In a previous ovarian cancer case control study that examined polymorphisms in 26 glycosylation-associated genes, we found strong statistical evidence (P = 0.00017) that women who inherited two copies of a single-nucleotide polymorphism in the UDP-N-acetylgalactosamine:polypeptide N-acetylgalactosaminyltransferase, GALNT1, had decreased ovarian cancer risk. The current study attempted to replicate this observation. The GALNT1 single-nucleotide polymorphism rs17647532 was genotyped in 6,965 cases and 8,377 controls from 14 studies forming the Ovarian Cancer Association Consortium. The fixed effects estimate per rs17647532 allele was null (odds ratio, 0.99; 95% confidence interval, 0.92-1.07). When a recessive model was fit, the results were unchanged. Test for heterogeneity of the odds ratios revealed consistency across the 14 replication sites but significant differences compared with the original study population (P = 0.03). This study underscores the need for replication of putative findings in genetic association studies.

Genetic Association Between the COMT Genotype and Urinary Levels of Tea Polyphenols and Their Metabolites among Daily Green Tea Drinkers

Available in vitro and animal studies have shown cancer protective effects of tea polyphenols. Recent study suggests a greater protective effect of green tea intake on breast cancer risk among women possessing the low-activity associated genotype of the catechol-O-methyltransferase (COMT), which may modulate the metabolism and excretion of tea polyphenols through urine. To determine the effect of COMT genotype on urinary excretion of tea polyphenol metabolites of daily green tea drinkers, a cross-sectional analysis was performed within the Shanghai Cohort Study, a population-based, prospective investigation of diet and cancer in 18,244men. In addition to an in-person interview, each participant provided a blood and urine sample at baseline. In the present study, COMT genotype (rs4680) and five urinary metabolites of tea polyphenols were determined in 660 cohort subjects who self-identified as daily drinkers of green tea. All urinary tea polyphenol measurements were expressed in units of urinary creatinine. Men possessing the homozygous low-activity associated COMT genotype (LL) exhibited statistically significantly lower urinary levels of individual as well as all of the five tea polyphenol metabolites under study relative to individuals possessing the wild type high-activity associated COMT genotype (HH)or the heterozygous variant genotype (HL). Levels of urinary tea polyphenol metabolites were comparable between men possessing the HH and HL genotypes. The present study demonstrated that men carrying low-activity associated COMT genotype excreted less tea polyphenols from urine, which suggests that they may retain more tea polyphenols in their bodies and derive greater health benefits from green tea intake.

A genome-wide association study identifies a new ovarian cancer susceptibility locus on 9p22.2

Epithelial ovarian cancer has a major heritable component, but the known susceptibility genes explain less than half the excess familial risk. We performed a genome-wide association study (GWAS) to identify common ovarian cancer susceptibility alleles. We evaluated 507,094 SNPs genotyped in 1,817 cases and 2,353 controls from the UK and approximately 2 million imputed SNPs. We genotyped the 22,790 top ranked SNPs in 4,274 cases and 4,809 controls of European ancestry from Europe, USA and Australia. We identified 12 SNPs at 9p22 associated with disease risk (P < 10(-8)). The most significant SNP (rs3814113; P = 2.5 x 10(-17)) was genotyped in a further 2,670 ovarian cancer cases and 4,668 controls, confirming its association (combined data odds ratio (OR) = 0.82, 95% confidence interval (CI) 0.79-0.86, P(trend) = 5.1 x 10(-19)). The association differs by histological subtype, being strongest for serous ovarian cancers (OR 0.77, 95% CI 0.73-0.81, P(trend) = 4.1 x 10(-21)).

Association between invasive ovarian cancer susceptibility and 11 best candidate SNPs from breast cancer genome-wide association study

Because both ovarian and breast cancer are hormone-related and are known to have some predisposition genes in common, we evaluated 11 of the most significant hits (six with confirmed associations with breast cancer) from the breast cancer genome-wide association study for association with invasive ovarian cancer. Eleven SNPs were initially genotyped in 2927 invasive ovarian cancer cases and 4143 controls from six ovarian cancer case-control studies. Genotype frequencies in cases and controls were compared using a likelihood ratio test in a logistic regression model stratified by study. Initially, three SNPs (rs2107425 in MRPL23, rs7313833 in PTHLH, rs3803662 in TNRC9) were weakly associated with ovarian cancer risk and one SNP (rs4954956 in NXPH2) was associated with serous ovarian cancer in non-Hispanic white subjects (P-trend < 0.1). These four SNPs were then genotyped in an additional 4060 cases and 6308 controls from eight independent studies. Only rs4954956 was significantly associated with ovarian cancer risk both in the replication study and in combined analyses. This association was stronger for the serous histological subtype [per minor allele odds ratio (OR) 1.07 95% CI 1.01-1.13, P-trend = 0.02 for all types of ovarian cancer and OR 1.14 95% CI 1.07-1.22, P-trend = 0.00017 for serous ovarian cancer]. In conclusion, we found that rs4954956 was associated with increased ovarian cancer risk, particularly for serous ovarian cancer. However, none of the six confirmed breast cancer susceptibility variants we tested was associated with ovarian cancer risk. Further work will be needed to identify the causal variant associated with rs4954956 or elucidate its function.

Green tea intake, MTHFR/TYMS genotype and breast cancer risk: the Singapore Chinese Health Study

The tea polyphenol (-)-epigallocatechin-3-gallate (EGCG) has been reported to act as a cancer preventive agent through folate pathway inhibition in experimental studies. We hypothesized that if folate pathway inhibition is the mechanism of cancer preventive activities of EGCG, then the protective effect against breast cancer would be stronger among women with low dietary folate intake and the high-activity methylenetetrahydrofolate reductase (MTHFR) and thymidylate synthase (TYMS) genotypes. In a nested case-control study of 380 women with incident breast cancer and 662 controls within the Singapore Chinese Health Study, we found no association between either green tea intake or gene polymorphisms of MTHFR (C677T and A1298C) and TYMS (1494 ins/del) and breast cancer risk. However, among women with low folate intake (<133.4 microg/day), weekly/daily green tea intake was inversely associated with breast cancer risk compared with less green tea intake [odds ratio (OR) = 0.45, 95% confidence interval (CI) = 0.26-0.79, P for interaction = 0.02]. Among women with high folate intake (>or=133.4 microg/day), green tea intake was not associated with breast cancer. Similarly, among women possessing the high-activity MTHFR/TYMS genotypes (0-1 variant allele), weekly/daily versus less frequent green tea intake was associated with lower breast cancer risk (OR = 0.66, 95% CI = 0.45-0.98), which was observed even more strongly among those who also had low folate intake (OR = 0.44, 95% CI = 0.22-0.89) than high folate intake (OR = 0.92, 95% CI = 0.55-1.54). This association was not observed among women possessing the low-activity genotypes (2-4 variant alleles). Our findings suggest that folate pathway inhibition may be one mechanism through which green tea protects against breast cancer in humans.

Consortium analysis of 7 candidate SNPs for ovarian cancer

The Ovarian Cancer Association Consortium selected 7 candidate single nucleotide polymorphisms (SNPs), for which there is evidence from previous studies of an association with variation in ovarian cancer or breast cancer risks. The SNPs selected for analysis were F31I (rs2273535) in AURKA, N372H (rs144848) in BRCA2, rs2854344 in intron 17 of RB1, rs2811712 5' flanking CDKN2A, rs523349 in the 3' UTR of SRD5A2, D302H (rs1045485) in CASP8 and L10P (rs1982073) in TGFB1. Fourteen studies genotyped 4,624 invasive epithelial ovarian cancer cases and 8,113 controls of white non-Hispanic origin. A marginally significant association was found for RB1 when all studies were included [ordinal odds ratio (OR) 0.88 (95% confidence interval (CI) 0.79-1.00) p = 0.041 and dominant OR 0.87 (95% CI 0.76-0.98) p = 0.025]; when the studies that originally suggested an association were excluded, the result was suggestive although no longer statistically significant (ordinal OR 0.92, 95% CI 0.79-1.06). This SNP has also been shown to have an association with decreased risk in breast cancer. There was a suggestion of an association for AURKA, when one study that caused significant study heterogeneity was excluded [ordinal OR 1.10 (95% CI 1.01-1.20) p = 0.027; dominant OR 1.12 (95% CI 1.01-1.24) p = 0.03]. The other 5 SNPs in BRCA2, CDKN2A, SRD5A2, CASP8 and TGFB1 showed no association with ovarian cancer risk; given the large sample size, these results can also be considered to be informative. These null results for SNPs identified from relatively large initial studies shows the importance of replicating associations by a consortium approach.

No association between the SRD5A2 gene A49T missense variant and prostate cancer risk: lessons learned

The steroid 5-alpha reductase type II gene (SRD5A2) encodes the enzyme which converts testosterone (T) to the more active androgen dihydrotestosterone. A non-synonymous single-nucleotide polymorphism, A49T (rs9282858), in SRD5A2 has been implicated in prostate cancer risk; however, results have been inconsistent. In 1999, we reported a strong association between the A49T variant and prostate cancer risk among African-Americans and Latinos in the Hawaii-Los Angeles Multiethnic Cohort (MEC). We report here an updated analysis of MEC data including the five major ethnic groups of the MEC, an increased sample size, improved genotyping technology and a comprehensive meta-analysis of the published literature. We found a non-statistically significant positive association between prostate cancer risk and carrying either the AT or TT genotype [odds ratio (OR) = 1.16, 95% confidence interval (CI) 0.79-1.69] in the MEC. This finding is in contrast to our previous results of ORs of 3.28 and 2.50 for the association between prostate cancer risk and the variant in African-American and Latino men, respectively; this can be accounted for by genotyping error in our earlier study. Meta-analysis of the published literature, including the current MEC data, shows a summary OR of 1.13 (95% CI 0.95-1.34) for the A49T variant with prostate cancer risk among sporadic, unselected cases. After evaluating more than 6000 cases and 6000 controls, there is little evidence of a role for the SRD5A2 A49T variant in prostate cancer risk. Overall, this report highlights the importance of rigorous genotyping quality control measures and replication efforts in genetic association studies.

Interleukin-6-Related Genotypes, Body Mass Index, and Risk of Multiple Myeloma and Plasmacytoma

Interleukin-6 (IL-6) promotes normal plasma cell development and proliferation of myeloma cells in culture. We evaluated IL-6 genotypes and body mass index (BMI) in a case-control study of multiple myeloma and plasmacytoma. DNA samples and questionnaires were obtained from incident cases of multiple myeloma (n = 134) and plasmacytoma (n = 16; plasma cell neoplasms) ascertained from the Los Angeles County population-based cancer registry and from siblings or cousins of cases (family controls, n = 112) and population controls (n = 126). Genotypes evaluated included IL-6 promoter gene single nucleotide polymorphisms (SNP) at positions -174, -572, and -597; one variable number of tandem repeats (-373 A(n)T(n)); and one SNP in the IL-6 receptor (IL-6ralpha) gene at position -358. The variant allele of the IL-6 promoter SNP -572 was associated with a roughly 2-fold increased risk of plasma cell neoplasms when cases were compared with family [odds ratio (OR), 1.8; 95% confidence interval (95% CI), 0.7-4.7] or population controls (OR, 2.4; 95% CI, 1.2-4.7). The -373 9A/9A genotype was associated with a decreased risk compared with the most common genotype (OR for cases versus family controls, 0.4; 95% CI, 0.1-1.7; OR for cases versus population controls, 0.3; 95% CI, 0.1-0.9). No other SNPs were associated with risk. Obesity (BMI >or= 30 kg/m(2)) increased risk nonsignificantly by 40% and 80% when cases were compared with family controls or population controls, respectively, relative to persons with a BMI of <25 kg/m(2). These results suggest that IL-6 promoter genotypes may be associated with increased risk of plasma cell neoplasms.