The debranching enzyme Dbr1 regulates lariat turnover and intron splicing.

The majority of genic transcription is intronic. Introns are removed by splicing as branched lariat RNAs which require rapid recycling. The branch site is recognized during splicing catalysis and later debranched by Dbr1 in the rate-limiting step of lariat turnover. Through generation of a viable DBR1 knockout cell line, we find the predominantly nuclear Dbr1 enzyme to encode the sole debranching activity in human cells. Dbr1 preferentially debranches substrates that contain canonical U2 binding motifs, suggesting that branchsites discovered through sequencing do not necessarily represent those favored by the spliceosome. We find that Dbr1 also exhibits specificity for particular 5' splice site sequences. We identify Dbr1 interactors through co-immunoprecipitation mass spectrometry. We present a mechanistic model for Dbr1 recruitment to the branchpoint through the intron-binding protein AQR. In addition to a 20-fold increase in lariats, Dbr1 depletion increases exon skipping. Using ADAR fusions to timestamp lariats, we demonstrate a defect in spliceosome recycling. In the absence of Dbr1, spliceosomal components remain associated with the lariat for a longer period of time. As splicing is co-transcriptional, slower recycling increases the likelihood that downstream exons will be available for exon skipping.

The debranching enzyme Dbr1 regulates lariat turnover and intron splicing.

The majority of genic transcription is intronic. Introns are removed by splicing as branched lariat RNAs which require rapid recycling. The branch site is recognized during splicing catalysis and later debranched by Dbr1 in the rate-limiting step of lariat turnover. Through generation of the first viable DBR1 knockout cell line, we find the predominantly nuclear Dbr1 enzyme to encode the sole debranching activity in human cells. Dbr1 preferentially debranches substrates that contain canonical U2 binding motifs, suggesting that branchsites discovered through sequencing do not necessarily represent those favored by the spliceosome. We find that Dbr1 also exhibits specificity for particular 5' splice site sequences. We identify Dbr1 interactors through co-immunoprecipitation mass spectroscopy. We present a mechanistic model for Dbr1 recruitment to the branchpoint through the intron-binding protein AQR. In addition to a 20-fold increase in lariats, Dbr1 depletion increases exon skipping. Using ADAR fusions to timestamp lariats, we demonstrate a defect in spliceosome recycling. In the absence of Dbr1, spliceosomal components remain associated with the lariat for a longer period of time. As splicing is co-transcriptional, slower recycling increases the likelihood that downstream exons will be available for exon skipping.

Hereditary cancer genes are highly susceptible to splicing mutations.

Substitutions that disrupt pre-mRNA splicing are a common cause of genetic disease. On average, 13.4% of all hereditary disease alleles are classified as splicing mutations mapping to the canonical 5' and 3' splice sites. However, splicing mutations present in exons and deeper intronic positions are vastly underreported. A recent re-analysis of coding mutations in exon 10 of the Lynch Syndrome gene, MLH1, revealed an extremely high rate (77%) of mutations that lead to defective splicing. This finding is confirmed by extending the sampling to five other exons in the MLH1 gene. Further analysis suggests a more general phenomenon of defective splicing driving Lynch Syndrome. Of the 36 mutations tested, 11 disrupted splicing. Furthermore, analyzing past reports suggest that MLH1 mutations in canonical splice sites also occupy a much higher fraction (36%) of total mutations than expected. When performing a comprehensive analysis of splicing mutations in human disease genes, we found that three main causal genes of Lynch Syndrome, MLH1, MSH2, and PMS2, belonged to a class of 86 disease genes which are enriched for splicing mutations. Other cancer genes were also enriched in the 86 susceptible genes. The enrichment of splicing mutations in hereditary cancers strongly argues for additional priority in interpreting clinical sequencing data in relation to cancer and splicing.

Defective splicing of the RB1 transcript is the dominant cause of retinoblastomas.

Defective splicing is a common cause of genetic diseases. On average, 13.4% of all hereditary disease alleles are classified as splicing mutations with most mapping to the critical GT or AG nucleotides within the 5' and 3' splice sites. However, splicing mutations are underreported and the fraction of splicing mutations that compose all disease alleles varies greatly across disease gene. For example, there is a great excess (46%; ~threefold) of hereditary disease alleles that map to splice sites in RB1 that cause retinoblastoma. Furthermore, mutations in the exons and deeper intronic position may also affect splicing. We recently developed a high-throughput method that assays reported disease mutations for their ability to disrupt pre-mRNA splicing. Surprisingly, 27% of RB1-coding mutations tested also disrupt splicing. High-throughput in vitro spliceosomal assembly assay reveals heterogeneity in which stage of spliceosomal assembly is affected by splicing mutations. 58% of exonic splicing mutations were primarily blocked at the A complex in transition to the B complex and 33% were blocked at the B complex. Several mutants appear to reduce more than one step in the assembly. As RB1 splicing mutants are enriched in retinoblastoma disease alleles, additional priority should be allocated to this class of allele while interpreting clinical sequencing experiments. Analysis of the spectrum of RB1 variants observed in 60,706 exomes identifies 197 variants that have enough potential to disrupt splicing to warrant further consideration.

The effects of structure on pre-mRNA processing and stability.

Pre-mRNA molecules can form a variety of structures, and both secondary and tertiary structures have important effects on processing, function and stability of these molecules. The prediction of RNA secondary structure is a challenging problem and various algorithms that use minimum free energy, maximum expected accuracy and comparative evolutionary based methods have been developed to predict secondary structures. However, these tools are not perfect, and this remains an active area of research. The secondary structure of pre-mRNA molecules can have an enhancing or inhibitory effect on pre-mRNA splicing. An example of enhancing structure can be found in a novel class of introns in zebrafish. About 10% of zebrafish genes contain a structured intron that forms a bridging hairpin that enforces correct splice site pairing. Negative examples of splicing include local structures around splice sites that decrease splicing efficiency and potentially cause mis-splicing leading to disease. Splicing mutations are a frequent cause of hereditary disease. The transcripts of disease genes are significantly more structured around the splice sites, and point mutations that increase the local structure often cause splicing disruptions. Post-splicing, RNA secondary structure can also affect the stability of the spliced intron and regulatory RNA interference pathway intermediates, such as pre-microRNAs. Additionally, RNA secondary structure has important roles in the innate immune defense against viruses. Finally, tertiary structure can also play a large role in pre-mRNA splicing. One example is the G-quadruplex structure, which, similar to secondary structure, can either enhance or inhibit splicing through mechanisms such as creating or obscuring RNA binding protein sites.

Pathogenic variants that alter protein code often disrupt splicing.

The lack of tools to identify causative variants from sequencing data greatly limits the promise of precision medicine. Previous studies suggest that one-third of disease-associated alleles alter splicing. We discovered that the alleles causing splicing defects cluster in disease-associated genes (for example, haploinsufficient genes). We analyzed 4,964 published disease-causing exonic mutations using a massively parallel splicing assay (MaPSy), which showed an 81% concordance rate with splicing in patient tissue. Approximately 10% of exonic mutations altered splicing, mostly by disrupting multiple stages of spliceosome assembly. We present a large-scale characterization of exonic splicing mutations using a new technology that facilitates variant classification and keeps pace with variant discovery.

Chromosomal Imbalances in Patients with Congenital Cardiac Defects: A Meta-analysis Reveals Novel Potential Critical Regions Involved in Heart Development.

OBJECTIVE: Congenital cardiac defects represent the most common group of birth defects, affecting an estimated six per 1000 births. Genetic characterization of patients and families with cardiac defects has identified a number of genes required for heart development. Yet, despite the rapid pace of these advances, mutations affecting known genes still account for only a small fraction of congenital heart defects suggesting that many more genes and developmental mechanisms remain to be identified. DESIGN: In this study, we reviewed 1694 described cases of patients with cardiac defects who were determined to have a significant chromosomal imbalance (a deletion or duplication). The cases were collected from publicly available databases (DECIPHER, ISCA, and CHDWiki) and from recent publications. An additional 68 nonredundant cases were included from the University of Michigan. Cases with multiple chromosomal or whole chromosome defects (trisomy 13, 18, 21) were excluded, and cases with overlapping deletions and/or insertions were grouped to identify regions potentially involved in heart development. RESULTS: Seventy-nine chromosomal regions were identified in which 5 or more patients had overlapping imbalances. Regions of overlap were used to determine minimal critical domains most likely to contain genes or regulatory elements involved in heart development. This approach was used to refine the critical regions responsible for cardiac defects associated with chromosomal imbalances involving 1q24.2, 2q31.1, 15q26.3, and 22q11.2. CONCLUSIONS: The pattern of chromosomal imbalances in patients with congenital cardiac defects suggests that many loci may be involved in normal heart development, some with very strong and direct effects and others with less direct effects. Chromosomal duplication/deletion mapping will provide an important roadmap for genome-wide sequencing and genetic mapping strategies to identify novel genes critical for heart development.

Genetic variation and RNA binding proteins: tools and techniques to detect functional polymorphisms.

At its most fundamental level the goal of genetics is to connect genotype to phenotype. This question is asked at a basic level evaluating the role of genes and pathways in genetic model organism. Increasingly, this question is being asked in the clinic. Genomes of individuals and populations are being sequenced and compared. The challenge often comes at the stage of analysis. The variant positions are analyzed with the hope of understanding human disease. However after a genome or exome has been sequenced, the researcher is often deluged with hundreds of potentially relevant variations. Traditionally, amino-acid changing mutations were considered the tractable class of disease-causing mutations; however, mutations that disrupt noncoding elements are the subject of growing interest. These noncoding changes are a major avenue of disease (e.g., one in three hereditary disease alleles are predicted to affect splicing). Here, we review some current practices of medical genetics, the basic theory behind biochemical binding and functional assays, and then explore technical advances in how variations that alter RNA protein recognition events are detected and studied. These advances are advances in scale-high-throughput implementations of traditional biochemical assays that are feasible to perform in any molecular biology laboratory. This chapter utilizes a case study approach to illustrate some methods for analyzing polymorphisms. The first characterizes a functional intronic SNP that deletes a high affinity PTB site using traditional low-throughput biochemical and functional assays. From here we demonstrate the utility of high-throughput splicing and spliceosome assembly assays for screening large sets of SNPs and disease alleles for allelic differences in gene expression. Finally we perform three pilot drug screens with small molecules (G418, tetracycline, and valproic acid) that illustrate how compounds that rescue specific instances of differential pre-mRNA processing can be discovered.

Functionally significant, rare transcription factor variants in tetralogy of Fallot.

OBJECTIVE: Rare variants in certain transcription factors involved in cardiac development cause Mendelian forms of congenital heart disease. The purpose of this study was to systematically assess the frequency of rare transcription factor variants in sporadic patients with the cardiac outflow tract malformation tetralogy of Fallot (TOF). METHODS AND RESULTS: We sequenced the coding, 5'UTR, and 3'UTR regions of twelve transcription factor genes implicated in cardiac outflow tract development (NKX2.5, GATA4, ISL1, TBX20, MEF2C, BOP/SMYD1, HAND2, FOXC1, FOXC2, FOXH, FOXA2 and TBX1) in 93 non-syndromic, non-Mendelian TOF cases. We also analysed Illumina Human 660W-Quad SNP Array data for copy number variants in these genes; none were detected. Four of the rare variants detected have previously been shown to affect transactivation in in vitro reporter assays: FOXC1 p.P297S, FOXC2 p.Q444R, FOXH1 p.S113T and TBX1 p.P43_G61del PPPPRYDPCAAAAPGAPGP. Two further rare variants, HAND2 p.A25_A26insAA and FOXC1 p.G378_G380delGGG, A488_491delAAAA, affected transactivation in in vitro reporter assays. Each of these six functionally significant variants was present in a single patient in the heterozygous state; each of the four for which parental samples were available were maternally inherited. Thus in the 93 TOF cases we identified six functionally significant mutations in the secondary heart field transcriptional network. SIGNIFICANCE: This study indicates that rare genetic variants in the secondary heart field transcriptional network with functional effects on protein function occur in 3-13% of patients with TOF. This is the first report of a functionally significant HAND2 mutation in a patient with congenital heart disease.

An international effort towards developing standards for best practices in analysis, interpretation and reporting of clinical genome sequencing results in the CLARITY Challenge.

BACKGROUND: There is tremendous potential for genome sequencing to improve clinical diagnosis and care once it becomes routinely accessible, but this will require formalizing research methods into clinical best practices in the areas of sequence data generation, analysis, interpretation and reporting. The CLARITY Challenge was designed to spur convergence in methods for diagnosing genetic disease starting from clinical case history and genome sequencing data. DNA samples were obtained from three families with heritable genetic disorders and genomic sequence data were donated by sequencing platform vendors. The challenge was to analyze and interpret these data with the goals of identifying disease-causing variants and reporting the findings in a clinically useful format. Participating contestant groups were solicited broadly, and an independent panel of judges evaluated their performance. RESULTS: A total of 30 international groups were engaged. The entries reveal a general convergence of practices on most elements of the analysis and interpretation process. However, even given this commonality of approach, only two groups identified the consensus candidate variants in all disease cases, demonstrating a need for consistent fine-tuning of the generally accepted methods. There was greater diversity of the final clinical report content and in the patient consenting process, demonstrating that these areas require additional exploration and standardization. CONCLUSIONS: The CLARITY Challenge provides a comprehensive assessment of current practices for using genome sequencing to diagnose and report genetic diseases. There is remarkable convergence in bioinformatic techniques, but medical interpretation and reporting are areas that require further development by many groups.

Recessive truncating titin gene, TTN, mutations presenting as centronuclear myopathy.

OBJECTIVE: To identify causative genes for centronuclear myopathies (CNM), a heterogeneous group of rare inherited muscle disorders that often present in infancy or early life with weakness and hypotonia, using next-generation sequencing of whole exomes and genomes. METHODS: Whole-exome or -genome sequencing was performed in a cohort of 29 unrelated patients with clinicopathologic diagnoses of CNM or related myopathy depleted for cases with mutations of MTM1, DNM2, and BIN1. Immunofluorescence analyses on muscle biopsies, splicing assays, and gel electrophoresis of patient muscle proteins were performed to determine the molecular consequences of mutations of interest. RESULTS: Autosomal recessive compound heterozygous truncating mutations of the titin gene, TTN, were identified in 5 individuals. Biochemical analyses demonstrated increased titin degradation and truncated titin proteins in patient muscles, establishing the impact of the mutations. CONCLUSIONS: Our study identifies truncating TTN mutations as a cause of congenital myopathy that is reported as CNM. Unlike the classic CNM genes that are all involved in excitation-contraction coupling at the triad, TTN encodes the giant sarcomeric protein titin, which forms a myofibrillar backbone for the components of the contractile machinery. This study expands the phenotypic spectrum associated with TTN mutations and indicates that TTN mutation analysis should be considered in cases of possible CNM without mutations in the classic CNM genes.

Genome-wide association study of multiple congenital heart disease phenotypes identifies a susceptibility locus for atrial septal defect at chromosome 4p16.

We carried out a genome-wide association study (GWAS) of congenital heart disease (CHD). Our discovery cohort comprised 1,995 CHD cases and 5,159 controls and included affected individuals from each of the 3 major clinical CHD categories (with septal, obstructive and cyanotic defects). When all CHD phenotypes were considered together, no region achieved genome-wide significant association. However, a region on chromosome 4p16, adjacent to the MSX1 and STX18 genes, was associated (P = 9.5 × 10⁻7) with the risk of ostium secundum atrial septal defect (ASD) in the discovery cohort (N = 340 cases), and this association was replicated in a further 417 ASD cases and 2,520 controls (replication P = 5.0 × 10⁻5; odds ratio (OR) in replication cohort = 1.40, 95% confidence interval (CI) = 1.19-1.65; combined P = 2.6 × 10⁻¹°). Genotype accounted for ~9% of the population-attributable risk of ASD.

Genome-wide association study identifies loci on 12q24 and 13q32 associated with tetralogy of Fallot.

We conducted a genome-wide association study to search for risk alleles associated with Tetralogy of Fallot (TOF), using a northern European discovery set of 835 cases and 5159 controls. A region on chromosome 12q24 was associated (P = 1.4 × 10(-7)) and replicated convincingly (P = 3.9 × 10(-5)) in 798 cases and 2931 controls [per allele odds ratio (OR) = 1.27 in replication cohort, P = 7.7 × 10(-11) in combined populations]. Single nucleotide polymorphisms in the glypican 5 gene on chromosome 13q32 were also associated (P = 1.7 × 10(-7)) and replicated convincingly (P = 1.2 × 10(-5)) in 789 cases and 2927 controls (per allele OR = 1.31 in replication cohort, P = 3.03 × 10(-11) in combined populations). Four additional regions on chromosomes 10, 15 and 16 showed suggestive association accompanied by nominal replication. This study, the first genome-wide association study of a congenital heart malformation phenotype, provides evidence that common genetic variation influences the risk of TOF.

Contribution of global rare copy-number variants to the risk of sporadic congenital heart disease.

Previous studies have shown that copy-number variants (CNVs) contribute to the risk of complex developmental phenotypes. However, the contribution of global CNV burden to the risk of sporadic congenital heart disease (CHD) remains incompletely defined. We generated genome-wide CNV data by using Illumina 660W-Quad SNP arrays in 2,256 individuals with CHD, 283 trio CHD-affected families, and 1,538 controls. We found association of rare genic deletions with CHD risk (odds ratio [OR] = 1.8, p = 0.0008). Rare deletions in study participants with CHD had higher gene content (p = 0.001) with higher haploinsufficiency scores (p = 0.03) than they did in controls, and they were enriched with Wnt-signaling genes (p = 1 × 10(-5)). Recurrent 15q11.2 deletions were associated with CHD risk (OR = 8.2, p = 0.02). Rare de novo CNVs were observed in ~5% of CHD trios; 10 out of 11 occurred on the paternally transmitted chromosome (p = 0.01). Some of the rare de novo CNVs spanned genes known to be involved in heart development (e.g., HAND2 and GJA5). Rare genic deletions contribute ~4% of the population-attributable risk of sporadic CHD. Second to previously described CNVs at 1q21.1, deletions at 15q11.2 and those implicating Wnt signaling are the most significant contributors to the risk of sporadic CHD. Rare de novo CNVs identified in CHD trios exhibit paternal origin bias.

Phenotype-specific effect of chromosome 1q21.1 rearrangements and GJA5 duplications in 2436 congenital heart disease patients and 6760 controls.

Recurrent rearrangements of chromosome 1q21.1 that occur via non-allelic homologous recombination have been associated with variable phenotypes exhibiting incomplete penetrance, including congenital heart disease (CHD). However, the gene or genes within the ~1 Mb critical region responsible for each of the associated phenotypes remains unknown. We examined the 1q21.1 locus in 948 patients with tetralogy of Fallot (TOF), 1488 patients with other forms of CHD and 6760 ethnically matched controls using single nucleotide polymorphism genotyping arrays (Illumina 660W and Affymetrix 6.0) and multiplex ligation-dependent probe amplification. We found that duplication of 1q21.1 was more common in cases of TOF than in controls [odds ratio (OR) 30.9, 95% confidence interval (CI) 8.9-107.6); P = 2.2 × 10(-7)], but deletion was not. In contrast, deletion of 1q21.1 was more common in cases of non-TOF CHD than in controls [OR 5.5 (95% CI 1.4-22.0); P = 0.04] while duplication was not. We also detected rare (n = 3) 100-200 kb duplications within the critical region of 1q21.1 in cases of TOF. These small duplications encompassed a single gene in common, GJA5, and were enriched in cases of TOF in comparison to controls [OR = 10.7 (95% CI 1.8-64.3), P = 0.01]. These findings show that duplication and deletion at chromosome 1q21.1 exhibit a degree of phenotypic specificity in CHD, and implicate GJA5 as the gene responsible for the CHD phenotypes observed with copy number imbalances at this locus.

Programmed death ligand 1 (PD-L1) gene variants contribute to autoimmune Addison's disease and Graves' disease susceptibility.

CONTEXT: Despite much investigation, a substantial amount of the genetic susceptibility to autoimmune diseases remains unaccounted for. Recently, a single-nucleotide polymorphism (SNP) in the programmed death ligand 1 (PD-L1) gene has been associated with Graves' disease (GD) in a Japanese patient cohort. Our aim was to determine whether variants in PD-L1 are also associated with autoimmune Addison's disease (AAD) and to replicate the previous association in patients with GD from the United Kingdom. DESIGN AND PATIENTS: We analyzed eight SNPs within PD-L1 in a United Kingdom cohort of 315 AAD subjects and 316 healthy controls. We then replicated our experiment in a cohort of 342 Norwegian AAD cases and 379 controls and in 496 United Kingdom GD subjects. RESULTS: Three of the eight SNPs studied, part of a haplotype block in the PD-L1 gene, showed modest association with both AAD and GD in the United Kingdom cohort, with maximum evidence at the marker RS1411262 [United Kingdom AAD odds ratio 1.33 (5-95% confidence interval 1.02-1.73), P(genotype) = 0.028; GD odds ratio 1.36 (5-95% confidence interval 1.07-1.72), P(genotype) = 0.033]. Association with genotypes at the same three markers was confirmed in the Norwegian AAD cohort [P(genotype) = 0.011-0.020]. A recessive effect at the most associated alleles was observed in both the AAD and GD cohorts. CONCLUSIONS: We confirm the role of PD-L1 variants in GD susceptibility and extend these findings to demonstrate association in two Northern European patient cohorts with AAD. PD-L1 joins the growing number of known susceptibility loci exerting modest effects in these autoimmune disorders.