TITINdb-a computational tool to assess titin's role as a disease gene

TITINdb2-expanding annotation and structural information for protein variants in the giant sarcomeric protein titin

Summary

We present TITINdb2, an update to the TITINdb database previously constructed to facilitate the identification of pathogenic missense variants in the giant protein titin, which are associated with a variety of skeletal and cardiac myopathies. The database and web portal have been substantially revised and include the following new features: (i) an increase in computational annotation from 4 to 20 variant impact predictors, available through a new custom data table dialogue; (ii) through structural coverage of single domains with AlphaFold2 predicted models; (iii) newly predicted domain-domain interface annotations; (iv) an expanded in silico saturation mutagenesis incorporating four variant impact predictors; (v) a comprehensive overhaul of available data, including population data sources and variants reported pathogenic in the literature; and (vi) a curated mapping of existing protein, transcript, and chromosomal sequence positions and a new variant conversion tool to translate variants in one format to any other format.

Availability and implementation

The database is accessible via titindb.kcl.ac.uk/TITINdb/.

High-coverage whole-genome sequencing of a Jakun individual from the "Orang Asli" Proto-Malay subtribe from Peninsular Malaysia

Jakun, a Proto-Malay subtribe from Peninsular Malaysia, is believed to have inhabited the Malay Archipelago during the period of agricultural expansion approximately 4 thousand years ago (kya). However, their genetic structure and population history remain inconclusive. In this study, we report the genome structure of a Jakun female, based on whole-genome sequencing, which yielded an average coverage of 35.97-fold. We identified approximately 3.6 million single-nucleotide variations (SNVs) and 517,784 small insertions/deletions (indels). Of these, 39,916 SNVs were novel (referencing dbSNP151), and 10,167 were nonsynonymous (nsSNVs), spanning 5674 genes. Principal Component Analysis (PCA) revealed that the Jakun genome sequence closely clustered with the genomes of the Cambodians (CAM) and the Metropolitan Malays from Singapore (SG_MAS). The ADMIXTURE analysis further revealed potential admixture from the EA and North Borneo populations, as corroborated by the results from the F3, F4, and TreeMix analyses. Mitochondrial DNA analysis revealed that the Jakun genome carried the N21a haplogroup (estimated to have occurred ~19 kya), which is commonly found among Malays from Malaysia and Indonesia. From the whole-genome sequence data, we identified 825 damaging and deleterious nonsynonymous single-nucleotide polymorphisms (nsSNVs) affecting 720 genes. Some of these variants are associated with age-related macular degeneration, atrial fibrillation, and HDL cholesterol level. Additionally, we located a total of 3310 variants on 32 core adsorption, distribution, metabolism, and elimination (ADME) genes. Of these, 193 variants are listed in PharmGKB, and 21 are nsSNVs. In summary, the genetic structure identified in the Jakun individual could enhance the mapping of genetic variants for disease-based population studies and further our understanding of the human migration history in Southeast Asia.

Variant Impact Predictor database (VIPdb), version 2: trends from three decades of genetic variant impact predictors

BACKGROUND

Variant interpretation is essential for identifying patients' disease-causing genetic variants amongst the millions detected in their genomes. Hundreds of Variant Impact Predictors (VIPs), also known as Variant Effect Predictors (VEPs), have been developed for this purpose, with a variety of methodologies and goals. To facilitate the exploration of available VIP options, we have created the Variant Impact Predictor database (VIPdb).

RESULTS

The Variant Impact Predictor database (VIPdb) version 2 presents a collection of VIPs developed over the past three decades, summarizing their characteristics, ClinGen calibrated scores, CAGI assessment results, publication details, access information, and citation patterns. We previously summarized 217 VIPs and their features in VIPdb in 2019. Building upon this foundation, we identified and categorized an additional 190 VIPs, resulting in a total of 407 VIPs in VIPdb version 2. The majority of the VIPs have the capacity to predict the impacts of single nucleotide variants and nonsynonymous variants. More VIPs tailored to predict the impacts of insertions and deletions have been developed since the 2010s. In contrast, relatively few VIPs are dedicated to the prediction of splicing, structural, synonymous, and regulatory variants. The increasing rate of citations to VIPs reflects the ongoing growth in their use, and the evolving trends in citations reveal development in the field and individual methods.

CONCLUSIONS

VIPdb version 2 summarizes 407 VIPs and their features, potentially facilitating VIP exploration for various variant interpretation applications. VIPdb is available at  https://genomeinterpretation.org/vipdb.

Variant Impact Predictor database (VIPdb), version 2: Trends from 25 years of genetic variant impact predictors

Background

Variant interpretation is essential for identifying patients' disease-causing genetic variants amongst the millions detected in their genomes. Hundreds of Variant Impact Predictors (VIPs), also known as Variant Effect Predictors (VEPs), have been developed for this purpose, with a variety of methodologies and goals. To facilitate the exploration of available VIP options, we have created the Variant Impact Predictor database (VIPdb).

Results

The Variant Impact Predictor database (VIPdb) version 2 presents a collection of VIPs developed over the past 25 years, summarizing their characteristics, ClinGen calibrated scores, CAGI assessment results, publication details, access information, and citation patterns. We previously summarized 217 VIPs and their features in VIPdb in 2019. Building upon this foundation, we identified and categorized an additional 186 VIPs, resulting in a total of 403 VIPs in VIPdb version 2. The majority of the VIPs have the capacity to predict the impacts of single nucleotide variants and nonsynonymous variants. More VIPs tailored to predict the impacts of insertions and deletions have been developed since the 2010s. In contrast, relatively few VIPs are dedicated to the prediction of splicing, structural, synonymous, and regulatory variants. The increasing rate of citations to VIPs reflects the ongoing growth in their use, and the evolving trends in citations reveal development in the field and individual methods.

Conclusions

VIPdb version 2 summarizes 403 VIPs and their features, potentially facilitating VIP exploration for various variant interpretation applications.

Availability

VIPdb version 2 is available at https://genomeinterpretation.org/vipdb.

A latent variable model for evaluating mutual exclusivity and co-occurrence between driver mutations in cancer

A key challenge in cancer genomics is understanding the functional relationships and dependencies between combinations of somatic mutations that drive cancer development. Such driver mutations frequently exhibit patterns of mutual exclusivity or co-occurrence across tumors, and many methods have been developed to identify such dependency patterns from bulk DNA sequencing data of a cohort of patients. However, while mutual exclusivity and co-occurrence are described as properties of driver mutations, existing methods do not explicitly disentangle functional, driver mutations from neutral, passenger mutations. In particular, nearly all existing methods evaluate mutual exclusivity or co-occurrence at the gene level, marking a gene as mutated if any mutation - driver or passenger - is present. Since some genes have a large number of passenger mutations, existing methods either restrict their analyses to a small subset of suspected driver genes - limiting their ability to identify novel dependencies - or make spurious inferences of mutual exclusivity and co-occurrence involving genes with many passenger mutations. We introduce DIALECT, an algorithm to identify dependencies between pairs of driver mutations from somatic mutation counts. We derive a latent variable mixture model for drivers and passengers that combines existing probabilistic models of passenger mutation rates with a latent variable describing the unknown status of a mutation as a driver or passenger. We use an expectation maximization (EM) algorithm to estimate the parameters of our model, including the rates of mutually exclusivity and co-occurrence between drivers. We demonstrate that DIALECT more accurately infers mutual exclusivity and co-occurrence between driver mutations compared to existing methods on both simulated mutation data and somatic mutation data from 5 cancer types in The Cancer Genome Atlas (TCGA).

Walking with giants: The challenges of variant impact assessment in the giant sarcomeric protein titin

Titin, the so-called "third filament" of the sarcomere, represents a difficult challenge for the determination of damaging genetic variants. A single titin molecule extends across half the length of a sarcomere in striated muscle, fulfilling a variety of vital structural and signaling roles, and has been linked to an equally varied range of myopathies, resulting in a significant burden on individuals and healthcare systems alike. While the consequences of truncating variants of titin are well-documented, the ramifications of the missense variants prevalent in the general population are less so. We here present a compendium of titin missense variants-those that result in a single amino-acid substitution in coding regions-reported to be pathogenic and discuss these in light of the nature of titin and the variant position within the sarcomere and their domain, the structural, pathological, and biophysical characteristics that define them, and the methods used for characterization. Finally, we discuss the current knowledge and integration of the multiple fields that have contributed to our understanding of titin-related pathology and offer suggestions as to how these concurrent methodologies may aid the further development in our understanding of titin and hopefully extend to other, less well-studied giant proteins. This article is categorized under: Cardiovascular Diseases > Genetics/Genomics/Epigenetics Congenital Diseases > Genetics/Genomics/Epigenetics Congenital Diseases > Molecular and Cellular Physiology.

Missense mutations in the central domains of cardiac myosin binding protein-C and their potential contribution to hypertrophic cardiomyopathy

Myosin binding protein-C (MyBP-C) is a multidomain protein that regulates muscle contraction. Mutations in MYBPC3, the gene encoding for the cardiac variant (henceforth called cMyBP-C), are amongst the most frequent causes of hypertrophic cardiomyopathy. Most mutations lead to a truncated version of cMyBP-C, which is most likely unstable. However, missense mutations have also been reported, which tend to cluster in the central domains of the cMyBP-C molecule. This suggests that these central domains are more than just a passive spacer between the better characterized N- and C-terminal domains. Here, we investigated the potential impact of four different missense mutations, E542Q, G596R, N755K, and R820Q, which are spread over the domains C3 to C6, on the function of MyBP-C on both the isolated protein level and in cardiomyocytes in vitro. Effect on domain stability, interaction with thin filaments, binding to myosin, and subcellular localization behavior were assessed. Our studies show that these missense mutations result in slightly different phenotypes at the molecular level, which are mutation specific. The expected functional readout of each mutation provides a valid explanation for why cMyBP-C fails to work as a brake in the regulation of muscle contraction, which eventually results in a hypertrophic cardiomyopathy phenotype. We conclude that missense mutations in cMyBP-C must be evaluated in context of their domain localization, their effect on interaction with thin filaments and myosin, and their effect on protein stability to explain how they lead to disease.

Aberrant mRNA processing caused by splicing mutations in TTN-related neuromuscular disorders

Mutations in the TTN gene have been reported to be responsible for a range of neuromuscular disorders, including recessive distal myopathy and congenital myopathy (CM). Only five splicing mutations have been identified to induce aberrant mRNA splicing in TTN-related neuromuscular disorders. In our study, we described detailed clinical characteristics, muscle pathology and genetic analysis of two probands with TTN-related autosomal recessive neuromuscular disorders. Besides, we identified two novel intronic mutations, c.107377+1 G > C in intron 362 and c.19994-2 A > G in intron 68, in the two probands. Through cDNA analysis, we revealed the c.107377+1 G > C mutation induced retention of the entire intron 362, and the c.19994-2 A > G mutation triggered skipping of the first 11 bp of exon 69. Our study broadens the aberrant splicing spectrum of neuromuscular disorders caused by splicing mutations in the TTN gene.

Successful heart transplant in a child with congenital core myopathy and delayed-onset restrictive cardiomyopathy due to recessive mutations in the titin (TTN) gene

BACKGROUND

Mutations in the TTN gene, encoding the muscle filament titin, are a major cause of inherited dilated cardiomyopathy. Early-onset skeletal muscle disorders due to recessive TTN mutations have recently been described, sometimes associated with cardiomyopathies.

CASE DESCRIPTION

We report the case of a boy with congenital core myopathy due to compound heterozygosity for TTN variants. He presented in infancy with rapidly evolving restrictive cardiomyopathy, requiring heart transplantation at the age of 5 years with favorable long-term cardiac and neuromuscular outcome.

CONCLUSION

Heart transplantation may have a role in selected patients with TTN-related congenital myopathy with disproportionally severe cardiac presentation compared to skeletal and respiratory muscle involvement.

Structure determination and analysis of titin A-band fibronectin type III domains provides insights for disease-linked variants and protein oligomerisation

Titin is the largest protein found in nature and spans half a sarcomere in vertebrate striated muscle. The protein has multiple functions, including in the organisation of the thick filament and acting as a molecular spring during the muscle contraction cycle. Missense variants in titin have been linked to both cardiac and skeletal myopathies. Titin is primarily composed of tandem repeats of immunoglobulin and fibronectin type III (Fn3) domains in a variety of repeat patterns; however, the vast majority of these domains have not had their high-resolution structure determined experimentally. Here, we present the crystal structures of seven wild type titin Fn3 domains and two harbouring rare missense variants reported in hypertrophic cardiomyopathy (HCM) patients. All domains present the typical Fn3 fold, with the domains harbouring variants reported in HCM patients retaining the wild-type conformation. The effect on domain folding and stability were assessed for five rare missense variants found in HCM patients: four caused thermal destabilization of between 7 and 13 °C and one prevented the folding of its domain. The structures also allowed us to locate the positions of residues whose mutations have been linked to congenital myopathies and rationalise how they convey their deleterious effects. We find no evidence of physiological homodimer formation, excluding one hypothesised mechanism as to how titin variants could exert pathological effects.

Titin related myopathy with ophthalmoplegia. A novel phenotype

Titin-related myopathy is an emerging genetic neuromuscular disorder with a wide spectrum of clinical phenotypes. To date, there have not been reports of patients with this disease that presented with extraocular muscle involvement. Here we discuss a 19-year-old male with congenital weakness, complete ophthalmoplegia, thoracolumbar scoliosis, and obstructive sleep apnea. Muscle magnetic resonance imaging revealed severe involvement of the gluteal and anterior compartment muscles, and clear adductor sparing, while muscle biopsy of the right vastus lateralis showed distinctive cap-like structures. Trio Whole Exome Sequencing (WES) showed compound heterozygous likely pathologic variants in the TTN gene. (c.82541_82544dup (p.Arg27515Serfs*2) in exon 327 (NM_001267550.2) and c.31846+1G>A (p.?) in exon 123 (NM_001267550.2). To our knowledge, this is the first report of a TTN-related disorder associated with ophthalmoplegia.

Structural and signaling proteins in the Z-disk and their role in cardiomyopathies

The sarcomere is the smallest functional unit of muscle contraction. It is delineated by a protein-rich structure known as the Z-disk, alternating with M-bands. The Z-disk anchors the actin-rich thin filaments and plays a crucial role in maintaining the mechanical stability of the cardiac muscle. A multitude of proteins interact with each other at the Z-disk and they regulate the mechanical properties of the thin filaments. Over the past 2 decades, the role of the Z-disk in cardiac muscle contraction has been assessed widely, however, the impact of genetic variants in Z-disk proteins has still not been fully elucidated. This review discusses the various Z-disk proteins (alpha-actinin, filamin C, titin, muscle LIM protein, telethonin, myopalladin, nebulette, and nexilin) and Z-disk-associated proteins (desmin, and obscurin) and their role in cardiac structural stability and intracellular signaling. This review further explores how genetic variants of Z-disk proteins are linked to inherited cardiac conditions termed cardiomyopathies.

How Functional Genomics Can Keep Pace With VUS Identification

Over the last two decades, an exponentially expanding number of genetic variants have been identified associated with inherited cardiac conditions. These tremendous gains also present challenges in deciphering the clinical relevance of unclassified variants or variants of uncertain significance (VUS). This review provides an overview of the advancements (and challenges) in functional and computational approaches to characterize variants and help keep pace with VUS identification related to inherited heart diseases.

Pairwise sequence similarity mapping with PaSiMap: Reclassification of immunoglobulin domains from titin as case study

•

A novel multidimensional scaling pipeline for sequence analysis.

•

A simple way to distinguish between unique and shared sequence features.

•

Titin domains were reclassified, improving upon earlier analysis.

Sequence comparison is critical for the functional assignment of newly identified protein genes. As uncharacterized protein sequences accumulate, there is an increasing need for sensitive tools for their classification. Here, we present a novel multidimensional scaling pipeline, PaSiMap, which creates a map of pairwise sequence similarities. Uniquely, PaSiMap distinguishes between unique and shared features, allowing for a distinct view of protein-sequence relationships. We demonstrate PaSiMap’s efficiency in detecting sequence groups and outliers using titin’s 169 immunoglobulin (Ig) domains. We show that Ig domain similarity is hierarchical, being firstly determined by chain location, then by the loop features of the Ig fold and, finally, by super-repeat position. The existence of a previously unidentified domain repeat in the distal, constitutive I-band is revealed. Prototypic Igs, plus notable outliers, are identified and thereby domain classification improved. This re-classification can now guide future molecular research. In summary, we demonstrate that PaSiMap is a sensitive tool for the classification of protein sequences, which adds a new perspective in the understanding of inter-protein relationships. PaSiMap is applicable to any biological system defined by a linear sequence, including polynucleotide chains.

In Silico Analysis of Novel Titin Non-Synonymous Missense Variants Detected by Targeted Next-Generation Sequencing in a Cohort of Romanian Index Patients with Hypertrophic Cardiomyopathy

Background and aim

Most of detected variants in cardiogenetic panels are still classified as variants of unknown significance, requiring supplementary analyses for a definite classification. Performing further in-depth studies on such vast number of candidates is unfeasible. We sought to prioritise the novel nonsynonymous missense variants identified in titin gene (TTN) in a cohort of Romanian index cases with hypertrophic cardiomyopathy (HCM).

Methods

45 unrelated probands with HCM were screened by targeted next generation sequencing (NGS) covering all TTN exons. A stepwise strategy was used to select and prioritize the candidate variants for subsequent investigation.

Results

Using rigorous bioinformatic filtering, 7 novel TTN nonsynonymous missense variants were identified and were the subject of in silico sequential analysis. 4 of the 7 variants were predicted to be possibly pathogenic by the Mendelian Clinically Applicable Pathogenicity (M-CAP) algorithm. Of these, three sequence variants (c.30392G>T, c.2518G>T, and c.49G>T) were also predicted to be destabilizing according to the second computational tool (TITINdb) and were designated as likely function-impacting.

Conclusions

Herein we presented our strategy to hand-pick the novel TTN missense variants to be considered for further experimental studies. By applying various in silico tools, we restricted the list of sequence variants to be investigated to those most likely to be disease-associated, and thus reducing the need to perform expensive and time-consuming additional studies.

Solution NMR Structure of Titin N2A Region Ig Domain I83 and Its Interaction with Metal Ions

Titin, the largest single chain protein known so far, has long been known to play a critical role in passive muscle function but recent studies have highlighted titin's role in active muscle function. One of the key elements in this role is the Ca²⁺-dependent interaction between titin's N2A region and the thin filament. An important element in this interaction is I83, the terminal immunoglobulin domain in the N2A region. There is limited structural information about this domain, but experimental evidence suggests that it plays a critical role in the N2A-actin binding interaction. We now report the solution NMR structure of I83 and characterize its dynamics and metal binding properties in detail. Its structure shows interesting relationships to other I-band Ig domains. Metal binding and dynamics data point towards the way the domain is evolutionarily optimized to interact with neighbouring domains. We also identify a calcium binding site on the N-terminal side of I83, which is expected to impact the interdomain interaction with the I82 domain. Together these results provide a first step towards a better understanding of the physiological effects associated with deletion of most of the I83 domain, as occurs in the mdm mouse model, as well as for future investigations of the N2A region.

Is Gene-Size an Issue for the Diagnosis of Skeletal Muscle Disorders?

Human genes have a variable length. Those having a coding sequence of extraordinary length and a high number of exons were almost impossible to sequence using the traditional Sanger-based gene-by-gene approach. High-throughput sequencing has partly overcome the size-related technical issues, enabling a straightforward, rapid and relatively inexpensive analysis of large genes.

Several large genes (e.g. TTN, NEB, RYR1, DMD) are recognized as disease-causing in patients with skeletal muscle diseases. However, because of their sheer size, the clinical interpretation of variants in these genes is probably the most challenging aspect of the high-throughput genetic investigation in the field of skeletal muscle diseases.

The main aim of this review is to discuss the technical and interpretative issues related to the diagnostic investigation of large genes and to reflect upon the current state of the art and the future advancements in the field.

Making sense of missense variants in TTN-related congenital myopathies

Mutations in the sarcomeric protein titin, encoded by TTN, are emerging as a common cause of myopathies. The diagnosis of a TTN-related myopathy is, however, often not straightforward due to clinico-pathological overlap with other myopathies and the prevalence of TTN variants in control populations. Here, we present a combined clinico-pathological, genetic and biophysical approach to the diagnosis of TTN-related myopathies and the pathogenicity ascertainment of TTN missense variants. We identified 30 patients with a primary TTN-related congenital myopathy (CM) and two truncating variants, or one truncating and one missense TTN variant, or homozygous for one TTN missense variant. We found that TTN-related myopathies show considerable overlap with other myopathies but are strongly suggested by a combination of certain clinico-pathological features. Presentation was typically at birth with the clinical course characterized by variable progression of weakness, contractures, scoliosis and respiratory symptoms but sparing of extraocular muscles. Cardiac involvement depended on the variant position. Our biophysical analyses demonstrated that missense mutations associated with CMs are strongly destabilizing and exert their effect when expressed on a truncating background or in homozygosity. We hypothesise that destabilizing TTN missense mutations phenocopy truncating variants and are a key pathogenic feature of recessive titinopathies that might be amenable to therapeutic intervention.

Improved Criteria for the Classification of Titin Variants in Inherited Skeletal Myopathies

BACKGROUND

Extensive genetic screening results in the identification of thousands of rare variants that are difficult to interpret. Because of its sheer size, rare variants in the titin gene (TTN) are detected frequently in any individual. Unambiguous interpretation of molecular findings is almost impossible in many patients with myopathies or cardiomyopathies.

OBJECTIVE

To refine the current classification framework for TTN-associated skeletal muscle disorders and standardize the interpretation of TTN variants.

METHODS

We used the guidelines issued by the American College of Medical Genetics and Genomics (ACMG) and the Association for Molecular Pathology (AMP) to re-analyze TTN genetic findings from our patient cohort.

RESULTS

We identified in the classification guidelines three rules that are not applicable to titin-related skeletal muscle disorders; six rules that require disease-/gene-specific adjustments and four rules requiring quantitative thresholds for a proper use. In three cases, the rule strength need to be modified.

CONCLUSIONS

We suggest adjustments are made to the guidelines. We provide frequency thresholds to facilitate filtering of candidate causative variants and guidance for the use and interpretation of functional data and co-segregation evidence. We expect that the variant classification framework for TTN-related skeletal muscle disorders will be further improved along with a better understanding of these diseases.

The importance of an integrated genotype-phenotype strategy to unravel the molecular bases of titinopathies

Next generation sequencing (NGS) has allowed the titin gene (TTN) to be identified as a major contributor to neuromuscular disorders, with high clinical heterogeneity. The mechanisms underlying the phenotypic variability and the dominant or recessive pattern of inheritance are unclear. Titin is involved in the formation and stability of the sarcomeres. The effects of the different TTN variants can be harmless or pathogenic (recessive or dominant) but the interpretation is tricky because the current bioinformatics tools can not predict their functional impact effectively. Moreover, TTN variants are very frequent in the general population. The combination of deep phenotyping associated with RNA molecular analyses, western blot (WB) and functional studies is often essential for the interpretation of genetic variants in patients suspected of titinopathy. In line with the current guidelines and suggestions, we implemented for patients with skeletal myopathy and with potentially disease causing TTN variant(s) an integrated genotype-transcripts-protein-phenotype approach, associated with phenotype and variants segregation studies in relatives and confrontation with published data on titinopathies to evaluate pathogenic effects of TTN variants (even truncating ones) on titin transcripts, amount, size and functionality. We illustrate this integrated approach in four patients with recessive congenital myopathy.

Genotype-phenotype correlations in recessive titinopathies

PURPOSE

High throughput sequencing analysis has facilitated the rapid analysis of the entire titin (TTN) coding sequence. This has resulted in the identification of a growing number of recessive titinopathy patients. The aim of this study was to (1) characterize the causative genetic variants and clinical features of the largest cohort of recessive titinopathy patients reported to date and (2) to evaluate genotype-phenotype correlations in this cohort.

METHODS

We analyzed clinical and genetic data in a cohort of patients with biallelic pathogenic or likely pathogenic TTN variants. The cohort included both previously reported cases (100 patients from 81 unrelated families) and unreported cases (23 patients from 20 unrelated families).

RESULTS

Overall, 132 causative variants were identified in cohort members. More than half of the cases had hypotonia at birth or muscle weakness and a delayed motor development within the first 12 months of life (congenital myopathy) with causative variants located along the entire gene. The remaining patients had a distal or proximal phenotype and a childhood or later (noncongenital) onset. All noncongenital cases had at least one pathogenic variant in one of the final three TTN exons (362-364).

CONCLUSION

Our findings suggest a novel association between the location of nonsense variants and the clinical severity of the disease.

The Axial Alignment of Titin on the Muscle Thick Filament Supports Its Role as a Molecular Ruler

The giant protein titin is expressed in vertebrate striated muscle where it spans half a sarcomere from the Z-disc to the M-band and is essential for muscle organisation, activity and health. The C-terminal portion of titin is closely associated with the thick, myosin-containing filament and exhibits a complex pattern of immunoglobulin and fibronectin domains. This pattern reflects features of the filament organisation suggesting that it acts as a molecular ruler and template, but the exact axial disposition of the molecule has not been determined. Here, we present data that allow us to precisely locate titin domains axially along the thick filament from its tip to the edge of the bare zone. We find that the domains are regularly distributed along the filament at 4-nm intervals and we can determine the domains that associate with features of the filament, such as the 11 stripes of accessory proteins. We confirm that the nine stripes ascribed to myosin binding protein-C are not related to the titin sequence previously assumed; rather, they relate to positions approximately 18 domains further towards the C terminus along titin. This disposition also allows a subgroup of titin domains comprising two or three fibronectin domains to associate with each of the 49 levels of myosin heads in each half filament. The results strongly support the role of titin as a blueprint for the thick filament and the arrangement of the myosin motor domains.

A long-read RNA-seq approach to identify novel transcripts of very large genes

RNA-seq is widely used for studying gene expression, but commonly used sequencing platforms produce short reads that only span up to two exon junctions per read. This makes it difficult to accurately determine the composition and phasing of exons within transcripts. Although long-read sequencing improves this issue, it is not amenable to precise quantitation, which limits its utility for differential expression studies. We used long-read isoform sequencing combined with a novel analysis approach to compare alternative splicing of large, repetitive structural genes in muscles. Analysis of muscle structural genes that produce medium (Nrap: 5 kb), large (Neb: 22 kb), and very large (Ttn: 106 kb) transcripts in cardiac muscle, and fast and slow skeletal muscles identified unannotated exons for each of these ubiquitous muscle genes. This also identified differential exon usage and phasing for these genes between the different muscle types. By mapping the in-phase transcript structures to known annotations, we also identified and quantified previously unannotated transcripts. Results were confirmed by endpoint PCR and Sanger sequencing, which revealed muscle-type-specific differential expression of these novel transcripts. The improved transcript identification and quantification shown by our approach removes previous impediments to studies aimed at quantitative differential expression of ultralong transcripts.

The importance of chain context in assessing small nucleotide variants in titin: in silico case study of the I10-I11 tandem and its arrhythmic right ventricular cardiomyopathy linked position T2580

Non-synonymous small nucleotide variations (nsSNVs) in the giant muscle protein, titin, have key roles in the development of several myopathologies. Although there is considerable motive to screen at-risk individuals for nsSNVs, to identify patients in early disease stages while therapeutic intervention is still possible, the clinical significance of most titin variations remains unclear. Therefore, there is a growing need to establish methods to classify nsSNVs in a simple, economic and rapid manner. Due to its strong correlation to arrhythmogenic right ventricular cardiomyopathy (ARVC), one particular mutation in titin-T2580I, located in the I10 immunoglobulin domain-has received considerable attention. Here, we use the I10-I11 tandem as a case study to explore the possible benefits of considering the titin chain context-i.e. domain interfaces-in the assessment of titin nsSNVs. Specifically, we investigate which exchanges mimic the conformational molecular phenotype of the T2580I mutation at the I10-I11 domain interface. Then, we computed a residue stability landscape for domains alone and in tandem to define a Domain Interface Score (DIS) which identifies several hotspot residues. Our findings suggest that the T2580 position is highly sensitive to exchange and that any variant found in this position should be considered with care. Furthermore, we conclude that the consideration of the higher order structure of the titin chain is important to gain accurate insights into the vulnerability of positions in linker regions and that titin nsSNV prediction benefits from a contextual analysis. Communicated by Ramaswamy H. Sarma.

A new congenital multicore titinopathy associated with fast myosin heavy chain deficiency

Congenital titinopathies are myopathies with variable phenotypes and inheritance modes. Here, we fully characterized, using an integrated approach (deep phenotyping, muscle morphology, mRNA and protein evaluation in muscle biopsies), two siblings with congenital multicore myopathy harboring three TTN variants predicted to affect titin stability and titin-myosin interactions. Muscle biopsies showed multicores, type 1 fiber uniformity and sarcomeric structure disruption with some thick filament loss. Immunohistochemistry and Western blotting revealed a marked reduction of fast myosin heavy chain isoforms. This is the first observation of a titinopathy suggesting that titin defect leads to secondary loss of fast myosin heavy chain isoforms.

Structural diversity in the atomic resolution 3D fingerprint of the titin M-band segment

In striated muscles, molecular filaments are largely composed of long protein chains with extensive arrays of identically folded domains, referred to as “beads-on-a-string”. It remains a largely unresolved question how these domains have developed a unique molecular profile such that each carries out a distinct function without false-positive readout. This study focuses on the M-band segment of the sarcomeric protein titin, which comprises ten identically folded immunoglobulin domains. Comparative analysis of high-resolution structures of six of these domains ‒ M1, M3, M4, M5, M7, and M10 ‒ reveals considerable structural diversity within three distinct loops and a non-conserved pattern of exposed cysteines. Our data allow to structurally interpreting distinct pathological readouts that result from titinopathy-associated variants. Our findings support general principles that could be used to identify individual structural/functional profiles of hundreds of identically folded protein domains within the sarcomere and other densely crowded cellular environments.

Novel TTN mutations and muscle imaging characteristics in congenital titinopathy

Objective

We present clinical features, muscle imaging findings, and genetic characteristics of five unrelated Chinese patients with congenital titinopathy, emphasizing the diagnostic role of muscle MRI.

Methods

Five patients who recessive titinopathies were recruited. All patients received muscle biopsies. Mutations were detected by panel massively parallel sequencing and confirmed by Sanger sequencing. Western blotting of muscle proteins was performed. Leg muscle MRIs were performed in four patients.

Results

Four patients aged 1–4 years old showed delayed motor development from early infancy, while a 17‐year‐old boy showed only a 1‐year history of exercise intolerance. Physical examination showed proximal weakness in three patients. Muscle biopsies demonstrated multiple myopathological changes, including increased internalized nuclei, multicores, central cores, and dystrophic changes. Genetic sequencing revealed compound heterozygous or homozygous novel TTN mutations, including six frameshift mutations, one nonsense mutation, two missense mutations, one splicing mutation, and one small nonframeshift deletion. Protein analyses revealed significant decrease of full‐length titin in all patients. Thigh muscle MRIs in four patients showed prominent fatty infiltration in the upper portion of semitendinosus and the peripheral portion of gluteus medius, while the sartorius and gracilis were relatively preserved.

Interpretation

These cases provided further evidence that TTN mutations are likely responsible for an increasing proportion of congenital myopathies than currently recognized. The novel mutations reported expand the mutation spectrum of the TTN gene. There is a characteristic pattern of muscle involvement in congenital titinopathy regardless of clinical or pathological phenotype, providing valuable clues for guiding a genetic diagnosis workup.

The giant titin: how to evaluate its role in cardiomyopathies

Titin, the largest protein known, has attracted a lot of interest in the cardiovascular field in recent years, since the discovery that truncating variants in titin are commonly found in patients with dilated cardiomyopathy. This review will discuss the contribution of variants in titin to inherited cardiac conditions (cardiomyopathies) and how model systems, such as animals and cellular systems, can help to provide insights into underlying disease mechanisms. It will also give an outlook onto exciting technological developments, such as in the field of CRISPR, which may facilitate future research on titin variants and their contributions to cardiomyopathies.

Fine mapping titin's C-zone: Matching cardiac myosin-binding protein C stripes with titin's super-repeats

Titin is largely comprised of serially-linked immunoglobulin (Ig) and fibronectin type-III (Fn3) domains. Many of these domains are arranged in an 11 domain super-repeat pattern that is repeated 11 times, forming the so-named titin C-zone in the A-band region of the sarcomere. Each super-repeat is thought to provide binding sites for thick filament proteins, such as cMyBP-C (cardiac myosin-binding protein C). However, it remains to be established which of titin's 11 C-zone super-repeats anchor cMyBP-C as titin contains 11 super-repeats and cMyBP-C is found in 9 stripes only. To study the layout of titin's C-zone in relation to MyBP-C, immunolabeling studies were performed on mouse skinned myocardium with antibodies to titin and cMyBP-C, using both immuno-electron microscopy and super-resolution optical microscopy. Results indicate that cMyBP-C locates near the interface between titin's C-zone super-repeats. Studies on a mouse model in which two of titin's C-zone repeats have been genetically deleted support that the first Ig domain of a super-repeat is important for anchoring cMyBP-C but also Fn3 domains located at the end of the preceding repeat. Furthermore, not all super-repeat interfaces are equal as the interface between super-repeat 1 and 2 (close to titin's D-zone) does not contain cMyBP-C. Finally, titin's C-zone does not extend all the way to the bare zone but instead terminates at the level of the second myosin crown. This study enhances insights in the molecular layout of the C-zone of titin, its relation to cMyBP-C, and its possible roles in cardiomyopathies.

Titin mutations and muscle disease

The introduction of next-generation sequencing technology has revealed that mutations in the gene that encodes titin (TTN) are linked to multiple skeletal and cardiac myopathies. The most prominent of these myopathies is dilated cardiomyopathy (DCM). Over 60 genes are linked to the etiology of DCM, but by far, the leading cause of DCM is mutations in TTN with truncating variants in TTN (TTNtvs) associated with familial DCM in ∼ 20% of the cases. Titin is a large (3-4 MDa) and abundant protein that forms the third myofilament type of striated muscle where it spans half the sarcomere, from the Z-disk to the M-line. The underlying mechanisms by which titin mutations induce disease are poorly understood and targeted therapies are not available. Here, we review what is known about TTN mutations in muscle disease, with a major focus on DCM. We highlight that exon skipping might provide a possible therapeutic avenue to address diseases that arise from TTNtvs.

[Protein-protein interacting networks, their structures and disease-related mutations]

In recent years, the comparison of protein interactomes has identified conserved modules, that could represent functional nuclei with a common ancestry. Within this context, recent analyses of protein-protein interacting networks have led to a debate on the influence of the experimental method on the quality and biological pertinence of these data. It is crucial to understand the measure in which divergence between networks of different species reflect sampling biases in respective experimental methods, as opposed to topological features dictated by biological functionality. This aspect requires novel, precise and practical mathematical tools, to quantify and compare high resolution networks. To this end, we have studied the relationship between pools of aleatory graphs and real biological signalization networks, while stressing the number of graph cycles in the networks, which represent complexes in experimental protein interactomes. By combining methods for graph and algorithm dynamics to count the loops, we evaluate the relative importance of the loops in biological networks in comparison with network analyses.

When signalling goes wrong: pathogenic variants in structural and signalling proteins causing cardiomyopathies

Cardiomyopathies are a diverse group of cardiac disorders with distinct phenotypes, depending on the proteins and pathways affected. A substantial proportion of cardiomyopathies are inherited and those will be the focus of this review article. With the wide application of high-throughput sequencing in the practice of clinical genetics, the roles of novel genes in cardiomyopathies are recognised. Here, we focus on a subgroup of cardiomyopathy genes [TTN, FHL1, CSRP3, FLNC and PLN, coding for Titin, Four and a Half LIM domain 1, Muscle LIM Protein, Filamin C and Phospholamban, respectively], which, despite their diverse biological functions, all have important signalling functions in the heart, suggesting that disturbances in signalling networks can contribute to cardiomyopathies.