Search

Targeted hybridization enrichment prior to next-generation sequencing is a widespread method for characterizing sequence variation in a research setting, and is being adopted by diagnostic laboratories. However, the number of variants identified can overwhelm clinical laboratories with strict time constraints, the final interpretation of likely pathogenicity being a particular bottleneck. To address this, we have developed an approach in which, after automatic variant calling on a standard unix pipeline, subsequent variant filtering is performed interactively, using AgileExomeFilter and AgilePindelFilter (http://dna.leeds.ac.uk/agile), tools designed for clinical scientists with standard desktop computers. To demonstrate the method's diagnostic efficacy, we tested 128 patients using (1) a targeted capture of 36 cancer-predisposing genes or (2) whole-exome capture for diagnosis of the genetically heterogeneous disorder primary ciliary dyskinesia (PCD). In the cancer cohort, complete concordance with previous diagnostic data was achieved across 793 variant genotypes. A high yield (42%) was also achieved for exome-based PCD diagnosis, underscoring the scalability of our method. Simple adjustments to the variant filtering parameters further allowed the identification of a homozygous truncating mutation in a presumptive new PCD gene, DNAH8. These tools should allow diagnostic laboratories to expand their testing portfolios flexibly, using a standard set of reagents and techniques., (© 2013 The Authors. *Human Mutation published by Wiley Periodicals, Inc.)

Massively parallel ("next generation") DNA sequencing (NGS) has quickly become the method of choice for seeking pathogenic mutations in rare uncharacterized monogenic diseases. Typically, before DNA sequencing, protein-coding regions are enriched from patient genomic DNA, representing either the entire genome ("exome sequencing") or selected mapped candidate loci. Sequence variants, identified as differences between the patient's and the human genome reference sequences, are then filtered according to various quality parameters. Changes are screened against datasets of known polymorphisms, such as dbSNP and the 1000 Genomes Project, in the effort to narrow the list of candidate causative variants. An increasing number of commercial services now offer to both generate and align NGS data to a reference genome. This potentially allows small groups with limited computing infrastructure and informatics skills to utilize this technology. However, the capability to effectively filter and assess sequence variants is still an important bottleneck in the identification of deleterious sequence variants in both research and diagnostic settings. We have developed an approach to this problem comprising a user-friendly suite of programs that can interactively analyze, filter and screen data from enrichment-capture NGS data. These programs ("Agile Suite") are particularly suitable for small-scale gene discovery or for diagnostic analysis., (© 2013 WILEY PERIODICALS, INC.)

Heterozygous mutations in DNA mismatch repair (MMR) genes result in predisposition to colorectal cancer (hereditary nonpolyposis colorectal cancer or Lynch syndrome). Patients with biallelic mutations in these genes, however, present earlier, with constitutional mismatch repair deficiency cancer syndrome (CMMRD), which is characterized by a spectrum of rare childhood malignancies and café-au-lait skin patches. The hallmark of MMR deficiency, microsatellite instability (MSI), is readily detectable in tumor DNA in Lynch syndrome, but is also present in constitutional DNA of CMMRD patients. However, detection of constitutional or germline MSI (gMSI) has hitherto relied on technically difficult assays that are not routinely applicable for clinical diagnosis. Consequently, we have developed a simple high-throughput screening methodology to detect gMSI in CMMRD patients based on the presence of stutter peaks flanking a dinucleotide repeat allele when amplified from patient blood DNA samples. Using the three different microsatellite markers, the gMSI ratio was determined in a cohort of normal individuals and 10 CMMRD patients, with biallelic germline mutations in PMS2 (seven patients), MSH2 (one patient), or MSH6 (two patients). Subjects with either PMS2 or MSH2 mutations were easily identified; however, this measure was not altered in patients with CMMRD due to MSH6 mutation., (© 2013 Wiley Periodicals, Inc.)

Pachydermoperiostosis, or primary hypertrophic osteoarthropathy (PHO), is an inherited multisystem disorder, whose features closely mimic the reactive osteoarthropathy that commonly accompanies neoplastic and inflammatory pathologies. We previously described deficiency of the prostaglandin-degrading enzyme 15-hydroxyprostaglandin dehydrogenase (HPGD) as a cause of this condition, implicating elevated circulating prostaglandin E(2) (PGE(2)) as causative of PHO, and perhaps also as the principal mediator of secondary HO. However, PHO is genetically heterogeneous. Here, we use whole-exome sequencing to identify recessive mutations of the prostaglandin transporter SLCO2A1, in individuals lacking HPGD mutations. We performed exome sequencing of four probands with severe PHO, followed by conventional mutation analysis of SLCO2A1 in nine others. Biallelic SLCO2A1 mutations were identified in 12 of the 13 families. Affected individuals had elevated urinary PGE(2), but unlike HPGD-deficient patients, also excreted considerable quantities of the PGE(2) metabolite, PGE-M. Clinical differences between the two groups were also identified, notably that SLCO2A1-deficient individuals have a high frequency of severe anemia due to myelofibrosis. These findings reinforce the key role of systemic or local prostaglandin excess as the stimulus to HO. They also suggest that the induction or maintenance of hematopoietic stem cells by prostaglandin may depend upon transporter activity., (© 2012 Wiley Periodicals, Inc.)

Autozygosity mapping has been a powerful method for the identification of autosomal recessive disease genes. However, the approach is limited by the availability of suitable consanguineous pedigrees. While rare autosomal recessive diseases are overrepresented in consanguineous families, a significant proportion of affected patients nonetheless originate in families where the parents are apparently unrelated. However, due to their relative rarity and the heterogeneity of disease alleles, it has proved difficult to use these patients to identify disease loci. Therefore, we developed "Phaser," a computer application that is able to infer the phase of SNP alleles and so haplotype entire chromosomes in small nuclear families (http://dna.leeds.ac.uk/Phaser). Once the index case's chromosomes have been haplotyped, it is then possible to deduce those of the parents and subsequently identify the parental origin of all the siblings' DNA. By combining information from a small number of nuclear families, it may then be possible to identify linkage to the recessive disease locus, in both in-bred and out-bred families. We have illustrated the program's utility by using it to correctly identify both the cystic fibrosis locus (using two unrelated compound heterozygous CEPH families) and a new gene mutated in early-onset myopathy with respiratory distress and dysphagia locus in a single consanguineous pedigree., (© 2011 Wiley Periodicals, Inc.)