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. ^

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.

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.

Heterozygosity and diversity analysis using mapped single nucelotide polymorphisms in a faba bean inbreeding programme

The aim of this study was to convert existing faba bean (Vicia faba L.) single nucleotide polymorphism (SNP) markers from cleaved amplification polymorphic sequences and SNaPshot® formats, which are expensive and time-consuming, to the more convenient KBiosciences competitive allele‐specific PCR (KASP) assay format. Out of 80 assays designed, 75 were validated, though a core set of 67 of the most robust markers is recommended for further use. The 67 best KASP SNP assays were used across two generations of single seed descent to detect unintended outcrossing and to track and quantify loss of heterozygosity, a capability that will significantly increase the efficiency and performance of pure line production and maintenance. This same set of assays was also used to examine genetic relationships between the 67 members of the partly inbred panel, and should prove useful for line identification and diversity studies in the future.

An interleukin-1 beta (IL-1 beta) single-nucleotide polymorphism at position 3954 and red complex periodontopathogens independently and additively modulate the levels of IL-1 beta in diseased periodontal tissues

Inflammatory cytokines such as interieukin-1 beta (IL-1 beta) are involved in the pathogenesis of periodontal diseases. A high individual variation in the levels of IL-10 mRNA has been verified, which is possibly determined by genetic polymorphisms and/or by the presence of periodontopathogens such as Porphyromonas gingivalis, Tannerella forsythia, Treponema denticola, and Aggregatibacter actinomycetemcomitans. In this study, we investigated the role of an IL-10 promoter single-nucleotide polymorphism at position 3954 [IL-1 beta(3954) SNP] and the presence of the periodontopathogens in the determination of the IL-1 beta levels in the periodontal tissues of nonsmoking chronic periodontitis (CP) patients (n = 117) and control (C) subjects in = 175) and the possible correlations with the clinical parameters of the disease. IL-1 beta(3954) SNP was investigated by restriction fragment length polymorphism, while the IL-1 beta levels and the presence of the periodontopathogens were determined by real-time PCR. Similar frequencies of IL-1 beta(3954) SNP were found in the C and CP groups, in spite of a trend toward a higher incidence of T alleles in the CP group. The IL-1 beta (3954) SNP CT and TT genotypes, as well as P. gingivalis, T. forsythia, and T. denticola, were associated with higher IL-1 beta levels and with higher values of the clinical parameters of disease severity. Concomitant analyses demonstrate that IL-1 beta(3954) and the red complex periodontopathogens were found to independently and additively modulate the levels of IL-1 beta in periodontal tissues. Similarly, the concurrent presence of both factors was associated with increased scores of disease severity. IL-1 beta(3954) genotypes and red complex periodontopathogens, individually and additively, modulate the levels of IL-1 beta in the diseased tissues of nonsmoking CP patients and, consequently, are potentially involved in the determination of the disease outcome.

Tumor necrosis factor-alpha-308G/A single nucleotide polymorphism and red-complex periodontopathogens are independently associated with increased levels of tumor necrosis factor-alpha in diseased periodontal tissues

Background and Objective: Inflammatory cytokines such as tumor necrosis factor-alpha are involved in the pathogenesis of periodontal diseases. A high between-subject variation in the level of tumor necrosis factor-alpha mRNA has been verified, which may be a result of genetic polymorphisms and/or the presence of periodontopathogens such as Porphyromonas gingivalis, Tannerella forsythia, Treponema denticola (called the red complex) and Aggregatibacter actinomycetemcomitans. In this study, we investigated the effect of the tumor necrosis factor-alpha (TNFA) -308G/A gene polymorphism and of periodontopathogens on the tumor necrosis factor-alpha levels in the periodontal tissues of nonsmoking patients with chronic periodontitis (n = 127) and in control subjects (n = 177). Material and Methods: The TNFA-308G/A single nucleotide polymorphism was investigated using polymerase chain reaction-restriction fragment length polymorphism analysis, whereas the tumor necrosis factor-alpha levels and the periodontopathogen load were determined using real-time polymerase chain reaction. Results: No statistically significant differences were found in the frequency of the TNFA-308 single nucleotide polymorphism in control and chronic periodontitis groups, in spite of the higher frequency of the A allele in the chronic periodontitis group. The concomitant analyses of genotypes and periodontopathogens demonstrated that TNFA-308 GA/AA genotypes and the red-complex periodontopathogens were independently associated with increased levels of tumor necrosis factor-alpha in periodontal tissues, and no additive effect was seen when both factors were present. P. gingivalis, T. forsythia and T. denticola counts were positively correlated with the level of tumor necrosis factor-alpha. TNFA-308 genotypes were not associated with the periodontopathogen detection odds or with the bacterial load. Conclusion: Our results demonstrate that the TNFA-308 A allele and red-complex periodontopathogens are independently associated with increased levels of tumor necrosis factor-alpha in diseased tissues of nonsmoking chronic periodontitis patients and consequently are potentially involved in determining the disease outcome.

Effects of a single nucleotide polymorphism in the leptin gene on the productive traits of dairy buffaloes (Bubalus bubalis)

The gene responsible for coding the leptin hormone has been associated with productive and reproductive traits in cattle. In dairy cattle, different polymorphisms found in the leptin gene have been associated with several traits of economic interest, such as energy balance, milk yield and composition, live weight, fertility and dry matter consumption. The aim of this study was to detect genetic variability in the leptin gene of buffaloes and to test possible associations with milk yield, fat and protein percentages, age at first calving and first calving interval. Three genotypes (AA, AG and GG) were identified by polymerase chain reaction-restriction fragment length polymorphism, which presented genotypic frequencies of 0.30, 0.54 and 0.16, respectively. The allele frequencies were 0.57 for the A allele and 0.43 for the G allele. No significant effects were found in the present study, but there is an indicative that leptin gene affects lipid metabolism. © 2013 Springer Science+Business Media Dordrecht.

Single-nucleotide polymorphism in two representative multidrug-resistant Mycobacterium bovis isolates collected from patients in a Spanish hospital harboring a human infection outbreak

Mycobacterium bovis is the etiological agent of tuberculosis in domestic and wild animals. Its involvement as a human pathogen has been highlighted again with the recent descriptions of transmission through dairy products (18), reactivation or primary infection in human immunodeficiency virus-infected patients (5), and association with meat industry workers, animal keepers, or hunters (3). Strains resistant to antituberculous drugs (M. bovis is naturally resistant to pyrazinamide) pose an additional risk (2). Several studies have demonstrated that mutations in target genes are associated with resistance to antituberculous drugs (4, 7, 10, 11, 16). However, most of them have been developed in Mycobacterium tuberculosis strains and limited data are available regarding M. bovis isolates. The aim of this study was to characterize by sequencing the main genes involved in antibiotic resistance in two multidrug-resistant (MDR) M. bovis isolates in a human outbreak detected in a hospital in Madrid that subsequently spread to several countries (5, 6, 15). The isolates were resistant to 11 drugs, but only their rpoB and katG genes have been analyzed so far (1, 14). We studied the first (93/R1) and last (95/R4) M. bovis isolates of this nosocomial outbreak, characterized by spoligotyping as SB0426 (hexacode 63-5F-5E-7F-FF-60 in the database at www.mbovis.org) (1, 13). Several genes involved in resistance to isoniazid (katG, ahpC, inhA, and the oxyR-ahpC intergenic region), rifampin (rpoB), streptomycin (rrs, rpsL), ethambutol (embB), and quinolones (gyrA) were studied. These genes, or fragments of genes, were amplified and sequenced as previously described (12). The sequence analysis revealed polymorphisms in five (ahpC, rpoB, rpsL, embB, and gyrA) out of nine analyzed genes (Table 1). Nucleotide substitutions in four genes cause a change in the encoded amino acid. Two additional synonymous mutations in ahpC and rpsL differentiated the first and last isolates from the outbreak.

Highly discriminatory single-nucleotide polymorphism interrogation of Escherichia coli by use of allele-specific real-time PCR and eBURST analysis

In total, 782 Escherichia coli strains originating from various host sources have been analyzed in this study by using a highly discriminatory single-nucleotide polymorphism (SNP) approach. A set of eight SNPs, with a discrimination value (Simpson's index of diversity [D]) of 0.96, was determined using the Minimum SNPs software, based on sequences of housekeeping genes from the E. coli multilocus sequence typing (MLST) database. Allele-specific real-time PCR was used to screen 114 E. coli isolates from various fecal sources in Southeast Queensland (SEQ). The combined analysis of both the MLST database and SEQ E. coli isolates using eight high-D SNPs resolved the isolates into 74 SNP profiles. The data obtained suggest that SNP typing is a promising approach for the discrimination of host-specific groups and allows for the identification of human-specific E. coli in environmental samples. However, a more diverse E. coli collection is required to determine animal- and environment-specific E. coli SNP profiles due to the abundance of human E. coli strains (56%) in the MLST database.

Association between polymorphisms in the progesterone receptor gene and endometriosis

The progesterone receptor (PR) is a candidate gene for the development of endometriosis, a complex disease with strong hormonal features, common in women of reproductive age. We typed the 306 base pair Alu insertion (AluIns) polymorphism in intron G of PR in 101 individuals, estimated linkage disequilibrium (LD) between five single-nucleotide polymorphisms (SNPs) across the PR locus in 980 Australian triads (endometriosis case and two parents) and used transmission disequilibrium testing (TDT) for association with endometriosis. The five SNPs showed strong pairwise LD, and the AluIns was highly correlated with proximal SNPs rs1042839 (Δ2 = 0.877, D9 = 1.00, P < 0.0001) and rs500760 (Δ2 = 0.438, D9 = 0.942, P < 0.0001). TDT showed weak evidence of allelic association between endometriosis and rs500760 (P = 0.027) but not in the expected direction. We identified a common susceptibility haplotype GGGCA across the five SNPs (P = 0.0167) in the whole sample, but likelihood ratio testing of haplotype transmission and non-transmission of the AluIns and flanking SNPs showed no significant pattern. Further, analysis of our results pooled with those from two previous studies suggested that neither the T2 allele of the AluIns nor the T1/T2 genotype was associated with endometriosis.

Human-specific E. coli single nucleotide polymorphism (SNP) genotypes detected in a South-East Queensland waterway, Australia

The World Health Organization recommends that the majority of water monitoring laboratories in the world should test for E. coli daily since thermotolerant coliforms and E. coli are key indicators for risk assessment of recreational waters. Recently, we developed a new SNP method for typing E. coli strains, by which human-specific genotypes were identified. Here, we report the presence of these previously described specific SNP profiles in environmental water, sourced from the Coomera River, located on South East Queensland, Australia, over a period of two years. This study tested for the presence of human-specific E. coli to ascertain whether hydrologic and anthropogenic activity plays a key role in the pollution of the investigated watershed or whether the pollution is from other sources. We found six human-specific SNP profiles and one animal-specific SNP profile consistently across sampling sites and times. We have demonstrated that our SNP genotyping method is able to rapidly identify and characterise human- and animal-specific E. coli isolates in water sources.

Genetic polymorphisms in the human tissue kallikrein (KLK) locus and their implication in various malignant and non-malignant diseases

The Kallikrein (KLK) gene locus encodes a family of serine proteases and is the largest contiguous cluster of protease-encoding genes attributed an evolutionary age of 330 million years. The KLK locus has been implicated as a high susceptibility risk loci in numerous cancer studies through the last decade. The KLK3 gene already has established clinical relevance as a biomarker in prostate cancer prognosis through its encoded protein, prostate-specific antigen. Data mined through genome-wide association studies (GWAS) and next-generation sequencing point to many important candidate single nucleotide polymorphisms (SNPs) in KLK3 and other KLK genes. SNPs in the KLK locus have been found to be associated with several diseases including cancer, hypertension, cardiovascular disease and atopic dermatitis. Moreover, introducing a model incorporating SNPs to improve the efficiency of prostate-specific antigen in detecting malignant states of prostate cancer has been recently suggested. Establishing the functional relevance of these newly-discovered SNPs, and their interactions with each other, through in silico investigations followed by experimental validation, can accelerate the discovery of diagnostic and prognostic biomarkers. In this review, we discuss the various genetic association studies on the KLK loci identified either through candidate gene association studies or at the GWAS and post-GWAS front to aid researchers in streamlining their search for the most significant, relevant and therapeutically promising candidate KLK gene and/or SNP for future investigations.

No role for estrogen receptor 1 gene intron 1 Pvu II and exon 4 C325G polymorphisms in migraine susceptibility

Background We have previously reported an association between the estrogen receptor 1 (ESR1) gene exon 8 G594A polymorphism and migraine susceptibility in two independent Australian cohorts. In this paper we report results of analysis of two further single nucleotide polymorphisms (SNPs) in the ESR1 gene in the same study group, the T/C Pvu II SNP in intron 1 and the C325G SNP in exon 4, as well as results of linkage disequilibrium (LD) analysis on these markers. Methods We investigated these variants by case-control association analysis in a cohort of 240 migraineurs and 240 matched controls. The SNPs were genotyped using specific restriction enzyme assays. Results were analysed using contingency table methods incorporating the chi-squared statistic. LD results are presented as D' statistics with associated P values. Results We found no evidence for association of the Pvu II T/C polymorphism and the C325G polymorphism and migraine susceptibility and no evidence for LD between these two SNPs and the previously implicated exon 8 G594A marker. Conclusion We have found no role for the polymorphisms in intron 1 and exon 4 with migraine susceptibility. To further investigate our previously implicated exon 8 marker, we suggest the need for studies with a high density of polymorphisms be undertaken, with particular focus on markers in LD with the exon 8 marker.

Genetic polymorphisms in DPF3 associated with risk of breast cancer and lymph node metastases

Background Several studies have identified rare genetic variations responsible for many cases of familial breast cancer but their contribution to total breast cancer incidence is relatively small. More common genetic variations with low penetrance have been postulated to account for a higher proportion of the population risk of breast cancer. Methods and Results In an effort to identify genes that influence non-familial breast cancer risk, we tested over 25,000 single nucleotide polymorphisms (SNPs) located within approximately 14,000 genes in a large-scale case-control study in 254 German women with breast cancer and 268 age-matched women without malignant disease. We identified a marker on chromosome 14q24.3-q31.1 that was marginally associated with breast cancer status (OR = 1.5, P = 0.07). Genotypes for this SNP were also significantly associated with indicators of breast cancer severity, including presence of lymph node metastases ( P = 0.006) and earlier age of onset ( P = 0.01). The association with breast cancer status was replicated in two independent samples (OR = 1.35, P = 0.05). High-density association fine mapping showed that the association spanned about 80 kb of the zinc-finger gene DPF3 (also known as CERD4 ). One SNP in intron 1 was found to be more strongly associated with breast cancer status in all three sample collections (OR = 1.6, P = 0.003) as well as with increased lymph node metastases ( P = 0.01) and tumor size ( P = 0.01). Conclusion Polymorphisms in the 5' region of DPF3 were associated with increased risk of breast cancer development, lymph node metastases, age of onset, and tumor size in women of European ancestry. This large-scale association study suggests that genetic variation in DPF3 contributes to breast cancer susceptibility and severity.