Candidate genes and sequence variants for susceptibility to mycobacterial infection identified by whole-exome sequencing.

Inborn errors of immunity are known to influence susceptibility to mycobacterial infections. The aim of this study was to characterize the genetic profile of nine patients with mycobacterial infections (eight with BCGitis and one with disseminated tuberculosis) from the Republic of Moldova using whole-exome sequencing. In total, 12 variants in eight genes known to be associated with Mendelian Susceptibility to Mycobacterial Disease (MSMD) were detected in six out of nine patients examined. In particular, a novel splice site mutation c.373-2A>C in STAT1 gene was found and functionally confirmed in a patient with disseminated tuberculosis. Trio analysis was possible for seven out of nine patients, and resulted in 23 candidate variants in 15 novel genes. Four of these genes - GBP2, HEATR3, PPP1R9B and KDM6A were further prioritized, considering their elevated expression in immune-related tissues. Compound heterozygosity was found in GBP2 in a single patient, comprising a maternally inherited missense variant c.412G>A/p.(Ala138Thr) predicted to be deleterious and a paternally inherited intronic mutation c.1149+14T>C. Functional studies demonstrated that the intronic mutation affects splicing and the level of transcript. Finally, we analyzed pathogenicity of variant combinations in gene pairs and identified five patients with putative oligogenic inheritance. In summary, our study expands the spectrum of genetic variation contributing to susceptibility to mycobacterial infections in children and provides insight into the complex/oligogenic disease-causing mode.

DeepCLIP: predicting the effect of mutations on protein-RNA binding with deep learning.

Nucleotide variants can cause functional changes by altering protein-RNA binding in various ways that are not easy to predict. This can affect processes such as splicing, nuclear shuttling, and stability of the transcript. Therefore, correct modeling of protein-RNA binding is critical when predicting the effects of sequence variations. Many RNA-binding proteins recognize a diverse set of motifs and binding is typically also dependent on the genomic context, making this task particularly challenging. Here, we present DeepCLIP, the first method for context-aware modeling and predicting protein binding to RNA nucleic acids using exclusively sequence data as input. We show that DeepCLIP outperforms existing methods for modeling RNA-protein binding. Importantly, we demonstrate that DeepCLIP predictions correlate with the functional outcomes of nucleotide variants in independent wet lab experiments. Furthermore, we show how DeepCLIP binding profiles can be used in the design of therapeutically relevant antisense oligonucleotides, and to uncover possible position-dependent regulation in a tissue-specific manner. DeepCLIP is freely available as a stand-alone application and as a webtool at http://deepclip.compbio.sdu.dk.

Down-regulation of CK2α correlates with decreased expression levels of DNA replication minichromosome maintenance protein complex (MCM) genes.

Protein kinase CK2 is a serine/threonine kinase composed of two catalytic subunits (CK2α and/or CK2α') and two regulatory subunits (CK2ß). It is implicated in every stage of the cell cycle and in the regulation of various intracellular pathways associated with health and disease states. The catalytic subunits have similar biochemical activity, however, their functions may differ significantly in cells and in vivo. In this regard, homozygous deletion of CK2α leads to embryonic lethality in mid-gestation potentially due to severely impaired cell proliferation. To determine the CK2α-dependent molecular mechanisms that control cell proliferation, we established a myoblast-derived cell line with inducible silencing of CK2α and carried out a comprehensive RNA-Seq analysis of gene expression. We report evidence that CK2α depletion causes delayed cell cycle progression through the S-phase and defective response to replication stress. Differential gene expression analysis revealed that the down-regulated genes were enriched in pathways implicated in cell cycle regulation, DNA replication and DNA damage repair. Interestingly, the genes coding for the minichromosome maintenance proteins (MCMs), which constitute the core of the replication origin recognition complex, were among the most significantly down-regulated genes. These findings were validated in cells and whole mouse embryos. Taken together, our study provides new evidence for a critical role of protein kinase CK2 in controlling DNA replication initiation and the expression levels of replicative DNA helicases, which ensure maintenance of proliferative potential and genome integrity in eukaryotic cells.

A pathogenic haplotype, common in Europeans, causes autosomal recessive albinism and uncovers missing heritability in OCA1.

Oculocutaneous albinism (OCA) is a genetically heterogeneous disorder. Six genes are associated with autosomal recessive OCA (TYR, OCA2, TYRP1, SLC45A2, SLC24A5 and LRMDA), and one gene, GPR143, is associated with X-linked ocular albinism (OA). Molecular genetic analysis provides a genetic diagnosis in approximately 60% of individuals with clinical OA/OCA. A considerably number of the remaining 40% are heterozygous for a causative sequence variation in TYR. To identify missing causative sequence variants in these, we used a NGS based approach, genotyping and segregation analysis. We report two putative pathogenic haplotypes which only differ by two extremely rare SNVs, indicating that the haplotypes have a common derivation. Both haplotypes segregate consistent with an autosomal recessive inheritance pattern and include the allele p.S192Y-p.R402Q. An explanation for the pathogenicity of the haplotypes could be the combination of p.S192Y and p.R402Q. Homozygosity for the pathogenic haplotypes causes a partial albinism phenotype. In our cohort, 15% of affected individuals had a molecular genetic diagnosis involving the pathogenic haplotype. Consequently, the prevalence of albinism seems to be substantially underestimated, and children with unexplained bilateral subnormal vision and/or nystagmus should be analysed clinically and molecularly for albinism.

Blocking of an intronic splicing silencer completely rescues IKBKAP exon 20 splicing in familial dysautonomia patient cells.

Familial dysautonomia (FD) is a severe genetic disorder causing sensory and autonomic dysfunction. It is predominantly caused by a c.2204+6T>C mutation in the IKBKAP gene. This mutation decreases the 5' splice site strength of IKBKAP exon 20 leading to exon 20 skipping and decreased amounts of full-length IKAP protein. We identified a binding site for the splicing regulatory protein hnRNP A1 downstream of the IKBKAP exon 20 5'-splice site. We show that hnRNP A1 binds to this splicing regulatory element (SRE) and that two previously described inhibitory SREs inside IKBKAP exon 20 are also bound by hnRNP A1. Knockdown of hnRNP A1 in FD patient fibroblasts increases IKBKAP exon 20 inclusion demonstrating that hnRNP A1 is a negative regulator of IKBKAP exon 20 splicing. Furthermore, by mutating the SREs in an IKBKAP minigene we show that all three SREs cause hnRNP A1-mediated exon repression. We designed splice switching oligonucleotides (SSO) that blocks the intronic hnRNP A1 binding site, and demonstrate that this completely rescues splicing of IKBKAP exon 20 in FD patient fibroblasts and increases the amounts of IKAP protein. We propose that this may be developed into a potential new specific treatment of FD.

An intronic variation in SLC52A1 causes exon skipping and transient riboflavin-responsive multiple acyl-CoA dehydrogenation deficiency.

Vitamin B2, riboflavin is essential for cellular function, as it participates in a diversity of redox reactions central to human metabolism, through its role as precursor for the cofactors flavin mononucleotide (FMN) and flavin adenine dinucleotide (FAD), which are electron carriers. The electron transfer flavoprotein (ETF) and its dehydrogenase (ETFDH), uses FAD as cofactor. The ETF and ETFDH are forming the electron transport pathway for many mitochondrial flavoprotein dehydrogenases involved in fatty acid, amino acid and choline metabolism. A variation in either ETF or ETFDH causes multiple acyl-CoA dehydrogenation deficiency (MADD), but genetic variations in the riboflavin metabolism or transportation of riboflavin can also cause MADD. The most common variations are located in the riboflavin transporter 2 (RFVT2) and 3 (RFVT3), that are highly expressed in brain and intestinal tissues, respectively. Deficiency of riboflavin transporter 1 (RFVT1), encoded by the SLC52A1 gene, highly expressed in the placenta, has only been reported once. We here report a case of transient MADD, caused by a heterozygous intronic variation, c.1134+11G>A, in the SLC52A1 gene encoding RFVT1. This variation creates a binding site for the splice inhibitory hnRNP A1 protein and causes exon 4 skipping. Riboflavin deficiency and maternal malnutrition during pregnancy might have been the determining factor in the outcome of this case.

Corrigendum: Splicing factor 1 modulates dietary restriction and TORC1 pathway longevity in C. elegans.

This corrects the article DOI: 10.1038/nature20789.

Splicing factor 1 modulates dietary restriction and TORC1 pathway longevity in C. elegans.

Ageing is driven by a loss of transcriptional and protein homeostasis and is the key risk factor for multiple chronic diseases. Interventions that attenuate or reverse systemic dysfunction associated with age therefore have the potential to reduce overall disease risk in the elderly. Precursor mRNA (pre-mRNA) splicing is a fundamental link between gene expression and the proteome, and deregulation of the splicing machinery is linked to several age-related chronic illnesses. However, the role of splicing homeostasis in healthy ageing remains unclear. Here we demonstrate that pre-mRNA splicing homeostasis is a biomarker and predictor of life expectancy in Caenorhabditis elegans. Using transcriptomics and in-depth splicing analysis in young and old animals fed ad libitum or subjected to dietary restriction, we find defects in global pre-mRNA splicing with age that are reduced by dietary restriction via splicing factor 1 (SFA-1; the C. elegans homologue of SF1, also known as branchpoint binding protein, BBP). We show that SFA-1 is specifically required for lifespan extension by dietary restriction and by modulation of the TORC1 pathway components AMPK, RAGA-1 and RSKS-1/S6 kinase. We also demonstrate that overexpression of SFA-1 is sufficient to extend lifespan. Together, these data demonstrate a role for RNA splicing homeostasis in dietary restriction longevity and suggest that modulation of specific spliceosome components may prolong healthy ageing.

RNA-sequencing of a mouse-model of spinal muscular atrophy reveals tissue-wide changes in splicing of U12-dependent introns.

Spinal Muscular Atrophy (SMA) is a neuromuscular disorder caused by insufficient levels of the Survival of Motor Neuron (SMN) protein. SMN is expressed ubiquitously and functions in RNA processing pathways that include trafficking of mRNA and assembly of snRNP complexes. Importantly, SMA severity is correlated with decreased snRNP assembly activity. In particular, the minor spliceosomal snRNPs are affected, and some U12-dependent introns have been reported to be aberrantly spliced in patient cells and animal models. SMA is characterized by loss of motor neurons, but the underlying mechanism is largely unknown. It is likely that aberrant splicing of genes expressed in motor neurons is involved in SMA pathogenesis, but increasing evidence indicates that pathologies also exist in other tissues. We present here a comprehensive RNA-seq study that covers multiple tissues in an SMA mouse model. We show elevated U12-intron retention in all examined tissues from SMA mice, and that U12-dependent intron retention is induced upon siRNA knock-down of SMN in HeLa cells. Furthermore, we show that retention of U12-dependent introns is mitigated by ASO treatment of SMA mice and that many transcriptional changes are reversed. Finally, we report on missplicing of several Ca2+ channel genes that may explain disrupted Ca2+ homeostasis in SMA and activation of Cdk5.

The prevalent deep intronic c. 639+919 G>A GLA mutation causes pseudoexon activation and Fabry disease by abolishing the binding of hnRNPA1 and hnRNP A2/B1 to a splicing silencer.

Fabry disease is an X-linked recessive inborn disorder of the glycosphingolipid metabolism, caused by total or partial deficiency of the lysosomal α-galactosidase A enzyme due to mutations in the GLA gene. The prevalent c.639+919 G>A mutation in GLA leads to pathogenic insertion of a 57bp pseudoexon sequence from intron 4, which is responsible for the cardiac variant phenotype. In this study we investigate the splicing regulatory mechanism leading to GLA pseudoexon activation. Splicing analysis of GLA minigenes revealed that pseudoexon activation is influenced by cell-type. We demonstrate that the wild-type sequence harbors an hnRNP A1 and hnRNP A2/B1-binding exonic splicing silencer (ESS) overlapping the 5'splice site (5'ss) that prevents pseudoexon inclusion. The c.639+919 G>A mutation disrupts this ESS allowing U1 snRNP recognition of the 5'ss. We show that the wild-type GLA 5'ss motif with the ESS is also able to inhibit inclusion of an unrelated pseudoexon in the FGB gene, and that also in the FGB context inactivation of the ESS by the c.639+919 G>A mutation causes pseudoexon activation, underscoring the universal nature of the ESS. Finally, we demonstrate that splice switching oligonucleotide (SSO) mediated blocking of the pseudoexon 3'ss and 5'ss effectively restores normal GLA splicing. This indicates that SSO based splicing correction may be a therapeutic alternative in the treatment of Fabry disease.

Global identification of hnRNP A1 binding sites for SSO-based splicing modulation.

BACKGROUND: Many pathogenic genetic variants have been shown to disrupt mRNA splicing. Besides splice mutations in the well-conserved splice sites, mutations in splicing regulatory elements (SREs) may deregulate splicing and cause disease. A promising therapeutic approach is to compensate for this deregulation by blocking other SREs with splice-switching oligonucleotides (SSOs). However, the location and sequence of most SREs are not well known. RESULTS: Here, we used individual-nucleotide resolution crosslinking immunoprecipitation (iCLIP) to establish an in vivo binding map for the key splicing regulatory factor hnRNP A1 and to generate an hnRNP A1 consensus binding motif. We find that hnRNP A1 binding in proximal introns may be important for repressing exons. We show that inclusion of the alternative cassette exon 3 in SKA2 can be significantly increased by SSO-based treatment which blocks an iCLIP-identified hnRNP A1 binding site immediately downstream of the 5' splice site. Because pseudoexons are well suited as models for constitutive exons which have been inactivated by pathogenic mutations in SREs, we used a pseudoexon in MTRR as a model and showed that an iCLIP-identified hnRNP A1 binding site downstream of the 5' splice site can be blocked by SSOs to activate the exon. CONCLUSIONS: The hnRNP A1 binding map can be used to identify potential targets for SSO-based therapy. Moreover, together with the hnRNP A1 consensus binding motif, the binding map may be used to predict whether disease-associated mutations and SNPs affect hnRNP A1 binding and eventually mRNA splicing.

Splice-shifting oligonucleotide (SSO) mediated blocking of an exonic splicing enhancer (ESE) created by the prevalent c.903+469T>C MTRR mutation corrects splicing and restores enzyme activity in patient cells.

The prevalent c.903+469T>C mutation in MTRR causes the cblE type of homocystinuria by strengthening an SRSF1 binding site in an ESE leading to activation of a pseudoexon. We hypothesized that other splicing regulatory elements (SREs) are also critical for MTRR pseudoexon inclusion. We demonstrate that the MTRR pseudoexon is on the verge of being recognized and is therefore vulnerable to several point mutations that disrupt a fine-tuned balance between the different SREs. Normally, pseudoexon inclusion is suppressed by a hnRNP A1 binding exonic splicing silencer (ESS). When the c.903+469T>C mutation is present two ESEs abrogate the activity of the ESS and promote pseudoexon inclusion. Blocking the 3'splice site or the ESEs by SSOs is effective in restoring normal splicing of minigenes and endogenous MTRR transcripts in patient cells. By employing an SSO complementary to both ESEs, we were able to rescue MTRR enzymatic activity in patient cells to approximately 50% of that in controls. We show that several point mutations, individually, can activate a pseudoexon, illustrating that this mechanism can occur more frequently than previously expected. Moreover, we demonstrate that SSO blocking of critical ESEs is a promising strategy to treat the increasing number of activated pseudoexons.

The ETFDH c.158A>G variation disrupts the balanced interplay of ESE- and ESS-binding proteins thereby causing missplicing and multiple Acyl-CoA dehydrogenation deficiency.

Multiple acyl-CoA dehydrogenation deficiency is a disorder of fatty acid and amino acid oxidation caused by defects of electron transfer flavoprotein (ETF) or its dehydrogenase (ETFDH). A clear relationship between genotype and phenotype makes genotyping of patients important not only diagnostically but also for prognosis and for assessment of treatment. In the present study, we show that a predicted benign ETFDH missense variation (c.158A>G/p.Lys53Arg) in exon 2 causes exon skipping and degradation of ETFDH protein in patient samples. Using splicing reporter minigenes and RNA pull-down of nuclear proteins, we show that the c.158A>G variation increases the strength of a preexisting exonic splicing silencer (ESS) motif UAGGGA. This ESS motif binds splice inhibitory hnRNP A1, hnRNP A2/B1, and hnRNP H proteins. Binding of these inhibitory proteins prevents binding of the positive splicing regulatory SRSF1 and SRSF5 proteins to nearby and overlapping exonic splicing enhancer elements and this causes exon skipping. We further suggest that binding of hnRNP proteins to UAGGGA is increased by triggering synergistic hnRNP H binding to GGG triplets located upstream and downsteam of the UAGGGA motif. A number of disease-causing exonic elements that induce exon skipping in other genes have a similar architecture as the one in ETFDH exon 2.

A synonymous polymorphic variation in ACADM exon 11 affects splicing efficiency and may affect fatty acid oxidation.

In recent studies combining genome-wide association and tandem-MS based metabolic profiling, a single-nucleotide polymorphism (SNP), rs211718C>T, located far upstream of the MCAD gene (ACADM) was found to be associated with serum concentrations of medium-chain acylcarnitines indicating improved beta-oxidation of medium-chain fatty acids. We examined the functional basis for this association and identified linkage between rs211718 and the intragenic synonymous polymorphic variant c.1161A>G in ACADM exon 11 (rs1061337). Employing minigene studies we show that the c.1161A allele is associated with exon 11 missplicing, and that the c.1161G allele corrects this missplicing. This may result in production of more full length MCAD protein from the c.1161G allele. Our analysis suggests that the improved splicing of the c.1161G allele is due to changes in the relative binding of splicing regulatory proteins SRSF1 and hnRNP A1. Using publicly available pre-aligned RNA-seq data, we find that the ACADM c.1161G allele is expressed at significantly higher levels than the c.1161A allele across different tissues. This supports that c.1161A>G is a functional SNP, which leads to higher MCAD expression, perhaps due to improved splicing. This study is a proof of principle that synonymous SNPs are not neutral. By changing the binding sites for splicing regulatory proteins they can have significant effects on pre-mRNA splicing and thus protein function. In addition, this study shows that for a sequence variation to have an effect, it might need to change the balance in the relative binding of positive and negative splicing factors.

Influence of polymorphic OATP1B-type carriers on the disposition of docetaxel.

PURPOSE: Docetaxel is extensively metabolized by CYP3A4 in the liver but mechanisms by which the drug is taken up into hepatocytes remain poorly understood. We hypothesized that (i) liver uptake of docetaxel is mediated by the polymorphic solute carriers OATP1B1 and OATP1B3 and (ii) inherited genetic defects in this process may impair systemic drug elimination. EXPERIMENTAL DESIGN: Transport of docetaxel was studied in vitro using various cell lines stably transfected with OATP1B1*1A (wild-type), OATP1B1*5 [c.521T>C (V174A); rs4149056], OATP1B3, or the mouse transporter Oatp1b2. Docetaxel clearance was evaluated in wild-type and Oatp1b2-knockout mice as well as in two cohorts of patients with multiple variant transporter genotypes (n = 213). RESULTS: Docetaxel was found to be a substrate for OATP1B1, OATP1B3, and Oatp1b2 but was not transported by OATP1B1*5. Deficiency of Oatp1b2 in mice was associated with an 18-fold decrease in docetaxel clearance (P = 0.0099), which was unrelated to changes in intrinsic metabolic capacity in mouse liver microsomes. In patients, however, none of the studied common reduced function variants in OATP1B1 or OATP1B3 were associated with docetaxel clearance (P > 0.05). CONCLUSIONS: The existence of at least two potentially redundant uptake transporters in the human liver with similar affinity for docetaxel supports the possibility that functional defects in both of these proteins may be required to confer substantially altered disposition phenotypes. In view of the established exposure-toxicity relationships for docetaxel, we suggest that caution is warranted if docetaxel has to be administered together with agents that potently inhibit both OATP1B1 and OATP1B3.

Adeno-associated virus-delivered polycistronic microRNA-clusters for knockdown of vascular endothelial growth factor in vivo.

BACKGROUND: Vascular endothelial growth factor (VEGF) is an angiogenic growth factor that plays a critical role in several diseases, including cancer, rheumatoid arthritis and diseases of the eye. Persistent regulation of VEGF by expression of small interfering RNAs targeting VEGF represents a potential future strategy for treatment of such diseases. As a step toward this goal, the present study combines the potency of VEGF-targeted miRNA mimics, produced from a miRNA cluster, with delivery by adeno-associated virus (AAV)-based vectors. METHODS: Nine different engineered tri-cistronic miRNA clusters encoding anti-VEGF effectors were generated and tested in adult human retinal pigment epithelial (ARPE-19) cells using Renilla luciferase screening, quantitative reverse transcriptase-polymerase chain reaction (RT-PCR), western blotting and immunostaining analysis. In vivo efficacy was tested by the injection of scAAV2/8 vectors expressing the most effective miRNA cluster into murine hindlimb muscles, followed by quantitative RT-PCR. RESULTS: Plasmids containing anti-VEGF miRNA clusters showed efficient silencing of VEGF and demonstrated a combined gene silencing effect for miRNA clusters composed of multiple miRNA-mimicked RNA interference effectors. The most potent molecule, miR-5,10,7, resulted in a knockdown of VEGF by approximately 75%. Injection of scAAV2/8 vectors expressing miR-5,10,7 into murine hindlimb muscles, resulted in a 44% reduction of endogenous VEGF. CONCLUSIONS: We have developed miRNA clusters encoding anti-VEGF effectors and shown, in a mouse model, that VEGF is efficiently down-regulated by scAAV2/8-delivered miRNA clusters, allowing potent attenuation of VEGF. These findings may contribute to the development of gene therapy based on AAV-mediated delivery of miRNA clusters.

Rare versus common variants in pharmacogenetics: SLCO1B1 variation and methotrexate disposition.

Methotrexate is used to treat autoimmune diseases and malignancies, including acute lymphoblastic leukemia (ALL). Inter-individual variation in clearance of methotrexate results in heterogeneous systemic exposure, clinical efficacy, and toxicity. In a genome-wide association study of children with ALL, we identified SLCO1B1 as harboring multiple common polymorphisms associated with methotrexate clearance. The extent of influence of rare versus common variants on pharmacogenomic phenotypes remains largely unexplored. We tested the hypothesis that rare variants in SLCO1B1 could affect methotrexate clearance and compared the influence of common versus rare variants in addition to clinical covariates on clearance. From deep resequencing of SLCO1B1 exons in 699 children, we identified 93 SNPs, 15 of which were non-synonymous (NS). Three of these NS SNPs were common, with a minor allele frequency (MAF) >5%, one had low frequency (MAF 1%-5%), and 11 were rare (MAF <1%). NS SNPs (common or rare) predicted to be functionally damaging were more likely to be found among patients with the lowest methotrexate clearance than patients with high clearance. We verified lower function in vitro of four SLCO1B1 haplotypes that were associated with reduced methotrexate clearance. In a multivariate stepwise regression analysis adjusting for other genetic and non-genetic covariates, SLCO1B1 variants accounted for 10.7% of the population variability in clearance. Of that variability, common NS variants accounted for the majority, but rare damaging NS variants constituted 17.8% of SLCO1B1's effects (1.9% of total variation) and had larger effect sizes than common NS variants. Our results show that rare variants are likely to have an important effect on pharmacogenetic phenotypes.

Treatment of deep caries lesions in adults: randomized clinical trials comparing stepwise vs. direct complete excavation, and direct pulp capping vs. partial pulpotomy.

Less invasive excavation methods have been suggested for deep caries lesions. We tested the effects of stepwise vs. direct complete excavation, 1 yr after the procedure had been carried out, in 314 adults (from six centres) who had received treatment of a tooth with deep caries. The teeth had caries lesions involving 75% or more of the dentin and were centrally randomized to stepwise or direct complete excavation. Stepwise excavation resulted in fewer pulp exposures compared with direct complete excavation [difference: 11.4%, 95% confidence interval (CI) (1.2; 21.3)]. At 1 yr of follow-up, there was a statistically significantly higher success rate with stepwise excavation, with success being defined as an unexposed pulp with sustained pulp vitality without apical radiolucency [difference: 11.7%, 95% CI (0.5; 22.5)]. In a subsequent nested trial, 58 patients with exposed pulps were randomized to direct capping or partial pulpotomy. We found no significant difference in pulp vitality without apical radiolucency between the two capping procedures after more than 1 yr [31.8% and 34.5%; difference: 2.7%, 95% CI (-22.7; 26.6)]. In conclusion, stepwise excavation decreases the risk of pulp exposure compared with direct complete excavation. In view of the poor prognosis of vital pulp treatment, a stepwise excavation approach for managing deep caries lesions is recommended.