A significant proportion of mammalian genomes corresponds to genes that transcribe long non-coding RNAs (lncRNAs). Throughout the last decade, the number of studies concerning the roles played by ...lncRNAs in different biological processes has increased considerably. This intense interest in lncRNAs has produced a major shift in our understanding of gene and genome regulation and structure. It became apparent that lncRNAs regulate gene expression through several mechanisms. These RNAs function as transcriptional or post-transcriptional regulators through binding to histone-modifying complexes, to DNA, to transcription factors and other DNA binding proteins, to RNA polymerase II, to mRNA, or through the modulation of microRNA or enzyme function. Often, the lncRNA transcription itself rather than the lncRNA product appears to be regulatory. In this review, we highlight studies identifying lncRNAs in the homeostasis of various cell and tissue types or demonstrating their effects in the expression of protein-coding or other non-coding RNA genes.
Abstract
Osteoarthritis is a prevalent joint disease and a major cause of disability worldwide with no curative therapy. Development of disease-modifying therapies requires a better understanding of ...the molecular mechanisms underpinning disease. A hallmark of osteoarthritis is cartilage degradation. To define molecular events characterizing osteoarthritis at the whole transcriptome level, we performed deep RNA sequencing in paired samples of low- and high-osteoarthritis grade knee cartilage derived from 124 patients undergoing total joint replacement. We detected differential expression between low- and high-osteoarthritis grade articular cartilage for 365 genes and identified a 38-gene signature in osteoarthritis cartilage by replicating our findings in an independent dataset. We also found differential expression for 25 novel long non-coding RNA genes (lncRNAs) and identified potential lncRNA interactions with RNA-binding proteins in osteoarthritis. We assessed alterations in the relative usage of individual gene transcripts and identified differential transcript usage for 82 genes, including ABI3BP, coding for an extracellular matrix protein, AKT1S1, a negative regulator of the mTOR pathway and TPRM4, coding for a transient receptor potential channel. We further assessed genome-wide differential splicing, for the first time in osteoarthritis, and detected differential splicing for 209 genes, which were enriched for extracellular matrix, proteoglycans and integrin surface interactions terms. In the largest study of its kind in osteoarthritis, we find that isoform and splicing changes, in addition to extensive differences in both coding and non-coding sequence expression, are associated with disease and demonstrate a novel layer of genomic complexity to osteoarthritis pathogenesis.
Objective
To identify robustly differentially expressed long noncoding RNAs (lncRNAs) with osteoarthritis (OA) pathophysiology in cartilage and to explore potential target messenger RNA (mRNA) by ...establishing coexpression networks, followed by functional validation.
Methods
RNA sequencing was performed on macroscopically lesioned and preserved OA cartilage from patients who underwent joint replacement surgery due to OA (n = 98). Differential expression analysis was performed on lncRNAs that were annotated in GENCODE and Ensembl databases. To identify potential interactions, correlations were calculated between the identified differentially expressed lncRNAs and the previously reported differentially expressed protein‐coding genes in the same samples. Modulation of chondrocyte lncRNA expression was achieved using locked nucleic acid GapmeRs.
Results
By applying our in‐house pipeline, we identified 5,053 lncRNAs that were robustly expressed, of which 191 were significantly differentially expressed (according to false discovery rate) between lesioned and preserved OA cartilage. Upon integrating mRNA sequencing data, we showed that intergenic and antisense differentially expressed lncRNAs demonstrate high, positive correlations with their respective flanking sense genes. To functionally validate this observation, we selected P3H2‐AS1, which was down‐regulated in primary chondrocytes, resulting in the down‐regulation of P3H2 gene expression levels. As such, we can confirm that P3H2‐AS1 regulates its sense gene P3H2.
Conclusion
By applying an improved detection strategy, robustly differentially expressed lncRNAs in OA cartilage were detected. Integration of these lncRNAs with differential mRNA expression levels in the same samples provided insight into their regulatory networks. Our data indicates that intergenic and antisense lncRNAs play an important role in regulating the pathophysiology of OA.
To uncover the microRNA (miRNA) interactome of the osteoarthritis (OA) pathophysiological process in the cartilage.
We performed RNA sequencing in 130 samples (n=35 and n=30 pairs for messenger RNA ...(mRNA) and miRNA, respectively) on macroscopically preserved and lesioned OA cartilage from the same patient and performed differential expression (DE) analysis of miRNA and mRNAs. To build an OA-specific miRNA interactome, a prioritisation scheme was applied based on inverse Pearson's correlations and inverse DE of miRNAs and mRNAs. Subsequently, these were filtered by those present in predicted (TargetScan/microT-CDS) and/or experimentally validated (miRTarBase/TarBase) public databases. Pathway enrichment analysis was applied to elucidate OA-related pathways likely mediated by miRNA regulatory mechanisms.
We found 142 miRNAs and 2387 mRNAs to be differentially expressed between lesioned and preserved OA articular cartilage. After applying prioritisation towards likely miRNA-mRNA targets, a regulatory network of 62 miRNAs targeting 238 mRNAs was created. Subsequent pathway enrichment analysis of these mRNAs (or genes) elucidated that genes within the 'nervous system development' are likely mediated by miRNA regulatory mechanisms (familywise error=8.4×10
). Herein
encodes neurotrophin-3, which controls survival and differentiation of neurons and which is closely related to the nerve growth factor.
By an integrated approach of miRNA and mRNA sequencing data of OA cartilage, an OA miRNA interactome and related pathways were elucidated. Our functional data demonstrated interacting levels at which miRNA affects expression of genes in the cartilage and exemplified the complexity of functionally validating a network of genes that may be targeted by multiple miRNAs.
Osteoarthritis is a complex degenerative joint disease. Here, we investigate matched genotype and methylation profiles of primary chondrocytes from macroscopically intact (low-grade) and degraded ...(high-grade) osteoarthritis cartilage and from synoviocytes collected from 98 osteoarthritis-affected individuals undergoing knee replacement surgery. We perform an epigenome-wide association study of knee cartilage degeneration and report robustly replicating methylation markers, which reveal an etiologic mechanism linked to the migration of epithelial cells. Using machine learning, we derive methylation models of cartilage degeneration, which we validate with 82% accuracy in independent data. We report a genome-wide methylation quantitative trait locus (mQTL) map of articular cartilage and synovium and identify 18 disease-grade-specific mQTLs in osteoarthritis cartilage. We resolve osteoarthritis GWAS loci through causal inference and colocalization analyses and decipher the epigenetic mechanisms that mediate the effect of genotype on disease risk. Together, our findings provide enhanced insights into epigenetic mechanisms underlying osteoarthritis in primary tissues.
Abstract
Objective
To identify OA subtypes based on cartilage transcriptomic data in cartilage tissue and characterize their underlying pathophysiological processes and/or clinically relevant ...characteristics.
Methods
This study includes n = 66 primary OA patients (41 knees and 25 hips), who underwent a joint replacement surgery, from which macroscopically unaffected (preserved, n = 56) and lesioned (n = 45) OA articular cartilage were collected Research Arthritis and Articular Cartilage (RAAK) study. Unsupervised hierarchical clustering analysis on preserved cartilage transcriptome followed by clinical data integration was performed. Protein–protein interaction (PPI) followed by pathway enrichment analysis were done for genes significant differentially expressed between subgroups with interactions in the PPI network.
Results
Analysis of preserved samples (n = 56) resulted in two OA subtypes with n = 41 (cluster A) and n = 15 (cluster B) patients. The transcriptomic profile of cluster B cartilage, relative to cluster A (DE-AB genes) showed among others a pronounced upregulation of multiple genes involved in chemokine pathways. Nevertheless, upon investigating the OA pathophysiology in cluster B patients as reflected by differentially expressed genes between preserved and lesioned OA cartilage (DE-OA-B genes), the chemokine genes were significantly downregulated with OA pathophysiology. Upon integrating radiographic OA data, we showed that the OA phenotype among cluster B patients, relative to cluster A, may be characterized by higher joint space narrowing (JSN) scores and low osteophyte (OP) scores.
Conclusion
Based on whole-transcriptome profiling, we identified two robust OA subtypes characterized by unique OA, pathophysiological processes in cartilage as well as a clinical phenotype. We advocate that further characterization, confirmation and clinical data integration is a prerequisite to allow for development of treatments towards personalized care with concurrently more effective treatment response.
Summary
Pemphigus foliaceus (PF) is a blistering autoimmune skin disease rare in most of the world but endemic in certain regions of Brazil. PF is characterized by the detachment of epidermal cells ...and the presence of autoantibodies against desmoglein 1. In previous studies, we have shown that genetic polymorphisms and variable expression levels of certain leucocyte receptor complex (LRC) genes were associated with PF. However, the role of the LRC on PF susceptibility remained to be investigated. Here, we analysed 527 tag single nucleotide polymorphisms (SNPs) distributed within the 1·5 Mb LRC. After quality control, a total of 176 SNPs were analysed in 229 patients with PF and 194 controls. Three SNPs were associated with differential susceptibility to PF. The intergenic variant rs465169 odds ratio (OR) = 1·50; P = 0·004 is located in a region that might regulate several immune‐related genes, including VSTM1, LILRB1/2, LAIR1/2, LILRA3/4 and LENG8. The rs35336528 (OR = 3·44; P = 0·009) and rs1865097 (OR = 0·57; P = 0·005) SNPs in LENG8 and FCAR genes, respectively, were also associated with PF. Moreover, we found four haplotypes with SNPs within the KIR3DL2/3, LAIR2 and LILRB1 genes associated with PF (P < 0·05), which corroborate previously reported associations. Thus, our results confirm the importance of the LRC for differential susceptibility to PF and reveal new markers that might influence expression levels of several LRC genes, as well as candidates for further functional studies.
Leucocyte receptor complex single nucleotide polymorphisms associated with pemphigus foliaceus are potentially implicated in gene regulation.
It has been found that the majority of disease-associated genetic variants identified by genome-wide association studies are located outside of protein-coding regions, where they seem to affect ...regions that control transcription (promoters, enhancers) and non-coding RNAs that also can influence gene expression. In this review, we focus on two classes of non-coding RNAs that are currently a major focus of interest: micro-RNAs and long non-coding RNAs. We describe their biogenesis, suggested mechanism of action, and discuss how these non-coding RNAs might be affected by disease-associated genetic alterations. The discovery of these alterations has already contributed to a better understanding of the etiopathology of human diseases and yielded insight into the function of these non-coding RNAs. We also provide an overview of available databases, bioinformatics tools, and high-throughput techniques that can be used to study the mechanism of action of individual non-coding RNAs. This article is part of a Special Issue entitled: From Genome to Function.
•Background on micro- and long non-coding RNAs•Genetic variants affecting ncRNA function are implicated in human disease.•Outline of bioinformatics tools, as well as high-throughput techniques to study ncRNAs
Abstract
Objectives
To present an unbiased approach to identify positional transcript single nucleotide polymorphisms (SNPs) of osteoarthritis (OA) risk loci by allelic expression imbalance (AEI) ...analyses using RNA sequencing of articular cartilage and subchondral bone from OA patients.
Methods
RNA sequencing from 65 articular cartilage and 24 subchondral bone from OA patients was used for AEI analysis. AEI was determined for all genes present in the 100 regions reported by the genome-wide association studies (GWAS) catalog that were also expressed in cartilage or bone. The count fraction of the alternative allele (φ) was calculated for each heterozygous individual with the risk SNP or with the SNP in linkage disequilibrium (LD) with it (r2 > 0.6). Furthermore, a meta-analysis was performed to generate a meta-φ (null hypothesis median φ = 0.49) and P-value for each SNP.
Results
We identified 30 transcript SNPs (28 in cartilage and two in subchondral bone) subject to AEI in 29 genes. Notably, 10 transcript SNPs were located in genes not previously reported in the GWAS catalog, including two long intergenic non-coding RNAs (lincRNAs), MALAT1 (meta-φ = 0.54, FDR = 1.7×10−4) and ILF3-DT (meta-φ = 0.6, FDR = 1.75×10−5). Moreover, 12 drugs were interacting with seven genes displaying AEI, of which seven drugs have been already approved.
Conclusions
By prioritizing proxy transcript SNPs that mark AEI in cartilage and/or subchondral bone at loci harbouring GWAS signals, we present an unbiased approach to identify the most likely functional OA risk-SNP and gene. We identified 10 new potential OA risk genes ready for further translation towards underlying biological mechanisms.
Celiac disease (CeD) is triggered by gluten and results in inflammation and villous atrophy of the small intestine. We aimed to explore the role of miRNA-mediated deregulation of the transcriptome in ...CeD. Duodenal biopsies of CeD patients (
= 33) and control subjects (
= 10) were available for miRNA-sequencing, with RNA-sequencing also available for controls (
= 5) and CeD (
= 6). Differential expression analysis was performed to select CeD-associated miRNAs and genes. MiRNA‒target transcript pairs selected from public databases that also displayed a strong negative expression correlation in the current dataset (R < -0.7) were used to construct a CeD miRNA‒target transcript interaction network. The network includes 2030 miRNA‒target transcript interactions, including 423 experimentally validated pairs. Pathway analysis found that interactions are involved in immune-related pathways (e.g., interferon signaling) and metabolic pathways (e.g., lipid metabolism). The network includes 13 genes previously prioritized to be causally deregulated by CeD-associated genomic variants, including STAT1. CeD-associated miRNAs might play a role in promoting inflammation and decreasing lipid metabolism in the small intestine, thereby contributing unbalanced cell turnover in the intestinal crypt. Some CeD-associated miRNAs deregulate genes that are also affected by genomic CeD-risk variants, adding an additional layer of complexity to the deregulated transcriptome in CeD.