Christopher M. Watson

Large‐scale analysis of association between GDF5 and FRZB variants and osteoarthritis of the hip, knee, and hand

OBJECTIVE GDF5 and FRZB have been proposed as genetic loci conferring susceptibility to osteoarthritis (OA); however, the results of several studies investigating the association of OA with the rs143383 polymorphism of the GDF5 gene or the rs7775 and rs288326 polymorphisms of the FRZB gene have been conflicting or inconclusive. To examine these associations, we performed a large-scale meta-analysis of individual-level data. METHODS Fourteen teams contributed data on polymorphisms and knee, hip, and hand OA. For rs143383, the total number of cases and controls, respectively, was 5,789 and 7,850 for hip OA, 5,085 and 8,135 for knee OA, and 4,040 and 4,792 for hand OA. For rs7775, the respective sample sizes were 4,352 and 10,843 for hip OA, 3,545 and 6,085 for knee OA, and 4,010 and 5,151 for hand OA, and for rs288326, they were 4,346 and 8,034 for hip OA, 3,595 and 6,106 for knee OA, and 3,982 and 5,152 for hand OA. For each individual study, sex-specific odds ratios (ORs) were calculated for each OA phenotype that had been investigated. The ORs for each phenotype were synthesized using both fixed-effects and random-effects models for allele-based effects, and also for haplotype effects for FRZB. RESULTS A significant random-effects summary OR for knee OA was demonstrated for rs143383 (1.15 [95% confidence interval 1.09-1.22]) (P=9.4x10(-7)), with no significant between-study heterogeneity. Estimates of effect sizes for hip and hand OA were similar, but a large between-study heterogeneity was observed, and statistical significance was borderline (for OA of the hip [P=0.016]) or absent (for OA of the hand [P=0.19]). Analyses for FRZB polymorphisms and haplotypes did not reveal any statistically significant signals, except for a borderline association of rs288326 with hip OA (P=0.019). CONCLUSION Evidence of an association between the GDF5 rs143383 polymorphism and OA is substantially strong, but the genetic effects are consistent across different populations only for knee OA. Findings of this collaborative analysis do not support the notion that FRZB rs7775 or rs288326 has any sizable genetic effect on OA phenotypes.

CCDC151 mutations cause primary ciliary dyskinesia by disruption of the outer dynein arm docking complex formation.

A diverse family of cytoskeletal dynein motors powers various cellular transport systems, including axonemal dyneins generating the force for ciliary and flagellar beating essential to movement of extracellular fluids and of cells through fluid. Multisubunit outer dynein arm (ODA) motor complexes, produced and preassembled in the cytosol, are transported to the ciliary or flagellar compartment and anchored into the axonemal microtubular scaffold via the ODA docking complex (ODA-DC) system. In humans, defects in ODA assembly are the major cause of primary ciliary dyskinesia (PCD), an inherited disorder of ciliary and flagellar dysmotility characterized by chronic upper and lower respiratory infections and defects in laterality. Here, by combined high-throughput mapping and sequencing, we identified CCDC151 loss-of-function mutations in five affected individuals from three independent families whose cilia showed a complete loss of ODAs and severely impaired ciliary beating. Consistent with the laterality defects observed in these individuals, we found Ccdc151 expressed in vertebrate left-right organizers. Homozygous zebrafish ccdc151ts272a and mouse Ccdc151Snbl mutants display a spectrum of situs defects associated with complex heart defects. We demonstrate that CCDC151 encodes an axonemal coiled coil protein, mutations in which abolish assembly of CCDC151 into respiratory cilia and cause a failure in axonemal assembly of the ODA component DNAH5 and the ODA-DC-associated components CCDC114 and ARMC4. CCDC151-deficient zebrafish, planaria, and mice also display ciliary dysmotility accompanied by ODA loss. Furthermore, CCDC151 coimmunoprecipitates CCDC114 and thus appears to be a highly evolutionarily conserved ODA-DC-related protein involved in mediating assembly of both ODAs and their axonemal docking machinery onto ciliary microtubules.

Concordant association of insulin degrading enzyme gene (IDE) variants with IDE mRNA, Abeta, and Alzheimer's disease.

Background The insulin-degrading enzyme gene (IDE) is a strong functional and positional candidate for late onset Alzheimers disease (LOAD). Methodology/Principal Findings We examined conserved regions of IDE and its 10 kb flanks in 269 AD cases and 252 controls thereby identifying 17 putative functional polymorphisms. These variants formed eleven haplotypes that were tagged with ten variants. Four of these showed significant association with IDE transcript levels in samples from 194 LOAD cerebella. The strongest, rs6583817, which has not previously been reported, showed unequivocal association (p = 1.5×10−8, fold-increase = 2.12,); the eleven haplotypes were also significantly associated with transcript levels (global p = 0.003). Using an in vitro dual luciferase reporter assay, we found that rs6583817 increases reporter gene expression in Be(2)-C (p = 0.006) and HepG2 (p = 0.02) cell lines. Furthermore, using data from a recent genome-wide association study of two Croatian isolated populations (n = 1,879), we identified a proxy for rs6583817 that associated significantly with decreased plasma Aβ40 levels (ß = −0.124, p = 0.011) and total measured plasma Aβ levels (b = −0.130, p = 0.009). Finally, rs6583817 was associated with decreased risk of LOAD in 3,891 AD cases and 3,605 controls. (OR = 0.87, p = 0.03), and the eleven IDE haplotypes (global p = 0.02) also showed significant association. Conclusions Thus, a previously unreported variant unequivocally associated with increased IDE expression was also associated with reduced plasma Aß40 and decreased LOAD susceptibility. Genetic association between LOAD and IDE has been difficult to replicate. Our findings suggest that targeted testing of expression SNPs (eSNPs) strongly associated with altered transcript levels in autopsy brain samples may be a powerful way to identify genetic associations with LOAD that would otherwise be difficult to detect.

Diagnosis of copy number variation by Illumina next generation sequencing is comparable in performance to oligonucleotide array comparative genomic hybridisation.

Array comparative genomic hybridisation (aCGH) profiling is currently the gold standard for genetic diagnosis of copy number. Next generation sequencing technologies provide an alternative and adaptable method of detecting copy number by comparing the number of sequence reads in non-overlapping windows between patient and control samples. Detection of copy number using the BlueGnome 8×60k oligonucleotide aCGH platform was compared with low resolution next generation sequencing using the Illumina GAIIx on 39 patients with developmental delay and/or learning difficulties who were referred to the Leeds Clinical Cytogenetics Laboratory. Sensitivity and workflow of the two platforms were compared. Customised copy number algorithms assessed sequence counts and detected changes in copy number. Imbalances detected on both platforms were compared. Of the thirty-nine patients analysed, all eleven imbalances detected by array CGH and confirmed by FISH or Q-PCR were also detected by CNV-seq. In addition, CNV-seq reported one purported pathogenic copy number variant that was not detected by array CGH. Non-pathogenic, unconfirmed copy number calls were detected by both platforms; however few were concordant between the two. CNV-seq offers an alternative to array CGH for copy number analysis with resolution and future costs comparable to conventional array CGH platforms and with less stringent sample requirements.

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.

Structure of Ubiquitin-fold Modifier 1-specific Protease UfSP2

Ubiquitin-fold modifier 1 (Ufm1)-specific protease 2 (UfSP2) is a cysteine protease that is responsible for the release of Ufm1 from Ufm1-conjugated cellular proteins, as well as for the generation of mature Ufm1 from its precursor. The 2.6 Å resolution crystal structure of mouse UfSP2 reveals that it is composed of two domains. The C-terminal catalytic domain is similar to UfSP1 with Cys294, Asp418, His420, Tyr282, and a regulatory loop participating in catalysis. The novel N-terminal domain shows a unique structure and plays a role in the recognition of its cellular substrate C20orf116 and thus in the recruitment of UfSP2 to the endoplasmic reticulum, where C20orf116 predominantly localizes. Mutagenesis studies were carried out to provide the structural basis for understanding the loss of catalytic activity observed in a recently identified UfSP2 mutation that is associated with an autosomal dominant form of hip dysplasia.

Biallelic Mutations in PDE10A Lead to Loss of Striatal PDE10A and a Hyperkinetic Movement Disorder with Onset in Infancy

Deficits in the basal ganglia pathways modulating cortical motor activity underlie both Parkinson disease (PD) and Huntington disease (HD). Phosphodiesterase 10A (PDE10A) is enriched in the striatum, and animal data suggest that it is a key regulator of this circuitry. Here, we report on germline PDE10A mutations in eight individuals from two families affected by a hyperkinetic movement disorder due to homozygous mutations c.320A>G (p.Tyr107Cys) and c.346G>C (p.Ala116Pro). Both mutations lead to a reduction in PDE10A levels in recombinant cellular systems, and critically, positron-emission-tomography (PET) studies with a specific PDE10A ligand confirmed that the p.Tyr107Cys variant also reduced striatal PDE10A levels in one of the affected individuals. A knock-in mouse model carrying the homologous p.Tyr97Cys variant had decreased striatal PDE10A and also displayed motor abnormalities. Striatal preparations from this animal had an impaired capacity to degrade cyclic adenosine monophosphate (cAMP) and a blunted pharmacological response to PDE10A inhibitors. These observations highlight the critical role of PDE10A in motor control across species.

Association between missense mutations in the BBS2 gene and nonsyndromic retinitis pigmentosa.

IMPORTANCE A large number of genes can cause inherited retinal degenerations when mutated. It is important to identify the cause of disease for a better disease prognosis and a possible gene-specific therapeutic intervention. OBJECTIVE To identify the cause of disease in families with nonsyndromic retinitis pigmentosa. DESIGN, SETTING, AND PARTICIPANTS Patients and family members were recruited for the study and underwent clinical evaluation and genetic analyses. MAIN OUTCOMES AND MEASURES Identification of sequence variants in genes using next-generation sequencing. RESULTS We performed exome sequencing for 4 families, which was followed by Sanger sequencing of the identified mutations in 120 ethnicity-matched patients. In total, we identified 4 BBS2 missense mutations that cause nonsyndromic retinitis pigmentosa. Three siblings of Moroccan Jewish ancestry were compound heterozygotes for p.A33D and p.P134R, and 6 patients belonging to 4 families of Ashkenazi Jewish ancestry were homozygous for either p.D104A or p.R632P, or compound heterozygous for these 2 mutations. The mutations cosegregated with retinitis pigmentosa in the studied families, and the affected amino acid residues are evolutionarily conserved. CONCLUSIONS AND RELEVANCE Our study shows that BBS2 mutations can cause nonsyndromic retinitis pigmentosa and highlights yet another candidate for this genetically heterogeneous condition.

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.