Christopher K. Edlund

Genome-wide association study of Tourette's syndrome

Tourettes syndrome (TS) is a developmental disorder that has one of the highest familial recurrence rates among neuropsychiatric diseases with complex inheritance. However, the identification of definitive TS susceptibility genes remains elusive. Here, we report the first genome-wide association study (GWAS) of TS in 1285 cases and 4964 ancestry-matched controls of European ancestry, including two European-derived population isolates, Ashkenazi Jews from North America and Israel and French Canadians from Quebec, Canada. In a primary meta-analysis of GWAS data from these European ancestry samples, no markers achieved a genome-wide threshold of significance (P<5 × 10−8); the top signal was found in rs7868992 on chromosome 9q32 within COL27A1 (P=1.85 × 10−6). A secondary analysis including an additional 211 cases and 285 controls from two closely related Latin American population isolates from the Central Valley of Costa Rica and Antioquia, Colombia also identified rs7868992 as the top signal (P=3.6 × 10−7 for the combined sample of 1496 cases and 5249 controls following imputation with 1000 Genomes data). This study lays the groundwork for the eventual identification of common TS susceptibility variants in larger cohorts and helps to provide a more complete understanding of the full genetic architecture of this disorder.

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.

Nicotinic Acetylcholine Receptor Variation and Response to Smoking Cessation Therapies

Objective To evaluate the association of nicotinic acetylcholine receptor (nAChR) single nucleotide polymorphism (SNP) with 7-day point prevalence abstinence (abstinence) in randomized clinical trials of smoking cessation therapies in individuals grouped by pharmacotherapy randomization to inform the development of personalized smoking cessation therapy. Materials and methods We quantified association of four SNPs at three nAChRs with abstinence in eight randomized clinical trials. Participants were 2633 outpatient treatment-seeking, self-identified European ancestry individuals smoking at least 10 cigarettes/day, recruited through advertisement, prescribed pharmacotherapy, and provided with behavioral therapy. Interventions included nicotine replacement therapy (NRT), bupropion, varenicline, placebo (PLA), or combined NRT and bupropion, and five modes of group and individual behavioral therapy. Outcome measures tested in multivariate logistic regression were end of treatment and 6 month (6MO) abstinence, with demographic, behavioral, and genetic covariates. Results ‘Risk’ alleles previously associated with smoking heaviness were significantly (P<0.05) associated with reduced abstinence in the PLA pharmacotherapy group (PG) at 6MO [for rs588765, odds ratio (95% confidence interval) 0.41 (0.17–0.99)], and at end of treatment and at 6MO [for rs1051730, 0.42 (0.19–0.93) and 0.31 (0.12–0.80)], and with increased abstinence in the NRT PG at 6MO [for rs588765, 2.07 (1.11–3.87) and for rs1051730, 2.54 (1.29–4.99)]. We observed significant heterogeneity in rs1051730 effects (F=2.48, P=0.021) between PGs. Conclusion chr15q25.1 nAChR SNP risk alleles for smoking heaviness significantly increase relapse with PLA treatment and significantly increase abstinence with NRT. These SNP–PG associations require replication in independent samples for validation, and testing in larger sample sizes to evaluate whether similar effects occur in other PGs.

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.

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.

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.

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. Cancer Epidemiol Biomarkers Prev; 20(5); 758–66. ©2011 AACR.

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

BackgroundThere 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.ResultsWe 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.ConclusionSnagger 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.

Genome-wide association study of colorectal cancer identifies six new susceptibility loci

Genetic susceptibility to colorectal cancer is caused by rare pathogenic mutations and common genetic variants that contribute to familial risk. Here we report the results of a two-stage association study with 18,299 cases of colorectal cancer and 19,656 controls, with follow-up of the most statistically significant genetic loci in 4,725 cases and 9,969 controls from two Asian consortia. We describe six new susceptibility loci reaching a genome-wide threshold of P<5.0E-08. These findings provide additional insight into the underlying biological mechanisms of colorectal cancer and demonstrate the scientific value of large consortia-based genetic epidemiology studies.