Evaluation of a target region capture sequencing platform using monogenic diabetes as a study-model
© Gao et al.; licensee BioMed Central Ltd. 2014
Received: 6 June 2013
Accepted: 11 December 2013
Published: 29 January 2014
Monogenic diabetes is a genetic disease often caused by mutations in genes involved in beta-cell function. Correct sub-categorization of the disease is a prerequisite for appropriate treatment and genetic counseling. Target-region capture sequencing is a combination of genomic region enrichment and next generation sequencing which might be used as an efficient way to diagnose various genetic disorders. We aimed to develop a target-region capture sequencing platform to screen 117 selected candidate genes involved in metabolism for mutations and to evaluate its performance using monogenic diabetes as a study-model.
The performance of the assay was evaluated in 70 patients carrying known disease causing mutations previously identified in HNF4A, GCK, HNF1A, HNF1B, INS, or KCNJ11. Target regions with a less than 20-fold sequencing depth were either introns or UTRs. When only considering translated regions, the coverage was 100% with a 50-fold minimum depth. Among the 70 analyzed samples, 63 small size single nucleotide polymorphisms and indels as well as 7 large deletions and duplications were identified as being the pathogenic variants. The mutations identified by the present technique were identical with those previously identified through Sanger sequencing and Multiplex Ligation-dependent Probe Amplification.
We hereby demonstrated that the established platform as an accurate and high-throughput gene testing method which might be useful in the clinical diagnosis of monogenic diabetes.
Since 1992, several genetic subtypes of monogenetic diabetes have been described in which gene mutations result in diabetes, primarily by causing beta-cell dysfunction . Thus, patients previously categorized clinically as having maturity-onset diabetes of the young (MODY), permanent neonatal diabetes mellitus (NDM), or transient NDM, can now be classified by genetic sub grouping. Identification of the genetic subgroups is a prerequisite for appropriate treatment, genetic counseling, and prognostic information.
MODY, defined as dominantly inherited, young-onset (at least one family member with onset before 25 years of age), non-autoimmune and non-insulin dependent diabetes (no insulin treatment needed three or more years after diabetes diagnosis or measurable circulating C-peptide) [1, 2], is the most common form of monogenetic diabetes. Molecular genetic findings have shown that MODY is very heterogeneous. To date, the MODY phenotype has been reported to be linked with mutations within 13 different genes: HNF4A (MODY1), GCK (MODY2), HNF1A (MODY3), PDX1 (MODY4), HNF1B (MODY5), NEUROD1 (MODY6), KLF11 (MODY7), CEL (MODY8), PAX4 (MODY9), INS (MODY10), BLK (MODY11) and very recently ABCC8 (MODY12) and KCNJ11 (MODY13) [3–5]. Approximately 70-80% of clinically diagnosed MODY patients carry mutations in one of these 13 genes, however, the remaining 20-30% are likely to carry disease causing mutations in yet unidentified genes [1, 6–8]. Thus, the term MODY encompasses a group of clinically and genetically heterogeneous forms of beta cell dysfunction and in a few instances also insulin resistance which at a molecular level is defined by mutations in various genes [2, 5].
NDM, a form of monogenic diabetes which is usually diagnosed within the first six months of life, has a rare incidence of one in 100,000 to 400,000 live births. Clinically, NDM can be divided into two subtypes: transient neonatal diabetes mellitus (TNDM) and permanent neonatal diabetes mellitus (PNDM). The most frequent cause of PNDM is due to mutations in KCNJ11 or ABCC8 which encode a subunit and a modulator of the ATP-sensitive potassium (KATP) channels in the pancreatic beta-cell, respectively . Patients with PNDM due to KCNJ11 or ABCC8 mutations may be treated with oral sulfonylurea instead of insulin .
The monogenic subset of the diabetic population is often mis-classified as type 1 diabetes (T1D) or type 2 diabetes (T2D) and thus is not offered optimal treatment regimens . In Europe and USA, approximately 2-5% of T2D and 10% of familial T1D have previously been estimated to be misdiagnosed MODY cases [11, 12], however, the exact incidence and prevalence of monogenic diabetes is not known. According to recent estimates, approximately 92 million Chinese are recognized diabetes patients, making China one of the world’s leading countries with regards to diabetes prevalence . In the latest report, 27.7 million Chinese children and 334 million Chinese adults are estimated to be pre-diabetic or diabetic . With this high diabetes prevalence and the known high prevalence of specific MODY forms in young onset diabetic patients in China [15, 16], the prevalence of monogenic diabetes cases in China is likely to be accordingly high, and the importance of clinical differentiation should thus not be ignored. A high-throughput method with high accuracy to diagnose monogenic forms of diabetes is needed not only for a scientific purpose but also for medical application. Due to its high accuracy, Sanger sequencing is still the gold standard in terms of gene variant detection. However, when it comes to heterogeneous disorders, the method requires accurate clinical assessments and gene selection. If the “wrong gene” is tested, a (false) negative result is obtained and the testing procedure has to be repeated. In this case, the method is time consuming and expensive.
Here, we present the use of a gene testing panel (117 genes) based on a target region capture system coupled with next-generation sequencing (NGS) technology. The performance of the capture probe was tested by use of in-house YanHuang (YH) DNA . The assay performance was validated by investigating six of the known monogenic diabetes genes (HNF4A, GCK, HNF1A, HNF1B, INS and KCNJ11) in 70 Danish patients carrying previously identified causative variants. All 70 causative variants were correctly identified. We therefore believe the established platform to be a high-throughput gene testing method which may be applied in clinical diagnosis of monogenic diabetes.
As a means to test the use of the platform, six YH samples were contributed from the Clinical Laboratory Center at Beijing Genomic Institute (BGI) , Shenzhen. For the validation of the method in relation to detection of known diabetes causing mutations, a total of 70 samples from unrelated Danish MODY and NDM probands collected at Steno Diabetes Center, Gentofte, Denmark, were examined. Prior to the participation in the study, informed consent was obtained from all subjects. The study was approved by the Ethical Committee of Copenhagen and was in accordance with the principles of the Declaration of Helsinki II.
Target region array design
A customized oligonucletide probe was designed to capture whole genome regions (including exons, introns, and untranslated regions (UTRs)) of 117 genes for monogenic diabetes [3–5] and selected candidate genes referred to type 2 diabetes mellitus, metabolism of glucose and oral hypoglycemic agents, and obesity [2, 18, 19] (Additional file 1: Table S1). Hence, the total size of probe set was 4962226 bp. The probes (GenCap™ Enrichment, MyGenostics, USA) were approximately 60 bp in length and with a 5-10 bp overlap in restricted regions.
DNA extraction, target region capture and next-generation sequencing
Methods for DNA extraction, target region capture, and NGS have previously been extensively described  and a brief experimental workflow is included in Additional file 1: Figure S1. In brief, genomic DNA was extracted from peripheral blood lymphocytes by standard procedures using QIAamp DNA Bloodmini kits (Qiagen, Germany). Next, 1 ug genomic DNA was fragmented by Covaris sonicator (Covaris S2, USA) to sizes of 150-300 bp and then purified. The blunt ends of the purified DNA fragments were then repaired, and A-tailing was added. The fragments were ligated overnight using standard Illumina paired-end (PE) adapter. The ligated products were then amplified through 4-cycle polymerase chain reactions (PCRs) using PE primers containing 8 bp index tags. The purified PCR products containing 3 ug DNA were hybridized to the GenCap™ probe (in solution) at 65°C for 22 hours using a PCR machine. The products were bound to a rotator for 1 hour at room temperature using Dynal Myone Streptavidin C1 magnetic beads (Invitrogen, USA), which had been activated beforehand, and the products were then washed with buffer according to the kit manual. The captured DNA libraries were amplified using 15-cycle PCRs, purified, and subsequently eluted in a 30 ul volume and subjected to Agilent 2100 Bioanalyzer and quantitative PCR to estimate the magnitude of enrichment. The final captured DNA libraries were sequenced using the Illumina HiSeq2000 Analyzers as PE 90 bp reads (following the manufacturer’s standard cluster generation and sequencing protocols), providing an average coverage depth for each sample of at least 100-fold.
Data filtering and analysis
Image analysis, error estimation, and base calling were performed using the Illumina pipeline (version 1.3.4) with default parameters. Indexed primers were used to identify the different samples in the primary data. All unqualified reads (defined as reads either polluted by adapter, containing more than 10% nucleotides out of read length, having an average quality of less than 10, or having 50% bases with a quality value less than 5) were removed using a local dynamic programming algorithm. The remaining reads were aligned to the reference human genome (UCSC hg19) using Burrows-Wheeler Alignment Tool (BWA-0.5.9). Next, SNPs and indels were identified using SOAPsnp software 2.0 and SAMtools v1.4 while using the recommended parameters [21, 22].
The accuracy and precision of the targeted region capture sequencing were assessed by comparing SNPs identified in one of the 6 YH samples with those reported in the online YH database (http://yh.genomics.org.cn/). The remaining 5 YH samples were used as inter-control samples.
According to the YH database reference, an in-house pipeline was used to estimate the true positive (TP) ratio, true negative (TN) ratio, false positive (FP) ratio, and false negative (FN) ratio of all SNPs . The accuracy and precision was calculated by (TP + TN)/(TP + FN + TN + FP) and TP/(TP + FN), respectively. In this calculation, only high quality SNPs fulfilling the following criteria were included: 1) SOAPsnp score ≥ 20; 2) Depth ≥ 20; 3) Percentage of reads supporting variation ≥ 28% .
Functional annotation of genetic variants
The variants were functionally annotated using an in-house pipeline as well as the reported frequencies available from public databases (dbSNP 135, HapMap database, 1000 genome variants database, and a local control database) and categorized into either missense, nonsense, splice-site, insertion, deletion, synonymous or noncoding mutations. For all variants, the results were filtered using a quality value of single base sequencing ≥ 20. The variants were filtered to potential mutation candidates through: 1). Functional variants (insertion/deletion: in CDS and splicing region, SNP: nonsense, splice site and missense) and 2). Variants with an allele frequency below 0.01 in either of the public databases mentioned above. To validate the pre-screened samples, variants were initially filtered by six known disease-causing genes (HNF4A, GCK, HNF1A, HNF1B, INS and KCNJ11).
The identification of known pathogenic variants was based on mutations previously reported to cause MODY in the literature, in Locus Specific Mutation Databases (LOVD, http://grenada.lumc.nl/LSDB_list/lsdbs) or previously identified to segregate with diabetes in Danish MODY families. Novel variants considered to be pathogenic were either: 1) stop/frameshift variants; 2) missense mutations positioned in the amino acid conservative region across species; 3) splice-site variations fulfilling the GT-AT rules; or 4) predicted to be possibly damaging or disease-causing by more than two of the bioinformatic programs (SIFT, Sorting Intolerant From Tolerant, http://sift.bii.a-star.edu.sg/; PolyPhen-2, http://genetics.bwh.harvard.edu/pph2/; Mutation Taster, http://www.mutationtaster.org/; BDGP, Berkeley Drosophila Genome Project, http://www.fruitfly.org/seq_tools/splice.html) .
Large deletions/duplications analysis
Formula.1 z-score is calculated for large deletions/duplications analysis. Mean Depth rate-region=; Mean Depth rate-total=, other region means other regions of the samples operated in the same run.
The large deletions and duplications were identified using a predefined cut-off point (±3) of derived z-score of each captured gene region. We used the cut-off value of 3 for absolute z-score, as it represents the 99.9th percentile of the normal samples set for one tailed region. Any region with a z-score above 3 was defined as either a deletion (<-3) or a duplication (> 3).
General performance of the overall sequencing
Statistics of the performance
Target region capture statistics
Target region (bp):
Raw data yield (Mb)
Data mapped to target region (Mb):
Reads mapped to genome:
Reads unique mapped to genome:
Reads mapped to target region:
Mean depth of target region:
Coverage of target region > =1X (%):
Coverage of target region > =20X (%):
Average read length (bp):
Capture specificity (%):
Mean depth of flanking region:
Coverage of 200 bp flanking region (%):
Statistics of the SNP count
Total number of SNPs
Statistics of the indel count
Total number of indels
Accuracy and precision of target region captured sequencing
The accuracy and precision assessment
SOAPsnp score ≥ 20; Depth ≥ 20;
Percentage of reads supporting ≥ 28%
Total number of SNPs:
Total number of target SNPs:
YH target in dbSNP:
True positive (TP ratio):
True negative (TN ratio):
False negative (FN ratio):
False positive (FP ratio):
Accuracy (%)((TP + TN)/(TP + FN + TN + FP)):
Precision (%)(TP/(TP + FN)):
Identified variants in monogenic diabetes samples
To validate whether the method can accurately detect causative variants in patients with known disease causing mutations in HNF4A, GCK, HNF1A, HNF1B, INS, or KCNJ11, we applied the method to 70 samples from Danish patients.
The clinical characteristics of Danish patients of monogenic diabetes mellitus
Age at examination (years) [NA]
Age at diagnosis (years) [NA]
Height (cm) [NA]
Weight (kg) [NA]
BMI (kg/m2) [NA]
Treatment (ins/ins + OHA/OHA/diet) [NA]
23.84 (1.4-73) 
16.26 (1.4-35) 
155.55 (80-176) 
55.53 (12-95.5) 
21.63 (13.46-37.59) 
6.56 (5.9-8.1) 
35.14 (8-60) 
23.24 (8-49) 
171.45 (155-188.5) 
71.48 (50.8-106.6) 
24.26 (18.82-37.77) 
7.42 (5.5-9.4) 
29.5 (26-33) 
29 (26-32) 
167 (162-172) 
58.25 (55.4-61.1) 
21.00 (18.73-23.28) 
6.05 (5.6-6.5) 
22.33 (4-35) 
15.67 (4-34) 
157.19 (96.5-183) 
57.34 (15.4-81.3) 
22.17 (16.54-25.59) 
7.09 (5.5-8.4) 
39.5 (35-44) 
10.6 (0.2-21) 
171.5 (169-174) 
72.2 (71-73.4) 
24.55 (24.24-24.86) 
6.8 (5.9-7.7) 
6.83 (3.5-13) 
0.33 (0.2-0.5) 
131.5 (108-155) 
29.9 (20.8-39) 
17.03 (16.23-17.83) 
7.37 (6.4-8.1) 
28.11 (1.4-73) 
18.38 (0.2-49) 
161.69 (80-188.5) 
61.26 (12-106.6) 
22.58 (13.46-37.77) 
6.96 (5.5-9.4) 
In the present study, target region enrichment combined with NGS technology proved to be a useful tool for the identification of genetic variants involved in monogenic diabetes. The hybridizing procedure to enrich target region is easy to use, and the duration of the total process is approximately 30 hours which is considerably less than performing directional sequencing of each gene. Combined with the applied sample pooling strategy, the target region capture procedure proved to be very effective. With the use of this method, 117 genes potentially implicated in metabolic disorders could be analyzed in one run. Each sample generated an average of 2005.13 Mb data mapped to the target region. Though the capture specificity of the applied chip is relatively low (53.72% compared to ~70% on other commercial chips, e.g. Nimblegen EZ Choice), the compatible performance in capture coverage (reducing the FN rate) as well as the lower cost (20% more sequencing cost is also considered) makes the present capture chip a good choice for target-region capture sequencing. For the entire testing panel, the coverage of the 117 genes was 99.38%, and the coverage was only slightly diminished (96.95%) at 20-fold minimum depth.
To assess the accuracy of the method, in-house YH DNA was initially tested. We found that after filtering of the low quality variants, the accuracy was 95%, and both the FP and FN ratio reached the lowest value. Nevertheless, comparison with the reference data revealed that the FN ratio was still quite high (9.04%). The relatively high FN ratio may be due to errors in the YH reference sequence; such false positives in the reference sequence will cause an inflated FN rate in the examined sequence. This is supported by the fact that an FN ratio of 5-10% has also been observed in other previous studies using the YH reference sequence .
By applying this method, a total of 70 pre-screened samples were tested, and in 100% of the samples, we were able to detect known mutations. With the use of this method, small variants (both SNPs and indels) as well as larger deletions and duplications could be identified in one analysis. In order to apply the method in clinical settings, several improvements could be implemented. Firstly, the gene panel could be re-designed to include only translated regions of phenotype specific genes (known causative genes, candidate genes for MODYX, other disease functional related genes), plus regulatory regions with reported causative mutations, for example the P2 promoter of HNF4A . Secondly, any low coverage regions should be supplemented by Sanger sequencing methods in order to improve the mutation detection rate. Thirdly, our assay at its current form has a turnaround time of seven weeks, which is relatively long for routine clinical testing and reporting. Therefore, transferring the sequencing procedure from Hiseq to Miseq or Proton might shorten turnaround time while maintaining the high accuracy and throughput. Last but not least, clinical variation interpretation can be a great challenge. Thus, we suggest a “phenotype-specific candidate-genes panel” approach i.e. only genes of relevance to the individuals’ phenotype are analysed in order to reduce the number of variants identified to facilitate the interpretation of the pathogenicity of the variants . A similar approach was recently published targeting the coding regions of 29 genes in which mutations have been reported to cause neonatal diabetes, MODY, maternally inherited diabetes and deafness (MIDD) or familial partial lipodystrophy (FPLD) . However, as the present platform includes a much larger number of genes, includes non-coding region, and has the ability to identify larger structural variations, we believe that the present platform has a greater usability.
We have established a straight-forward target region sequencing method for examining a large number of genes. The method was successfully evaluated using monogenic diabetes as a study-model. The turnaround time was seven weeks, and the coverage of the target gene coding region reached 100%. The method can be used to detect small size mutations and larger deletions and/or duplications in the same run. In terms of the 70 pre-screened samples, we correctly identified all the causative variants. We believe our method could be a useful tool in large-scale monogenic diabetes gene testing.
Beijing Genomic Institute
Body mass index
Locus specific mutation databases
Matrix-Assisted Laser Desorption/Ionization time-of-flight mass spectrometry
Multiplex ligation-dependent probe amplification
Maturity-Onset Diabetes of the Young
Neonatal diabetes mellitus
Next generation sequencing
Polymerase chain reaction
Persistent hyperinsulinemic hypoglycemia of infancy
Permanent neonatal diabetes mellitus
Oral hypoglycemic agents
Single nucleotide polymorphism
Type 1 diabetes
Transient neonatal diabetes mellitus
This study was supported by BGI-Shenzhen Research (China) (The Guangdong Enterprise Key Laboratory of Human Disease Genomics; ShenZhen Engineering Laboratory for Clinical molecular diagnostic), Steno Diabetes Center and The Novo Nordisk Foundation Center for Basic Metabolic Research (Denmark). We thank all the staff at BGI and at Metabolic Genetics Section, The Novo Nordisk Foundation Center for Basic Metabolic Research. The Novo Nordisk Foundation Center for Basic Metabolic Research is an independent Research Center at the University of Copenhagen partially funded by an unrestricted donation from the Novo Nordisk Foundation (http://www.metabol.ku.dk).
- Murphy R, Ellard S, Hattersley AT: Clinical implications of a molecular genetic classification of monogenic beta-cell diabetes. Nat Clin Pract Endocrinol Metab. 2008, 4 (4): 200-213. 10.1038/ncpendmet0778.PubMedView ArticleGoogle Scholar
- Fajans SS, Bell GI, Polonsky KS: Molecular mechanisms and clinical pathophysiology of maturity-onset diabetes of the young. N Engl J Med. 2001, 345 (13): 971-980. 10.1056/NEJMra002168.PubMedView ArticleGoogle Scholar
- Flannick J, Beer NL, Bick AG, Agarwala V, Molnes J, Gupta N, Burtt NP, Florez JC, Meigs JB, Taylor H, Lyssenko V, Irgens H, Fox E, Burslem F, Johansson S, Brosnan MJ, Trimmer JK, Newton-Cheh C, Tuomi T, Molven A, Wilson JG, O'Donnell CJ, Kathiresan S, Hirschhorn JN, Njølstad PR, Rolph T, Seidman JG, Gabriel S, Cox DR, Seidman CE, et al: Assessing the phenotypic effects in the general population of rare variants in genes for a dominant Mendelian form of diabetes. Nature Genetics. 2013, 45 (11): 1380-1385. 10.1038/ng.2794.PubMedPubMed CentralView ArticleGoogle Scholar
- Bowman P, Flanagan SE, Edghill EL, Damhuis A, Shepherd MH, Paisey R, Hattersley AT, Ellard S: Heterozygous ABCC8 mutations are a cause of MODY. Diabetologia. 2012, 55 (1): 123-127. 10.1007/s00125-011-2319-x.PubMedView ArticleGoogle Scholar
- Bonnefond A, Philippe J, Durand E, Dechaume A, Huyvaert M, Montagne L, Marre M, Balkau B, Fajardy I, Vambergue A, Vatin V, Delplanque J, Le Guilcher D, De Graeve F, Lecoeur C, Sand O, Vaxillaire M, Froguel P: Whole-exome sequencing and high throughput genotyping identified KCNJ11 as the thirteenth MODY gene. PLoS One. 2012, 7 (6): e37423-10.1371/journal.pone.0037423.PubMedPubMed CentralView ArticleGoogle Scholar
- Craig ME, Hattersley A, Donaghue KC: Definition, epidemiology and classification of diabetes in children and adolescents. Pediatr Diabetes. 2009, 10 (Suppl 12): 3-12.PubMedView ArticleGoogle Scholar
- Johansen A, Ek J, Mortensen HB, Pedersen O, Hansen T: Half of clinically defined maturity-onset diabetes of the young patients in Denmark do not have mutations in HNF4A, GCK, and TCF1. J Clin Endocrinol Metab. 2005, 90 (8): 4607-4614. 10.1210/jc.2005-0196.PubMedView ArticleGoogle Scholar
- McCarthy MI, Hattersley AT: Learning from molecular genetics: novel insights arising from the definition of genes for monogenic and type 2 diabetes. Diabetes. 2008, 57 (11): 2889-2898. 10.2337/db08-0343.PubMedPubMed CentralView ArticleGoogle Scholar
- Gloyn AL, Pearson ER, Antcliff JF, Proks P, Bruining GJ, Slingerland AS, Howard N, Srinivasan S, Silva JM, Molnes J, Edghill EL, Frayling TM, Temple IK, Mackay D, Shield JP, Sumnik Z, van Rhijn A, Wales JK, Clark P, Gorman S, Aisenberg J, Ellard S, Njølstad PR, Ashcroft FM, Hattersley AT: Activating mutations in the gene encoding the ATP-sensitive potassium-channel subunit Kir6.2 and permanent neonatal diabetes. N Engl J Med. 2004, 350 (18): 1838-1849. 10.1056/NEJMoa032922.PubMedView ArticleGoogle Scholar
- Shields BM, Hicks S, Shepherd MH, Colclough K, Hattersley AT, Ellard S: Maturity-onset diabetes of the young (MODY): how many cases are we missing?. Diabetologia. 2010, 53 (12): 2504-2508. 10.1007/s00125-010-1799-4.PubMedView ArticleGoogle Scholar
- Velho G, Robert JJ: Maturity-onset diabetes of the young (MODY): genetic and clinical characteristics. Horm Res. 2002, 57 (Suppl 1): 29-33.PubMedView ArticleGoogle Scholar
- Moller AM, Dalgaard LT, Pociot F, Nerup J, Hansen T, Pedersen O: Mutations in the hepatocyte nuclear factor-1alpha gene in Caucasian families originally classified as having type I diabetes. Diabetologia. 1998, 41 (12): 1528-1531. 10.1007/s001250051101.PubMedView ArticleGoogle Scholar
- Yang SH, Dou KF, Song WJ: Prevalence of diabetes among men and women in China. N Engl J Med. 2010, 362 (25): 2425-2426. author reply 2426PubMedView ArticleGoogle Scholar
- Yan S, Li J, Li S, Zhang B, Du S, Gordon-Larsen P, Adair L, Popkin B: The expanding burden of cardiometabolic risk in China: the China health and nutrition survey. Obes Rev. 2012, 13 (9): 810-821. 10.1111/j.1467-789X.2012.01016.x.PubMedPubMed CentralView ArticleGoogle Scholar
- Xu JY, Chan V, Zhang WY, Wat NM, Lam KS: Mutations in the hepatocyte nuclear factor-1alpha gene in Chinese MODY families: prevalence and functional analysis. Diabetologia. 2002, 45 (5): 744-746. 10.1007/s00125-002-0814-9.PubMedView ArticleGoogle Scholar
- Xu JY, Dan QH, Chan V, Wat NM, Tam S, Tiu SC, Lee KF, Siu SC, Tsang MW, Fung LM, Xu JY, Dan QH, Chan V, Wat NM, Tam S, Tiu SC, Lee KF, Siu SC, Tsang MW, Fung LM, Chan KW, Lam KS: Genetic and clinical characteristics of maturity-onset diabetes of the young in Chinese patients. Eur J Hum Genet. 2005, 13 (4): 422-427. 10.1038/sj.ejhg.5201347.PubMedView ArticleGoogle Scholar
- Wang J, Wang W, Li R, Li Y, Tian G, Goodman L, Fan W, Zhang J, Li J, Guo Y, Feng B, Li H, Lu Y, Fang X, Liang H, Du Z, Li D, Zhao Y, Hu Y, Yang Z, Zheng H, Hellmann I, Inouye M, Pool J, Yi X, Zhao J, Duan J, Zhou Y, Qin J, Ma L, et al: The diploid genome sequence of an Asian individual. Nature. 2008, 456 (7218): 60-65. 10.1038/nature07484.PubMedPubMed CentralView ArticleGoogle Scholar
- Ramachandrappa S, Farooqi IS: Genetic approaches to understanding human obesity. J Clin Invest. 2011, 121 (6): 2080-2086. 10.1172/JCI46044.PubMedPubMed CentralView ArticleGoogle Scholar
- Morris AP, Voight BF, Teslovich TM, Ferreira T, Segre AV, Steinthorsdottir V, Strawbridge RJ, Khan H, Grallert H, Mahajan A, Feng B, Li H, Lu Y, Fang X, Liang H, Du Z, Li D, Zhao Y, Hu Y, Yang Z, Zheng H, Hellmann I, Inouye M, Pool J, Yi X, Zhao J, Duan J, Zhou Y, Qin J: Large-scale association analysis provides insights into the genetic architecture and pathophysiology of type 2 diabetes. Nat Genet. 2012, 44 (9): 981-990. 10.1038/ng.2383.PubMedPubMed CentralView ArticleGoogle Scholar
- Bettegowda C, Agrawal N, Jiao Y, Sausen M, Wood LD, Hruban RH, Rodriguez FJ, Cahill DP, McLendon R, Riggins G, Velculescu VE, Oba-Shinjo SM, Nagahashi Marie SK, Vogelstein B, Bigner D, Yan H, Papadopoulos N, Kinzler KW: Mutations in CIC and FUBP1 contribute to human oligodendroglioma. Science. 2011, 333 (6048): 1453-1455. 10.1126/science.1210557.PubMedPubMed CentralView ArticleGoogle Scholar
- Li R, Li Y, Fang X, Yang H, Wang J, Kristiansen K: SNP detection for massively parallel whole-genome resequencing. Genome Res. 2009, 19 (6): 1124-1132. 10.1101/gr.088013.108.PubMedPubMed CentralView ArticleGoogle Scholar
- Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R: The sequence alignment/map format and SAMtools. Bioinformatics. 2009, 25 (16): 2078-2079. 10.1093/bioinformatics/btp352.PubMedPubMed CentralView ArticleGoogle Scholar
- Wei X, Ju X, Yi X, Zhu Q, Qu N, Liu T, Chen Y, Jiang H, Yang G, Zhen R, Lan Z, Qi M, Wang J, Yang Y, Chu Y, Li X, Guang Y, Huang J: Identification of sequence variants in genetic disease-causing genes using targeted next-generation sequencing. PLoS One. 2011, 6 (12): e29500-10.1371/journal.pone.0029500.PubMedPubMed CentralView ArticleGoogle Scholar
- Berg JS, Evans JP, Leigh MW, Omran H, Bizon C, Mane K, Knowles MR, Weck KE, Zariwala MA: Next generation massively parallel sequencing of targeted exomes to identify genetic mutations in primary ciliary dyskinesia: implications for application to clinical testing. Genet Med. 2011, 13 (3): 218-229. 10.1097/GIM.0b013e318203cff2.PubMedPubMed CentralView ArticleGoogle Scholar
- Audo I, Bujakowska KM, Leveillard T, Mohand-Said S, Lancelot ME, Germain A, Antonio A, Michiels C, Saraiva JP, Letexier M, Sahel J-A, Bhattacharya SS, Zeitz C: Development and application of a next-generation-sequencing (NGS) approach to detect known and novel gene defects underlying retinal diseases. Orphanet J Rare Dis. 2012, 7: 8-10.1186/1750-1172-7-8.PubMedPubMed CentralView ArticleGoogle Scholar
- Hansen SK, Parrizas M, Jensen ML, Pruhova S, Ek J, Boj SF, Johansen A, Maestro MA, Rivera F, Eiberg H, Andel M, Lebl J, Pedersen O, Ferrer J, Hansen T: Genetic evidence that HNF-1alpha-dependent transcriptional control of HNF-4alpha is essential for human pancreatic beta cell function. J Clin Invest. 2002, 110 (6): 827-833. 10.1172/JCI0215085.PubMedPubMed CentralView ArticleGoogle Scholar
- Rossor AM, Polke JM, Houlden H, Reilly MM: Clinical implications of genetic advances in Charcot-Marie-tooth disease. Nat Rev Neurol. 2013, 9 (10): 562-571. 10.1038/nrneurol.2013.179.PubMedView ArticleGoogle Scholar
- Ellard S, Lango Allen H, De Franco E, Flanagan SE, Hysenaj G, Colclough K, Houghton JA, Shepherd M, Hattersley AT, Weedon MN, Caswell R: Improved genetic testing for monogenic diabetes using targeted next-generation sequencing. Diabetologia. 2013, 56 (9): 1958-1963. 10.1007/s00125-013-2962-5.PubMedPubMed CentralView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.