Exonic splicing code and coordination of divalent metals in proteins.

Exonic sequences contain both protein-coding and RNA splicing information but the interplay of the protein and splicing code is complex and poorly understood. Here, we have studied traditional and auxiliary splicing codes of human exons that encode residues coordinating two essential divalent metals at the opposite ends of the Irving-Williams series, a universal order of relative stabilities of metal-organic complexes. We show that exons encoding Zn2+-coordinating amino acids are supported much less by the auxiliary splicing motifs than exons coordinating Ca2+. The handicap of the former is compensated by stronger splice sites and uridine-richer polypyrimidine tracts, except for position -3 relative to 3' splice junctions. However, both Ca2+ and Zn2+ exons exhibit close-to-constitutive splicing in multiple tissues, consistent with their critical importance for metalloprotein function and a relatively small fraction of expendable, alternatively spliced exons. These results indicate that constraints imposed by metal coordination spheres on RNA splicing have been efficiently overcome by the plasticity of exon-intron architecture to ensure adequate metalloprotein expression.

Exonic splicing code and protein binding sites for calcium.

Auxilliary splicing sequences in exons, known as enhancers (ESEs) and silencers (ESSs), have been subject to strong selection pressures at the RNA and protein level. The protein component of this splicing code is substantial, recently estimated at ∼50% of the total information within ESEs, but remains poorly understood. The ESE/ESS profiles were previously associated with the Irving-Williams (I-W) stability series for divalent metals, suggesting that the ESE/ESS evolution was shaped by metal binding sites. Here, we have examined splicing activities of exonic sequences that encode protein binding sites for Ca2+, a weak binder in the I-W affinity order. We found that predicted exon inclusion levels for the EF-hand motifs and for Ca2+-binding residues in nonEF-hand proteins were higher than for average exons. For canonical EF-hands, the increase was centred on the EF-hand chelation loop and, in particular, on Ca2+-coordinating residues, with a 1>12>3∼5>9 hierarchy in the 12-codon loop consensus and usage bias at codons 1 and 12. The same hierarchy but a lower increase was observed for noncanonical EF-hands, except for S100 proteins. EF-hand loops preferentially accumulated exon splits in two clusters, one located in their N-terminal halves and the other around codon 12. Using splicing assays and published crosslinking and immunoprecipitation data, we identify candidate trans-acting factors that preferentially bind conserved GA-rich motifs encoding negatively charged amino acids in the loops. Together, these data provide evidence for the high capacity of codons for Ca2+-coordinating residues to be retained in mature transcripts, facilitating their exon-level expansion during eukaryotic evolution.

A CRISPR and high-content imaging assay compliant with ACMG/AMP guidelines for clinical variant interpretation in ciliopathies.

Ciliopathies are a broad range of inherited developmental and degenerative diseases associated with structural or functional defects in motile or primary non-motile cilia. There are around 200 known ciliopathy disease genes and whilst genetic testing can provide an accurate diagnosis, 24-60% of ciliopathy patients who undergo genetic testing do not receive a genetic diagnosis. This is partly because following current guidelines from the American College of Medical Genetics and the Association for Molecular Pathology, it is difficult to provide a confident clinical diagnosis of disease caused by missense or non-coding variants, which account for more than one-third of cases of disease. Mutations in PRPF31 are the second most common cause of the degenerative retinal ciliopathy autosomal dominant retinitis pigmentosa. Here, we present a high-throughput high-content imaging assay providing quantitative measure of effect of missense variants in PRPF31 which meets the recently published criteria for a baseline standard in vitro test for clinical variant interpretation. This assay utilizes a new PRPF31+/- human retinal cell line generated using CRISPR gene editing to provide a stable cell line with significantly fewer cilia in which novel missense variants are expressed and characterised. We show that high-content imaging of cells expressing missense variants in a ciliopathy gene on a null background can allow characterisation of variants according to the cilia phenotype. We hope that this will be a useful tool for clinical characterisation of PRPF31 variants of uncertain significance, and can be extended to variant classification in other ciliopathies.

Sequencing era methods for identifying signatures of selection in the genome.

Insights into genetic loci which are under selection and their functional roles contribute to increased understanding of the patterns of phenotypic variation we observe today. The availability of whole-genome sequence data, for humans and other species, provides opportunities to investigate adaptation and evolution at unprecedented resolution. Many analytical methods have been developed to interrogate these large data sets and characterize signatures of selection in the genome. We review here recently developed methods and consider the impact of increased computing power and data availability on the detection of selection signatures. Consideration of demography, recombination and other confounding factors is important, and use of a range of methods in combination is a powerful route to resolving different forms of selection in genome sequence data. Overall, a substantial improvement in methods for application to whole-genome sequencing is evident, although further work is required to develop robust and computationally efficient approaches which may increase reproducibility across studies.

Understanding the disease genome: gene essentiality and the interplay of selection, recombination and mutation.

Despite the identification of many genetic variants contributing to human disease (the 'disease genome'), establishing reliable molecular diagnoses remain challenging in many cases. The ability to sequence the genomes of patients has been transformative, but difficulty in interpretation of voluminous genetic variation often confounds recognition of underlying causal variants. There are numerous predictors of pathogenicity for individual DNA variants, but their utility is reduced because many plausibly pathogenic variants are probably neutral. The rapidly increasing quantity and quality of information on the properties of genes suggests that gene-specific information might be useful for prediction of causal variation when used alongside variant-specific predictors of pathogenicity. The key to understanding the role of genes in disease relates in part to gene essentiality, which has recently been approximated, for example, by quantifying the degree of intolerance of individual genes to loss-of-function variation. Increasing understanding of the interplay between genetic recombination, selection and mutation and their relationship to gene essentiality suggests that gene-specific information may be useful for the interpretation of sequenced genomes. Considered alongside additional distinctive properties of the disease genome, such as the timing of the evolutionary emergence of genes and the roles of their products in protein networks, the case for using gene-specific measures to guide filtering of sequenced genomes seems strong.

Linkage disequilibrium maps to guide contig ordering for genome assembly.

MOTIVATION: Efforts to establish reference genome sequences by de novo sequence assembly have to address the difficulty of linking relatively short sequence contigs to form much larger chromosome assemblies. Efficient strategies are required to span gaps and establish contig order and relative orientation. We consider here the use of linkage disequilibrium (LD) maps of sequenced contigs and the utility of LD for ordering, orienting and positioning linked sequences. LD maps are readily constructed from population data and have at least an order of magnitude higher resolution than linkage maps providing the potential to resolve difficult areas in assemblies. We empirically evaluate a linkage disequilibrium map-based method using single nucleotide polymorphism genotype data in a 216 kilobase region of human 6p21.3 from which three shorter contigs are formed. RESULTS: LD map length is most informative about the correct order and orientation and is suggested by the shortest LD map where the residual error variance is close to one. For regions in strong LD this method may be less informative for correcting inverted contigs than for identifying correct contig orders. For positioning two contigs in linkage disequilibrium with each other the inter-contig distances may be roughly estimated by this method. AVAILABILITY AND IMPLEMENTATION: The LDMAP program is written in C for a linux platform and is available at https://www.soton.ac.uk/genomicinformatics/research/ld.page. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Gene-dense autosomal chromosomes show evidence for increased selection.

Purifying selection tends to reduce nucleotide and haplotype diversity leading to increased linkage disequilibrium. However, detection of evidence for selection is difficult as the signature is confounded by wide variation in the recombination rate which has a complex relationship with selection. The effective bottleneck time (the ratio of the linkage disequilibrium map to the genetic map in Morgans) controls for variability in the recombination rate. Reduced effective bottleneck times indicate stronger residual linkage disequilibrium, consistent with increased selection. Using whole genome sequence data from one European and three Sub-Saharan African human populations we find, in the African samples, strong correlations between high gene densities and reduced effective bottleneck time for autosomal chromosomes. This suggests that gene-dense autosomes have been subject to increased purifying selection reducing effective bottleneck times compared to gene-poor autosomes. Although previous studies have shown unusually strong linkage disequilibrium for the sex chromosomes variation within the autosomes has not been recognised. The strongest relationship is between effective bottleneck time and the density of essential genes, which are likely targets of greater selective pressure (p = 0.006, for the 22 autosomes). The magnitude of the reduction in chromosome-specific effective bottleneck times from the least to the most gene-dense autosomes is ~17-21% for Sub-Saharan African populations. The effect size is greater in Sub-Saharan African populations, compared to a European sample, consistent with increased efficiency of selection in populations with larger effective population sizes which have not been subject to intense population bottlenecks as experienced by populations of European ancestry. The findings highlight the value of deeper analyses of selection within Sub-Saharan African populations.

Autosomal dominant tubulointerstitial kidney disease-UMOD is the most frequent non polycystic genetic kidney disease.

BACKGROUND: Autosomal dominant tubulointerstitial kidney disease (ADTKD) caused by mutations in the UMOD gene (ADTKD-UMOD) is considered rare and often remains unrecognised. We aimed to establish the prevalence of genetic kidney diseases, ADTKD and ADTKD-UMOD in adult chronic kidney disease (CKD) patients, and to investigate characteristic features. METHODS: We sent questionnaires on family history to all patients with CKD stages 3-5 in our tertiary renal centre to identify patients with inherited renal disease. Details on clinical and family history were obtained from patient interviews and clinical records. Sanger sequencing of the UMOD gene was performed from blood or saliva samples. RESULTS: 2027 of 3770 sent questionnaires were returned. 459 patients reported a family history, which was consistent with inherited kidney disease in 217 patients. 182 non-responders with inherited kidney diseases were identified through a database search. Of these 399 individuals, 252 had autosomal dominant polycystic kidney disease (ADPKD), 28 had ADTKD, 25 had Alports, and 44 were unknown, resulting in 11% of CKD 3-5 patients and 19% of end-stage renal disease patients with genetic kidney diseases. Of the unknown, 40 were genotyped, of whom 31 had findings consistent with ADTKD. 30% of unknowns and 39% of unknowns with ADTKD had UMOD mutations. Altogether, 35 individuals from 18 families were found to have ten distinct UMOD mutations (three novel), making up 1% of patients with CKD 3-5, 2% of patients with end-stage renal disease, 9% of inherited kidney diseases and 56% with ADTKD. ADTKD-UMOD was the most common genetic kidney disease after ADPKD with a population prevalence of 9 per million. Less proteinuria and haematuria, but not hyperuricaemia or gout were predictive of ADTKD-UMOD. The main limitations of the study are the single-centre design and a predominantly Caucasian population. CONCLUSIONS: The prevalence of genetic kidney diseases and ADTKD-UMOD is significantly higher than previously described. Clinical features poorly predicted ADTKD-UMOD, highlighting the need for genetic testing guided by family history alone.

Exome sequence read depth methods for identifying copy number changes.

Copy number variants (CNVs) play important roles in a number of human diseases and in pharmacogenetics. Powerful methods exist for CNV detection in whole genome sequencing (WGS) data, but such data are costly to obtain. Many disease causal CNVs span or are found in genome coding regions (exons), which makes CNV detection using whole exome sequencing (WES) data attractive. If reliably validated against WGS-based CNVs, exome-derived CNVs have potential applications in a clinical setting. Several algorithms have been developed to exploit exome data for CNV detection and comparisons made to find the most suitable methods for particular data samples. The results are not consistent across studies. Here, we review some of the exome CNV detection methods based on depth of coverage profiles and examine their performance to identify problems contributing to discrepancies in published results. We also present a streamlined strategy that uses a single metric, the likelihood ratio, to compare exome methods, and we demonstrated its utility using the VarScan 2 and eXome Hidden Markov Model (XHMM) programs using paired normal and tumour exome data from chronic lymphocytic leukaemia patients. We use array-based somatic CNV (SCNV) calls as a reference standard to compute prevalence-independent statistics, such as sensitivity, specificity and likelihood ratio, for validation of the exome-derived SCNVs. We also account for factors known to influence the performance of exome read depth methods, such as CNV size and frequency, while comparing our findings with published results.

Mutations specific to the Rac-GEF domain of TRIO cause intellectual disability and microcephaly.

BACKGROUND: Neurodevelopmental disorders have challenged clinical genetics for decades, with over 700 genes implicated and many whose function remains unknown. The application of whole-exome sequencing is proving pivotal in closing the genotype/phenotype gap through the discovery of new genes and variants that help to unravel the pathogenic mechanisms driving neuropathogenesis. One such discovery includes TRIO, a gene recently implicated in neurodevelopmental delay. Trio is a Dbl family guanine nucleotide exchange factor (GEF) and a major regulator of neuronal development, controlling actin cytoskeleton dynamics by activating the GTPase Rac1. METHODS: Whole-exome sequencing was undertaken on a family presenting with global developmental delay, microcephaly and mild dysmorphism. Father/daughter exome analysis was performed, followed by confirmatory Sanger sequencing and segregation analysis on four individuals. Three further patients were recruited through the deciphering developmental disorders (DDD) study. Functional studies were undertaken using patient-specific Trio protein mutations. RESULTS: We identified a frameshift deletion in TRIO that segregated autosomal dominantly. By scrutinising data from DDD, we further identified three unrelated children with a similar phenotype who harboured de novo missense mutations in TRIO. Biochemical studies demonstrated that in three out of four families, the Trio mutations led to a markedly reduced Rac1 activation. CONCLUSIONS: We describe an inherited global developmental delay phenotype associated with a frameshift deletion in TRIO. Additionally, we identify pathogenic de novo missense mutations in TRIO associated with the same consistent phenotype, intellectual disability, microcephaly and dysmorphism with striking digital features. We further functionally validate the importance of the GEF domain in Trio protein function. Our study demonstrates how genomic technologies are yet again proving prolific in diagnosing and advancing the understanding of neurodevelopmental disorders.

Collagen (COL4A) mutations are the most frequent mutations underlying adult focal segmental glomerulosclerosis.

BACKGROUND: Multiple genes underlying focal segmental glomerulosclerosis (FSGS) and/or steroid-resistant nephrotic syndrome (SRNS) have been identified, with the recent inclusion of collagen IV mutations responsible for Alport disease (AD) or thin basement membrane nephropathy (TBMN). We aimed to investigate the distribution of gene mutations in adult patients with primary FSGS/SRNS by targeted next generation sequencing (NGS). METHODS: Eighty-one adults from 76 families were recruited; 24 families had a history of renal disease. A targeted NGS panel was designed and applied, covering 39 genes implicated in FSGS/SRNS including COL4A3-5. RESULTS: Confirmed pathogenic mutations were found in 10 patients (6 with family history) from 9 families (diagnostic rate 12%). Probably pathogenic mutations were identified in an additional six patients (combined diagnostic rate 20%). Definitely pathogenic mutations were identified in 22% of patients with family history and 10% without. Mutations in COL4A3-5 were present in eight patients from six families, representing 56% of definitely pathogenic mutations, and establishing a diagnosis of AD in six patients and TBMN in two patients. Collagen mutations were identified in 38% of families with familial FSGS, and 3% with sporadic FSGS, with over half the mutations occurring in COL4A5. Patients with collagen mutations were younger at presentation and more likely to have family history, haematuria and glomerular basement membrane abnormalities. CONCLUSIONS: We show that collagen IV mutations, including COL4A5, frequently underlie FSGS and should be considered, particularly with a positive family history. Targeted NGS improves diagnostic efficiency by investigating many candidate genes in parallel.

Whole genome sequences are required to fully resolve the linkage disequilibrium structure of human populations.

BACKGROUND: An understanding of linkage disequilibrium (LD) structures in the human genome underpins much of medical genetics and provides a basis for disease gene mapping and investigating biological mechanisms such as recombination and selection. Whole genome sequencing (WGS) provides the opportunity to determine LD structures at maximal resolution. RESULTS: We compare LD maps constructed from WGS data with LD maps produced from the array-based HapMap dataset, for representative European and African populations. WGS provides up to 5.7-fold greater SNP density than array-based data and achieves much greater resolution of LD structure, allowing for identification of up to 2.8-fold more regions of intense recombination. The absence of ascertainment bias in variant genotyping improves the population representativeness of the WGS maps, and highlights the extent of uncaptured variation using array genotyping methodologies. The complete capture of LD patterns using WGS allows for higher genome-wide association study (GWAS) power compared to array-based GWAS, with WGS also allowing for the analysis of rare variation. The impact of marker ascertainment issues in arrays has been greatest for Sub-Saharan African populations where larger sample sizes and substantially higher marker densities are required to fully resolve the LD structure. CONCLUSIONS: WGS provides the best possible resource for LD mapping due to the maximal marker density and lack of ascertainment bias. WGS LD maps provide a rich resource for medical and population genetics studies. The increasing availability of WGS data for large populations will allow for improved research utilising LD, such as GWAS and recombination biology studies.

Rapid identification of Saccharomyces eubayanus and its hybrids.

The Saccharomyces species Saccharomyces eubayanus was recently discovered in Patagonia. Genomic, genetic and phylogenetic data all suggest that it is one of the two parents of the hybrid yeast S. pastorianus (S. cerevisiae being the other). Saccharomyces eubayanus genomes can also be found in strains of the hybrid species S. bayanus. Here, we describe a novel pair of PCR primers targeted against the S. eubayanus FSY1 gene that will specifically detect S. eubayanus and hybrids containing this allele of the gene. The primer pair can be used to identify the species using a rapid, inexpensive colony PCR method suitable for a preliminary identification of wild isolates.

Bridging the language gap - A call for the wider use of Human phenotype ontology by non-geneticist clinicians when requesting genomic tests.

Advances in genomic technology including the development of next-generation sequencing (NGS) have enabled the identification of thousands of variations at a time, allowing the discovery of novel genetic diseases. Given the volume of data generated by these investigations, attention is drawn towards reporting relevant clinical features by clinicians to guide the diagnosis and management of their patients. The Human Phenotype Ontology (HPO) developed in 2008, revolutionized the semantic vocabulary of phenotypic descriptions in genomic medicine allowing researchers, laboratories and clinical geneticists to better understand each other. In this era of personalized medicine where genetic tests are becoming more accessible, non-geneticist clinicians are expected to be more involved than ever in the process of ordering genetic tests and interpreting genetic reports. It is therefore essential that they understand and adequately apply HPO nomenclature to integrate the patient care chain and seize the opportunity offered by this tailored language. The current article highlights the importance of using HPO vocabularies in clinical practice and advocates for its wider use by non-geneticist clinicians. Correct use of HPO will reduce misunderstandings between healthcare professionals and ultimately improve the healthcare system.

Pharmacological EZH2 inhibition combined with retinoic acid treatment promotes differentiation and apoptosis in rhabdomyosarcoma cells.

BACKGROUND: Rhabdomyosarcomas (RMS) are predominantly paediatric sarcomas thought to originate from muscle precursor cells due to impaired myogenic differentiation. Despite intensive treatment, 5-year survival for patients with advanced disease remains low (< 30%), highlighting a need for novel therapies to improve outcomes. Differentiation therapeutics are agents that induce differentiation of cancer cells from malignant to benign. The histone methyltransferase, Enhancer of Zeste Homolog 2 (EZH2) suppresses normal skeletal muscle differentiation and is highly expressed in RMS tumours. RESULTS: We demonstrate combining inhibition of the epigenetic modulator EZH2 with the differentiating agent retinoic acid (RA) is more effective at reducing cell proliferation in RMS cell lines than single agents alone. In PAX3-FOXO1 positive RMS cells this is due to an RA-driven induction of the interferon pathway resulting in apoptosis. In fusion negative RMS, combination therapy led to an EZH2i-driven upregulation of myogenic signalling resulting in differentiation. In both subtypes, EZH2 is significantly associated with enrichment of trimethylated lysine 27 on histone 3 (H3K27me3) in genes that are downregulated in untreated RMS cells and upregulated with EZH2 inhibitor treatment. These results provide insight into the mechanism that drives the anti-cancer effect of the EZH2/RA single agent and combination treatment and indicate that the reduction of EZH2 activity combined with the induction of RA signalling represents a potential novel therapeutic strategy to treat both subtypes of RMS. CONCLUSIONS: The results of this study demonstrate the potential utility of combining EZH2 inhibitors with differentiation agents for the treatment of paediatric rhabdomyosarcomas. As EZH2 inhibitors are currently undergoing clinical trials for adult and paediatric solid tumours and retinoic acid differentiation agents are already in clinical use this presents a readily translatable potential therapeutic strategy. Moreover, as inhibition of EZH2 in the poor prognosis FPRMS subtype results in an inflammatory response, it is conceivable that this strategy may also synergise with immunotherapies for a more effective treatment in these patients.

The Palestinian primary ciliary dyskinesia population: first results of the diagnostic and genetic spectrum.

Background: Diagnostic testing for primary ciliary dyskinesia (PCD) started in 2013 in Palestine. We aimed to describe the diagnostic, genetic and clinical spectrum of the Palestinian PCD population. Methods: Individuals with symptoms suggestive of PCD were opportunistically considered for diagnostic testing: nasal nitric oxide (nNO) measurement, transmission electron microscopy (TEM) and/or PCD genetic panel or whole-exome testing. Clinical characteristics of those with a positive diagnosis were collected close to testing including forced expiratory volume in 1 s (FEV1) Global Lung Index z-scores and body mass index z-scores. Results: 68 individuals had a definite positive PCD diagnosis, 31 confirmed by genetic and TEM results, 23 by TEM results alone, and 14 by genetic variants alone. 45 individuals from 40 families had 17 clinically actionable variants and four had variants of unknown significance in 14 PCD genes. CCDC39, DNAH11 and DNAAF11 were the most commonly mutated genes. 100% of variants were homozygous. Patients had a median age of 10.0 years at diagnosis, were highly consanguineous (93%) and 100% were of Arabic descent. Clinical features included persistent wet cough (99%), neonatal respiratory distress (84%) and situs inversus (43%). Lung function at diagnosis was already impaired (FEV1 z-score median -1.90 (-5.0-1.32)) and growth was mostly within the normal range (z-score mean -0.36 (-3.03-2.57). 19% individuals had finger clubbing. Conclusions: Despite limited local resources in Palestine, detailed geno- and phenotyping forms the basis of one of the largest national PCD populations globally. There was notable familial homozygosity within the context of significant population heterogeneity.

Comparison of Mendeliome exome capture kits for use in clinical diagnostics.

Next generation sequencing has disrupted genetic testing, allowing far more scope in the tests applied. The appropriate sections of the genome to be tested can now be readily selected, from single mutations to whole-genome sequencing. One product offering within this spectrum are focused exomes, targeting ~5,000 genes know to be implicated in human disease. These are designed to offer a flexible platform offering high diagnostic yield with a reduction in sequencing requirement compared to whole exome sequencing. Here, we have undertaken sequencing of control DNA samples and compare two kits, the Illumina TruSight One and the Agilent SureSelect Focused Exome. Characteristics of the kits are comprehensively evaluated. Despite the larger design region of the Agilent kit, we find that the Illumina kit performs better in terms of gene coverage, as well as coverage of clinically relevant loci. We provide exhaustive coverage statistics for each kit to aid the assessment of their suitability and provide read data for control DNA samples to allow for bioinformatic benchmarking by users developing pipelines for these data.

Cold-induced urticarial autoinflammatory syndrome related to factor XII activation.

Hereditary autoinflammatory diseases are caused by gene mutations of the innate immune pathway, e.g. nucleotide receptor protein 3 (NLRP3). Here, we report a four-generation family with cold-induced urticarial rash, arthralgia, chills, headache and malaise associated with an autosomal-dominant inheritance. Genetic studies identify a substitution mutation in gene F12 (T859A, resulting in p.W268R) which encodes coagulation factor XII (FXII). Functional analysis reveals enhanced autocatalytic cleavage of the mutated protein and spontaneous FXII activation in patient plasma and in supernatant of transfected HEK293 cells expressing recombinant W268R-mutated proteins. Furthermore, we observe reduced plasma prekallikrein, cleaved high molecular weight kininogen and elevated plasma bradykinin. Neutrophils are identified as a local source of FXII. Interleukin-1ß (IL-1ß) is upregulated in lesional skin and mononuclear donor cells exposed to recombinant mutant proteins. Treatment with icatibant (bradykinin-B2-antagonist) or anakinra (interleukin-1-antagonist) reduces disease activity in patients. In conclusion, our findings provide a link between contact system activation and cytokine-mediated inflammation.