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.

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 FEV1, FVC, and FEV1/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 FEV1, FVC, or FEV1/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.

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) .

Genome-wide interaction studies reveal sex-specific asthma risk alleles.

Asthma is a complex disease with sex-specific differences in prevalence. Candidate gene studies have suggested that genotype-by-sex interaction effects on asthma risk exist, but this has not yet been explored at a genome-wide level. We aimed to identify sex-specific asthma risk alleles by performing a genome-wide scan for genotype-by-sex interactions in the ethnically diverse participants in the EVE Asthma Genetics Consortium. We performed male- and female-specific genome-wide association studies in 2653 male asthma cases, 2566 female asthma cases and 3830 non-asthma controls from European American, African American, African Caribbean and Latino populations. Association tests were conducted in each study sample, and the results were combined in ancestry-specific and cross-ancestry meta-analyses. Six sex-specific asthma risk loci had P-values < 1 × 10(-6), of which two were male specific and four were female specific; all were ancestry specific. The most significant sex-specific association in European Americans was at the interferon regulatory factor 1 (IRF1) locus on 5q31.1. We also identify a Latino female-specific association in RAP1GAP2. Both of these loci included single-nucleotide polymorphisms that are known expression quantitative trait loci and have been associated with asthma in independent studies. The IRF1 locus is a strong candidate region for male-specific asthma susceptibility due to the association and validation we demonstrate here, the known role of IRF1 in asthma-relevant immune pathways and prior reports of sex-specific differences in interferon responses.

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 meta-analysis of nodular sclerosing Hodgkin lymphoma identifies risk loci at 6p21.32.

Nodular sclerosing Hodgkin lymphoma (NSHL) is a distinct, highly heritable Hodgkin lymphoma subtype. We undertook a genome-wide meta-analysis of 393 European-origin adolescent/young adult NSHL patients and 3315 controls using the Illumina Human610-Quad Beadchip and Affymetrix Genome-Wide Human SNP Array 6.0. We identified 3 single nucleotide polymorphisms (SNPs) on chromosome 6p21.32 that were significantly associated with NSHL risk: rs9268542 (P = 5.35 × 10(-10)), rs204999 (P = 1.44 × 10(-9)), and rs2858870 (P = 1.69 × 10(-8)). We also confirmed a previously reported association in the same region, rs6903608 (P = 3.52 × 10(-10)). rs204999 and rs2858870 were weakly correlated (r(2) = 0.257), and the remaining pairs of SNPs were not correlated (r(2) < 0.1). In an independent set of 113 NSHL cases and 214 controls, 2 SNPs were significantly associated with NSHL and a third showed a comparable odds ratio (OR). These SNPs are found on 2 haplotypes associated with NSHL risk (rs204999-rs9268528-rs9268542-rs6903608-rs2858870; AGGCT, OR = 1.7, P = 1.71 × 10(-6); GAATC, OR = 0.4, P = 1.16 × 10(-4)). All individuals with the GAATC haplotype also carried the HLA class II DRB1*0701 allele. In a separate analysis, the DRB1*0701 allele was associated with a decreased risk of NSHL (OR = 0.5, 95% confidence interval = 0.4, 0.7). These data support the importance of the HLA class II region in NSHL etiology.

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.

A genome-wide association study of breast cancer in women of African ancestry.

Genome-wide association studies (GWAS) in diverse populations are needed to reveal variants that are more common and/or limited to defined populations. We conducted a GWAS of breast cancer in women of African ancestry, with genotyping of >1,000,000 SNPs in 3,153 African American cases and 2,831 controls, and replication testing of the top 66 associations in an additional 3,607 breast cancer cases and 11,330 controls of African ancestry. Two of the 66 SNPs replicated (p < 0.05) in stage 2, which reached statistical significance levels of 10(-6) and 10(-5) in the stage 1 and 2 combined analysis (rs4322600 at chromosome 14q31: OR = 1.18, p = 4.3 × 10(-6); rs10510333 at chromosome 3p26: OR = 1.15, p = 1.5 × 10(-5)). These suggestive risk loci have not been identified in previous GWAS in other populations and will need to be examined in additional samples. Identification of novel risk variants for breast cancer in women of African ancestry will demand testing of a substantially larger set of markers from stage 1 in a larger replication sample.

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.

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.

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 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.

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.

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.

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.

A meta-analysis of genome-wide association studies of multiple myeloma among men and women of African ancestry.

Persons of African ancestry (AA) have a twofold higher risk for multiple myeloma (MM) compared with persons of European ancestry (EA). Genome-wide association studies (GWASs) support a genetic contribution to MM etiology in individuals of EA. Little is known about genetic risk factors for MM in individuals of AA. We performed a meta-analysis of 2 GWASs of MM in 1813 cases and 8871 controls and conducted an admixture mapping scan to identify risk alleles. We 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. In GWAS meta-analysis, we identified 2 suggestive novel loci located at 9p24.3 and 9p13.1 at P < 1 × 10-6; however, 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. Of the 23 known EA risk variants, 20 showed directional consistency, and 9 replicated at P < .05 in AA individuals. In 8 regions, we identified markers that better capture MM risk in persons with AA. AA individuals with a PRS in the top 10% had a 1.82-fold (95% confidence interval, 1.56-2.11) increased MM risk compared with 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). Our study shows that common genetic variation contributes to MM risk in individuals with AA.

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.

Snagger: a user-friendly program for incorporating additional information for tagSNP selection.

BACKGROUND: There has been considerable effort focused on developing efficient programs for tagging single-nucleotide polymorphisms (SNPs). Many of these programs do not account for potential reduced genomic coverage resulting from genotyping failures nor do they preferentially select SNPs based on functionality, which may be more likely to be biologically important. RESULTS: We have developed a user-friendly and efficient software program, Snagger, as an extension to the existing open-source software, Haploview, which uses pairwise r2 linkage disequilibrium between single nucleotide polymorphisms (SNPs) to select tagSNPs. Snagger distinguishes itself from existing SNP selection algorithms, including Tagger, by providing user options that allow for: (1) prioritization of tagSNPs based on certain characteristics, including platform-specific design scores, functionality (i.e., coding status), and chromosomal position, (2) efficient selection of SNPs across multiple populations, (3) selection of tagSNPs outside defined genomic regions to improve coverage and genotyping success, and (4) picking of surrogate tagSNPs that serve as backups for tagSNPs whose failure would result in a significant loss of data. Using HapMap genotype data from ten ENCODE regions and design scores for the Illumina platform, we show similar coverage and design score distribution and fewer total tagSNPs selected by Snagger compared to the web server Tagger. CONCLUSION: Snagger improves upon current available tagSNP software packages by providing a means for researchers to select tagSNPs that reliably capture genetic variation across multiple populations while accounting for significant genotyping failure risk and prioritizing on SNP-specific characteristics.

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.