Identification of polymorphisms and balancing selection in the male infertility candidate gene, ornithine decarboxylase antizyme 3
© Christensen et al; licensee BioMed Central Ltd. 2006
Received: 09 December 2005
Accepted: 16 March 2006
Published: 16 March 2006
The antizyme family is a group of small proteins that play a role in cell growth and division by regulating the biosynthesis of polyamines (putrescine, spermidine, spermine). Antizymes regulate polyamine levels primarily through binding ornithine decarboxylase (ODC), an enzyme key to polyamine production, and targeting ODC for destruction by the 26S proteosome. Ornithine decarboxylase antizyme 3 (OAZ3) is a testis-specific antizyme paralog and the only antizyme expressed in the mid to late stages of spermatogenesis.
To see if mutations in the OAZ3 gene are responsible for some cases of male infertility, we sequenced and evaluated the genomic DNA of 192 infertile men, 48 men of known paternity, and 34 African aborigines from the Mbuti tribe in the Democratic Republic of the Congo. The coding sequence of OAZ3 was further screened for polymorphisms by SSCP analysis in the infertile group and an additional 250 general population controls. Identified polymorphisms in the OAZ3 gene were further subjected to a haplotype analysis using PHASE 2.02 and Arlequin 2.0 software programs.
A total of 23 polymorphisms were identified in the promoter, exons or intronic regions of OAZ3. The majority of these fell within a region of less than two kilobases. Two of the polymorphisms, -239 A/G in the promoter and 4280 C/T, a missense polymorphism in exon 5, may show evidence of association with male infertility. Haplotype analysis identified 15 different haplotypes, which can be separated into two divergent clusters.
Mutations in the OAZ3 gene are not a common cause of male infertility. However, the presence of the two divergent haplotypes at high frequencies in all three of our subsamples (infertile, control, African) suggests that they have been maintained in the genome by balancing selection, which was supported by a test of Tajima's D statistic. Evidence for natural selection in this region implies that these haplotypes may be associated with a trait other than infertility. This trait may be related to another function of OAZ3 or a region in tight linkage disequilibrium to the gene.
Polyamines are small, ubiquitous organic molecules derived from arginine or methionine. They carry a positive charge on each nitrogen atom, which allows them to interact with polyanionic macromolecules like DNA, RNA and some proteins . Often, these interactions are very specific, and play critical roles in normal cell growth and function . Abnormal expression of polyamines has been associated with both apoptosis and tumor growth . Consequently, polyamine concentrations are tightly regulated at transcriptional, translational, and post-translational levels .
Ornithine decarboxylase (ODC), which converts ornithine to the first polyamine, putrescine, is the rate-limiting step in polyamine biosynthesis and the key point at which polyamine levels are regulated . In eukaryotic cells, the most important regulator of ODC is the protein ornithine decarboxylase antizyme (antizyme) . When polyamine levels rise to a critical level they induce a +1 ribosomal frameshift in antizyme mRNA, producing a functional antizyme protein . The antizyme protein binds ODC, preventing formation of the enzymatically-active ODC homodimer and targeting ODC for ubiquitin-independent proteolysis by the 26S proteosome . In addition to inhibiting the production of intracellular polyamines, antizyme can also inhibit the import of extracellular polyamines .
Antizyme genes have been identified from yeast to mammals. Most invertebrates have a single antizyme gene, while there are at least three independently conserved antizyme isoforms in vertebrates . Three main antizymes have been described in mammals. Antizyme 1 (AZ1) and antizyme 2 (AZ2) are generally thought to be present in all tissues except haploid male germ cells, with AZ1 mRNA being present in a 10- to 20- fold excess of AZ2 . AZ1 is more active in the degradation of ODC than AZ2, and both contribute to the negative regulation of polyamine transport .
Antizyme 3 (OAZ3) is expressed predominanty in round and elongating spermatids. Though AZ1 is expressed in early, diploid testicular germ cells, OAZ3 is the only antizyme present in post-meiotic, male germ cells [11–13]. Expression of OAZ3 follows closely behind ODC expression, which is tightly regulated in the testis, with a sharp spike in expression occurring during the late spermatocyte and early spermatid stages . Transgenic male mice that overexpress ODC, resulting in increased synthesis of polyamines are infertile, with an almost complete lack of mature spermatozoa [14, 15] This suggests that OAZ3 might play an important role in maintaining spermatogenesis by controlling ODC and polyamine levels. A recent report indicates that OAZ3 may also interact with at least one additional protein in the testis, gametogenetin 1, , in an unspecified manner. The function of gametogenetin 1 is unknown, but its testis-specific expression is similar to OAZ3 both spatially and temporally.
Infertility affects approximately 10% of couples, and in about half of these cases a male factor problem can be identified . The etiology of male factor infertility is poorly understood. While some can be explained due to Y-chromosome microdeletions, endocrine disruptions, developmental abnormalities or environmental insults, the majority remains idiopathic and potentially genetic in origin . Multiple studies, primarily through transgenic animal models, have helped develop a list of candidate genes that may contribute to human, male infertility [19, 20]. Given its testis-specific expression and the clear link between deregulation of polyamine expression and infertility , OAZ3 is a clear choice for a male infertility candidate gene.
To determine if mutations in OAZ3 account for some cases of male infertility we screened 192 patients and 334 fertile or general population controls by either single-strand conformational polymorphism analysis (SSCP) or directly sequencing the exons and flanking intronic sequence of OAZ3. In addition to screening the OAZ3 gene for possible mutations associated with infertility we have also examined the frequency of haplotypes that emerged and the evolutionary relationships among the haplotypes. Finally, to examine the relationship between variation in this gene and disease phenotype, we tested for associations between haplotype and phenotype.
Patient and control DNA samples
After obtaining Institutional Review Board approval, a total of 192 men presenting with primary infertility were enrolled at the University of Utah School of Medicine (Salt Lake City, UT) or the Weill-Cornell Medical Center (New York, NY). The patients were either azoospermic, having no sperm in the ejaculate (n = 143) or severely oligozoospermic, having less than five million sperm/ml. of semen (n = 49). These patients were selected because their phenotypes match those observed in mice overexpressing ODC, which exhibit both azoospermia, primarily due to maturation arrest, and oligozoospermia [14, 15]. To increase the chance of identifying a mutation in the OAZ3 gene potentially responsible for male infertility or subfertility, and to decrease possible confounding factors, patients were excluded if they had any suspected or known causes of male infertility. For example, patients with a known Y-chromosome microdeletion, cystic fibrosis, varicocele, Klinefelter's, exposure to chemotherapeutics or radiation etc., were not included. Venous blood was drawn from each patient using standard techniques and DNA extracted using the Puregene DNA extraction kit (Genzyme, Minneapolis, MN).
In addition to the infertile group, all or part of OAZ3 was screened in three different control groups. The first group (n = 48), representing men with established paternity, was obtained from the Utah Genetic Reference Project (UGRP). These individuals are primarily of Northern European ancestry and are genetically very similar to the infertile group. The second group was a sample of central African aborigines (n = 34), all members of the Mbuti tribe from the Democratic Republic of the Congo. The group is predominantly male, though there are a few females in the group. It is unknown which, if any of the group, are infertile. They were included in the study because they represent a group of genetically diverse individuals, geographically isolated from the infertile and fertile study groups, which is useful in determining how haplotypes may have evolved. The infertile and African groups were screened identically to the patient group by direct sequencing. A third and final general population group of men (n = 250), which included 90 African and 80 Asian males, was screened by SSCP analysis, specifically for the Pro164Ser change identified in one of the patient samples, but also for mutations in the coding regions in the rest of the gene. To assess evolutionary relationships with non-human primates, DNA was obtained and sequenced for both the common chimpanzee (Pan troglodytes) and gorilla (Gorilla gorilla).
PCR and sequencing primers for the OAZ3 gene. An "F" in the name designates a forward primer, an "R" designates a reverse primer. The position given for each primer is that of the first nucleotide, in relation to the first nucleotide of the ATG start codon.
PCR Primers 5'-3'
Sequencing Primers 5'-3'
EXON 4–5, 3'UTR
(PCR fragment includes 487 bp of the 3'UTR)
Sequence trace files generated from the ABI 3700s were analyzed using the Phred, Phrap, and Consed software programs . Phred assigns a quantitative value to the quality of each sequenced base. This base quality provides a probabilistic estimate of the correctness of the base call. The sequences were assembled, and a consensus sequence generated from the most common base calls using the Phrap program, and potential mutations identified using Consed, which has the ability to search for high quality base discrepancies, based on the Phred values, in the assembled sequence. A visual analysis was also conducted of the trace files to confirm identified polymorphisms and potential mutations. In addition to direct sequencing, the five coding exons of the 192 patient samples and exon 5 of the 250-member panel were screened for variants by SSCP. This was done according to methods described previously .
Haplotypes were inferred using PHASE 2.02 . Evolutionary relationships among haplotypes were inferred using Arlequin 2.0  to generate minimum spanning trees. Population genetic diversity was measured using two standard statistics, S and pi. S is the number of variable nucleotide positions in a sample, pi is the mean pairwise difference (per nucleotide) between sampled sequences. Pairwise differences among sequences were calculated by counting the number of nucleotide positions distinguishing each possible combination of the 2*N sequences in each sample, where N is the number of individuals in the sample. Tests of evolutionary neutrality were performed using Tajima's D statistic . Tajima's D compares S and pi to identify departures from neutrality and constant population size. Because much evidence suggests that modern human populations are the product of a major population expansion in the Upper Pleistocene [27–29], we took population size change into account in D tests using the method of Wooding et al. .
Mutational analysis of OAZ3
Nucleotide variants detected in the human OAZ3 gene
Amino acid change
-239 A – G
1079 A – G
2992 G – C*
3043 A – G
3068 G – A
3079 C – 0
3153 A – G
See Figure 1
3280 T – A
3420 T – G
3427 T – G
3457 T – C
3595 G – T
3619 G – A
3710 A – G
3735 C – T
3802 C – T
3934 T – G
4023 C – A
4280 C – T
4394 A – C
4487 G – C
4531 A – 0
Ten of the identified variants (2, 3, 5, 7, 12, 14, 17–19, 23) (Table 2) had been reported previously in dbSNPs . Five of the remaining 13 nucleotide variants, numbers 1, 9, 11, 20, and 22 were unique to the infertile group while four variants, numbers 13, 15, 16 and 21 were unique to the aborigines. No significant differences were observed in allele frequencies between the infertile and control group. Significant differences in allele frequencies were found between the African group and both of the other two groups for all polymorphisms identified that were present in more than one individual (p < .001).
Assuming OAZ3 plays an integral role in spermatogenesis, At least two of the identified variants, variants 1 and 20, may be candidates for involvement in male infertility or subfertility. Variant 1, -239A→G, identified in a heterozygous state in the 5' untranslated region (5'UTR) of an azoospermic patient, generates a new potential translation initiation ATG codon within a Kozak consensus sequence and occurs two basepairs downstream from a putative Inr binding site. Though translation initiated at this position encounters a stop codon just prior to the normal translation initiation site, it could potentially reduce the amount of OAZ3 mRNA available to translate. In this way it may reduce the total amount of OAZ3 in the testis and affect fertility.
Of the 23 variants identified by direct sequencing, 16 were used for haplotype analysis. The 7 variants not used consisted of the 3 insertion/deletion polymorphisms (6, 8, 23) and 4 variants (1–4) for which complete data was not available in all samples. Of the 16 variable nucleotides used, separate analyses of the affected, control, and African subsamples revealed 12, 7, and 12 variable nucleotide positions, respectively.
When haplotypes 5, 9, and 12 were aligned, it was apparent that haplotype 9 was derived from haplotypes 5 and 12. This appears to have taken place as the result of a gene conversion event coupled with a recombination event (Figure 1). No intermediary haplotypes between 5 and 9 or 12 and 9 were detected, suggesting that the gene conversion and recombination event from which 9 was derived must have occurred very close together historically.
The population of infertile patients we screened displayed phenotypes consistent with what we expected to observe if they harbored debilitating mutations in OAZ3. Our results indicate that mutations in the OAZ3 gene are not strongly associated with human male infertility, though they may contribute to infertility in isolated instances. It is also possible that mutations in OAZ3 result in a phenotype different than the one we selected for. If that is the case, information from the OAZ3 knockout mouse, currently in progress, may provide a clue as to which human phenotype should be screened. One patient had a heterozygous change of a single, nonpolar proline residue to a charged, polar serine residue in exon five, Pro164Ser. It is not known if the proline residue is part of a functional domain or how its alteration might affect the function of the OAZ3 protein. It is conserved in all of the other OAZ3 sequences examined (Gorilla, chimpanzee, rat, mouse). We are aware of at least four proteins that interact with OAZ3 (unpublished data) and the proline residue may be important for correct protein interaction. The variation was not identified in any of the additional 191 infertile or 332 control individuals screened, indicating it is very rare. Currently, there is no direct evidence that the Pro164Ser mutation affects OAZ3 function. Additional studies will need to be conducted to determine how this mutation affects OAZ3 and spermatogenesis. Variant 1, a -239A→G change in the 5' UTR of one infertile patient may also be of significance, as it was not identified in the fertile control group and might lead to reduced levels of OAZ3 being translated. The -239A→G variant also occurs two basepairs downstream from a potential Inr transcription factor binding site , which may possibly attenuate the transcription factor's interaction with the binding sequence.
Based on the NCBI SNPs database, at least 10 additional SNPs, which we did not identify, are located in the region we screened. Only one of these, rs 5777784, appears to have potential clinical significance. It shows the deletion of a guanidine in exon three, which would cause a frameshift mutation in OAZ3 and a truncated protein. As we did not detect this SNP in any of our samples, we cannot verify its existence.
The presence of two divergent haplotype clusters in our sample suggested that either balancing natural selection or population subdivision has been active in our study population. Both of these factors can maintain divergent alleles for extended periods of time . Balancing natural selection, for example, can maintain divergent alleles when heterozygotes have a fitness advantage, while population subdivision can maintain divergent alleles in different populations. However, the presence of the two divergent haplotypes at high frequencies in all three of our subsamples (infertile, control, and African) shows that subdivision cannot account for the presence of these divergent lineages. A better explanation is that balancing natural selection has historically been active in this region of the genome, or in some region in tight linkage disequilibrium with it.
To test whether balancing selection might be occurring in a region connected to OAZ3, we reviewed SNPs identified in OAZ3s nearest neighbors that might have functional implications for the protein. Only one SNP, an A→G that converts a histidine residue to arginine in the final exon of the tudor and KH domain containing (TDRKH) gene, was identified. The 3' end of TDRKH, located on the opposite strand of OAZ3, overlaps with the 3' end of OAZ3, and the SNP in question lies in close proximity to the region from which the haplotypes were derived. The region containing the TDRKH SNP was amplified by PCR in 12 samples, six from each of the two major haplotypes and sequenced. While the SNP was identified in several of the patients, it appeared in both major haplotypes, indicating a different allele in the region is the one being maintained by balancing selection.
Evidence that balancing natural selection may have been active in the OAZ3 region is further supported by tests of Tajima's D statistic . Tajima's D compares estimates of diversity based on S and pi, which are affected differently by natural selection. Positive natural selection, for instance, which drives a new favored variant rapidly to high frequency, tends to reduce pi relative to S, causing negative values of D. In contrast, balancing natural selection tends to increase pi relative to S, causing positive values of D.
Initial tests of Tajima's D statistic showed that the D values observed in our sample, while positive, were not significantly different from neutral expectation (Figure 5). However, these initial tests were performed under the assumption that human population sizes have been constant. This assumption is likely inappropriate, because archaeological and genetic evidence suggest that human populations expanded greatly – more than 100-fold – beginning roughly 100,000 years ago [33, 34]. For this reason, we retested Tajima's D statistic under the conservative assumption that the human population increased 100-fold, from an initial effective population size of 10,000 (a typical estimate; ), 100,000 years ago. Under this assumption, the observed values of D departed from expectation in every subsample, and in the sample as a whole. Thus, patterns of variation in this region are consistent with the hypothesis that balancing natural selection has been active.
One possible explanation for balancing selection is heterozygote advantage, which can result from increased homozygote mortality or increased heterozygote reproductive success. To test for this possibility we performed Chi-squared tests of Hardy-Weinberg equilibrium on each SNP in our sample. Only two variants, variant 18 in the infertile group and variant 7 in the African group were significantly different than expected (p < .01). These tests show no evidence of a heterozygote excess and suggest that balancing selection in OAZ3 would be better explained by differences in pure reproductive success than by differences in survivorship. In the case of favored reproductive success, no excess of heterozygotes is expected, as random matings will return the population to Hardy-Weinberg equilibrium.
Signatures of natural selection can only arise when some sort of variation with phenotypic consequences, upon which selection can act, is present. Thus, evidence that balancing natural selection has been active in OAZ3 suggests that some kind of functional variation is present in or near this region. This variation is likely in linkage disequilibrium with the genetic variants that distinguish the two major haplotypes, which seem to have been maintained by selection. While the similar frequencies of these divergent haplotypes in our affected and control samples suggests that they are not likely associated with the phenotype examined here, they could be associated with other phenotypes related to OAZ3.
Polymorphisms in OAZ3 are not strongly associated with human male infertility, though some of the identified SNPs, like Pro164Ser in exon 5, may contribute to isolated cases of infertility. The presence of two divergent haplotype clusters in OAZ3, in all subsamples tested (infertile, control, African), implies the haplotypes are being maintained by balancing natural selection. Natural selection in this area indicates these haplotypes may be associated with a trait other than infertility. This trait may be related to either a different function of OAZ3 or a region in tight linkage disequilibrium to the OAZ3 gene.
The African genomic DNA, 250 member panel used for SSCP and the chimpanzee and gorilla DNA were kindly provided by Lynn Jorde, Ph.D. Control samples obtained from the Utah Genetic Reference Project (UGRP) were collected with support from a Public Health Services research grant to the Huntsman General Clinical Research Center, number M01-RR00064, from the National Center for Research Resources. The UGRP was also supported by generous gifts from the W.M. Keck and Delores Dore Eccles Foundations. We'd like to extend our sincere thanks to all family members who participated in the UGRP and the individuals who consented to DNA collection for the infertile panel.
- Wallace HM, Fraser AV, Hughes A: A perspective of polyamine metabolism. Biochem J. 2003, 376: 1-14. 10.1042/BJ20031327.View ArticlePubMedPubMed CentralGoogle Scholar
- Thomas T, Thomas TJ: Polyamine metabolism and cancer. J Cell Mol Med. 2003, 7: 113-126.View ArticlePubMedGoogle Scholar
- Gerner EW, Meyskens FLJ: Polyamines and cancer: old molecules, new understanding. Nat Rev Cancer. 2004, 4: 781-792. 10.1038/nrc1454.View ArticlePubMedGoogle Scholar
- Mangold U, Leberer E: Regulation of all members of the antizyme family by antizyme inhibitor. Biochem J. 2005, 385: 21-28. 10.1042/BJ20040547.View ArticlePubMedGoogle Scholar
- Hayashi S, Murakami Y: Rapid and regulated degradation of ornithine decarboxylase. Biochem J. 1995, 306 ( Pt 1): 1-10.View ArticleGoogle Scholar
- Hayashi S, Murakami Y, Matsufuji S: Ornithine decarboxylase antizyme: a novel type of regulatory protein. Trends Biochem Sci. 1996, 21: 27-30. 10.1016/0968-0004(96)80882-6.View ArticlePubMedGoogle Scholar
- Matsufuji S, Matsufuji T, Miyazaki Y, Murakami Y, Atkins JF, Gesteland RF, Hayashi S: Autoregulatory frameshifting in decoding mammalian ornithine decarboxylase antizyme. Cell. 1995, 80: 51-60. 10.1016/0092-8674(95)90450-6.View ArticlePubMedGoogle Scholar
- Mitchell JL, Judd GG, Bareyal-Leyser A, Ling SY: Feedback repression of polyamine transport is mediated by antizyme in mammalian tissue-culture cells. Biochem J. 1994, 299 ( Pt 1): 19-22.View ArticleGoogle Scholar
- Coffino P: Regulation of cellular polyamines by antizyme. Nat Rev Mol Cell Biol. 2001, 2: 188-194. 10.1038/35056508.View ArticlePubMedGoogle Scholar
- Ivanov IP, Gesteland RF, Atkins JF: A second mammalian antizyme: conservation of programmed ribosomal frameshifting. Genomics. 1998, 52: 119-129. 10.1006/geno.1998.5434.View ArticlePubMedGoogle Scholar
- Ivanov IP, Rohrwasser A, Terreros DA, Gesteland RF, Atkins JF: Discovery of a spermatogenesis stage-specific ornithine decarboxylase antizyme: antizyme 3. Proc Natl Acad Sci U S A. 2000, 97: 4808-4813. 10.1073/pnas.070055897.View ArticlePubMedPubMed CentralGoogle Scholar
- Tosaka Y, Tanaka H, Yano Y, Masai K, Nozaki M, Yomogida K, Otani S, Nojima H, Nishimune Y: Identification and characterization of testis specific ornithine decarboxylase antizyme (OAZ-t) gene: expression in haploid germ cells and polyamine-induced frameshifting. Genes Cells. 2000, 5: 265-276. 10.1046/j.1365-2443.2000.00324.x.View ArticlePubMedGoogle Scholar
- Ike A, Yamada S, Tanaka H, Nishimune Y, Nozaki M: Structure and promoter activity of the gene encoding ornithine decarboxylase antizyme expressed exclusively in haploid germ cells in testis (OAZt/Oaz3). Gene. 2002, 298: 183-193. 10.1016/S0378-1119(02)00978-2.View ArticlePubMedGoogle Scholar
- Hakovirta H, Keiski A, Toppari J, Halmekyto M, Alhonen L, Janne J, Parvinen M: Polyamines and regulation of spermatogenesis: selective stimulation of late spermatogonia in transgenic mice overexpressing the human ornithine decarboxylase gene. Mol Endocrinol. 1993, 7: 1430-1436. 10.1210/me.7.11.1430.PubMedGoogle Scholar
- Halmekyto M, Hyttinen JM, Sinervirta R, Utriainen M, Myohanen S, Voipio HM, Wahlfors J, Syrjanen S, Syrjanen K, Alhonen L, et al: Transgenic mice aberrantly expressing human ornithine decarboxylase gene. J Biol Chem. 1991, 266: 19746-19751.PubMedGoogle Scholar
- Zhang J, Wang Y, Zhou Y, Cao Z, Huang P, Lu B: Yeast two-hybrid screens imply that GGNBP1, GGNBP2 and OAZ3 are potential interaction partners of testicular germ cell-specific protein GGN1. FEBS Lett. 2005, 579: 559-566. 10.1016/j.febslet.2004.10.112.View ArticlePubMedGoogle Scholar
- Okabe M, Ikawa M, Ashkenas J: Male infertility and the genetics of spermatogenesis. Am J Hum Genet. 1998, 62: 1274-1281. 10.1086/301895.View ArticlePubMedPubMed CentralGoogle Scholar
- Shah K, Sivapalan G, Gibbons N, Tempest H, Griffin DK: The genetic basis of infertility. Reproduction. 2003, 126: 13-25. 10.1530/rep.0.1260013.View ArticlePubMedGoogle Scholar
- Matzuk MM, Lamb DJ: Genetic dissection of mammalian fertility pathways. Nat Cell Biol. 2002, 4 Suppl: s41-9.PubMedGoogle Scholar
- Christensen GL, Carrell DT: Animal models of genetic causes of male infertility. Asian J Androl. 2002, 4: 213-219.PubMedGoogle Scholar
- Pegg AE, Feith DJ, Fong LY, Coleman CS, O'Brien TG, Shantz LM: Transgenic mouse models for studies of the role of polyamines in normal, hypertrophic and neoplastic growth. Biochem Soc Trans. 2003, 31: 356-360. 10.1042/BST0310356.View ArticlePubMedGoogle Scholar
- Laboratory of Phil Green. [http://www.phrap.org]
- Splawski I, Timothy KW, Vincent GW, Atkinson MT: Molecular basis of the long-QT syndrome associated with deafness. N Engl J Med. 1997, 336 (22): 1562-1567. 10.1056/NEJM199705293362204.View ArticlePubMedGoogle Scholar
- Stephens M, Smith NJ, Donnelly P: A new statistical method for haplotype reconstruction from population data. Am J Hum Genet. 2001, 68: 978-989. 10.1086/319501.View ArticlePubMedPubMed CentralGoogle Scholar
- Schneider S, Roessli D, Excoffier L: ARLEQUIN version 2.000: a software for population genetics data analysis. Genetics and Biometry Laboratory, Department of Anthroplogy, University of Geneva. Geneva Switzerland. 2000Google Scholar
- Tajima F: Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989, 123: 585-595.PubMedPubMed CentralGoogle Scholar
- Harpending H, Rogers A: Genetic perspectives on human origins and differentiation. Annu Rev Genomics Hum Genet. 2000, 1: 361-385. 10.1146/annurev.genom.1.1.361.View ArticlePubMedGoogle Scholar
- Excoffier L: Human demographic history: refining the recent African origin model. Curr Opin Genet Dev. 2002, 12: 675-682. 10.1016/S0959-437X(02)00350-7.View ArticlePubMedGoogle Scholar
- Tishkoff SA, Verrelli BC: Patterns of human genetic diversity: implications for human evolutionary history and disease. Annu Rev Genomics Hum Genet. 2003, 4: 293-340. 10.1146/annurev.genom.4.070802.110226.View ArticlePubMedGoogle Scholar
- Wooding S, Kim UK, Bamshad MJ, Larsen J, Jorde LB, Drayna D: Natural selection and molecular evolution in PTC, a bitter-taste receptor gene. Am J Hum Genet. 2004, 74: 637-646. 10.1086/383092.View ArticlePubMedPubMed CentralGoogle Scholar
- National Center for Biotechnology Information. [http://www.ncbi.nih.gov]
- Bamshad M, Wooding SP: Signatures of natural selection in the human genome. Nat Rev Genet. 2003, 4: 99-111. 10.1038/nrg999.View ArticlePubMedGoogle Scholar
- Klein EK, Austerlitz F, Laredo C: Some statistical improvements for estimating population size and mutation rate from segregating sites in DNA sequences. Theor Popul Biol. 1999, 55: 235-247. 10.1006/tpbi.1998.1401.View ArticlePubMedGoogle Scholar
- Rogers A: Population forecasting: do simple models outperform complex models?. Math Popul Stud. 1995, 5: 187-202, 291.View ArticlePubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2350/7/27/prepub
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.