Transcriptome profiling identifies regulators of pathogenesis in collagen VI related muscular dystrophy

Transcriptome profiling of skeletal muscles from Korean patients with Bethlem myopathy

Bethlem myopathy is one of the collagens VI-related muscular dystrophies caused by mutations in the collagen VI genes. The study was designed to analyze the gene expression profiles in the skeletal muscle of patients with Bethlem myopathy. Six skeletal muscle samples from 3 patients with Bethlem myopathy and 3 control subjects were analyzed by RNA-sequencing. 187 transcripts were significantly differentially expressed, with 157 upregulated and 30 downregulated transcripts in the Bethlem group. Particularly, 1 (microRNA-133b) was considerably upregulated, and 4 long intergenic non-protein coding RNAs, LINC01854, MBNL1-AS1, LINC02609, and LOC728975, were significantly downregulated. We categorized differentially expressed gene using Gene Ontology and showed that Bethlem myopathy is strongly associated with the organization of extracellular matrix (ECM). Kyoto Encyclopedia of Genes and Genomes pathway enrichment reflected themes with significant enrichment of the ECM-receptor interaction (hsa04512), complement and coagulation cascades (hsa04610), and focal adhesion (hsa04510). We confirmed that Bethlem myopathy is strongly associated with the organization of ECM and the wound healing process. Our results demonstrate transcriptome profiling of Bethlem myopathy, and provide new insights into the path mechanism of Bethlem myopathy associated with non-protein coding RNAs.

COLLAGEN VI-RELATED MYOPATHIES: CLINICAL VARIABILITY, PHENOTYPE-GENOTYPE CORRELATION AND EXPLORATORY TRANSCRIPTOME STUDY

Collagen VI-related myopathies are a group of disorders that cause muscle weakness and joint contractures with significant variability in disease severity among patients. Here we report the clinical and genetic characteristics of 13 Chinese patients. Detailed histological, radiological and muscle transcriptomic evaluations were also conducted for selected representative patients. Across the cohort, fifteen putative disease causal variants were identified in three genes encoding collagen VI subunits, COL6A1 (n=6), COL6A2 (n=5), and COL6A3 (n=4). Most of these variants (12/15, 80%) were dominant negative and occurred at the triple helical domain. The rest (3/15, 20%) were located at the C-terminus. Two previously unreported variants, an in-frame mutation (COL6A1:c.1084_1092del) and a missense mutation (COL6A2:c.811G>C), were also noted. The transcriptome data from the muscle biopsies of two patients in the study with dominant negative mutations [COL6A2:c.811G>C and COL6A1:c.930+189C>T] supports the accepted aetiology of Collagen VI myopathy as dysfunction of the extracellular matrix. It also suggests there are perturbations to skeletal muscle differentiation and skeletal system development. It should be noted that although the phenotypes of patients can be mostly explained by the position and dominant-negative effect of the variants, exceptions and variability still exist and have to be reckoned with. This study provides valuable data explaining the varying severity of phenotypes among ethnically Chinese patients.

Reversal of viral and epigenetic HLA class I repression in Merkel cell carcinoma

Cancers avoid immune surveillance through an array of mechanisms, including perturbation of HLA class I antigen presentation. Merkel cell carcinoma (MCC) is an aggressive, HLA-I-low, neuroendocrine carcinoma of the skin often caused by the Merkel cell polyomavirus (MCPyV). Through the characterization of 11 newly generated MCC patient-derived cell lines, we identified transcriptional suppression of several class I antigen presentation genes. To systematically identify regulators of HLA-I loss in MCC, we performed parallel, genome-scale, gain- and loss-of-function screens in a patient-derived MCPyV-positive cell line and identified MYCL and the non-canonical Polycomb repressive complex 1.1 (PRC1.1) as HLA-I repressors. We observed physical interaction of MYCL with the MCPyV small T viral antigen, supporting a mechanism of virally mediated HLA-I suppression. We further identify the PRC1.1 component USP7 as a pharmacologic target to restore HLA-I expression in MCC.

Collagen misfolding mutations: the contribution of the unfolded protein response to the molecular pathology

Mutations in collagen genes cause a broad range of connective tissue pathologies. Structural mutations that impact procollagen assembly or triple helix formation and stability are a common and important mutation class. How misfolded procollagens engage with the cellular proteostasis machinery and whether they can elicit a cytotoxic unfolded protein response (UPR) is a topic of considerable research interest. Such interest is well justified since modulating the UPR could offer a new approach to treat collagenopathies for which there are no current disease mechanism-targeting therapies. This review scrutinizes the evidence underpinning the view that endoplasmic reticulum stress and chronic UPR activation contributes significantly to the pathophysiology of the collagenopathies. While there is strong evidence that the UPR contributes to the pathology for collagen X misfolding mutations, the evidence that misfolding mutations in other collagen types induce a canonical, cytotoxic UPR is incomplete. To gain a more comprehensive understanding about how the UPR amplifies to pathology, and thus what types of manipulations of the UPR might have therapeutic relevance, much more information is needed about how specific misfolding mutation types engage differentially with the UPR and downstream signaling responses. Most importantly, since the capacity of the proteostasis machinery to respond to collagen misfolding is likely to vary between cell types, reflecting their functional roles in collagen and extracellular matrix biosynthesis, detailed studies on the UPR should focus as much as possible on the actual target cells involved in the collagen pathologies.

Caspase-9 inhibition decreases expression of Mmp9 during chondrogenesis

Besides cell death, caspase-9 participates in non-apoptotic events, including cell differentiation. To evaluate a possible impact on the expression of chondrogenic/osteogenic factors, a caspase-9 inhibitor was tested in vitro. For this purpose, mouse forelimb-derived micromass cultures, the most common chondrogenic in vitro model, were used. The following analyses were performed based on polymerase chain reaction (PCR) arrays and real-time PCR. The expression of several chondrogenesis-related genes was shown to be altered, some of which may impact chondrogenic differentiation (Bmp4, Bmp7, Sp7, Gli1), mineral deposition (Alp, Itgam) or the remodelling of the extracellular matrix (Col1a2, Mmp9) related to endochondral ossification. From the cluster of genes with altered expression, Mmp9 showed the most significant decrease in expression, of more than 50-fold. Additionally, we determined the possible impact of caspase-9 downregulation on the expression of other Mmp genes. A mild increase in Mmp14 was observed, but there was no change in the expression of other studied Mmp genes (-2, -3, -8, -10, -12, -13). Interestingly, inhibition of Mmp9 in micromasses led to decreased expression of some chondrogenic markers related to caspase-9. These samples also showed a decreased expression of caspase-9 itself, suggesting a bidirectional regulation of these two enzymes. These results indicate a specific impact of caspase-9 inhibition on the expression of Mmp9. The localisation of these two enzymes overlaps in resting, proliferative and pre-hypertrophic chondrocytes during in vivo development, which supports their multiple functions, either apoptotic or non-apoptotic. Notably, a coincidental expression pattern was identified in Pik3cg, a possible candidate for Mmp9 regulation.

Transcriptome analysis of collagen VI-related muscular dystrophy muscle biopsies

Objective

To define the transcriptomic changes responsible for the histologic alterations in skeletal muscle and their progression in collagen VI‐related muscular dystrophy (COL6‐RD).

Methods

COL6‐RD patient muscle biopsies were stratified into three groups based on the overall level of pathologic severity considering degrees of fibrosis, muscle fiber atrophy, and fatty replacement of muscle tissue. Using microarray and RNA‐Seq, we then performed global gene expression profiling on the same muscle biopsies and compared their transcriptome with age‐ and sex‐matched controls.

Results

COL6‐RD muscle biopsy transcriptomes as a group revealed prominent upregulation of muscle extracellular matrix component genes and the downregulation of skeletal muscle and mitochondrion‐specific genes. Upregulation of the TGFβ pathway was the most conspicuous change across all biopsies and was fully evident even in the mildest/earliest histological group. There was no difference in the overall transcriptional signature between the different histologic groups but polyserial analysis identified relative changes along with COL6‐RD histological severity.

Interpretation

Overall, our study establishes the prominent dysregulation of extracellular matrix genes, TGFβ signaling, and its downstream cellular pathways at the transcriptomic level in COL6‐RD muscle.

A new single nucleotide polymorphism affects the predisposition to thoracic ossification of the posterior longitudinal ligament

Background

Thoracic ossification of the posterior longitudinal ligament (T-OPLL) is one of the common factors that cause thoracic spinal stenosis, which results in intractable myelopathy and radiculopathy. Our previous study first reported rs201153092A site mutation in the collagen 6A1 (COL6A1) gene as a potentially pathogenic locus for T-OPLL. We aimed to determine whether the rs201153092A site mutation causes abnormal expression of the COL6A1 in Han Chinese patients with T-OPLL and whether this locus is also associated with cervical-OPLL.

Methods

Peripheral blood was collected from a total of 60 patients with T-OPLL disease (30 patients carrying the rs201153092A site mutation in COL6A1 and 30 wild-type patients) and 400 northern Chinese individuals (200 cervical-OPLL patients and 200 control subjects) using the Sequenom system. The expression of COL6A1 was analyzed by enzyme-linked immunosorbent assay, reverse transcription-quantitative polymerase chain reaction, and Western blotting.

Results

rs201153092A mutation resulted in markedly increased COL6A1 gene expression levels in peripheral blood samples. The allele frequency and genotype frequency results showed that this locus is no difference between cervical-OPLL patients and controls.

Conclusions

The rs201153092A site mutation of COL6A1 can significantly increase the expression of COL6A1. The COL6A1 gene rs201153092A site polymorphism is a potential pathogenic mutation in T-OPLL disease, which may be only associated with the occurrence of T-OPLL.

AAV-Mediated Gene Transfer Restores a Normal Muscle Transcriptome in a Canine Model of X-Linked Myotubular Myopathy

Multiple clinical trials employing recombinant adeno-associated viral (rAAV) vectors have been initiated for neuromuscular disorders, including Duchenne and limb-girdle muscular dystrophies, spinal muscular atrophy, and recently X-linked myotubular myopathy (XLMTM). Our previous work on a canine model of XLMTM showed that a single rAAV8-cMTM1 systemic infusion corrected structural abnormalities within the muscle and restored contractile function, with affected dogs surviving more than 4 years post injection. This remarkable therapeutic efficacy presents a unique opportunity to identify the downstream molecular drivers of XLMTM pathology and to what extent the whole muscle transcriptome is restored to normal after gene transfer. Herein, RNA-sequencing was used to examine the transcriptomes of the Biceps femoris and Vastus lateralis in a previously described canine cohort that showed dose-dependent clinical improvements after rAAV8-cMTM1 gene transfer. Our analysis confirmed several dysregulated genes previously observed in XLMTM mice but also identified transcripts linked to XLMTM pathology. We demonstrated XLMTM transcriptome remodeling and dose-dependent normalization of gene expression after gene transfer and created metrics to pinpoint potential biomarkers of disease progression and correction.

Etiopathogenetic Mechanisms in Diverticular Disease of the Colon

This article reviews epidemiological evidence of heritability and putative mechanisms in diverticular disease, with greatest attention to 3 recent studies of genetic associations with diverticular disease based on genome-wide or whole-genome sequencing studies in large patient cohorts. We provide an analysis of the biological plausibility of the significant associations with gene variants reported and highlight the relevance of ANO1, CPI-17 (aka PPP1R14A), COLQ6, COL6A1, CALCB or CALCA, COL6A1, ARHGAP15, and S100A10 to colonic neuromuscular function and tissue properties that may result in altered compliance and predispose to the development of diverticular disease. Such studies also identify candidate genes for future studies.

Transcriptome Profiling of Primary Skin Fibroblasts Reveal Distinct Molecular Features Between PLOD1- and FKBP14-Kyphoscoliotic Ehlers-Danlos Syndrome

Kyphoscoliotic Ehlers-Danlos Syndrome (kEDS) is a rare genetic heterogeneous disease clinically characterized by congenital muscle hypotonia, kyphoscoliosis, and joint hypermobility. kEDS is caused by biallelic pathogenic variants in either PLOD1 or FKBP14. PLOD1 encodes the lysyl hydroxylase 1 enzyme responsible for hydroxylating lysyl residues in the collagen helix, which undergo glycosylation and form crosslinks in the extracellular matrix thus contributing to collagen fibril strength. FKBP14 encodes a peptidyl-prolyl cis-trans isomerase that catalyzes collagen folding and acts as a chaperone for types III, VI, and X collagen. Despite genetic heterogeneity, affected patients with mutations in either PLOD1 or FKBP14 are clinically indistinguishable. We aim to better understand the pathomechanism of kEDS to characterize distinguishing and overlapping molecular features underlying PLOD1-kEDS and FKBP14-kEDS, and to identify novel molecular targets that may expand treatment strategies. Transcriptome profiling by RNA sequencing of patient-derived skin fibroblasts revealed differential expression of genes encoding extracellular matrix components that are unique between PLOD1-kEDS and FKBP14-kEDS. Furthermore, we identified genes involved in inner ear development, vascular remodeling, endoplasmic reticulum (ER) stress, and protein trafficking that were differentially expressed in patient fibroblasts compared to controls. Overall, our study presents the first transcriptomics data in kEDS revealing distinct molecular features between PLOD1-kEDS and FKBP14-kEDS, and serves as a tool to better understand the disease.

Novel defects in collagen XII and VI expand the mixed myopathy/Ehlers-Danlos syndrome spectrum and lead to variant-specific alterations in the extracellular matrix

PURPOSE

To date, heterozygous or homozygous COL12A1 variants have been reported in 13 patients presenting with a clinical phenotype overlapping with collagen VI-related myopathies and Ehlers-Danlos syndrome (EDS). The small number of reported patients limits thorough investigation of this newly identified syndrome, currently coined as myopathic EDS.

METHODS

DNA from 78 genetically unresolved patients fulfilling the clinical criteria for myopathic EDS was sequenced using a next-generation panel of COL12A1, COL6A1, COL6A2, and COL6A3.

RESULTS

Among this cohort, we identified four pathogenic heterozygous in-frame exon skipping (∆) defects in COL12A1, clustering to the thrombospondin N-terminal region and the adjacent collagenous domain (Δ52, Δ53, Δ54, and Δ56 respectively), one heterozygous COL12A1 arginine-to-cysteine substitution of unclear significance (p.(Arg1863Cys)), and compound heterozygous pathogenic COL6A1 variants (c.[98-6G>A];[301C>T]) in one proband. Variant-specific intracellular accumulation of collagen XII chains, extracellular overmodification of the long isoform and near-absence of the short isoform of collagen XII, and extracellular decrease of decorin and tenascin-X were observed for the COL12A1 variants. In contrast, the COL6A1 variants abolished collagen VI and V deposition and increased tenascin-X levels.

CONCLUSION

Our data further support the significant clinical overlap between myopathic EDS and collagen VI-related myopathies, and emphasize the variant-specific consequences of collagen XII defects.

Expanding the Boundaries of RNA Sequencing as a Diagnostic Tool for Rare Mendelian Disease

Gene-panel and whole-exome analyses are now standard methodologies for mutation detection in Mendelian disease. However, the diagnostic yield achieved is at best 50%, leaving the genetic basis for disease unsolved in many individuals. New approaches are thus needed to narrow the diagnostic gap. Whole-genome sequencing is one potential strategy, but it currently has variant-interpretation challenges, particularly for non-coding changes. In this study we focus on transcriptome analysis, specifically total RNA sequencing (RNA-seq), by using monogenetic neuromuscular disorders as proof of principle. We examined a cohort of 25 exome and/or panel "negative" cases and provided genetic resolution in 36% (9/25). Causative mutations were identified in coding and non-coding exons, as well as in intronic regions, and the mutational pathomechanisms included transcriptional repression, exon skipping, and intron inclusion. We address a key barrier of transcriptome-based diagnostics: the need for source material with disease-representative expression patterns. We establish that blood-based RNA-seq is not adequate for neuromuscular diagnostics, whereas myotubes generated by transdifferentiation from an individual's fibroblasts accurately reflect the muscle transcriptome and faithfully reveal disease-causing mutations. Our work confirms that RNA-seq can greatly improve diagnostic yield in genetically unresolved cases of Mendelian disease, defines strengths and challenges of the technology, and demonstrates the suitability of cell models for RNA-based diagnostics. Our data set the stage for development of RNA-seq as a powerful clinical diagnostic tool that can be applied to the large population of individuals with undiagnosed, rare diseases and provide a framework for establishing minimally invasive strategies for doing so.