Two-stage designs in case-control association analysis

DNA pooling is a cost-effective approach for collecting information on marker allele frequency in genetic studies. It is often suggested as a screening tool to identify a subset of candidate markers from a very large number of markers to be followed up by more accurate and informative individual genotyping. In this article, we investigate several statistical properties and design issues related to this two-stage design, including the selection of the candidate markers for second-stage analysis, statistical power of this design, and the probability that truly disease-associated markers are ranked among the top after second-stage analysis. We have derived analytical results on the proportion of markers to be selected for second-stage analysis. For example, to detect disease-associated markers with an allele frequency difference of 0.05 between the cases and controls through an initial sample of 1000 cases and 1000 controls, our results suggest that when the measurement errors are small (0.005), approximately 3% of the markers should be selected. For the statistical power to identify disease-associated markers, we find that the measurement errors associated with DNA pooling have little effect on its power. This is in contrast to the one-stage pooling scheme where measurement errors may have large effect on statistical power. As for the probability that the disease-associated markers are ranked among the top in the second stage, we show that there is a high probability that at least one disease-associated marker is ranked among the top when the allele frequency differences between the cases and controls are not <0.05 for reasonably large sample sizes, even though the errors associated with DNA pooling in the first stage are not small. Therefore, the two-stage design with DNA pooling as a screening tool offers an efficient strategy in genomewide association studies, even when the measurement errors associated with DNA pooling are nonnegligible. For any disease model, we find that all the statistical results essentially depend on the population allele frequency and the allele frequency differences between the cases and controls at the disease-associated markers. The general conclusions hold whether the second stage uses an entirely independent sample or includes both the samples used in the first stage and an independent set of samples.     

Asymptotic variances of QTL estimators with selective DNA pooling

Investigation on QTL-marker linkage usually requires a great number of observed recombinations, inferred from combined analysis of phenotypes and genotypes. To avoid costly individual genotyping, inferences on QTL position and effects can instead make use of marker allele frequencies. DNA pooling of selected samples makes allele frequency estimation feasible for studies involving large sample sizes. Linkage studies in outbred populations have traditionally exploited half-sib family designs; within the animal production context, half-sibships provide large families that are highly suitable for DNA pooling. Estimators for QTL position and effect have been proposed that make use of information from flanking markers. We present formulas derived by the delta method for the asymptotic variance of these estimators.     

Kinetic FP-TDI assay for SNP allele frequency determination

Strategies for identifying genetic risk factors in complex diseases by association studies require the comparison of allele frequencies of numerous SNPs between affected and control populations. Theoretically, hundreds of thousands of SNP markers across the genome will have to be genotyped in these studies. Genotyping SNPs one sample at a time is extremely costly and time consuming. To streamline whole genome association studies, some have proposed to screen SNPs by pooling the DNA samples initially for allele frequency determination and perform individual genotyping only when there is a significant discrepancy in allele frequencies between the affected and control populations. Here we describe a new method for determining the allele frequency of SNPs in pooled DNA samples using a two-color primer extension assay with real-time monitoring of fluorescence polarization (named kinetic FP-TDI assay). By comparing the ratio of the rate of incorporation of the two allele-specific dye-terminators, one can calculate the relative amounts of each allele in the pooled sample. The accuracy of allele frequency determination with pooled samples is within 3.3 +/- 0.8% of that determined by genotyping individual samples that make up the pool.     

Development of a Genotyping Microarray for Studying the Role of Gene-Environment Interactions in Risk for Lung Cancer

A microarray (LungCaGxE), based on Illumina BeadChip technology, was developed for high-resolution genotyping of genes that are candidates for involvement in environmentally driven aspects of lung cancer oncogenesis and/or tumor growth. The iterative array design process illustrates techniques for managing large panels of candidate genes and optimizing marker selection, aided by a new bioinformatics pipeline component, Tagger Batch Assistant. The LungCaGxE platform targets 298 genes and the proximal genetic regions in which they are located, using ∼13,000 DNA single nucleotide polymorphisms (SNPs), which include haplotype linkage markers with a minimum allele frequency of 1% and additional specifically targeted SNPs, for which published reports have indicated functional consequences or associations with lung cancer or other smoking-related diseases. The overall assay conversion rate was 98.9%; 99.0% of markers with a minimum Illumina design score of 0.6 successfully generated allele calls using genomic DNA from a study population of 1873 lung-cancer patients and controls.     

基于DNA池测序法筛选奶牛高信息量SNP标记的可行性

首先选择139个牛SNP标记, 利用DNA池测序法, 根据测序峰图中不同碱基信号峰高的比值确定了92个SNP为高信息量标记(比值>1/2); 为了进一步验证筛选的准确性, 对其中59个标记采用基质辅助激光解析电离飞行时间质谱(Matrix-assisted laser desorption/ionisation time-of-flight mass spectrometry, MALDI-TOF MS)技术检测了122头荷斯坦牛的基因型。结果显示, 检出率高于85%的标记有56个, 其平均最小等位基因频率(Minor allele frequency, MAF)为0.41, 最小值为0.27, 最大值为0.5; MAF>0.3的标记有54个, 占96.4%(54/56)。文章结果表明, 采用DNA池测序法筛选高信息量SNP标记是可行和可信的。     

Novel linkage mapping approach using DNA pooling in human and animal genetics. II. Detection of quantitative traits loci in dairy cattle

Selective DNA pooling is an advanced methodology for linkage mapping of quantitative trait loci (QTL) in farm animals. The principle is based on densitometric estimates of marker allele frequency in pooled DNA samples of phenotypically extreme individuals from half-sib, backcross and F(2) experimental designs in farm animals. This methodology provides a rapid and efficient analysis of a large number of individuals with short tandem repeat markers that are essential to detect QTL through the genome - wide searching approach. Several strategies involving whole genome scanning with a high statistical power have been developed for systematic search to detect the quantitative traits loci and linked loci of complex traits. In recent studies, greater success has been achieved in mapping several QTLs in Israel-Holstein cattle using selective DNA pooling. This paper outlines the currently emerged novel strategies of linkage mapping to identify QTL based on selective DNA pooling with more emphasis on its theoretical pre-requisite to detect linked QTLs, applications, a general theory for experimental half-sib designs, the power of statistics and its feasibility to identify genetic markers linked QTL in dairy cattle. The study reveals that the application of selective DNA pooling in dairy cattle can be best exploited in the genome-wide detection of linked loci with small and large QTL effects and applied to a moderately sized half-sib family of about 500 animals.     

Estimation of single nucleotide polymorphism allele frequency in DNA pools by using Pyrosequencing

Positional cloning of genes underlying complex diseases, such as type 2 diabetes mellitus (T2DM), typically follows a two-tiered process in which a chromosomal region is first identified by genome-wide linkage scanning, followed by association analyses using densely spaced single nucleotide polymorphic markers to identify the causal variant(s). The success of genome-wide single nucleotide polymorphism (SNP) detection has resulted in a vast number of potential markers available for use in the construction of such dense SNP maps. However, the cost of genotyping large numbers of SNPs in appropriately sized samples is nearly prohibitive. We have explored pooled DNA genotyping as a means of identifying differences in allele frequency between pools of individuals with T2DM and unaffected controls by using Pyrosequencing technology. We found that allele frequencies in pooled DNA were strongly correlated with those in individuals (r=0.99, P<0.0001) across a wide range of allele frequencies (0.02-0.50). We further investigated the sensitivity of this method to detect allele frequency differences between contrived pools, also over a wide range of allele frequencies. We found that Pyrosequencing was able to detect an allele frequency difference of less than 2% between pools, indicating that this method may be sensitive enough for use in association studies involving complex diseases where a small difference in allele frequency between cases and controls is expected.     

Localization of the gene for autosomal recessive congenital hereditary endothelial dystrophy (CHED2) to chromosome 20 by homozygosity mapping.

Congenital hereditary endothelial dystrophy (CHED) is a corneal disorder that presents with diffuse bilateral corneal clouding. Vision may be severely impaired, and many patients require corneal transplantation. Both autosomal dominant (AD) and autosomal recessive (AR) forms of the disorder have been described. The gene responsible for AD CHED (HGMW-approved symbol CHED1) has been mapped to the pericentromeric region of chromosome 20. Investigating a large, consanguineous Irish pedigree with autosomal recessive CHED, we have previously excluded linkage to this AD CHED locus. We now describe a genome-wide search using homozygosity mapping and DNA pooling. Evidence of linkage to chromosome 20p was demonstrated with a maximum lod score of 9.30 at a recombination fraction of 0.0 using microsatellite marker D20S482. A region of homozygosity in all affected individuals was identified, narrowing the disease gene locus to an 8-cM region flanked by markers D20S113 and D20S882. This AR CHED (HGMW-approved symbol CHED2) disease gene locus is physically and genetically distinct from the AD CHED locus.     

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.     

Applying polynomial standard curve method to correct bias encountered in estimating allele frequencies using DNA pooling strategy

DNA pooling approach is a cost-saving strategy which is crucial for multiple-SNP association study and particularly for laboratories with limited budget. However, the biased allele frequency estimates cannot be completely abolished by κ correction. Using the SNaPshot™, we systematically examined the relations between actual minor allele frequencies (AMiAFs) levels and estimates obtained from the pooling process for all six types of SNPs. We applied principle of polynomial standard curves method (PSCM) to produce allele frequency estimates in pooled DNA samples and compared it with the κ method. The results showed that estimates derived from the PSCM were in general closer to AMiAFs than those from the κ method, particularly for C/G and G/T polymorphisms at the range of AMiAF between 20–40%. We demonstrated that applying PSCM in the SNaPshot™ platform is suitable for multiple-SNP association study using pooling strategy, due to its cost effectiveness and estimation accuracy.     

Allele frequency of two intragenic microsatellite loci of SEL1L gene in Northern Italian population

Two cytosine-adenine (CA) repeats CAR/CAL and RepIN20 occur in the human SEL1L gene, which is regarded as a candidate gene for insulin-dependent diabetes mellitus (IDDM) and Grave's disease. We have characterized these repeats to determine if they might serve as effective microsatellite markers for linkage analysis to clarify whether SEL1L gene plays a role in the pathogenesis of these autoimmune diseases. The allele frequencies and average heterozygosity of the microsatellite repeats were analysed in 94 DNA samples from peripheral blood mononuclear (PBMC) cells from adults of Northern Italy. The average heterozygosity was 0.68 for CAR/CAL polymorphism and 0.85 for RepIN20. The size of PCR fragments of CAR/CAL ranged from 207–225 bp and the most frequent allele was 207 bp (40.4%). The size of the fragments of RepIN20 ranged from 237–255 bp and the most frequent allele was 249 bp (30.8%). In the light of the highly polymorphic nature of both microsatellites and their intragenic location in SEL1L gene, we suggest that they could provide a means for linkage analysis to clarify the potential role of SEL1L in conferring susceptibility to IDDM or Grave's disease.     

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.     

Genotyping pooled DNA using 100K SNP microarrays: a step towards genomewide association scans

The identification of quantitative trait loci (QTLs) of small effect size that underlie complex traits poses a particular challenge for geneticists due to the large sample sizes and large numbers of genetic markers required for genomewide association scans. An efficient solution for screening purposes is to combine single nucleotide polymorphism (SNP) microarrays and DNA pooling (SNP-MaP), an approach that has been shown to be valid, reliable and accurate in deriving relative allele frequency estimates from pooled DNA for groups such as cases and controls for 10K SNP microarrays. However, in order to conduct a genomewide association study many more SNP markers are needed. To this end, we assessed the validity and reliability of the SNP-MaP method using Affymetrix GeneChip^® Mapping 100K Array set. Interpretable results emerged for 95% of the SNPs (nearly 110000 SNPs). We found that SNP-MaP allele frequency estimates correlated 0.939 with allele frequencies for 97605 SNPs that were genotyped individually in an independent population; the correlation was 0.971 for 26 SNPs that were genotyped individually for the 1028 individuals used to construct the DNA pools. We conclude that extending the SNP-MaP method to the Affymetrix GeneChip^® Mapping 100K Array set provides a useful screen of >100000 SNP markers for QTL association scans.     

Evidence for association of D1S249 locus on human chromosome 1 with the susceptibility to essential hypertension in Han Chinese

Essential hypertension (EH) is thought to result from theinteraction of environmental and genetic factors. The molecular genetics of EH has witnessed considerable progress during the past few years. However, the number of genes involved, their chromosomal location and the magnitude of their effect on EH susceptibility are unknown. We conducted the present study to screen susceptibility genes to essential hypertension using a genome-wide scanning method in a group of Han people from Fangshan district located in the southwest of Beijing. A case-control study and affected sibpair were performed. Genotyping was carried out using a fluorescence-based semiautomated technique on automated DNA sequencer (ABI 377, PE). The basis for the genome-screen was the ABI prism linkage mapping sets of 400 microsatellite markers (version 2, PE, Co.). PCR for amplification of markers was carried out as multiplex reactions with Ampli Taq gold (PE, Co.) following protocols developed in our laboratory. Data were exported as a text file from genotyper for subsequent two-point affected sibpair linkage analysis. The data from case-control association study showed a linkage disequilibrium between EH and marker D1S249 locus (X2 = 14.6, P = 0.002). There are 12 alleles in the D1S249 locus. The frequency of A9 allele in hypertension was higher than in normotensives, (13.6% v.s. 2.7%, X2 = 6.30, p = 0.01, OR = 4.57, 95%CI = 1.24-25.4). The data from two-point affected sibpair linkage analysis demonstrated a linkage between EH and A9 allele, P＜0.05. It suggested that microsatellite marker D1S249 locus might be associated with the genetic susceptibility to essential hypertension in Han Chinese.     

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.     

An autosomal recessive nonsyndromic-hearing-loss locus identified by DNA pooling using two inbred Bedouin kindreds.

Autosomal recessive nonsyndromic hearing loss (ARNSHL) is the most common form of severe inherited childhood deafness. We present the linkage analysis of two inbred Bedouin kindreds from Israel that are affected with ARNSHL. A rapid genomewide screen for markers linked to the disease was performed by using pooled DNA samples. This screen revealed evidence for linkage with markers D9S922 and D9S301 on chromosome 9q. Genotyping of individuals from both kindreds confirmed linkage to chromosome 9q and a maximum combined LOD score of 26.2 (recombination fraction [theta] .025) with marker D9S927. The disease locus was mapped to a 1.6-cM region of chromosome 9ql3-q2l, between markers D9S15 and D9S927. The disease segregates with a common haplotype in the two kindreds, at markers D9S927, D9S175, and D9S284 in the linked interval, supporting the hypothesis that both kindreds inherited the deafness gene from a common ancestor. Although this nonsyndromic-hearing-loss (NSHL) locus maps to the same cytogenetic interval as DFNB7, it does not overlap the currently defined DFNB7 interval and may represent (1) a novel form of NSHL in close proximity to DFNB7 or (2) a relocalization of the DFNB7 interval to a region telomeric to its reported location. This study further demonstrates that DNA pooling is an effective means of quickly identifying regions of linkage in inbred families with heterogeneous autosomal recessive disorders.     

SNP haplotype tagging from DNA pools of two individuals

Background  

DNA pooling is a technique to reduce genotyping effort while incurring only minor losses in accuracy of allele frequency estimates for single nucleotide polymorphism (SNP) markers.     

Identification of a RAPD marker for palmitic-acid concentration in the seed oil of spring turnip rape (Brassica rapa ssp. oleifera)

F₂ progeny (105 individuals) from the cross Jo4002 x Sv3402 were used to identify DNA markers associated with palmitic-acid content in spring turnip rape (Brassica rapa ssp. oleifera). QTL mapping and ANOVA analysis of 140 markers exposed one linkage group with a locus controlling palmitic-acid content (LOD score 27), and one RAPD (random amplified polymorphic DNA) marker, OPB-11a, closely linked (1.4 cM) to this locus. Palmitic-acid content in the 62 F₂ plants with the visible allele of marker OPB-11a was 8.45 ±3.15%, while that in the 24 plants without it was 4.59 ±0.97%. As oleic-acid concentration is affected by a locus on the same linkage group as the palmitic-acid locus, this locus probably controls the chain elongation from palmitic acid to oleic acid (through stearic acid). Marker OPB-11a may be used in future breeding programs of spring turnip rape to simplify and hasten the selection for palmitic-acid content.     

Metalloproteinase’s Activity and Oxidative Stress in Mild Cognitive Impairment and Alzheimer’s Disease

Matrix metalloproteinases (MMPs) and oxidative stress have been implicated in neurological diseases such as Alzheimer’s disease (AD). Plasma MMP-2 and MMP-9 activities were assessed in Mild Cognitive Impairment (MCI) and AD subjects compared with aged-matched controls, and subsequently analysed in relation to oxidative stress markers. Both MMP-2 and MMP-9 showed no significant changes versus control subjects. Plasma glutathione peroxidase Se-dependent (GPx-Se) activity and malondialdehyde (MDA) levels were higher in AD than in controls (P < 0.05), suggesting a role for GPx-Se in controlling oxidative stress in AD. Negative correlations were observed between MMPs and MDA in AD and MCI patients (P < 0.05). In conclusion, oxidative stress events did not include activation of MMPs and this similar pattern in AD and MCI suggests that both are biochemically equivalent.     

DNA 池结合DHPLC 和直接测序技术在江豚SNPs 检测中的应用

选取江豚基因组中的2 个已知单核苷酸多态性(single nucleotide polymorphisms,SNPs)位点,通过PCR 扩增,将PCR 产物按基因频率不同制备成0 ～ 50% 的11 个DNA 池(DNA pool),用于变性高效液相色谱(denaturing high performance liquid chromatography,DHPLC)和直接测序分析,以探讨DNA 池中基因频率的最低要求。结果显示,当稀有等位基因的基因频率不少于5% 时可在DHPLC 检测过程中明显分辨;而利用DNA 池进行直接测序时的基因频率则需达到10% 。这提示,为保证DHPLC 分析的准确性和可靠性,制备DNA 池时等摩尔DNA 混合的个体数最好不超过10 个。DNA 池结合DHPLC 技术的高效性与准确性可在大规模的SNPs 位点筛选中发挥作用。