Primate-specific evolution of noncoding element insertion into PLA2G4Cand human preterm birth
- Jevon Plunkett1, 2,
- Scott Doniger3,
- Thomas Morgan1, 4,
- Ritva Haataja5,
- Mikko Hallman5,
- Hilkka Puttonen6,
- Ramkumar Menon7, 8,
- Edward Kuczynski9,
- Errol Norwitz9,
- Victoria Snegovskikh9,
- Aarno Palotie10, 11, 12, 13,
- Leena Peltonen10, 11, 12,
- Vineta Fellman14, 15,
- Emily A DeFranco16,
- Bimal P Chaudhari17,
- John Oates18,
- Olivier Boutaud18,
- Tracy L McGregor1,
- Jude J McElroy1,
- Kari Teramo6,
- Ingrid Borecki19,
- Justin C Fay20 and
- Louis J Muglia1, 21Email author
© Plunkett et al; licensee BioMed Central Ltd. 2010
Received: 22 July 2010
Accepted: 24 December 2010
Published: 24 December 2010
The onset of birth in humans, like other apes, differs from non-primate mammals in its endocrine physiology. We hypothesize that higher primate-specific gene evolution may lead to these differences and target genes involved in human preterm birth, an area of global health significance.
We performed a comparative genomics screen of highly conserved noncoding elements and identified PLA2G4C, a phospholipase A isoform involved in prostaglandin biosynthesis as human accelerated. To examine whether this gene demonstrating primate-specific evolution was associated with birth timing, we genotyped and analyzed 8 common single nucleotide polymorphisms (SNPs) in PLA2G4C in US Hispanic (n = 73 preterm, 292 control), US White (n = 147 preterm, 157 control) and US Black (n = 79 preterm, 166 control) mothers.
Detailed structural and phylogenic analysis of PLA2G4C suggested a short genomic element within the gene duplicated from a paralogous highly conserved element on chromosome 1 specifically in primates. SNPs rs8110925 and rs2307276 in US Hispanics and rs11564620 in US Whites were significant after correcting for multiple tests (p < 0.006). Additionally, rs11564620 (Thr360Pro) was associated with increased metabolite levels of the prostaglandin thromboxane in healthy individuals (p = 0.02), suggesting this variant may affect PLA2G4C activity.
Our findings suggest that variation in PLA2G4C may influence preterm birth risk by increasing levels of prostaglandins, which are known to regulate labor.
A growing body of evidence supports genetic influences on preterm birth risk; however, few genes have been consistently associated with the disorder [1, 2]. Investigators have typically focused on candidate genes selected based on predicted parturition physiology; however, this approach may be limited by the divergence in physiological mechanisms between humans and model organisms that have been typically studied. For example, while a rapid decline in progesterone plays a prominent role in initiating parturition in rodents and sheep, this signal does not seem to precede human labor . Other parturition-related traits, such as placental morphology and source of progesterone, also differ importantly in humans compared to model organisms typically studied and may limit what generalizations can be made .
Differences in parturition physiology between apes, including humans, and other mammals may have developed in response to uniquely human adaptations including relatively large human head size and narrow birth canal cross-sectional area . Genes involved in parturition likely have evolved differentially along the human and/or higher primate phylogenetic lineages to decrease the length of gestation and alleviate the complications arising from such cephalopelvic constraints. As a result, the set of genes rapidly evolving on the human and/or higher primate lineage likely includes genes that play important roles in regulating parturition and potentially influence preterm birth risk. Consistent with our hypothesis, we identified FSHR as having rapidly evolved by nucleotide substitution and as being associated with preterm birth risk across independent populations ( and (Plunkett J, Doniger S, Orabona G, Morgan T, Haataja R, Hallman M, Puttonen H, Menon R, Kuczynski E, Norwitz E et al: Evolutionary history of FSHR in human predicts role in birth time, submitted)).
In addition to nucleotide substitution, genomic rearrangements account for a substantial portion of genomic divergence among species. For example, Frazer et al.  and Wetterbom et al.  observed insertions and deletions frequently when comparing genome sequences among humans, chimpanzees and other primate species. Moreover, such rearrangements may account for a larger fraction of genomic divergence than nucleotide substitutions . Rearrangements can lead to loss or acquisition of exons, splice sites and promoters, facilitating differences in expression patterns, such as those observed for transcript variants of CHRM3 and SFTPB with differing transposable element insertion events [8, 9]. Hence, genomic rearrangement may contribute to rapid evolution along the human and/or higher primate lineages in response to unique physiological constraints.
We hypothesized that genes with genomic rearrangements departing from the ancestral state and occurring on the human and/or higher primate lineages may play important roles in birth timing and preterm delivery. Thus, we investigated association with preterm birth for common variants in a gene, PLA2G4C, which is expressed in the uterus  and involved in prostaglandin synthesis, suggesting a potential role in parturition, and in which we have identified a primate-specific insertion.
Evolutionary history of a primate-specific PLA2G4Cnoncoding element
Association with preterm birth
Additionally, 2, 3 and 4 SNP haplotypes containing SNPs rs8110925 and rs2307276 were significant in the US Hispanics after correcting for 18 haplotype comparisons (p < 0.003), although not more significant than single SNP association findings (Additional file 3 Table S3, S4). 2 SNP haplotypes containing rs11564620 were moderately significant (p < 0.05) in US Whites (Additional file 3 Table S3, Additional file 4 Table S4). Linkage disequilibrium (LD) among SNPs rs8110925, rs2307276, and rs11564620 was very low (r2 < 0.1) in the three populations studied (Additional file 5 Figure S1), suggesting multiple independent associations were observed.
Association with prostaglandin concentrations
Comparative genomic analysis is an attractive method for identifying genetic variation among species that may correlate with inter-specific phenotypic variation in fundamental processes such as parturition. Through such an analysis, we have identified a noncoding element in intron 14 of PLA2G4C on chromosome 19 representing a primate-specific change involving amplification and subsequent divergence but without increased nucleotide substitution. Having identified PLA2G4C as a candidate gene, we proposed that this duplicated element represents a primate-specific change with a potential regulatory role in human parturition.
Case-control association results for significant and suggestive SNPs in the PLA2G4C gene region tested across 3 independent populations.
Allelic OR (95% CI)
(73 cases, 292 controls)
7.92 × 10 -4 a, b
5.66 × 10 -5 b
5.45 × 10 -3 b
(147 cases, 157 controls)
6.98 × 10 -3
1.03 × 10 -3 b
(79 cases, 166 controls)
PLA2G4C encodes cytosolic phospholipase A2 gamma, which hydrolizes phospholipids from the cellular membrane to form free arachidonic acid, from which prostaglandins, including prostaglandins D, E, F, I2 (also known as prostacyclin), and thromboxane A2 are generated. Prostaglandins play an important role in parturition. Pharmacologically, prostaglandins are used to induce abortion, for cervical ripening, and labor induction and drugs inhibiting prostaglandin synthesis are successful in preventing preterm labor . Levels of prostaglandins, including thromboxane A2, are elevated in pregnant compared to non-pregnant women, and in late (36 weeks) compared to early (20, 30 weeks) gestation , suggesting a link between prostaglandin abundance and parturition timing. Prostaglandins may facilitate labor by several mechanisms. These hormones are known uterotonic agents and also promote luteolysis and the onset of labor in species that exhibit progesterone withdrawal prior to birth . Prostaglandins may also facilitate delivery by affecting placenta function, since thromboxane A2 induces platelet aggregation and acts as a vasoconstrictor . Hence, higher prostaglandin levels than expected may initiate parturition prematurely and lead to preterm delivery.
The PLA2G4C enzyme is the only cytosolic phospholipase A2 family member that is constitutively associated with the cellular membrane, the site of prostaglandin synthesis, rather than translocating to the membrane in response to calcium signaling . Hence, dysregulation of PLA2G4C may alter prostaglandins levels independent of other parturition signals, such as oxytocin , that act via intracellular calcium signaling. For example, rs11564620 may contribute to a conformational change in the enzyme's active site, rendering it more active than usual and leading to increased synthesis of prostaglandins, as demonstrated by our observation of higher levels of thromboxane A2 in minor allele carriers for this polymorphism. Moreover, multiple splice isoforms of PLA2G4C exist, differing in transcript length, presence of certain exons and overlapping exons with different boundaries (AceView, NCBI, http://www.ncbi.nlm.nih.gov/IEB/Research/Acembly/). As a result, variation in PLA2G4C may contribute to differences in tissue-specific expression or relative abundance of various PLA2G4C isoforms, potentially altering function. Further study of the region encompassing these SNPs, including the genomic element of interest, is needed to examine the mechanism by which variation in PLA2G4C influences birth timing.
Specialization within multi-gene families, like the large phospholipase A2 gene family, can create individualized functions among paralogous genes. For example, PLA2G4C has a continuous association with the cellular membrane, unlike other phospholipase A2 genes, potentially differentiating its role in prostaglandin synthesis from those of other family members. Genomic variation, such as the element insertion observed in PLA2G4C, may contribute to gene specialization, as demonstrated by divergence in PLA2G4C expression patterns in humans versus mice, who lack the element insertion and express PLA2G4C only in ovary and oocytes. Specialized genes are potentially better therapeutic targets than gene products with multiple roles within cell, since pharmaceutically targeting such genes may lead to fewer side effects. As a result, PLA2G4C may be a useful target for designing novel therapies to prolong pregnancy and reduce the incidence of preterm birth.
A higher primate-specific noncoding element insertion into intron 14 of the phospholipase A2 gene PLA2G4C was identified, demonstrating the gene's rapid evolution along the higher primate lineage by genomic rearrangement. Results from our genetic analysis suggest common variation in PLA2G4C influences preterm birth risk. One of the variants associated with preterm birth is a nonsynonymous coding polymorphism, rs11564620, predicted to be potentially altering protein structure. This coding polymorphism is also associated with thromboxane levels, suggesting that genetic variation in PLA2G4C may increase risk for preterm birth by increasing levels of prostaglandins, which are known to regulate labor. By examining rapid evolution along human and higher primate lineages by genomic rearrangement, we have identified a novel gene associated with preterm birth. This approach can be readily applied to other traits differing among humans and/or higher primates and other species to aid in gene discovery.
Genomic alignments to investigate evolutionary history of PLA2G4C
Noting the deletion of PLA2G4C reported in the Taurine Cattle genome  contrasted with its presence in the 17-way MultiZ alignments  we used to identify the gene as rapidly evolving (analysis conducted Spring 2007 and presented in detail in (Plunkett J, Doniger S, Orabona G, Morgan T, Haataja R, Hallman M, Puttonen H, Menon R, Kuczynski E, Norwitz E et al: Evolutionary history of FSHR in human predicts role in birth time, submitted)), we examined the history of this region in greater depth. We extracted sequence surrounding the 130 bp highly conserved noncoding element (human chromosome 19: 48,560,500 -48,560,630; hg19 genome build) which largely contributed to our designation of PLA2G4C as rapid evolving along the human lineage. A BLASTN search of the element revealed highly identical conserved noncoding elements on human chromosomes 1 (87% identity) and 2 (85% identity) (Figure 1B). We compared the human chromosome 19 and chromosome 1 sequences to 31 eutherian mammalian genomes using Ensembl Genomic alignments (accessed September 2009), and ClustalW alignment, and to specific primate genomes using BLASTN searches of human, chimpanzee, gorilla, orangutan, macaque, and bushbaby genomes (accessed September 2009). We then reconstructed history of the element by creating phylogenies using maximum likelihood with sequences homologous to the human chromosome 19 element (Figure 2) and coding sequences homologous to human PLA2G4C (Figure 3).
Study subjects were enrolled for genetic analysis by methods approved by Institutional Review Boards/Ethics Committees at each participating institution. Informed consent was obtained for all participants. Mothers with preterm birth were included if the birth was spontaneous (non-iatrogenic), singleton, had no obvious precipitating stimulus (trauma, infection, drug use), and was less than 37 weeks (Yale University; New York University) or 36 weeks (Centennial Hospital, Nashville, TN) of completed gestation. Control mothers were included if they had delivered two or more children at 37 weeks or later spontaneously. Healthy volunteers were recruited at Vanderbilt University for studies of prostaglandin metabolism. DNA from blood or saliva was prepared by standard methods. Race/ethnicity was assigned by self-report. All specimens were linked with demographic and medical data abstracted from maternal/neonatal records. DNA from blood or saliva was prepared by standard methods. Maternal age did not differ between cases and controls in the different populations (27.3 y vs. 28.4 y, p = 0.10 US White; 25.3 y vs. 25.2 y, p = 0.88 US Black; 26.0 y vs. 25.0 y, p = 0.20 US Hispanic).
Prostaglandin metabolite levels
For individuals enrolled in the prostaglandin study, urine was collected by standard methods. Levels of the urinary metabolites of prostaglandin E (PGE), prostaglandin I (PGI) and thromboxane (11-DTXB2) were quantified by mass spectrometry and normalized to creatinine levels, an indicator of renal function, in 44 healthy control individuals of Black, Hispanic or White race (median age 29, 60% male, 77% White).
We genotyped 14 SNPs spanning the PLA2G4C gene region (Additional file 1 Table S1) on human chromosome 19 in cohorts of US Hispanics (n = 73 preterm, 292 control mothers), US Whites (n = 147 preterm, 157 control mothers) and US Blacks (n = 79 preterm, 166 control mothers). For SNP selection, data from the HapMap Release 27 CEU population was examined in the Haploview program , using tagger and haplotype block functions, to identify regions of high LD. We selected 1 SNP per haplotype block, defined using the D' confidence interval method , having the highest minor allele frequency (MAF) in the CEU population for genotyping. We also included coding SNPs and other noncoding SNPs to improve coverage of conserved elements contributing to the gene's designation as "rapidly evolving." This selection scheme resulted in approximately 35% coverage of the gene region at r2≥0.8. SNPs showing evidence of association in one or more cohort (p < 0.01; n = 4) were then genotyped in healthy individuals on whom data on their concentrations of several prostaglandin metabolites was available to examine potential functional effects of the variants. All SNPs were genotyped using the Sequenom iPLEX massARRAY technology (Sequenom, San Diego, CA).
Data cleaning and analysis was performed with Whole-genome Association Study Pipeline (WASP)  and PLINK . We excluded individuals based on genotyping quality (< 90% call rate) and SNPs based on the following criteria: not in Hardy-Weinberg Equilibrium in controls (p < 0.001 χ2 test), <90% genotype call rate, MAF < 0.01). Linkage disequilibrium among SNPs tested was determined using the Haploview program . We chose this Hardy-Weinberg threshold for two reasons. First, we hypothesize this locus is under selective pressure which could result in some deviation from HWE in the control population. Second, samples with p < 0.001 also had low genotype call rates, suggesting genotyping error, while those with p > 0.001 had high call rates. For significant SNPs, the Hardy-Weinberg deviation was greater than 0.05 unless otherwise indicated. We corrected for multiple testing using the simpleM method , which estimates the number of independent tests, given the LD relationships among SNPs, used to obtain a Bonferroni-corrected critical value.
Our analysis considered preterm birth affection status (i.e. delivery <37 weeks) as a binary trait, comparing frequencies between case and control groups of alleles and genotypes by χ2 test. Sliding windows of 2,3 and 4 SNP haplotypes also were compared between cases and controls . Meta-analysis of data for significant SNPs was done using the Mantel-Haenszel method, after successfully passing the test of homogeneity.
To test the potential functional effect of associated PLA2G4C variants on prostaglandin metabolism, we examined the levels of PGE, PGI, and 11-DTXB2, standardized to normal distributions (μ = 0, σ = 1), as quantitative traits. A Wald test was performed to compare the mean phenotype between different allele or genotype classes for associated SNPs. We also tested whether rs11564620 risk-allele carriers had higher prostaglandin levels than noncarriers, by comparing the 11-DTXB2 value distribution among genotype classes with box plots and one-sided Wilcoxon nonparametric test performed in R .
We thank Cara Sutcliffe and Rachel Wiseman in the DNA Resources Core at Vanderbilt University Medical Center for their assistance with genotyping. This work was supported by grants from the Children's Discovery Institute at Washington University School of Medicine and St. Louis Children's Hospital awarded to JF and LJM and from the March of Dimes awarded to LJM and EN. This research was also supported by T32 GM081739 from the National Institute of General Medical Science and the Mr. and Mrs. Spencer T. Olin Fellowship for Women in Graduate Study at Washington University in St. Louis awarded to JP, a grant from the Sigrid Juselius Foundation awarded to MH, a grant from the Signe and Ane Gyllenberg foundation to VF and grants from the Academy of Finland to RH and MH. TM was supported by a Turner-Hazinski grant award from Vanderbilt University.
- Chaudhari BP, Plunkett J, Ratajczak CK, Shen TT, DeFranco EA, Muglia LJ: The genetics of birth timing: insights into a fundamental component of human development. Clin Genet. 2008, 74 (6): 493-501. 10.1111/j.1399-0004.2008.01124.x.View ArticlePubMedGoogle Scholar
- Plunkett J, Muglia LJ: Genetic contributions to preterm birth: Implications from epidemiological and genetic association studies. Ann Med. 2008, 40 (3): 167-179. 10.1080/07853890701806181.View ArticlePubMedGoogle Scholar
- Mitchell BF, Taggart MJ: Are animal models relevant to key aspects of human parturition?. Am J Physiol Regul Integr Comp Physiol. 2009, 297 (3): R525-545.View ArticlePubMedGoogle Scholar
- Rosenberg K, Trevathan W: Birth, obstetrics and human evolution. BJOG. 2002, 109 (11): 1199-1206. 10.1046/j.1471-0528.2002.00010.x.View ArticlePubMedGoogle Scholar
- Fay JC, Wyckoff GJ, Wu CI: Positive and negative selection on the human genome. Genetics. 2001, 158 (3): 1227-1234.PubMedPubMed CentralGoogle Scholar
- Frazer KA, Chen X, Hinds DA, Pant PV, Patil N, Cox DR: Genomic DNA insertions and deletions occur frequently between humans and nonhuman primates. Genome Res. 2003, 13 (3): 341-346. 10.1101/gr.554603.View ArticlePubMedPubMed CentralGoogle Scholar
- Wetterbom A, Sevov M, Cavelier L, Bergstrom TF: Comparative genomic analysis of human and chimpanzee indicates a key role for indels in primate evolution. J Mol Evol. 2006, 63 (5): 682-690. 10.1007/s00239-006-0045-7.View ArticlePubMedGoogle Scholar
- Huh JW, Kim YH, Lee SR, Kim H, Kim DS, Kim HS, Kang HS, Chang KT: Gain of new exons and promoters by lineage-specific transposable elements-integration and conservation event on CHRM3 gene. Mol Cells. 2009, 28 (2): 111-117. 10.1007/s10059-009-0106-z.View ArticlePubMedGoogle Scholar
- Lee JR, Huh JW, Kim DS, Ha HS, Ahn K, Kim YJ, Chang KT, Kim HS: Lineage specific evolutionary events on SFTPB gene: Alu recombination-mediated deletion (ARMD), exonization, and alternative splicing events. Gene. 2009, 435 (1-2): 29-35. 10.1016/j.gene.2009.01.008.View ArticlePubMedGoogle Scholar
- Tithof PK, Roberts MP, Guan W, Elgayyar M, Godkin JD: Distinct phospholipase A2 enzymes regulate prostaglandin E2 and F2alpha production by bovine endometrial epithelial cells. Reprod Biol Endocrinol. 2007, 5: 16-10.1186/1477-7827-5-16.View ArticlePubMedPubMed CentralGoogle Scholar
- Committee on Understanding Premature Birth and Assuring Healthy Outcomes: Preterm Birth: Causes, Consequences, and Prevention. 2006, Washington, D.C.: The National Academies PressGoogle Scholar
- Elsik CG, Tellam RL, Worley KC, Gibbs RA, Muzny DM, Weinstock GM, Adelson DL, Eichler EE, Elnitski L, Guigo R, et al: The genome sequence of taurine cattle: a window to ruminant biology and evolution. Science. 2009, 324 (5926): 522-528. 10.1126/science.1169588.View ArticlePubMedPubMed CentralGoogle Scholar
- Blanchette M, Kent WJ, Riemer C, Elnitski L, Smit AF, Roskin KM, Baertsch R, Rosenbloom K, Clawson H, Green ED, et al: Aligning multiple genomic sequences with the threaded blockset aligner. Genome Res. 2004, 14 (4): 708-715. 10.1101/gr.1933104.View ArticlePubMedPubMed CentralGoogle Scholar
- Kistka ZA, DeFranco EA, Ligthart L, Willemsen G, Plunkett J, Muglia LJ, Boomsma DI: Heritability of parturition timing: an extended twin design analysis. Am J Obstet Gynecol. 2008, 199 (1): 10.1016/j.ajog.2007.12.014. 43 e41 45Google Scholar
- Plunkett J, Feitosa MF, Trusgnich M, Wangler MF, Palomar L, Kistka ZA, DeFranco EA, Shen TT, Stormo AE, Puttonen H, et al: Mother's genome or maternally-inherited genes acting in the fetus influence gestational age in familial preterm birth. Hum Hered. 2009, 68 (3): 209-219. 10.1159/000224641.View ArticlePubMedPubMed CentralGoogle Scholar
- Wilcox AJ, Skaerven R, Lie RT: Familial patterns of preterm delivery: maternal and fetal contributions. Am J Epidemiol. 2008, 167 (4): 474-479. 10.1093/aje/kwm319.View ArticlePubMedGoogle Scholar
- The International HapMap Consortium: A second generation human haplotype map of over 3.1 million SNPs. Nature. 2007, 449 (7164): 851-861. 10.1038/nature06258.View ArticlePubMed CentralGoogle Scholar
- Ramensky V, Bork P, Sunyaev S: Human non-synonymous SNPs: server and survey. Nucleic Acids Research. 2002, 30 (17): 3894-3900. 10.1093/nar/gkf493.View ArticlePubMedPubMed CentralGoogle Scholar
- Ng PC, Henikoff S: SIFT: predicting amino acid changes that affect protein function. Nucleic Acids Research. 2003, 31 (13): 3812-3814. 10.1093/nar/gkg509.View ArticlePubMedPubMed CentralGoogle Scholar
- Underwood KW, Song C, Kriz RW, Chang XJ, Knopf JL, Lin LL: A novel calcium-independent phospholipase A2, cPLA2-gamma, that is prenylated and contains homology to cPLA2. J Biol Chem. 1998, 273 (34): 21926-21932. 10.1074/jbc.273.34.21926.View ArticlePubMedGoogle Scholar
- Khan AH, Carson RJ, Nelson SM: Prostaglandins in labor--a translational approach. Front Biosci. 2008, 13: 5794-5809. 10.2741/3117.View ArticlePubMedGoogle Scholar
- Hayashi M, Inoue T, Hoshimoto K, Hirabayashi H, Negishi H, Ohkura T: The levels of five markers of hemostasis and endothelial status at different stages of normotensive pregnancy. Acta Obstet Gynecol Scand. 2002, 81 (3): 208-213.View ArticlePubMedGoogle Scholar
- Gross G, Imamura T, Vogt SK, Wozniak DF, Nelson DM, Sadovsky Y, Muglia LJ: Inhibition of cyclooxygenase-2 prevents inflammation-mediated preterm labor in the mouse. Am J Physiol Regul Integr Comp Physiol. 2000, 278 (6): R1415-1423.PubMedGoogle Scholar
- Burns PD, Mendes JO, Yemm RS, Clay CM, Nelson SE, Hayes SH, Silvia WJ: Cellular mechanisms by which oxytocin mediates ovine endometrial prostaglandin F2alpha synthesis: role of G(i) proteins and mitogen-activated protein kinases. Biol Reprod. 2001, 65 (4): 1150-1155. 10.1095/biolreprod65.4.1150.View ArticlePubMedGoogle Scholar
- Vitale A, Perlin J, Leonelli L, Herr J, Wright P, Digilio L, Coonrod S: Mouse cPLA2[gamma], a novel oocyte and early embryo-abundant phospholipase A2 gamma-like protein, is targeted to the nuclear envelope during germinal vesicle breakdown. Developmental Biology. 2005, 282 (2): 374-384. 10.1016/j.ydbio.2005.03.018.View ArticlePubMedGoogle Scholar
- Barrett JC, Fry B, Maller J, Daly MJ: Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005, 21 (2): 263-265. 10.1093/bioinformatics/bth457.View ArticlePubMedGoogle Scholar
- Gabriel SB, Schaffner SF, Nguyen H, Moore JM, Roy J, Blumenstiel B, Higgins J, DeFelice M, Lochner A, Faggart M, et al: The structure of haplotype blocks in the human genome. Science. 2002, 296 (5576): 2225-2229. 10.1126/science.1069424.View ArticlePubMedGoogle Scholar
- Hafler DA, Compston A, Sawcer S, Lander ES, Daly MJ, De Jager PL, de Bakker PI, Gabriel SB, Mirel DB, Ivinson AJ, et al: Risk alleles for multiple sclerosis identified by a genomewide study. N Engl J Med. 2007, 357 (9): 851-862. 10.1056/NEJMoa073493.View ArticlePubMedGoogle Scholar
- Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, Maller J, Sklar P, de Bakker PI, Daly MJ, et al: PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007, 81 (3): 559-575. 10.1086/519795.View ArticlePubMedPubMed CentralGoogle Scholar
- Gao X, Becker LC, Becker DM, Starmer JD, Province MA: Avoiding the high Bonferroni penalty in genome-wide association studies. Genet Epidemiol. 2009, 34 (1): 100-105.Google Scholar
- R Development Core Team: R: A language and environment for statistical computing. 2009, R Foundation for Statistical ComputingGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1755-8794/3/62/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 (<url>http://creativecommons.org/licenses/by/2.0</url>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.