Single nucleotide polymorphism (SNP) selection for genotype-phenotype association studies

With hundreds of single nucleotide polymorphisms (SNPs) in a candidate gene and millions of SNPs across the genome, selecting an informative subset of SNPs to maximize the ability to detect genotype-phenotype association is of great interest and importance. In addition, with a large number of SNPs, analytic methods are needed that allow investigators to control the false positive rate resulting from large numbers of SNP genotype-phenotype analyses. This dissertation uses simulated data to explore methods for selecting SNPs for genotype-phenotype association studies. I examined the pattern of linkage disequilibrium (LD) across a candidate gene region and used this pattern to aid in localizing a disease-influencing mutation. The results indicate that the r2 measure of linkage disequilibrium is preferred over the common D′ measure for use in genotype-phenotype association studies. Using step-wise linear regression, the best predictor of the quantitative trait was not usually the single functional mutation. Rather it was a SNP that was in high linkage disequilibrium with the functional mutation. Next, I compared three strategies for selecting SNPs for application to phenotype association studies: based on measures of linkage disequilibrium, based on a measure of haplotype diversity, and random selection. The results demonstrate that SNPs selected based on maximum haplotype diversity are more informative and yield higher power than randomly selected SNPs or SNPs selected based on low pair-wise LD. The data also indicate that for genes with small contribution to the phenotype, it is more prudent for investigators to increase their sample size than to continuously increase the number of SNPs in order to improve statistical power. When typing large numbers of SNPs, researchers are faced with the challenge of utilizing an appropriate statistical method that controls the type I error rate while maintaining adequate power. We show that an empirical genotype based multi-locus global test that uses permutation testing to investigate the null distribution of the maximum test statistic maintains a desired overall type I error rate while not overly sacrificing statistical power. The results also show that when the penetrance model is simple the multi-locus global test does as well or better than the haplotype analysis. However, for more complex models, haplotype analyses offer advantages. The results of this dissertation will be of utility to human geneticists designing large-scale multi-locus genotype-phenotype association studies. ^

The human XPG gene: gene architecture, alternative splicing and single nucleotide polymorphisms

Defects in the XPG DNA repair endonuclease gene can result in the cancer-prone disorders xeroderma pigmentosum (XP) or the XP–Cockayne syndrome complex. While the XPG cDNA sequence was known, determination of the genomic sequence was required to understand its different functions. In cells from normal donors, we found that the genomic sequence of the human XPG gene spans 30 kb, contains 15 exons that range from 61 to 1074 bp and 14 introns that range from 250 to 5763 bp. Analysis of the splice donor and acceptor sites using an information theory-based approach revealed three splice sites with low information content, which are components of the minor (U12) spliceosome. We identified six alternatively spliced XPG mRNA isoforms in cells from normal donors and from XPG patients: partial deletion of exon 8, partial retention of intron 8, two with alternative exons (in introns 1 and 6) and two that retained complete introns (introns 3 and 9). The amount of alternatively spliced XPG mRNA isoforms varied in different tissues. Most alternative splice donor and acceptor sites had a relatively high information content, but one has the U12 spliceosome sequence. A single nucleotide polymorphism has allele frequencies of 0.74 for 3507G and 0.26 for 3507C in 91 donors. The human XPG gene contains multiple splice sites with low information content in association with multiple alternatively spliced isoforms of XPG mRNA.

The Genetics of Success: How Single-Nucleotide Polymorphisms Associated With Educational Attainment Relate to Life-Course Development.

A previous genome-wide association study (GWAS) of more than 100,000 individuals identified molecular-genetic predictors of educational attainment. We undertook in-depth life-course investigation of the polygenic score derived from this GWAS using the four-decade Dunedin Study (N = 918). There were five main findings. First, polygenic scores predicted adult economic outcomes even after accounting for educational attainments. Second, genes and environments were correlated: Children with higher polygenic scores were born into better-off homes. Third, children's polygenic scores predicted their adult outcomes even when analyses accounted for their social-class origins; social-mobility analysis showed that children with higher polygenic scores were more upwardly mobile than children with lower scores. Fourth, polygenic scores predicted behavior across the life course, from early acquisition of speech and reading skills through geographic mobility and mate choice and on to financial planning for retirement. Fifth, polygenic-score associations were mediated by psychological characteristics, including intelligence, self-control, and interpersonal skill. Effect sizes were small. Factors connecting DNA sequence with life outcomes may provide targets for interventions to promote population-wide positive development.

Development of single nucleotide polymorphism (SNP) markers from the mango (Mangifera indica) transcriptome for mapping and estimation of genetic diversity

The development of molecular markers for genomic studies in Mangifera indica (mango) will allow marker-assisted selection and identification of genetically diverse germplasm, greatly aiding mango breeding programs. We report here our identification of thousands of unambiguous molecular markers that can be easily assayed across genotypes of the species. With origin centered in Southeast Asia, mangos are grown throughout the tropics and subtropics as a nutritious fruit that exhibits remarkable intraspecific phenotypic diversity. With the goal of building a high density genetic map, we have undertaken discovery of sequence variation in expressed genes across a broad range of mango cultivars. A transcriptome sequence reference was built de novo from extensive sequencing and assembly of RNA from cultivar 'Tommy Atkins'. Single nucleotide polymorphisms (SNPs) in protein coding transcripts were determined from alignment of RNA reads from 24 mango cultivars of diverse origins: 'Amin Abrahimpur' (India), 'Aroemanis' (Indonesia), 'Burma' (Burma), 'CAC' (Hawaii), 'Duncan' (Florida), 'Edward' (Florida), 'Everbearing' (Florida), 'Gary' (Florida), 'Hodson' (Florida), 'Itamaraca' (Brazil), 'Jakarata' (Florida), 'Long' (Jamaica), 'M. Casturi Purple' (Borneo), 'Malindi' (Kenya), 'Mulgoba' (India), 'Neelum' (India), 'Peach' (unknown), 'Prieto' (Cuba), 'Sandersha' (India), 'Tete Nene' (Puerto Rico), 'Thai Everbearing' (Thailand), 'Toledo' (Cuba), 'Tommy Atkins' (Florida) and 'Turpentine' (West Indies). SNPs in a selected subset of protein coding transcripts are currently being converted into Fluidigm assays for genotyping of mapping populations and germplasm collections. Using an alternate approach, SNPs (144) discovered by sequencing of candidate genes in 'Kensington Pride' have already been converted and used for genotyping.

Single nucleotide polymorphisms in the IS900 sequence of Mycobacterium avium subsp. paratuberculosis are strain type specific

Insertion sequence IS900 is used as a target for the identification of Mycobacterium avium subsp. paratuberculosis. Previous reports have revealed single nucleotide polymorphisms within IS900. This study, which analyzed the IS900 sequences of a panel of isolates representing M. avium subsp. paratuberculosis strain types I, II, and III, revealed conserved type-specific polymorphisms that could be utilized as a tool for diagnostic and epidemiological purposes.

A single nucleotide polymorphism in NEUROD1 is associated with production traits in nelore beef cattle.

Feed efficiency and carcass characteristics are late-measured traits. The detection of molecular markers associated with them can help breeding programs to select animals early in life, and to predict breeding values with high accuracy. The objective of this study was to identify polymorphisms in the functional and positional candidate gene NEUROD1 (neurogenic differentiation 1), and investigate their associations with production traits in reference families of Nelore cattle. A total of 585 steers were used, from 34 sires chosen to represent the variability of this breed. By sequencing 14 animals with extreme residual feed intake (RFI) values, seven single nucleotide polymorphisms (SNPs) in NEUROD1 were identified. The investigation of marker effects on the target traits RFI, backfat thickness (BFT), ribeye area (REA), average body weight (ABW), and metabolic body weight (MBW) was performed with a mixed model using the restricted maximum likelihood method. SNP1062, which changes cytosine for guanine, had no significant association with RFI or REA. However, we found an additive effect on ABW (P ≤ 0.05) and MBW (P ≤ 0.05), with an estimated allele substitution effect of -1.59 and -0.93 kg0.75, respectively. A dominant effect of this SNP for BFT was also found (P ≤ 0.010). Our results are the first that identify NEUROD1 as a candidate that affects BFT, ABW, and MBW. Once confirmed, the inclusion of this SNP in dense panels may improve the accuracy of genomic selection for these traits in Nelore beef cattle as this SNP is not currently represented on SNP chips.

Identificação de SNPs (Single Nucleiotide Polymorphisms) no gene colina monooxigenase relacionado ao metabolismo da glicina betaína em Eucalyptus

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

A single multiplex PCR and SNaPshot minisequencing reaction of 42 SNPs to classify admixture populations into mitochondrial DNA haplogroups

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

A single multiplex SNaPshot reaction of 42 SNPs to classify admixture populations into mitochondrial DNA haplogroups

SNaPshot minisequencing reaction is in increasing use because of its fast detection of many polymorphisms in a single assay. In this work we described a highly sensitive single nucleotide polymorphisms (SNPs) typing method with detection of 42 mitochondrial DNA (mtDNA) SNPs in a single PCR and SNaPshot multiplex reaction in order to allow haplogroup classification in Latin American admixture population. We validated the panel typing 160 Brazilian individuals. DNA was extracted from blood spotted on filter paper using Chelex protocol. Forty SNPs were selected targeting haplogroup-specific mutations in Europeans, Africans and Asians (only precursors of Native Americans haplogroups A2, B2, C1, and D1) and two non-coding SNPs were chosen to increase the power of discrimination between individuals (SNPs positions 16,519 and 16,362). It was done using a modified version of a previously published multiplex SNaPshot minisequencing reaction established to resolve European haplogroups, adding SNPs targeting Africans (L0, L1, L2, L3, and L*) and Asians (A, B, C, and D) haplogroups based on SNPs described at PhyloTree.org build 2. PCR primers were designed using PerlPrimer software and checked with the Autodimer program. Thirty-three primer-pairs were used to amplify 42 SNPs. Using this panel, we were able to successfully classify 160 individuals into their correct haplogroups. Complete SNP profiles were obtained from 10. pg of total DNA. We conclude that it is possible to build and genotype more than 40 mtDNA SNPs in a single multiplex PCR and SNaPshot reaction, with sensitivity and reliability, resolving haplogroup classification in admixture populations. © 2011.

Association of SNPs on CAPN1 and CAST genes with tenderness in Nellore cattle

We examined whether single-nucleotide polymorphisms (SNPs) in the calpain (CAPN) and calpastatin (CAST) genes, described from Bos primigenius taurus, are polymorphic in Nellore cattle. We also looked for a possible association of linkage disequilibrium of this polymorphism with tenderness of the longissimus dorsi muscle after 7, 14 and 21 days of postmortem aging in 638 purebred Nellore bulls. Meat tenderness was measured as Warner-Bratzler shear force. Additive and dominance effects were tested for SNPs of the three genotypic classes; the substitution effect was tested for SNPs with missing genotypic classes. Genotypic and gene frequencies were also calculated for the different SNPs. An increase in tenderness was observed from 7 to 21 days; the average values for shear force at 7, 14 and 21 days of aging were 5.92 +/- 0.06, 4.92 +/- 0.05, and 4.38 +/- 0.04 kg, respectively. All markers showed polymorphism, but there was no CC genotype for CAPN316, and few animals showed the AA genotype for CAPN530. The alleles CAPN4751, UOGCAST1, and WSUCAST were found to have additive and dominance effects for shear force at 7, 14 and 21 days, while CAPN316 showed a substitution effect for shear force at 7 and 21 days. An additive-by-additive epistatic interaction was observed between CAPN4751 and markers on the CAST gene. In conclusion, these markers should be considered for use in breeding programs.

CYP1A2*1C, CYP2E1*5B, and GSTM1 polymorphisms are predictors of risk and poor outcome in head and neck squamous cell carcinoma patients

Head and neck squamous cell carcinoma (HNSCC) is associated with environmental factors, especially tobacco and alcohol consumption. Most of the carcinogens present in tobacco smoke are converted into DNA-reactive metabolites by cytochrome P450 (CYPs) enzymes and detoxification of these substances is performed by glutathione S-transferases (GSTs). It has been suggested that genetic alterations, such as polymorphisms, play an important role in tumorigenesis and HNSCC progression. The aim of this study was to investigate CYP1A1, CYP1A2, CYP2E1, GSTM1, and GSTT1 polymorphisms as risk factors in HNSCC and their association with clinicopathologic data. The patients comprised 153 individuals with HNSCC (cases) and 145 with no current or previous diagnosis of cancer (controls). Genotyping of the single nucleotide polymorphisms (SNPs) of the CYP1A1, CYP1A2, and CYP2E1 genes was performed by PCR-RFLP and the GSTM1 and GSTT1 copy number polymorphisms (CNPs) were analyzed by PCR-multiplex. As expected, a significant difference was detected for tobacco and alcohol consumption between cases and controls (P < 0.001). It was observed that the CYP1A2*1D (OR = 16.24) variant and GSTM1 null alleles (OR = 0.02) confer increased risk of HNSCC development (P < 0.001). In addition, head and neck cancer alcohol consumers were more frequently associated with the CYP2E1*5B variant allele than control alcohol users (P < 0.0001, OR = 190.6). The CYP1A2*1C polymorphism was associated with tumor recurrence (log-rank test, P = 0.0161). The CYP2E1*5B and GSTM1 null alleles were significantly associated with advanced clinical stages (T3 + T4; P = 0.022 and P = 0.028, respectively). Overall, the findings suggested that the genetic polymorphisms studied are predictors of risk and are also associated with tumor recurrence, since they are important for determining the parameters associated with tumor progression and poor outcomes in HNSCC. (C) 2009 Elsevier Ltd. All rights reserved.

DNA polymorphisms as tools for spinal cord injury research

Study Design: Data mining of single nucleotide polymorphisms (SNPs) in gene pathways related to spinal cord injury (SCI). Objectives: To identify gene polymorphisms putatively implicated with neuronal damage evolution pathways, potentially useful to SCI study. Setting: Departments of Psychiatry and Orthopedics, Faculdade de Medicina, Universidade de Sao Paulo, Brazil. Methods: Genes involved with processes related to SCI, such as apoptosis, inflammatory response, axonogenesis, peripheral nervous system development and axon ensheathment, were determined by evaluating the `Biological Process` annotation of Gene Ontology (GO). Each gene of these pathways was mapped using MapViewer, and gene coordinates were used to identify their polymorphisms in the SNP database. As a proof of concept, the frequency of subset of SNPs, located in four genes (ALOX12, APOE, BDNF and NINJ1) was evaluated in the DNA of a group of 28 SCI patients and 38 individuals with no SC lesions. Results: We could identify a total of 95 276 SNPs in a set of 588 genes associated with the selected GO terms, including 3912 nucleotide alterations located in coding regions of genes. The five non-synonymous SNPs genotyped in our small group of patients, showed a significant frequency, reinforcing their potential use for the investigation of SCI evolution. Conclusion: Despite the importance of SNPs in many aspects of gene expression and protein activity, these gene alterations have not been explored in SCI research. Here we describe a set of potentially useful SNPs, some of which could underlie the genetic mechanisms involved in the post trauma spinal cord damage.