Search

We report a patient case of pseudomembranous colitis associated with a monotoxin-producing Clostridioides difficile belonging to the very rarely diagnosed polymerase chain reaction (PCR) ribotype (RT) 151. To understand why this isolate was not identified using a routine commercial test, we performed a genomic analysis of RT151.Illumina short-read sequencing was performed on n = 11 RT151s from various geographical regions to study their genomic characteristics and relatedness. Subsequently, we used PacBio circular consensus sequencing to determine the complete genome sequence of isolates belonging to cryptic clades C-I and C-II, which includes the peatient isolate.We found that 1) RT151s are polyphyletic with isolates falling into clades 1 and cryptic clades C-I and C-II; 2) RT151 contains both nontoxigenic and toxigenic isolates and 3) RT151 C-II isolates contained monotoxin pathogenicity loci. The isolate from our patient case report contains a novel-pathogenicity loci insertion site, lacked tcdA and had a divergent tcdB sequence that might explain the failure of the diagnostic test.This study shows that RT151 encompasses both typical and cryptic clades and provides conclusive evidence for C. difficile infection due to clade C-II isolates that was hitherto lacking. Vigilance towards C. difficile infection as a result of cryptic clade isolates is warranted.

Duchenne and Becker muscular dystrophy are X-linked recessive inherited disorders characterized by progressive weakness due to skeletal muscle degeneration. Different mutations in the DMD gene, which encodes for dystrophin protein, are responsible for these disorders. The aim of our study was to investigate the relationship between type, size, and location of the mutation that occurs in the DMD gene and their effect on dystrophin protein expression in a cohort of 40 male dystrophinopathy patients and nine females, possible carriers. We evaluated the expression of dystrophin by immunofluorescence and immunoblotting. The mutational spectrum of the DMD gene was established by MLPA for large copy number variants, followed by HRM analysis for point mutations and sequencing of samples with an abnormal melting profile. MLPA revealed 30 deletions (75%) and three duplications (7.5%). HRM analysis accounted for seven-point mutations (17.5%). We also report four novel small mutations (c. 8507G>T, c.3021delG, c.9563_9563+1insAGCATGTTTATGATACAGCA, c.7661-60T>A) in DMD gene. Our work shows that the DNA translational open reading frame and the location of the mutation both influence the expression of dystrophin and disease severity phenotype. The proposed algorithm used in this study demonstrates its accuracy for the characterization of dystrophinopathy patients.

Pharmacogenomics is a key component of personalized medicine that promises safer and more effective drug treatment by individualizing drug choice and dose based on genetic profiles. In clinical practice, genetic biomarkers are used to categorize patients into *-alleles to predict CYP450 enzyme activity and adjust drug dosages accordingly. However, this approach leaves a large part of variability in drug response unexplained. Here, we present a proof-of-concept approach that uses continuous-scale (instead of categorical) assignments to predict enzyme activity. We used full CYP2D6 gene sequences obtained with long-read amplicon-based sequencing and cytochrome P450 (CYP) 2D6-mediated tamoxifen metabolism data from a prospective study of 561 patients with breast cancer to train a neural network. The model explained 79% of interindividual variability in CYP2D6 activity compared to 54% with the conventional *-allele approach, assigned enzyme activities to known alleles with previously reported effects, and predicted the activity of previously uncharacterized combinations of variants. The results were replicated in an independent cohort of tamoxifen-treated patients (model R-2 adjusted = 0.66 versus *-allele R-2 adjusted = 0.35) and a cohort of patients treated with the CYP2D6 substrate venlafaxine (model R2 adjusted = 0.64 versus *-allele R-2 adjusted = 0.55). Human embryonic kidney cells were used to confirm the effect of five genetic variants on metabolism of the CYP2D6 substrate bufuralol in vitro. These results demonstrate the advantage of a continuous scale and a completely phased genotype for prediction of CYP2D6 enzyme activity and could potentially enable more accurate prediction of individual drug response.

Pharmacogenomics is a key component of personalized medicine that promises safer and more effective drug treatment by individualizing drug choice and dose based on genetic profiles. In clinical practice, genetic biomarkers are used to categorize patients into *-alleles to predict CYP450 enzyme activity and adjust drug dosages accordingly. However, this approach leaves a large part of variability in drug response unexplained. Here, we present a proof-of-concept approach that uses continuous-scale (instead of categorical) assignments to predict enzyme activity. We used full

Pharmacogenomics is a key component of personalized medicine. It promises a safer and more effective drug treatment by individualizing the choice of drug and dose based on an individual’s genetic profile1,2. The majority of commonly prescribed drugs are metabolized by a small set of Cytochrome P450 (CYP) enzymes3. In clinical practice, genetic biomarkers are being used to categorize patients into predefined *-alleles to predict CYP450 enzyme activity and adjust drug dosages accordingly. Yet, this approach has important limitations as it leaves a large part of variability in drug response unexplained4,5. Here, we present a novel approach and introduce a continuous scale (instead of categorical) assignments to predict metabolic enzyme activity. The proposed strategy uses full gene sequencing data, a neural network model and CYP2D6 mediated tamoxifen metabolism from a prospective study of 561 breast cancer patients. The model explains 79% of the interindividual variability in CYP2D6 activity compared to 54% with the conventional approach. It is capable of assigning accurate enzyme activity to alleles containing previously uncharacterized combinations of variants and were replicated in an independent cohort of tamoxifen treated patients, a cohort of Venlafaxine users as well as in vitro functional assays using HEK cells. These results demonstrate the advantage of a continuous scale and a completely phased genotype for prediction of CYP450 enzyme activity and thereby enables more accurate prediction of individual drug response.

The use of blood-circulating cell-free DNA (cfDNA) as a “liquid biopsy” in oncology is being explored for its potential as a cancer biomarker. Mitochondria contain their own circular genomic entity (mitochondrial DNA, mtDNA), up to even thousands of copies per cell. The mutation rate of mtDNA is several orders of magnitude higher than that of the nuclear DNA. Tumor-specific variants have been identified in tumors along the entire mtDNA, and their number varies among and within tumors. The high mtDNA copy number per cell and the high mtDNA mutation rate make it worthwhile to explore the potential of tumor-specific cf-mtDNA variants as cancer marker in the blood of cancer patients. We used single-molecule real-time (SMRT) sequencing to profile the entire mtDNA of 19 tissue specimens (primary tumor and/or metastatic sites, and tumor-adjacent normal tissue) and 9 cfDNA samples, originating from 8 cancer patients (5 breast, 3 colon). For each patient, tumor-specific mtDNA variants were detected and traced in cfDNA by SMRT sequencing and/or digital PCR to explore their feasibility as cancer biomarker. As a reference, we measured other blood-circulating biomarkers for these patients, including driver mutations in nuclear-encoded cfDNA and cancer-antigen levels or circulating tumor cells. Four of the 24 (17%) tumor-specific mtDNA variants were detected in cfDNA, however at much lower allele frequencies compared to mutations in nuclear-encoded driver genes in the same samples. Also, extensive heterogeneity was observed among the heteroplasmic mtDNA variants present in an individual. We conclude that there is limited value in tracing tumor-specific mtDNA variants in blood-circulating cfDNA with the current methods available.

Cytochrome P450 2D6 (CYP2D6) is among the most important genes involved in drug metabolism. Specific variants are associated with changes in the enzyme's amount and activity. Multiple technologies exist to determine these variants, like the AmpliChip CYP450 test, Taqman qPCR, or Second‐Generation Sequencing, however, sequence homology between cytochrome P450 genes and pseudogene CYP2D7 impairs reliable CYP2D6 genotyping, and variant phasing cannot accurately be determined using these assays. To circumvent this, we sequenced CYP2D6 using the Pacific Biosciences RSII and obtained high‐quality, full‐length, phased CYP2D6 sequences, enabling accurate variant calling and haplotyping of the entire gene‐locus including exonic, intronic, and upstream and downstream regions. Unphased diplotypes (Roche AmpliChip CYP450 test) were confirmed for 24 of the 25 samples, including gene duplications. Cases with gene deletions required additional specific assays to resolve. In total, 61 unique variants were detected, including variants that had not previously been associated with specific haplotypes. To further aid genomic analysis using standard reference sequences, we have established an LOVD‐powered CYP2D6 gene‐variant database, and added all reference haplotypes and data reported here. We conclude that our CYP2D6 genotyping approach produces reliable CYP2D6 diplotypes and reveals information about additional variants, including phasing and copy‐number variation.

BackgroundRead binning is a key step in proper and accurate analysis of metagenomics data. Typically, this is performed by comparing metagenomics reads to known microbial sequences. However, microbial communities usually contain mixtures of hundreds to thousands of unknown bacteria. This restricts the accuracy and completeness of alignment-based approaches. The possibility of reference-free deconvolution of environmental sequencing data could benefit the field of metagenomics, contributing to the estimation of metagenome complexity, improving the metagenome assembly, and enabling the investigation of new bacterial species that are not visible using standard laboratory or alignment-based bioinformatics techniques.ResultsHere, we apply an alignment-free method that leverages on k-mer frequencies to classify reads within a single long read metagenomic dataset. In addition to a series of simulated metagenomic datasets, we generated sequencing data from a bioreactor microbiome using the PacBio RSII single-molecule real-time sequencing platform. We show that distances obtained after the comparison of k-mer profiles can reveal relationships between reads within a single metagenome, leading to a clustering per species.ConclusionsIn this study, we demonstrated the possibility to detect substructures within a single metagenome operating only with the information derived from the sequencing reads. The obtained results are highly important as they establish a principle that might potentially expand the toolkit for the detection and investigation of previously unknow microorganisms.

Background Theory predicts that dependency within host-endosymbiont interactions results in endosymbiont genome size reduction. Unexpectedly, the largest Wolbachia genome was found in the obligate, parthenogenesis-associated wFol. In this study, we investigate possible processes underlying this genome expansion by comparing a re-annotated wFol genome to other Wolbachia genomes. In addition, we also search for candidate genes related to parthenogenesis induction (PI). Results Within wFol, we found five phage WO regions representing 25.4% of the complete genome, few pseudogenized genes, and an expansion of DNA-repair genes in comparison to other Wolbachia. These signs of genome conservation were mirrored in the wFol host, the springtail F. candida, which also had an expanded DNA-repair gene family and many horizontally transferred genes. Across all Wolbachia genomes, there was a strong correlation between gene numbers of Wolbachia strains and their hosts. In order to identify genes with a potential link to PI, we assembled the genome of an additional PI strain, wLcla. Comparisons between four PI Wolbachia, including wFol and wLcla, and fourteen non-PI Wolbachia yielded a small set of potential candidate genes for further investigation. Conclusions The strong similarities in genome content of wFol and its host, as well as the correlation between host and Wolbachia gene numbers suggest that there may be some form of convergent evolution between endosymbiont and host genomes. If such convergent evolution would be strong enough to overcome the evolutionary forces causing genome reduction, it would enable expanded genomes within long-term obligate endosymbionts. Electronic supplementary material The online version of this article (10.1186/s12864-019-5492-9) contains supplementary material, which is available to authorized users.

Background Investigating how epigenetic information is transmitted through the mammalian germline is the key to understanding how this information impacts on health and disease susceptibility in offspring. EED is essential for regulating the repressive histone modification, histone 3 lysine 27 tri-methylation (H3K27me3) at many developmental genes. Results In this study, we used oocyte-specific Zp3-Cre recombinase (Zp3Cre) to delete Eed specifically in mouse growing oocytes, permitting the study of EED function in oocytes and the impact of depleting EED in oocytes on outcomes in offspring. As EED deletion occurred only in growing oocytes and females were mated to normal wild type males, this model allowed the study of oocyte programming without confounding factors such as altered in utero environment. Loss of EED from growing oocytes resulted in a significant overgrowth phenotype that persisted into adult life. Significantly, this involved increased adiposity (total fat) and bone mineral density in offspring. Similar overgrowth occurs in humans with Cohen-Gibson (OMIM 617561) and Weaver (OMIM 277590) syndromes, that result from de novo germline mutations in EED or its co-factor EZH2, respectively. Consistent with a role for EZH2 in human oocytes, we demonstrate that de novo germline mutations in EZH2 occurred in the maternal germline in some cases of Weaver syndrome. However, deletion of Ezh2 in mouse oocytes resulted in a distinct phenotype compared to that resulting from oocyte-specific deletion of Eed. Conclusions This study provides novel evidence that altering EED-dependent oocyte programming leads to compromised offspring growth and development in the next generation. Electronic supplementary material The online version of this article (10.1186/s13148-018-0526-8) contains supplementary material, which is available to authorized users.

Large variation exists in mitochondrial DNA (mtDNA) not only between but also within individuals. Also in human cancer, tumor-specific mtDNA variation exists. In this work, we describe the comparison of four methods to extract mtDNA as pure as possible from frozen tumor tissue. Also, three state-of-the-art methods for sensitive detection of mtDNA variants were evaluated. The main aim was to develop a procedure to detect low-frequent single-nucleotide mtDNA-specific variants in frozen tumor tissue. We show that of the methods evaluated, DNA extracted from cytosol fractions following exonuclease treatment results in highest mtDNA yield and purity from frozen tumor tissue (270-fold mtDNA enrichment). Next, we demonstrate the sensitivity of detection of low-frequent single-nucleotide mtDNA variants (≤1% allele frequency) in breast cancer cell lines MDA-MB-231 and MCF-7 by single-molecule real-time (SMRT) sequencing, UltraSEEK chemistry based mass spectrometry, and digital PCR. We also show de novo detection and allelic phasing of variants by SMRT sequencing. We conclude that our sensitive procedure to detect low-frequent single-nucleotide mtDNA variants from frozen tumor tissue is based on extraction of DNA from cytosol fractions followed by exonuclease treatment to obtain high mtDNA purity, and subsequent SMRT sequencing for (de novo) detection and allelic phasing of variants.

Anaerobic ammonium-oxidizing (anammox) bacteria are a group of strictly anaerobic chemolithoautotrophic microorganisms. They are capable of oxidizing ammonium to nitrogen gas using nitrite as a terminal electron acceptor, thereby facilitating the release of fixed nitrogen into the atmosphere. The anammox process is thought to exert a profound impact on the global nitrogen cycle and has been harnessed as an environment-friendly method for nitrogen removal from wastewater. In this study, we present the first closed genome sequence of an anammox bacterium, Kuenenia stuttgartiensis MBR1. It was obtained through Single-Molecule Real-Time (SMRT) sequencing of an enrichment culture constituting a mixture of at least two highly similar Kuenenia strains. The genome of the novel MBR1 strain is different from the previously reported Kuenenia KUST reference genome as it contains numerous structural variations and unique genomic regions. We find new proteins, such as a type 3b (sulf)hydrogenase and an additional copy of the hydrazine synthase gene cluster. Moreover, multiple copies of ammonium transporters and proteins regulating nitrogen uptake were identified, suggesting functional differences in metabolism. This assembly, including the genome-wide methylation profile, provides a new foundation for comparative and functional studies aiming to elucidate the biochemical and metabolic processes of these organisms.

BACKGROUND: A substantial fraction of familial colorectal cancer (CRC) and polyposis heritability remains unexplained. This study aimed to identify predisposing loci in patients with these disorders.METHODS: Homozygosity mapping was performed using 222 563 SNPs in 302 index patients with various colorectal neoplasms and 3367 controls. Linkage analysis, exome and whole-genome sequencing were performed in a family affected by microsatellite stable CRCs. Candidate variants were genotyped in 10 554 cases and 21 480 controls. Gene expression was assessed at the mRNA and protein level.RESULTS: Homozygosity mapping revealed a disease-associated region at 1q32.3 which was part of the linkage region 1q32.2-42.2 identified in the CRC family. This includes a region previously associated with risk of CRC. Sequencing identified the p.Asp1432Glu variant in the MIA3 gene (known as TANGO1 or TANGO) and 472 additional rare, shared variants within the linkage region. In both cases and controls the population frequency was 0.02% for this MIA3 variant. The MIA3 mutant allele showed predominant mRNA expression in normal, cancer and precancerous tissues. Furthermore, immunohistochemistry revealed increased expression of MIA3 in adenomatous tissues.CONCLUSIONS: Taken together, our two independent strategies associate genetic variations in chromosome 1q loci and predisposition to familial CRC and polyps, which warrants further investigation.British Journal of Cancer advance online publication: 25 July 2017; doi:10.1038/bjc.2017.240 www.bjcancer.com.

A genetic diagnosis of autosomal‐dominant polycystic kidney disease (ADPKD) is challenging due to allelic heterogeneity, high GC content, and homology of the PKD1 gene with six pseudogenes. Short‐read next‐generation sequencing approaches, such as whole‐genome sequencing and whole‐exome sequencing, often fail at reliably characterizing complex regions such as PKD1. However, long‐read single‐molecule sequencing has been shown to be an alternative strategy that could overcome PKD1 complexities and discriminate between homologous regions of PKD1 and its pseudogenes. In this study, we present the increased power of resolution for complex regions using long‐read sequencing to characterize a cohort of 19 patients with ADPKD. Our approach provided high sensitivity in identifying PKD1 pathogenic variants, diagnosing 94.7% of the patients. We show that reliable screening of ADPKD patients in a single test without interference of PKD1 homologous sequences, commonly introduced by residual amplification of PKD1 pseudogenes, by direct long‐read sequencing is now possible. This strategy can be implemented in diagnostics and is highly suitable to sequence and resolve complex genomic regions that are of clinical relevance.

Motivation: Advances in sequencing technologies and computational algorithms have enabled the study of genomic variants to dissect their functional consequence. Despite this unprecedented progress, current tools fail to reliably detect and characterize more complex allelic variants, such as short tandem repeats (STRs). We developed TSSV as an efficient and sensitive tool to specifically profile all allelic variants present in targeted loci. Based on its design, requiring only two short flanking sequences, TSSV can work without the use of a complete reference sequence to reliably profile highly polymorphic, repetitive or uncharacterized regions. Results: We show that TSSV can accurately determine allelic STR structures in mixtures with 10% representation of minor alleles or complex mixtures in which a single STR allele is shared. Furthermore, we show the universal utility of TSSV in two other independent studies: characterizing de novo mutations introduced by transcription activator-like effector nucleases (TALENs) and profiling the noise and systematic errors in an IonTorrent sequencing experiment. TSSV complements the existing tools by aiding the study of highly polymorphic and complex regions and provides a high-resolution map that can be used in a wide range of applications, from personal genomics to forensic analysis and clinical diagnostics. Availability and implementation: We have implemented TSSV as a Python package that can be installed through the command-line using pip install TSSV command. Its source code and documentation are available at https://pypi.python.org/pypi/tssv and http://www.lgtc.nl/tssv. Contact: S.Y.Anvar@lumc.nl Supplementary information: Supplementary data are available at Bioinformatics online.

Droplet digital PCR (ddPCR) is based on the isolated amplification of thousands of individual DNA molecules simultaneously, with each molecule compartmentalized in a droplet. The presence of amplified product in each droplet is indicated by a fluorescent signal, and the proportion of positive droplets allows the precise quantification of a given sequence. In this chapter we briefly outline the basis of ddPCR, and describe two different applications using the Bio-Rad QX200 system: genotyping copy number variation and quantification of Illumina sequencing libraries.

Conventional mitochondrial-DNA (MT DNA) sequencing approaches use Sanger sequencing of 20-40 partially overlapping PCR fragments per individual, which is a time- and resource-consuming process. We have developed a high-throughput, accurate, fast, and cost-effective human MT DNA sequencing approach. In this setup we first generate long-range PCR products for two partially overlapping 7.7 and 9.2 kb MT DNA-specific amplicons, add sample-specific barcodes, and sequence these on the PacBio RSII system to obtain full-length MT DNA sequences for genotyping/haplotyping purposes.

High-resolution melting analysis (HRMA) is a simple, quick, and effective method to scan and screen PCR amplicons for sequence variants. HRMA is a nondestructive closed tube assay; after PCR, DNA melting can directly be performed on the amplified samples without any purification or separation steps. For single SNP genotyping, HRMA is an attractive alternative to Sanger sequencing, restriction enzyme analysis, and hydrolysis probes.

High-resolution melting analysis (HRMA) is a simple, quick, and effective method to scan and screen PCR amplicons for sequence variants. HRMA is a nondestructive closed tube assay; after PCR, DNA melting can directly be performed on the amplified samples without any purification or separation steps. For single SNP genotyping, HRMA is an attractive alternative to Sanger sequencing, restriction enzyme analysis, and hydrolysis probes.

Droplet digital PCR (ddPCR) is based on the isolated amplification of thousands of individual DNA molecules simultaneously, with each molecule compartmentalized in a droplet. The presence of amplified product in each droplet is indicated by a fluorescent signal, and the proportion of positive droplets allows the precise quantification of a given sequence. In this chapter we briefly outline the basis of ddPCR, and describe two different applications using the Bio-Rad QX200 system: genotyping copy number variation and quantification of Illumina sequencing libraries.

Conventional mitochondrial-DNA (MT DNA) sequencing approaches use Sanger sequencing of 20-40 partially overlapping PCR fragments per individual, which is a time- and resource-consuming process. We have developed a high-throughput, accurate, fast, and cost-effective human MT DNA sequencing approach. In this setup we first generate long-range PCR products for two partially overlapping 7.7 and 9.2 kb MT DNA-specific amplicons, add sample-specific barcodes, and sequence these on the PacBio RSII system to obtain full-length MT DNA sequences for genotyping/haplotyping purposes.

AIMS/HYPOTHESIS: High levels of circulating mannan-binding lectin (MBL) are associated with the development of diabetic nephropathy and hyperglycaemia-induced vasculopathy. Here, we aimed to assess the effect of glycaemic control on circulating levels of MBL and the relationship of these levels with vascular damage. METHODS: We assessed MBL levels and corresponding MBL2 genotype, together with vascular endothelial growth factor (VEGF) levels as a marker of vascular damage, in type 1 diabetes patients with diabetic nephropathy before and after simultaneous pancreas-kidney (SPK) transplantation. We included diabetic nephropathy patients (n = 21), SPK patients (n = 37), healthy controls (n = 19), type 1 diabetes patients (n = 15) and diabetic nephropathy patients receiving only a kidney transplant (n = 15). Fourteen diabetic nephropathy patients were followed up for 12 months after SPK. RESULTS: We found elevated circulating MBL levels in diabetic nephropathy patients, and a trend towards elevated circulating MBL levels in type 1 diabetes patients, compared with healthy control individuals. MBL levels in SPK patients completely normalised and our data indicate that this predominantly occurs in patients with a polymorphism in the MBL2 gene. By contrast, MBL levels in kidney transplant only patients remained elevated, suggesting that glycaemic control but not reversal of renal failure is associated with decreased MBL levels. In line, levels of glucose and HbA1c, but not creatinine levels and estimated GFR, were correlated with MBL levels. VEGF levels were associated with levels of MBL and HbA1c in an MBL-polymorphism-dependent manner. CONCLUSIONS/INTERPRETATION: Taken together, circulating MBL levels are associated with diabetic nephropathy and are dependent on glycaemic control, possibly in an MBL2-genotype-dependent manner.

High Resolution Melting Analysis (HRMA) is a rapid and sensitive method for single nucleotide polymorphism (SNP) analysis. In the present study we present a novel HRMA assay to detect three SNPs in close proximity of each other in the first exon of the gene encoding mannosebinding lectin (MBL), a key molecule of innate immunity. These SNPs have been selected for their known biological and clinical relevance. The three SNPs in MBL2 were simultaneously determined in sixty- nine human DNA samples using HRMA and a single non- fluorescent melting probe, without any post- PCR processing of samples. Combining analyses from amplicon melting and probe melting, we have been able to discriminate ten exon 1 MBL2 genotypes with HRMA, making it a suitable tool for MBL genotyping. A second HRMA assay is presented to detect a relevant polymorphism (Y/X SNP) in the MBL2 promoter region. In conclusion, HRMA is a closed tube assay that is easy to setup and lends itself perfectly for high throughput genotyping of MBL2 variants. The present study thereby facilitates further clinical studies into the role of MBL in inflammatory and infectious disease. (C) 2010 Wiley-Liss, Inc.

Kirsten RAS (KRAS) is a small GTPase that plays a key role in Ras/mitogen-activated protein kinase signaling; somatic mutations in KRAS are frequently found in many cancers. The most common KRAS mutations result in a constitutively active protein. Accurate detection of KRAS mutations is pivotal to the molecular diagnosis of cancer and may guide proper treatment selection. Here, we describe a two-step KRAS mutation screening protocol that combines whole-genome amplification (WGA), high-resolution melting analysis (HRM) as a prescreen method for mutation carrying samples, and direct Sanger sequencing of DNA from formalin-fixed, paraffin-embedded (FFPE) tissue, from which limited amounts of DNA are available. We developed target-specific primers, thereby avoiding amplification of homologous KRAS sequences. The addition of herring sperm DNA facilitated WGA in DNA samples isolated from as few as 100 cells. KRAS mutation screening using high-resolution melting analysis on wgaDNA from formalin-fixed, paraffin-embedded tissue is highly sensitive and specific; additionally, this method is feasible for screening of clinical specimens, as illustrated by our analysis of pancreatic cancers. Furthermore, PCR on wgaDNA does not introduce genotypic changes, as opposed to unamplified genomic DNA. This method can, after validation, be applied to virtually any potentially mutated region in the genome.

Transition of the double-stranded DNA molecule to its two single strands, DNA denaturation or melting, has been used for many years to study DNA structure and composition. Recent technological advances have improved the potential of this technology, especially to detect variants in the DNA sequence. Sensitivity and specificity were increased significantly by the development of so-called saturating DNA dyes and by improvements in the instrumentation to measure the melting behavior (improved temperature precision combined with increased measurements per time unit and drop in temperature). Melt analysis using these new instruments has been designated high-resolution melting curve analysis (HRM or HRMA). Based on its ease of use, simplicity, flexibility, low cost, nondestructive nature, superb sensitivity, and specificity, HRMA is quickly becoming the tool of choice to screen patients for pathogenic variants. Here we will briefly discuss the latest developments in HRMA and review in particular other applications that have thus far received less attention, including presequence screening, single nucleotide polymorphism (SNP) typing, methylation analysis, quantification (copy number variants and mosaicism), an alternative to gel-electrophoresis and clone characterization. Together, these diverse applications make HRMA a multipurpose technology and a standard tool that should be present in any laboratory studying nucleic acids.

BackgroundMinimally differentiated acute myeloid leukemia is heterogeneous in karyotype and is defined by immature morphological and molecular characteristics. This originally French-American-British classification is still used in the new World Health Organization classification when other criteria are not met. Apart from RUNX1 mutation, no characteristic molecular aberrations are recognized.Design and MethodsWe performed whole genome single nucleotide polymorphism analysis and extensive molecular analysis in a cohort of 52 patients with minimally differentiated acute myeloid leukemia.ResultsMany recurring and potentially relevant regions of loss of heterozygosity were revealed. These point towards a variety of candidate genes that could contribute to the pathogenesis of minimally differentiated acute myeloid leukemia, including the tumor suppressor genes TP53 and NF1, and reinforced the importance of RUNX1 in this leukemia . Furthermore, for the first time in this minimally differentiated form of leukemia we detected mutations in the transactivation domain of RUNX1. Mutations in other acute myeloid leukemia associated transcriptions factors were infrequent. In contrast, FLT3, RAS, PTPN11 and JAK2 were often mutated. Irrespective of the RUNX1 mutation status, our results show that RAS signaling is the most important pathway for proliferation in minimally differentiated acute myeloid leukemia. Importantly, we found that high terminal deoxynucleotidyl transferase expression is closely associated with RUNX1 mutation, which could allow an easier diagnosis of RUNX1 mutation in this hematologic malignancy.ConclusionsOur results suggest that in patients without RUNX1 mutation, several other molecular aberrations, separately or in combination, contribute to a common minimally differentiated phenotype.

The detection of quantitative changes in genomic DNA, i.e. deletions and duplications or Copy Number Variants (CNVs), has recently gained considerable interest. First, detailed analysis of the human genome showed a surprising amount of CNVs, involving thousands of genes. Second, it was realised that the detection of CNVs as a cause of genetic disease was often neglected, but should be an essential part of a complete screening strategy. In both cases new efficient CNV screening methods, covering the entire range from specific loci to genome-wide, were behind these developments. This paper will briefly review the methods that are available to detect CNVs, discuss their strong and weak points, show some new developments and look ahead. Methods covered include microscopy, fluorescence in situ hybridization (including fiber-FISH), Southern blotting, PCR-based methods (including MLPA), array technology and massive parallel sequencing. In addition, we will show some new developments, including a 1400-plex CNV bead assay, fast-MLPA (from DNA to result in ∼6 h) and a simple Melting Curve Analysis assay to confirm potential CNVs. Using the 1400-plex CNV bead assay, targeting selected chromosomal regions only, we detected confirmed rearrangements in 9% of 320 mental retardation patients studied.

Objectives To validate the use of Real Time PCR, a widely used technique that can detect very low levels of Y chromosomal sequence, and to assess the use of a highly sensitive PCR technique, pyrophosphorolysis-activated polymerisation (PAP), for fetal sex determination using free fetal DNA (ffDNA). Methods The fetal sex was determined by Real Time PCR in 58 pregnancies using ffDNA isolated from maternal plasma. In parallel with the Real Time PCR experiments, the presence of Y chromosome sequence was also determined using PAP on 54 isolated ffDNA samples. Results Both techniques detected Y chromosome sequence at very low levels with 98% specificity and 100% sensitivity (Real Time n = 44, PAP n = 54). Furthermore, the PAP technique was shown to be more robust than the Real Time PCR as none of the samples tested failed to meet the acceptance criteria. Combining the two techniques for male fetal sex detection from maternal blood plasma increases the sensitivity and specificity to 100% in this series. Conclusions This study shows that both Real Time PCR and PAP can be used for Y chromosome detection on ffDNA. Furthermore, by using PAP in combination with Real Time PCR more reliable early prenatal sexing can be performed using ffDNA. Copyright © 2007 John Wiley & Sons, Ltd.

Familial adenomatous polyposis is most frequently caused by pathogenic variants in either the APC gene or the MUTYH gene. The detection rate of pathogenic variants depends on the severity of the phenotype and sensitivity of the screening method, including sensitivity for mosaic variants. For 171 patients with multiple colorectal polyps without previously detectable pathogenic variant, APC was reanalyzed in leukocyte DNA by one uniform technique: high-resolution melting (HRM) analysis. Serial dilution of heterozygous DNA resulted in a lowest detectable allelic fraction of 6 % for the majority of variants. HRM analysis and subsequent sequencing detected pathogenic fully heterozygous APC variants in 10 (6 %) of the patients and pathogenic mosaic variants in 2 (1 %). All these variants were previously missed by various conventional scanning methods. In parallel, HRM APC scanning was applied to DNA isolated from polyp tissue of two additional patients with apparently sporadic polyposis and without detectable pathogenic APC variant in leukocyte DNA. In both patients a pathogenic mosaic APC variant was present in multiple polyps. The detection of pathogenic APC variants in 7 % of the patients, including mosaics, illustrates the usefulness of a complete APC gene reanalysis of previously tested patients, by a supplementary scanning method. HRM is a sensitive and fast pre-screening method for reliable detection of heterozygous and mosaic variants, which can be applied to leukocyte and polyp derived DNA. Electronic supplementary material The online version of this article (doi:10.1007/s10689-015-9780-5) contains supplementary material, which is available to authorized users.

We have designed a multiplex ligation-dependent probe amplification (MLPA) assay to simultaneously screen all 79 DMD gene exons for deletions and duplications in Duchenne and Becker muscular dystrophy (DMD/BMD) patients. We validated the assay by screening 123 unrelated patients from Serbia and Montenegro already screened using multiplex PCR. MLPA screening confirmed the presence of all previously detected deletions. In addition, we detected seven new deletions, nine duplications, one point mutation, and we were able to precisely determine the extent of all rearrangements. To facilitate MLPA-based screening in laboratories lacking specific equipment, we designed the assay such that it can also be performed using agarose gel analysis and ethidium bromide staining. The MLPA assay as described provides a simple and cheap method for deletion and duplication screening in DMD/BMD patients. The assay outperforms the Beggs and Chamberlain multiplex-PCR test, and should be considered as the method of choice for an initial DNA analysis of DMD/BMD patients.

The therapeutic potential of frame-restoring exon skipping by antisense oligonucleotides (AONs) has recently been demonstrated in cultured muscle cells from a series of Duchenne muscular dystrophy (DMD) patients. To facilitate clinical application, in vivo studies in animal models are required to develop safe and efficient AON-delivery methods. However, since exon skipping is a sequence-specific therapy, it is desirable to target the human DMD gene directly. We therefore set up human sequence-specific exon skipping in transgenic mice carrying the full-size human gene (hDMD). We initially compared the efficiency and toxicity of intramuscular AON injections using different delivery reagents in wild-type mice. At a dose of 3.6 nmol AON and using polyethylenimine, the skipping levels accumulated up to 3% in the second week postinjection and lasted for 4 weeks. We observed a correlation of this long-term effect with the intramuscular persistence of the AON. In regenerating myofibers higher efficiencies (up to 9%) could be obtained. Finally, using the optimized protocols in hDMD mice, we were able to induce the specific skipping of human DMD exons without affecting the endogenous mouse gene. These data highlight the high sequence specificity of this therapy and present the hDMD mouse as a unique model to optimize human-specific exon skipping in vivo.

Duchenne and Becker muscular dystrophy (DMD and BMD) are caused by mutations in the dystrophin gene. Large rearrangements in the gene are found in about two,thirds of DMD patients, with similar to60% carrying deletions and 5-10% carrying duplications. Most of the remaining 30-35% of patients are expected to have small nucleotide substitutions, insertions, or deletions. To detect these subtle changes within the coding and splice site determining sequences of the dystrophin gene, we established a semiautomated denaturing gradient gel electrophoresis (DGGE) mutation scanning system. The DGGE scan covers the dystrophin gene with 95 amplicons, PCRed either individually or in a multiplex setup. PCR and pooling were performed semiautomatically, using a pipetting robot and 384-well plates, enabling concurrent amplification of DNA of four patients in one run. Amplification of individual fragments was performed using one PCR program. The products were pooled just before gel loading; DGGE requires only a single gel condition. Validation was performed using DNA samples harboring 39 known DMD variants, all of which could be readily detected. DGGE mutation scanning was applied to analyze 135 DMD/BMD patients and potential DMD carriers without large deletions or duplications. In DNA from 25 out of 44 DMD patients (57%) and from 5 out of 39 BMD patients (13%), we identified clear pathogenic changes. All mutations were different, with the exception of one DMD mutation, which occurred twice. In DNA from 10 out of 44 potential DMD carriers, including four obligate carriers, we detected causative changes, including one pathogenic change in every obligate carrier. In addition to these pathogenic changes, we detected 15 unique unclassified variants, i.e., changes for which a pathogenic nature is uncertain. (C) 2003 Wiley-Liss, Inc.

Hemoglobinopathies, such as sickle cell disease (SCD) and beta-thalassemia major (TM), are severe diseases and the most common autosomal recessive condition worldwide and in particular in Oman. Early screening and diagnosis of carriers are the key for primary prevention. Once a country-wide population screening program is mandated by law, a sequencing technology that can rapidly confirm or identify disease-causing mutations for a large number of patients in a short period of time will be necessary. While Sanger sequencing is the standard protocol for molecular diagnosis, next generation sequencing starts to become available to reference laboratories. Using the Ion Torrent PGM sequencer, we have analyzed a cohort of 297 unrelated Omani cases and reliably identified mutations in the beta-globin (HBB) gene. Our model study has shown that Ion Torrent PGM can rapidly sequence such a small gene in a large number of samples using a barcoded uni-directional or bi-directional sequence methodology, reducing cost, workload and providing accurate diagnosis. Based on our results we believe that the Ion Torrent PGM sequencing platform, able to analyze hundreds of patients simultaneously for a single disease gene can be a valid molecular screening alternative to ABI sequencing in the diagnosis of hemoglobinopathies and other genetic disorders in the near future.

A funding source was omitted from the published funding statement. The authors would like to add: The research leading to these results has received funding from the European Union's Seventh Framework Programme (FP7/2007-2011) under grant agreement n° 259679 (IDEAL).

Antisense oligonucleotide (AON)-mediated exon skipping aimed at restoring the reading frame is a promising therapeutic approach for Duchenne muscular dystrophy that is currently tested in clinical trials. Numerous AONs have been tested in (patient-derived) cultured muscle cells and the mdx mouse model. The main outcome to measure AON efficiency is usually the exon-skipping percentage, though different groups use different methods to assess these percentages. Here, we compare a series of techniques to quantify exon skipping levels in AON-treated mdx mouse muscle. We compared densitometry of RT-PCR products on ethidium bromide-stained agarose gels, primary and nested RT-PCR followed by bioanalyzer analysis and melting curve analysis. The digital array system (Fluidigm) allows absolute quantification of skipped vs non-skipped transcripts and was used as a reference. Digital array results show that 1 ng of mdx gastrocnemius muscle-derived mRNA contains B1100 dystrophin transcripts and that 665 transcripts are sufficient to determine exon-skipping levels. Quantification using bioanalyzer or densitometric analysis of primary PCR products resulted in values close to those obtained with digital array. The use of the same technique allows comparison between different groups working on exon skipping in the mdx mouse model. Laboratory Investigation (2010) 90, 1396-1402; doi:10.1038/labinvest.2010.98; published online 10 May 2010

Keratosis Follicularis Spinulosa Decalvans (KFSD) is a rare genetic disorder characterized by development of hyperkeratotic follicular papules on the scalp followed by progressive alopecia of the scalp, eyelashes, and eyebrows. Associated eye findings include photophobia in childhood and corneal dystrophy. Due to the genetic and clinical heterogeneity of similar disorders, a definitive diagnosis of KFSD is often challenging. Toward identification of the causative gene we reanalyzed a large Dutch KFSD family. SNP arrays (1 M) redefined the locus to a 2.9-Mb region at Xp22.12-Xp22.11. Screening of all 14 genes in the candidate region identified MBTPS2 as the candidate gene carrying a c.1523A>G (p.Asn508Ser) missense mutation. The variant was also identified in two unrelated X-linked KFSD families and cosegregated with KFSD in all families. In symptomatic female carriers, skewed X-inactivation of the normal allele matched with increased severity of symptoms. MBTPS2 is required for cleavage of sterol regulatory element-binding proteins (SREBPs). In vitro functional expression studies of the c. 1523A>G mutation showed that sterol responsiveness was reduced by half. Other missense mutations in MBTPS2 have recently been identified in patients with IFAP syndrome. We postulate that both phenotypes are in the spectrum of one genetic disorder with a partially overlapping phenotype. Hum Mutat 31:1125-1133, 2010. (C) 2010 Wiley-Liss, Inc.

Background Methodologies like phage display selection, in vitro mutagenesis and the determination of allelic expression differences include steps where large numbers of clones need to be compared and characterised. In the current study we show that high-resolution melt curve analysis (HRMA) is a simple, cost-saving tool to quickly study clonal variation without prior nucleotide sequence knowledge. Results HRMA results nicely matched those obtained with ELISA and compared favourably to DNA fingerprinting of restriction digested clone insert-PCR. DNA sequence analysis confirmed that HRMA-clustered clones contained identical inserts. Conclusion Using HRMA, analysis of up to 384 samples can be done simultaneously and will take approximately 30 minutes. Clustering of clones can be largely automated using the system's software within 2 hours. Applied to the analysis of clones obtained after phage display antibody selection, HRMA facilitated a quick overview of the overall success as well as the identification of identical clones. Our approach can be used to characterize any clone set prior to sequencing, thereby reducing sequencing costs significantly.

The protein truncation test detects mutations at the protein level that lead to premature translation termination. The method has evolved considerably since is original publication in this manual. This thoroughly revised unit describes what is now the preferred method for performing the protein truncation test. Transcription and translation are performed in separate reactions; during translation, biotin-labeled or N-terminally tagged proteins are synthesized. The translation products are detected on immunoblots via chemiluminescence. An Alternate Protocol using coupled in vitro transcription/translation and radiolabeled proteins is also presented.

We report the generation of mice with an intact and functional copy of the 2.3-megabase human dystrophin gene (hDMD), the largest functional stretch of human DNA thus far integrated into a mouse chromosome. Yeast spheroplasts containing an artificial chromosome with the full-length hDMD gene were fused with mouse embryonic stem cells and were subsequently injected into mouse blastocysts to produce transgenic hDMD mice. Human-specific PCR, Southern blotting, and fluorescent in situ hybridization techniques demonstrated the intactness and stable chromosomal integration of the hDMD gene on mouse chromosome 5. Expression of the transgene was confirmed by RT-PCR and Western blotting. The tissue-specific expression pattern of the different DMD transcripts was maintained. However, the human Dp427p and Dp427m transcripts were expressed at 2-fold higher levels and human Dp427c and Dp260 transcripts were expressed at 2- and 4-fold lower levels than their endogenous counterparts. Ultimate functional proof of the hDMD transgene was obtained by crossing of hDMD mice with dystrophin-deficient mdx mice and dystrophin and utrophin-deficient mdx x Utrn-/- mice. The hDMD transgene rescued the lethal dystrophic phenotype of the mdx x Utrn-/- mice. All signs of muscular dystrophy disappeared in the rescued mice, as demonstrated by histological staining of muscle sections and gene expression profiling experiments. Currently, hDMD mice are extensively used for preclinical testing of sequence-specific therapeutics for the treatment of Duchenne muscular dystrophy. In addition, the hDMD mouse can be used to study the influence of the genomic context on deletion and recombination frequencies, genome stability, and gene expression regulation.

BackgroundThe comparability of gene expression data generated with different microarray platforms is still a matter of concern. Here we address the performance and the overlap in the detection of differentially expressed genes for five different microarray platforms in a challenging biological context where differences in gene expression are few and subtle.ResultsGene expression profiles in the hippocampus of five wild-type and five transgenic δC-doublecortin-like kinase mice were evaluated with five microarray platforms: Applied Biosystems, Affymetrix, Agilent, Illumina, LGTC home-spotted arrays. Using a fixed false discovery rate of 10% we detected surprising differences between the number of differentially expressed genes per platform. Four genes were selected by ABI, 130 by Affymetrix, 3,051 by Agilent, 54 by Illumina, and 13 by LGTC. Two genes were found significantly differentially expressed by all platforms and the four genes identified by the ABI platform were found by at least three other platforms. Quantitative RT-PCR analysis confirmed 20 out of 28 of the genes detected by two or more platforms and 8 out of 15 of the genes detected by Agilent only. We observed improved correlations between platforms when ranking the genes based on the significance level than with a fixed statistical cut-off. We demonstrate significant overlap in the affected gene sets identified by the different platforms, although biological processes were represented by only partially overlapping sets of genes. Aberrances in GABA-ergic signalling in the transgenic mice were consistently found by all platforms.ConclusionThe different microarray platforms give partially complementary views on biological processes affected. Our data indicate that when analyzing samples with only subtle differences in gene expression the use of two different platforms might be more attractive than increasing the number of replicates. Commercial two-color platforms seem to have higher power for finding differentially expressed genes between groups with small differences in expression.

Mannose-binding lectin (MBL) is a key molecule of innate immunity. Binding of MBL to carbohydrates present on pathogens activates the lectin pathway of complement activation, resulting into opsonization and anti-microbial protection. Three frequently occurring single nucleotide polymorphisms (SNPs) are described in the coding region of the MBL2 gene that are associated with abnormal polymerization of the MBL molecule, decreased serum concentrations of high molecular weight MBL, and strongly impaired function. Clinical studies have shown that these MBL SNPs are associated with increased susceptibility to infections, especially in immune-compromised persons, as well as with accelerated progression of chronic diseases. The present study describes a novel method to detect the three major MBL SNPs by pyrosequencing. The close proximity of these SNPs allows their detection in one single pyrosequencing reaction, resulting in clearly distinguishable patterns for each allele combination described until now. This method can be used for the easy and reliable detection of MBL SNPs to identify the basis of functional MBL deficiency in clinical diagnostics and research.

Germline mutations in the adenomatous polyposis coli gene cause familial adenomatous polyposis, a colon cancer predisposition syndrome. More than 95% of the identified mutations result in the generation of stop codons or reading frame shifts and encode a truncated gene product, a mutation profile also found in other tumor predisposition genes such as the breast cancer or the hereditary non-polyposis coli. Therefore the protein truncation test is ideally suited for screening of mutations in these genes, starting from simple blood samples. Gene segments of interest are amplified from genomic DNA or mRNA, thereby incorporating a T7 promoter at the 5′-end. After in vitro transcription and translation of the PCR products, the resulting protein is analysed by gel electrophoresis. Truncated translation products indicate the presence of a stop mutation. We have developed a non-radioactive protein truncation test that uses a biotinylated Lys-t-RNA to label the translation products and allows a chemiluminescent detection instead of the standard radioactive method. This generic protein truncation test kit was then used to develop a parameter-specific protein truncation test for adenomatous polyposis coli. The adenomatous polyposis coli gene was divided in 5 overlapping segments, and primers were optimized to produce distinct bands with very low background in the protein truncation test. The assay was tested on 20 familial adenomatous polyposis patient samples, where 18 mutations were found, demonstrating the efficiency of this method.

Facioscapulohumeral muscular dystrophy (FSHD) has recently been shown to be associated with deletions that are detectable using probe p13E-11 (D4F104S1). Although these deletions reside within large, highly polymorphic restriction fragments (20-300 kb), the "mutant" fragment is usually shorter than 28 kb and can routinely be detected using conventional agarose gel electrophoresis. Yet, the complete visualization of the alleles requires pulsed-field gel electrophoresis (PFGE). Family studies showed that p13E-11 detects two nonallelic loci in this size range, only one of which originates from chromosome 4q35. We have assigned the other p13E-11 locus to chromosome 10qter by linkage analysis in CEPH pedigrees. Knowing the location of both loci improves the diagnostic reliability, as the exact origin of "small" EcoRI fragments can be determined by haplotyping. Since FSHD shows genetic heterogeneity, this 10qter locus became an interesting candidate to be the second FSHD locus. However, analysis of a large chromosome 4-unlinked FSHD family did not provide evidence for linkage on chromosome 10qter.

Histidine-rich glycoprotein (HRG) belongs to the cystatin superfamily and appears to be a potential risk factor for thrombosis. An increased prevalence of elevated HRG plasma levels in patients with venous thrombosis and families with thrombophilia has been reported. It is interesting to note that the genes of four different members of the cystatin superfamily are located on the distal section of the long arm of chromosome 3: Stefin A (STF1) on 3q21, Kininogen (KNG) on 3q26-qter, [alpha]-2-HS-glycoprotein (AHSG) on 3q27-q28, and HRG on 3q21-qter. To further investigate the evolutionary relationship between HRG and members of the cystatin superfamily, the authors isolated a cosmid that was used to refine the chromosomal localization of HRG by in situ hybridization. In addition, they used a dinucleotide repeat polymorphism to localize HRG on the linkage map of chromosome 3q. 10 refs., 2 figs.

The hippocampal expression profiles of wild-type mice and mice transgenic for deltaC-doublecortin-like kinase were compared with Solexa/Illumina deep sequencing technology and five different microarray platforms. With Illumina's digital gene expression assay, we obtained approximately 2.4 million sequence tags per sample, their abundance spanning four orders of magnitude. Results were highly reproducible, even across laboratories. With a dedicated Bayesian model, we found differential expression of 3179 transcripts with an estimated false-discovery rate of 8.5%. This is a much higher figure than found for microarrays. The overlap in differentially expressed transcripts found with deep sequencing and microarrays was most significant for Affymetrix. The changes in expression observed by deep sequencing were larger than observed by microarrays or quantitative PCR. Relevant processes such as calmodulin-dependent protein kinase activity and vesicle transport along microtubules were found affected by deep sequencing but not by microarrays. While undetectable by microarrays, antisense transcription was found for 51% of all genes and alternative polyadenylation for 47%. We conclude that deep sequencing provides a major advance in robustness, comparability and richness of expression profiling data and is expected to boost collaborative, comparative and integrative genomics studies.