PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of scirepAboutEditorial BoardFor AuthorsScientific Reports
 
Sci Rep. 2017; 7: 7051.
Published online 2017 August 1. doi:  10.1038/s41598-017-07629-3
PMCID: PMC5539152

The phenotypic variability of HK1-associated retinal dystrophy

Abstract

Inherited retinal dystrophies (IRDs) are a clinically and genetically heterogeneous group of Mendelian disorders primarily affecting photoreceptor cells. The same IRD-causing variant may lead to different retinal symptoms, demonstrating pleiotropic phenotype traits influenced by both underlying genetic and environmental factors. In the present study, we identified four unrelated IRD families with the HK1 p.E851K variant, which was previously reported to cause autosomal dominant retinitis pigmentosa (RP), and described their detailed clinical phenotypes. Interestingly, we found that in addition to RP, this particular variant can also cause dominant macular dystrophy and cone-rod dystrophy, which primarily affect cone photoreceptors instead of rods. Our results identified pleiotropic effects for an IRD-causing variant and provide more insights into the involvement of a hexokinase in retinal pathogenesis.

Introduction

Inherited retinal dystrophies (IRDs) are a group of Mendelian diseases caused by photoreceptor dysfunction and cell death. IRDs demonstrate both extensive clinical and genetic heterogeneity. To date, nearly 300 genes have been identified that associate with various forms of IRDs (https://sph.uth.edu/retnet/) (accessed Jun 23rd, 2017) and they are involved in a wide range of biological pathways such as photo-transduction1, 2, ciliogenesis3, 4, pre-mRNA splicing5, 6, transcription regulation7 and protein transport/quality control8, 9. Nevertheless, novel disease-causing genes and pathways are still being identified by whole exome sequencing or whole genome sequencing, adding to the pathological mechanisms of retinal dystrophy.

Different types of IRDs can be caused by mutations in the same gene. For example, mutations in ABCA4 can lead to Stargardt disease, cone-rod dystrophy (CRD), or retinitis pigmentosa (RP), which differentially affects cone and/or rod photoreceptors1012. Furthermore, one specific disease-causing variant may also lead to different clinical presentations of IRD. The PRPH2 p.P210R mutation, which was identified in autosomal dominant disease families, displays a broad phenotypic variability ranging from macular to peripheral retinal dystrophy13, 14. This phenomenon demonstrated the complexity of IRD genotype-phenotype correlations and that modifier alleles or environmental factors may play important roles in determining the phenotype of a patient.

One missense mutation affecting the human hexokinase 1 (HK1, NM_033497.2, c.2551 G > A, p.E851K) was previously identified in multiple families to associate with autosomal dominant retinitis pigmentosa15, 16. Here, by targeted capture sequencing of known IRD-associated genes, we report four additional IRD cases with this variant, further confirming the disease association and supplementing the described phenotypic variability of HK1-related retinal dystrophy.

Results

Genetic findings

A total of four unrelated IRD probands (one Asian and three Caucasians) were found to possess the previously reported HK1 p.E851K variant by target capture panel sequencing. High quality next-generation sequencing data were obtained for each sample (mean coverage: 128.7; coverage  10×: 97.5%) and no other plausible disease-causing variants were identified. Additional rare protein-coding variants identified are summarized in Supplementary Table 1. Sanger sequencing was performed to confirm this variant and genotype-phenotype co-segregation. The pedigrees of these four families are shown in Fig. 1.

Figure 1
The pedigrees of the four IRD families in this study. “M” indicates the HK1 p.E851K mutation. Arrow indicates the proband. Grey coloring in the Family 2 indicates an unknown retinal phenotype.

We then explored the disease-associated haplotype of these four families and the previous published HK1 disease family16. As shown in Supplementary Table 2, the Han Chinese Family 1 in this study possesses a haplotype clearly different from the previously reported Caucasian family outside the 16 kb window between rs749105 and rs14006. This result strongly suggested that the HK1 E851K variant is bona fide the disease-causing variant.

Clinical findings

The clinical features of all the patients are summarized in Table 1. Family 1 (F1) is of the Han Chinese ethnicity. There are four patients in this family presenting a dominant inheritance pattern. The proband F1:II-1 is a 54 year-old female. Since the initial diagnosis of macular dystrophy at the age of 24, she has been followed up for 30 years. Her central visual acuity has progressively worsened to its current state of hand motion OU (oculus uterque, both eyes). Recently, detailed ophthalmologic examinations were performed in all three available affected family members and representative clinical phenotypes are shown in Fig. 2. The proband’s color vision is impaired and her fundus images showed features typical of macular dystrophy features including posterior pole atrophy with scattered pigment deposits throughout the retina of both eyes. In accordance with the ophthalmoscopy result, optic coherence tomography (OCT) revealed the disappearance of the ellipsoid and interdigitation zones in both eyes. Bilateral large central scotomas were present in her 30 degree visual field exams. Electroretinogram (ERG) results showed rod (scotopic 0.01) and rod-cone combined (scotopic 3.0) responses were normal, and cone (photopic 3.0 and 30 Hz flicker) responses were mildly reduced in both eyes. However, different from the macular dystrophy phenotype in the proband, the patient F1:II-2 was diagnosed as CRD and F1:II-3 was diagnosed as RP. F1:II-2 is a 50 year-old female who has complained about poor darkness adaptation for 15 years. The diagnosis of CRD was confirmed by the ERG results with mildly reduced rod and rod-cone combined responses, as well as moderately reduced cone responses. The best corrected vision acuity (BCVA) of patient F1:II-2 is 20/25 for OD (oculus dexter, the right eye) and 20/30 for OS (oculus sinister, the left eye). Patient F1:II-3 is a 48-year-old male. Since his initial onset of RP, he has suffered from night blindness for about 20 years and reduced visual acuity for 4 years. His best corrected vision acuity is 20/30 (OU). Apparent retinal atrophy and bone-spicule pigmentary deposits along and beyond the arcade with a normal fovea were observed in his ophthalmoscopy. OCT images revealed the loss of the ellipsoid and interdigitation zone except for in the foveal area, with increased foveal thickness in both eyes. His visual field test demonstrated a constricted visual field. The ERG results showed bilateral remarkably reduced rod and rod-cone combined responses, and moderately reduced cone responses. The anterior segments of all three patients was normal. As for the deceased patient F1:I-1, an early-onset macular dystrophy was reported.

Table 1
Ophthalmological findings in the patients.
Figure 2
The clinical phenotypes of the patients. The phenotypes of the three affected individuals in Family 1 and the proband in Family 4. OCT, optical coherence tomography; FC, ophthalmoscopy; FAF, fundus autofluorescence; VF, visual field. OD, the right eye; ...

The proband in Family 2 is a 38 year-old Caucasian female. She was diagnosed with RP at the age of 18. Her 63-year-old mother and 34-year-old sister were reported to have “spots” in their retina but the disease phenotype was not confirmed. By age 29, ERG results revealed less than half of normal responses under all testing conditions. At the age of 30, she was aware of progressive constriction of her visual fields and an increasing nyctalopia. Her BCVA was 20/25 OU (OD: −2.50 + 2.00 × 10; OS: −1.75 + 1.25 × 175). Her ocular fundus showed typical features of RP including a mild diffuse optic pallor, mild vascular attenuation throughout, bone-spicule migration and pigmented epithelial atrophy in four quadrants OU. Formal Goldmann visual field tests showed a modest concentric constriction, even with the V4e targets, and large predominantly nasal dense scotomata in each eye. At 34 years, visual field tests demonstrated that a dense circumferential scotoma (III4e) essentially completely surrounds her central 5 degree island in each eye. Her central acuity decreased to 20/30 OU.

The proband in Family 3 is a 37-year-old Caucasian female diagnosed with RP in her late teenage years. She has no known family history of RP. Biomicroscopy showed no abnormality in the anterior segment but notable vitreous syneresis with pigment and cells in each eye, more anteriorly than posteriorly. Her BCVA was stable at 20/30 OD, 20/25 OS, with a mild myopic refractive error: OD, −2.75 + 1.75 × 180; OS, −2.25 + 1.00 × 175. Ophthalmoscopy revealed an atrophic, degenerating retina with apparent vascular attenuation, a mild diffuse pallor of each optic nerve and scattered bone spicules OU.

Family 4 is a dominant RP family of Caucasian ethnicity and the proband is a 14 year-old boy. In the proband, biomicroscopy showed early RP features in the retina and a mild cystoid macular edema in the left eye. Fundus photos revealed very mild vessel attenuation but pallor in the optic nerve of both eyes. Bilateral thinning of the retina and a significant vitreomacular adhesion in the left eye with a cystoid macular edema were striking in his OCT and autofluorescence images. Autofluorescence also showed a hyperfluorescent stressed ring in both eyes (Figure 2). All these observations were in consistent with the diagnosis of RP.

Discussion

HK1 encodes hexokinase 1, an important enzyme in the glycolysis pathway17. Recessive HK1 mutations are known to lead to nonspherocytic hemolytic anemia, in which cases hexokinase deficiency is the primary cause18, 19 and the phenotype is completely different from IRD. Biochemical studies have shown that the HK1 p.E851K variant does not affect the hexokinase activity or stability of the HK1 protein16. There are also >20 individuals with heterozygous HK1 protein-truncating alleles in the ExAC database, suggesting that HK1-associated dominant IRD probably results from either a gain-of-function or dominant-negative mechanism instead of haploinsufficiency. The disease-causing variant p.E851K is extremely rare but still observed once in the ExAC database (minor allele frequency 1/121,278), which only includes presumably “healthy” controls. This is consistent with our observation of incomplete penetrance in Family 3 and in previous studies15, 16. Incomplete penetrance is frequently observed in dominant IRD cases2022. With the availability of genomic variant databases with significantly larger sample size23, we are very likely to encounter dominant IRD-causing variants in these control databases at a very low frequency, like PRPF31 loss-of-function variants and the HK1 variant in the present study. The variant occurrence, its associated population group, and the disease prevalence will be important for estimating the penetrance of certain dominant IRD-causing variants, thus leading to a more accurate molecular diagnosis. In addition, incomplete penetrance in IRD strongly suggests the effect of additional genetic modifiers, as exemplified by a previous report identifying CNOT3 in PRPF31-associated dominant RP24. It is likely that the incomplete penetrance of HK1-assoicated IRD also arises from genetic modifiers, which can be further revealed by the recruitment of additional large HK1-associated IRD families.

The retina features a high energy consumption, with the particular need for a highly regulated metabolic homeostasis25, probably explaining the involvement of metabolic enzymes in retinopathies. The association between glucose metabolism and retinal function has long been studied, and specifically, light versus darkness is known to lead to different rates of glycolysis in the retina26, 27. Interestingly, we observed variable phenotypes caused by the HK1 p.E851K variant including both macular and peripheral retinal dystrophy, suggesting the universality of HK1 in controlling retinal functions. Since HK1 is involved in glycolysis, a ubiquitous biological process rather than a restricted retinal pathway, it is possible that HK1 genetic defects can be more widely modified by both genetic and non-genetic factors, thus resulting in variable clinical outcomes. Here, by perusing the list of additional protein-coding variants identified by capture panel sequencing, we did not identify variants clearly associated with macular dystrophy phenotype, suggesting that the genetic variation in other genes or intronic variants may be responsible for the phenotypic variability.

In summary, our study identified additional IRD cases with the HK1 p.E851K mutation, confirming the disease association of this variant. We also report that this variant can be associated not only RP but also macular dystrophy, thus expanding the HK1-related phenotypic spectrum. The incomplete penetrance, phenotypic variability, and the probable existence of genetic modifiers associated with this particular IRD-causing variant would provide additional information and guidance on the molecular diagnosis, disease counseling, prediction, and management in the future.

Methods

Patient recruitment

The four families investigated in this study were collected from Peking Union Medical College Hospital, Beijing, China or Baylor College of Medicine, Houston, USA. Detailed medical family histories were recorded. Pedigrees of all the patients are drawn based on interview. Detailed ophthalmologic examinations including best visual acuities, color vision tests (D-15 color plates), slit lamp examination, fundus photography, visual field tests, fundus autofluorescence (AF), optical coherence tomography (OCT, Topcon, Tokyo, Japan) (Carl Zeiss, Germany) and Electroretinograms (ERGs, RetiPort ERG System, Roland Consult, Wiesbaden, Germany) were performed on all the patients. ERGs were performed using corneal “ERGjet” contact lens electrodes. The ERG protocol followed the standards published by the International Society for Clinical Electrophysiology of Vision. Genomic DNA was extracted from peripheral blood leukocytes using a commercial kit (QIAamp DNA Blood Midi Kit, Qiagen, Hilden, Germany) according to the manufacturer’s protocol. Written informed consent was obtained from all participating individuals or their guardians. All experimental methods were approved by the Institutional Review Boards of Peking Union Medical College Hospital and Baylor College of Medicine, and were performed in accordance with relative guidelines and regulations.

Next generation sequencing and data analysis

Pre-capture Illumina libraries were generated as in previous literature8, 28, 29. The targeted DNA was captured, washed and recovered using Agilent Hybridization and Wash Kits (Agilent Technologies). The panel sequencing was performed by capturing DNA with probes targeting the exons of known retinal disease-causing genes30. Captured DNA libraries were sequenced on Illumina HiSeq. 2000 (Illumina, Inc.). After sequencing, the reads were aligned to human hg19 genome using BWA version 0.6.131. Base quality recalibration and local realignment was performed by the Genome Analysis Tool Kit version 3.632. Atlas-SNP2 and Atlas-Indel2 were used for variant calling33. Variant frequency data were obtained from public and internal control databases including the Exome Aggregation Consortium (ExAC) database23, CHARGE consortium34, ESP-650035 and 1000 Genome Project36. Variants with a minor allele frequency higher than the set threshold (1/200 for recessive; 1/10,000 for dominant) were filtered out. Synonymous and deep intronic (distance > 10 bp from exon-intron junctions) variants were also excluded from further analysis. ANNOVAR37 (version 06/17/2015) and dbNSFP38 (version 2.9, includes SIFT, PolyPhen-2, LRT, MutationTaster, MutationAssessor, etc.) were used to annotate protein-altering effects. Reported retinal disease-causing variants were detected based on the HGMD professional database (version 08/15/2016).

Electronic supplementary material

Acknowledgements

We would like to thank the probands and their family members for participating in this study. We thank Andrew Smith for editing the manuscript. This work is supported by grants to R.S. from the Ministry of Science and Technology of the People’s Republic of China (Grant No.: 2010DFB33430), the Foundation Fighting Blindness USA (CD-CL-0214-0631-PUMCH), the National Natural Science Foundation of China (81470669), Beijing Natural Science Foundation (7152116) and CAMS Innovation Fund for Medical Sciences (CIFMS 2016-12M-1-002). This work is also supported by the National Eye Institute [R01EY022356, R01EY018571, vision core grant P30EY002520], the Retina Research Foundation, the Foundation Fighting Blindness [grant number BR-GE-0613-0618-BCM] to R.C.

Author Contributions

Author Contributions

M.X. and R.C. conceived the project. Z.Y., E.Y.C., H.L., L.Y., S.W., L.-J.C.W., R.A.L. and R.S. collected clinical data, Y.L. performed next-generation sequencing. Z.Y., B.L., M.X. and Z.T.S. analyzed genetic data. M.X. prepared Figure 1. Z.Y., E.Y.C., and R.S. prepared Figure 2. Z.Y., B.L. and M.X. wrote the manuscript. All authors reviewed and revised the manuscript.

Notes

Competing Interests

The authors declare that they have no competing interests.

Footnotes

Zhisheng Yuan and Baiyu Li contributed equally to this work.

Electronic supplementary material

Supplementary information accompanies this paper at doi:10.1038/s41598-017-07629-3

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Contributor Information

Ruifang Sui, moc.361@iusfrh.

Rui Chen, ude.mcb@nehciur.

References

1. Janecke AR, et al. Mutations in RDH12 encoding a photoreceptor cell retinol dehydrogenase cause childhood-onset severe retinal dystrophy. Nat. Genet. 2004;36:850–854. doi: 10.1038/ng1394. [PubMed] [Cross Ref]
2. Morimura H, et al. Mutations in the RPE65 gene in patients with autosomal recessive retinitis pigmentosa or leber congenital amaurosis. Proc. Natl. Acad. Sci. U. S. A. 1998;95:3088–3093. doi: 10.1073/pnas.95.6.3088. [PubMed] [Cross Ref]
3. Bujakowska KM, et al. Mutations in IFT172 cause isolated retinal degeneration and Bardet-Biedl syndrome. Hum. Mol. Genet. 2015;24:230–242. doi: 10.1093/hmg/ddu441. [PMC free article] [PubMed] [Cross Ref]
4. Xu M, et al. Mutations in human IFT140 cause non-syndromic retinal degeneration. Hum. Genet. 2015;134:1069–1078. doi: 10.1007/s00439-015-1586-x. [PMC free article] [PubMed] [Cross Ref]
5. Mordes D, et al. Pre-mRNA splicing and retinitis pigmentosa. Mol. Vis. 2006;12:1259–1271. [PMC free article] [PubMed]
6. Xu M, et al. Mutations in the Spliceosome Component CWC27 Cause Retinal Degeneration with or without Additional Developmental Anomalies. Am. J. Hum. Genet. 2017;100:592–604. doi: 10.1016/j.ajhg.2017.02.008. [PubMed] [Cross Ref]
7. Wang F, et al. A homozygous missense mutation in NEUROD1 is associated with nonsyndromic autosomal recessive retinitis pigmentosa. Invest. Ophthalmol. Vis. Sci. 2015;56:150–155. doi: 10.1167/iovs.14-15382. [PMC free article] [PubMed] [Cross Ref]
8. Xu M, et al. ATF6 Is Mutated in Early Onset Photoreceptor Degeneration With Macular Involvement. Invest. Ophthalmol. Vis. Sci. 2015;56:3889–3895. doi: 10.1167/iovs.15-16778. [PMC free article] [PubMed] [Cross Ref]
9. Xu M, et al. Mutations in POMGNT1 cause non-syndromic retinitis pigmentosa. Hum. Mol. Genet. 2016;25:1479–1488. doi: 10.1093/hmg/ddw022. [PMC free article] [PubMed] [Cross Ref]
10. Allikmets R, et al. A photoreceptor cell-specific ATP-binding transporter gene (ABCR) is mutated in recessive Stargardt macular dystrophy. Nat. Genet. 1997;15:236–246. doi: 10.1038/ng0397-236. [PubMed] [Cross Ref]
11. Martinez-Mir A, et al. Retinitis pigmentosa caused by a homozygous mutation in the Stargardt disease gene ABCR. Nat. Genet. 1998;18:11–12. doi: 10.1038/ng0198-11. [PubMed] [Cross Ref]
12. Maugeri A, et al. Mutations in the ABCA4 (ABCR) gene are the major cause of autosomal recessive cone-rod dystrophy. Am. J. Hum. Genet. 2000;67:960–966. doi: 10.1086/303079. [PubMed] [Cross Ref]
13. Feist RM, White MF, Jr., Skalka H, Stone EM. Choroidal neovascularization in a patient with adult foveomacular dystrophy and a mutation in the retinal degeneration slow gene (Pro 210 Arg) Am. J. Ophthalmol. 1994;118:259–260. doi: 10.1016/S0002-9394(14)72913-7. [PubMed] [Cross Ref]
14. Gorin MB, et al. A peripherin/retinal degeneration slow mutation (Pro-210-Arg) associated with macular and peripheral retinal degeneration. Ophthalmology. 1995;102:246–255. doi: 10.1016/S0161-6420(95)31029-9. [PubMed] [Cross Ref]
15. Sullivan LS, et al. A dominant mutation in hexokinase 1 (HK1) causes retinitis pigmentosa. Invest. Ophthalmol. Vis. Sci. 2014;55:7147–7158. doi: 10.1167/iovs.14-15419. [PMC free article] [PubMed] [Cross Ref]
16. Wang F, et al. A missense mutation in HK1 leads to autosomal dominant retinitis pigmentosa. Invest. Ophthalmol. Vis. Sci. 2014;55:7159–7164. doi: 10.1167/iovs.14-15520. [PMC free article] [PubMed] [Cross Ref]
17. Aleshin AE, et al. The mechanism of regulation of hexokinase: new insights from the crystal structure of recombinant human brain hexokinase complexed with glucose and glucose-6-phosphate. Structure. 1998;6:39–50. doi: 10.1016/S0969-2126(98)00006-9. [PubMed] [Cross Ref]
18. Bianchi M, Magnani M. Hexokinase mutations that produce nonspherocytic hemolytic anemia. Blood Cells Mol. Dis. 1995;21:2–8. doi: 10.1006/bcmd.1995.0002. [PubMed] [Cross Ref]
19. van Wijk R, Rijksen G, Huizinga EG, Nieuwenhuis HK, van Solinge WW. HK Utrecht: missense mutation in the active site of human hexokinase associated with hexokinase deficiency and severe nonspherocytic hemolytic anemia. Blood. 2003;101:345–347. doi: 10.1182/blood-2002-06-1851. [PubMed] [Cross Ref]
20. Saini S, Robinson PN, Singh JR, Vanita V. A novel 7 bp deletion in PRPF31 associated with autosomal dominant retinitis pigmentosa with incomplete penetrance in an Indian family. Exp. Eye Res. 2012;104:82–88. doi: 10.1016/j.exer.2012.09.010. [PubMed] [Cross Ref]
21. Rose AM, Bhattacharya SS. Variant haploinsufficiency and phenotypic non-penetrance in PRPF31-associated retinitis pigmentosa. Clin. Genet. 2016;90:118–126. doi: 10.1111/cge.12758. [PubMed] [Cross Ref]
22. Moore AT, et al. Autosomal dominant retinitis pigmentosa with apparent incomplete penetrance: a clinical, electrophysiological, psychophysical, and molecular genetic study. Br. J. Ophthalmol. 1993;77:473–479. doi: 10.1136/bjo.77.8.473. [PMC free article] [PubMed] [Cross Ref]
23. Lek M, et al. Analysis of protein-coding genetic variation in 60,706 humans. Nature. 2016;536:285–291. doi: 10.1038/nature19057. [PMC free article] [PubMed] [Cross Ref]
24. Venturini G, Rose AM, Shah AZ, Bhattacharya SS, Rivolta C. CNOT3 is a modifier of PRPF31 mutations in retinitis pigmentosa with incomplete penetrance. PLoS genetics. 2012;8:e1003040. doi: 10.1371/journal.pgen.1003040. [PMC free article] [PubMed] [Cross Ref]
25. Wong-Riley MT. Energy metabolism of the visual system. Eye and brain. 2010;2:99–116. doi: 10.2147/EB.S9078. [PMC free article] [PubMed] [Cross Ref]
26. Wang L, Kondo M, Bill A. Glucose metabolism in cat outer retina. Effects of light and hyperoxia. Invest. Ophthalmol. Vis. Sci. 1997;38:48–55. [PubMed]
27. Wang L, Tornquist P, Bill A. Glucose metabolism in pig outer retina in light and darkness. Acta Physiol. Scand. 1997;160:75–81. doi: 10.1046/j.1365-201X.1997.00131.x. [PubMed] [Cross Ref]
28. Salvo J, et al. Next-generation sequencing and novel variant determination in a cohort of 92 familial exudative vitreoretinopathy patients. Invest. Ophthalmol. Vis. Sci. 2015;56:1937–1946. doi: 10.1167/iovs.14-16065. [PMC free article] [PubMed] [Cross Ref]
29. Tajiguli A, et al. Next-generation sequencing-based molecular diagnosis of 12 inherited retinal disease probands of Uyghur ethnicity. Sci. Rep. 2016;6:21384. doi: 10.1038/srep21384. [PMC free article] [PubMed] [Cross Ref]
30. Xu M, et al. ADIPOR1 Is Mutated in Syndromic Retinitis Pigmentosa. Hum. Mutat. 2016;37:246–249. doi: 10.1002/humu.22940. [PMC free article] [PubMed] [Cross Ref]
31. Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009;25:1754–1760. doi: 10.1093/bioinformatics/btp324. [PMC free article] [PubMed] [Cross Ref]
32. McKenna A, et al. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20:1297–1303. doi: 10.1101/gr.107524.110. [PubMed] [Cross Ref]
33. Challis D, et al. An integrative variant analysis suite for whole exome next-generation sequencing data. BMC Bioinformatics. 2012;13:8. doi: 10.1186/1471-2105-13-8. [PMC free article] [PubMed] [Cross Ref]
34. Psaty BM, et al. Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium: Design of prospective meta-analyses of genome-wide association studies from 5 cohorts. Circ. Cardiovasc. Genet. 2009;2:73–80. doi: 10.1161/CIRCGENETICS.108.829747. [PMC free article] [PubMed] [Cross Ref]
35. Tennessen JA, et al. Evolution and functional impact of rare coding variation from deep sequencing of human exomes. Science. 2012;337:64–69. doi: 10.1126/science.1219240. [PMC free article] [PubMed] [Cross Ref]
36. Abecasis GR, et al. A map of human genome variation from population-scale sequencing. Nature. 2010;467:1061–1073. doi: 10.1038/nature09534. [PMC free article] [PubMed] [Cross Ref]
37. Wang K, Li M, Hakonarson H. ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010;38:e164. doi: 10.1093/nar/gkq603. [PMC free article] [PubMed] [Cross Ref]
38. Liu X, Jian X, Boerwinkle E. dbNSFP: a lightweight database of human nonsynonymous SNPs and their functional predictions. Hum. Mutat. 2011;32:894–899. doi: 10.1002/humu.21517. [PMC free article] [PubMed] [Cross Ref]

Articles from Scientific Reports are provided here courtesy of Nature Publishing Group