QTL for white spot syndrome virus resistance and the sex-determining locus in the Indian black tiger shrimp (Penaeus monodon)

BMC Genomics

Alphis G Ponniah¹  Thomas Gitterle⁴  Pitchaiyappan Ravichandran¹  Chavali Gopal¹  Sarangapani Jothivel¹  Jayakani Shanmugakarthik¹  Mudagandur S Shekhar¹  Vinaya Kumar Katneni¹  Matthew Baranski²  Gopalapillay Gopikrishna¹  Nicholas A Robinson³ 

[1] Central Institute of Brackishwater Aquaculture, Raja Annamalai Puram, Chennai 600028, Tamil Nadu, India;Nofima, PO Box 210, 1431 Ås, Norway;Flinders University, Sturt Road, Bedford Park, SA 5042, Australia;Akvaforsk Genetics Centre, Sjølseng N-6600, Sunndalsøra, Norway

关键词: Quantitative trait loci;    Sex-linked gene markers;    Disease resistance;    White spot syndrome virus;    Single nucleotide polymorphism;    Penaeus monodon;   

Others  :  1141182
DOI  :  10.1186/1471-2164-15-731

 received in 2014-04-08, accepted in 2014-08-20,  发布年份 2014

【 摘 要 】

Background

Shrimp culture is a fast growing aquaculture sector, but in recent years there has been a shift away from tiger shrimp Penaeus monodon to other species. This is largely due to the susceptibility of P. monodon to white spot syndrome virus disease (Whispovirus sp.) which has impacted production around the world. As female penaeid shrimp grow more rapidly than males, mono-sex production would be advantageous, however little is known about genes controlling or markers associated with sex determination in shrimp. In this study, a mapped set of 3959 transcribed single nucleotide polymorphisms were used to scan the P. monodon genome for loci associated with resistance to white-spot syndrome virus and sex in seven full-sibling tiger shrimp families challenged with white spot syndrome virus.

Results

Linkage groups 2, 3, 5, 6, 17, 18, 19, 22, 27 and 43 were found to contain quantitative trait loci significantly associated with hours of survival after white spot syndrome virus infection (P < 0.05 after Bonferroni correction). Nine QTL were significantly associated with hours of survival. Of the SNPs mapping to these and other regions with suggestive associations, many were found to occur in transcripts showing homology to genes with putative immune functions of interest, including genes affecting the action of the ubiquitin-proteasome pathway, lymphocyte-cell function, heat shock proteins, the TOLL pathway, protein kinase signal transduction pathways, mRNA binding proteins, lectins and genes affecting the development and differentiation of the immune system (eg. RUNT protein 1A). Several SNPs significantly associated with sex were mapped to linkage group 30, the strongest associations (P < 0.001 after Bonferroni correction) for 3 SNPs located in a 0.8 cM stretch between positions 43.5 and 44.3 cM where the feminisation gene (FEM-1, affecting sexual differentiation in Caenorhabditis elegans) mapped.

Conclusions

The markers for disease resistance and sexual differentiation identified by this study could be useful for marker assisted selection to improve resistance to WSSV and for identifying homogametic female individuals for mono-sex (all female) production. The genes with putative functions affecting immunity and sexual differentiation that were found to closely map to these loci provide leads about the mechanisms affecting these important economic traits in shrimp.

【 授权许可】

   
2014 Robinson et al.; licensee BioMed Central Ltd.

【 图 表 】

Figure 1.

Figure 2.

Figure 3.

Figure 4.

Figure 5.

【 参考文献 】

• [1]Food and Agriculture Organisation of the United Nations: The State of World Fisheries and Aquaculture. Rome: FAO; 2010.
• [2]Escobedo-Bonilla CM, Alday-Sanz V, Wille M, Sorgeloos P, Pensaert MB, Nauwynck HJ: A review on the morphology, molecular characterization, morphogenesis and pathogenesis of white spot syndrome virus. J Fish Dis 2008, 31(1):1-18.
• [3]Sanchez-Martinez JG, Aguirre-Guzman G, Mejia-Ruiz H: White spot syndrome virus in cultured shrimp: a review. Aquac Res 2007, 38(13):1339-1354.
• [4]Lo CF, Kou GH: Virus-associated white spot syndrome of shrimp in Taiwan: a review. Fish Pathol 1998, 33(4):365-371.
• [5]Gjedrem T: Selection and Breeding Programs in Aquaculture. Netherlands: Springer; 2005.
• [6]Cock J, Gitterle T, Salazar M, Rye M: Breeding for disease resistance of Penaeid shrimps. Aquaculture 2009, 286(1–2):1-11.
• [7]Gitterle T: Genetic Analyses for Resistance to White Spot Syndrome Virus (WSSV), Harvest Body Weight and Pond Survival in the Pacific White Shrimp Penaeus (Litopenaeus) Vannamei. (Dr. Scientiarum Thesis 2005:9). Aas Norway: Agricultural University of Norway, Department of Animal and Aquacultural Sciences; 2005.
• [8]Gitterle T, Gjerde B, Cock J, Salazar M, Rye M, Vidal O, Lozano C, Erazo C, Salte R: Optimization of experimental infection protocols for the estimation of genetic parameters of resistance to White Spot Syndrome Virus (WSSV) in Penaeus (Litopenaeus) vannamei. Aquaculture 2006, 261(2):501-509.
• [9]Gitterle T, Ødegård J, Gjerde B, Rye M, Salte R: Genetic parameters and accuracy of selection for resistance to White Spot Syndrome Virus (WSSV) in Penaeus (Litopenaeus) vannamei using different statistical models. Aquaculture 2006, 251(2/4):210-218.
• [10]Gitterle T, Rye M, Salte R, Cock J, Johansen H, Lozano C, Suarez JA, Gjerde B: Genetic (co)variation in harvest body weight and survival in Penaeus (Litopenaeus) vannamei under standard commercial conditions. Aquaculture 2005, 243(1/4):83-92.
• [11]Gitterle T, Salte R, Gjerde B, Cock J, Johansen H, Salazar M, Lozano C, Rye M: Genetic (co)variation in resistance to White Spot Syndrome Virus (WSSV) and harvest weight in Penaeus (Litopenaeus) vannamei. Aquaculture 2005, 246(1/4):139-149.
• [12]Hayes BJ, Gitterle T, Gopikrishna G, Gopal C, Krishna G, Jahageerdar S, Lozano C, Alavandi S, Paulpandi S, Ravichandran P, Rye M: Limited evidence for genetic variation for resistance to the white spot syndrome virus in Indian populations of Penaeus monodon. Aquac Res 2010, 41(11):872-877.
• [13]Luo K, Kong J, Luan S, Meng X-H, Zhang T-S, Wang Q-Y: Effect of inbreeding on survival, WSSV tolerance and growth at the postlarval stage of experimental full-sibling inbred populations of the Chinese shrimp Fenneropenaeus chinensis. Aquaculture 2014, 420–421:32-38.
• [14]Arala-Chaves M, Sequeira T: Is there any kind of adaptive immunity in invertebrates? Aquaculture 2000, 191:247-258.
• [15]Fagutao FF, Yasuike M, Caipang CM, Kondo H, Hirono I, Takahashi Y, Aoki T: Gene expression profile of hemocytes of kuruma shrimp, Marsupenaeus japonicus following peptidoglycan stimulation. Mar Biotechnol 2008, 10(6):731-740.
• [16]Swapna AP, Rosamma P, Valsamma J, Singh ISB: Anti-lipopolysaccharide factor and Crustin-III, the anti-white spot virus peptides in Penaeus monodon: control of viral infection by up-regulation. Aquaculture 2011, 319:11-17.
• [17]Garcia JC, Reyes A, Salazar M, Granja CB: Differential gene expression in White Spot Syndrome Virus (WSSV)-infected naive and previously challenged Pacific white shrimp Penaeus (Litopenaeus) vannamei. Aquaculture 2009, 289(3–4):253-258.
• [18]Jeswin J, Anju A, Thomas PC, Paulton MP, Vijayan KK: Survivability of Penaeus monodon during white spot syndrome virus infection and its correlation with immune related genes. Aquaculture 2013, 380:84-90.
• [19]Woramongkolchai N, Supungul P, Tassanakajon A: The possible role of penaeidin5 from the black tiger shrimp, Penaeus monodon, in protection against viral infection. Dev Comp Immunol 2011, 35(5):530-536.
• [20]Huang PY, Kang ST, Chen WY, Hsu TC, Lo CF, Li KF, Chen LL: Identification of the small heat shock protein, HSP21, of shrimp Penaeus monodon and the gene expression of HSP21 is inactivated after white spot syndrome virus (WSSV) infection. Fish Shellfish Immunol 2008, 25(3):250-257.
• [21]Mai WJ, Wang WN: Protection of blue shrimp (Litopenaeus stylirostris) against the White Spot Syndrome Virus (WSSV) when injected with shrimp lysozyme. Fish Shellfish Immunol 2010, 28(4):727-733.
• [22]Youtong W, Deachamag P, Phongdara A, Chotigeat W: WSSV: VP26 binding protein and its biological activity. Fish Shellfish Immunol 2011, 30(1):77-83.
• [23]Tonganunt M, Nupan B, Saengsakda M, Suklour S, Wanna W, Senapin S, Chotigeat W, Phongdara A: The role of Pm-fortilin in protecting shrimp from white spot syndrome virus (WSSV) infection. Fish Shellfish Immunol 2008, 25(5):633-637.
• [24]Graidist P, Fujise K, Wanna W, Sritunyalucksana K, Phongdara A: Establishing a role for shrimp fortilin in preventing cell death. Aquaculture 2006, 255(1/4):157-164.
• [25]Ruan Y-H, Kuo C-M, Lo C-F, Lee M-H, Lian J-L, Hsieh S-L: Ferritin administration effectively enhances immunity, physiological responses, and survival of Pacific white shrimp (Litopenaeus vannamei) challenged with white spot syndrome virus. Fish Shellfish Immunol 2010, 28(4):542-548.
• [26]Sarathi M, Simon MC, Venkatesan C, Hameed ASS: Oral administration of bacterially expressed VP28dsRNA to protect Penaeus monodon from white spot syndrome virus. Mar Biotechnol 2008, 10(3):242-249.
• [27]Kumar SR, Ahamed VPI, Sarathi M, Basha AN, Hameed ASS: Immunological responses of Penaeus monodon to DNA vaccine and its efficacy to protect shrimp against white spot syndrome virus (WSSV). Fish Shellfish Immunol 2008, 24(4):467-478.
• [28]Vaseeharan B, Anand TP, Murugan T, Chen JC: Shrimp vaccination trials with the VP292 protein of white spot syndrome virus. Lett Appl Microbiol 2006, 43(2):137-142.
• [29]Sun YM, Li FH, Chi YH, Xiang JH: Enhanced resistance of marine shrimp Exopalamon carincauda Holthuis to WSSV by injecting live VP28-recombinant bacteria. Acta Oceanol Sin 2013, 32(2):52-58.
• [30]Rodriguez J, Espinosa Y, Echeverria F, Cardenas G, Roman R, Stern S: Exposure to probiotics and beta-1,3/1,6-glucans in larviculture modifies the immune response of Penaeus vannamei juveniles and both the survival to White Spot Syndrome Virus challenge and pond culture. Aquaculture 2007, 273(4):405-415.
• [31]Amar EC, Faisan JP: Induction of Immunity and Resistance to White Spot Syndrome Virus (WSSV) in Shrimp Penaeus monodon (Fabricius) by Synthetic Oligodeoxynucleotide and Bacterial DNA. Philipp Agric Sci 2012, 95(3):267-277.
• [32]Robalino J, Browdy CL, Prior S, Metz A, Parnell P, Gross P, Warr G: Induction of antiviral immunity by double-stranded RNA in a marine invertebrate. J Virol 2004, 78(19):10442-10448.
• [33]Luo T, Zhang X, Shao Z, Xu X: PmAV, a novel gene involved in virus resistance of shrimp Penaeus monodon. FEBS Lett 2003, 551(1/3):53-57.
• [34]Moss DR, Henning OL, Moss SM: Sexual growth dimorphism in penaeid shrimp. Potential for all female culture? Global Aquaculture Advocate 2002, 5:60-61.
• [35]Zhang LS, Yang CJ, Zhang Y, Li L, Zhang XM, Zhang QL, Xiang JH: A genetic linkage map of Pacific white shrimp (Litopenaeus vannamei): sex-linked microsatellite markers and high recombination rates. Genetica 2007, 131(1):37-49.
• [36]Benzie JAH, Kenway M, Ballment E: Growth of Penaeus monodon X Penaeus esculentus tiger prawn hybrids relative to the parental species. Aquaculture 2001, 193(3–4):227-237.
• [37]Li YT, Byrne K, Miggiano E, Whan V, Moore S, Keys S, Crocos P, Preston N, Lehnert S: Genetic mapping of the kuruma prawn Penaeus japonicus using AFLP markers. Aquaculture 2003, 219(1–4):143-156.
• [38]Staelens J, Rombaut D, Vercauteren I, Argue B, Benzie J, Vuylsteke M: High-density linkage maps and sex-linked markers for the black tiger shrimp (Penaeus monodon). Genetics 2008, 179(2):917-925.
• [39]Baranski M, Gopikrishna G, Robinson N, Vinaya Kumar K, Shekhar MS, Shanmugakarthik J, Jothivel S, Gopal C, Ravichandran P, Kent M, Arnyasi M, Ponniah AG: The development of a high density linkage map for black tiger shrimp (Penaeus monodon) based on cSNPs. PLoS One 2014, 9(1):e85413.
• [40]Kemper KE, Daetwyler HD, Visscher PM, Goddard ME: Comparing linkage and association analyses in sheep points to a better way of doing GWAS. Genet Res 2012, 94(4):191-203.
• [41]Risch N, Merikangas K: The future of genetic studies of complex human diseases. Science 1996, 273:1516-1517.
• [42]Amin N, van Duijn CM, Aulchenko YS: A genomic back-ground based method for association analysis in related individuals. PLoS One 2007, 2:e1274.
• [43]Danwattananusorn T, Fagutao FF, Shitara A, Kondo H, Aoki T, Nozaki R, Hirono I: Molecular characterization and expression analysis of heat shock proteins 40, 70 and 90 from kuruma shrimp Marsupenaeus japonicus. Fish Sci 2011, 77(6):929-937.
• [44]Zeng Y, Lu CP: Identification of differentially expressed genes in haemocytes of the crayfish (Procambarus clarkii) infected with white spot syndrome virus by suppression subtractive hybridization and cDNA microarrays. Fish Shellfish Immunol 2009, 26(4):646-650.
• [45]Shi H, Yan XF, Xu X, Ruan LW: Molecular cloning and characterization of a cDNA encoding extracellular signal-regulated kinase from Litopenaeus vannamei. Fish Shellfish Immunol 2012, 33(4):813-820.
• [46]Jeena K, Prasad KP, Pathan MK, Bab PG: Expression Profiling of WSSV ORF 199 and Shrimp Ubiquitin Conjugating Enzyme in WSSV Infected Penaeus monodon. Asian Australas J Anim Sci 2012, 25(8):1184-1189.
• [47]He F, Fenner BJ, Godwin AK, Kwang J: White spot syndrome virus open reading frame 222 encodes a viral E3 ligase and mediates degradation of a host tumor suppressor via ubiquitination. J Virol 2006, 80(8):3884-3892.
• [48]Mallery DL, McEwan WA, Bidgood SR, Towers GJ, Johnson CM, James LC: Antibodies mediate intracellular immunity through tripartite motif-containing 21 (TRIM21). Proc Natl Acad Sci 2010, 107:19985-19990.
• [49]He F, Kwang J: Identification and characterization of a new E3 ubiquitin ligase in white spot syndrome virus involved in virus latency. Virol J 2009, 5:151-151.
• [50]He F, Syed SM, Hameed ASS, Kwang J: Viral ubiquitin ligase WSSV222 is required for efficient white spot syndrome virus replication in shrimp. J Gen Virol 2009, 90:1483-1490.
• [51]Li F, Li MY, Ke W, Ji YC, Bian XF, Yan XM: Identification of the immediate-early genes of white spot syndrome virus. Virology 2009, 385(1):267-274.
• [52]Wang ZL, Chua HK, Gusti A, He F, Fenner B, Manopo I, Wang H, Kwang J: RING-H2 protein WSSV249 from white spot syndrome virus sequesters a shrimp ubiquitin-conjugating enzyme, PvUbc, for viral pathogenesis. J Virol 2005, 79(14):8764-8772.
• [53]Chen AJ, Wang SA, Zhao XF, Yu XQ, Wang JX: Enzyme E2 from Chinese white shrimp inhibits replication of white spot syndrome virus and ubiquitinates its RING domain proteins. J Virol 2011, 85(16):8069-8079.
• [54]Wang B, Li FH, Dong B, Zhang XJ, Zhang CS, Xiang JH: Discovery of the genes in response to white spot syndrome virus (WSSV) infection in Fenneropenaeus chinensis through cDNA microarray. Mar Biotechnol 2006, 8(5):491-500.
• [55]Roberts RJ, Agius C, Saliba C, Bossier P, Sung YY: Heat shock proteins (chaperones) in fish and shellfish and their potential role in relation to fish health: a review. J Fish Dis 2010, 33(10):789-801.
• [56]Moseley P: Stress proteins and the immune response. Immunopharmacology 2000, 48(3):299-302.
• [57]Shaulian E, Karin M: AP-1 in cell proliferation and survival. Oncogene 2001, 20(19):2390-2400.
• [58]Hsu HY, Chiu SL, Wen MH, Chen KY, Hua KF: Ligands of macrophage scavenger receptor induce cytokine expression via differential modulation of protein kinase signaling pathways. J Biol Chem 2001, 276(31):28719-28730.
• [59]Hsu HY, Hua KF, Lin CC, Lin CH, Hsu J, Wong CH: Extract of Reishi polysaccharides induces cytokine expression via TLR4-modulated protein kinase signaling pathways. J Immunol 2004, 173(10):5989-5999.
• [60]Hsu HY, Wen MH: Lipopolysaccharide-mediated reactive oxygen species and signal transduction in the regulation of interleukin-1 gene expression. J Biol Chem 2002, 277(25):22131-22139.
• [61]Dinarello CA: Biologic basis for interleukin-1 in disease. Blood 1996, 87(6):2095-2147.
• [62]Stohr N, Lederer M, Reinke C, Meyer S, Hatzfeld M, Singer RH, Huttelmaier S: ZBP1 regulates mRNA stability during cellular stress. J Cell Biol 2006, 175(4):527-534.
• [63]Marques MRF, Barracco MA: Lectins, as non-self-recognition factors, in crustaceans. Aquaculture 2000, 191(1–3):23-44.
• [64]Zhao Z-Y, Yin Z-X, Weng S-P, Guan H-J, Li S-D, Xing K, Chan S-M, He J-G: Profiling of differentially expressed genes in hepatopancreas of white spot syndrome virus-resistant shrimp (Litopenaeus vannamei) by suppression subtractive hybridisation. Fish Shellfish Immunol 2007, 22(5):520-534.
• [65]He NH, Qin QW, Xu X: Differential profile of genes expressed in hemocytes of White Spot Syndrome Virus-resistant shrimp (Penaeus japonicus) by combining suppression subtractive hybridization and differential hybridization. Antivir Res 2005, 66(1):39-45.
• [66]Pan D, He NH, Yang ZY, Liu HP, Xu X: Differential gene expression profile in hepatopancreas of WSSV-resistant shrimp (Penaeus japonicus) by suppression subtractive hybridization. Dev Comp Immunol 2005, 29(2):103-112.
• [67]Soderhall I, Bangyeekhun E, Mayo S, Soderhall K: Hemocyte production and maturation in an invertebrate animal; proliferation and gene expression in hematopoietic stem cells of Pacifastacus leniusculus. Dev Comp Immunol 2003, 27(8):661-672.
• [68]Puig-Kroger A, Aguilera-Montilla N, Martinez-Nunez R, Dominguez-Soto A, Sanchez-Cabo F, Martin-Gayo E, Zaballos A, Toribio ML, Groner Y, Ito Y, Dopazo A, Corcuera MT, Martin MJA, Vega MA, Corbi AL: The novel RUNX3/p33 isoform is induced upon monocyte-derived dendritic cell maturation and downregulates IL-8 expression. Immunobiology 2010, 215(9–10):812-820.
• [69]Alarcon-Riquelme ME: A RUNX trio with a taste for autoimmunity. Nat Genet 2003, 35(4):299-300.
• [70]Hernroth B, Baden SP, Holm K, Andre T, Soderhall I: Manganese induced immune suppression of the lobster. Nephrops norvegicus. Aquat Toxicol 2004, 70(3):223-231.
• [71]Hodgkin J: Exploring the envelope: systematic alteration in the sex-determination system of the nematode Caeraorhabditis elegans. Genetics 2002, 162(2):767-780.
• [72]Jager S, Schwartz HT, Horvitz HR, Conradt B: The Caenorhabditis elegans F-box orotein SEL-10 promotes female development and may target FEM-1 and FEM-3 for degradation by the proteasome. Proc Natl Acad Sci U S A 2004, 101(34):12549-12554.
• [73]Kuwabara PE: A complex solution to a sexual dilemma. Dev Cell 2007, 13(1):6-8.
• [74]Starostina NG, Lim JM, Schvarzstein M, Wells L, Spence AM, Kipreos ET: A CUL-2 ubiquitin ligase containing three FEM proteins degrades TRA-1 to regulate C-elegans sex determination. Dev Cell 2007, 13(1):127-139.
• [75]Johnson CL, Spence AM: Epigenetic licensing of germline gene expression by maternal RNA in C. elegans. Science 2011, 333(6047):1311-1314.
• [76]Meuwissen TH, Hayes BJ, Goddard ME: Prediction of total genetic value using genome-wide dense marker maps. Genetics 2001, 157(4):1819-1829.
• [77]Bovenhuis H, Spelman RJ: Selective genotyping to detect quantitative trait loci for multiple traits in outbred populations. J Dairy Sci 2000, 83(1):173-180.
• [78]Ødegård J, Meuwissen THE: Identity-by-descent genomic selection using selective and sparse genotyping. Genet Sel Evol 2014, 46:3-3. BioMed Central Full Text
• [79]Grobet L, Poncelet D, Royo LJ, Brouwers B, Pirottin D, Michaux C, Menissier F, Zanotti M, Dunner S, Georges M: Molecular definition of an allelic series of mutations disrupting the myostatin function and causing double-muscling in cattle. Mamm Genome 1998, 9(3):210-213.
• [80]Grisart B, Coppieters W, Farnir F, Karim L, Ford C, Berzi P, Cambisano N, Mni M, Reid S, Simon P, Spelman R, Georges M, Snell R: Positional candidate cloning of a QTL in dairy cattle: identification of a missense mutation in the bovine DGAT1 gene with major effect on milk yield and composition. Genome Res 2002, 12(2):222-231.
• [81]Grisart B, Farnir F, Karim L, Cambisano N, Kim JJ, Kvasz A, Mni M, Simon P, Frere JM, Coppieters W, Georges M: Genetic and functional confirmation of the causality of the DGAT1 K232A quantitative trait nucleotide in affecting milk yield and composition. Proc Natl Acad Sci U S A 2004, 101(8):2398-2403.
• [82]Du H, Fu L, Xu Y, Kil Z, Xu Z: Improvement in a simple method for isolating white spot syndrome virus (WSSV) from the crayfish Procambarus clarkia. Aquaculture 2007, 262:532-534.
• [83]Sambrook J, Fritsch EF, Maniatis T: Molecular cloning: a laboratory manual. 2nd edition. New York: Cold Spring Harbor Laboratory Press; 1989.
• [84]Hadfield JD: MCMC methods for multi-response generalised linear mixed models: the MCMCglmm R package. J Stat Softw 2010, 33:1-22.
• [85]Wilson AJ, Reale D, Clements MN, Morrissey MM, Postma E, Walling CA, Kruuk LEB, Nussey DH: An ecologists guide to the animal model. J Anim Ecol 2009, 79:13-26.
• [86]Seaton G, Hernandez J, Grunchec JA, White I, Allen J, De Koning DJ, Wei W, Berry D, Haley C, Knott S: GridQTL: a grid portal for QTL mapping of compute intensive datasets. In Proceedings of the 8th World Congress on Genetics Applied to Livestock Production: 2006. Brazil: Belo Horizonte; 2006.
• [87]Benjamini Y, Hochberg Y: Controlling false discovery rate: a practical and powerful approach to multiple testing. J R Statist Soc B 1995, 57(1):289-300.
• [88]Chen W-M, Abecasis GR: Family-based association tests for genomewide association scans. Am J Hum Genet 2007, 81:913-926.
• [89]Devlin B, Roeder K: Genomic control for association studies. Biometrics 1999, 55:997-1004.
• [90]Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MAR, Bender D, Maller J, Sklar P, de Bakker PIW, Daly MJ, Sham PC: PLINK: a toolset for whole-genome association and population-based linkage analysis. Am J Hum Genet 2007, 81:559-575.