Candidate Functional Single Nucleotide Polymorphisms Research Articles

219 Identification of Functional Single Nucleotide Polymorphisms Responsible for Alternative Splicing Affecting Gastrointestinal Nematodes Resistance in Grazing Sheep

Abstract Sheep infected with gastrointestinal nematodes (GIN) can show weight loss, anemia, diarrhea, decrease appetite, and hypoproteinemia, which can sometimes lead to death. These signs of infection decrease productivity and cause economic losses. Economic losses can also be associated with treatment without efficacy due to GIN resistance to most anthelmintic classes. Combined approaches need to be used to reduce the reliance on anthelmintic as the only control method to GIN. The aim of this study is the identification of potential functional single nucleotide polymorphisms (SNPs) responsible for alternate splicing in Rideau X Dorset crossbred sheep with different immune profiles using RNA-Sequencing (RNA-Seq). Animals were exposed to GIN while grazing, after that those with the greatest and least gastrointestinal tract parasite burdens were selected for liver RNA extracting. These animals can be separated into groups of high (H; n = 5) and medium (M; n = 6) innate immune responses. In addition, liver RNA was extracted for GIN-unexposed control lambs (U; n = 4). RNA-Seq analysis was performed using the CLC genomics workbench software, version 20.0.4 (CLC Bio, Aarhus, Denmark). Quality control analyses, including guanine-cytosine (GC) content, ambiguous base content, Phred score, base coverage, nucleotide contributions, and over-represented sequences parameters, were performed on fastq files. Sequence reads were aligned to the annotated Oar_rambouillet_v.1 (release 109) ovine reference genome. After filtering steps, mRNA isoforms with multiple (³ 2) expressed mRNA isoforms with at least one mRNA isoform differently expressed (DE) between H- and M-innate immune response and GIN-unexposed control groups (FDR < 0.05 and |FC| > 2) were keep. A total of 1, 75, and 58 DE mRNA isoforms were identified for M versus H, U versus M, and U versus H comparisons, respectively. Gene ontology (GO) analysis including the biological process, molecular function, and cellular component GO categories and metabolic pathways analysis were performed using STRING-db software. Enriched biological processes and metabolic pathways related to the immune system, as well as cholesterol and other sterol metabolism (e.g., sterol biosynthesis, steroid metabolism, cholesterol biosynthesis, cholesterol metabolism, ascorbate and aldarate metabolism, drug metabolism, and terpenoid biosynthesis) were identified (FDR < 0.05). Next, the list of DE mRNA isoforms identified will be deeply studied by performing SNP discovery analysis to validate if specific functional SNPs are responsible for the alternative splicing affecting GIN resistance in grazing sheep. The complete results will help to better understand the biological mechanisms that may be regulating the immune response to GIN. Candidate functional SNPs have the potential to be included in a commercial SNP panel for selection of more resistant animals to improve GIN resistance in sheep flocks.

Identification of a non‐coding SNP associated with risk for non‐syndromic orofacial clefting with allele‐specific effects on IRF6 expression in vitro

Introduction Genome wide association studies (GWAS) have identified loci associated with risk for non-syndromic cleft lip with or without cleft palate (NSCL/P). At every locus there are multiple single nucleotide polymorphisms (SNPs) associated with disease and it is a challenge to distinguish functional SNPs from those merely in linkage disequilibrium with functional SNPs. Objective At many of the associated loci, the presumed risk-relevant gene is expressed in oral epithelium. We hypothesize that at such loci functional SNPs have allele-specific effects on the activity of oral epithelium enhancers. Our objective is to identify the candidate functional SNPs for NSCL/P at each of the loci identified by GWAS. Materials and Methods To test SNPs for their effects on enhancer activity we carry out massively parallel reporter assays (MPRAs) on 890 SNPs from 8 loci in the GMSM-K human fetal oral epithelium cell line. To test the SNP-target gene association we use the Activity-by-Contact method, which incorporates open chromatin, H3K27Ac, and HiC data. To confirm MPRA results on top-scoring SNPs we perform luciferase reporter assays. We use CRISPR-mediated homology directed repair to engineer GMSM-K cells to be homozygous for risk or non-risk alleles for promising SNPs and assess the expression level of the target gene in engineered cells of each genotype. Assays that pass these assays are tested in allele-specific reporter assays in zebrafish and mouse embryos. Results Using the described methods we tested 608 NSCL/P-associated SNPs in the IRF6 locus and 11 such SNPs in the FOXE1 locus. Two SNPs in the first group and one in the second had consistent allele-specific effects in the MPRA and in luciferase assays. Importantly, for one SNP in the IRF6 locus, we have engineered cells to be homozygous for the non-risk-associated or risk-associated allele and found expression of IRF6 is lower in the latter. In vivo enhancer reporter assays in mouse and zebrafish embryos are ongoing. Conclusion Out of more than 600 NSCL/P-risk-associated SNPs in the IRF6 locus, we have evidence that a particular one directly contributes to pathogenesis for this disorder. Significance This study illustrates a method to screen among multiple non-coding SNPs identified by GWAS to identify candidates for those that are functional. Moreover, identifying the functional SNP near IRF6 is a step towards illuminating how a common variant contributes risk for this disorder.

Comprehensive Integrative Analyses Identify TIGD5 rs75547282 as a Risk Variant for Autism Spectrum Disorder.

Although recent genome-wide association studies have identified risk loci that strongly associates with autism spectrum disorder (ASD), how to pinpoint the causal genes remains a challenge. We aimed to pinpoint the potential causal genes and explore the possible susceptibility and mechanism. A convergent functional genomics (CFG) method was used to prioritize the candidate genes by combining lines of evidence, including Sherlock analysis, spatio-temporal expression patterns, expression analysis, protein-protein interactions, co-expression and association with brain structure. A higher score in the CFG approach suggested that more evidence supported this gene as an ASD risk gene. We screened genes with higher CFG scores for candidate functional single nucleotide polymorphisms (SNPs). A genotyping experiment (602 ASD children and 604 healthy sex-matched children) and the dual-luciferase reporter gene assay were followed to validate the effects of SNPs. We identified three genes (MAPT, ZNF285, and TIGD5) as candidate causal genes using the CFG approach. The genotyping experiment showed that TIGD5 rs75547282 was associated with an increased risk of ASD under the dominant model (OR = 1.37, 95% CI = 1.09-1.72, P = 0.006) though the statistical power was limited (5.2%). The T allele of rs75547282 activated the expression of TIGD5 compared with the C allele in the dual-luciferase reporter assay. Our study indicates that such comprehensive integrative analyses may be an effective way to explore promising ASD susceptibility variants and needs to be further investigated in future research. Genotyping experiments should, however, be based on a larger population sample to increase statistical power. LAY SUMMARY: We set out to pinpoint the potential causal genes of ASD and explore the possible susceptibility and mechanism by combining lines of evidence from different analyses. Our results show that TIGD5 rs75547282 is associated with the risk of ASD in the Han Chinese population. In addition, a similar framework to seek promising ASD risk variants could be further investigated in future research Autism Res 2021, 14: 631-644. © 2021 International Society for Autism Research and Wiley Periodicals LLC.

A functional variant on 20q13.33 related to glioma risk alters enhancer activity and modulates expression of multiple genes.

Genome‐wide association studies (GWAS) have identified single‐nucleotide polymorphisms (SNPs) associated with glioma risk on 20q13.33, but the biological mechanisms underlying this association are unknown. We tested the hypothesis that a functional SNP on 20q13.33 impacted the activity of an enhancer, leading to an altered expression of nearby genes. To identify candidate functional SNPs, we identified all SNPs in linkage disequilibrium with the risk‐associated SNP rs2297440 that mapped to putative enhancers. Putative enhancers containing candidate functional SNPs were tested for allele‐specific effects in luciferase enhancer activity assays against glioblastoma multiforme (GBM) cell lines. An enhancer containing SNP rs3761124 exhibited allele‐specific effects on activity. Deletion of this enhancer by CRISPR‐Cas9 editing in GBM cell lines correlated with an altered expression of multiple genes, including STMN3, RTEL1, RTEL1‐TNFRSF6B, GMEB2, and SRMS. Expression quantitative trait loci (eQTL) analyses using nondiseased brain samples, isocitrate dehydrogenase 1 (IDH1) wild‐type glioma, and neurodevelopmental tissues showed STMN3 to be a consistent significant eQTL with rs3761124. RTEL1 and GMEB2 were also significant eQTLs in the context of early CNS development and/or in IDH1 wild‐type glioma. We provide evidence that rs3761124 is a functional variant on 20q13.33 related to glioma/GBM risk that modulates the expression of STMN3 and potentially other genes across diverse cellular contexts.

Functional Genetic Polymorphisms in the IL1RL1–IL18R1 Region Confer Risk for Ocular Behçet’s Disease in a Chinese Han Population

Single nucleotide polymorphisms (SNPs) in the IL1RL1–IL18R1 region are associated with various immune-mediated diseases. This study was carried out to investigate the causal variant for ocular Behçet’s disease (BD) and elucidate its target genes in the IL1RL1–IL18R1 region. Nine candidate functional SNPs were prioritized with bioinformatics analysis, followed by a two-stage association study in 694 ocular BD patients and 1,458 unaffected controls. Functional studies were performed in the peripheral blood mononuclear cells (PBMCs) of 45 healthy men and 16 active male BD patients. Genotyping was performed using the MassARRAY System. The mRNA expressions of IL1RL1, IL18R1, IL18RAP, and SLC9A4 were assayed by real-time PCR and secretion of cytokines was examined by ELISA. Significantly lower frequencies of the rs12987977 GG genotype/G allele (Pc = 8.93 × 10–7, OR = 0.39; Pc = 2.60 × 10–3, OR = 0.77, respectively), rs12999364 TT genotype/T allele (Pc = 3.15 × 10–4, OR = 0.51; Pc = 1.13 × 10–2, OR = 0.80, respectively), and rs4851569 AA genotype/A allele (Pc = 3.29 × 10–4, OR = 0.52; Pc = 9.72 × 10–3, OR = 0.80, respectively) were observed in BD patients compared with the controls. Functional experiments revealed a downregulation of IL1RL1, IL18R1, and SLC9A4 and a decreased secretion of IFN-γ in the anti-CD3/CD28 antibody-treated PBMCs as well as a decreased production of TNF-α in the lipopolysaccharide (LPS)-stimulated PBMCs in carriers of the protective homozygous rs12987977/GG genotype compared with the TT genotype. Our findings show that functional SNPs—rs12987977, rs12999364, and rs4851569—in the IL1RL1–IL18R1 region confer susceptibility to ocular BD in a Chinese Han population. And IL1RL1, IL18R1, and SLC9A4 may be the target genes of rs12987977.

Analysis of 47 Non-MHC Ankylosing Spondylitis Susceptibility Loci Regarding Associated Variants across Whites and Han Chinese.

To present a systematic evaluation of 47 non-MHC ankylosing spondylitis (AS) susceptibility loci that have been initially discovered through white genome-wide association studies in Han Chinese. Originally, 10,743 samples representing north and south Chinese in 4 datasets were obtained. After data quality control and imputation, metaanalysis results of 94,621 variants within 47 loci were extracted. Four ERAP1 single-nucleotide polymorphisms (SNP) and HLA-B27 tag SNP rs13202464 were used for interaction analysis. Population-attributable risk percentages of AS-associated variants were compared. Functional annotations of AS-associated variants were conducted using HaploReg, RegulomeDB, and rVarBase databases. We revealed 16 AS-associated variants with nominal evidence in Han Chinese, including rs10865331 (p = 6.30 × 10-10), rs10050860 (p = 4.09 × 10-5) and rs8070463 (p = 1.03 × 10-4). Potential susceptible SNP within these 47 loci were also identified, such as rs13024541 (2p15), rs17401719 (5q15), and rs62074054 (17q21). Epistatic interactions between 3 ERAP1 SNP (rs17401719, rs30187, and rs10050860) and HLA-B27 were confirmed. Among the 16 AS-associated variants, rs30187 showed weaker risk effect, while rs10050860 and rs12504282 seemed to attribute more risk in Han Chinese than in whites. Further genomic annotation pinpointed 35 candidate functional SNP, especially in the 2p15, ERAP1, and NPEPPS-TBKBP1 regions. Our results provided a detailed spectrum of all the reported non-MHC AS susceptibility loci in Han Chinese, which comprehensively exhibited the ethnic heterogeneity of AS susceptibility and highlighted that 2p15, ERAP1, and NPEPPS-TBKBP1 regions may play a critical role in AS pathogenesis across diverse populations.

Genetic variants in SLC22A3 contribute to the susceptibility to colorectal cancer.

Previous a genome‐wide association study (GWAS) of colorectal cancer in Japanese population has identified a risk region at the chromosome 6q26‐q27 associated with colorectal cancer risk. However, the causal gene at this locus remained unclear. In our study, we enrolled a total of 14 candidate functional single nucleotide polymorphisms (SNPs) at 6q26‐q27 (318 kb), and then genotyped them by TaqMan method in a Chinese population including 1,147 colorectal cancer cases and 1,203 controls. Among that, 5 SNPs were identified statistical association with colorectal cancer risk by logistic regression analysis. Of which, SNP rs420038 G > A in SLC22A3 was related to decreased risk of colorectal cancer (adjusted odds ratio (OR) = 0.79, 95% confidence interval (CI) = 0.67–0.94, p = 0.007), and also associated with lower expression of SLC22A3 (p = 0.040) using expression quantitative trait loci (eQTL) analysis. Moreover, by the luciferase assays, we found that compared to the G allele of rs420038, the A allele could suppress the activity of the promoter in SLC22A3. Furthermore, the expression of SLC22A3 was significantly higher in colorectal cancer tissues than that in paired normal tissues (p < 0.001). Meanwhile, the phenotypes of proliferation, migration, invasion, cell cycle and apoptosis of colorectal cancer cell were significantly affected by SLC22A3 in vitro. Our results revealed a novel susceptible locus, rs420038 in SLC22A3, which may be involved in colorectal cancer development and progression.

Whole-genome sequencing of eight goat populations for the detection of selection signatures underlying production and adaptive traits

The goat (Capra hircus) is one of the first farm animals that have undergone domestication and extensive natural and artificial selection by adapting to various environments, which in turn has resulted in its high level of phenotypic diversity. Here, we generated medium-coverage (9–13×) sequences from eight domesticated goat breeds, representing morphologically or geographically specific populations, to identify genomic regions representing selection signatures. We discovered ~10 million single nucleotide polymorphisms (SNPs) for each breed. By combining two approaches, ZHp and di values, we identified 22 genomic regions that may have contributed to the phenotypes in coat color patterns, body size, cashmere traits, as well as high altitude adaptation in goat populations. Candidate genes underlying strong selection signatures including coloration (ASIP, KITLG, HTT, GNA11, and OSTM1), body size (TBX15, DGCR8, CDC25A, and RDH16), cashmere traits (LHX2, FGF9, and WNT2), and hypoxia adaptation (CDK2, SOCS2, NOXA1, and ENPEP) were identified. We also identified candidate functional SNPs within selected genes that may be important for each trait. Our results demonstrated the potential of using sequence data in identifying genomic regions that are responsible for agriculturally significant phenotypes in goats, which in turn can be used in the selection of goat breeds for environmental adaptation and domestication.

Fine-mapping of two differentiated thyroid carcinoma susceptibility loci at 9q22.33 and 14q13.3 detects novel candidate functional SNPs in Europeans from metropolitan France and Melanesians from New Caledonia.

Incidence of differentiated thyroid carcinoma varies considerably between countries and ethnic groups, with particularly high incidence rates in Melanesians of New Caledonia. Differentiated thyroid cancer (DTC) has a familial relative risk higher than other cancers, highlighting the contribution of inherited factors to the disease. Recently, genome-wide association studies (GWAS) identified several DTC susceptibility loci. The most robust associations were reported at loci 9q22 (rs965513 and rs1867277) and 14q13 (rs944289 and rs116909734). In this study, we performed a fine-mapping study of the two gene regions among Europeans and Melanesians from Metropolitan France and New Caledonia. We examined 81 single nucleotide polymorphisms (SNPs) at 9q22 and 561 SNPs at 14q13 in Europeans (625 cases/776 controls) and in Melanesians (244 cases/189 controls). The association with the four SNPs previously identified in GWAS was replicated in Europeans while only rs944289 was replicated in Melanesians. Among Europeans, we found that the two SNPs previously reported at 9q22 were not independently associated to DTC and that rs965513 was the predominant signal; at 14q13, we showed that the haplotype rs944289[C]-rs116909374[C]-rs999460[T] was significantly associated with DTC risk and that the association with rs116909374 differed by smoking status (p-interaction = 0.03). Among Melanesians, a new independent signal was observed at 14q13 for rs1755774 which is strongly correlated to rs2787423; this latter is potentially a functional variant. Significant interactions with parity (p < 0.05) and body mass index were observed for rs1755774 and rs2787423. This study contributed to a better characterization of the DTC loci 9q22 and 14q13 in Europeans and in Melanesians and has identified novel variants to be prioritized for further functional studies.

CAUSEL: an epigenome- and genome-editing pipeline for establishing function of noncoding GWAS variants.

The vast majority of disease-associated single nucleotide polymorphisms (SNPs) mapped by genome-wide association studies (GWAS) are located in the non-protein coding genome, but establishing the functional and mechanistic roles of these sequence variants has proven challenging. Here, we describe a general pipeline in which candidate functional SNPs are first evaluated by fine-mapping, epigenomic profiling, and epigenome editing and then interrogated for causal function by using genome editing to create isogenic cell lines. To validate this approach, we analyzed the 6q22.1 prostate cancer risk locus and identified rs339331 as the top scoring SNP. Epigenome editing confirmed that rs339331 possessed regulatory potential. Using transcription activator-like effector nuclease (TALEN)-mediated genome-editing, we created a panel of isogenic 22Rv1 prostate cancer cell lines representing all three genotypes (TT, TC, CC) at rs339331. Introduction of the “T” risk allele increased transcription of the RFX6 gene, increased HOXB13 binding at the rs339331 region, and increased deposition of the enhancer-associated H3K4me2 histone mark at the rs339331 region. The cell lines also differed in cellular morphology and adhesion, and pathway analysis of differentially expressed genes suggested an influence of androgens. In summary, we have developed and validated a widely accessible approach to establish functional causality for non-coding sequence variants identified by GWAS.

Brain tumor risk according to germ-line variation in the MLLT10 locus.

Genome-wide association studies have recently identified a cancer susceptibility locus at 10p12 mapping to MLLT10 associated with the onset of diverse tumors. We genotyped two tightly linked single-nucleotide polymorphisms (SNPs) at MLLT10 associated with meningioma (rs12770228) or ovarian cancer (rs1243180), and tested for associations among 295 meningioma cases, 606 glioma cases and 646 noncancer controls, all of European descent. The variant 'A' allele in MLLT10 rs12770228 was associated with an increased risk of meningioma (per allele odds ratio: 1.25; 95% confidence interval: 1.02, 1.53; P=0.031). Similar associations were observed for rs1243180. MLLT10 variants were unrelated to glioma. Functional investigation identified 22 candidate functional SNPs mapping to this region. The present study further validates 10p12 as a meningioma risk locus.

Abstract IA19: The GAME-ON (Genetic Associations and Mechanisms in Oncology) Consortium: Assessing the functional relevance of genetic variation.

Abstract Recently, significant advances in genetics, in particular in genome-wide association studies (GWAS), have uncovered a large number of chromosomal loci associated with cancer risk. Despite these advances the determination of the molecular mechanism behind the detected associations has proven to be a significant challenge. The design of most GWAS genotyping chips take advantage of well characterized tagging single nucleotide polymorphisms (SNPs). These “tags” represent specific proxy markers for a series of other SNPs in high linkage disequilibrium with the tagging SNP in a defined chromosomal region. While the use of tagging SNPs result in reproducible and robust signals during genotyping and minimize the number of SNPs to be tested, the expectation is that they may not necessarily constitute the causal SNP, the nucleotide change that result in a relevant biological activity responsible for cancer predisposition. Identification of predisposition loci in prostate, ovary, and breast have been accelerated by the COGS (Collaborative Oncologic Gene-Environment Study) project which typed 200,000 SNPs and genotyped an unprecedented number of cancer cases and controls (˜200,000). The GAME-ON (Genetic Associations and Mechanisms in Oncology) consortium has poured over GWAS data for prostate, ovary, breast, lung, and colorectal cancers to develop systematic procedures to dissect functional contribution of SNPs and their target genes. Here I describe our experience in conducting functional analyses for these loci, focusing on ovarian cancer susceptibility loci. The overall approach can be loosely divided in three stages: in the first we use the most significant SNP found in GWAS analyses or meta-analyses as a starting point to interrogate other candidate functional SNPs in the region. The second stage explores potential target genes in the region. Finally, we focus on tying the candidate functional SNP to the candidate target gene. This final stage may also focus on, depending on the availability of in vivo models, attempts to uncover direct evidence for the participation of the SNP/target gene pairs in the mechanism of oncogenesis. Although in some cases the analytic approaches and tools might be specific to the cancer in question, we feel that the main framework can be utilized for different cancers. Most importantly, perhaps, are the lessons learned during this process as well as identifying the limitations and challenges that lay ahead. Citation Format: Alvaro N.A. Monteiro. The GAME-ON (Genetic Associations and Mechanisms in Oncology) Consortium: Assessing the functional relevance of genetic variation. [abstract]. In: Proceedings of the AACR Special Conference on Post-GWAS Horizons in Molecular Epidemiology: Digging Deeper into the Environment; 2012 Nov 11-14; Hollywood, FL. Philadelphia (PA): AACR; Cancer Epidemiol Biomarkers Prev 2012;21(11 Suppl):Abstract nr IA19.

FunciSNP: an R/bioconductor tool integrating functional non-coding data sets with genetic association studies to identify candidate regulatory SNPs.

Single nucleotide polymorphisms (SNPs) are increasingly used to tag genetic loci associated with phenotypes such as risk of complex diseases. Technically, this is done genome-wide without prior restriction or knowledge of biological feasibility in scans referred to as genome-wide association studies (GWAS). Depending on the linkage disequilibrium (LD) structure at a particular locus, such tagSNPs may be surrogates for many thousands of other SNPs, and it is difficult to distinguish those that may play a functional role in the phenotype from those simply genetically linked. Because a large proportion of tagSNPs have been identified within non-coding regions of the genome, distinguishing functional from non-functional SNPs has been an even greater challenge. A strategy was recently proposed that prioritizes surrogate SNPs based on non-coding chromatin and epigenomic mapping techniques that have become feasible with the advent of massively parallel sequencing. Here, we introduce an R/Bioconductor software package that enables the identification of candidate functional SNPs by integrating information from tagSNP locations, lists of linked SNPs from the 1000 genomes project and locations of chromatin features which may have functional significance. Availability: FunciSNP is available from Bioconductor (bioconductor.org).

Application of statistical and functional methodologies for the investigation of genetic determinants of coronary heart disease biomarkers: lipoprotein lipase genotype and plasma triglycerides as an exemplar

Genome-wide association studies have proved very successful in identifying novel single-nucleotide polymorphisms (SNPs) associated with disease or traits, but the related, functional SNP is usually unknown. In this paper, we describe a methodology to locate and validate candidate functional SNPs using lipoprotein lipase (LPL), a gene previously associated with triglyceride levels, as an exemplar. Two thousand seven hundred and eighty-six healthy middle-aged men from the NPHSII UK prospective study (with up to six measures of plasma lipid levels) were genotyped for 20 LPL tagging (t)SNPs using Illumina Bead technology. Using model-selection procedures and haplotypes, we identified eight SNPs that consistently maximized the fit of the model to the phenotype. Fifteen SNPs in high linkage disequilibrium with these were identified, and functional assays were carried out on all 23 SNPs. Electrophoretic mobility shift assay (EMSA) was used to identify SNPs that had the potential to alter DNA–protein interactions, reducing the number to eight possible candidate SNPs. These were examined for ability to alter expression using a luciferase reporter assay, and two regulatory SNPs, showing genotype differences, rs327 and rs3289, were identified. Finally, multiplexed-competitor-EMSA (MC-EMSA) and supershift EMSA identified FOXA2 to rs327T, and CREB-binding protein (CBP) and CCAAT displacement protein (CDP) to rs3289C as the factors responsible for transcription binding. We have identified two novel candidate functional SNPs in LPL and presented a procedure aimed to efficiently detect SNPs potentially causal to genetic association. We believe that this methodology could be successfully applied to future re-sequencing data.

Variation in genes required for normal mitosis and risk of breast cancer

The down-regulation of genes involved in normal cell division can cause aberrant mitoses and increased cell death. Surviving cells exhibit aneuploidy and/or polyploidy. Since mitotic disruption has been linked with tumor development and progression, alterations in the expression or activity of these mitotic regulators may contribute to breast tumor formation. We evaluated associations between common inherited variation in these genes and breast cancer risk. Two hundred and five tagging and candidate functional single nucleotide polymorphisms in 30 genes required for normal cell division were genotyped in 798 breast cancer cases and 843 controls from the Mayo Clinic breast cancer study. Two variants in EIF3A (rs10787899 and rs3824830; P < 0.01) and four variants in SART1 (rs660118, rs679581, rs754532, and rs735942; P(trend) < or = 0.02) were significantly associated with an altered risk of breast cancer along with single variants in RRM2, PSCD3, C11orf51, CDC16, SNW1, MFAP1, and CDC2 (P < 0.05). Variation in both SART1 (P = 0.009) and EIF3A (P = 0.02) was also significant at the gene level. Analyses suggested that SART1 SNPs rs660118 and rs679581 accounted for the majority of the association of that gene with breast cancer. The observed associations between breast cancer risk and genetic variation in the SART1 and EIF3A genes that are required for maintenance of normal mitosis suggest a direct role for these genes in the development of breast cancer.

LS-SNP: large-scale annotation of coding non-synonymous SNPs based on multiple information sources

The NCBI dbSNP database lists over 9 million single nucleotide polymorphisms (SNPs) in the human genome, but currently contains limited annotation information. SNPs that result in amino acid residue changes (nsSNPs) are of critical importance in variation between individuals, including disease and drug sensitivity. We have developed LS-SNP, a genomic scale software pipeline to annotate nsSNPs. LS-SNP comprehensively maps nsSNPs onto protein sequences, functional pathways and comparative protein structure models, and predicts positions where nsSNPs destabilize proteins, interfere with the formation of domain-domain interfaces, have an effect on protein-ligand binding or severely impact human health. It currently annotates 28,043 validated SNPs that produce amino acid residue substitutions in human proteins from the SwissProt/TrEMBL database. Annotations can be viewed via a web interface either in the context of a genomic region or by selecting sets of SNPs, genes, proteins or pathways. These results are useful for identifying candidate functional SNPs within a gene, haplotype or pathway and in probing molecular mechanisms responsible for functional impacts of nsSNPs. http://www.salilab.org/LS-SNP CONTACT: rachelk@salilab.org http://salilab.org/LS-SNP/supp-info.pdf.

A single-nucleotide polymorphism tagging set for human drug metabolism and transport

Interindividual variability in drug response, ranging from no therapeutic benefit to life-threatening adverse reactions, is influenced by variation in genes that control the absorption, distribution, metabolism and excretion of drugs. We genotyped 904 single-nucleotide polymorphisms (SNPs) from 55 such genes in two population samples (European and Japanese) and identified a set of tagging SNPs that represents the common variation in these genes, both known and unknown. Extensive empirical evaluations, including a direct assessment of association with candidate functional SNPs in a new, larger population sample, validated the performance of these tagging SNPs and confirmed their utility for linkage-disequilibrium mapping in pharmacogenetics. The analyses also suggest that rare variation is not amenable to tagging strategies.