Susan Land

A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3

We describe a new computer program, SnpEff, for rapidly categorizing the effects of variants in genome sequences. Once a genome is sequenced, SnpEff annotates variants based on their genomic locations and predicts coding effects. Annotated genomic locations include intronic, untranslated region, upstream, downstream, splice site, or intergenic regions. Coding effects such as synonymous or non-synonymous amino acid replacement, start codon gains or losses, stop codon gains or losses, or frame shifts can be predicted. Here the use of SnpEff is illustrated by annotating ~356,660 candidate SNPs in ~117 Mb unique sequences, representing a substitution rate of ~1/305 nucleotides, between the Drosophila melanogaster w1118; iso-2; iso-3 strain and the reference y1; cn1 bw1 sp1 strain. We show that ~15,842 SNPs are synonymous and ~4,467 SNPs are non-synonymous (N/S ~0.28). The remaining SNPs are in other categories, such as stop codon gains (38 SNPs), stop codon losses (8 SNPs), and start codon gains (297 SNPs) in the 5′UTR. We found, as expected, that the SNP frequency is proportional to the recombination frequency (i.e., highest in the middle of chromosome arms). We also found that start-gain or stop-lost SNPs in Drosophila melanogaster often result in additions of N-terminal or C-terminal amino acids that are conserved in other Drosophila species. It appears that the 5′ and 3′ UTRs are reservoirs for genetic variations that changes the termini of proteins during evolution of the Drosophila genus. As genome sequencing is becoming inexpensive and routine, SnpEff enables rapid analyses of whole-genome sequencing data to be performed by an individual laboratory.

Using Drosophila melanogaster as a Model for Genotoxic Chemical Mutational Studies with a New Program, SnpSift

This paper describes a new program SnpSift for filtering differential DNA sequence variants between two or more experimental genomes after genotoxic chemical exposure. Here, we illustrate how SnpSift can be used to identify candidate phenotype-relevant variants including single nucleotide polymorphisms, multiple nucleotide polymorphisms, insertions, and deletions (InDels) in mutant strains isolated from genome-wide chemical mutagenesis of Drosophila melanogaster. First, the genomes of two independently isolated mutant fly strains that are allelic for a novel recessive male-sterile locus generated by genotoxic chemical exposure were sequenced using the Illumina next-generation DNA sequencer to obtain 20- to 29-fold coverage of the euchromatic sequences. The sequencing reads were processed and variants were called using standard bioinformatic tools. Next, SnpEff was used to annotate all sequence variants and their potential mutational effects on associated genes. Then, SnpSift was used to filter and select differential variants that potentially disrupt a common gene in the two allelic mutant strains. The potential causative DNA lesions were partially validated by capillary sequencing of polymerase chain reaction-amplified DNA in the genetic interval as defined by meiotic mapping and deletions that remove defined regions of the chromosome. Of the five candidate genes located in the genetic interval, the Pka-like gene CG12069 was found to carry a separate pre-mature stop codon mutation in each of the two allelic mutants whereas the other four candidate genes within the interval have wild-type sequences. The Pka-like gene is therefore a strong candidate gene for the male-sterile locus. These results demonstrate that combining SnpEff and SnpSift can expedite the identification of candidate phenotype-causative mutations in chemically mutagenized Drosophila strains. This technique can also be used to characterize the variety of mutations generated by genotoxic chemicals.

Genetic and vascular modifiers of age-sensitive cognitive skills: effects of COMT, BDNF, ApoE, and hypertension.

Several single nucleotide polymorphisms have been linked to neural and cognitive variation in healthy adults. We examined contribution of three polymorphisms frequently associated with individual differences in cognition (Catechol-O-Methyl-Transferase Val158Met, Brain-Derived-Neurotrophic-Factor Val66Met, and Apolipoprotein E epsilon4) and a vascular risk factor (hypertension) in a sample of 189 volunteers (age 18-82). Genotypes were determined from buccal culture samples, and cognitive performance was assessed in 4 age-sensitive domains?fluid intelligence, executive function (inhibition), associative memory, and processing speed. We found that younger age and COMT Met/Met genotype, associated with low COMT activity and higher prefrontal dopamine content, were independently linked to better performance in most of the tested domains. Homozygotes for Val allele of BDNF polymorphism exhibited better associative memory and faster speed of processing than the Met allele carriers, with greater effect for women and persons with hypertension. Carriers of ApoE epsilon4 allele evidenced steeper age-related increase in costs of Stroop color interference, but showed no negative effects on memory. The findings indicate that age-related cognitive performance is differentially affected by distinct genetic factors and their interactions with vascular health status.

Epigenetic silencing of multiple interferon pathway genes after cellular immortalization

Abrogating cellular senescence is a necessary step in the formation of a cancer cell. Promoter hypermethylation is an epigenetic mechanism of gene regulation known to silence gene expression in carcinogenesis. Treatment of spontaneously immortal Li–Fraumeni fibroblasts with 5-aza-2′-deoxycytidine (5AZA-dC), an inhibitor of DNA methyltransferase (DNMT), induces a senescence-like state. We used microarrays containing 12 558 genes to determine the gene expression profile associated with cellular immortalization and also regulated by 5AZA-dC. Remarkably, among 85 genes with methylation-dependent downregulation (silencing) after immortalization, 39 (46%) are known to be regulated during interferon signaling, a known growth-suppressive pathway. This work indicates that gene silencing may be associated with an early event in carcinogenesis, cellular immortalization.

Gene expression profiles predict survival and progression of pleural mesothelioma

Purpose: Clinical outcomes for malignant pleural mesothelioma (MPM) patients having surgery are imprecisely predicted by histopathology and intraoperative staging. We hypothesized that gene expression profiles could predict time to progression and survival in surgically cytoreduced pleural mesothelioma of all stages. Experimental Design: Gene expression analyses from 21 MPM patients having cytoreductions and identical postoperative adjuvant therapy were performed using the U95 Affymetrix gene chip. Using both dChip and SAM, neural networks constructed a common 27 gene classifier, which was associated with either the high-risk and low-risk group of patients. Data were validated using real-time PCR and immunohistochemical staining. The 27 gene classifier was also used for validation in a separate set of 17 MPM patients from another institution. Results: The groups predicted by the gene classifier recapitulated the actual time to progression and survival of the test set with 95.2% accuracy using 10-fold cross-validation. Clinical outcomes were independent of histology, and heterogeneity of progression and survival in early stage patients was defined by the classifier. The gene classifier had a 76% accuracy in the separate validation set of MPMs. Conclusions: These data suggest that pretherapy gene expression analysis of mesothelioma biopsies may predict which patients may benefit from a surgical approach.

Tobacco and estrogen metabolic polymorphisms and risk of non-small cell lung cancer in women

To explore the potential role for estrogen in lung cancer susceptibility, candidate single-nucleotide polymorphism (SNPs) in tobacco and estrogen metabolism genes were evaluated. Population-based cases (n = 504) included women aged 18-74, diagnosed with NSCLC in metropolitan Detroit between November 2001 and October 2005. Population-based controls (n = 527) were identified through random digit dialing and matched on race and age. Eleven SNPs in 10 different genes were examined in relation to risk: CYP1A1 Msp1, CYP1A1 Ile462Val, CYP1B1 Leu432Val, CYP17, CYP19A1, XRCC1 Gln399Arg, COMT Val158Met, NQO1 Pro187Ser, GSTM1, GSTT1 and GSTP1 Ile105Val. Lung cancer risk associated with individual SNPs was seen for GSTP1 [A allele; odds ratio (OR) = 1.85; 95% confidence interval (CI), 1.04-3.27] and XRCC1 (A/A genotype; OR = 1.68; 95% CI, 1.01-2.79) in white women and CYP1B1 (G allele; OR = 11.1; 95% CI, 1.18-104) in black women smokers. White women smokers carrying two risk genotypes at the following loci were at increased risk of lung cancer compared with individuals not carrying risk alleles at these loci: CYP17 and GSTM1, COMT and GSTM1, CYP17 and GSTT1, XRCC1 and GSTP1, CYP1B1 and XRCC1 and COMT and XRCC1. The most parsimonious model of lung cancer risk in white smoking women included age, family history of lung cancer, history of chronic lung disease, pack-years, body mass index, XRCC1 A/A genotype, GSTM1 null and COMT A/G or G/G genotype. These findings support the need for continued study of estrogen in relation to lung cancer risk. Polymorphisms in the tobacco metabolism, estrogen metabolism and DNA repair pathways will be useful in developing more predictive models of individual risk.

Volume of white matter hyperintensities in healthy adults: Contribution of age, vascular risk factors, and inflammation-related genetic variants

Aging is associated with appearance of white matter hyperintensities (WMH) on MRI scans. Vascular risk and inflammation, which increase with age, may contribute to white matter deterioration and proliferation of WMH. We investigated whether circulating biomarkers and genetic variants associated with elevated vascular risk and inflammation are associated with WMH volume in healthy adults (144 volunteers, 44-77 years of age). We examined association of WMH volume with age, sex, hypertension, circulating levels of total plasma homocysteine (tHcy), cholesterol (low-density lipoprotein), and C-reactive protein (CRP), and four polymorphisms related to vascular risk and inflammation: Apolipoprotein ε (ApoE ε2,3,4), Angiotensin-Converting Enzyme insertion/deletion (ACE I/D), methylenetetrahydrofolate reductase (MTHFR) C677T, C-reactive protein (CRP)-286C>A>T, and interleukin-1β (IL-1β) C-511T. We found that larger WMH volume was associated with advanced age, hypertension, and elevated levels of homocysteine and CRP but not with low-density lipoprotein levels. Homozygotes for IL-1β-511T allele and carriers of CRP-286T allele that are associated with increased inflammatory response had larger WMH than the other allelic combinations. Carriers of the APOE ε2 allele had larger frontal WMH than ε3 homozygotes and ε4 carriers did. Thus, in healthy adults, who are free of neurological and vascular disease, genetic variants that promote inflammation and elevated levels of vascular risk biomarkers can contribute to brain abnormalities. This article is part of a Special Issue entitled: Imaging Brain Aging and Neurodegenerative disease.

Genes Associated with Prostate Cancer Are Differentially Expressed in African American and European American Men

Background: Despite more aggressive screening across all demographics and gradual declines in mortality related to prostate cancer (PCa) in the United States, disparities among populations persist. A substantial proportion of African American men (AAM) have a higher overall incidence, earlier age of onset, increased proportion of clinically advanced disease, and increased bone metastases and mortality from PCa compared to European American men (EAM). Limited early evidence indicates that underlying causes for disparities may be observed in tumor-specific gene expression programs. Methods: This study used microarray-based methods to measure expression levels for 517 genes that were previously associated with PCa in archived formalin-fixed paraffin embedded (FFPE) specimens; testing the hypothesis that gene expression features of functional consequence to cancer distinguish PCa from AAM and EAM. A t test was conducted comparing AAM to EAM expression levels for each probe on the array. Results: Analysis of 639 tumor samples (270 AAM, 369 EAM) showed that 95 genes were overexpressed specifically in PCa from AAM relative to EAM and 132 were overexpressed in PCa from EAM relative to AAM. Furthermore, systems-level analyses highlight the relevant signaling pathways and functions associated with the EAM- or AAM-specific overexpressed gene sets, for example, inflammation and lipid metabolism. Conclusions: Results here bring further understanding to the potential for molecular differences for PCa in AAM versus EAM. Impact: The results support the notion that therapeutic benefits will be realized when targeted treatments are designed to acknowledge and address a greater spectrum of PCa subtypes and molecular distinctions. Cancer Epidemiol Biomarkers Prev; 22(5); 891–7. ©2013 AACR.

Racial Differences in the Association Between SNPs on 15q25.1, Smoking Behavior, and Risk of Non-small Cell Lung Cancer

Introduction: Three genome-wide association studies identified a region on chromosome 15q25.1 associated with lung cancer and measures of nicotine addiction. This region includes nicotinic acetylcholine receptor subunit genes CHRNA3 and CHRNA5. These studies were conducted in European or European American populations and do not provide risk estimates for African Americans. The goal of this study was to determine whether recently identified genetic variation in 3 SNPs (rs1051730, rs931794, rs8034191) on chromosome 15q25.1 contributes to risk of lung cancer in African Americans. Methods: Data were derived from three case-control studies. Participants included 1058 population-based non-small cell lung cancer cases selected from the Detroit area SEER registry and 1314 controls matched within study by age, race, and sex. Thirty-nine percent of participants were African American. Results: Risk associated with rs1051730 (odds ratio 1.59; 95% confidence interval 1.16-2.19) and rs931794 (odds ratio 1.39; 95% confidence interval 1.09-1.78) increased in ever smoking African Americans adjusting for cigarettes smoked per day. Among white cases, the number of cigarettes smoked varied by genotype at all three SNPs, and when smoking quantity was included in the models, risk was not significantly associated with any of the three SNPs. Conclusions: These findings suggest that SNPs in the CHRNA3 and CHRNA5 region contribute to lung cancer risk, and while variant alleles are less frequent in African Americans, risk in this group may be greater than in whites and less likely to reflect an indirect effect on lung cancer risk through nicotine dependence.

Cytokine SNPs: Comparison of allele frequencies by race and implications for future studies.

The role of inflammation is being considered in chronic diseases. Previous studies have examined SNPs in a few key inflammatory genes and have included small numbers of African American participants. Variation in the frequencies of inflammatory pathway SNPs may help to explain racial disparities in disease risk. Through a population-based study of 103 African American and 380 Caucasian unrelated, healthy women, we examined the relationships between race and allele frequencies of 70 cytokine and cytokine receptor SNPs. The associations between genotypic and haplotype frequencies and race were also analyzed. Allelic frequencies for 52 out of the 70 SNPs meeting criteria for analysis differed significantly by race. Of the 32 pro-inflammatory and 20 anti-inflammatory SNPs for which the allele frequencies varied significantly by race, variant allele frequency differences between Caucasians and African Americans ranged between 6-37% and 7-53% for pro-inflammatory SNPs and anti-inflammatory SNPs, respectively. Our findings suggest that while allele frequencies do vary by race, racial groups are not simplistically represented by a pro-inflammatory or anti-inflammatory genetic profile. Given the racial variability in allele frequencies in inflammatory gene SNPs, studies examining the association between these SNPs and disease should at least incorporate self-reported race in their analyses.