Developing genomic resources for the common bottlenose dolphin (Tursiops truncatus): isolation and characterization of 153 single nucleotide polymorphisms and 53 genotyping assays

Although single nucleotide polymorphisms (SNPs) are commonly used in human genetics, they have only recently been incorporated into genetic studies of non‐model organisms, including cetaceans. SNPs have several advantages over other molecular markers for studies of population genetics: they are quicker and more straightforward to score, cross‐laboratory comparisons of data are less complicated, and they can be used successfully with low‐quality DNA. We screened portions of the genome of one of the most abundant cetaceans in U.S. waters, the common bottlenose dolphin (Tursiops truncatus), and identified 153 SNPs resulting in an overall average of one SNP every 463 base pairs. Custom TaqMan^® Assays were designed for 53 of these SNPs, and their performance was tested by genotyping a set of bottlenose dolphin samples, including some with low‐quality DNA. We found that in 19% of the loci examined, the minor allele frequency (MAF) estimated during initial SNP ascertainment using a DNA pool of 10 individuals differed significantly from the final MAF after genotyping over 100 individuals, suggesting caution when making inferences about MAF values based on small data sets. For two assays, we also characterized the basis for unusual clustering patterns to determine whether their data could still be utilized for further genetic studies. Overall results support the use of these SNPs for accurate analysis of both poor and good‐quality DNA. We report the first SNP markers and genotyping assays for use in population and conservation genetic studies of bottlenose dolphins.     

Large SNP arrays for genotyping in crop plants

Genotyping with large numbers of molecular markers is now an indispensable tool within plant genetics and breeding. Especially through the identification of large numbers of single nucleotide polymorphism (SNP) markers using the novel high-throughput sequencing technologies, it is now possible to reliably identify many thousands of SNPs at many different loci in a given plant genome. For a number of important crop plants, SNP markers are now being used to design genotyping arrays containing thousands of markers spread over the entire genome and to analyse large numbers of samples. In this article, we discuss aspects that should be considered during the design of such large genotyping arrays and the analysis of individuals. The fact that crop plants are also often autopolyploid or allopolyploid is given due consideration. Furthermore, we outline some potential applications of large genotyping arrays including high-density genetic mapping, characterization (fingerprinting) of genetic material and breeding-related aspects such as association studies and genomic selection.     

Investigating the utility of combining Φ29 whole genome amplification and highly multiplexed single nucleotide polymorphism BeadArray™ genotyping

Background

Sustainable DNA resources and reliable high-throughput genotyping methods are required for large-scale, long-term genetic association studies. In the genetic dissection of common disease it is now recognised that thousands of samples and hundreds of thousands of markers, mostly single nucleotide polymorphisms (SNPs), will have to be analysed. In order to achieve these aims, both an ability to boost quantities of archived DNA and to genotype at low costs are highly desirable. We have investigated Φ29 polymerase Multiple Displacement Amplification (MDA)-generated DNA product (MDA product), in combination with highly multiplexed BeadArray? genotyping technology. As part of a large-scale BeadArray genotyping experiment we made a direct comparison of genotyping data generated from MDA product with that from genomic DNA (gDNA) templates.

Results

Eighty-six MDA product and the corresponding 86 gDNA samples were genotyped at 345 SNPs and a concordance rate of 98.8% was achieved. The BeadArray sample exclusion rate, blind to sample type, was 10.5% for MDA product compared to 5.8% for gDNA.

Conclusions

We conclude that the BeadArray technology successfully produces high quality genotyping data from MDA product. The combination of these technologies improves the feasibility and efficiency of mapping common disease susceptibility genes despite limited stocks of gDNA samples.     

Risk conferring genes in multiple sclerosis

Multiple sclerosis (MS) is characterized by inflammation, axonal and oligodendrocyte pathology and progressive neurological disability. Epidemiologic data indicate that MS may be caused by interplay of genetic and environmental factors. Large samples collected in cooperative efforts and new technologies such as high throughput single nucleotide polymorphism (SNP) genotyping allowed recently to discover non-HLA genes associated with MS susceptibility that are mostly involved in the immune response. In addition, several studies indicate an effect of genetic variations on disease onset, progression and response to therapy. However, the polymorphisms discovered so far explain the genetic variation in MS only in part and are mostly common variants that have only low impact on MS susceptibility. Functional studies are required to validate the importance of the newly identified SNPs. Taking into account the interplay of genetic and environmental factors a combination of genome wide genotyping including HLA-typing and genome wide expression profiling as well as a collection on relevant or putatively relevant environmental factors in patients well characterized clinically and by MRI is a promising way to identify new disease relevant biomarkers.     

Linkage disequilibrium-based quality control for large-scale genetic studies

Quality control (QC) is a critical step in large-scale studies of genetic variation. While, on average, high-throughput single nucleotide polymorphism (SNP) genotyping assays are now very accurate, the errors that remain tend to cluster into a small percentage of "problem" SNPs, which exhibit unusually high error rates. Because most large-scale studies of genetic variation are searching for phenomena that are rare (e.g., SNPs associated with a phenotype), even this small percentage of problem SNPs can cause important practical problems. Here we describe and illustrate how patterns of linkage disequilibrium (LD) can be used to improve QC in large-scale, population-based studies. This approach has the advantage over existing filters (e.g., HWE or call rate) that it can actually reduce genotyping error rates by automatically correcting some genotyping errors. Applying this LD-based QC procedure to data from The International HapMap Project, we identify over 1,500 SNPs that likely have high error rates in the CHB and JPT samples and estimate corrected genotypes. Our method is implemented in the software package fastPHASE, available from the Stephens Lab website (http://stephenslab.uchicago.edu/software.html).     

Effective algorithms for tag SNP selection

Single nucleotide polymorphisms (SNPs), due to their abundance and low mutation rate, are very useful genetic markers for genetic association studies. However, the current genotyping technology cannot afford to genotype all common SNPs in all the genes. By making use of linkage disequilibrium, we can reduce the experiment cost by genotyping a subset of SNPs, called Tag SNPs, which have a strong association with the ungenotyped SNPs, while are as independent from each other as possible. The problem of selecting Tag SNPs is NP-complete; when there are large number of SNPs, in order to avoid extremely long computational time, most of the existing Tag SNP selection methods first partition the SNPs into blocks based on certain block definitions, then Tag SNPs are selected in each block by brute-force search. The size of the Tag SNP set obtained in this way may usually be reduced further due to the inter-dependency among blocks. This paper proposes two algorithms, TSSA and TSSD, to tackle the block-independent Tag SNP selection problem. TSSA is based on A* search algorithm, and TSSD is a heuristic algorithm. Experiments show that TSSA can find the optimal solutions for medium-sized problems in reasonable time, while TSSD can handle very large problems and report approximate solutions very close to the optimal ones.     

SNP-based pathway enrichment analysis for genome-wide association studies

Background  

Recently we have witnessed a surge of interest in using genome-wide association studies (GWAS) to discover the genetic basis of complex diseases. Many genetic variations, mostly in the form of single nucleotide polymorphisms (SNPs), have been identified in a wide spectrum of diseases, including diabetes, cancer, and psychiatric diseases. A common theme arising from these studies is that the genetic variations discovered by GWAS can only explain a small fraction of the genetic risks associated with the complex diseases. New strategies and statistical approaches are needed to address this lack of explanation. One such approach is the pathway analysis, which considers the genetic variations underlying a biological pathway, rather than separately as in the traditional GWAS studies. A critical challenge in the pathway analysis is how to combine evidences of association over multiple SNPs within a gene and multiple genes within a pathway. Most current methods choose the most significant SNP from each gene as a representative, ignoring the joint action of multiple SNPs within a gene. This approach leads to preferential identification of genes with a greater number of SNPs.     

A pipeline for high throughput detection and mapping of SNPs from EST databases

Single nucleotide polymorphisms (SNPs) represent the most abundant type of genetic variation that can be used as molecular markers. The SNPs that are hidden in sequence databases can be unlocked using bioinformatic tools. For efficient application of these SNPs, the sequence set should be error-free as much as possible, targeting single loci and suitable for the SNP scoring platform of choice. We have developed a pipeline to effectively mine SNPs from public EST databases with or without quality information using QualitySNP software, select reliable SNP and prepare the loci for analysis on the Illumina GoldenGate genotyping platform. The applicability of the pipeline was demonstrated using publicly available potato EST data, genotyping individuals from two diploid mapping populations and subsequently mapping the SNP markers (putative genes) in both populations. Over 7000 reliable SNPs were identified that met the criteria for genotyping on the GoldenGate platform. Of the 384 SNPs on the SNP array approximately 12% dropped out. For the two potato mapping populations 165 and 185 SNPs segregating SNP loci could be mapped on the respective genetic maps, illustrating the effectiveness of our pipeline for SNP selection and validation.     

Estimation of relative allele frequencies of single-nucleotide polymorphisms in different populations by microarray hybridization of pooled DNA

Single-nucleotide polymorphisms (SNPs) are considered useful polymorphic markers for genetic studies of polygenic traits. A new practical approach to high-throughput genotyping of SNPs in a large number of individuals is needed in association study and other studies on relationships between genes and diseases. We have developed an accurate and high-throughput method for determining the allele frequencies by pooling the DNA samples and applying a DNA microarray hybridization analysis. In this method, the combination of the microarray, DNA pooling, probe pair hybridization, and fluorescent ratio analysis solves the dual problems of parallel multiple sample analysis, and parallel multiplex SNP genotyping for association study. Multiple DNA samples are immobilized on a slide and a single hybridization is performed with a pool of allele-specific oligonucleotide probes. The results of this study show that hybridization of microarray from pooled DNA samples can accurately obtain estimates of absolute allele frequencies in a sample pool. This method can also be used to identify differences in allele frequencies in distinct populations. It is amenable to automation and is suitable for immediate utilization for high-throughput genotyping of SNP.     

How Many Genetic Variants Remain to Be Discovered?

A great majority of genetic markers discovered in recent genome-wide association studies have small effect sizes, and they explain only a small fraction of the genetic contribution to the diseases. How many more variants can we expect to discover and what study sizes are needed? We derive the connection between the cumulative risk of the SNP variants to the latent genetic risk model and heritability of the disease. We determine the sample size required for case-control studies in order to achieve a certain expected number of discoveries in a collection of most significant SNPs. Assuming similar allele frequencies and effect sizes of the currently validated SNPs, complex phenotypes such as type-2 diabetes would need approximately 800 variants to explain its 40% heritability. Much smaller numbers of variants are needed if we assume rare-variants but higher penetrance models. We estimate that up to 50,000 cases and an equal number of controls are needed to discover 800 common low-penetrant variants among the top 5000 SNPs. Under common and rare low-penetrance models, the very large studies required to discover the numerous variants are probably at the limit of practical feasibility. Under rare-variant with medium- to high-penetrance models (odds-ratios between 1.6 and 4.0), studies comparable in size to many existing studies are adequate provided the genotyping technology can interrogate more and rarer variants.     

SYBR green dye-based probe-free SNP genotyping: Introduction of T-Plex real-time PCR assay

Single-nucleotide polymorphism (SNP) genotyping is widely used in genetic association studies to characterize genetic factors underlying inherited traits. Despite many recent advances in high-throughput SNP genotyping, inexpensive and flexible methods with reasonable throughput levels are still needed. Real-time PCR methods for discovering and genotyping SNPs are becoming increasingly important in various fields of biology. In this study, we introduce a new, single-tube strategy that combines the tetra-primer ARMS PCR assay, SYBR Green I-based real-time PCR, and melting-point analysis with primer design strategies to detect the SNP of interest. This assay, T-Plex real-time PCR, is based on the T_m discrimination of the amplified allele-specific amplicons in a single tube. The specificity, sensitivity, and robustness of the assay were evaluated for common mutations in the FV, PII, MTHFR, and FGFR3 genes. We believe that T-Plex real-time PCR would be a useful alternative for either individual genotyping requests or large epidemiological studies.     

DNA sequence variation and selection of tag single-nucleotide polymorphisms at candidate genes for drought-stress response in Pinus taeda L

Genetic association studies are rapidly becoming the experimental approach of choice to dissect complex traits, including tolerance to drought stress, which is the most common cause of mortality and yield losses in forest trees. Optimization of association mapping requires knowledge of the patterns of nucleotide diversity and linkage disequilibrium and the selection of suitable polymorphisms for genotyping. Moreover, standard neutrality tests applied to DNA sequence variation data can be used to select candidate genes or amino acid sites that are putatively under selection for association mapping. In this article, we study the pattern of polymorphism of 18 candidate genes for drought-stress response in Pinus taeda L., an important tree crop. Data analyses based on a set of 21 putatively neutral nuclear microsatellites did not show population genetic structure or genomewide departures from neutrality. Candidate genes had moderate average nucleotide diversity at silent sites (pi(sil) = 0.00853), varying 100-fold among single genes. The level of within-gene LD was low, with an average pairwise r2 of 0.30, decaying rapidly from approximately 0.50 to approximately 0.20 at 800 bp. No apparent LD among genes was found. A selective sweep may have occurred at the early-response-to-drought-3 (erd3) gene, although population expansion can also explain our results and evidence for selection was not conclusive. One other gene, ccoaomt-1, a methylating enzyme involved in lignification, showed dimorphism (i.e., two highly divergent haplotype lineages at equal frequency), which is commonly associated with the long-term action of balancing selection. Finally, a set of haplotype-tagging SNPs (htSNPs) was selected. Using htSNPs, a reduction of genotyping effort of approximately 30-40%, while sampling most common allelic variants, can be gained in our ongoing association studies for drought tolerance in pine.     

Clinical significance of germline copy number variation in susceptibility of human diseases

Germline copy number variation (CNV) is considered to be an important form of human genetic polymorphisms. Previous studies have identified amounts of CNVs in human genome by advanced technologies, such as comparative genomic hybridization, single nucleotide genotyping, and high-throughput sequencing. CNV is speculated to be derived from multiple mechanisms, such as nonallelic homologous recombination (NAHR) and nonhomologous end-joining (NHEJ). CNVs cover a much larger genome scale than single nucleotide polymorphisms (SNPs), and may alter gene expression levels by means of gene dosage, gene fusion, gene disruption, and long-range regulation effects, thus affecting individual phenotypes and playing crucial roles in human pathogenesis. The number of studies linking CNVs with common complex diseases has increased dramatically in recent years. Here, we provide a comprehensive review of the current understanding of germline CNVs, and summarize the association of germline CNVs with the susceptibility to a wide variety of human diseases that were identified in recent years. We also propose potential issues that should be addressed in future studies.     

Analysis of 17,576 potentially functional SNPs in three case-control studies of myocardial infarction

Myocardial infarction (MI) is a common complex disease with a genetic component. While several single nucleotide polymorphisms (SNPs) have been reported to be associated with risk of MI, they do not fully explain the observed genetic component of MI. We have been investigating the association between MI and SNPs that are located in genes and have the potential to affect gene function or expression. We have previously published studies that tested about 12,000 SNPs for association with risk of MI, early-onset MI, or coronary stenosis. In the current study we tested 17,576 SNPs that could affect gene function or expression. In order to use genotyping resources efficiently, we staged the testing of these SNPs in three case-control studies of MI. In the first study (762 cases, 857 controls) we tested 17,576 SNPs and found 1,949 SNPs that were associated with MI (P<0.05). We tested these 1,949 SNPs in a second study (579 cases and 1159 controls) and found that 24 SNPs were associated with MI (1-sided P<0.05) and had the same risk alleles in the first and second study. Finally, we tested these 24 SNPs in a third study (475 cases and 619 controls) and found that 5 SNPs in 4 genes (ENO1, FXN (2 SNPs), HLA-DPB2, and LPA) were associated with MI in the third study (1-sided P<0.05), and had the same risk alleles in all three studies. The false discovery rate for this group of 5 SNPs was 0.23. Thus, we have identified 5 SNPs that merit further examination for their potential association with MI. One of these SNPs (in LPA), has been previously shown to be associated with risk of cardiovascular disease in other studies.     

SNP analysis to dissect human traits

The analysis of complex human diseases has been spurred by the number of published genomic sequence variants - many identified in the course of sequencing the human genome. But, to be useful for genetic analysis, variants have to be mapped accurately, their frequencies in various populations determined, and automated high-throughput assay techniques developed. Recently proposed methods address these issues: the use of 'reduced representation shotgun' methods for more efficient detection of single nucleotide polymorphisms (SNPs), the employment of high-throughput genotyping techniques, the development of SNP maps that incorporate information about linkage disequilibrium, and the use of SNPs in identifying susceptibility genes for common illnesses.     

Shifting paradigm of association studies: value of rare single-nucleotide polymorphisms

Currently, single-nucleotide polymorphisms (SNPs) with minor allele frequency (MAF) of >5% are preferentially used in case-control association studies of common human diseases. Recent technological developments enable inexpensive and accurate genotyping of a large number of SNPs in thousands of cases and controls, which can provide adequate statistical power to analyze SNPs with MAF <5%. Our purpose was to determine whether evaluating rare SNPs in case-control association studies could help identify causal SNPs for common diseases. We suggest that slightly deleterious SNPs (sdSNPs) subjected to weak purifying selection are major players in genetic control of susceptibility to common diseases. We compared the distribution of MAFs of synonymous SNPs with that of nonsynonymous SNPs (1) predicted to be benign, (2) predicted to be possibly damaging, and (3) predicted to be probably damaging by PolyPhen. Our sources of data were the International HapMap Project, ENCODE, and the SeattleSNPs project. We found that the MAF distribution of possibly and probably damaging SNPs was shifted toward rare SNPs compared with the MAF distribution of benign and synonymous SNPs that are not likely to be functional. We also found an inverse relationship between MAF and the proportion of nsSNPs predicted to be protein disturbing. On the basis of this relationship, we estimated the joint probability that a SNP is functional and would be detected as significant in a case-control study. Our analysis suggests that including rare SNPs in genotyping platforms will advance identification of causal SNPs in case-control association studies, particularly as sample sizes increase.     

Single step genome-wide association studies based on genotyping by sequence data reveals novel loci for the litter traits of domestic pigs

In this study, data genotyping by sequence (GBS) was used to perform single step GWAS (ssGWAS) to identify SNPs associated with the litter traits in domestic pigs and search for candidate genes in the region of significant SNPs. After quality control, 167,355 high-quality SNPs from 532 pigs were obtained. Phenotypic traits on 2112 gilt litters from 532 pigs were recorded including total number born (TNB), number born alive (NBA), and litter weight born alive (LWB). A single-step genomic BLUP approach (ssGBLUP) was used to implement the genome-wide association analysis at a 5% genome-wide significance level. A total of 8, 23 and 20 significant SNPs were associated with TNB, NBA, and LWB, respectively, and these significant SNPs accounted for 62.78%, 79.75%, and 58.79% of genetic variance. Furthermore, 1 (SSC14: 16314857), 4 (SSC1: 81986236, SSC1: 66599775, SSC1: 161999013, and SSC1: 267883107), and 5 (SSC9: 29030061, SSC2: 32368561, SSC5: 110375350, SSC13: 45619882 and SSC13: 45647829) significant SNPs for TNB, NBA, and LWB were inferred to be novel loci. At SSC1, the AIM1 and FOXO3 genes were found to be associated with NBA; these genes increase ovarian reproductive capacity and follicle number and decrease gonadotropin levels. The genes SLC36A4 and INTU are involved in cell growth, cytogenesis and development were found to be associated with LWB. These significant SNPs can be used as an indication for regions in the Sus scrofa genome for variability in litter traits, but further studies are expected to confirm causative mutations.     

Prioritize and select SNPs for association studies with multi-stage designs.

Large-scale whole genome association studies are increasingly common, due in large part to recent advances in genotyping technology. With this change in paradigm for genetic studies of complex diseases, it is vital to develop valid, powerful, and efficient statistical tools and approaches to evaluate such data. Despite a dramatic drop in genotyping costs, it is still expensive to genotype thousands of individuals for hundreds of thousands single nucleotide polymorphisms (SNPs) for large-scale whole genome association studies. A multi-stage (or two-stage) design has been a promising alternative: in the first stage, only a fraction of samples are genotyped and tested using a dense set of SNPs, and only a small subset of markers that show moderate associations with the disease will be genotyped in later stages. Multi-stage designs have also been used in candidate gene association studies, usually in regions that have shown strong signals by linkage studies. To decide which set of SNPs to be genotyped in the next stage, a common practice is to utilize a simple test (such as a chi2 test for case-control data) and a liberal significance level without corrections for multiple testing, to ensure that no true signals will be filtered out. In this paper, I have developed a novel SNP selection procedure within the framework of multi-stage designs. Based on data from stage 1, the method explicitly explores correlations (linkage disequilibrium) among SNPs and their possible interactions in determining the disease phenotype. Comparing with a regular multi-stage design, the approach can select a much reduced set of SNPs with high discriminative power for later stages. Therefore, not only does it reduce the genotyping cost in later stages, it also increases the statistical power by reducing the number of tests. Combined analysis is proposed to further improve power, and the theoretical significance level of the combined statistic is derived. Extensive simulations have been performed, and results have shown that the procedure can reduce the number of SNPs required in later stages, with improved power to detect associations. The procedure has also been applied to a real data set from a genome-wide association study of the sporadic amyotrophic lateral sclerosis (ALS) disease, and an interesting set of candidate SNPs has been identified.     

Evaluating coverage of exons by HapMap SNPs

Genome-wide association (GWA) studies are currently one of the most powerful tools in identifying disease-associated genes or variants. In typical GWA studies, single-nucleotide polymorphisms (SNPs) are often used as genetic makers. Therefore, it is critical to estimate the percentage of genetic variations which can be covered by SNPs through linkage disequilibrium (LD). In this study, we use the concept of haplotype blocks to evaluate the coverage of five SNP sets including the HapMap and four commercial arrays, for every exon in the human genome. We show that although some Chips can reach similar coverage as the HapMap, only about 50% of exons are completely covered by haplotype blocks of HapMap SNPs. We suggest further high-resolution genotyping methods are required, to provide adequate genome-wide power for identifying variants.     

SNPselector: a web tool for selecting SNPs for genetic association studies

SUMMARY: Single nucleotide polymorphisms (SNPs) are commonly used for association studies to find genes responsible for complex genetic diseases. With the recent advance of SNP technology, researchers are able to assay thousands of SNPs in a single experiment. But the process of manually choosing thousands of genotyping SNPs for tens or hundreds of genes is time consuming. We have developed a web-based program, SNPselector, to automate the process. SNPselector takes a list of gene names or a list of genomic regions as input and searches the Ensembl genes or genomic regions for available SNPs. It prioritizes these SNPs on their tagging for linkage disequilibrium, SNP allele frequencies and source, function, regulatory potential and repeat status. SNPselector outputs result in compressed Excel spreadsheet files for review by the user. AVAILABILITY: SNPselector is freely available at http://primer.duhs.duke.edu/