Fine-mapping analysis including over 254,000 East Asian and European descendants identifies 136 putative colorectal cancer susceptibility genes

Genome-wide association studies (GWAS) have identified more than 200 common genetic variants independently associated with colorectal cancer (CRC) risk, but the causal variants and target genes are mostly unknown. We sought to fine-map all known CRC risk loci using GWAS data from 100,204 cases and 154,587 controls of East Asian and European ancestry. Our stepwise conditional analyses revealed 238 independent association signals of CRC risk, each with a set of credible causal variants (CCVs), of which 28 signals had a single CCV. Our cis-eQTL/mQTL and colocalization analyses using colorectal tissue-specific transcriptome and methylome data separately from 1299 and 321 individuals, along with functional genomic investigation, uncovered 136 putative CRC susceptibility genes, including 56 genes not previously reported. Analyses of single-cell RNA-seq data from colorectal tissues revealed 17 putative CRC susceptibility genes with distinct expression patterns in specific cell types. Analyses of whole exome sequencing data provided additional support for several target genes identified in this study as CRC susceptibility genes. Enrichment analyses of the 136 genes uncover pathways not previously linked to CRC risk. Our study substantially expanded association signals for CRC and provided additional insight into the biological mechanisms underlying CRC development.

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.

Rare Variants in the DNA Repair Pathway and the Risk of Colorectal Cancer

BACKGROUND

Inherited susceptibility is an important contributor to colorectal cancer risk, and rare variants in key genes or pathways could account in part for the missing proportion of colorectal cancer heritability.

METHODS

We conducted an exome-wide association study including 2,327 cases and 2,966 controls of European ancestry from three large epidemiologic studies. Single variant associations were tested using logistic regression models, adjusting for appropriate study-specific covariates. In addition, we examined the aggregate effects of rare coding variation at the gene and pathway levels using Bayesian model uncertainty techniques.

RESULTS

In an exome-wide gene-level analysis, we identified ST6GALNAC2 as the top associated gene based on the Bayesian risk index (BRI) method [summary Bayes factor (BF)_BRI = 2604.23]. A rare coding variant in this gene, rs139401613, was the top associated variant (P = 1.01 × 10^-6) in an exome-wide single variant analysis. Pathway-level association analyses based on the integrative BRI (iBRI) method found extreme evidence of association with the DNA repair pathway (BF_iBRI = 17852.4), specifically with the nonhomologous end joining (BF_iBRI = 437.95) and nucleotide excision repair (BF_iBRI = 36.96) subpathways. The iBRI method also identified RPA2, PRKDC, ERCC5, and ERCC8 as the top associated DNA repair genes (summary BF_iBRI ≥ 10), with rs28988897, rs8178232, rs141369732, and rs201642761 being the most likely associated variants in these genes, respectively.

CONCLUSIONS

We identified novel variants and genes associated with colorectal cancer risk and provided additional evidence for a role of DNA repair in colorectal cancer tumorigenesis.

IMPACT

This study provides new insights into the genetic predisposition to colorectal cancer, which has potential for translation into improved risk prediction.

Cumulative Burden of Colorectal Cancer-Associated Genetic Variants Is More Strongly Associated With Early-Onset vs Late-Onset Cancer

BACKGROUND & AIMS

Early-onset colorectal cancer (CRC, in persons younger than 50 years old) is increasing in incidence; yet, in the absence of a family history of CRC, this population lacks harmonized recommendations for prevention. We aimed to determine whether a polygenic risk score (PRS) developed from 95 CRC-associated common genetic risk variants was associated with risk for early-onset CRC.

METHODS

We studied risk for CRC associated with a weighted PRS in 12,197 participants younger than 50 years old vs 95,865 participants 50 years or older. PRS was calculated based on single nucleotide polymorphisms associated with CRC in a large-scale genome-wide association study as of January 2019. Participants were pooled from 3 large consortia that provided clinical and genotyping data: the Colon Cancer Family Registry, the Colorectal Transdisciplinary Study, and the Genetics and Epidemiology of Colorectal Cancer Consortium and were all of genetically defined European descent. Findings were replicated in an independent cohort of 72,573 participants.

RESULTS

Overall associations with CRC per standard deviation of PRS were significant for early-onset cancer, and were stronger compared with late-onset cancer (P for interaction = .01); when we compared the highest PRS quartile with the lowest, risk increased 3.7-fold for early-onset CRC (95% CI 3.28-4.24) vs 2.9-fold for late-onset CRC (95% CI 2.80-3.04). This association was strongest for participants without a first-degree family history of CRC (P for interaction = 5.61 × 10-5). When we compared the highest with the lowest quartiles in this group, risk increased 4.3-fold for early-onset CRC (95% CI 3.61-5.01) vs 2.9-fold for late-onset CRC (95% CI 2.70-3.00). Sensitivity analyses were consistent with these findings.

CONCLUSIONS

In an analysis of associations with CRC per standard deviation of PRS, we found the cumulative burden of CRC-associated common genetic variants to associate with early-onset cancer, and to be more strongly associated with early-onset than late-onset cancer, particularly in the absence of CRC family history. Analyses of PRS, along with environmental and lifestyle risk factors, might identify younger individuals who would benefit from preventive measures.

Publisher Correction: Shared heritability and functional enrichment across six solid cancers

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

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.

Shared heritability and functional enrichment across six solid cancers

Quantifying the genetic correlation between cancers can provide important insights into the mechanisms driving cancer etiology. Using genome-wide association study summary statistics across six cancer types based on a total of 296,215 cases and 301,319 controls of European ancestry, here we estimate the pair-wise genetic correlations between breast, colorectal, head/neck, lung, ovary and prostate cancer, and between cancers and 38 other diseases. We observed statistically significant genetic correlations between lung and head/neck cancer (rg = 0.57, p = 4.6 × 10-8), breast and ovarian cancer (rg = 0.24, p = 7 × 10-5), breast and lung cancer (rg = 0.18, p =1.5 × 10-6) and breast and colorectal cancer (rg = 0.15, p = 1.1 × 10-4). We also found that multiple cancers are genetically correlated with non-cancer traits including smoking, psychiatric diseases and metabolic characteristics. Functional enrichment analysis revealed a significant excess contribution of conserved and regulatory regions to cancer heritability. Our comprehensive analysis of cross-cancer heritability suggests that solid tumors arising across tissues share in part a common germline genetic basis.

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.

Abstract 491: Elucidating de novo PATRR-mediated t(3;8) balanced translocation and clear cell renal cell carcinoma

Palindromic AT-Rich Repeat (PATRR) - mediated translocations entail exchange across chromosomes at sites enriched in palindromic repeats of the nucleotides adenine (A) and thymine (T). Their precise embryologic origin and associated pathobiology with clear cell renal cell carcinoma (ccRCC) remain incompletely described. In the present study we document an individual with familial non-VHL ccRCC (7 primary renal tumors) who harbors a germline de novo PATRR mediated balanced translocation involving chromosomes 3 and 8 [t(3;8)] validated by spectral karyotyping (SKY).

Using translocation specific PCR and DNA sequencing we determined the chromosome 3 breakpoint to be located in an AT-rich palindromic sequence in the third intron of the FHIT gene (chr3p14.2) and the chromosome 8 breakpoint in the first intron of the TRC8 gene (chr8q24.1). Genotyping analysis, using a high density custom exomechip array by Illumina, revealed a loss of the entire aberrant chromosome 8 carrying the 3p segment [der(8)] in all renal tumors tested. We also determined that the [46, XY, t(3;8), (p14.2, q24.1)] translocation was paternally derived by performing a genotypic assessment of the regions that differ between the parental alleles and then establishing which haplotypes are associated with the translocation. The somatic mutational landscape was assessed by Whole Exome sequencing of the renal tumors and the proband’s germline DNA. No germline or somatic deleterious mutations were detected in VHL gene suggesting that the ccRCC phenotype is not associated with Von Hippel-Lindau disease.

Furthermore, we measured the transcriptomic profiles of the renal tumors and matched normal tissues employing barcoded probe hybridization technology (NanoString) and RNASeq to fully characterize differential gene expression and define gene ontology networks that are dysregulated. Most significantly, no difference in expression of the VHL gene was detected between tumors and normal counterparts. Aberrant regulation was detected in members of WNT signaling pathway, TGF-beta and MAPK pathways, in addition to NOTCH signaling pathway genes and members of the VEGF family. In summary, in this study we employed advanced genomic and transcriptomic techniques, to delineate the embryologic origin and deleterious biological consequences of the PATRR-mediated t(3;8) balanced translocation associated with clear cell renal cell carcinoma.

Citation Format: Marilena Melas, Kevin J. McDonnell, Christopher K. Edlund, Sarah J. Tash, Duveen Y. Sturgeon, Chenxu Qu, Charalampos Lazaris, Peter J. Gruber, Thomas W. Glover, Beverly S. Emanuel, Stephen B. Gruber. Elucidating de novo PATRR-mediated t(3;8) balanced translocation and clear cell renal cell carcinoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 491. doi:10.1158/1538-7445.AM2017-491

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.

Cross-Cancer Genome-Wide Analysis of Lung, Ovary, Breast, Prostate, and Colorectal Cancer Reveals Novel Pleiotropic Associations

Identifying genetic variants with pleiotropic associations can uncover common pathways influencing multiple cancers. We took a two-stage approach to conduct genome-wide association studies for lung, ovary, breast, prostate, and colorectal cancer from the GAME-ON/GECCO Network (61,851 cases, 61,820 controls) to identify pleiotropic loci. Findings were replicated in independent association studies (55,789 cases, 330,490 controls). We identified a novel pleiotropic association at 1q22 involving breast and lung squamous cell carcinoma, with eQTL analysis showing an association with ADAM15/THBS3 gene expression in lung. We also identified a known breast cancer locus CASP8/ALS2CR12 associated with prostate cancer, a known cancer locus at CDKN2B-AS1 with different variants associated with lung adenocarcinoma and prostate cancer, and confirmed the associations of a breast BRCA2 locus with lung and serous ovarian cancer. This is the largest study to date examining pleiotropy across multiple cancer-associated loci, identifying common mechanisms of cancer development and progression. Cancer Res; 76(17); 5103-14. ©2016 AACR.

Abstract LB-365: Novel susceptibility loci for colorectal cancer: Findings from the colorectal transdisciplinary (CORECT) study OncoArray analysis

Genome-wide association studies (GWAS) have identified 58 susceptibility alleles across 37 regions associated with the incidence of colorectal cancer (CRC) with P&lt;5E-08, yet much of the disease's familial risk remains to be explained. In its first phase, the Colorectal Transdisciplinary (CORECT) study, in collaboration with the Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO), the Colon Cancer Family Registry (CCFR), the Asian Colorectal Cancer Consortium (ACCC) and the Multiethnic Cohort (MEC), identified 6 of these loci based on 18,299 cases and 19,656 controls. The CORECT consortium recently completed genotyping on approximately 26,000 cases and 18,000 controls using the Illumina Infinium® OncoArray platform with the goals of 1) identifying novel overall and ethnic-specific susceptibility alleles and 2) fine-mapping known CRC risk loci. Newly genotyped samples represent multiple ethnic groups and include individuals primarily of European (75%), Asian (19%) and Hispanic (2.5%) origins. Combining both phases of CORECT, this study constitutes the largest GWAS meta-analysis of CRC to date. Here, we present European-specific results derived from genetically-defined Europeans genotyped on the OncoArray (16,456 cases and 10,442 controls), CORECT Phase I (5,584 cases and 5,329 controls), and CCFR/GECCO (12,715 cases and 14,327 controls), for a total of 34,755 cases and 30,098 controls. OncoArray genotype data were imputed to the 1000 Genomes Phase III reference panel, and summary results from logistic regression adjusting for age, sex, global ancestry, and study-specific covariates were combined in a fixed-effects inverse variance-weighted meta-analysis. Preliminary results indicate at least 10 new low-penetrance risk variants that reach genome-wide significance (P&lt;5E-08) and that are independent of known risk loci. Further, fine-mapping of well-established CRC susceptibility regions is underway. This investigation provides additional insight into the etiology of CRC and informs future risk modeling efforts.

Citation Format: Stephanie L. Schmit, Fredrick R. Schumacher, Christopher K. Edlund, Jian Gong, Gad Rennert, Wei Zheng, Loic Le Marchand, Ulrike Peters, Graham Casey, Li Hsu, Stephen B. Gruber, David V. Conti, CORECT, GECCO, CCFR, ACCC, and MEC. Novel susceptibility loci for colorectal cancer: Findings from the colorectal transdisciplinary (CORECT) study OncoArray analysis. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr LB-365.

Genome-Wide Association of the Laboratory-Based Nicotine Metabolite Ratio in Three Ancestries

Introduction:

Metabolic enzyme variation and other patient and environmental characteristics influence smoking behaviors, treatment success, and risk of related disease. Population-specific variation in metabolic genes contributes to challenges in developing and optimizing pharmacogenetic interventions. We applied a custom genome-wide genotyping array for addiction research (Smokescreen), to three laboratory-based studies of nicotine metabolism with oral or venous administration of labeled nicotine and cotinine, to model nicotine metabolism in multiple populations. The trans-3′-hydroxycotinine/cotinine ratio, the nicotine metabolite ratio (NMR), was the nicotine metabolism measure analyzed.

Methods:

Three hundred twelve individuals of self-identified European, African, and Asian American ancestry were genotyped and included in ancestry-specific genome-wide association scans (GWAS) and a meta-GWAS analysis of the NMR. We modeled natural-log transformed NMR with covariates: principal components of genetic ancestry, age, sex, body mass index, and smoking status.

Results:

African and Asian American NMRs were statistically significantly (P values ≤ 5E-5) lower than European American NMRs. Meta-GWAS analysis identified 36 genome-wide significant variants over a 43 kilobase pair region at CYP2A6 with minimum P = 2.46E-18 at rs12459249, proximal to CYP2A6. Additional minima were located in intron 4 (rs56113850, P = 6.61E-18) and in the CYP2A6-CYP2A7 intergenic region (rs34226463, P = 1.45E-12). Most (34/36) genome-wide significant variants suggested reduced CYP2A6 activity; functional mechanisms were identified and tested in knowledge-bases. Conditional analysis resulted in intergenic variants of possible interest (P values < 5E-5).

Conclusions:

This meta-GWAS of the NMR identifies CYP2A6 variants, replicates the top-ranked single nucleotide polymorphism from a recent Finnish meta-GWAS of the NMR, identifies functional mechanisms, and provides pan-continental population biomarkers for nicotine metabolism.

Implications:

This multiple ancestry meta-GWAS of the laboratory study-based NMR provides novel evidence and replication for genome-wide association of CYP2A6 single nucleotide and insertion–deletion polymorphisms. We identify three regions of genome-wide significance: proximal, intronic, and distal to CYP2A6. We replicate the top-ranking single nucleotide polymorphism from a recent GWAS of the NMR in Finnish smokers, identify a functional mechanism for this intronic variant from in silico analyses of RNA-seq data that is consistent with CYP2A6 expression measured in postmortem lung and liver, and provide additional support for the intergenic region between CYP2A6 and CYP2A7.

Genome-wide association study of colorectal cancer in Hispanics

This manuscript describes the first large-scale genome-wide association study of colorectal cancer in Hispanics and Latinos. Our results demonstrate the broad replication of known susceptibility regions and the importance of fine-mapping in ethnic minority populations.

Genome-wide association studies (GWAS) have identified 58 susceptibility alleles across 37 regions associated with the risk of colorectal cancer (CRC) with P < 5×10⁻⁸. Most studies have been conducted in non-Hispanic whites and East Asians; however, the generalizability of these findings and the potential for ethnic-specific risk variation in Hispanic and Latino (HL) individuals have been largely understudied. We describe the first GWAS of common genetic variation contributing to CRC risk in HL (1611 CRC cases and 4330 controls). We also examine known susceptibility alleles and implement imputation-based fine-mapping to identify potential ethnicity-specific association signals in known risk regions. We discovered 17 variants across 4 independent regions that merit further investigation due to suggestive CRC associations (P < 1×10⁻⁶) at 1p34.3 (rs7528276; Odds Ratio (OR) = 1.86 [95% confidence interval (CI): 1.47–2.36); P = 2.5×10⁻⁷], 2q23.3 (rs1367374; OR = 1.37 (95% CI: 1.21–1.55); P = 4.0×10⁻⁷), 14q24.2 (rs143046984; OR = 1.65 (95% CI: 1.36–2.01); P = 4.1×10⁻⁷) and 16q12.2 [rs142319636; OR = 1.69 (95% CI: 1.37–2.08); P=7.8×10⁻⁷]. Among the 57 previously published CRC susceptibility alleles with minor allele frequency ≥1%, 76.5% of SNPs had a consistent direction of effect and 19 (33.3%) were nominally statistically significant (P < 0.05). Further, rs185423955 and rs60892987 were identified as novel secondary susceptibility variants at 3q26.2 (P = 5.3×10^–5) and 11q12.2 (P = 6.8×10⁻⁵), respectively. Our findings demonstrate the importance of fine mapping in HL. These results are informative for variant prioritization in functional studies and future risk prediction modeling in minority populations.

GWASeq: targeted re-sequencing follow up to GWAS

BACKGROUND

For the last decade the conceptual framework of the Genome-Wide Association Study (GWAS) has dominated the investigation of human disease and other complex traits. While GWAS have been successful in identifying a large number of variants associated with various phenotypes, the overall amount of heritability explained by these variants remains small. This raises the question of how best to follow up on a GWAS, localize causal variants accounting for GWAS hits, and as a consequence explain more of the so-called "missing" heritability. Advances in high throughput sequencing technologies now allow for the efficient and cost-effective collection of vast amounts of fine-scale genomic data to complement GWAS.

RESULTS

We investigate these issues using a colon cancer dataset. After QC, our data consisted of 1993 cases, 899 controls. Using marginal tests of associations, we identify 10 variants distributed among six targeted regions that are significantly associated with colorectal cancer, with eight of the variants being novel to this study. Additionally, we perform so-called 'SNP-set' tests of association and identify two sets of variants that implicate both common and rare variants in the etiology of colorectal cancer.

CONCLUSIONS

Here we present a large-scale targeted re-sequencing resource focusing on genomic regions implicated in colorectal cancer susceptibility previously identified in several GWAS, which aims to 1) provide fine-scale targeted sequencing data for fine-mapping and 2) provide data resources to address methodological questions regarding the design of sequencing-based follow-up studies to GWAS. Additionally, we show that this strategy successfully identifies novel variants associated with colorectal cancer susceptibility and can implicate both common and rare variants.

Smokescreen: a targeted genotyping array for addiction research

BACKGROUND

Addictive disorders are a class of chronic, relapsing mental disorders that are responsible for increased risk of mental and medical disorders and represent the largest, potentially modifiable cause of death. Tobacco dependence is associated with increased risk of disease and premature death. While tobacco control efforts and therapeutic interventions have made good progress in reducing smoking prevalence, challenges remain in optimizing their effectiveness based on patient characteristics, including genetic variation. In order to maximize collaborative efforts to advance addiction research, we have developed a genotyping array called Smokescreen. This custom array builds upon previous work in the analyses of human genetic variation, the genetics of addiction, drug metabolism, and response to therapy, with an emphasis on smoking and nicotine addiction.

RESULTS

The Smokescreen genotyping array includes 646,247 markers in 23 categories. The array design covers genome-wide common variation (65.67, 82.37, and 90.72% in African (YRI), East Asian (ASN), and European (EUR) respectively); most of the variation with a minor allele frequency ≥ 0.01 in 1014 addiction genes (85.16, 89.51, and 90.49% for YRI, ASN, and EUR respectively); and nearly all variation from the 1000 Genomes Project Phase 1, NHLBI GO Exome Sequencing Project and HapMap databases in the regions related to smoking behavior and nicotine metabolism: CHRNA5-CHRNA3-CHRNB4 and CYP2A6-CYP2B6. Of the 636 pilot DNA samples derived from blood or cell line biospecimens that were genotyped on the array, 622 (97.80%) passed quality control. In passing samples, 90.08% of markers passed quality control. The genotype reproducibility in 25 replicate pairs was 99.94%. For 137 samples that overlapped with HapMap2 release 24, the genotype concordance was 99.76%. In a genome-wide association analysis of the nicotine metabolite ratio in 315 individuals participating in nicotine metabolism laboratory studies, we identified genome-wide significant variants in the CYP2A6 region (min p = 9.10E-15).

CONCLUSIONS

We developed a comprehensive genotyping array for addiction research and demonstrated its analytic validity and utility through pilot genotyping of HapMap and study samples. This array allows researchers to perform genome-wide, candidate gene, and pathway-based association analyses of addiction, tobacco-use, treatment response, comorbidities, and associated diseases in a standardized, high-throughput platform.

Germline Mutations in the BRIP1, BARD1, PALB2, and NBN Genes in Women With Ovarian Cancer

BACKGROUND

Epithelial ovarian cancer (EOC) is the most lethal gynecological malignancy, responsible for 13 000 deaths per year in the United States. Risk prediction based on identifying germline mutations in ovarian cancer susceptibility genes could have a clinically significant impact on reducing disease mortality.

METHODS

Next generation sequencing was used to identify germline mutations in the coding regions of four candidate susceptibility genes-BRIP1, BARD1, PALB2 and NBN-in 3236 invasive EOC case patients and 3431 control patients of European origin, and in 2000 unaffected high-risk women from a clinical screening trial of ovarian cancer (UKFOCSS). For each gene, we estimated the prevalence and EOC risks and evaluated associations between germline variant status and clinical and epidemiological risk factor information. All statistical tests were two-sided.

RESULTS

We found an increased frequency of deleterious mutations in BRIP1 in case patients (0.9%) and in the UKFOCSS participants (0.6%) compared with control patients (0.09%) (P = 1 x 10(-4) and 8 x 10(-4), respectively), but no differences for BARD1 (P = .39), NBN1 ( P = .61), or PALB2 (P = .08). There was also a difference in the frequency of rare missense variants in BRIP1 between case patients and control patients (P = 5.5 x 10(-4)). The relative risks associated with BRIP1 mutations were 11.22 for invasive EOC (95% confidence interval [CI] = 3.22 to 34.10, P = 1 x 10(-4)) and 14.09 for high-grade serous disease (95% CI = 4.04 to 45.02, P = 2 x 10(-5)). Segregation analysis in families estimated the average relative risks in BRIP1 mutation carriers compared with the general population to be 3.41 (95% CI = 2.12 to 5.54, P = 7×10(-7)).

CONCLUSIONS

Deleterious germline mutations in BRIP1 are associated with a moderate increase in EOC risk. These data have clinical implications for risk prediction and prevention approaches for ovarian cancer and emphasize the critical need for risk estimates based on very large sample sizes before genes of moderate penetrance have clinical utility in cancer prevention.

Contribution of Germline Mutations in the RAD51B, RAD51C, and RAD51D Genes to Ovarian Cancer in the Population

PURPOSE

The aim of this study was to estimate the contribution of deleterious mutations in the RAD51B, RAD51C, and RAD51D genes to invasive epithelial ovarian cancer (EOC) in the population and in a screening trial of individuals at high risk of ovarian cancer.

PATIENTS AND METHODS

The coding sequence and splice site boundaries of the three RAD51 genes were sequenced and analyzed in germline DNA from a case-control study of 3,429 patients with invasive EOC and 2,772 controls as well as in 2,000 unaffected women who were BRCA1/BRCA2 negative from the United Kingdom Familial Ovarian Cancer Screening Study (UK_FOCSS) after quality-control analysis.

RESULTS

In the case-control study, we identified predicted deleterious mutations in 28 EOC cases (0.82%) compared with three controls (0.11%; P < .001). Mutations in EOC cases were more frequent in RAD51C (14 occurrences, 0.41%) and RAD51D (12 occurrences, 0.35%) than in RAD51B (two occurrences, 0.06%). RAD51C mutations were associated with an odds ratio of 5.2 (95% CI, 1.1 to 24; P = .035), and RAD51D mutations conferred an odds ratio of 12 (95% CI, 1.5 to 90; P = .019). We identified 13 RAD51 mutations (0.65%) in unaffected UK_FOCSS participants (RAD51C, n = 7; RAD51D, n = 5; and RAD51B, n = 1), which was a significantly greater rate than in controls (P < .001); furthermore, RAD51 mutation carriers were more likely than noncarriers to have a family history of ovarian cancer (P < .001).

CONCLUSION

These results confirm that RAD51C and RAD51D are moderate ovarian cancer susceptibility genes and suggest that they confer levels of risk of EOC that may warrant their use alongside BRCA1 and BRCA2 in routine clinical genetic testing.

Identification of a Novel Mucin Gene HCG22 Associated With Steroid-Induced Ocular Hypertension

PURPOSE

The pathophysiology of ocular hypertension (OH) leading to primary open-angle glaucoma shares many features with a secondary form of OH caused by treatment with glucocorticoids, but also exhibits distinct differences. In this study, a pharmacogenomics approach was taken to discover candidate genes for this disorder.

METHODS

A genome-wide association study was performed, followed by an independent candidate gene study, using a cohort enrolled from patients treated with off-label intravitreal triamcinolone, and handling change in IOP as a quantitative trait.

RESULTS

An intergenic quantitative trait locus (QTL) was identified at chromosome 6p21.33 near the 5' end of HCG22 that attained the accepted statistical threshold for genome-level significance. The HCG22 transcript, encoding a novel mucin protein, was expressed in trabecular meshwork cells, and expression was stimulated by IL-1, and inhibited by triamcinolone acetate and TGF-β. Bioinformatic analysis defined the QTL as an approximately 4 kilobase (kb) linkage disequilibrium block containing 10 common single nucleotide polymorphisms (SNPs). Four of these SNPs were identified in the National Center for Biotechnology Information (NCBI) GTEx eQTL browser as modifiers of HCG22 expression. Most are predicted to disrupt or improve motifs for transcription factor binding, the most relevant being disruption of the glucocorticoid receptor binding motif. A second QTL was identified within the predicted signal peptide of the HCG22 encoded protein that could affect its secretion. Translation, O-glycosylation, and secretion of the predicted HCG22 protein was verified in cultured trabecular meshwork cells.

CONCLUSIONS

Identification of two independent QTLs that could affect expression of the HCG22 mucin gene product via two different mechanisms (transcription or secretion) is highly suggestive of a role in steroid-induced OH.

15q12 variants, sputum gene promoter hypermethylation, and lung cancer risk: a GWAS in smokers

BACKGROUND

Lung cancer is the leading cause of cancer-related mortality worldwide. Detection of promoter hypermethylation of tumor suppressor genes in exfoliated cells from the lung provides an assessment of field cancerization that in turn predicts lung cancer. The identification of genetic determinants for this validated cancer biomarker should provide novel insights into mechanisms underlying epigenetic reprogramming during lung carcinogenesis.

METHODS

A genome-wide association study using generalized estimating equations and logistic regression models was conducted in two geographically independent smoker cohorts to identify loci affecting the propensity for cancer-related gene methylation that was assessed by a 12-gene panel interrogated in sputum. All statistical tests were two-sided.

RESULTS

Two single nucleotide polymorphisms (SNPs) at 15q12 (rs73371737 and rs7179575) that drove gene methylation were discovered and replicated with rs73371737 reaching genome-wide significance (P = 3.3×10(-8)). A haplotype carrying risk alleles from the two 15q12 SNPs conferred 57% increased risk for gene methylation (P = 2.5×10(-9)). Rs73371737 reduced GABRB3 expression in lung cells and increased risk for smoking-induced chronic mucous hypersecretion. Furthermore, subjects with variant homozygote of rs73371737 had a two-fold increase in risk for lung cancer (P = .0043). Pathway analysis identified DNA double-strand break repair by homologous recombination (DSBR-HR) as a major pathway affecting susceptibility for gene methylation that was validated by measuring chromatid breaks in lymphocytes challenged by bleomycin.

CONCLUSIONS

A functional 15q12 variant was identified as a risk factor for gene methylation and lung cancer. The associations could be mediated by GABAergic signaling that drives the smoking-induced mucous cell metaplasia. Our findings also substantiate DSBR-HR as a critical pathway driving epigenetic gene silencing.

Ethnic-specific associations of rare and low-frequency DNA sequence variants with asthma

Common variants at many loci have been robustly associated with asthma but explain little of the overall genetic risk. Here we investigate the role of rare (<1%) and low-frequency (1-5%) variants using the Illumina HumanExome BeadChip array in 4,794 asthma cases, 4,707 non-asthmatic controls and 590 case-parent trios representing European Americans, African Americans/African Caribbeans and Latinos. Our study reveals one low-frequency missense mutation in the GRASP gene that is associated with asthma in the Latino sample (P=4.31 × 10(-6); OR=1.25; MAF=1.21%) and two genes harbouring functional variants that are associated with asthma in a gene-based analysis: GSDMB at the 17q12-21 asthma locus in the Latino and combined samples (P=7.81 × 10(-8) and 4.09 × 10(-8), respectively) and MTHFR in the African ancestry sample (P=1.72 × 10(-6)). Our results suggest that associations with rare and low-frequency variants are ethnic specific and not likely to explain a significant proportion of the 'missing heritability' of asthma.

Cross-disorder genome-wide analyses suggest a complex genetic relationship between Tourette's syndrome and OCD

OBJECTIVE

Obsessive-compulsive disorder (OCD) and Tourette's syndrome are highly heritable neurodevelopmental disorders that are thought to share genetic risk factors. However, the identification of definitive susceptibility genes for these etiologically complex disorders remains elusive. The authors report a combined genome-wide association study (GWAS) of Tourette's syndrome and OCD.

METHOD

The authors conducted a GWAS in 2,723 cases (1,310 with OCD, 834 with Tourette's syndrome, 579 with OCD plus Tourette's syndrome/chronic tics), 5,667 ancestry-matched controls, and 290 OCD parent-child trios. GWAS summary statistics were examined for enrichment of functional variants associated with gene expression levels in brain regions. Polygenic score analyses were conducted to investigate the genetic architecture within and across the two disorders.

RESULTS

Although no individual single-nucleotide polymorphisms (SNPs) achieved genome-wide significance, the GWAS signals were enriched for SNPs strongly associated with variations in brain gene expression levels (expression quantitative loci, or eQTLs), suggesting the presence of true functional variants that contribute to risk of these disorders. Polygenic score analyses identified a significant polygenic component for OCD (p=2×10(-4)), predicting 3.2% of the phenotypic variance in an independent data set. In contrast, Tourette's syndrome had a smaller, nonsignificant polygenic component, predicting only 0.6% of the phenotypic variance (p=0.06). No significant polygenic signal was detected across the two disorders, although the sample is likely underpowered to detect a modest shared signal. Furthermore, the OCD polygenic signal was significantly attenuated when cases with both OCD and co-occurring Tourette's syndrome/chronic tics were included in the analysis (p=0.01).

CONCLUSIONS

Previous work has shown that Tourette's syndrome and OCD have some degree of shared genetic variation. However, the data from this study suggest that there are also distinct components to the genetic architectures of these two disorders. Furthermore, OCD with co-occurring Tourette's syndrome/chronic tics may have different underlying genetic susceptibility compared with OCD alone.

Multiple functional risk variants in a SMAD7 enhancer implicate a colorectal cancer risk haplotype

Genome-wide association studies (GWAS) of colorectal cancer (CRC) have led to the identification of a number of common variants associated with modest risk. Several risk variants map within the vicinity of TGFβ/BMP signaling pathway genes, including rs4939827 within an intron of SMAD7 at 18q21.1. A previous study implicated a novel SNP (novel 1 or rs58920878) as a functional variant within an enhancer element in SMAD7 intron 4. In this study, we show that four SNPs including novel 1 (rs6507874, rs6507875, rs8085824, and rs58920878) in linkage disequilibrium (LD) with the index SNP rs4939827 demonstrate allele-specific enhancer effects in a large, multi-component enhancer of SMAD7. All four SNPs demonstrate allele-specific protein binding to nuclear extracts of CRC cell lines. Furthermore, some of the risk-associated alleles correlate with increased expression of SMAD7 in normal colon tissues. Finally, we show that the enhancer is responsive to BMP4 stimulation. Taken together, we propose that the associated CRC risk at 18q21.1 is due to four functional variants that regulate SMAD7 expression and potentially perturb a BMP negative feedback loop in TGFβ/BMP signaling pathways.

Organic cation transporter variation and response to smoking cessation therapies

INTRODUCTION

We evaluated chr6q25.3 organic cation transporter gene (SLC22A1, SLC22A2, SLC22A3) variation and response to smoking cessation therapies. The corresponding proteins are low-affinity transporters of choline, acetylcholine and monoamines, and smoking cessation pharmacotherapies expressed in multiple tissues.

METHODS

We selected 7 common polymorphisms for mega-regression analysis. We assessed additive model association of polymorphisms with 7-day point prevalence abstinence overall and by assigned pharmacotherapy at end of treatment and at 6 months among European-ancestry participants of 7 randomized controlled trials adjusted for demographic, population genetic, and trial covariates.

RESULTS

Initial results were obtained in 6 trials with 1,839 participants. Nominally statistically significant associations of 2 SLC22A2 polymorphisms were observed: (1) with rs316019 at 6 months, overall ([c.808T>G; p.Ser270Ala], OR = 1.306, 95% CI = 1.034-1.649, p = .025), and among those randomized to nicotine replacement therapy (NRT) (OR = 1.784, 95% CI = 1.072-2.970, p = .026); and (2) with rs316006 (c.1502-529A>T) among those randomized to varenicline (OR = 1.420, 95% CI = 1.038-1.944, p = .028, OR = 1.362, 95% CI = 1.001-1.853, p = .04) at end of treatment and 6 months. Individuals randomized to NRT from a seventh trial were genotyped for rs316019; rs316019 was associated with a nominally statistically significant effect on abstinence overall at 6 months among 2,233 participants (OR = 1.249, 95% CI = 1.007-1.550, p = .043).

CONCLUSIONS

The functional OCT2 Ser270Ala polymorphism is nominally statistically significantly associated with abstinence among European-ancestry treatment-seeking smokers after adjustments for pharmacotherapy, demographics, population genetics, and without adjustment for multiple testing of 7 SNPs. Replication of these preliminary findings in additional randomized controlled trials of smoking cessation therapies and from multiple continental populations would describe another pharmacogenetic role for SLC22A2/OCT2.

A novel colorectal cancer risk locus at 4q32.2 identified from an international genome-wide association study

Only a fraction of colorectal cancer heritability is explained by known risk-conferring genetic variation. This study was designed to identify novel risk alleles in Europeans. We conducted a genome-wide association study (GWAS) meta-analysis of colorectal cancer in participants from a population-based case-control study in Israel (n = 1616 cases, 1329 controls) and a consortium study from the Colon Cancer Family Registry (n = 1977 cases, 999 controls). We used a two-stage (discovery-replication) GWAS design, followed by a joint meta-analysis. A combined analysis identified a novel susceptibility locus that reached genome-wide significance on chromosome 4q32.2 [rs35509282, risk allele = A (minor allele frequency = 0.09); odds ratio (OR) per risk allele = 1.53; P value = 8.2 × 10(-9); nearest gene = FSTL5]. The direction of the association was consistent across studies. In addition, we confirmed that 14 of 29 previously identified susceptibility variants were significantly associated with risk of colorectal cancer in this study. Genetic variation on chromosome 4q32.2 is significantly associated with risk of colorectal cancer in Ashkenazi Jews and Europeans in this study.

A meta-analysis of Hodgkin lymphoma reveals 19p13.3 TCF3 as a novel susceptibility locus

Recent genome-wide association studies (GWAS) of Hodgkin lymphoma (HL) have identified associations with genetic variation at both HLA and non-HLA loci; however, much of heritable HL susceptibility remains unexplained. Here we perform a meta-analysis of three HL GWAS totaling 1,816 cases and 7,877 controls followed by replication in an independent set of 1,281 cases and 3,218 controls to find novel risk loci. We identify a novel variant at 19p13.3 associated with HL (rs1860661; odds ratio (OR)=0.81, 95% confidence interval (95% CI) = 0.76-0.86, P(combined) = 3.5 × 10(-10)), located in intron 2 of TCF3 (also known as E2A), a regulator of B- and T-cell lineage commitment known to be involved in HL pathogenesis. This meta-analysis also notes associations between previously published loci at 2p16, 5q31, 6p31, 8q24 and 10p14 and HL subtypes. We conclude that our data suggest a link between the 19p13.3 locus, including TCF3, and HL risk.

Comprehensive analyses of DNA repair pathways, smoking and bladder cancer risk in Los Angeles and Shanghai

Tobacco smoking is a bladder cancer risk factor and a source of carcinogens that induce DNA damage to urothelial cells. Using data and samples from 988 cases and 1,004 controls enrolled in the Los Angeles County Bladder Cancer Study and the Shanghai Bladder Cancer Study, we investigated associations between bladder cancer risk and 632 tagSNPs that comprehensively capture genetic variation in 28 DNA repair genes from four DNA repair pathways: base excision repair (BER), nucleotide excision repair (NER), non-homologous end-joining (NHEJ) and homologous recombination repair (HHR). Odds ratios (ORs) and 95% confidence intervals (CIs) for each tagSNP were corrected for multiple testing for all SNPs within each gene using pACT and for genes within each pathway and across pathways with Bonferroni. Gene and pathway summary estimates were obtained using ARTP. We observed an association between bladder cancer and POLB rs7832529 (BER) (pACT = 0.003; ppathway = 0.021) among all, and SNPs in XPC (NER) and OGG1 (BER) among Chinese men and women, respectively. The NER pathway showed an overall association with risk among Chinese males (ARTP NER p = 0.034). The XRCC6 SNP rs2284082 (NHEJ), also in LD with SREBF2, showed an interaction with smoking (smoking status interaction pgene = 0.001, ppathway = 0.008, poverall = 0.034). Our findings support a role in bladder carcinogenesis for regions that map close to or within BER (POLB, OGG1) and NER genes (XPC). A SNP that tags both the XRCC6 and SREBF2 genes strongly modifies the association between bladder cancer risk and smoking.

Identification and characterization of functional risk variants for colorectal cancer mapping to chromosome 11q23.1

Genome-wide association studies of colorectal cancer (CRC) have identified a number of common variants associated with modest risk, including rs3802842 at chromosome 11q23.1. Several genes map to this region but rs3802842 does not map to any known transcribed or regulatory sequences. We reasoned, therefore, that rs3802842 is not the functional single-nucleotide polymorphism (SNP), but is in linkage disequilibrium (LD) with a functional SNP(s). We performed ChIP-seq for histone modifications in SW480 and HCT-116 CRC cells, and incorporated ChIP-seq and DNase I hypersensitivity data available through ENCODE within a 137-kb genomic region containing rs3802842 on 11q23.1. We identified SNP rs10891246 in LD with rs3802842 that mapped within a bidirectional promoter region of genes C11orf92 and C11orf93. Following mutagenesis to the risk allele, the promoter demonstrated lower levels of reporter gene expression. A second SNP rs7130173 was identified in LD with rs3802842 that mapped to a candidate enhancer region, which showed strong unidirectional activity in both HCT-116 and SW480 CRC cells. The risk allele of rs7130173 demonstrated reduced enhancer activity compared with the common allele, and reduced nuclear protein binding affinity in electromobility shift assays compared with the common allele suggesting differential transcription factor (TF) binding. SNPs rs10891246 and rs7130173 are on the same haplotype, and expression quantitative trait loci (eQTL) analyses of neighboring genes implicate C11orf53, C11orf92 and C11orf93 as candidate target genes. These data imply that rs10891246 and rs7130173 are functional SNPs mapping to 11q23.1 and that C11orf53, C11orf92 and C11orf93 represent novel candidate target genes involved in CRC etiology.

The sex hormone system in carriers of BRCA1/2 mutations: a case-control study

BACKGROUND

Penetrance for breast cancer, ovarian cancer, or both in carriers of BRCA1/BRCA2 mutations is disproportionately high. Sex hormone dysregulation and altered end-organ hormone sensitivity might explain this organ-specific penetrance. We sought to identify differences in hormone regulation between carriers of BRCA1/2 and women who are negative for BRCA1/2 mutations.

METHODS

We assessed endometrial thickness for each menstrual cycle day (as an index of hormone regulation) in 393 scans from 228 women in the UK Familial Ovarian Cancer Screening Study (UK FOCSS) known to carry either mutation and 1573 scans from 754 women known to be negative for the mutations. To quantify differences in endometrial thickness we focused on days 10-14 and days 21-26, and calculated the area under the curve. We then compared serum oestradiol and progesterone titres during these days of the menstrual cycle in the same groups. Follicular and luteal oestradiol and progesterone serum titres were grouped into quartiles and odds ratios were calculated with logistic regression.

FINDINGS

Follicular phase endometrial thickness of carriers of the mutations adjusted for age and day of the menstrual cycle was higher (odds ratio [OR] 1·11, 95% CI 1·03-1·20; p=0·0063) and luteal phase endometrial thickness lower (0·90, 0·83-0·98; p=0·027) than for women negative for the mutations. Median luteal phase titres of progesterone were 121% higher (p=0·00037) in carriers than in women negative for the mutations, and for oestradiol were 33% higher (p=0·007)-ie, 59% of carriers had concentrations of serum progesterone that would have been in the top quartile of concentrations in the control group (OR 8·0, 95% CI 2·1-52·57; p=0·008).

INTERPRETATION

Carriers of BRCA1/BRCA2 mutations are exposed to higher titres of oestradiol and progesterone-known risk-factors for breast cancer. Higher titres of oestradiol in carriers are compatible with this hormone having a role in ovarian carcinogenesis in such women. Our findings could not be explained by differential contraceptive pill use.

Partitioning the heritability of Tourette syndrome and obsessive compulsive disorder reveals differences in genetic architecture

The direct estimation of heritability from genome-wide common variant data as implemented in the program Genome-wide Complex Trait Analysis (GCTA) has provided a means to quantify heritability attributable to all interrogated variants. We have quantified the variance in liability to disease explained by all SNPs for two phenotypically-related neurobehavioral disorders, obsessive-compulsive disorder (OCD) and Tourette Syndrome (TS), using GCTA. Our analysis yielded a heritability point estimate of 0.58 (se = 0.09, p = 5.64e-12) for TS, and 0.37 (se = 0.07, p = 1.5e-07) for OCD. In addition, we conducted multiple genomic partitioning analyses to identify genomic elements that concentrate this heritability. We examined genomic architectures of TS and OCD by chromosome, MAF bin, and functional annotations. In addition, we assessed heritability for early onset and adult onset OCD. Among other notable results, we found that SNPs with a minor allele frequency of less than 5% accounted for 21% of the TS heritability and 0% of the OCD heritability. Additionally, we identified a significant contribution to TS and OCD heritability by variants significantly associated with gene expression in two regions of the brain (parietal cortex and cerebellum) for which we had available expression quantitative trait loci (eQTLs). Finally we analyzed the genetic correlation between TS and OCD, revealing a genetic correlation of 0.41 (se = 0.15, p = 0.002). These results are very close to previous heritability estimates for TS and OCD based on twin and family studies, suggesting that very little, if any, heritability is truly missing (i.e., unassayed) from TS and OCD GWAS studies of common variation. The results also indicate that there is some genetic overlap between these two phenotypically-related neuropsychiatric disorders, but suggest that the two disorders have distinct genetic architectures.

Genome-wide association study of obsessive-compulsive disorder

Obsessive-compulsive disorder (OCD) is a common, debilitating neuropsychiatric illness with complex genetic etiology. The International OCD Foundation Genetics Collaborative (IOCDF-GC) is a multi-national collaboration established to discover the genetic variation predisposing to OCD. A set of individuals affected with DSM-IV OCD, a subset of their parents, and unselected controls, were genotyped with several different Illumina SNP microarrays. After extensive data cleaning, 1465 cases, 5557 ancestry-matched controls and 400 complete trios remained, with a common set of 469,410 autosomal and 9657 X-chromosome single nucleotide polymorphisms (SNPs). Ancestry-stratified case-control association analyses were conducted for three genetically-defined subpopulations and combined in two meta-analyses, with and without the trio-based analysis. In the case-control analysis, the lowest two P-values were located within DLGAP1 (P=2.49 × 10(-6) and P=3.44 × 10(-6)), a member of the neuronal postsynaptic density complex. In the trio analysis, rs6131295, near BTBD3, exceeded the genome-wide significance threshold with a P-value=3.84 × 10(-8). However, when trios were meta-analyzed with the case-control samples, the P-value for this variant was 3.62 × 10(-5), losing genome-wide significance. Although no SNPs were identified to be associated with OCD at a genome-wide significant level in the combined trio-case-control sample, a significant enrichment of methylation QTLs (P<0.001) and frontal lobe expression quantitative trait loci (eQTLs) (P=0.001) was observed within the top-ranked SNPs (P<0.01) from the trio-case-control analysis, suggesting these top signals may have a broad role in gene expression in the brain, and possibly in the etiology of OCD.

Genetic variation in insulin pathway genes and distal colorectal adenoma risk

BACKGROUND

Insulin, glucose, and other insulin-related proteins that mediate insulin signaling are associated with colorectal neoplasia risk, but associations with common genetic variation in insulin axis genes are less clear. In this study, we used a comprehensive tag single-nucleotide polymorphisms (SNPs) approach to define genetic variation in six insulin axis genes (IGF1, IGF2, IGFBP1, IGFBP3, IRS1, and IRS2) and three genes associated with estrogen signaling (ESR1, ESR2, and PGR).

METHODS

We assessed associations between SNPs and distal colorectal adenoma (CRA) risk in a case-control study of 1,351 subjects. Cases were individuals with one or more adenomas diagnosed during sigmoidoscopy, and controls were individuals with no adenomas at the sigmoidoscopy exam. We used unconditional logistic regression assuming an additive model to assess SNP-specific risks adjusting for multiple comparisons with P (act).

RESULTS

Distal adenoma risk was significantly increased for one SNP in IGF2 [per minor allele OR = 1.41; 95 % confidence interval (CI) = 1.16, 1.67; P (act) = 0.005] and decreased for an ESR2 SNP (per minor allele OR = 0.78; 95 % CI = 0.66, 0.91; P (act) = 0.041). There was no statistically significant heterogeneity of these associations by race, sex, BMI, physical activity, or, in women, hormone replacement therapy use. Risk estimates did not differ in the colon versus rectum or for smaller (<1 cm) versus larger (>1 cm) adenomas.

CONCLUSIONS

These data suggest that selected genetic variability in IGF2 and ESR2 may be modifiers of CRA risk.

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

Genetic determinants for promoter hypermethylation in the lungs of smokers: a candidate gene-based study

The detection of tumor suppressor gene promoter methylation in sputum-derived exfoliated cells predicts early lung cancer. Here, we identified genetic determinants for this epigenetic process and examined their biologic effects on gene regulation. A two-stage approach involving discovery and replication was used to assess the association between promoter hypermethylation of a 12-gene panel and common variation in 40 genes involved in carcinogen metabolism, regulation of methylation, and DNA damage response in members of the Lovelace Smokers Cohort (N = 1,434). Molecular validation of three identified variants was conducted using primary bronchial epithelial cells. Association of study-wide significance (P < 8.2 × 10(-5)) was identified for rs1641511, rs3730859, and rs1883264 in TP53, LIG1, and BIK, respectively. These single-nucleotide polymorphisms (SNP) were significantly associated with altered expression of the corresponding genes in primary bronchial epithelial cells. In addition, rs3730859 in LIG1 was also moderately associated with increased risk for lung cancer among Caucasian smokers. Together, our findings suggest that genetic variation in DNA replication and apoptosis pathways impacts the propensity for gene promoter hypermethylation in the aerodigestive tract of smokers. The incorporation of genetic biomarkers for gene promoter hypermethylation with clinical and somatic markers may improve risk assessment models for lung cancer.

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.

Genotype-environment interactions in microsatellite stable/microsatellite instability-low colorectal cancer: results from a genome-wide association study

BACKGROUND

Genome-wide association studies (GWAS) have led to the identification of a number of common susceptibility loci for colorectal cancer (CRC); however, none of these GWAS have considered gene-environment (G × E) interactions. Therefore, it is unclear whether current hits are modified by environmental exposures or whether there are additional hits whose effects are dependent on environmental exposures.

METHODS

We conducted a systematic search for G × E interactions using genome wide data from the Colon Cancer Family Registry that included 1,191 cases of microsatellite stable (MSS) or microsatellite instability-low (MSI-L) CRC and 999 controls genotyped using either the Illumina Human1M or Human1M-Duo BeadChip. We tested for interactions between genotypes and 14 environmental factors using 3 methods: a traditional case-control test, a case-only test, and the recently proposed 2-step method by Murcray and colleagues. All potentially significant findings were replicated in the ARCTIC Study.

RESULTS

No G × E interactions were identified that reached genome-wide significance by any of the 3 methods. When analyzing previously reported susceptibility loci, 7 significant G × E interactions were found at a 5% significance level. We investigated these 7 interactions in an independent sample and none of the interactions were replicated.

CONCLUSIONS

Identifying G × E interactions will present challenges in a GWAS setting. Our power calculations illustrate the need for larger sample sizes; however, as CRC is a heterogeneous disease, a tradeoff between increasing sample size and heterogeneity needs to be considered.

IMPACT

The results from this first genome-wide analysis of G × E in CRC identify several challenges, which may be addressed by large consortium efforts.

Variation in folate pathway genes and distal colorectal adenoma risk: a sigmoidoscopy-based case-control study

Background

Folate-associated one-carbon metabolism (FOCM) is an important pathway in colorectal neoplasia risk but data on genetic variation in this pathway are largely limited to studies of single SNPs in selected genes.

Methods

We used a comprehensive tagSNP approach to study the association between genetic variation in 11 genes in the FOCM pathway and risk of incident distal colorectal adenomas in a sigmoidoscopy-based case–control study. We included 655 cases (one or more adenomas) and 695 controls (no adenomas) recruited from one of two Kaiser Permanente clinics between 1991 and 1995. We assessed a total of 159 tagSNPs selected using Haploview Tagger as well as selected non-synonymous SNPs. We used unconditional logistic regression to model the association between SNPs and risk of distal adenomas, assuming a log-additive model.

Results

Five SNPs in the SLC19A1 (RFC1) gene: rs1051266 (G80A), rs283895, rs2236484, rs12482346, and rs2838958 were associated with adenoma risk after correction for multiple testing (all corrected p values ≤0.043). The non-synonymous SLC19A1 SNP G80A interacted significantly with the MTHFRC677T genotype (interaction p value = 0.018).

Conclusion

Our data suggest that genetic variation in SLC19A1 may modify the risk of distal colorectal adenoma.

Electronic supplementary material

The online version of this article (doi:10.1007/s10552-011-9726-7) contains supplementary material, which is available to authorized users.

Validating genetic risk associations for ovarian cancer through the international Ovarian Cancer Association Consortium

The search for genetic variants associated with ovarian cancer risk has focused on pathways including sex steroid hormones, DNA repair, and cell cycle control. The Ovarian Cancer Association Consortium (OCAC) identified 10 single-nucleotide polymorphisms (SNPs) in genes in these pathways, which had been genotyped by Consortium members and a pooled analysis of these data was conducted. Three of the 10 SNPs showed evidence of an association with ovarian cancer at P< or =0.10 in a log-additive model: rs2740574 in CYP3A4 (P=0.011), rs1805386 in LIG4 (P=0.007), and rs3218536 in XRCC2 (P=0.095). Additional genotyping in other OCAC studies was undertaken and only the variant in CYP3A4, rs2740574, continued to show an association in the replication data among homozygous carriers: OR(homozygous(hom))=2.50 (95% CI 0.54-11.57, P=0.24) with 1406 cases and 2827 controls. Overall, in the combined data the odds ratio was 2.81 among carriers of two copies of the minor allele (95% CI 1.20-6.56, P=0.017, p(het) across studies=0.42) with 1969 cases and 3491 controls. There was no association among heterozygous carriers. CYP3A4 encodes a key enzyme in oestrogen metabolism and our finding between rs2740574 and risk of ovarian cancer suggests that this pathway may be involved in ovarian carcinogenesis. Additional follow-up is warranted.

Polymorphism in the IL18 gene and epithelial ovarian cancer in non-Hispanic white women

Over 22,000 cases of ovarian cancer were diagnosed in 2007 in the United States, but only a fraction of them can be attributed to mutations in highly penetrant genes such as BRCA1. To determine whether low-penetrance genetic variants contribute to ovarian cancer risk, we genotyped 1,536 single nucleotide polymorphisms (SNP) in several candidate gene pathways in 848 epithelial ovarian cancer cases and 798 controls in the North Carolina Ovarian Cancer Study (NCO) using a customized Illumina array. The inflammation gene interleukin-18 (IL18) showed the strongest evidence for association with epithelial ovarian cancer in a gene-by-gene analysis (P = 0.002) with a <25% chance of being a false-positive finding (q value = 0.240). Using a multivariate model search algorithm over 11 IL18 tagging SNPs, we found that the association was best modeled by rs1834481. Further, this SNP uniquely tagged a significantly associated IL18 haplotype and there was an increased risk of epithelial ovarian cancer per rs1834481 allele (odds ratio, 1.24; 95% confidence interval, 1.06-1.45). In a replication stage, 12 independent studies from the Ovarian Cancer Association Consortium (OCAC) genotyped rs1834481 in an additional 5,877 cases and 7,791 controls. The fixed effects estimate per rs1834481 allele was null (odds ratio, 0.99; 95% confidence interval, 0.94-1.05) when data from the 12 OCAC studies were combined. The effect estimate remained unchanged with the addition of the initial North Carolina Ovarian Cancer Study data. This analysis shows the importance of consortia, like the OCAC, in either confirming or refuting the validity of putative findings in studies with smaller sample sizes. (Cancer Epidemiol Biomarkers Prev 2008;17(12):3567-72).

Lift/transfer and technical aids for persons with severe acquired brain injury. An inventory of problems

A questionnaire about assistive devices for lift and transfer was distributed to all individuals evaluated at the Solberga Project outpatient service unit, a regional centre for long-term follow-up of severely brain injured persons in the Greater Stockholm area. The target group was 60 severely brain injured adults, of whom 57 (30 women and 27 men) answered a questionnaire created by one of the authors. The causes of injury in the study population were trauma (n = 27), cerebrovascular accident (n = 19), anoxia (n = 10) and other (n = 1). Thirty-two persons were quadriplegic after the injury and 21 were hemiplegic. More than half (33/57) reported problems with all five defined lift and transfer situations. Most (42/44) reported problems getting in and out of their wheelchairs. Technical aids were seldom used; 24 persons did not use any aids at all. The most commonly used aid was an adjustable bed. Most of the technical aids were used when the individual also had personal assistance. The persons who managed lift and transfer by themselves used few aids. It is reasonable to assume that severely brain injured individuals would be involved in more activities if they used technical aids more readily.

Properties of person hoist spreader bars and their influence on sitting/lifting position

The purpose of this study is to analyse mechanical lifts with regard to how the spreader bar design, e.g. frontal length, height and sagittal depth, affects sitting positions and safety during lift and transfer. The seven spreader bars available on the Swedish market had two-point to four-point suspensions. Twelve subjects were photographed individually sitting in each of three types of sling: divided, one-piece with divided leg support, and one-piece. All the slings had significant effects on all specific sitting aspects. However, the spreader bar design only affected the sitting position in the hoist. A long (0.62 m) spreader bar gave a relatively upright sitting position. A deep (0.18 m) spreader bar gave a backward inclination. For all specific sitting aspects, there was a significant interaction between slings and spreader bars. This study demonstrates that the design of the spreader bar, as well as that of the sling, is of decisive importance for the sitting position during lifting.

Hoists for locomotor disability--Properties of the person-carrying section and biomechanics of sitting positions

The purpose of this study was to make an inventory of the slings available for person hoists, to describe them and to carry out biomechanical analyses of the sitting positions in different slings. On the Swedish market 18 different slings were available. Nine were photographed with a subject sitting in them. The sagittal lengths of the supporting parts of the slings varied between 0.10 m and 1.0 m. The risk of falling through was considered largest in divided slings, but even in one-piece slings with divided leg supports, the person lifted can most probably slip through inadvertently. Almost all the slings gave a backward-inclined sitting position. However, in one of the slings the trunk inclination was small and the seat was slightly forward-downward inclined, giving a risk of falling/sliding forwards. Some of the slings entailed a trunk thigh angle smaller than 90 degrees, perhaps a problem for patients with hip pain and reduced hip flexion. Given the great variety of slings available, it seems important to work for relevant instructions and information about their use.