Laura A. Crinnion

Robust Diagnostic Genetic Testing Using Solution Capture Enrichment and a Novel Variant‐Filtering Interface

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 methods 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.

Identification of a mutation in the ubiquitin-fold modifier 1-specific peptidase 2 gene, UFSP2, in an extended South African family with Beukes hip dysplasia.

BACKGROUND Beukes hip dysplasia (BHD) is an autosomal dominant disorder of variable penetrance that was originally identified in a large South African family of European origin. BHD is characterised by bilateral dysmorphism of the proximal femur, which results in severe degenerative osteoarthropathy. Previous studies mapped the disorder to a 3.34 Mb region on chromosome 4q35. OBJECTIVE To fine-map the BHD locus and identify the disease-causing mutation by direct sequencing. RESULTS The linked BHD allele was refined to 1.33 Mb, reducing the number of candidate genes from 25 to 16. Analysis of protein coding and invariant splice-site sequences in three distantly related individuals identified a single-candidate disease-causing variant c.868T>C within exon 8 of the ubiquitin-fold modifier 1 (Ufm1)-specific peptidase 2 gene, UFSP2. The presence of this unique mutation was confirmed in all 17 affected members of the BHD family who were genotyped. The mutation segregated with the BHD phenotype in the extended family with a two-point (single marker) LOD score of 10.4 (θ=0.0 and 80% penetrance). The mutation predicts the substitution of a highly conserved amino acid, p.Tyr290His, in the encoded protein. In vitro functional assays performed using purified recombinant wild-type and mutant UFSP2 protein demonstrated that the BHD mutation abolishes UFSP2-mediated C-terminal cleavage of its substrate, Ufm1. CONCLUSION We report a unique UFSP2 mutation that segregates with the BHD phenotype. The predicted amino acid substitution inactivates UFSP2 proteolytic function, thus implicating the ubiquitin-fold modifier 1 cascade in this form of severe hip osteoarthropathy. The facile polymerase chain reaction-based assay we describe could be used to confirm the diagnosis of BHD, or for presymptomatic testing of members of the extended BHD family.

GeneTIER: prioritization of candidate disease genes using tissue-specific gene expression profiles

Motivation: In attempts to determine the genetic causes of human disease, researchers are often faced with a large number of candidate genes. Linkage studies can point to a genomic region containing hundreds of genes, while the high-throughput sequencing approach will often identify a great number of non-synonymous genetic variants. Since systematic experimental verification of each such candidate gene is not feasible, a method is needed to decide which genes are worth investigating further. Computational gene prioritization presents itself as a solution to this problem, systematically analyzing and sorting each gene from the most to least likely to be the disease-causing gene, in a fraction of the time it would take a researcher to perform such queries manually. Results: Here, we present Gene TIssue Expression Ranker (GeneTIER), a new web-based application for candidate gene prioritization. GeneTIER replaces knowledge-based inference traditionally used in candidate disease gene prioritization applications with experimental data from tissue-specific gene expression datasets and thus largely overcomes the bias toward the better characterized genes/diseases that commonly afflict other methods. We show that our approach is capable of accurate candidate gene prioritization and illustrate its strengths and weaknesses using case study examples. Availability and Implementation: Freely available on the web at http://dna.leeds.ac.uk/GeneTIER/. Contact: umaan@leeds.ac.uk Supplementary information: Supplementary data are available at Bioinformatics online.

Diagnostic whole genome sequencing and split‐read mapping for nucleotide resolution breakpoint identification in CNTNAP2 deficiency syndrome

Whole genome sequencing (WGS) has the potential to report on all types of genetic abnormality, thus converging diagnostic testing on a single methodology. Although WGS at sufficient depth for robust detection of point mutations is still some way from being affordable for diagnostic purposes, low‐coverage WGS is already an excellent method for detecting copy number variants (“CNVseq”). We report on a family in which individuals presented with a presumed autosomal recessive syndrome of severe intellectual disability and epilepsy. Array comparative genomic hybridization (CGH) analysis had revealed a homozygous deletion apparently lying within intron 3 of CNTNAP2. Since this was too small for confirmation by FISH, CNVseq was used, refining the extent of this mutation to approximately 76.8 kb, encompassing CNTNAP2 exon 3 (an out‐of‐frame deletion). To characterize the precise breakpoints and provide a rapid molecular diagnostic test, we resequenced the CNVseq library at medium coverage and performed split read mapping. This yielded information for a multiplex polymerase chain reaction (PCR) assay, used for cascade screening and/or prenatal diagnosis in this family. This example demonstrates a rapid, low‐cost approach to converting molecular cytogenetic findings into robust PCR‐based tests.

OVA: integrating molecular and physical phenotype data from multiple biomedical domain ontologies with variant filtering for enhanced variant prioritization.

Motivation: Exome sequencing has become a de facto standard method for Mendelian disease gene discovery in recent years, yet identifying disease-causing mutations among thousands of candidate variants remains a non-trivial task. Results: Here we describe a new variant prioritization tool, OVA (ontology variant analysis), in which user-provided phenotypic information is exploited to infer deeper biological context. OVA combines a knowledge-based approach with a variant-filtering framework. It reduces the number of candidate variants by considering genotype and predicted effect on protein sequence, and scores the remainder on biological relevance to the query phenotype. We take advantage of several ontologies in order to bridge knowledge across multiple biomedical domains and facilitate computational analysis of annotations pertaining to genes, diseases, phenotypes, tissues and pathways. In this way, OVA combines information regarding molecular and physical phenotypes and integrates both human and model organism data to effectively prioritize variants. By assessing performance on both known and novel disease mutations, we show that OVA performs biologically meaningful candidate variant prioritization and can be more accurate than another recently published candidate variant prioritization tool. Availability and implementation: OVA is freely accessible at http://dna2.leeds.ac.uk:8080/OVA/index.jsp Supplementary information: Supplementary data are available at Bioinformatics online. Contact: umaan@leeds.ac.uk

Rapid Detection of Rare Deleterious Variants by Next Generation Sequencing with Optional Microarray SNP Genotype Data

Autozygosity mapping is a powerful technique for the identification of rare, autosomal recessive, disease‐causing genes. The ease with which this category of disease gene can be identified has greatly increased through the availability of genome‐wide SNP genotyping microarrays and subsequently of exome sequencing. Although these methods have simplified the generation of experimental data, its analysis, particularly when disparate data types must be integrated, remains time consuming. Moreover, the huge volume of sequence variant data generated from next generation sequencing experiments opens up the possibility of using these data instead of microarray genotype data to identify disease loci. To allow these two types of data to be used in an integrated fashion, we have developed AgileVCFMapper, a program that performs both the mapping of disease loci by SNP genotyping and the analysis of potentially deleterious variants using exome sequence variant data, in a single step. This method does not require microarray SNP genotype data, although analysis with a combination of microarray and exome genotype data enables more precise delineation of disease loci, due to superior marker density and distribution.

A tubulin alpha 8 mouse knockout model indicates a likely role in spermatogenesis but not in brain development.

Tubulin alpha 8 (Tuba8) is the most divergent member of the highly conserved alpha tubulin family, and uniquely lacks two key post-translational modification sites. It is abundantly expressed in testis and muscle, with lower levels in the brain. We previously identified homozygous hypomorphic TUBA8 mutations in human subjects with a polymicrogyria (PMG) syndrome, suggesting its involvement in development of the cerebral cortex. We have now generated and characterized a Tuba8 knockout mouse model. Homozygous mice were confirmed to lack Tuba8 protein in the testis, but did not display PMG and appeared to be neurologically normal. In response to this finding, we re-analyzed the human PMG subjects using whole exome sequencing. This resulted in identification of an additional homozygous loss-of-function mutation in SNAP29, suggesting that SNAP29 deficiency, rather than TUBA8 deficiency, may underlie most or all of the neurodevelopmental anomalies in these subjects. Nonetheless, in the mouse brain, Tuba8 specifically localised to the cerebellar Purkinje cells, suggesting that the human mutations may affect or modify motor control. In the testis, Tuba8 localisation was cell-type specific. It was restricted to spermiogenesis with a strong acrosomal localization that was gradually replaced by cytoplasmic distribution and was absent from spermatozoa. Although the knockout mice were fertile, the localisation pattern indicated that Tuba8 may have a role in spermatid development during spermatogenesis, rather than as a component of the mature microtubule-rich flagellum itself.

A distinctive oral phenotype points to FAM20A mutations not identified by Sanger sequencing.

Biallelic FAM20A mutations cause two conditions where Amelogenesis Imperfecta (AI) is the presenting feature: Amelogenesis Imperfecta and Gingival Fibromatosis Syndrome; and Enamel Renal Syndrome. A distinctive oral phenotype is shared in both conditions. On Sanger sequencing of FAM20A in cases with that phenotype, we identified two probands with single, likely pathogenic heterozygous mutations. Given the recessive inheritance pattern seen in all previous FAM20A mutation‐positive families and the potential for renal disease, further screening was carried out to look for a second pathogenic allele. Reverse transcriptase‐PCR on cDNA was used to determine transcript levels. CNVseq was used to screen for genomic insertions and deletions. In one family, FAM20A cDNA screening revealed only a single mutated FAM20A allele with the wild‐type allele not transcribed. In the second family, CNV detection by whole genome sequencing (CNVseq) revealed a heterozygous 54.7 kb duplication encompassing exons 1 to 4 of FAM20A. This study confirms the link between biallelic FAM20A mutations and the characteristic oral phenotype. It highlights for the first time examples of FAM20A mutations missed by the most commonly used mutation screening techniques. This information informed renal assessment and ongoing clinical care.

m6aViewer: software for the detection, analysis, and visualization of N6-methyladenosine peaks from m6A-seq/ME-RIP sequencing data

Recent methods for transcriptome-wide N6-methyladenosine (m6A) profiling have facilitated investigations into the RNA methylome and established m6A as a dynamic modification that has critical regulatory roles in gene expression and may play a role in human disease. However, bioinformatics resources available for the analysis of m6A sequencing data are still limited. Here, we describe m6aViewer-a cross-platform application for analysis and visualization of m6A peaks from sequencing data. m6aViewer implements a novel m6A peak-calling algorithm that identifies high-confidence methylated residues with more precision than previously described approaches. The application enables data analysis through a graphical user interface, and thus, in contrast to other currently available tools, does not require the user to be skilled in computer programming. m6aViewer and test data can be downloaded here: http://dna2.leeds.ac.uk/m6a.

Deficiency of the myogenic factor MyoD causes a perinatally lethal fetal akinesia

Background Lethal fetal akinesia deformation sequence (FADS) describes a clinically and genetically heterogeneous phenotype that includes fetal akinesia, intrauterine growth retardation, arthrogryposis and developmental anomalies. Affected babies die as a result of pulmonary hypoplasia. We aimed to identify the underlying genetic cause of this disorder in a family in which there were three affected individuals from two sibships. Methods Autosomal-recessive inheritance was suggested by a family history of consanguinity and by recurrence of the phenotype between the two sibships. We performed exome sequencing of the affected individuals and their unaffected mother, followed by autozygosity mapping and variant filtering to identify the causative gene. Results Five autozygous regions were identified, spanning 31.7 Mb of genomic sequence and including 211 genes. Using standard variant filtering criteria, we excluded all variants as being the likely pathogenic cause, apart from a single novel nonsense mutation, c.188C>A p.(Ser63*) (NM_002478.4), in MYOD1. This gene encodes an extensively studied transcription factor involved in muscle development, which has nonetheless not hitherto been associated with a hereditary human disease phenotype. Conclusions We provide the first description of a human phenotype that appears to result from MYOD1 mutation. The presentation with FADS is consistent with a large body of data demonstrating that in the mouse, MyoD is a major controller of precursor cell commitment to the myogenic differentiation programme.