PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
 
Nat Genet. Author manuscript; available in PMC 2010 September 1.
Published in final edited form as:
Published online 2010 January 31. doi:  10.1038/ng.526
PMCID: PMC2835984
NIHMSID: NIHMS167247

Mutations in PNKP cause microcephaly, seizures and defects in DNA repair

Maintenance of DNA integrity is critical for all cell types, but neurons are particularly sensitive to mutations in DNA repair genes, which lead to both abnormal development and neurodegeneration 1. We describe a previously unknown autosomal recessive disease characterized by microcephaly, early-onset, intractable seizures and developmental delay (MCSZ). Using genome-wide linkage analysis in consanguineous families, we mapped the disease locus to chromosome 19q13.33 and identified multiple mutations in PNKP that results in severe neurological disease, whereas a splicing mutation is associated with more moderate symptoms. Surprisingly, while patient cells are sensitive to radiation and other DNA damaging agents, no patient has yet developed cancer or immunodeficiency. Unlike other DNA repair defects that affect humans, PNKP mutations universally cause severe seizures. The neurological abnormalities in MCSZ patients may reflect a role for PNKP in several DNA repair pathways.

We identified the autosomal recessive disorder MCSZ with the following features: microcephaly, infantile onset seizures, developmental delay, and variable behavioral problems especially hyperactivity. MCSZ was observed in multiple pedigrees of Middle Eastern and European origin (Figure 1). The first three pedigrees were Arabic Palestinians living in Jordan and the USA. Three other families with similar manifestations were Arabic (Kingdom of Saudi Arabia), Turkish, and of mixed European ancestry (USA). We later found less severely affected patients from Family 7, which is also American of mixed European heritage (USA). Brain MRIs consistently show microcephaly with preserved brain structures, without apparent neuronal migration or other structural abnormalities, and with no evidence of degeneration (Figure 2). The patients did not develop ataxia or other neurological symptoms. Routine clinical genetic and metabolic screening showed no abnormalities. Despite careful inquiry, MCSZ patients did not have a higher frequency of common or uncommon infections, offering no clinical evidence of immunodeficiency. Cells from one patient showed sensitivity to irradiation in a standard colony survival assay2,3. However, no patient has developed cancer by age 21 and heterozygous carriers have not developed early onset cancer or any sign of immunodeficiency (see Clinical Information and Summary in Supplementary Note).

Figure 1
Pedigrees of MCSZ families
Figure 2
Brain MRIs of MCSZ patients

Genome-wide linkage screens suggested a single common homozygous region on chromosome 19q in all six MCSZ patients in the consanguineous Families 1-3. This same locus showed homozygosity in the two patients from Family 4 and the single patient from Family 5 whose parents were not known to be related. No other regions of linkage were observed including known primary microcephaly loci. Interestingly, all six patients in the three Palestinian families (Families 1-3) were homozygous for the same 3cM/1.5Mb haplotype (between markers D19S879 and D19S907) at chromosome 19q13.33, suggesting a common ancestor (Supplementary Information Figure 1). Further investigation revealed a common town of origin for these families. Linkage analysis of Families 1-4 and 6 generated a combined maximum two point LOD score of 5.60 at D19S867 with θ =0 (Supplementary Information Table 2) and a maximum multipoint LOD score of 7.12 (Supplementary Information Figure 2).

We sequenced 41 genes within the minimal region (see Supplementary Information) and only the PNKP gene contained mutations. Homozygous mutations were present in each of the first five families and compound heterozygous mutations were found in Families 6 and 7 (Figure 3 and Supplementary Information Figure 3 for detailed analyses of all mutations). The three Palestinian Families (Families 1-3) shared a homozygous base pair substitution in exon 11 (975G>A), resulting in the non-conservative amino acid change E326K. The families from the Kingdom of Saudi Arabia (Family 4) and Turkey (Family 5) had the same homozygous 17 bp duplication (17 bp dup) in exon 14 (1250_1266dup), resulting in a frame-shift, T424GfsX48. Family 6 (European) had two heterozygous mutations: the same 17 bp dup in exon 14 as the Saudi and Turkish families and a point mutation in exon 5 (526C>T) resulting in L176F. The moderately affected members of Family 7 (European) displayed compound heterozygosity, carrying the 17 bp dup mutation in exon 14 and a 17 bp deletion in intron 15 that disrupts proper mRNA splicing.

Figure 3
PNKP mutations in MCSZ patients

To confirm that these mutations were pathogenic, we analyzed patient-derived, Epstein Bar Virus (EBV)-transformed, lymphocytes from Families 3 (E326K) and 7 (17 bp dup/intron 15 del). The samples from the severely affected patient (Family 3) and the mildly affected patients (Family 7) had greatly diminished PNKP protein compared to unaffected or heterozygous family members (Figure 3b). RT-PCR analysis showed that the intron 15 deletion in Family 7 disrupts mRNA splicing and causes skipping of exon 15; a barely detectable level of properly spliced mRNA remains (Figure 3c). The low levels of PNKP protein and the indistinguishable phenotype among severely affected individuals from Families 1-6 suggest that all of these mutations impair PNKP severely or completely, whereas patients in Family 7 that exhibit a slightly milder phenotype and carry one non-coding mutation may retain some PNKP activity.

Surprisingly, despite their diverse ancestries (Saudi, Turkish, and mixed European), all chromosomes with the 17 bp dup mutation appear to share the same haplotype for 18 SNPs in or near the PNKP locus (Supplementary Figure 4a). The low frequency of this mutant haplotype (Supplementary Figure 4b), suggests that the 17 bp dup in exon 14 is likely of shared origin and very old in all or most of these families. Nonetheless, the mutation shows an extremely low carrier frequency in the general population, as the 17 bp dup mutation was absent in 1080 Middle Eastern or Caucasian control chromosomes, while all other mutations were absent in all of at least 280 control chromosomes screened.

PNKP has been implicated in repair of both double strand breaks (DSBs) and single strand breaks (SSBs) 4,5, since its phosphatase domain removes 3′ phosphates and the kinase domain phosphorylates 5′ hydroxyl groups, required for DNA ligation 6. PNKP's forkhead domain mediates interaction with the non-homologous end joining (NHEJ) complex via XRCC4, or the SSB and base excision repair (BER) pathways via XRCC1 (Supplemental Figure 6 a and b)4-8. PNKP has been further implicated in the repair pathway disrupted in an ataxic neurodegenerative disease, spinocerebellar ataxia with axonal neuropathy, SCAN1 (TDP1) (see Supplementary Information Figure 6 for additional details). The abnormalities in a clinical irradiation sensitivity test already reflect deficiencies in NHEJ required to repair radiation induced DSBs. Because PNKP potentially plays roles in additional repair pathways, we tested cells in response to other DNA damaging agents.

The response to free radical damage from hydrogen peroxide, predominately requiring BER, and camptothecin, requiring TDP1 via its activity on topoisomerase I, was examined via alkaline comet assay which detects both SSBs and DSBs by quantitating the amount of DNA that moves from a nucleus after electrophoresis. Typical nuclei with various amounts of DNA damage are shown in Figure 4a. The relative amount of DNA damage was determined by measuring the % Tail DNA (DNA that has moved out of the nucleus). The relative amount of DNA damage repair was quantitated as % Tail DNA after toxin exposure at various times of recovery (see Supplementary Figure 5 for all data) and normalized to the maximum and minimum damage for each cell line in each experiment. MCSZ derived cell lines were significantly impaired in their ability to repair hydrogen peroxide induced damage (Figure 4b) and were also delayed in repairing camptothecin induced damage (Figure 4c). However, MCSZ derived cells were eventually able to repair the camptothecin induced damage after 45 minutes while they were not able to repair all hydrogen peroxide induced damage even after 90 min. These data suggest that camptothecin induced damage may be easier to repair than free radical damage in the setting of PNKP mutations. Although MCSZ patients do not develop the ataxia that characterizes TDP1 mutations that also have camptothecin sensitivity 9,10, our patients are perhaps too young (oldest 21 years) or maybe there is a compensatory mechanism for the loss of PNKP.

Figure 4
Patient derived lymphocytes show abnormal DNA repair

Microcephaly can result from failure to produce enough neurons during development (primary microcephaly) or degeneration after normal development. MCSZ patients have no evidence of brain atrophy or clinical regression, making degeneration unlikely (see Supplemental Note, Clinical Information for further details). In situ hybridization indicated that human and mouse PNKP mRNAs are expressed in both dividing neuronal precursors in the cerebral cortical ventricular zone (Figure 5 a and c VZ) as well as post-mitotic neurons of the cortical plate (Figure 5d CP) consistent with potential roles in both dividing and post-mitotic neurons. In addition, we found that when levels of Pnkp were reduced in dissociated mouse neurons in vitro via RNAi, there was a small but statistically significant increase in apoptosis in both neuronal precursors (Supplemental Information Figure 7c) and differentiated neurons (7d) compared to cells transfected with control plasmids. This indicates that microcephaly could result from apoptosis of precursors, differentiated neurons or both cell types.

Figure 5
PNKP in situ hybridization

Although disruption of genes encoding NHEJ repair proteins leads to microcephaly in human and/or mouse in LIG4 syndrome, Lig411-13, severe combined immunodeficiency with microcephaly, NHEJ1 (cernunnos) 14, Xrcc4 15, Xrcc6 (Ku70)/ Xrcc5 (Ku80)16, and Prkdc (DNA-PKcs) 17, disruption of NHEJ function is not known to cause seizures, which are a prominent clinical feature of MCSZ. In contrast, disruption of the BER/SSB pathways by deletion of Xrcc1 in mice produced normal brain size but the mutant mice developed ataxia, loss of interneurons within the cerebellum, and seizure-like behavior 18. The intriguing similarity of phenotypes between patients with mutations in PNKP and mice with targeted deletion of Xrcc1 does suggest a potential common mechanism demonstrating the requirement for the BER/SSB pathways in the prevention of seizures, potentially via an interneuron specific requirement of the BER/SSB pathway. Therefore, the unique pattern of neurological symptoms of MCSZ may reflect a requirement for PNKP activity in multiple DNA repair pathways.

Online Methods

Genetic screening

Family 1 underwent a genome-wide linkage screen using about 400 microsatellite markers in the ABI linkage mapping set MD v2.5 at 10 cM average density (Applied Biosystems) at the Children's Hospital in Boston (Boston, MA, Genotyping Core facility). Genome-wide screens for Families 3 and 4 were performed at the Center for Inherited Disease Research using microsatellite markers also at 10 cM average spacing. Fine mapping was done using polymorphic microsatellite markers from the ABI linkage mapping set HD v2.5 at 5cM average density (Applied Biosystems) along with additional microsatellite markers identified using the UCSC Human Genome Browser 19 and synthesized primers (Sigma-Genosys). Two point and multipoint LOD scores calculated using Allegro 20 assumed recessive inheritance, full penetrance, and a disease allele frequency of 0.0001. All nucleotide numbers are in reference to cDNA where A (of the ATG start site) is +1 except for the intronic deletion which is in reference to the genomic sequence where A of ATG of the translational start site is +1 (all from UCSC genome browser, NCBI Build 36.1).

PNKP Western blot

For PNKP, EBV transformed lymphocytes were grown in DMEM with 15% FCS plus normocin. 2×105 cells were lysed in sodium dodecyl sulfate (SDS) loading buffer at 90°C and fractionated by SDS-polyacrylamide gel electrophoresis (PAGE). Proteins were transferred to Hybond-C Extra Nitrocellulose (GE Healthcare), stained with Ponceau S solution (Sigma), and washed in TBST before immunoblotting with anti- PNK (SK3195) 21 at a dilution of 1 in 1000 in 5% milk. Blots were washed and probed with goat anti-rabbit IgG HRP secondary antibody (DakoCytomation) and detected with ECL detection reagents (GE Healthcare). The blot washed in TBST and re-probed with anti-β-actin monoclonal antibody (Sigma) using rabbit anti-mouse IgG HRP (DakoCytomation) as a secondary antibody.

PNKP RT-PCR

EBV transformed lymphocytes were grown as above. RNA was isolated via RNeasy Mini Kit (Qiagen). 5 μg of total RNA was used for first strand synthesis via oligo dT primers via SuperScript III First-Strand Synthesis SuperMix (Invitrogen Corporation). 1 μl of the product of the RT reaction was used in PCR using primers from Exon 10 to Exon 17.

Comet Assays

Alkaline comet assays were performed via the protocol from CometAssayTM ES unit (Trevigen, Inc.) Briefly, patient derived EBV transformed lymphocytes were grown as described. Cells were obtained prior to treatment or exposed to 100 μM hydrogen peroxide (Sigma) or 10 μM camptothecin (Sigma) for 30 and 60 minutes respectively at 37°C. After exposure, cells were immediately collected for 0 min recovery time point or washed once, resuspended in growth media and incubated at 37°C. Cells were subsequently collected at times described. Cells were embedded in low melt agarose, plated upon microscope slides, lysed, treated with alkaline solution, and slides electrophoresed in alkaline solution at 1 Volt/cm (21 V) with ~300 mAmps for 30 minutes. Slides were washed, dried, and DNA stained with SYBR green. Images of nuclei and tails were taken with a Nikon TE2000-E fluorescent microscope with CCD camera and % Tail DNA determined with CometScore 1.5 software (TriTek, Corp.).

In situ hybridization

Human in situ hybridizations were performed at the Human Developmental Biology Resource (Institute of Human Genetics, International Centre for Life, Newcastle upon Tyne, NE1 3BZ) with probes to anti-sense human PNKP cDNA (+488-1500). Mouse in situs were done as previously published 22, with anti-sense probe from mouse Pnkp cDNA (+386-1566).

Statistics were done within Microsoft Excel software unless otherwise described. The % DNA Repair in Figure 4 was determined by measuring mean of % Tail DNA in 3 separate wells (raw data in Supplementary Information Figure 5) and using combined % Tail DNA at 0 minute recovery (maximum damage) and minimum damage measured as full recovery for each individual cell line. P-values were determined by comparing the baseline adjusted DNA repair level from 3 separate wells using T-Test comparing each MCSZ deficient line to 2 normal control cell lines with two-tailed distribution and homoscedastic test.

PNKP mutation analysis

PNKP protein alignments were performed with MegAlign (Lasergene). Three dimensional structure of mouse Pnkp was examined in MacPyMOL from previously published information 23. The Branch Point analysis was done with Human Splicing Finder version 2.3 (http://www.umd.be/HSF/).

Exon 14 17 bp dup haplotype determination

Sequenom SNP genotypes were determined in patient and controls samples at the Molecular Genetics Core Facility at Children's Hospital Boston (Boston, MA, USA). Haplotypes were determined by Mendelian inheritance patterns for parents and offspring and Phase haplotype determining software24-26 for control samples.

Pnkp RNAi studies

Silencer 1.0 GFP had RNAi oligos targeted to mouse Pnkp mRNA identified via the Broad Institute RNAi Consortium shRNA Library, ligated into EcoRI and ApaI sites (see Supplemental Information (Invitrogen)). Cerebral cortices from E13.5 day old Swiss-Webster mice were dissected and dissociated with Papain Dissociation System (Worthington Biochemical Corp.), cells were washed twice in HBSS buffer and transfected via the Amaxa Nucleofactor 96 well shuttle system (Amaxa Biosystems) with either pSilencer mouse Pnkp-RNAi GFP, pSilencer GFP, with or without pSport-human PNKP (Open Biosystems, Entrez clone ID BC057659). Cells were grown on poly-L-ornithine (Sigma-Aldrich) treated Lab-Tek chamber slides with Permox (Nalge Nunc International) in DMEM with 10% FBS for the Pax6 studies and for the NeuN studies cultures were changed to Neurobasal complete media after 4 hours. Cultures were fixed with 4% paraformaldehyde 24 hours post-transfection. Immunostaining was performed in PBS with 0.2% Triton X-100 and 10% goat serum with chick anti-GFP (Abcam), rabbit anti-cleaved caspase 3 (Asp175) (Cell Signaling Technology, Inc.), mouse anti-Pax6 (The Developmental Studies Hybridoma Bank) or mouse anti-NeuN (Chemicon). Primary antibodies were detected with anti-chicken IgG Alexa 488 (Invitrogen), and donkey anti-mouse IgG Cy3 or donkey anti-rabbit IgG Cy5 (Jackson Immunoresearch Inc.). Images were obtained by looking for GFP+ cells with a Nikon TE2000-E fluorescent microscope and Metamorph imaging software. Cell counts were made from digital images with the counter blinded as to whether the cultures were derived from Pnkp-RNAi or control transfections. P-value in Supplementary Figure 7 were calculations were calculated using G-test of goodness-of fit with 5 degrees of freedom, two-tailed P-value and no Williams correction27.

Mouse Pnkp-RNAi testing

Mo-Pnkp RNAi construct was tested against Mo-Pnkp-DsRed or Hu-PNKP-DsRed both expressed by pCAG-DsRed vector 28. RNAi control vector (empty vector)/mouse (Mo)-Pnkp-DsRed fusion protein, MoPnkp-RNAi/MoPnkp-DsRed, RNAi-control/human (Hu)-PNKP-DsRed fusion, and MoPnkp-RNAi/HuPNKP-DsRed were transfected into NIH-3T3 cells with Lipofectamine 2000 (Invitrogen) per manufacturer protocol. Cells were counted after 24 hours in culture. Random fields were chosen by expression of GFP, the number of high GFP expressing cells were counted followed by the number of Ds-Red fusion protein expressing cells. DsRed expression was nuclear as expected from PNKP localization (data not shown). For Western Blot analysis, NIH-3T3 cells were transfected as above. After 24 hours of culture, cells were lysed via the protocol described in the methods of the main body of the manuscript. After protein quantitation via BCA assay, 25 μg of protein per lane was separated via PAGE. Protein was transferred to membrane, blocked with Odyssey blocking buffer (Li-Cor Biosciences) and probed with Mo anti-Actin (AC-15), Rb anti-GFP (Molecular Probes) and Rb anti-DsRed (Clontech). Primary antibody was detected with anti-mouse labeled with IRDye 700 or anti-rabbit IRDye 800 (Li-Cor Biosciences) and detected with Li-Cor imaging system. In the bar graph protein amounts are normalized with RNAi Control conditions.

Supplementary Material

Acknowledgements

JS was supported by the Victoria and Stuart Quan Fellowship. Research was supported by grants from the National Institute of Neurological Disorders and Stroke to ECG (5K08NS059673-02) and to CAW (NSR01-35129), the Fogarty International Center (R21 NS061772), by the Manton Center for Orphan Disease Research, the Simons Foundation, and by the Dubai Harvard Foundation for Medical Research. ECG was also supported by a K12 Child Health Research Center awarded to Children's Hospital Boston (1 K12 HD052896-01A1) and a Research in Training Award from the Child Neurology Foundation. JR and KWC are supported by Medical Research Council (G0600776 & G0400959). CAW is an Investigator of the Howard Hughes Medical Institute.

We thank Sue Lindsay, and Steve Lisgo from the Medical Research Council-Wellcome Trust Human Developmental Biology Resource (HDBR), Institute of Human Genetics, International Centre for Life, Newcastle upon Tyne, UK for performing human in situ hybridizations.

We are grateful for genotyping services provided by the Center for Inherited Disease Research (CIDR). CIDR is fully funded through a federal contract from the National Institutes of Health to The Johns Hopkins University, contract number HHSN268200782096C and NIH N01-HG-65403. Genotyping at Children's Hospital Boston is supported by the IDDRC (CHB MRDDRC, P30 HD18655).

Footnotes

The genetic study was approved by Beth Israel Deaconess Medical Center and Children's Hospital Boston Institutional Review Boards. Appropriate informed consent was obtained from all involved human subjects. All animal work was approved by Harvard Medical School, Beth Israel Deaconess Medical Center, and Children's Hospital Boston Institutional Animal Care and Use Committees.

ELECTRONIC-DATABASE INFORMATION

UCSC Human Genome Browser: http://genome.ucsc.edu/

References

1. McKinnon PJ. DNA repair deficiency and neurological disease. Nat Rev Neurosci. 2009;10:100–12. [PMC free article] [PubMed]
2. Huo YK, et al. Radiosensitivity of ataxia-telangiectasia, X-linked agammaglobulinemia, and related syndromes using a modified colony survival assay. Cancer Res. 1994;54:2544–7. [PubMed]
3. Sun X, et al. Early diagnosis of ataxia-telangiectasia using radiosensitivity testing. J Pediatr. 2002;140:724–31. [PubMed]
4. Jilani A, et al. Molecular cloning of the human gene, PNKP, encoding a polynucleotide kinase 3′-phosphatase and evidence for its role in repair of DNA strand breaks caused by oxidative damage. J Biol Chem. 1999;274:24176–86. [PubMed]
5. Karimi-Busheri F, et al. Molecular characterization of a human DNA kinase. J Biol Chem. 1999;274:24187–94. [PubMed]
6. Chappell C, Hanakahi LA, Karimi-Busheri F, Weinfeld M, West SC. Involvement of human polynucleotide kinase in double-strand break repair by non-homologous end joining. EMBO J. 2002;21:2827–32. [PubMed]
7. Whitehouse CJ, et al. XRCC1 stimulates human polynucleotide kinase activity at damaged DNA termini and accelerates DNA single-strand break repair. Cell. 2001;104:107–17. [PubMed]
8. Ali AA, Jukes RM, Pearl LH, Oliver AW. Specific recognition of a multiply phosphorylated motif in the DNA repair scaffold XRCC1 by the FHA domain of human PNK. Nucleic Acids Res. 2009;37:1701–12. [PMC free article] [PubMed]
9. Takashima H, et al. Mutation of TDP1, encoding a topoisomerase I–dependent DNA damage repair enzyme, in spinocerebellar ataxia with axonal neuropathy. Nat Genet. 2002;32:267–272. [PubMed]
10. El-Khamisy SF, et al. Defective DNA single-strand break repair in spinocerebellar ataxia with axonal neuropathy-1. Nature. 2005;434:108–13. [PubMed]
11. Frank KM, et al. Late embryonic lethality and impaired V(D)J recombination in mice lacking DNA ligase IV. Nature. 1998;396:173–7. [PubMed]
12. Barnes DE, Stamp G, Rosewell I, Denzel A, Lindahl T. Targeted disruption of the gene encoding DNA ligase IV leads to lethality in embryonic mice. Curr Biol. 1998;8:1395–8. [PubMed]
13. O'Driscoll M, et al. DNA ligase IV mutations identified in patients exhibiting developmental delay and immunodeficiency. Mol Cell. 2001;8:1175–85. [PubMed]
14. Buck D, et al. Cernunnos, a novel nonhomologous end-joining factor, is mutated in human immunodeficiency with microcephaly. Cell. 2006;124:287–99. [PubMed]
15. Gao Y, et al. A critical role for DNA end-joining proteins in both lymphogenesis and neurogenesis. Cell. 1998;95:891–902. [PubMed]
16. Gu Y, et al. Defective embryonic neurogenesis in Ku-deficient but not DNA-dependent protein kinase catalytic subunit-deficient mice. Proc Natl Acad Sci USA. 2000;97:2668–73. [PubMed]
17. Vemuri MC, Schiller E, Naegele JR. Elevated DNA double strand breaks and apoptosis in the CNS of scid mutant mice. Cell Death Differ. 2001;8:245–55. [PubMed]
18. Lee Y, et al. The genesis of cerebellar interneurons and the prevention of neural DNA damage require XRCC1. Nat Neurosci. 2009;12:973–80. [PMC free article] [PubMed]
19. Karolchik D, et al. The UCSC Genome Browser Database: 2008 update. Nucleic Acids Res. 2008;36:D773–9. [PMC free article] [PubMed]
20. Gudbjartsson DF, Jonasson K, Frigge ML, Kong A. Allegro, a new computer program for multipoint linkage analysis. Nat Genet. 2000;25:12–3. [PubMed]
21. Breslin C, Caldecott KW. DNA 3′-phosphatase activity is critical for rapid global rates of single-strand break repair following oxidative stress. Mol Cell Biol. 2009;29:4653–62. [PMC free article] [PubMed]
22. Berger UV, Hediger MA. Differential distribution of the glutamate transporters GLT-1 and GLAST in tanycytes of the third ventricle. J Comp Neurol. 2001;433:101–14. [PubMed]
23. Bernstein NK, et al. The molecular architecture of the mammalian DNA repair enzyme, polynucleotide kinase. Mol Cell. 2005;17:657–70. [PubMed]
24. Stephens M, Smith NJ, Donnelly P. A new statistical method for haplotype reconstruction from population data. Am J Hum Genet. 2001;68:978–89. [PubMed]
25. Stephens M, Donnelly P. A comparison of bayesian methods for haplotype reconstruction from population genotype data. Am J Hum Genet. 2003;73:1162–9. [PubMed]
26. Stephens M, Scheet P. Accounting for decay of linkage disequilibrium in haplotype inference and missing-data imputation. Am J Hum Genet. 2005;76:449–62. [PubMed]
27. McDonald JH. Handbook of biological statistics. Vol. 293. Sparky House Publishing; Baltimore: 2008.
28. Matsuda T, Cepko CL. Electroporation and RNA interference in the rodent retina in vivo and in vitro. Proc Natl Acad Sci USA. 2004;101:16–22. [PubMed]