Andrew Crenshaw

Multiple loci identified in a genome-wide association study of prostate cancer

We followed our initial genome-wide association study (GWAS) of 527,869 SNPs on 1,172 individuals with prostate cancer and 1,157 controls of European origin—nested in the Prostate, Lung, Colorectal, and Ovarian (PLCO) Cancer Screening Trial prospective study—by testing 26,958 SNPs in four independent studies (total of 3,941 cases and 3,964 controls). In the combined joint analysis, we confirmed three previously reported loci (two independent SNPs at 8q24 and one in HNF1B (formerly known as TCF2 on 17q); P < 10−10). In addition, loci on chromosomes 7, 10 (two loci) and 11 were highly significant (between P < 7.31 × 10−13 and P < 2.14 × 10−6). Loci on chromosome 10 include MSMB, which encodes β-microseminoprotein, a primary constituent of semen and a proposed prostate cancer biomarker, and CTBP2, a gene with antiapoptotic activity; the locus on chromosome 7 is at JAZF1, a transcriptional repressor that is fused by chromosome translocation to SUZ12 in endometrial cancer. Of the nine loci that showed highly suggestive associations (P < 2.5 × 10−5), four best fit a recessive model and included candidate susceptibility genes: CPNE3, IL16 and CDH13. Our findings point to multiple loci with moderate effects associated with susceptibility to prostate cancer that, taken together, in the future may predict high risk in select individuals.

Detectable clonal mosaicism from birth to old age and its relationship to cancer

We detected clonal mosaicism for large chromosomal anomalies (duplications, deletions and uniparental disomy) using SNP microarray data from over 50,000 subjects recruited for genome-wide association studies. This detection method requires a relatively high frequency of cells with the same abnormal karyotype (>5–10%; presumably of clonal origin) in the presence of normal cells. The frequency of detectable clonal mosaicism in peripheral blood is low (<0.5%) from birth until 50 years of age, after which it rapidly rises to 2–3% in the elderly. Many of the mosaic anomalies are characteristic of those found in hematological cancers and identify common deleted regions with genes previously associated with these cancers. Although only 3% of subjects with detectable clonal mosaicism had any record of hematological cancer before DNA sampling, those without a previous diagnosis have an estimated tenfold higher risk of a subsequent hematological cancer (95% confidence interval = 6–18).

SNP500Cancer: a public resource for sequence validation, assay development, and frequency analysis for genetic variation in candidate genes

The SNP500Cancer database provides sequence and genotype assay information for candidate SNPs useful in mapping complex diseases, such as cancer. The database is an integral component of the NCI Cancer Genome Anatomy Project (). SNP500Cancer reports sequence analysis of anonymized control DNA samples (n = 102 Coriell samples representing four self-described ethnic groups: African/African-American, Caucasian, Hispanic and Pacific Rim). The website is searchable by gene, chromosome, gene ontology pathway, dbSNP ID and SNP500Cancer SNP ID. As of October 2005, the database contains >13 400 SNPs, 9124 of which have been sequenced in the SNP500Cancer population. For each analysed SNP, gene location and >200 bp of surrounding annotated sequence (including nearby SNPs) are provided, with frequency information in total and per subpopulation as well as calculation of Hardy–Weinberg equilibrium for each subpopulation. The website provides the conditions for validated sequencing and genotyping assays, as well as genotype results for the 102 samples, in both viewable and downloadable formats. A subset of sequence validated SNPs with minor allele frequency >5% are entered into a high-throughput pipeline for genotyping analysis to determine concordance for the same 102 samples. In addition, the results of genotype analysis for select validated SNP assays (defined as 100% concordance between sequence analysis and genotype results) are posted for an additional 280 samples drawn from the Human Diversity Panel (HDP). SNP500Cancer provides an invaluable resource for investigators to select SNPs for analysis, design genotyping assays using validated sequence data, choose selected assays already validated on one or more genotyping platforms, and select reference standards for genotyping assays. The SNP500Cancer database is freely accessible via the web page at .

SNP500Cancer: a public resource for sequence validation and assay development for genetic variation in candidate genes

The SNP500Cancer Database provides sequence and genotype assay information for candidate single nucleotide polymorphisms (SNPs) useful in mapping complex diseases, such as cancer. The database is an integral component of the NCIs Cancer Genome Anatomy Project. SNP500Cancer provides bi-directional sequencing information on a set of control DNA samples derived from anonymized subjects (102 Coriell samples representing four self-described ethnic groups: African/African-American, Caucasian, Hispanic and Pacific Rim). All SNPs are chosen from public databases and reports, and the choice of genes includes a bias towards non-synonymous and promoter SNPs in genes that have been implicated in one or more cancers. The web site is searchable by gene, chromosome, gene ontology pathway and by known dbSNP ID. As of July 2003, the database contains over 3400 SNPs, 2490 of which have been sequenced in the SNP500Cancer population. For each analyzed SNP, gene location and over 200 bp of surrounding annotated sequence (including nearby SNPs) are provided, with frequency information in total and per subpopulation, and calculation of Hardy-Weinberg Equilibrium (HWE) for each subpopulation. Sequence validated SNPs with minor allele frequency > 5% are entered into a high-throughput pipeline for genotyping analysis to determine concordance for the same 102 samples. The website provides the conditions for validated genotyping assays. SNP500Cancer provides an invaluable resource for investigators to select SNPs for analysis, design genotyping assays using validated sequence data, choose selected assays already validated on one or more genotyping platforms, and select reference standards for genotyping assays. The SNP500Cancer Database is freely accessible via the web page at http://snp500cancer.nci.nih.gov/.

Common Variants at 9p21 and 8q22 Are Associated with Increased Susceptibility to Optic Nerve Degeneration in Glaucoma

Optic nerve degeneration caused by glaucoma is a leading cause of blindness worldwide. Patients affected by the normal-pressure form of glaucoma are more likely to harbor risk alleles for glaucoma-related optic nerve disease. We have performed a meta-analysis of two independent genome-wide association studies for primary open angle glaucoma (POAG) followed by a normal-pressure glaucoma (NPG, defined by intraocular pressure (IOP) less than 22 mmHg) subgroup analysis. The single-nucleotide polymorphisms that showed the most significant associations were tested for association with a second form of glaucoma, exfoliation-syndrome glaucoma. The overall meta-analysis of the GLAUGEN and NEIGHBOR dataset results (3,146 cases and 3,487 controls) identified significant associations between two loci and POAG: the CDKN2BAS region on 9p21 (rs2157719 [G], OR = 0.69 [95%CI 0.63–0.75], p = 1.86×10−18), and the SIX1/SIX6 region on chromosome 14q23 (rs10483727 [A], OR = 1.32 [95%CI 1.21–1.43], p = 3.87×10−11). In sub-group analysis two loci were significantly associated with NPG: 9p21 containing the CDKN2BAS gene (rs2157719 [G], OR = 0.58 [95% CI 0.50–0.67], p = 1.17×10−12) and a probable regulatory region on 8q22 (rs284489 [G], OR = 0.62 [95% CI 0.53–0.72], p = 8.88×10−10). Both NPG loci were also nominally associated with a second type of glaucoma, exfoliation syndrome glaucoma (rs2157719 [G], OR = 0.59 [95% CI 0.41–0.87], p = 0.004 and rs284489 [G], OR = 0.76 [95% CI 0.54–1.06], p = 0.021), suggesting that these loci might contribute more generally to optic nerve degeneration in glaucoma. Because both loci influence transforming growth factor beta (TGF-beta) signaling, we performed a genomic pathway analysis that showed an association between the TGF-beta pathway and NPG (permuted p = 0.009). These results suggest that neuro-protective therapies targeting TGF-beta signaling could be effective for multiple forms of glaucoma.

Quality Control Procedures for Genome‐Wide Association Studies

Genome‐wide association studies (GWAS) are being conducted at an unprecedented rate in population‐based cohorts and have increased our understanding of the pathophysiology of complex disease. Regardless of context, the practical utility of this information will ultimately depend upon the quality of the original data. Quality control (QC) procedures for GWAS are computationally intensive, operationally challenging, and constantly evolving. Here we enumerate some of the challenges in QC of GWAS data and describe the approaches that the electronic MEdical Records and Genomics (eMERGE) network is using for quality assurance in GWAS data, thereby minimizing potential bias and error in GWAS results. We discuss common issues associated with QC of GWAS data, including data file formats, software packages for data manipulation and analysis, sex chromosome anomalies, sample identity, sample relatedness, population substructure, batch effects, and marker quality. We propose best practices and discuss areas of ongoing and future research. Curr. Protoc. Hum. Genet. 68:1.19.1‐1.19.18

Epidermal growth factor receptor inhibition attenuates liver fibrosis and development of hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is the most rapidly increasing cause of cancer‐related mortality in the United States. Because of the lack of viable treatment options for HCC, prevention in high‐risk patients has been proposed as an alternative strategy. The main risk factor for HCC is cirrhosis and several lines of evidence implicate epidermal growth factor (EGF) in the progression of cirrhosis and development of HCC. We therefore examined the effects of the EGF receptor (EGFR) inhibitor erlotinib on liver fibrogenesis and hepatocellular transformation in three different animal models of progressive cirrhosis: a rat model induced by repeated, low‐dose injections of diethylnitrosamine (DEN), a mouse model induced by carbon tetrachloride (CCl4), and a rat model induced by bile duct ligation (BDL). Erlotinib reduced EGFR phosphorylation in hepatic stellate cells (HSC) and reduced the total number of activated HSC. Erlotinib also decreased hepatocyte proliferation and liver injury. Consistent with all these findings, pharmacological inhibition of EGFR signaling effectively prevented the progression of cirrhosis and regressed fibrosis in some animals. Moreover, by alleviating the underlying liver disease, erlotinib blocked the development of HCC and its therapeutic efficacy could be monitored with a previously reported gene expression signature predictive of HCC risk in human cirrhosis patients. Conclusion: These data suggest that EGFR inhibition using Food and Drug Administration‐approved inhibitors provides a promising therapeutic approach for reduction of fibrogenesis and prevention of HCC in high‐risk cirrhosis patients who can be identified and monitored by gene expression signatures. (Hepatology 2014;59:1577‐1590)

Candidate gene association resource (CARe): design, methods, and proof of concept.

Background— The National Heart, Lung, and Blood Institutes Candidate Gene Association Resource (CARe), a planned cross-cohort analysis of genetic variation in cardiovascular, pulmonary, hematologic, and sleep-related traits, comprises >40 000 participants representing 4 ethnic groups in 9 community-based cohorts. The goals of CARe include the discovery of new variants associated with traits using a candidate gene approach and the discovery of new variants using the genome-wide association mapping approach specifically in African Americans. Methods and Results— CARe has assembled DNA samples for >40 000 individuals self-identified as European American, African American, Hispanic, or Chinese American, with accompanying data on hundreds of phenotypes that have been standardized and deposited in the CARe Phenotype Database. All participants were genotyped for 7 single-nucleotide polymorphisms (SNPs) selected based on prior association evidence. We performed association analyses relating each of these SNPs to lipid traits, stratified by sex and ethnicity, and adjusted for age and age squared. In at least 2 of the ethnic groups, SNPs near CETP , LIPC , and LPL strongly replicated for association with high-density lipoprotein cholesterol concentrations, PCSK9 with low-density lipoprotein cholesterol levels, and LPL and APOA5 with serum triglycerides. Notably, some SNPs showed varying effect sizes and significance of association in different ethnic groups. Conclusions— The CARe Pilot Study validates the operational framework for phenotype collection, SNP genotyping, and analytic pipeline of the CARe project and validates the planned candidate gene study of ≈2000 biological candidate loci in all participants and genome-wide association study in ≈8000 African American participants. CARe will serve as a valuable resource for the scientific community.Background—The National Heart, Lung, and Blood Institute’s Candidate Gene Association Resource (CARe), a planned cross-cohort analysis of genetic variation in cardiovascular, pulmonary, hematologic, and sleep-related traits, comprises >40 000 participants representing 4 ethnic groups in 9 community-based cohorts. The goals of CARe include the discovery of new variants associated with traits using a candidate gene approach and the discovery of new variants using the genome-wide association mapping approach specifically in African Americans. Methods and Results—CARe has assembled DNA samples for >40 000 individuals self-identified as European American, African American, Hispanic, or Chinese American, with accompanying data on hundreds of phenotypes that have been standardized and deposited in the CARe Phenotype Database. All participants were genotyped for 7 single-nucleotide polymorphisms (SNPs) selected based on prior association evidence. We performed association analyses relating each of these SNPs to lipid traits, stratified by sex and ethnicity, and adjusted for age and age squared. In at least 2 of the ethnic groups, SNPs near CETP, LIPC, and LPL strongly replicated for association with high-density lipoprotein cholesterol concentrations, PCSK9 with low-density lipoprotein cholesterol levels, and LPL and APOA5 with serum triglycerides. Notably, some SNPs showed varying effect sizes and significance of association in different ethnic groups. Conclusions—The CARe Pilot Study validates the operational framework for phenotype collection, SNP genotyping, and analytic pipeline of the CARe project and validates the planned candidate gene study of ≈2000 biological candidate loci in all participants and genome-wide association study in ≈8000 African American participants. CARe will serve as a valuable resource for the scientific community.

Common variants near CAV1 and CAV2 are associated with primary open-angle glaucoma in Caucasians from the USA

Primary open-angle glaucoma (POAG) is a genetically complex common disease characterized by progressive optic nerve degeneration that results in irreversible blindness. Recently, a genome-wide association study (GWAS) for POAG in an Icelandic population identified significant associations with single nucleotide polymorphisms (SNPs) between the CAV1 and CAV2 genes on chromosome 7q31. In this study, we confirm that the identified SNPs are associated with POAG in our Caucasian US population and that specific haplotypes located in the CAV1/CAV2 intergenic region are associated with the disease. We also present data suggesting that associations with several CAV1/CAV2 SNPs are significant mostly in women.

zCall: a rare variant caller for array-based genotyping

SUMMARY zCall is a variant caller specifically designed for calling rare single-nucleotide polymorphisms from array-based technology. This caller is implemented as a post-processing step after a default calling algorithm has been applied. The algorithm uses the intensity profile of the common allele homozygote cluster to define the location of the other two genotype clusters. We demonstrate improved detection of rare alleles when applying zCall to samples that have both Illumina Infinium HumanExome BeadChip and exome sequencing data available. AVAILABILITY http://atguweb.mgh.harvard.edu/apps/zcall. CONTACT bneale@broadinstitute.org SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.