A simple and accurate method for determination of microsatellite total allele content differences between DNA pools

DNA pooling is a potential tool for the efficient analysis of the large numbers of samples and DNA markers that are necessary for genome-wide association studies. A simple accurate method for measuring total allele differences in comparisons between two pools containing large numbers of DNA samples is presented. This method compares relative peak height differences between electrophoretograms for each allele of a microsatellite. The method was evaluated by the analysis of 11 microsatellite markers and DNA pooled sample sizes of 50, 100, and 200 individual DNA samples from the same number of different subjects. Pools were created from previously individually genotyped subjects and constructed so that the pool comparisons would provide real total allele differences varying from 0% to 55%. Calculated pool differences were then compared with the real total allele differences determined by individual genotyping results. Together over 200 comparisons demonstrated a correlation coefficient of 0.96, which compared favorably with other previous methods of analysis. This method could provide a rapid screen for total allele differences of greater than 10%, a threshold that should be applicable to detecting low relative risk genes in common diseases. Therefore, these studies suggest that DNA pooling could be a useful tool in association studies for the determination of candidate regions for a range of complex genetic diseases.     

Novel linkage mapping approach using DNA pooling in human and animal genetics. I. Detection of complex disease loci

DNA pooling is a potential methodology for genetic loci with small effect contributing to complex diseases and quantitative traits. This is accomplished by the rapid preliminary screening of the genome for the allelic association with the most common class of polymorphic short tandem repeat markers. The methodology assumes as a common founder for the linked disease locus of interest and searches for a region of a chromosome shared between affected individuals. The general theory of DNA pooling basically relies on the observed differences in the allelic distribution between pools from affected and unaffected individuals, including a reduction in the number of alleles in the affected pool, which indicate the sharing of a chromosomal region. The power of statistic for associated linkage mapping can be determined using two recently developed strategies, firstly, by measuring the differences of allelic image patterns produced by two DNA pools of extreme character and secondly, by measuring total allele content differences by comparing between two pools containing large numbers of DNA samples. These strategies have effectively been utilized to identify the shared chromosomal regions for linkage studies and to investigate the candidate disease loci for fine structure gene mapping using allelic association. This paper outlines the utilization of DNA pooling as a potential tool to locate the complex disease loci, statistical methods for accurate estimates of allelic frequencies from DNA pools, its advantages, drawbacks and significance in associate linkage mapping using pooled DNA samples.     

Determination of SNP allele frequencies in pooled DNAs by primer extension genotyping and denaturing high-performance liquid chromatography

By testing DNA pools rather than single samples the number of tests for a case-control association study can be decreased to only two for each marker: one on the patient and one on the control pool. A fundamental requirement is that each pool represents the frequency of the markers in the corresponding population beyond the influence of experimental errors. Consequently the latter must be carefully determined. To this aim, we prepared pools of different size (49-402 individuals) with accurately quantified DNAs, estimated the allelic frequencies in the pools of two SNPs by primer extension genotyping followed by DHPLC analysis and compared them with the real frequencies determined in the single samples. Our data show that (1) the method is highly reproducible: the standard deviation of repeated determinations was +/-0.014; (2) the experimental error (i.e., the discrepancy between the estimated and real frequencies) was +/-0.013 (95% C.I.: 0.0098-0.0165). The magnitude of this error was not correlated to the pool size or to the type of SNP. The effect of the observed experimental error on the power of the association test was evaluated. We conclude that this method constitutes an efficient tool for high-throughput association screenings provided that the experimental error is low. We therefore recommend that before a pool is used for extensive association studies, its quality, i.e., the experimental error, is verified by determining the difference between estimated and real frequencies for at least one marker.     

Evaluation of DNA pooling for the estimation of microsatellite allele frequencies: a case study using striped bass (Morone saxatilis)

Using striped bass (Morone saxatilis) and six multiplexed microsatellite markers, we evaluated procedures for estimating allele frequencies by pooling DNA from multiple individuals, a method suggested as cost-effective relative to individual genotyping. Using moment-based estimators, we estimated allele frequencies in experimental DNA pools and found that the three primary laboratory steps, DNA quantitation and pooling, PCR amplification, and electrophoresis, accounted for 23, 48, and 29%, respectively, of the technical variance of estimates in pools containing DNA from 2-24 individuals. Exact allele-frequency estimates could be made for pools of sizes 2-8, depending on the locus, by using an integer-valued estimator. Larger pools of size 12 and 24 tended to yield biased estimates; however, replicates of these estimates detected allele frequency differences among pools with different allelic compositions. We also derive an unbiased estimator of Hardy-Weinberg disequilibrium coefficients that uses multiple DNA pools and analyze the cost-efficiency of DNA pooling. DNA pooling yields the most potential cost savings when a large number of loci are employed using a large number of individuals, a situation becoming increasingly common as microsatellite loci are developed in increasing numbers of taxa.     

ExactFDR: exact computation of false discovery rate estimate in case-control association studies

Genome-wide association studies require accurate and fast statistical methods to identify relevant signals from the background noise generated by a huge number of simultaneously tested hypotheses. It is now commonly accepted that exact computations of association probability value (P-value) are preferred to chi(2) and permutation-based approximations. Following the same principle, the ExactFDR software package improves speed and accuracy of the permutation-based false discovery rate (FDR) estimation method by replacing the permutation-based estimation of the null distribution by the generalization of the algorithm used for computing individual exact P-values. It provides a quick and accurate non-conservative estimator of the proportion of false positives in a given selection of markers, and is therefore an efficient and pragmatic tool for the analysis of genome-wide association studies.     

Novel Jefferson salamander,Ambystoma jeffersonianum,microsatellite DNA markers detect population structure and hybrid complexes

Twenty polymorphic microsatellite DNA markers were isolated and characterized in Ambystoma jeffersonianum collected from three vernal pools in the mid‐Atlantic region of the U.S. These markers revealed a high degree of genetic diversity (7–23 alleles per locus), heterozygosity (46.7% to 100%), and allelic heterogeneity (96% of comparisons were statistically significant). Genetic distances were greatest in comparisons between collections, intermediate within collections, and least among sibling pairs. Six markers were trisomic in A. jeffersonianum‐A. laterale hybrids. These microsatellite DNA loci should allow delineation of genetic structure within and among populations of the diploid A. jeffersonianum and provide an effective method for identification of triploid hybrid individuals.     

Loss of genetic diversity as a signature of apricot domestication and diffusion into the Mediterranean Basin

ABSTRACT: BACKGROUND: Domestication generally implies a loss of diversity in crop species relative to their wild ancestors because of genetic drift through bottleneck effects. Compared to native Mediterranean fruit species like olive and grape, the loss of genetic diversity is expected to be more substantial for fruit species introduced into Mediterranean areas such as apricot (Prunus armeniaca L.), which was probably primarily domesticated in China. By comparing genetic diversity among regional apricot gene pools in several Mediterranean areas, we investigated the loss of genetic diversity associated with apricot selection and diffusion into the Mediterranean Basin. RESULTS: According to the geographic origin of apricots and using Bayesian clustering of genotypes, Mediterranean apricot (207 genotypes) was structured into three main gene pools: 'Irano-Caucasian', 'North Mediterranean Basin' and 'South Mediterranean Basin'. Among the 25 microsatellite markers used, only one displayed deviations from the frequencies expected under neutrality. Similar genetic diversity parameters were obtained within each of the three main clusters using both all SSR loci and only 24 SSR loci based on the assumption of neutrality. A significant loss of genetic diversity, as assessed by the allelic richness and private allelic richness, was revealed from the 'Irano-Caucasian' gene pool, considered as a secondary centre of diversification, to the northern and southwestern Mediterranean Basin. A substantial proportion of shared alleles was specifically detected when comparing gene pools from the 'North Mediterranean Basin' and 'South Mediterranean Basin' to the secondary centre of diversification. CONCLUSIONS: A marked domestication bottleneck was detected with microsatellite markers in the Mediterranean apricot material, depicting a global image of two diffusion routes from the 'Irano-Caucasian' gene pool: North Mediterranean and Southwest Mediterranean. This study generated genetic insight that will be useful for management of Mediterranean apricot germplasm as well as genetic selection programs related to adaptive traits.     

A genome-wide association scan for asthma in a general Australian population

To date, almost every chromosome has been implicated in genetic susceptibility to asthma to some degree. When compared with single nucleotide polymorphism, microsatellite markers exhibit high levels of heterozygosity and therefore provide higher statistical power in association. The objective of this study was to perform a genome-wide association study using 23,465 in-house microsatellite markers to detect asthma susceptibility regions in the Busselton population. In this study, three separate pooled DNA screenings yielded 18 markers with significantly different estimated frequencies in the three separate “case and control” pools: each pool consisting of 60 males and 60 females. These markers were evaluated by individual typing in 360 cases and 360 controls. Two markers showed significant differences between cases and controls (P = 0.001 and P = 0.003). Regions surrounding the two markers were subjected to high-density association mapping with a total of 14 additional markers. We were able to confirm and fine map the association in these two regions by typing 14 additional microsatellite markers (1805A09 (D18S0325i), P = 0.002; 1806D05 (D18S0181i), P = 0.001). Each region contains a predicted gene that showed strong associations with asthma. Further studies are underway to characterize the novel candidate asthma susceptibility genes identified in this genome-wide study.     

Detection of agriculturally important QTLs in chickens and analysis of the factors affecting genotyping strategy

Three single cross populations were generated in order to analyze factors affecting the ability to detect true linkage with minimum false positive or false negative associations, and to detect associations between markers and quantitative traits. The three populations are: (1) a broiler x broiler cross of a single sire and 34 dams, resulting in 266 progeny; (2) a broiler x broiler cross of a single sire and 41 dams resulting in 360 progeny; and (3) a broiler x layer cross of a single sire with 56 dams resulting in 1180 progeny. Based on these three resource populations we show that: a) gradient selective genotyping was more effective than the random selective genotyping; b) selective genotyping was significant at a selected proportion less than 62% of the cumulative truncation point; c) as few as 10% of selected individuals (5% of each of the two tails) were sufficient to show significant association between markers and phenotypes; d) a gradient slices approach was more powerful than using replicates of the extreme groups; and e) in resource populations resulting from crosses between lines of different backgrounds, most of the microsatellite markers used are polymorphic. We also used simulation to test factors affecting power to detect true associations between markers and traits that are hard to detect in experimental resource populations. Using defined populations in the simulation, we concluded that the following guidelines provide reliable detection of linked QTLs: 1) the resource population size should be larger than 100; 2) a QTL effect larger than 0.4 SD is detectable with a reasonable number of markers (>100) and resource population size (>200 subjects); 3) the DNA pool from each tail of the trait distribution should contain at least 10% of the resource family; 4) each of the two DNA pools should include more than 35 individuals. Some of these guidelines that were deduced from the simulation analysis have been confirmed in the experimental part of this study.     

Polymorphic microsatellite markers for Philippine tarsiers (<Emphasis Type="Italic">Tarsius syrichta</Emphasis>)

This study reports the isolation and characterization of eight microsatellite markers for the study of Philippine tarsiers (Tarsius syrichta), small primates endemic to this Southeast Asian archipelago. The markers were used to screen 14 Tarsius syrichta for allelic diversity. This suite of highly polymorphic microsatellites provides the first chance to genetically study parentage and dispersal patterns in Philippine tarsiers.     

Maximum-likelihood estimation of allelic dropout and false allele error rates from microsatellite genotypes in the absence of reference data

The importance of quantifying and accounting for stochastic genotyping errors when analyzing microsatellite data is increasingly being recognized. This awareness is motivating the development of data analysis methods that not only take errors into consideration but also recognize the difference between two distinct classes of error, allelic dropout and false alleles. Currently methods to estimate rates of allelic dropout and false alleles depend upon the availability of error-free reference genotypes or reliable pedigree data, which are often not available. We have developed a maximum-likelihood-based method for estimating these error rates from a single replication of a sample of genotypes. Simulations show it to be both accurate and robust to modest violations of its underlying assumptions. We have applied the method to estimating error rates in two microsatellite data sets. It is implemented in a computer program, Pedant, which estimates allelic dropout and false allele error rates with 95% confidence regions from microsatellite genotype data and performs power analysis. Pedant is freely available at http://www.stats.gla.ac.uk/ approximately paulj/pedant.html.     

Case-control association studies in mixed populations: correcting using genomic control

OBJECTIVE: Case-control association studies in mixed populations can result in spurious disease-marker associations if subpopulation disease prevalence and marker frequencies both differ. Genomic control (GC) uses neutral loci to correct for spurious association (due to population stratification), but how well this works remains undetermined. METHODS: We simulated and mixed populations with different disease and marker frequencies but without marker-disease association. We generated case-control datasets, calculated the chi2 for disease association with each marker, and applied two GC procedures, dividing by the mean chi2 or median-chi2/0.456. RESULTS: Corrections became conservative (false positive rate [FPR] <5%) with increasing subpopulation prevalence and marker differences. The mean correction resulted in FPRs close to 5% at average subpopulation allele frequency differences <0.26, but inclusion of just a few markers with large frequency differences resulted in conservative FPRs. FPRs from the median correction were mostly conservative but became anticonservative when a few markers with large frequency differences were included. CONCLUSION: GC can both lead to a notable loss of power to detect a true association (conservative) in many circumstances or may fail to eliminate the spurious associations (anticonservative). The mean correction factor is useful in certain situations to correct population stratification, but it is difficult to know when those situations exist.     

Isolation and characterization of polymorphic microsatellite markers in the gray,short‐tailed opossum (Monodelphis domestica)

Twenty‐six polymorphic microsatellite markers were isolated from (AC)_n and (AG)_n microsatellite‐enhanced genomic libraries of the gray, short‐tailed opossum Monodelphis domestica. All 26 loci showed high allelic diversity, with allele numbers ranging from five to 11 in a subset of 35 animals. Normal Mendelian inheritance was confirmed for 24 loci by analysing allelic segregation in 10, two‐generation, families. Non‐amplifying (null) alleles were detected at two loci, which we recommend be used only if pedigree data are available. We conclude that all of these microsatellite markers would be useful for quantitative trait locus mapping and population genetic studies.     

An Investigation on the Allelic Chi‐Square Test Used in Genetic Association Studies

Case‐control studies are primary study designs used in genetic association studies. Sasieni (Biometrics 1997, 53, 1253–1261) pointed out that the allelic chi‐square test used in genetic association studies is invalid when Hardy‐Weinberg equilibrium (HWE) is violated in a combined population. It is important to know how much type I error rate is deviated from the nominal level under violated HWE. We examine bounds of type I error rate of the allelic chi‐square test. We also investigate power of the goodness‐of‐fit test for HWE which can be used as a guideline for selecting an appropriate test between the allelic chi‐square test and the modified allelic chi‐square test, the latter of which was proposed for cases of violated HWE. In small samples, power is not large enough to detect the Wright's inbreeding model of small values of inbreeding coefficient. Therefore, when the null hypothesis of HWE is barely accepted, the modified test should be considered as an alternative method. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

A meta-analysis of kidney microarray datasets: investigation of cytokine gene detection and correlation with rt-PCR and detection thresholds

Background

Large-scale mutagenesis screens in the zebrafish employing the mutagen ENU have isolated several hundred mutant loci that represent putative developmental control genes. In order to realize the potential of such screens, systematic genetic mapping of the mutations is necessary. Here we report on a large-scale effort to map the mutations generated in mutagenesis screening at the Max Planck Institute for Developmental Biology by genome scanning with microsatellite markers.

Results

We have selected a set of microsatellite markers and developed methods and scoring criteria suitable for efficient, high-throughput genome scanning. We have used these methods to successfully obtain a rough map position for 319 mutant loci from the Tübingen I mutagenesis screen and subsequent screening of the mutant collection. For 277 of these the corresponding gene is not yet identified. Mapping was successful for 80 % of the tested loci. By comparing 21 mutation and gene positions of cloned mutations we have validated the correctness of our linkage group assignments and estimated the standard error of our map positions to be approximately 6 cM.

Conclusion

By obtaining rough map positions for over 300 zebrafish loci with developmental phenotypes, we have generated a dataset that will be useful not only for cloning of the affected genes, but also to suggest allelism of mutations with similar phenotypes that will be identified in future screens. Furthermore this work validates the usefulness of our methodology for rapid, systematic and inexpensive microsatellite mapping of zebrafish mutations.     

Genomewide association analysis of human narcolepsy and a new resistance gene

Human narcolepsy is a hypersomnia that is affected by multiple genetic and environmental factors. One genetic factor strongly associated with narcolepsy is the HLA-DRB1*1501-DQB1*0602 haplotype in the human leukocyte antigen region on chromosome 6, whereas the other genetic factors are not clear. To discover additional candidate regions for susceptibility or resistance to human narcolepsy, we performed a genomewide association study, using 23,244 microsatellite markers. Two rounds of screening with the use of pooled DNAs yielded 96 microsatellite markers (including 16 markers on chromosome 6) with significantly different estimated frequencies in case and control pools. Markers not located on chromosome 6 were evaluated by the individual typing of 95 cases and 95 controls; 30 markers still showed significant associations. A strong association was displayed by a marker on chromosome 21 (21q22.3). The surrounding region was subjected to high-density association mapping with 14 additional microsatellite markers and 74 SNPs. One microsatellite marker (D21S0012m) and two SNPs (rs13048981 and rs13046884) showed strong associations (P < .0005; odds ratios 0.19-0.33). These polymorphisms were in a strong linkage disequilibrium, and no other polymorphism in the region showed a stronger association with narcolepsy. The region contains three predicted genes--NLC1-A, NLC1-B, and NLC1-C--tentatively named "narcolepsy candidate-region 1 genes," and NLC1-A and NLC1-C were expressed in human hypothalamus. Reporter-gene assays showed that the marker D21S0012m in the promoter region and the SNP rs13046884 in the intron of NLC1-A significantly affected expression levels. Therefore, NLC1-A is considered to be a new resistance gene for human narcolepsy.     

Isolation and characterization of 12 microsatellite loci for population studies of Sulawesi tarsiers (Tarsius spp.)

This study reports the development and characterization of the first 12 microsatellite markers for tarsiers. Nine loci were isolated from Dian's tarsier, Tarsius dianae and three from the Philippine tarsier, Tarsius syrichta. The 12 markers were used to screen 40 individuals of Dian's tarsier and 40 individuals of the Lariang tarsier, Tarsius lariang for allelic diversity. This suite of highly polymorphic microsatellites provides the first chance to genetically study parentage patterns in tarsiers.     

Genetic Diversity and Population Structure of the Brachiaria brizantha Germplasm

Brachiaria brizantha (Hochst. ex A. Rich.) Stapf. (syn. Urochloa brizantha (Hochst. ex A. Rich.) R.D. Webster) is a species used primarily as forage in tropical America and Southeast Asia. B. brizantha has been extensively researched since the 1980s with the initiation of the Tropical Forages Breeding Program conducted by the Brazilian Agricultural Research Corporation (Empresa Brasileira de Pesquisa Agropecuária; EMBRAPA), holding one of the largest germplasm collections in the world. This work has identified 15 new microsatellite markers for this species, which have been used in addition to five previously reported markers, to estimate the genetic similarities among 172 accessions and six cultivars of this species. Similarity index values ranged from 0.40 to 1.00. Two duplications were found in the germplasm. A Bayesian analysis performed using the STRUCTURE 2.3.3 program revealed the presence of three clusters with different allelic pools. This analysis is valuable for the performance of crosses to explore heterosis; however, the mode of reproduction of the accessions and ploidy barriers must be observed for effective exploration. A grouping analysis using the neighbor-joining method was consistent with the STRUCTURE analysis, and a combination approach suggested that this germplasm collection does not exhibit considerable genetic variability despite the presence of three distinct allelic pools. The lack of correlation between the genetic and geographic distances is also discussed.     

Low gene flow but high genetic diversity in the threatened Mallorcan midwife toad Alytes muletensis

We investigated fine-scale genetic structuring in the rare and vulnerable Mallorcan midwife toad Alytes muletensis using eight polymorphic microsatellite markers. The current range of this amphibian is restricted to some 19 sites of which six are derived from reintroductions, all located in the mountain ranges of Mallorca. We sampled tadpoles from 14 pools covering 10 natural sites and two reintroduction sites for microsatellite DNA analyses. Relatively high levels of genetic variation were found in most pools (H(E) = 0.38-0.71, allelic richness = 2.6-6.2). Only at one pool has the population recently gone through a bottleneck. Dispersal between pools in different torrents does not occur whereas downstream dispersal between pools within the same torrent does happen at low frequencies. This occasional exchange of individuals does not lead to neighbouring pools in the same torrent being panmictic. This can be concluded because all F(ST) values (0.12-0.53) differ significantly from zero and STRUCTURE analyses identified neighbouring pools as separate populations. Furthermore, assignment and migration tests showed little exchange between neighbouring pools. If upstream locations or complete torrents go extinct, they are unlikely to be recolonized naturally. For conservation purposes, reintroductions of tadpoles to sites where local extinctions have occurred may therefore be advisable.     

First set of microsatellite markers for genetic characterization of the Eurasian beaver (Castor fiber) based on tissue and hair samples

Noninvasive genetic techniques have become indispensible tools in wildlife conservation and management. Here, we report the development of the first set of microsatellite markers for the Eurasian beaver (Castor fiber). All 15 loci show considerable variation within the sampled region in southwestern Germany, with number of alleles ranging from two to six alleles per locus. A comparison between tissue and hair samples revealed that amplification success was only slightly lower for hair samples, making their use in noninvasive monitoring feasible. Despite some evidence for false alleles and allelic dropout, 77% of all loci were genotyped successfully among all hair samples and loci tested. The developed markers will be used for subspecies differentiation and reconstruction of dispersal routes, following reintroductions in Central Europe.