The present and future of whole-exome sequencing in studying and treating human reproductive disorders

The causes of recurrent spontaneous abortion (RSA) and fetal malformations are multifactorial and unclear in most cases. Environmental, maternal, and genetic factors have been shown to contribute to these defects. Whole-exome sequencing (WES) is widely used to detect genetic variations associated with human diseases and has recently been successfully applied to unveil genetic causes of unexplained recurrent spontaneous abortion (URSA) and fetal malformations. Here, we review the current discovery and diagnosis strategies to identify the underlying pathogenic mutations of URSA and fetal malformations using WES technology and propose to further develop WES, both to advance our understanding of these diseases and to eventually lead to targeted therapies for reproductive disorders.     

Paired analysis of tumor mutation burden calculated by targeted deep sequencing panel and whole exome sequencing in non-small cell lung cancer

Owing to rapid advancements in NGS (next generation sequen-cing), genomic alteration is now considered an essential pre-dictive biomarkers that impact the treatment decision in many cases of cancer. Among the various predictive biomarkers, tumor mutation burden (TMB) was identified by NGS and was con-sidered to be useful in predicting a clinical response in cancer cases treated by immunotherapy. In this study, we directly com-pared the lab-developed-test (LDT) results by target sequencing panel, K-MASTER panel v3.0 and whole-exome sequencing (WES) to evaluate the concordance of TMB. As an initial step, the reference materials (n = 3) with known TMB status were used as an exploratory test. To validate and evaluate TMB, we used one hundred samples that were acquired from surgically resected tissues of non-small cell lung cancer (NSCLC) patients. The TMB of each sample was tested by using both LDT and WES methods, which extracted the DNA from samples at the same time. In addition, we evaluated the impact of capture re-gion, which might lead to different values of TMB; the evalu-ation of capture region was based on the size of NGS and target sequencing panels. In this pilot study, TMB was evalu-ated by LDT and WES by using duplicated reference samples; the results of TMB showed high concordance rate (R² = 0.887). This was also reflected in clinical samples (n = 100), which showed R² of 0.71. The difference between the coding sequence ratio (3.49%) and the ratio of mutations (4.8%) indicated that the LDT panel identified a relatively higher number of mutations. It was feasible to calculate TMB with LDT panel, which can be useful in clinical practice. Furthermore, a customized approach must be developed for calculating TMB, which differs according to cancer types and specific clinical settings.     

Detection of oncogenic mutations in paired circulating tumor DNA and circulating tumor cells in patients with hepatocellular carcinoma

Background and aimsCirculating tumor cells (CTCs) or circulating tumor DNA (ctDNA) may be used for diagnostic or prognostic purposes in patients with hepatocellular carcinoma (HCC). We aim to determine whether CTCs or ctDNA are suitable to determine oncogenic mutations in HCC patients.MethodsTwenty-six mostly advanced HCC patients were enrolled. 30 mL peripheral blood from each patient was obtained. CellSearch system was used for CTC detection. A sequencing panel covering 14 cancer-relevant genes was used to identify oncogenic mutations. TERT promoter C228T and C250T mutations were determined by droplet digital PCR.ResultsCTCs were detected in 27% (7/26) of subjects but at low numbers (median: 2 cells, range: 1–15 cells) and ctDNA in 77% (20/26) of patients. Mutations in ctDNA were identified in several genes: TERT promoter C228T (77%, 20/26), TP53 (23%, 6/26), CTNNB1 (12%, 3/26), PIK3CA (12%, 3/26) and NRAS (4%, 1/26). The TERT C228T mutation was present in all patients with one or more ctDNA mutations, or detectable CTCs. The TERT C228T and TP53 mutations detected in ctDNA were present at higher levels in matched primary HCC tumor tissue. The maximal variant allele frequency (VAF) of ctDNA was linearly correlated with largest tumor size and AFP level (Log10). CtDNA (or TERT C228T) positivity was associated with macrovascular invasion, and positivity of ctDNA (or TERT C228T) or CTCs (≥ 2) correlated with poor patient survival.ConclusionsOncogenic mutations could be detected in ctDNA from advanced HCC patients. CtDNA analysis may serve as a promising liquid biopsy to identify druggable mutations.     

NBN and XRCC3 genetic variants in childhood acute lymphoblastic leukaemia

Nibrin and DNA repair protein XRCC3 are involved in DNA double-strand break repair. We genotyped seven tagging SNPs in these genes (rs1805794, rs709816; rs1063054; rs7141928, rs1799794, rs861530, rs861539) with the aim to analyse their association with acute lymphoblastic leukaemia (ALL), a disease, that is characterised by elevated genetic instability. Study consisted of 460 paediatric ALL cases and 552 healthy controls. For selection of DNA sequence variants we employed SNP-tagging approach, incorporating the HAPMAP CEU reference panel data.We did not find association of analysed and tagged SNPs and derived haplotypes with the ALL risk thus did not confirm the hypothesis that analysed DNA recombination repair variants account for increased susceptibility to ALL.     

Accurate detection of subclonal single nucleotide variants in whole genome amplified and pooled cancer samples using HaloPlex target enrichment

Background

Target enrichment and resequencing is a widely used approach for identification of cancer genes and genetic variants associated with diseases. Although cost effective compared to whole genome sequencing, analysis of many samples constitutes a significant cost, which could be reduced by pooling samples before capture. Another limitation to the number of cancer samples that can be analyzed is often the amount of available tumor DNA. We evaluated the performance of whole genome amplified DNA and the power to detect subclonal somatic single nucleotide variants in non-indexed pools of cancer samples using the HaloPlex technology for target enrichment and next generation sequencing.

Results

We captured a set of 1528 putative somatic single nucleotide variants and germline SNPs, which were identified by whole genome sequencing, with the HaloPlex technology and sequenced to a depth of 792–1752. We found that the allele fractions of the analyzed variants are well preserved during whole genome amplification and that capture specificity or variant calling is not affected. We detected a large majority of the known single nucleotide variants present uniquely in one sample with allele fractions as low as 0.1 in non-indexed pools of up to ten samples. We also identified and experimentally validated six novel variants in the samples included in the pools.

Conclusion

Our work demonstrates that whole genome amplified DNA can be used for target enrichment equally well as genomic DNA and that accurate variant detection is possible in non-indexed pools of cancer samples. These findings show that analysis of a large number of samples is feasible at low cost, even when only small amounts of DNA is available, and thereby significantly increases the chances of indentifying recurrent mutations in cancer samples.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-14-856) contains supplementary material, which is available to authorized users.     

A whole genome analyses of genetic variants in two Kelantan Malay individuals

BackgroundThe sequencing of two members of the Royal Kelantan Malay family genomes will provide insights on the Kelantan Malay whole genome sequences. The two Kelantan Malay genomes were analyzed for the SNP markers associated with thalassemia and Helicobacter pylori infection. Helicobacter pylori infection was reported to be low prevalence in the north-east as compared to the west coast of the Peninsular Malaysia and beta-thalassemia was known to be one of the most common inherited and genetic disorder in Malaysia.ResultBy combining SNP information from literatures, GWAS study and NCBI ClinVar, 18 unique SNPs were selected for further analysis. From these 18 SNPs, 10 SNPs came from previous study of Helicobacter pylori infection among Malay patients, 6 SNPs were from NCBI ClinVar and 2 SNPs from GWAS studies. The analysis reveals that both Royal Kelantan Malay genomes shared all the 10 SNPs identified by Maran (Single Nucleotide Polymorphims (SNPs) genotypic profiling of Malay patients with and without Helicobacter pylori infection in Kelantan, 2011) and one SNP from GWAS study. In addition, the analysis also reveals that both Royal Kelantan Malay genomes shared 3 SNP markers; HBG1 (rs1061234), HBB (rs1609812) and BCL11A (rs766432) where all three markers were associated with beta-thalassemia.ConclusionsOur findings suggest that the Royal Kelantan Malays carry the SNPs which are associated with protection to Helicobacter pylori infection. In addition they also carry SNPs which are associated with beta-thalassemia. These findings are in line with the findings by other researchers who conducted studies on thalassemia and Helicobacter pylori infection in the non-royal Malay population.     

The prognostic value of circulating tumor DNA in patients with melanoma: A systematic review and meta-analysis

BackgroundCirculating tumor DNA (ctDNA) has been investigated as a potential prognostic biomarker to evaluate the therapeutic efficacy and disease progression in melanoma patients, yet results remain inconclusive. The purpose of this study was to illustrate the prognostic value of ctDNA in melanoma.ObjectivesTo describe the clinical prognostic value of ctDNA for melanoma patients.MethodsSearched for eligible articles from Pubmed, Web of Science and Embase. Pooled hazard ratios (HRs) and 95% confidence intervals (CIs) were calculated to evaluate the association between ctDNA at baseline or during treatment and overall survival (OS) and progression-free survival (PFS).ResultsA total of 9 articles were obtained, involving 617 melanoma patients. The pooled HRs revealed that compared with baseline undetectable ctDNA patients, detectable ctDNA was highly correlated with poor OS (HR 2.91, 95% CI: 2.22–3.82; p < 0.001) and PFS (HR 2.75, 95% CI: 1.98–3.83; p < 0.001). A meta-analysis of these adjusted HRs was performed and confirmed that ctDNA collected at baseline was associated with poorer OS/PFS (OS: HR 3.00, 95% CI 2.19–4.11, p < 0.001/PFS: HR 2.68, 95% CI 1.77–4.06, p < 0.001). During treatment, a significant association was shown between ctDNA and poorer OS/PFS (OS: HR 6.26, 95% CI 2.48–15.80, p < 0.001; PFS: HR 4.93, 95% CI 2.36–10.33, p < 0.001).ConclusionInvestigation and application of ctDNA will improve "liquid biopsy" and play a role in early prediction, monitoring disease progression and precise adjusting treatment strategies in melanoma patients.     

Targeted genotyping by sequencing: a new way to genome profile the cat

Targeted GBS is a recent approach for obtaining an effective characterization for hundreds to thousands of markers. The high throughput of next‐generation sequencing technologies, moreover, allows sample multiplexing. The aims of this study were to (i) define a panel of single nucleotide polymorphisms (SNPs) in the cat, (ii) use GBS for profiling 16 cats, and (iii) evaluate the performance with respect to the inference using standard approaches at different coverage thresholds, thereby providing useful information for designing similar experiments. Probes for sequencing 230 variants were designed based on the Felis_catus_8.0. 8.0 genome. The regions comprised anonymous and non‐anonymous SNPs. Sixteen cat samples were analysed, some of which had already been genotyped in a large group of loci and one having been whole‐genome sequenced in the 99_Lives Cat Genome Sequencing Project. The accuracy of the method was assessed by comparing the GBS results with the genotypes already available. Overall, GBS achieved good performance, with 92–96% correct assignments, depending on the coverage threshold used to define the set of trustable genotypes. Analyses confirmed that (i) the reliability of the inference of each genotype depends on the coverage at that locus and (ii) the fraction of target loci whose genotype can be inferred correctly is a function of the total coverage. GBS proves to be a valid alternative to other methods. Data suggested a depth of less than 11× is required for greater than 95% accuracy. However, sequencing depth must be adapted to the total size of the targets to ensure proper genotype inference.     

利用DNA池和测序技术快速筛查SNPs及估算基因频率

选取产蛋性能具有明显差异的4个鸡品种(莱航鸡、阳山鸡、丝羽乌骨鸡和隐性白洛克鸡)构建品种DNA池,采用测序的方法研究鸡催乳素基因5′侧翼调控区远端序列(1028bp)的多态性,快速筛查到8个可能与产蛋性能相关的SNPs(C-2402T、T-2192C、C-2161G、C-2134G、C-2062G、G-2040A、A-1944G和C-1884A)。进一步利用测序图中SNP等位基因峰高的比值估算各鸡品种等位基因的频率,其中C-2402T、C-2161G、C-1884A和C-2062G、G-2040A位点等位基因频率的估算结果分别被PCR-RFLP、PCR-SSCP所验证,说明测序峰高比值估算等位基因频率的方法具有一定的可行性。     

Circulating tumor cell-derived organoids: Current challenges and promises in medical research and precision medicine

Traditional 2D cell cultures do not accurately recapitulate tumor heterogeneity, and insufficient human cell lines are available. Patient-derived xenograft (PDX) models more closely mimic clinical tumor heterogeneity, but are not useful for high-throughput drug screening. Recently, patient-derived organoid cultures have emerged as a novel technique to fill this critical need. Organoids maintain tumor tissue heterogeneity and drug-resistance responses, and thus are useful for high-throughput drug screening. Among various biological tissues used to produce organoid cultures, circulating tumor cells (CTCs) are promising, due to relative ease of ascertainment. CTC-derived organoids could help to acquire relevant genetic and epigenetic information about tumors in real time, and screen and test promising drugs. This could reduce the need for tissue biopsies, which are painful and may be difficult depending on the tumor location. In this review, we have focused on advances in CTC isolation and organoid culture methods, and their potential applications in disease modeling and precision medicine.     

Ecological genomics of local adaptation in Cornus florida L. by genotyping by sequencing

Discovering local adaptation, its genetic underpinnings, and environmental drivers is important for conserving forest species. Ecological genomic approaches coupled with next‐generation sequencing are useful means to detect local adaptation and uncover its underlying genetic basis in nonmodel species. We report results from a study on flowering dogwood trees (Cornus florida L.) using genotyping by sequencing (GBS). This species is ecologically important to eastern US forests but is severely threatened by fungal diseases. We analyzed subpopulations in divergent ecological habitats within North Carolina to uncover loci under local selection and associated with environmental–functional traits or disease infection. At this scale, we tested the effect of incorporating additional sequencing before scaling for a broader examination of the entire range. To test for biases of GBS, we sequenced two similarly sampled libraries independently from six populations of three ecological habitats. We obtained environmental–functional traits for each subpopulation to identify associations with genotypes via latent factor mixed modeling (LFMM) and gradient forests analysis. To test whether heterogeneity of abiotic pressures resulted in genetic differentiation indicative of local adaptation, we evaluated F_st per locus while accounting for genetic differentiation between coastal subpopulations and Piedmont‐Mountain subpopulations. Of the 54 candidate loci with sufficient evidence of being under selection among both libraries, 28–39 were Arlequin–BayeScan F_st outliers. For LFMM, 45 candidates were associated with climate (of 54), 30 were associated with soil properties, and four were associated with plant health. Reanalysis of combined libraries showed that 42 candidate loci still showed evidence of being under selection. We conclude environment‐driven selection on specific loci has resulted in local adaptation in response to potassium deficiencies, temperature, precipitation, and (to a marginal extent) disease. High allele turnover along ecological gradients further supports the adaptive significance of loci speculated to be under selection.     

Structural analysis of new compound heterozygous variants in PEPD gene identified in a patient with Prolidase Deficiency diagnosed by exome sequencing

Prolidase Deficiency (PD) is an autosomal recessive rare disorder caused by loss or reduction of prolidase enzymatic activity due to variants in the PEPD gene. PD clinical features vary among affected individuals: skin ulcerations, recurrent infections, and developmental delay are common. In this study, we describe a 16-year-old boy with a mild PD phenotype comprising chronic eczema, recurrent infections and elevated IgE. Whole exome sequencing analysis revealed three PEPD variants: c.575T>C p.(Leu192Pro) inherited from the mother, and c.692_694del p.(Tyr231del) and c.1409G>A p.(Arg470His), both inherited from the father. The variant p.(Tyr231del) has been previously characterized by high-resolution X-ray structure analysis as altering protein dynamics/flexibility. In order to study the effects of the other two prolidase variants, we performed site directed mutagenesis purification and crystallization studies. A high-resolution X-ray structure could only be obtained for the p.(Arg470His) variant, which showed no significant structural differences in comparison to WT prolidase. On the other hand, the p.(Leu192Pro) variant led to significant protein destabilization. Hence, we conclude that the maternal p.(Leu192Pro) variant was likely causally associated with the proband´s disease, together with the known pathogenic paternal variant p.(Tyr231del). Our results demonstrated the utility of exome sequencing to perform diagnosis in PD cases with mild phenotype.     

Mobile element insertions are frequent in oesophageal adenocarcinomas and can mislead paired-end sequencing analysis

Background

Mobile elements are active in the human genome, both in the germline and cancers, where they can mutate driver genes.

Results

While analysing whole genome paired-end sequencing of oesophageal adenocarcinomas to find genomic rearrangements, we identified three ways in which new mobile element insertions appear in the data, resembling translocation or insertion junctions: inserts where unique sequence has been transduced by an L1 (Long interspersed element 1) mobile element; novel inserts that are confidently, but often incorrectly, mapped by alignment software to L1s or polyA tracts in the reference sequence; and a combination of these two ways, where different sequences within one insert are mapped to different loci. We identified nine unique sequences that were transduced by neighbouring L1s, both L1s in the reference genome and L1s not present in the reference. Many of the resulting inserts were small fragments that include little or no recognisable mobile element sequence. We found 6 loci in the reference genome to which sequence reads from inserts were frequently mapped, probably erroneously, by alignment software: these were either L1 sequence or particularly long polyA runs. Inserts identified from such apparent rearrangement junctions averaged 16 inserts/tumour, range 0–153 insertions in 43 tumours. However, many inserts would not be detected by mapping the sequences to the reference genome, because they do not include sufficient mappable sequence. To estimate total somatic inserts we searched for polyA sequences that were not present in the matched normal or other normals from the same tumour batch, and were not associated with known polymorphisms. Samples of these candidate inserts were verified by sequencing across them or manual inspection of surrounding reads: at least 85 % were somatic and resembled L1-mediated events, most including L1Hs sequence. Approximately 100 such inserts were detected per tumour on average (range zero to approximately 700).

Conclusions

Somatic mobile elements insertions are abundant in these tumours, with over 75 % of cases having a number of novel inserts detected. The inserts create a variety of problems for the interpretation of paired-end sequencing data.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1685-z) contains supplementary material, which is available to authorized users.     

Distinct functional consequences of MUTYH variants associated with colorectal cancer: Damaged DNA affinity,glycosylase activity and interaction with PCNA and Hus1

MUTYH is a base excision repair (BER) enzyme that prevents mutations in DNA associated with 8-oxoguanine (OG) by catalyzing the removal of adenine from inappropriately formed OG:A base-pairs. Germline mutations in the MUTYH gene are linked to colorectal polyposis and a high risk of colorectal cancer, a syndrome referred to as MUTYH-associated polyposis (MAP). There are over 300 different MUTYH mutations associated with MAP and a large fraction of these gene changes code for missense MUTYH variants. Herein, the adenine glycosylase activity, mismatch recognition properties, and interaction with relevant protein partners of human MUTYH and five MAP variants (R295C, P281L, Q324H, P502L, and R520Q) were examined. P281L MUTYH was found to be severely compromised both in DNA binding and base excision activity, consistent with the location of this variation in the iron-sulfur cluster (FCL) DNA binding motif of MUTYH. Both R295C and R520Q MUTYH were found to have low fractions of active enzyme, compromised affinity for damaged DNA, and reduced rates for adenine excision. In contrast, both Q324H and P502L MUTYH function relatively similarly to WT MUTYH in both binding and glycosylase assays. However, P502L and R520Q exhibited reduced affinity for PCNA (proliferation cell nuclear antigen), consistent with their location in the PCNA-binding motif of MUTYH. Whereas, only Q324H, and not R295C, was found to have reduced affinity for Hus1 of the Rad9–Hus1–Rad1 complex, despite both being localized to the same region implicated for interaction with Hus1. These results underscore the diversity of functional consequences due to MUTYH variants that may impact the progression of MAP.     

Correlation between the number of false positive variants and the quality of results using Ion Torrent PGM™ sequencing to screen BRCA genes

Introduction: Next Generation Sequencing (NGS) is cost-effective and a faster method to study genes, but its protocol is challenging.Objective: To analyze different adjustments to the protocol for screening the BRCA genes using Ion Torrent PGM sequencing and correlate the results with the number of false positive (FP) variants.Material and methods: We conducted a library preparation process and analyzed the number of FP InDels, the library concentration, the number of cycles in the target amplification step, the purity of the nucleic acid, the input, and the number of samples/Ion 314 chips in association with the results obtained by NGS.Results: We carried out 51 reactions and nine adjustments of protocols and observed eight FP InDels in homopolymer regions. No FP Single-Nucleotide Polymorphism variant was observed; 67.5% of protocol variables were jointly associated with the quality of the results obtained (p<0.05). The number of FP InDels decreased when the quality of results increased.Conclusion: The Ion AmpliSeq BRCA1/BRCA2 Community Panel had a better performance using four samples per Ion-314 chip instead of eight and the optimum number of cycles in the amplification step, even when using high-quality DNA, was 23. We observed better results with the manual equalization process and not using the Ion Library Equalizer kit. These adjustments provided a higher coverage of the variants and fewer artifacts (6.7-fold). Laboratories must perform internal validation because FP InDel variants can vary according to the quality of results while the NGS assay should be validated with Sanger.     

DNA Methylation and SNPs in VCX are Correlated with Sex Differences in the Response to Chronic Hepatitis B

The study was conducted to explore the mechanisms of sex differences in the response to chronic hepatitis B(CHB) in terms of DNA methylation, SNP genotype, and gene expression. Genomic DNA was isolated from peripheral blood mononuclear cells(PBMCs) of CHB patients and healthy controls and evaluated using the Human Methylation 450 K Assay. The DNA methylation level at hg37 chromosome(CHR) X: 7810800 was further validated using pyrosequencing.SNP genotypes, VCX m RNA expression of PBMCs, and plasma VCX protein concentration were further examined using SNaPshot, RT-qPCR, and Western blot, respectively. Results showed that a total of 5529 CpG loci were differentially methylated between male and female CHB patients. DNA methylation level and CC+CT frequency at CHR X: 7810800,VCX mRNA expression of PBMCs, and plasma VCX protein concentration were higher in female than in male CHB patients. The CHR X: 7810800 locus was hypermethylated in CHB patients with CC+CT genotypes in comparison with those with the TT genotype. In cases of CC+CT genotypes, VCX mRNA expression was negatively correlated with the DNA methylation level. CHB patients with higher levels of HBV DNA, AST, and GGT or higher GPRI scores exhibited lower VCX expression. In conclusion, SNPs and DNA methylation at the CHR X: 7810800 locus cooperatively regulate VCX expression in CHB. The upregulated VCX expression in female CHB patients might represent a mechanism of protection from more severe liver dysfunction and extensive fibrosis, as observed in male CHB patients.     

Simple and accurate single base resolution analysis of 5-hydroxymethylcytosine by catalytic oxidative bisulfite sequencing using micelle incarcerated oxidants

Oxidation of 5-methylcytosine (5mC) is catalyzed by ten-eleven translocation (TET) enzymes to produce 5-hydroxymethylcytosine (5hmC) and following oxidative products. The oxidized nucleotides were shown to be the intermediates for DNA demethylation, as the nucleotides are removed by base excision repair system initiated by thymine DNA glycosylase. A simple and accurate method to determine initial oxidation product 5hmC at single base resolution in genomic DNA is necessary to understand demethylation mechanism. Recently, we have developed a new catalytic oxidation reaction using micelle-incarcerated oxidants to oxidize 5hmC to form 5-formylcytosine (5fC), and subsequent bisulfite sequencing can determine the positions of 5hmC in DNA. In the present study, we described the optimization of the catalytic oxidative bisulfite sequencing (coBS-seq), and its application to the analysis of 5hmC in genomic DNA at single base resolution in a quantitative manner. As the oxidation step showed quite low damage on genomic DNA, the method allows us to down scale the sample to be analyzed.     

Identification of high utility SNPs for population assignment and traceability purposes in the pig using high-throughput sequencing

The objectives of this study were to develop breed-specific single nucleotide polymorphisms (SNPs) in five pig breeds sequenced with Illumina's Genome Analyzer and to investigate their usefulness for breed assignment purposes. DNA pools were prepared for Duroc, Landrace, Large White, Pietrain and Wild Boar. The total number of animals used for sequencing was 153. SNP discovery was performed by aligning the filtered reads against Build 7 of the pig genome. A total of 313,964 high confidence SNPs were identified and analysed for the presence of breed-specific SNPs (defined in this context as SNPs for which one of the alleles was detected in only one breed). There were 29,146 putative breed-specific SNPs identified, of which 4441 were included in the PorcineSNP60 beadchip. Upon re-examining the genotypes obtained using the beadchip, 193 SNPs were confirmed as being breed specific. These 193 SNPs were subsequently used to assign an additional 490 individuals from the same breeds, using the sequenced individuals as reference populations. In total, four breed assignment tests were performed. Results showed that for all methods tested 99% of the animals were correctly assigned, with an average probability of assignment of at least 99.2%, indicating the high utility of breed-specific markers for breed assignment and traceability. This study provides a blueprint for the way next-generation sequencing technologies can be used for the identification of breed-specific SNPs, as well as evidence that these SNPs may be a powerful tool for breed assignment and traceability of animal products to their breeds of origin.     

Mutations in glucokinase and other genes detected in neonatal and type 1B diabetes patient using whole exome sequencing may lead to disease-causing changes in protein activity

Monogenic diabetes is caused by mutations that reduce β-cell function. While Sanger sequencing is the standard method used to detect mutated genes. Next-generation sequencing techniques, such as whole exome sequencing (WES), can be used to find multiple gene mutations in one assay. We used WES to detect genetic mutations in both permanent neonatal (PND) and type 1B diabetes (T1BD).A total of five PND and nine T1BD patients were enrolled in this study. WES variants were assessed using VarioWatch, excluding those identified previously. Sanger sequencing was used to confirm the mutations, and their pathogenicity was established via the literature or bioinformatic/functional analysis. The PND and T1BD patients were diagnosed at 0.1–0.5 and 0.8–2.7?years of age, respectively. Diabetic ketoacidosis was present at diagnosis in 60% of PND patients and 44.4% of T1BD patients. We found five novel mutations in five different genes. Notably, patient 602 had a novel homozygous missense mutation c.1295C?>?A (T432?K) in the glucokinase (GCK) gene. Compared to the wild-type recombinant protein, the mutant protein had significantly lower enzymatic activity (2.5%, p?=?0.0002) and V_max (1.23?±?0.019 vs. 0.33?±?0.016, respectively; p?=?0.005). WES is a robust technique that can be used to unravel the etiologies of genetically heterogeneous forms of diabetes. Homozygous inactivating mutations of the GCK gene may have a significant role in PND pathogenesis.