High-throughput genotyping of a common deletion polymorphism disrupting the TRY6 gene and its association with breast cancer risk
© Wagner et al; licensee BioMed Central Ltd. 2007
Received: 03 April 2007
Accepted: 29 June 2007
Published: 29 June 2007
Copy number polymorphisms caused by genomic rearrangements like deletions, make a significant contribution to the genomic differences between two individuals and may add to disease predisposition. Therefore, genotyping of such deletion polymorphisms in case-control studies could give important insights into risk associations.
We mapped the breakpoints and developed a fluorescent fragment analysis for a deletion disrupting the TRY6 gene to exemplify a quick and cheap genotyping approach for such structural variants. We showed that the deletion is larger than predicted and encompasses also the pseudogene TRY5. We performed a case-control study to test an association of the TRY6 deletion polymorphism with breast cancer using a single nucleotide polymorphism which is in 100% linkage disequilibrium with the deletion. We did not observe an effect of the deletion on breast cancer risk (OR 1.05, 95% CI 0.71–1.56).
Although we did not observe an association between the TRY6 deletion polymorphism and breast cancer risk, the identification and investigation of further deletions using the present approach may help to elucidate their effect on disease susceptibility.
An unexpectedly large amount of structural variation in the human genome has recently become apparent . The identification of genomic deletions, insertions and inversions, ranging from 1 kb to several Mb in size, has broadened our knowledge about genetic variation between individuals, which may lead to new insights into phenotypic variation and clinical outcome . The rearrangements can encompass exons or even entire genes and their regulatory regions. Such structural variants, called copy number polymorphisms (CNPs), have been identified by array-based methods [3–6]. Interestingly, a high percentage of copy number variants has been found to be in close proximity to segmental duplications, suggesting a recombination mechanism between the repetitive sequences [4, 7, 8]. Owing to the high sequence similarity between duplications and their consequent high recombination potential, segmental duplications are frequently found at breakpoints of both disease-associated and evolutionary rearrangements . The HapMap data on family trios have been used to identify CNPs [10, 11]. As deletions normally are not detected by standard single nucleotide polymorphism (SNP) genotyping methods, hemizygotes are usually misclassified as homozygotes. When a deletion is transmitted from parent to offspring, the child will show a null genotype or a genotype violating the rules of mendelian inheritance . However, the exact borders of the structural variants still remain to be ascertained, which makes it difficult to perform high-throughput genotyping of CNPs. Furthermore, linkage disequilibrium (LD) between deletions and SNPs has been demonstrated [11, 12], suggesting that these variants show a similar evolutionary history . The use of SNP genotyping may also be an alternative for the genotyping of the deletions. However, recent reports have raised concern over the use of SNPs as surrogates for CNPs in structurally dynamic regions [7, 13].
We present an approach to specify the breakpoints of the common deletion polymorphism encompassing the TRY6 gene and then describe a genotyping assay, using fluorescent multiplex PCR, to investigate whether the deletion is associated with breast cancer risk. Alternatively, the deletion can be genotyped by analysing a neighbouring SNP which is in 100% LD with the deletion polymorphism. These methods can be adapted for other deletion polymorphisms.
Identification of TRY6 breakpoints
Fragment analysis for a case-control study
rs13230029 is not associated with breast cancer risk
Primer sequences used to map the breakpoints.
Product size (bp)
F: GAC ACA AAG ACG TGG GAG TC
R: GAG GCA GCC TGG CTG GGA
F: GTG AAA GAG GCT GGG AAG GTG
R: CCC TTC TTT CAC AGC TGG GGC
F: ACA GGT GAT AAA AGC CCG AGC C
R: GCC CTC AGA TGA TTA TTG GTG AAG T
4F † /R
F: TTG AGA GAT CTA GAA GCC CCA C
R: GCC CTC AGA TGA TTA TTG GTG AAG T
F: TAA GTC TCC TTT TAG ATG CCA CC
R: GGA GTT TTC ATG TAA GCA GCA GTG
F: CTC CAG AGC TAT AAA GAC GGG C
R: CGG CAG GGC ATA TGT CTG CCA
F: AAA AGA GAG AAG CAT TCA GTG GG
R: CGG TTG CTT CCT GCT AAT TAG AA
F: GAA TGA ACA GTT CAT CTA TGC G
R: CGG TGT GCT TCG TTC TGG AAA T
F: TCT GAG GGC TGT GAC ACC ATG
R: GAC ACA GGT GAC ATG AAG CCT G
F: CTG GGA AAG GAT CCC TCA AG
R: GTC AAC TGT GGC TGC CAG TG
F: ACA TGT TCA GGG ACA ACA CAG C
R: CAC ACC TCT CTG CTC ATG AAT AA
F: AAC TCT GAC ATG TGA TCA GGG G
R: CAT GTG CAT CCT GTC ATA GGT TC
F: AAG CAG CCA CAG GCT GGG AGC
R: CCT GGC TGG GAG TTT TGC AGT
F: GAT GAA TTG GGG TAT CAG GAA TG
R: GCT CAC TAA TCC ACT GTG TTT CAT
F: CTG TGA TTG TTT AAG GAA GAG CG
R: CAC AAC TCT CTG CTC ACA GAT AA
F: CCC AAG TGT TTT TGT GTG GCG AAT GG
R: CAG CTA CCC ATG CTC AAT TTG G
F: GGA AGG ACA GTG GTT TGG CAT G
R: GTC AAC TGT GGC TGC CAG TG
Allele and genotype frequencies, odds ratios, 95% confidence intervals and p-values for the SNP rs13230029.
OR (95% CI)
Large genomic variations show a potential at least as great as SNPs to have an impact on disease susceptibility, given that regulatory regions, exons and even whole genes can be deleted, duplicated or disrupted . Deletions result from non-allelic homologous recombination events that occur between blocks of duplicated sequences [4, 7]. Recent studies have provided a preliminary characterisation of deletions and paved the way to disease association studies [3–6, 10–12]. However, fast, cheap and reliable techniques are needed to perform case-control studies for the identification of disease correlation.
Real-time PCR can be used for the detection of individual deletions or duplications but it is not appropriate for multiplexing . Furthermore, the optimisation process takes time and the throughput is limited. For the amplification of multiple regions, quantitative multiplex PCR of short fluorescent fragments (QMPSF), multiplex amplifiable probe hybridisation (MAPH) and multiplex ligation-dependent probe amplification (MLPA) are more adequate . However, MAPH and MLPA require several handling steps, the commercially available kits are still expensive and pre-designed for only common sites and these methods are not suitable for the identification of heterozygous genotypes. QMPSF is a faster method, yet the optimisation, especially for heterozygotes, is time-consuming and tedious and one should have control samples with known genotypes for comparison of the results. All the above described methods could be useful for identifying breakpoints in homozygous samples where it is only important to know if the product is present or absent. A quicker and safer approach is a triplex PCR using primers that flank the deletion break points, as described here and elsewhere [3, 17, 18]. If the primers are labelled fluorescently, an automatic analysis ensures a fast evaluation of the results. However, it is essential to know the deletion borders to design the primers.
We refined the borders of the deletion by amplifying short fragments in intervals of 1500–2000 bp throughout the predicted and the extended region of the TRY6 gene deletion polymorphism in individuals with known homozygous genotypes obtained from the HapMap subjects . Homozygote carriers showed the absence of the PCR products, while the fragments were present in the wild type and the heterozygote samples. As an internal control for the sample quality, we simultaneously amplified a fragment of an independent gene. We mapped the breakpoints by sequencing the complete region between the two fragments flanking the deletion in the homozygous samples and by comparing the obtained sequence to the published one (NT_007914.14, NCBI Build 36). A triplex PCR with primers neighbouring the borders enabled discrimination between the three genotypes. Due to the fluorescence of the reverse primer, the products could be visualised on a sequencer and analysed automatically with the appropriate software. This method is fast, cheap (0.70€/genotype) and easy to perform.
TRY 6 is deleted in a common polymorphism [11, 14]. Although this locus appears to encode a protein similar to trypsinogen, the locus is thought to be a transcribed pseudogene . ESTs support its transcription, but expression of its predicted protein has not been observed. The predicted protein sequence of TRY 6 differs significantly from the known functional trypsinogens, including a different amino acid at the conserved residue 122 which is important for autolysis. However, it has been suggested that it is expressed in minute amounts in the thymus . A recent study has identified TRY6 to be up-regulated in metastatic non-small cell lung cancer leading to an enhanced cell migration . Patients with high expression of TRY6 also had a substantially worse prognosis than patients without TRY6 expression. We picked this deletion polymorphism to test our method because it was shown to be real, frequent and affecting the coding region . We did not find an association of the TRY6 deletion with breast cancer susceptibility. With our sample size, considering that the cases had familial or early onset breast cancer and that ~17% of the controls were homozygous for the deletion, we had an >80% power to detect an OR of 0.6. Thus, we can only exclude a major to moderate effect of the TRY6 deletion on breast cancer risk. The negative findings of our study may be explained by the fact that TRY6, and the deleted neighbouring gene, TRY5, are pseudogenes. Although TRY6 mRNA expression has been established, the function of the gene remains to be clarified.
Several reports have shown that deletion polymorphisms show strong LD with common SNPs [11, 12, 18]. We genotyped a SNP and confirmed its LD with the deletion. Genotyping a SNP is faster and cheaper than detecting the structural variant. However, recent studies have found only modest evidence of LD between CNPs and HapMap SNPs, suggesting that duplication rich regions are not suitable for the identification of linked SNPs due to technical restrictions [7, 13]. This raises the concern that several genomic regions most prone to rearrangements have inadequate SNP density to successfully map CNPs with the use of available SNP markers. Thus, additional methods, such as the fluorescent fragment analysis for deletion polymorphisms described by us, are needed to reliably genotype CNPs.
The comparability of the results generated from various array-based studies is currently the major obstacle in the field, contributing to confusion in data interpretation between different reports. Because accurate genotyping requires exact sequence data, refining the deletion breakpoints should rank first. Furthermore, the mapping of CNP breakpoints will also enable the identification of nearby SNPs, which would be in LD with the deletion. Thus, the description of all variants will help to design studies relating to disease, pharmacogenomics and clinical practice. They can be searched in the Database of Genomics Variants , which is regularly updated with the latest detected CNPs. The development of methods to identify other subtle variations such as inversions or insertions still remains a challenge.
We presented a simple high-throughput approach to genotype common deletion polymorphisms. Although we did not observe an association between the TRY6 deletion polymorphism and breast cancer risk, the identification and investigation of further deletions using the present approach may help to elucidate their effect on disease susceptibility.
Based on the analysis carried out on the Hapmap samples by McCarroll et al. [11, 21], we selected CEPH individuals homozygous for the TRY6 deletion polymorphism to define the breakpoints of the deletion. We amplified several fragments throughout the TRY6 gene region (Figure 1B). The fragments were first limited to the predicted region but then extended outwards on both directions. Homozygotes for the deletion should not give a PCR product within the deleted region, which was used in mapping the breakpoint. The primer sequences and product sizes are listed in Table 1. Amplification was performed with 5 ng genomic DNA in a 10 μl reaction volume using 1× PCR buffer, 1.5 mM MgCl2, 0.11 μM dNTP Mixture (Invitrogen, Paisley, UK), 0.15 μM of each primer (Thermo Electron Corporation, Ulm, Germany) and 0.3 U PlatinumTaq Polymerase (Invitrogen). The PCR was carried out in a GeneAmp 9700 PCR system (Applied Biosystems, Foster City, USA) using 94°C for 2 min, followed by 3 cycles of 94°C for 1 min, 67°C for 1 min and 72°C for 1 min and 32 cycles of 94°C for 30 s, 66°C for 30 s and 72°C for 30 s. The final extension was performed for 6 min at 72°C. The PCR products were visualised on 1.5% agarose gels stained with ethidium bromide.
The PCR product was cleaned-up using 0.75 μl ExoSapIT (USB Amersham, Uppsala, Sweden) for 40 min at 37°C followed by 15 min at 85°C. The sequencing reaction was carried out using the BigDye Terminator Cycle Sequencing Ready Reaction Kit (Applied Biosystems) on the Genetic Analyzer 3130XL as described earlier The original data was analysed by the Sequencing Analysis 5.2 and DNASTAR Lasergene 5.0 (DNASTAR Inc., Madison, USA) softwares.
The sequence of a sample with a homozygous deletion was aligned to the sequence surrounding the observed breakpoint on both sides. The reference sequence was derived from Genbank NT_007914.14, NCBI build 36. The alignment (Figure 2) was performed using MegAlign 5.05 (DNAStar Inc., USA).
The analysis of the TRY6 deletion was performed on genomic DNA from 397 Polish familial and early age breast cancer cases (mean age 46 years, range 26–81) and 454 regionally and ethnically matched unrelated female controls (mean age 42 years, range 16–79). The inclusion criteria for the cases were (i) at least 2 first-degree relatives with breast and/or ovarian cancer regardless of age, (ii) breast cancer diagnosed below the age of 35 without family history, (iii) bilateral breast cancer regardless of the family history, (iv) breast and ovarian cancer diagnosed in one patient regardless of the family history and (v) breast cancer diagnosed below 50 years of age regardless of family history The subjects corresponding to criteria (i)-(iv), 215 cases, were collected during the years 1997 to 2002 by the Chemotherapy Clinics and the Genetic Counselling Service (Gliwice, Poland) and the subjects corresponding to criteria (v), 182 cases, were collected during the time period December 2002 to March 2004 by the Surgery Clinics (Gliwice, Poland). The cases were tested for four founder mutations in BRCA1 and two in BRCA2, which account for ~90% of the mutations in the Polish population, and found to be negative . We included familial and early onset cases to our study because it has been shown that selection of cases based on the family history increases the power to detect low-penetrance variants [25, 26]. The study was approved by the ethical committee of the University of Heidelberg.
Fluorescent fragment analysis
In order to genotype individuals for the deletion polymorphism (Figure 3), we used a triplex PCR with the primers 4F, 4aR and 10aR (Table 1). The reverse primers were labelled with 6-FAM (Invitrogen). We performed the amplification with the conditions mentioned above. The resulting product was digested with 4 U Sch I (Fermentas, St. Leon-Roth, Germany) for 4 h at 37°C. 1 μl of the digestion product was added to 10 μl HIDI-Formamide/GeneScan ROX350 standard size marker (mixed according to the manufacturer's instruction, Applied Biosystems). The mixture was denatured for 5 min at 95°C, loaded onto the ABI PRISM 3130XL Genetic Analyser and analysed by the GeneMapper software version 4.0 (Applied Biosystems).
The TaqMan assay for SNP rs13230029 was obtained as custom assay from Applied Biosystems and the assay details are available on request. The reaction was performed in 5 μl using 5 ng of genomic DNA, 1× TaqMan Universal Master Mix (Applied Biosystems) and 0.6× Assay-Mix (40×) per reaction. PCR was performed at 50°C for 2 min, 95°C for 10 min and 35–45 cycles at 92°C for 15 s and 60°C for 1 min. PCR was performed in a GeneAmp PCR System 9700 thermocycler and the number of cycles was dependent on the genotype clustering. The samples were read and analysed in an ABI Prism 7900 HT sequence detection system using the SDS 1.2 software (Applied Biosystems).
The observed genotype frequencies in the breast cancer cases and controls were tested for Hardy-Weinberg equilibrium (HWE) and the difference between the observed and expected frequencies was tested for significance using the χ2-test. Odds ratios (ORs) and 95% confidence intervals (95% CIs) were calculated for associations between genotypes and breast cancer. The calculations were carried out using the HWE test tool offered by the Institute of Human Genetics, TU Munich . Power calculation was carried out using the PS software for power and sample size calculation [28, 29].
This study was supported by the grants from State Committee for Scientific Research (PBZ-KBN-090 P05/02 to EG) and a grant from EU (LSHC-CT-2004-503465).
- Feuk L, Marshall CR, Wintle RF, Scherer SW: Structural variants: changing the landscape of chromosomes and design of disease studies. Hum Mol Genet. 2006, 15 Spec No 1: R57-66. 10.1093/hmg/ddl057.View ArticlePubMedGoogle Scholar
- Feuk L, Carson AR, Scherer SW: Structural variation in the human genome. Nat Rev Genet. 2006, 7 (2): 85-97. 10.1038/nrg1767.View ArticlePubMedGoogle Scholar
- Tuzun E, Sharp AJ, Bailey JA, Kaul R, Morrison VA, Pertz LM, Haugen E, Hayden H, Albertson D, Pinkel D, Olson MV, Eichler EE: Fine-scale structural variation of the human genome. Nat Genet. 2005, 37 (7): 727-732. 10.1038/ng1562.View ArticlePubMedGoogle Scholar
- Sharp AJ, Locke DP, McGrath SD, Cheng Z, Bailey JA, Vallente RU, Pertz LM, Clark RA, Schwartz S, Segraves R, Oseroff VV, Albertson DG, Pinkel D, Eichler EE: Segmental duplications and copy-number variation in the human genome. Am J Hum Genet. 2005, 77 (1): 78-88. 10.1086/431652.PubMed CentralView ArticlePubMedGoogle Scholar
- Iafrate AJ, Feuk L, Rivera MN, Listewnik ML, Donahoe PK, Qi Y, Scherer SW, Lee C: Detection of large-scale variation in the human genome. Nat Genet. 2004, 36 (9): 949-951. 10.1038/ng1416.View ArticlePubMedGoogle Scholar
- Sebat J, Lakshmi B, Troge J, Alexander J, Young J, Lundin P, Maner S, Massa H, Walker M, Chi M, Navin N, Lucito R, Healy J, Hicks J, Ye K, Reiner A, Gilliam TC, Trask B, Patterson N, Zetterberg A, Wigler M: Large-scale copy number polymorphism in the human genome. Science. 2004, 305 (5683): 525-528. 10.1126/science.1098918.View ArticlePubMedGoogle Scholar
- Locke DP, Sharp AJ, McCarroll SA, McGrath SD, Newman TL, Cheng Z, Schwartz S, Albertson DG, Pinkel D, Altshuler DM, Eichler EE: Linkage disequilibrium and heritability of copy-number polymorphisms within duplicated regions of the human genome. Am J Hum Genet. 2006, 79 (2): 275-290. 10.1086/505653.PubMed CentralView ArticlePubMedGoogle Scholar
- Goidts V, Cooper DN, Armengol L, Schempp W, Conroy J, Estivill X, Nowak N, Hameister H, Kehrer-Sawatzki H: Complex patterns of copy number variation at sites of segmental duplications: an important category of structural variation in the human genome. Hum Genet. 2006, 120 (2): 270-284. 10.1007/s00439-006-0217-y.View ArticlePubMedGoogle Scholar
- Samonte RV, Eichler EE: Segmental duplications and the evolution of the primate genome. Nat Rev Genet. 2002, 3 (1): 65-72. 10.1038/nrg705.View ArticlePubMedGoogle Scholar
- Conrad DF, Andrews TD, Carter NP, Hurles ME, Pritchard JK: A high-resolution survey of deletion polymorphism in the human genome. Nat Genet. 2006, 38 (1): 75-81. 10.1038/ng1697.View ArticlePubMedGoogle Scholar
- McCarroll SA, Hadnott TN, Perry GH, Sabeti PC, Zody MC, Barrett JC, Dallaire S, Gabriel SB, Lee C, Daly MJ, Altshuler DM: Common deletion polymorphisms in the human genome. Nat Genet. 2006, 38 (1): 86-92. 10.1038/ng1696.View ArticlePubMedGoogle Scholar
- Hinds DA, Kloek AP, Jen M, Chen X, Frazer KA: Common deletions and SNPs are in linkage disequilibrium in the human genome. Nat Genet. 2006, 38 (1): 82-85. 10.1038/ng1695.View ArticlePubMedGoogle Scholar
- Redon R, Ishikawa S, Fitch KR, Feuk L, Perry GH, Andrews TD, Fiegler H, Shapero MH, Carson AR, Chen W, Cho EK, Dallaire S, Freeman JL, Gonzalez JR, Gratacos M, Huang J, Kalaitzopoulos D, Komura D, MacDonald JR, Marshall CR, Mei R, Montgomery L, Nishimura K, Okamura K, Shen F, Somerville MJ, Tchinda J, Valsesia A, Woodwark C, Yang F, Zhang J, Zerjal T, Zhang J, Armengol L, Conrad DF, Estivill X, Tyler-Smith C, Carter NP, Aburatani H, Lee C, Jones KW, Scherer SW, Hurles ME: Global variation in copy number in the human genome. Nature. 2006, 444 (7118): 444-454. 10.1038/nature05329.PubMed CentralView ArticlePubMedGoogle Scholar
- Rowen L, Koop BF, Hood L: The complete 685-kilobase DNA sequence of the human beta T cell receptor locus. Science. 1996, 272 (5269): 1755-1762. 10.1126/science.272.5269.1755.View ArticlePubMedGoogle Scholar
- Yamamoto H, Iku S, Adachi Y, Imsumran A, Taniguchi H, Nosho K, Min Y, Horiuchi S, Yoshida M, Itoh F, Imai K: Association of trypsin expression with tumour progression and matrilysin expression in human colorectal cancer. J Pathol. 2003, 199 (2): 176-184. 10.1002/path.1277.View ArticlePubMedGoogle Scholar
- Diederichs S, Bulk E, Steffen B, Ji P, Tickenbrock L, Lang K, Zanker KS, Metzger R, Schneider PM, Gerke V, Thomas M, Berdel WE, Serve H, Muller-Tidow C: S100 family members and trypsinogens are predictors of distant metastasis and survival in early-stage non-small cell lung cancer. Cancer Res. 2004, 64 (16): 5564-5569. 10.1158/0008-5472.CAN-04-2004.View ArticlePubMedGoogle Scholar
- Newman TL, Rieder MJ, Morrison VA, Sharp AJ, Smith JD, Sprague LJ, Kaul R, Carlson CS, Olson MV, Nickerson DA, Eichler EE: High-throughput genotyping of intermediate-size structural variation. Hum Mol Genet. 2006, 15 (7): 1159-1167. 10.1093/hmg/ddl031.View ArticlePubMedGoogle Scholar
- Wagner K, Hemminki K, Grzybowska E, Bermejo JL, Butkiewicz D, Pamula J, Pekala W, Forsti A: Polymorphisms in the growth hormone receptor: A case-control study in breast cancer. x. 2006, 118 (11): 2903-2906.Google Scholar
- Chen JM, Montier T, Ferec C: Molecular pathology and evolutionary and physiological implications of pancreatitis-associated cationic trypsinogen mutations. Hum Genet. 2001, 109 (3): 245-252. 10.1007/s004390100580.View ArticlePubMedGoogle Scholar
- Database of Genomics Variants. [http://projects.tcag.ca/variation/]
- A haplotype map of the human genome. Nature. 2005, 437 (7063): 1299-1320. 10.1038/nature04226.Google Scholar
- Wagner K, Hemminki K, Grzybowska E, Klaes R, Burwinkel B, Bugert P, Schmutzler RK, Wappenschmidt B, Butkiewicz D, Pamula J, Pekala W, Forsti A: Polymorphisms in genes involved in GH1 release and their association with breast cancer risk. Carcinogenesis. 2006, 27 (9): 1867-1875. 10.1093/carcin/bgl036.View ArticlePubMedGoogle Scholar
- Forsti A, Jin Q, Grzybowska E, Soderberg M, Zientek H, Sieminska M, Rogozinska-Szczepka J, Chmielik E, Utracka-Hutka B, Hemminki K: Sex hormone-binding globulin polymorphisms in familial and sporadic breast cancer. Carcinogenesis. 2002, 23 (8): 1315-1320. 10.1093/carcin/23.8.1315.View ArticlePubMedGoogle Scholar
- Gorski B, Jakubowska A, Huzarski T, Byrski T, Gronwald J, Grzybowska E, Mackiewicz A, Stawicka M, Bebenek M, Sorokin D, Fiszer-Maliszewska L, Haus O, Janiszewska H, Niepsuj S, Gozdz S, Zaremba L, Posmyk M, Pluzanska M, Kilar E, Czudowska D, Wasko B, Miturski R, Kowalczyk JR, Urbanski K, Szwiec M, Koc J, Debniak B, Rozmiarek A, Debniak T, Cybulski C, Kowalska E, Toloczko-Grabarek A, Zajaczek S, Menkiszak J, Medrek K, Masojc B, Mierzejewski M, Narod SA, Lubinski J: A high proportion of founder BRCA1 mutations in Polish breast cancer families. x. 2004, 110 (5): 683-686.Google Scholar
- Houlston RS, Peto J: The future of association studies of common cancers. Hum Genet. 2003, 112 (4): 434-435.PubMedGoogle Scholar
- Antoniou AC, Easton DF: Polygenic inheritance of breast cancer: Implications for design of association studies. Genet Epidemiol. 2003, 25 (3): 190-202. 10.1002/gepi.10261.View ArticlePubMedGoogle Scholar
- HWE test tool. [http://ihg.gsf.de/cgi-bin/hw/hwa1.pl]
- Power and Sample Size Calculation. [http://biostat.mc.vanderbilt.edu/twiki/bin/view/Main/PowerSampleSize]
- Dupont WD, Plummer WD: Power and sample size calculations for studies involving linear regression. Control Clin Trials. 1998, 19 (6): 589-601. 10.1016/S0197-2456(98)00037-3.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.