Expansion of GGC Repeat in GIPC1 Is Associated with Oculopharyngodistal Myopathy.

Oculopharyngodistal myopathy (OPDM) is an adult-onset inherited neuromuscular disorder characterized by progressive ptosis, external ophthalmoplegia, and weakness of the masseter, facial, pharyngeal, and distal limb muscles. The myopathological features are presence of rimmed vacuoles (RVs) in the muscle fibers and myopathic changes of differing severity. Inheritance is variable, with either putative autosomal-dominant or autosomal-recessive pattern. Here, using a comprehensive strategy combining whole-genome sequencing (WGS), long-read whole-genome sequencing (LRS), linkage analysis, repeat-primed polymerase chain reaction (RP-PCR), and fluorescence amplicon length analysis polymerase chain reaction (AL-PCR), we identified an abnormal GGC repeat expansion in the 5' UTR of GIPC1 in one out of four families and three sporadic case subjects from a Chinese OPDM cohort. Expanded GGC repeats were further confirmed as the cause of OPDM in an additional 2 out of 4 families and 6 out of 13 sporadic Chinese individuals with OPDM, as well as 7 out of 194 unrelated Japanese individuals with OPDM. Methylation, qRT-PCR, and western blot analysis indicated that GIPC1 mRNA levels were increased while protein levels were unaltered in OPDM-affected individuals. RNA sequencing indicated p53 signaling, vascular smooth muscle contraction, ubiquitin-mediated proteolysis, and ribosome pathways were involved in the pathogenic mechanisms of OPDM-affected individuals with GGC repeat expansion in GIPC1. This study provides further evidence that OPDM is associated with GGC repeat expansions in distinct genes and highly suggests that expanded GGC repeat units are essential in the pathogenesis of OPDM, regardless of the genes in which the expanded repeats are located.

Single-cell RNA-seq analysis of mouse preimplantation embryos by third-generation sequencing.

The development of next generation sequencing (NGS) platform-based single-cell RNA sequencing (scRNA-seq) techniques has tremendously changed biological researches, while there are still many questions that cannot be addressed by them due to their short read lengths. We developed a novel scRNA-seq technology based on third-generation sequencing (TGS) platform (single-cell amplification and sequencing of full-length RNAs by Nanopore platform, SCAN-seq). SCAN-seq exhibited high sensitivity and accuracy comparable to NGS platform-based scRNA-seq methods. Moreover, we captured thousands of unannotated transcripts of diverse types, with high verification rate by reverse transcription PCR (RT-PCR)-coupled Sanger sequencing in mouse embryonic stem cells (mESCs). Then, we used SCAN-seq to analyze the mouse preimplantation embryos. We could clearly distinguish cells at different developmental stages, and a total of 27,250 unannotated transcripts from 9,338 genes were identified, with many of which showed developmental stage-specific expression patterns. Finally, we showed that SCAN-seq exhibited high accuracy on determining allele-specific gene expression patterns within an individual cell. SCAN-seq makes a major breakthrough for single-cell transcriptome analysis field.

Genetic origin of sporadic cases and RNA toxicity in neuronal intranuclear inclusion disease.

BACKGROUND: GGC repeat expansion in NOTCH2NLC has been recently linked to neuronal intranuclear inclusion disease (NIID) via unknown disease mechanisms. Herein, we explore the genetic origin of the sporadic cases and toxic RNA gain-of-function mechanism in NIID. METHODS: Multiple genetic screenings were performed on NIID individuals and their available family members. Methylation status of blood DNA, NOTCH2NLC mRNA level from muscle biopsies and RNA foci from skin biopsies of NIID individuals or asymptomatic carriers were evaluated and compared. RESULTS: In two sporadic NIID families, we identified two clinically and pathologically asymptomatic fathers carrying large GGC repeat expansion, above 300 repeats, with offspring repeat numbers of 172 and 148, respectively. Further evaluation revealed that the GGC repeat numbers in the sperm from two asymptomatic fathers were only 63 and 98, respectively. The CpG island in NOTCH2NLC of the asymptomatic carriers was hypermethylated, and accordingly, the NOTCH2NLC mRNA levels were decreased in the asymptomatic fathers. GGC repeat expansion RNA formed RNA foci and sequestered RNA binding proteins into p62 positive intranuclear inclusions in NIID individuals but not in the control or asymptomatic carrier. CONCLUSION: Our study suggested the GGC repeat expansion in NOTCH2NLC might have a disease-causing number ranging from ~41 to ~300 repeats. The contraction of GGC repeat expansion in sperm could be a possible mechanism for the paternal-biased origin in some sporadic or recessive inherited NIID individuals. The toxic RNA gain-of-function mechanism was identified to be involved in the pathogenicity of this disease.

The GGC repeat expansion in NOTCH2NLC is associated with oculopharyngodistal myopathy type 3.

Oculopharyngodistal myopathy (OPDM) is an adult-onset neuromuscular disease characterized by progressive ocular, facial, pharyngeal and distal limb muscle involvement. Trinucleotide repeat expansions in LRP12 or GIPC1 were recently reported to be associated with OPDM. However, a significant portion of OPDM patients have unknown genetic causes. In this study, long-read whole-genome sequencing and repeat-primed PCR were performed and we identified GGC repeat expansions in the NOTCH2NLC gene in 16.7% (4/24) of a cohort of Chinese OPDM patients, designated as OPDM type 3 (OPDM3). Methylation analysis indicated that methylation levels of the NOTCH2NLC gene were unaltered in OPDM3 patients, but increased significantly in asymptomatic carriers. Quantitative real-time PCR analysis indicated that NOTCH2NLC mRNA levels were increased in muscle but not in blood of OPDM3 patients. Immunofluorescence on OPDM muscle samples and expressing mutant NOTCH2NLC with (GGC)69 repeat expansions in HEK293 cells indicated that mutant NOTCH2NLC-polyglycine protein might be a major component of intranuclear inclusions, and contribute to toxicity in cultured cells. In addition, two RNA-binding proteins, hnRNP A/B and MBNL1, were both co-localized with p62 in intranuclear inclusions in OPDM muscle samples. These results indicated that a toxic protein gain-of-function mechanism and RNA gain-of-function mechanism may both play a vital role in the pathogenic processes of OPDM3. This study extended the spectrum of NOTCH2NLC repeat expansion-related diseases to a predominant myopathy phenotype presenting as OPDM, and provided evidence for possible pathogenesis of these diseases.

Single-molecule optical mapping enables quantitative measurement of D4Z4 repeats in facioscapulohumeral muscular dystrophy (FSHD).

PURPOSE: Facioscapulohumeral muscular dystrophy (FSHD) is a common adult muscular dystrophy. Over 95% of FSHD cases are associated with contraction of the D4Z4 tandem repeat (~3.3 kb per unit) at 4q35 with a specific genomic configuration (haplotype) called 4qA. Molecular diagnosis of FSHD typically requires pulsed-field gel electrophoresis with Southern blotting. We aim to develop novel genomic and computational methods for characterising D4Z4 repeat numbers in FSHD. METHODS: We leveraged a single-molecule optical mapping platform that maps locations of restriction enzyme sites on high molecular weight (>150 kb) DNA molecules. We developed bioinformatics methods to address several challenges, including the differentiation of 4qA with 4qB alleles, the differentiation of 4q35 and 10q26 segmental duplications, the quantification of repeat numbers with different enzymes that may or may not have recognition sites within D4Z4 repeats. We evaluated the method on 25 human subjects (13 patients, 3 individual control subjects, 9 control subjects from 3 families) labelled by the Nb.BssSI and/or Nt.BspQI enzymes. RESULTS: We demonstrated that the method gave a direct quantitative measurement of repeat numbers on D4Z4 repeats with 4qA allelic configuration and the levels of postzygotic mosaicism. Our method had high concordance with Southern blots from several cohorts on two platforms (Bionano Saphyr and Bionano Irys), but with improved quantification of repeat numbers. CONCLUSION: While the study is limited by small sample size, our results demonstrated that single-molecule optical mapping is a viable approach for more refined analysis on genotype-phenotype relationships in FSHD, especially when postzygotic mosaicism is present.

Rapid prenatal diagnosis of Facioscapulohumeral Muscular Dystrophy 1 by combined Bionano optical mapping and karyomapping.

PURPOSE: To explore the feasibility of performing rapid prenatal diagnoses of FSHD1 using a combination of Bianano optical mapping and linkage-based karyomapping. METHODS: DNA specimens from a family that had been previously diagnosed with FSHD1 using Southern Blot analysis were used for this study. Genetic diagnosis of the proband, fetus chorionic amniotic fluid, and aborted fetal tissue was performed using Bianano optical mapping (BOM) together with linkage-based karyomapping. RESULTS: BOM analysis showed that the proband's 4q35.2 region contained four D4Z4 repeats and the 4qA permissible allele, consistent with the previous FSHD1 diagnosis obtained by Southern Blotting. BOM analysis of the fetus' 4q35.2 region was consistent with that of the proband. Karyomap analysis revealed that the fetus inherited the affected chromosome segment from the proband. After genetic counseling, the couple choose termination of pregnancy, and we performed gene diagnosis of the abortus tissue by BOM. CONCLUSIONS: Bianano optical mapping can determine the number of D4Z4 repeats and exclude interference of the 10q26.3 homologous region, and in combination with karyomapping, can be used for rapid and accurate prenatal diagnosis of FSHD1.

Long-read sequencing identified repeat expansions in the 5'UTR of the NOTCH2NLC gene from Chinese patients with neuronal intranuclear inclusion disease.

BACKGROUND: Neuronal intranuclear inclusion disease (NIID) is a heterogenous neurodegenerative disorder named after its pathological features. It has long been considered a disease of genetic origin. Recently, the GGC repeated expansion in the 5'-untranslated region (5'UTR) of the NOTCH2NLC gene has been found in adult-onset NIID in Japanese individuals. This study was aimed to investigate the causative mutations of NIID in Chinese patients. METHODS: Fifteen patients with NIID were identified from five academic neurological centres. Biopsied skin samples were analysed by histological staining, immunostaining and electron microscopic observation. Whole-genome sequencing (WGS) and long-read sequencing (LRS) were initially performed in three patients with NIID. Repeat-primed PCR was conducted to confirm the genetic variations in the three patients and the other 12 cases. RESULTS: Our patients included 14 adult-onset patients and 1 juvenile-onset patient characterised by degeneration of multiple nervous systems. All patients were identified with intranuclear inclusions in the nuclei of fibroblasts, fat cells and ductal epithelial cells of sweat glands. The WGS failed to find any likely pathogenic variations for NIID. The LRS successfully identified that three patients with adult-onset NIID showed abnormalities of GGC expansion in 5'UTR of the NOTCH2NLC gene. The GGC repeated expansion was further confirmed by repeat-primed PCR in seven familial cases and eight sporadic cases. CONCLUSION: Our findings provided evidence that confirmed the GGC repeated expansion in the 5'UTR of the NOTCH2NLC gene is associated with the pathogenesis of NIID. Additionally, the GGC expansion was not only responsible for adult-onset patients, but also responsible for juvenile-onset patients.

Drug-symptom networking: Linking drug-likeness screening to drug discovery.

Understanding the relationships between drugs and symptoms has broad medical consequences, yet a comprehensive description of the drug-symptom associations is currently lacking. Here, 1441 FDA-approved drugs were collected, and PCA was used to extract 122 descriptors which explained 91% of the variance. Then, a k-means++ method was employed to partition the drug dataset into 3 clusters, and 3 corresponding SVDD models (drug-likeness screening models) were constructed with an overall accuracy of up to 95.6%. Furthermore, 6878 herbal molecules from the TcmSP™ database were screened by the above 3 SVDD model to obtain 5309 candidate drug molecules with highly accept classification of 77.19%. To assess the accuracy of the SVDD models, 8559 herbal molecule-symptom co-occurrences were mined from Pubmed abstracts, involving 697 herbal molecules and 314 symptoms. Most of the 697 herbal molecules could be found in the accepted SVDD data (5309 molecules), showing the potential of the SVDD for the screening of drug candidates. Moreover, a herbal molecule-herbal molecule network and a herbal molecule-symptom were constructed. Overall, the results provided a new drug-likeness screening approach independent to abnormal training data, and the comprehensive collection of herbal molecule-symptom associations formed a new data resource for systematic characterization of the symptom-oriented medicines.

Identification and characterization of two DMD pedigrees with large inversion mutations based on a long-read sequencing pipeline.

Pathogenic large inversions are rarely reported on DMD gene due to the lack of effective detection methods. Here we report two DMD pedigrees and proposed a reliable pipeline to define large inversions in DMD patients. In the first pedigree, conventional approaches including multiplex ligation-dependent probe amplification, and whole-exome sequencing by next generation sequencing were failed to detect any pathologic variant. Then an advanced analysis pipeline which consists of RNA-seq, cDNA array capture sequencing, optical mapping, long-read sequencing was built. RNA-seq and cDNA capture sequencing showed a complete absence of transcripts of exons 3-55. Optical mapping identified a 55 Mb pericentric inversion between Xp21 and Xq21. Subsequently, long-read sequencing and Sanger sequencing determined the inversion breakpoints at 32,915,769 and 87,989,324 of the X chromosomes. In the second pedigree, long-read sequencing was directly conducted and Sanger sequencing was performed to verify the mutation. Long-read sequencing and Sanger sequencing found breakpoints at 32,581,576 and 127,797,236 on DMD gene directly. In conclusion, large inversion might be a rare but important mutation type in DMD gene. An effective pipeline was built in detecting large inversion mutations based on long-read sequencing platforms.

A large pedigree study confirmed the CGG repeat expansion of RILPL1 Is associated with oculopharyngodistal myopathy.

BACKGROUND: Oculopharyngodistal myopathy (OPDM) is an autosomal dominant adult-onset degenerative muscle disorder characterized by ptosis, ophthalmoplegia and weakness of the facial, pharyngeal and limb muscles. Trinucleotide repeat expansions in non-coding regions of LRP12, G1PC1, NOTCH2NLC and RILPL1 were reported to be the etiologies for OPDM. RESULTS: In this study, we performed long-read whole-genome sequencing in a large five-generation family of 156 individuals, including 21 patients diagnosed with typical OPDM. We identified CGG repeat expansions in 5'UTR of RILPL1 gene in all patients we tested while no CGG expansion in unaffected family members. Repeat-primed PCR and fluorescence amplicon length analysis PCR were further confirmed the segregation of CGG expansions in other family members and 1000 normal Chinese controls. Methylation analysis indicated that methylation levels of the RILPL1 gene were unaltered in OPDM patients, which was consistent with previous studies. Our findings provide evidence that RILPL1 is associated OPDM in this large pedigree. CONCLUSIONS: Our results identified RILPL1 is the associated the disease in this large pedigree.

Amniotes co-opt intrinsic genetic instability to protect germ-line genome integrity.

Unlike PIWI-interacting RNA (piRNA) in other species that mostly target transposable elements (TEs), >80% of piRNAs in adult mammalian testes lack obvious targets. However, mammalian piRNA sequences and piRNA-producing loci evolve more rapidly than the rest of the genome for unknown reasons. Here, through comparative studies of chickens, ducks, mice, and humans, as well as long-read nanopore sequencing on diverse chicken breeds, we find that piRNA loci across amniotes experience: (1) a high local mutation rate of structural variations (SVs, mutations ≥ 50 bp in size); (2) positive selection to suppress young and actively mobilizing TEs commencing at the pachytene stage of meiosis during germ cell development; and (3) negative selection to purge deleterious SV hotspots. Our results indicate that genetic instability at pachytene piRNA loci, while producing certain pathogenic SVs, also protects genome integrity against TE mobilization by driving the formation of rapid-evolving piRNA sequences.

Study of complex structural variations of X-linked deafness-2 based on single-molecule sequencing.

X-linked deafness-2 (DFNX2) is cochlear incomplete partition type III (IP-III), one of inner ear malformations characterized by an abnormally wide opening in the bone separating the basal turn of the cochlea from the internal auditory canal, fixation of the stapes and cerebrospinal fluid (CSF) gusher upon stapedectomy or cochleostomy. The causative gene of DFNX2 was POU3F4. To investigate the genetic causes of DFNX2 and compare the efficiency of different sequencing methods, 12 unrelated patients were enrolled in the present study. Targeted next-generation sequencing (NGS) and long-read sequencing were used to analyze the genetic etiology of DFNX2. Six variants of POU3F4 were identified in this cohort by NGS. Three patients with a negative diagnosis based on NGS were enrolled in further long-read sequencing. Two of them were all found to carry structural variations (SVs) on chromosome X, consisting of an 870-kb deletion (DEL) at upstream of POU3F4 and an 8-Mb inversion (INV). The 870-kb DEL may have been due to non-homologous end joining (NHEJ), while non-allelic homologous recombination (NAHR) within a single chromatid may have accounted for the 8-Mb INV. Common POU3F4 mutations in DFNX2 included point mutations, small insertions and deletions (INDELs), and exon mutations, which can be detected by Sanger sequencing and NGS. Single-molecule long-read sequencing constitutes an additional and valuable method for accurate detection of pathogenic SVs in IP-III patients with negative NGS results.

GGC repeat expansions in NOTCH2NLC causing a phenotype of distal motor neuropathy and myopathy.

BACKGROUND: The expansion of GGC repeat in the 5' untranslated region of the NOTCH2NLC has been associated with various neurogenerative disorders of the central nervous system and, more recently, oculopharyngodistal myopathy. This study aimed to report patients with distal weakness with both neuropathic and myopathic features on electrophysiology and pathology who present GGC repeat expansions in the NOTCH2NLC. METHODS: Whole-exome sequencing (WES) and long-read sequencing were implemented to identify the candidate genes. In addition, the available clinical data and the pathological changes associated with peripheral nerve and muscle biopsies were reviewed and studied. RESULTS: We identified and validated GGC repeat expansions of NOTCH2NLC in three unrelated patients who presented with progressive weakness predominantly affecting distal lower limb muscles, following negative results in an initial WES. We found intranuclear inclusions with multiple proteins deposits in the nuclei of both myofibers and Schwann cells. The clinical features of these patients are compatible with the diagnosis of distal motor neuropathy and rimmed vacuolar myopathy. INTERPRETATION: These phenotypes enrich the class of features associated with NOTCH2NLC-related repeat expansion disorders (NRED), and provide further evidence that the neurological symptoms of NRED include not only brain, spinal cord, and peripheral nerves damage, but also myopathy, and that overlapping symptoms might exist.

Insights from systems pharmacology into cardiovascular drug discovery and therapy.

BACKGROUND: Given the complex nature of cardiovascular disease (CVD), information derived from a systems-level will allow us to fully interrogate features of CVD to better understand disease pathogenesis and to identify new drug targets. RESULTS: Here, we describe a systematic assessment of the multi-layer interactions underlying cardiovascular drugs, targets, genes and disorders to reveal comprehensive insights into cardiovascular systems biology and pharmacology. We have identified 206 effect-mediating drug targets, which are modulated by 254 unique drugs, of which, 43% display activities across different protein families (sequence similarity < 30%), highlighting the fact that multitarget therapy is suitable for CVD. Although there is little overlap between cardiovascular protein targets and disease genes, the two groups have similar pleiotropy and intimate relationships in the human disease gene-gene and cellular networks, supporting their similar characteristics in disease development and response to therapy. We also characterize the relationships between different cardiovascular disorders, which reveal that they share more etiological commonalities with each other rooted in the global disease-disease networks. Furthermore, the disease modular analysis demonstrates apparent molecular connection between 227 cardiovascular disease pairs. CONCLUSIONS: All these provide important consensus as to the cause, prevention, and treatment of various CVD disorders from systems-level perspective.

TCMSP: a database of systems pharmacology for drug discovery from herbal medicines.

BACKGROUND: Modern medicine often clashes with traditional medicine such as Chinese herbal medicine because of the little understanding of the underlying mechanisms of action of the herbs. In an effort to promote integration of both sides and to accelerate the drug discovery from herbal medicines, an efficient systems pharmacology platform that represents ideal information convergence of pharmacochemistry, ADME properties, drug-likeness, drug targets, associated diseases and interaction networks, are urgently needed. DESCRIPTION: The traditional Chinese medicine systems pharmacology database and analysis platform (TCMSP) was built based on the framework of systems pharmacology for herbal medicines. It consists of all the 499 Chinese herbs registered in the Chinese pharmacopoeia with 29,384 ingredients, 3,311 targets and 837 associated diseases. Twelve important ADME-related properties like human oral bioavailability, half-life, drug-likeness, Caco-2 permeability, blood-brain barrier and Lipinski's rule of five are provided for drug screening and evaluation. TCMSP also provides drug targets and diseases of each active compound, which can automatically establish the compound-target and target-disease networks that let users view and analyze the drug action mechanisms. It is designed to fuel the development of herbal medicines and to promote integration of modern medicine and traditional medicine for drug discovery and development. CONCLUSIONS: The particular strengths of TCMSP are the composition of the large number of herbal entries, and the ability to identify drug-target networks and drug-disease networks, which will help revealing the mechanisms of action of Chinese herbs, uncovering the nature of TCM theory and developing new herb-oriented drugs. TCMSP is freely available at http://sm.nwsuaf.edu.cn/lsp/tcmsp.php.