The French National Registry of patients with Facioscapulohumeral muscular dystrophy.

BACKGROUND: Facioscapulohumeral muscular dystrophy is a rare inherited neuromuscular disease with an estimated prevalence of 1/20,000 and France therefore harbors about 3000 FSHD patients. With research progress and the development of targeted therapies, patients' identification through registries can facilitate and improve recruitment in clinical trials and studies. RESULTS: The French National Registry of FSHD patients was designed as a mixed model registry involving both patients and physicians, through self-report and clinical evaluation questionnaires respectively, to collect molecular and clinical data. Because of the limited number of patients, data quality is a major goal of the registry and various automatic data control features have been implemented in the bioinformatics system. In parallel, data are manually validated by molecular and clinical curators. Since its creation in 2013, data from 638 FSHD patients have been collected, representing about 21% of the French FSHD population. The mixed model strategy allowed to collect 59.1% of data from both patients and clinicians; 26 and 14.9% from respectively patients and clinicians only. With the identification of the FSHD1 and FSHD2 forms, specific questionnaires have been designed. Though FSHD2 patients are progressively included, FSHD1 patients still account for the majority (94.9%). The registry is compatible with the FAIR principles as data are Findable, Accessible and Interoperable. We thus used molecular standards and standardized clinical terms used by the FILNEMUS French network of reference centers for the diagnosis and follow-up of patients suffering from a rare neuromuscular disease. The implemented clinical terms mostly map to dictionaries and terminology systems such as SNOMED-CT (75% of terms), CTV3 (61.7%) and NCIt (53.3%). Because of the sensitive nature of data, they are not directly reusable and can only be accessed as aggregated data after evaluation and approval by the registry oversight committee. CONCLUSIONS: The French National Registry of FSHD patients belongs to a national effort to develop databases, which should now interact with other initiatives to build a European and/or an international FSHD virtual registry for the benefits of patients. It is accessible at www.fshd.fr and various useful information, links, and documents, including a video, are available for patients and professionals.

Identification of splicing defects caused by mutations in the dysferlin gene.

Missense, iso-semantic, and intronic mutations are challenging for interpretation, in particular for their impact in mRNA. Various tools such as the Human Splicing Finder (HSF) system could be used to predict the impact on splicing; however, no diagnosis result could rely on predictions alone, but requires functional testing. Here, we report an in vitro approach to study the impact of DYSF mutations on splicing. It was evaluated on a series of 45 DYSF mutations, both intronic and exonic. We confirmed splicing alterations for all intronic mutations localized in 5' or 3' splice sites. Then, we showed that DYSF missense mutations could also result in splicing defects: mutations c.463G>A and c.2641A>C abolished ESEs and led to exon skipping; mutations c.565C>G and c.1555G>A disrupted Exonic Splicing Enhancer (ESE), while concomitantly creating new 5' or 3' splice site leading to exonic out of frame deletions. We demonstrated that 20% of DYSF missense mutations have a strong impact on splicing. This minigene strategy is an efficient tool for the detection of splicing defects in dysferlinopathies, which could allow for a better comprehension of splicing defects due to mutations and could improve prediction tools evaluating splicing defects.

Clinical heterogeneity and a high proportion of novel mutations in a Chinese cohort of patients with dysferlinopathy.

BACKGROUND AND AIMS: Dysferlinopathies are a group of autosomal recessive muscular dystrophies caused by mutations in the dysferlin gene. This study presents clinical features and the mutational spectrum in the largest cohort of Chinese patients analyzed to date. PATIENTS AND METHODS: A total of 36 unrelated Chinese patients with diagnostic suspicion of dysferlinopathy were clinically and genetically characterized. RESULTS: Patients were divided into five phenotypes: 19 patients with limb girdle muscular dystrophy (LGMD) type 2B, 10 with Miyoshi myopathy (MM), 1 with distal anterior compartment myopathy (DACM), 3 with exercise intolerance, and 3 with asymptomatic hypercreatine phosphokinasemia (hyperCPKemia). Thirty-one patients showed an absence or drastic reduction of dysferlin expression by Westernblot. Forty-three mutations were identified in DYSF, including 31 novel. CONCLUSION: Our study underlines clinical heterogeneity and a high proportion of novel mutations in Chinese patients affected with dysferlinopathy.

A human skeletal muscle interactome centered on proteins involved in muscular dystrophies: LGMD interactome.

BACKGROUND: The complexity of the skeletal muscle and the identification of numerous human disease-causing mutations in its constitutive proteins make it an interesting tissue for proteomic studies aimed at understanding functional relationships of interacting proteins in both health and diseases. METHOD: We undertook a large-scale study using two-hybrid screens and a human skeletal-muscle cDNA library to establish a proteome-scale map of protein-protein interactions centered on proteins involved in limb-girdle muscular dystrophies (LGMD). LGMD is a group of more than 20 different neuromuscular disorders that principally affect the proximal pelvic and shoulder girdle muscles. RESULTS AND CONCLUSION: The interaction network we unraveled incorporates 1018 proteins connected by 1492 direct binary interactions and includes 1420 novel protein-protein interactions. Computational, experimental and literature-based analyses were performed to assess the overall quality of this network. Interestingly, LGMD proteins were shown to be highly interconnected, in particular indirectly through sarcomeric proteins. In-depth mining of the LGMD-centered interactome identified new candidate genes for orphan LGMDs and other neuromuscular disorders. The data also suggest the existence of functional links between LGMD2B/dysferlin and gene regulation, between LGMD2C/γ-sarcoglycan and energy control and between LGMD2G/telethonin and maintenance of genome integrity. This dataset represents a valuable resource for future functional investigations.

UMD-DYSF, a novel locus specific database for the compilation and interactive analysis of mutations in the dysferlin gene.

Mutations in the dysferlin gene (DYSF) lead to a complete or partial absence of the dysferlin protein in skeletal muscles and are at the origin of dysferlinopathies, a heterogeneous group of rare autosomal recessive inherited neuromuscular disorders. As a step towards a better understanding of the DYSF mutational spectrum, and towards possible inclusion of patients in future therapeutic clinical trials, we set up the Universal Mutation Database for Dysferlin (UMD-DYSF), a Locus-Specific Database developed with the UMD® software. The main objective of UMD-DYSF is to provide an updated compilation of mutational data and relevant interactive tools for the analysis of DYSF sequence variants, for diagnostic and research purposes. In particular, specific algorithms can facilitate the interpretation of newly identified intronic, missense- or isosemantic-exonic sequence variants, a problem encountered recurrently during genetic diagnosis in dysferlinopathies. UMD-DYSF v1.0 is freely accessible at www.umd.be/DYSF/. It contains a total of 742 mutational entries corresponding to 266 different disease-causing mutations identified in 558 patients worldwide diagnosed with dysferlinopathy. This article presents for the first time a comprehensive analysis of the dysferlin mutational spectrum based on all compiled DYSF disease-causing mutations reported in the literature to date, and using the main bioinformatics tools offered in UMD-DYSF.

Interactions with M-band titin and calpain 3 link myospryn (CMYA5) to tibial and limb-girdle muscular dystrophies.

Mutations in the C terminus of titin, situated at the M-band of the striated muscle sarcomere, cause tibial muscular dystrophy (TMD) and limb-girdle muscular dystrophy (LGMD) type 2J. Mutations in the protease calpain 3 (CAPN3), in turn, lead to LGMD2A, and secondary CAPN3 deficiency in LGMD2J suggests that the pathomechanisms of the diseases are linked. Yeast two-hybrid screens carried out to elucidate the molecular pathways of TMD/LGMD2J and LGMD2A resulted in the identification of myospryn (CMYA5, cardiomyopathy-associated 5) as a binding partner for both M-band titin and CAPN3. Additional yeast two-hybrid and coimmunoprecipitation studies confirmed both interactions. The interaction of myospryn and M-band titin was supported by localization of endogenous and transfected myospryn at the M-band level. Coexpression studies showed that myospryn is a proteolytic substrate for CAPN3 and suggested that myospryn may protect CAPN3 from autolysis. Myospryn is a muscle-specific protein of the tripartite motif superfamily, reported to function in vesicular trafficking and protein kinase A signaling and implicated in the pathogenesis of Duchenne muscular dystrophy. The novel interactions indicate a role for myospryn in the sarcomeric M-band and may be relevant for the molecular pathomechanisms of TMD/LGMD2J and LGMD2A.

The genome of the blood fluke Schistosoma mansoni.

Schistosoma mansoni is responsible for the neglected tropical disease schistosomiasis that affects 210 million people in 76 countries. Here we present analysis of the 363 megabase nuclear genome of the blood fluke. It encodes at least 11,809 genes, with an unusual intron size distribution, and new families of micro-exon genes that undergo frequent alternative splicing. As the first sequenced flatworm, and a representative of the Lophotrochozoa, it offers insights into early events in the evolution of the animals, including the development of a body pattern with bilateral symmetry, and the development of tissues into organs. Our analysis has been informed by the need to find new drug targets. The deficits in lipid metabolism that make schistosomes dependent on the host are revealed, and the identification of membrane receptors, ion channels and more than 300 proteases provide new insights into the biology of the life cycle and new targets. Bioinformatics approaches have identified metabolic chokepoints, and a chemogenomic screen has pinpointed schistosome proteins for which existing drugs may be active. The information generated provides an invaluable resource for the research community to develop much needed new control tools for the treatment and eradication of this important and neglected disease.

The Trypanosoma cruzi L1Tc and NARTc non-LTR retrotransposons show relative site specificity for insertion.

The trypanosomatid protozoan Trypanosoma cruzi contains long autonomous (L1Tc) and short nonautonomous (NARTc) non-long terminal repeat retrotransposons. NARTc (0.25 kb) probably derived from L1Tc (4.9 kb) by 3'-deletion. It has been proposed that their apparent random distribution in the genome is related to the L1Tc-encoded apurinic/apyrimidinic endonuclease (APE) activity, which repairs modified residues. To address this question we used the T. cruzi (CL-Brener strain) genome data to analyze the distribution of all the L1Tc/NARTc elements present in contigs larger than 10 kb. This data set, which represents 0.91x sequence coverage of the haploid nuclear genome ( approximately 55 Mb), contains 419 elements, including 112 full-length L1Tc elements (14 of which are potentially functional) and 84 full-length NARTc. Approximately half of the full-length elements are flanked by a target site duplication, most of them (87%) are 12 bp long. Statistical analyses of sequences flanking the full-length elements show the same highly conserved pattern upstream of both the L1Tc and NARTc retrotransposons. The two most conserved residues are a guanine and an adenine, which flank the site where first-strand cleavage is performed by the element-encoded endonuclease activity. This analysis clearly indicates that the L1Tc and NARTc elements display relative site specificity for insertion, which suggests that the APE activity is not responsible for first-strand cleavage of the target site.

Evolution of non-LTR retrotransposons in the trypanosomatid genomes: Leishmania major has lost the active elements.

The ingi and L1Tc non-LTR retrotransposons--which constitute the ingi clade--are abundant in the genome of the trypanosomatid species Trypanosoma brucei and Trypanosoma cruzi, respectively. The corresponding retroelements, however, are not present in the genome of a closely related trypanosomatid, Leishmania major. To study the evolution of non-LTR retrotransposons in trypanosomatids, we have analyzed all ingi/L1Tc elements and highly degenerate ingi/L1Tc-related sequences identified in the recently completed T. brucei, T. cruzi and L. major genomes. The coding sequences of 242 degenerate ingi/L1Tc-related elements (DIREs) in all three genomes were reconstituted by removing the numerous frame shifts. Three independent phylogenetic analyses conducted on the conserved domains encoded by these elements show that all DIREs, including the 52 L. major DIREs, form a monophyletic group belonging to the ingi clade. This indicates that the trypanosomatid ancestor contained active mobile elements that have been retained in the Trypanosoma species, but were lost from L. major genome, where only remnants (DIRE) are detectable. All 242 DIREs analyzed group together according to their species origin with the exception of 11 T. cruzi DIREs which are close to the T. brucei ingi/DIRE families. Considering the absence of known horizontal transfer between the African T. brucei and the South-American T. cruzi, this suggests that this group of elements evolved at a lower rate when compared to the other trypanosomatid elements. Interestingly, the only nucleotide sequence conserved between ingi and L1Tc (the first 79 residues) is also present at the 5'-extremity of all the full length DIREs and suggests a possible role for this conserved motif, as well as for DIREs.

Comparative genomics of trypanosomatid parasitic protozoa.

A comparison of gene content and genome architecture of Trypanosoma brucei, Trypanosoma cruzi, and Leishmania major, three related pathogens with different life cycles and disease pathology, revealed a conserved core proteome of about 6200 genes in large syntenic polycistronic gene clusters. Many species-specific genes, especially large surface antigen families, occur at nonsyntenic chromosome-internal and subtelomeric regions. Retroelements, structural RNAs, and gene family expansion are often associated with syntenic discontinuities that-along with gene divergence, acquisition and loss, and rearrangement within the syntenic regions-have shaped the genomes of each parasite. Contrary to recent reports, our analyses reveal no evidence that these species are descended from an ancestor that contained a photosynthetic endosymbiont.

The genome sequence of Trypanosoma cruzi, etiologic agent of Chagas disease.

Whole-genome sequencing of the protozoan pathogen Trypanosoma cruzi revealed that the diploid genome contains a predicted 22,570 proteins encoded by genes, of which 12,570 represent allelic pairs. Over 50% of the genome consists of repeated sequences, such as retrotransposons and genes for large families of surface molecules, which include trans-sialidases, mucins, gp63s, and a large novel family (>1300 copies) of mucin-associated surface protein (MASP) genes. Analyses of the T. cruzi, T. brucei, and Leishmania major (Tritryp) genomes imply differences from other eukaryotes in DNA repair and initiation of replication and reflect their unusual mitochondrial DNA. Although the Tritryp lack several classes of signaling molecules, their kinomes contain a large and diverse set of protein kinases and phosphatases; their size and diversity imply previously unknown interactions and regulatory processes, which may be targets for intervention.

The genome of the African trypanosome Trypanosoma brucei.

African trypanosomes cause human sleeping sickness and livestock trypanosomiasis in sub-Saharan Africa. We present the sequence and analysis of the 11 megabase-sized chromosomes of Trypanosoma brucei. The 26-megabase genome contains 9068 predicted genes, including approximately 900 pseudogenes and approximately 1700 T. brucei-specific genes. Large subtelomeric arrays contain an archive of 806 variant surface glycoprotein (VSG) genes used by the parasite to evade the mammalian immune system. Most VSG genes are pseudogenes, which may be used to generate expressed mosaic genes by ectopic recombination. Comparisons of the cytoskeleton and endocytic trafficking systems with those of humans and other eukaryotic organisms reveal major differences. A comparison of metabolic pathways encoded by the genomes of T. brucei, T. cruzi, and Leishmania major reveals the least overall metabolic capability in T. brucei and the greatest in L. major. Horizontal transfer of genes of bacterial origin has contributed to some of the metabolic differences in these parasites, and a number of novel potential drug targets have been identified.

The sequence and analysis of Trypanosoma brucei chromosome II.

We report here the sequence of chromosome II from Trypanosoma brucei, the causative agent of African sleeping sickness. The 1.2-Mb pairs encode about 470 predicted genes organised in 17 directional clusters on either strand, the largest cluster of which has 92 genes lined up over a 284-kb region. An analysis of the GC skew reveals strand compositional asymmetries that coincide with the distribution of protein-coding genes, suggesting these asymmetries may be the result of transcription-coupled repair on coding versus non-coding strand. A 5-cM genetic map of the chromosome reveals recombinational 'hot' and 'cold' regions, the latter of which is predicted to include the putative centromere. One end of the chromosome consists of a 250-kb region almost exclusively composed of RHS (pseudo)genes that belong to a newly characterised multigene family containing a hot spot of insertion for retroelements. Interspersed with the RHS genes are a few copies of truncated RNA polymerase pseudogenes as well as expression site associated (pseudo)genes (ESAGs) 3 and 4, and 76 bp repeats. These features are reminiscent of a vestigial variant surface glycoprotein (VSG) gene expression site. The other end of the chromosome contains a 30-kb array of VSG genes, the majority of which are pseudogenes, suggesting that this region may be a site for modular de novo construction of VSG gene diversity during transposition/gene conversion events.

Molecular evolution of eukaryotic genomes: hemiascomycetous yeast spliceosomal introns.

As part of the exploratory sequencing program Génolevures, visual scrutinisation and bioinformatic tools were used to detect spliceosomal introns in seven hemiascomycetous yeast species. A total of 153 putative novel introns were identified. Introns are rare in yeast nuclear genes (<5% have an intron), mainly located at the 5' end of ORFs, and not highly conserved in sequence. They all share a clear non-random vocabulary: conserved splice sites and conserved nucleotide contexts around splice sites. Homologues of metazoan snRNAs and putative homologues of SR splicing factors were identified, confirming that the spliceosomal machinery is highly conserved in eukaryotes. Several introns' features were tested as possible markers for phylogenetic analysis. We found that intron sizes vary widely within each genome, and according to the phylogenetic position of the yeast species. The evolutionary origin of spliceosomal introns was examined by analysing the degree of conservation of intron positions in homologous yeast genes. Most introns appeared to exist in the last common ancestor of present day yeast species, and then to have been differentially lost during speciation. However, in some cases, it is difficult to exclude a possible sliding event affecting a pre-existing intron or a gain of a novel intron. Taken together, our results indicate that the origin of spliceosomal introns is complex within a given genome, and that present day introns may have resulted from a dynamic flux between intron conservation, intron loss and intron gain during the evolution of hemiascomycetous yeasts.