Defects in cilia formation and function result in a range of human skeletal and visceral abnormalities. Mutations in several genes have been identified to cause a proportion of these disorders, some of which display genetic (locus) heterogeneity. Mouse models are valuable for dissecting the function of these genes, as well as for more detailed analysis of the underlying developmental defects. The short-rib polydactyly (SRP) group of disorders are among the most severe human phenotypes caused by cilia dysfunction. We mapped the disease locus from two siblings affected by a severe form of SRP to 2p24, where we identified an in-frame homozygous deletion of exon 5 in WDR35. We subsequently found compound heterozygous missense and nonsense mutations in WDR35 in an independent second case with a similar, severe SRP phenotype. In a mouse mutation screen for developmental phenotypes, we identified a mutation in Wdr35 as the cause of midgestation lethality, with abnormalities characteristic of defects in the Hedgehog signaling pathway. We show that endogenous WDR35 localizes to cilia and centrosomes throughout the developing embryo and that human and mouse fibroblasts lacking the protein fail to produce cilia. Through structural modeling, we show that WDR35 has strong homology to the COPI coatamers involved in vesicular trafficking and that human SRP mutations affect key structural elements in WDR35. Our report expands, and sheds new light on, the pathogenesis of the SRP spectrum of ciliopathies.

Gillespie syndrome (GS) is characterized by bilateral iris hypoplasia, congenital hypotonia, non-progressive ataxia, and progressive cerebellar atrophy. Trio-based exome sequencing identified de novo mutations in ITPR1 in three unrelated individuals with GS recruited to the Deciphering Developmental Disorders study. Whole-exome or targeted sequence analysis identified plausible disease-causing ITPR1 mutations in 10/10 additional GS-affected individuals. These ultra-rare protein-altering variants affected only three residues in ITPR1: Glu2094 missense (one de novo, one co-segregating), Gly2539 missense (five de novo, one inheritance uncertain), and Lys2596 in-frame deletion (four de novo). No clinical or radiological differences were evident between individuals with different mutations. ITPR1 encodes an inositol 1,4,5-triphosphate-responsive calcium channel. The homo-tetrameric structure has been solved by cryoelectron microscopy. Using estimations of the degree of structural change induced by known recessive- and dominant-negative mutations in other disease-associated multimeric channels, we developed a generalizable computational approach to indicate the likely mutational mechanism. This analysis supports a dominant-negative mechanism for GS variants in ITPR1. In GS-derived lymphoblastoid cell lines (LCLs), the proportion of ITPR1-positive cells using immunofluorescence was significantly higher in mutant than control LCLs, consistent with an abnormality of nuclear calcium signaling feedback control. Super-resolution imaging supports the existence of an ITPR1-lined nucleoplasmic reticulum. Mice with Itpr1 heterozygous null mutations showed no major iris defects. Purkinje cells of the cerebellum appear to be the most sensitive to impaired ITPR1 function in humans. Iris hypoplasia is likely to result from either complete loss of ITPR1 activity or structure-specific disruption of multimeric interactions.

The rapidly increasing amount of data on human genetic variation has resulted in a growing demand to identify pathogenic mutations computationally, as their experimental validation is currently beyond reach. Here we show that alpha helices and beta strands differ significantly in their ability to tolerate mutations: helices can accumulate more mutations than strands without change, due to the higher numbers of inter-residue contacts in helices. This results in two patterns: a) the same number of mutations causes less structural change in helices than in strands; b) helices diverge more rapidly in sequence than strands within the same domains. Additionally, both helices and strands are significantly more robust than coils. Based on this observation we show that human missense mutations that change secondary structure are more likely to be pathogenic than those that do not. Moreover, inclusion of predicted secondary structure changes shows significant utility for improving upon state-of-the-art pathogenicity predictions.

The sensing of nucleic acids by receptors of the innate immune system is a key component of antimicrobial immunity. RNA:DNA hybrids, as essential intracellular replication intermediates generated during infection, could therefore represent a class of previously uncharacterised pathogen-associated molecular patterns sensed by pattern recognition receptors. Here we establish that RNA:DNA hybrids containing viral-derived sequences efficiently induce pro-inflammatory cytokine and antiviral type I interferon production in dendritic cells. We demonstrate that MyD88-dependent signalling is essential for this cytokine response and identify TLR9 as a specific sensor of RNA:DNA hybrids. Hybrids therefore represent a novel molecular pattern sensed by the innate immune system and so could play an important role in host response to viruses and the pathogenesis of autoimmune disease.

Transposases are important tools in genome engineering, and there is considerable interest in engineering more efficient ones. Here, we seek to understand the factors determining their activity using the Sleeping Beauty transposase. Recent work suggests that protein coevolutionary information can be used to classify groups of physically connected, coevolving residues into elements called "sectors", which have proven useful for understanding the folding, allosteric interactions, and enzymatic activity of proteins. Using extensive mutagenesis data, protein modeling and analysis of folding energies, we show that (i) The Sleeping Beauty transposase contains two sectors, which span across conserved domains, and are enriched in DNA-binding residues, indicating that the DNA binding and endonuclease functions of the transposase coevolve; (ii) Sector residues are highly sensitive to mutations, and most mutations of these residues strongly reduce transposition rate; (iii) Mutations with a strong effect on free energy of folding in the DDE domain of the transposase significantly reduce transposition rate. (iv) Mutations that influence DNA and protein-protein interactions generally reduce transposition rate, although most hyperactive mutants are also located on the protein surface, including residues with protein-protein interactions. This suggests that hyperactivity results from the modification of protein interactions, rather than the stabilization of protein fold.

Late-onset retinal degeneration (L-ORD) is a rare autosomal dominant retinal dystrophy, characterised by extensive sub-retinal pigment epithelium (RPE) deposits, RPE atrophy, choroidal neovascularisation and photoreceptor cell death associated with severe visual loss. L-ORD shows striking phenotypic similarities to age-related macular degeneration (AMD), a common and genetically complex disorder, which can lead to misdiagnosis in the early stages. To date, a single missense mutation (S163R) in the C1QTNF5 gene, encoding C1q And Tumor Necrosis Factor Related Protein 5 (C1QTNF5) has been shown to cause L-ORD in a subset of affected families. Here, we describe the identification and characterisation of three novel pathogenic mutations in C1QTNF5 in order to elucidate disease mechanisms. In silico and in vitro characterisation show that these mutations perturb protein folding, assembly or polarity of secretion of C1QTNF5 and, importantly, all appear to destabilise the wildtype protein in co-transfection experiments in a human RPE cell line. This suggests that the heterozygous mutations in L-ORD show a dominant negative, rather than a haploinsufficient, disease mechanism. The function of C1QTNF5 remains unclear but this new insight into the pathogenetic basis of L-ORD has implications for future therapeutic strategies such as gene augmentation therapy.

Trio-based whole-exome sequence (WES) data have established confident genetic diagnoses in ∼40% of previously undiagnosed individuals recruited to the Deciphering Developmental Disorders (DDD) study. Here we aim to use the breadth of phenotypic information recorded in DDD to augment diagnosis and disease variant discovery in probands. Median Euclidean distances (mEuD) were employed as a simple measure of similarity of quantitative phenotypic data within sets of ≥10 individuals with plausibly causative de novo mutations (DNM) in 28 different developmental disorder genes. 13/28 (46.4%) showed significant similarity for growth or developmental milestone metrics, 10/28 (35.7%) showed similarity in HPO term usage, and 12/28 (43%) showed no phenotypic similarity. Pairwise comparisons of individuals with high-impact inherited variants to the 32 individuals with causative DNM in ANKRD11 using only growth z-scores highlighted 5 likely causative inherited variants and two unrecognized DNM resulting in an 18% diagnostic uplift for this gene. Using an independent approach, naive Bayes classification of growth and developmental data produced reasonably discriminative models for the 24 DNM genes with sufficiently complete data. An unsupervised naive Bayes classification of 6,993 probands with WES data and sufficient phenotypic information defined 23 in silico syndromes (ISSs) and was used to test a "phenotype first" approach to the discovery of causative genotypes using WES variants strictly filtered on allele frequency, mutation consequence, and evidence of constraint in humans. This highlighted heterozygous de novo nonsynonymous variants in SPTBN2 as causative in three DDD probands.

Approximately half of proteins with experimentally determined structures can interact with other copies of themselves and assemble into homomeric complexes, the overwhelming majority of which (>96%) are symmetric. Although homomerisation is often assumed to a functionally beneficial result of evolutionary selection, there has been little systematic analysis of the relationship between homomer structure and function. Here, utilizing the large numbers of structures and functional annotations now available, we have investigated how proteins that assemble into different types of homomers are associated with different biological functions. We observe that homomers from different symmetry groups are significantly enriched in distinct functions, and can often provide simple physical and geometrical explanations for these associations in regards to substrate recognition or physical environment. One of the strongest associations is the tendency for metabolic enzymes to form dihedral complexes, which we suggest is closely related to allosteric regulation. We provide a physical explanation for why allostery is related to dihedral complexes: it allows for efficient propagation of conformational changes across isologous (i.e. symmetric) interfaces. Overall we demonstrate a clear relationship between protein function and homomer symmetry that has important implications for understanding protein evolution, as well as for predicting protein function and quaternary structure.

High temperature requirement A1 (HtrA1) belongs to an ancient protein family that is linked to various human disorders. The precise role of exon 1-encoded N-terminal domains and how these influence the biological functions of human HtrA1 remain elusive. In this study, we traced the evolutionary origins of these N-terminal domains to a single gene fusion event in the most recent common ancestor of vertebrates. We hypothesized that human HtrA1 is implicated in unfolded protein response. In highly secretory cells of the retinal pigmented epithelia, endoplasmic reticulum (ER) stress upregulated HtrA1. HtrA1 co-localized with vimentin intermediate filaments in highly arborized fashion. Upon ER stress, HtrA1 tracked along intermediate filaments, which collapsed and bundled in an aggresome at the microtubule organizing center. Gene silencing of HtrA1 altered the schedule and amplitude of adaptive signaling and concomitantly resulted in apoptosis. Restoration of wild-type HtrA1, but not its protease inactive mutant, was necessary and sufficient to protect from apoptosis. A variant of HtrA1 that harbored exon 1 substitutions displayed reduced efficacy in rescuing cells from proteotoxicity. Our results illuminate the integration of HtrA1 in the toolkit of mammalian cells against protein misfolding and the implications of defects in HtrA1 in proteostasis.

DNA replication is fundamental for cell proliferation in all organisms. Nonetheless, components of the replisome have been implicated in human disease, and here we report PRIM1 encoding the catalytic subunit of DNA primase as a novel disease gene. Using a variant classification agnostic approach, biallelic mutations in PRIM1 were identified in five individuals. PRIM1 protein levels were markedly reduced in patient cells, accompanied by replication fork asymmetry, increased interorigin distances, replication stress, and prolonged S-phase duration. Consequently, cell proliferation was markedly impaired, explaining the patients' extreme growth failure. Notably, phenotypic features distinct from those previously reported with DNA polymerase genes were evident, highlighting differing developmental requirements for this core replisome component that warrant future investigation.

The structure of ligand-binding sites has been shown to profoundly influence the evolution of function in homomeric protein complexes. Complexes with multichain binding sites (MBSs) have more conserved quaternary structure, more similar binding sites and ligands between homologs, and evolve new functions slower than homomers with single-chain binding sites (SBSs). Here, using in silico analyses of protein dynamics, we investigate whether ligand-binding-site structure shapes allosteric signal transduction pathways, and whether the structural similarity of binding sites influences the evolution of allostery. Our analyses show that: 1) allostery is more frequent among MBS complexes than in SBS complexes, particularly in homomers; 2) in MBS homomers, semirigid communities and critical residues frequently connect interfaces and thus they are characterized by signal transduction pathways that cross protein-protein interfaces, whereas SBS homomers usually not; 3) ligand binding alters community structure differently in MBS and SBS homomers; and 4) except MBS homomers, allosteric proteins are more likely to have homologs with similar binding site than nonallosteric proteins, suggesting that binding site similarity is an important factor driving the evolution of allostery.

Protein complex assembly often occurs while subunits are being translated, resulting in complexes whose subunits were translated from the same mRNA in an allele-specific manner. It has thus been hypothesized that such cotranslational assembly may counter the assembly-mediated dominant-negative effect, whereby co-assembly of mutant and wild-type subunits "poisons" complex activity. Here, we show that cotranslationally assembling subunits are much less likely to be associated with autosomal dominant relative to recessive disorders, and that subunits with dominant-negative disease mutations are significantly depleted in cotranslational assembly compared to those associated with loss-of-function mutations. We also find that complexes with known dominant-negative effects tend to expose their interfaces late during translation, lessening the likelihood of cotranslational assembly. Finally, by combining complex properties with other features, we trained a computational model for predicting proteins likely to be associated with non-loss-of-function disease mechanisms, which we believe will be of considerable utility for protein variant interpretation.

Ligand binding site structure has profound consequences for the evolution of function of protein complexes, particularly in homomers-complexes comprising multiple copies of the same protein. Previously, we have shown that homomers with multichain binding sites (MBSs) are characterized by more conserved binding sites and quaternary structure, and qualitatively different allosteric pathways than homomers with single-chain binding sites (SBSs) or monomers. Here, using computational methods, we show that the folds of single-domain MBS and SBS homomers are different, and SBS homomers are likely to be folded cotranslationally, while MBS homomers are more likely to form post-translationally and rely on more advanced folding-assistance and quality control mechanisms, which include chaperonins. In addition, our findings demonstrate that MBS homomers are qualitatively different from monomers, while SBS homomers are much less distinct, supporting the hypothesis that the evolution of quaternary structure in SBS homomers is significantly influenced by stochastic processes.

RNase H2 is a holoenzyme, composed of 3 subunits (ribonuclease H2 subunits A, B, and C), that cleaves RNA:DNA hybrids and removes mis-incorporated ribonucleotides from genomic DNA through ribonucleotide excision repair. Ribonucleotide incorporation by eukaryotic DNA polymerases occurs during every round of genome duplication and produces the most frequent type of naturally occurring DNA lesion. We investigated whether intestinal epithelial proliferation requires RNase H2 function and whether RNase H2 activity is disrupted during intestinal carcinogenesis.

The process of resonance assignment is fundamental to most NMR studies of protein structure and dynamics. Unfortunately, the manual assignment of residues is tedious and time-consuming, and can represent a significant bottleneck for further characterization. Furthermore, while automated approaches have been developed, they are often limited in their accuracy, particularly for larger proteins. Here, we address this by introducing the software COMPASS, which, by combining automated resonance assignment with manual intervention, is able to achieve accuracy approaching that from manual assignments at greatly accelerated speeds. Moreover, by including the option to compensate for isotope shift effects in deuterated proteins, COMPASS is far more accurate for larger proteins than existing automated methods. COMPASS is an open-source project licensed under GNU General Public License and is available for download from http://www.liu.se/forskning/foass/tidigare-foass/patrik-lundstrom/software?l=en. Source code and binaries for Linux, Mac OS X and Microsoft Windows are available.

Genome integrity is continuously challenged by the DNA damage that arises during normal cell metabolism. Biallelic mutations in the genes encoding the genome surveillance enzyme ribonuclease H2 (RNase H2) cause Aicardi-Goutières syndrome (AGS), a pediatric disorder that shares features with the autoimmune disease systemic lupus erythematosus (SLE). Here we determined that heterozygous parents of AGS patients exhibit an intermediate autoimmune phenotype and demonstrated a genetic association between rare RNASEH2 sequence variants and SLE. Evaluation of patient cells revealed that SLE- and AGS-associated mutations impair RNase H2 function and result in accumulation of ribonucleotides in genomic DNA. The ensuing chronic low level of DNA damage triggered a DNA damage response characterized by constitutive p53 phosphorylation and senescence. Patient fibroblasts exhibited constitutive upregulation of IFN-stimulated genes and an enhanced type I IFN response to the immunostimulatory nucleic acid polyinosinic:polycytidylic acid and UV light irradiation, linking RNase H2 deficiency to potentiation of innate immune signaling. Moreover, UV-induced cyclobutane pyrimidine dimer formation was markedly enhanced in ribonucleotide-containing DNA, providing a mechanism for photosensitivity in RNase H2-associated SLE. Collectively, our findings implicate RNase H2 in the pathogenesis of SLE and suggest a role of DNA damage-associated pathways in the initiation of autoimmunity.

It has been suggested that the evolution of protein complexes is significantly influenced by stochastic, non-adaptive processes. Using ligand binding as a proxy of function, we show that the structure of ligand-binding sites significantly influences the evolution of protein complexes. We show that homomers with multi-chain binding sites (MBSs) evolve new functions slower than monomers or other homomers, and those binding cofactors and metals have more conserved quaternary structure than other homomers. Moreover, the ligands and ligand-binding pockets of homologous MBS homomers are more similar than monomers and other homomers. Our results suggest strong evolutionary selection for quaternary structure in cofactor-binding MBS homomers, whereas neutral processes are more important in complexes with single-chain binding sites. They also have pharmacological implications, suggesting that complexes with single-chain binding sites are better targets for selective drugs, whereas MBS homomers are good candidates for broad-spectrum antibiotic and multitarget drug design.

Protein complexes are assemblies of subunits that have co-evolved to execute one or many coordinated functions in the cellular environment. Functional annotation of mammalian protein complexes is critical to understanding biological processes, as well as disease mechanisms. Here, we used genetic co-essentiality derived from genome-scale RNAi- and CRISPR-Cas9-based fitness screens performed across hundreds of human cancer cell lines to assign measures of functional similarity. From these measures, we systematically built and characterized functional similarity networks that recapitulate known structural and functional features of well-studied protein complexes and resolve novel functional modules within complexes lacking structural resolution, such as the mammalian SWI/SNF complex. Finally, by integrating functional networks with large protein-protein interaction networks, we discovered novel protein complexes involving recently evolved genes of unknown function. Taken together, these findings demonstrate the utility of genetic perturbation screens alone, and in combination with large-scale biophysical data, to enhance our understanding of mammalian protein complexes in normal and disease states.

Ribonuclease H2 (RNase H2) is the major nuclear enzyme involved in the degradation of RNA/DNA hybrids and removal of ribonucleotides misincorporated in genomic DNA. Mutations in each of the three RNase H2 subunits have been implicated in a human auto-inflammatory disorder, Aicardi-Goutières Syndrome (AGS). To understand how mutations impact on RNase H2 function we determined the crystal structure of the human heterotrimer. In doing so, we correct several key regions of the previously reported murine RNase H2 atomic model and provide biochemical validation for our structural model. Our results provide new insights into how the subunits are arranged to form an enzymatically active complex. In particular, we establish that the RNASEH2A C terminus is a eukaryotic adaptation for binding the two accessory subunits, with residues within it required for enzymatic activity. This C-terminal extension interacts with the RNASEH2C C terminus and both are necessary to form a stable, enzymatically active heterotrimer. Disease mutations cluster at this interface between all three subunits, destabilizing the complex and/or impairing enzyme activity. Altogether, we locate 25 out of 29 residues mutated in AGS patients, establishing a firm basis for future investigations into disease pathogenesis and function of the RNase H2 enzyme.

Ribonuclease H2 is the major nuclear enzyme degrading cellular RNA/DNA hybrids in eukaryotes and the sole nuclease known to be able to hydrolyze ribonucleotides misincorporated during genomic replication. Mutation in RNASEH2 causes Aicardi-Goutières syndrome, an auto-inflammatory disorder that may arise from nucleic acid byproducts generated during DNA replication. Here, we report the crystal structures of Archaeoglobus fulgidus RNase HII in complex with PCNA, and human PCNA bound to a C-terminal peptide of RNASEH2B. In the archaeal structure, three binding modes are observed as the enzyme rotates about a flexible hinge while anchored to PCNA by its PIP-box motif. PCNA binding promotes RNase HII activity in a hinge-dependent manner. It enhances both cleavage of ribonucleotides misincorporated in DNA duplexes, and the comprehensive hydrolysis of RNA primers formed during Okazaki fragment maturation. In addition, PCNA imposes strand specificity on enzyme function, and by localizing RNase H2 and not RNase H1 to nuclear replication foci in vivo it ensures that RNase H2 is the dominant RNase H activity during nuclear replication. Our findings provide insights into how type 2 RNase H activity is directed during genome replication and repair, and suggest a mechanism by which RNase H2 may suppress generation of immunostimulatory nucleic acids.