Use of RNA sequencing to evaluate rheumatic disease patients

Transcriptomic analysis of the Non-Obstructive Azoospermia (NOA) to address gene expression regulation in human testis

There is a need to understand the molecular basis of testes under Non-Obstructive Azoospermia (NOA), a state of failed spermatogenesis. There has been a lack of attention to the transcriptome at the level of alternatively spliced mRNAs (iso-mRNAs) and the mechanism of gene expression regulation. Hence, we aimed to establish a reliable iso-mRNA profile of NOA-testes, and explore molecular mechanisms - especially those related to gene expression regulation. We sequenced mRNAs from testicular samples of donors with complete spermatogenesis (control samples) and a failure of spermatogenesis (NOA samples). We identified differentially expressed genes and their iso-mRNAs via standard NGS data analyses. We then listed these iso-mRNAs hierarchically based on the extent of consistency of differential quantities across samples and groups, and validated the lists via RT-qPCRs (for 80 iso-mRNAs). In addition, we performed extensive bioinformatic analysis of the splicing features, domains, interactions, and functions of differentially expressed genes and iso-mRNAs. Many top-ranking down-regulated genes and iso-mRNAs, i.e., those down-regulated more consistently across the NOA samples, are associated with mitosis, replication, meiosis, cilium, RNA regulation, and post-translational modifications such as ubiquitination and phosphorylation. Most down-regulated iso-mRNAs correspond to full-length proteins that include all expected domains. The predominance of alternative promoters and termination sites in these iso-mRNAs indicate their gene expression regulation via promoters and UTRs. We compiled a new, comprehensive list of human transcription factors (TFs) and used it to identify TF-'TF gene' interactions with potential significance in down-regulating genes under the NOA condition. The results indicate that RAD51 suppression by HSF4 prevents SP1-activation, and SP1, in turn, could regulate multiple TF genes. This potential regulatory axis and other TF interactions identified in this study could explain the down-regulation of multiple genes in NOA-testes. Such molecular interactions may also have key regulatory roles during normal human spermatogenesis.

Possible Roles of tRNA Fragments, as New Regulatory ncRNAs, in the Pathogenesis of Rheumatoid Arthritis

Understanding the pathophysiology of rheumatoid arthritis (RA) has led to the successful development of molecule-targeted drugs for the treatment of RA. However, some RA patients are refractory to these treatments, suggesting that the pathological mechanism of the disease is not entirely understood. Genome and transcriptome analysis is essential for understanding the unknown pathophysiology of human diseases. Rapid and more comprehensive gene analysis technologies have revealed notable changes in the expression of coding RNA and non-coding RNA in RA patients. This review focuses on the current state of non-coding RNA research in relation to RA, especially on tRNA fragments. Interestingly, it has been found that tRNA fragments repress translation and are antiapoptotic. The association between tRNA fragments and various diseases has been studied, and this article reviews the possible role of tRNA fragments in RA.

Identification of differentially expressed genes and pathways crosstalk analysis in Rheumatoid and Osteoarthritis using next-generation sequencing and protein-protein networks

Osteoarthritis occurs when protective cartilage of bones worn out. Similarlty, cartilage damage occurs mainly in the pannus cartilage in rheumatoid arthritis. It is a potentially debilitating condition, affecting women two to three times more often than men. The cause and prognosis of rheumatoid and osteoarthritis are still poorly known. However, advances in the study of disease pathogenesis have encouraged the creation of new therapeutics with improved outcomes. The purpose of this study is to investigate the differentially expressed genes potentially involved in dysregulated rheumatoid arthritis (RA) and their association to other types of arthritis, including osteoarthritis (OA). Complete RNAs were isolated for RNA expression profiling using next-generation sequencing from human primary cultured normal and RA chondrocytes. From RNA sequencing results 250 differentially expressed genes were identified using bioinformatics analysis, of which 32 were found to be significantly playing role in RA pathogenesis and its associated diseases. Molecular ontologies of the identified genes showed they are connected to Innate immune response, Protein phosphorylation, Transcription initiation from RNA polymerase II promoter, Immune response, Neoplasms of bones, as well as osteorthritis, and Rheumatoid arthritis. Among the identified genes, TRAF1, TRAF2, BAMP, STX11, MEOX2, AES, REL, FHL3, PNMA1, SGTA, LZTS2, SIAH2, PNMA1, and TFCP2 were found to be highly enriched in the protein-protein interaction network. The significant cross talks were found in Hypertrophic cardiomyopathy, Small cell lung cancer, Proteasome, p53 signaling pathway, Arrhythmogenic right ventricular cardiomyopathy, Small cell lung cancer, SNARE interactions in vesicular transport, RIG-I-like receptor signaling pathway, and Hypertrophic cardiomyopathy pathways. The results offer new opportunities for target gene control in RA and OA cartilage destruction.

CD14 Involvement in Third-degree Skin Burn-induced Myocardial Injury via the MAPK Signaling Pathway

This study investigated the potential genes and related pathways in burn-induced myocardial injury. Rat myocardial injury induced by third-degree burn and the histopathological structures, apoptosis, and cardiac injury markers were then identified using hematoxylin & eosin staining, terminal deoxynucleotidyl transferase-mediated dUTP nick end labeling staining, and enzyme-linked immunosorbent assay. Next, differentially expressed mRNAs were screened through next-generation sequencing (NGS), followed by functional annotation and key gene validation through quantitative reverse transcription-polymerase chain reaction. Subsequently, CD14 was screened out, and small interfering RNAs against CD14 were transfected to H9C2 cells to further verify the role of CD14 in burn-induced injury. The results showed that third-degree burn could markedly damage the structure of myocardial tissue, induce the apoptosis of myocardial cells, and increase the levels of myocardial injury-related markers, suggesting that burns could induce myocardial injury in rats. Besides, NGS data discovered that third-degree burn could result in 416 differentially upregulated mRNAs and 285 differentially downregulated mRNAs in myocardial tissue. It was also disclosed that differentially expressed mRNAs were mainly enriched in the phosphatidylinositol 3-kinase/Akt, mitogen-activated protein kinase (MAPK), and tumor necrosis factor signaling pathways. Furthermore, cell viability was significantly decreased in H9C2 cells treated with 10% rat burn serum. CD14 was significantly differentially expressed and screened out for further studies. Treatment with burn serum can significantly upregulate the phosphorylation level of extracellular signal-regulated kinase, p38, and c-Jun N-terminal kinase and the expression of cleaved caspase-3 and downregulate the expression of Bcl2 when compared with those in negative control of small interfering RNA transfected H9C2 cells, whereas interfering with CD14 expression reversed the effects of burn serum. The study demonstrated that burn serum treatment could activate the MAPK signaling pathway to promote cell apoptosis, and it can be reversed by interfering with the expression of CD14.

Nutrition and Rheumatoid Arthritis in the 'Omics' Era

Modern high-throughput 'omics' science tools (including genomics, transcriptomics, proteomics, metabolomics and microbiomics) are currently being applied to nutritional sciences to unravel the fundamental processes of health effects ascribed to particular nutrients in humans and to contribute to more precise nutritional advice. Diet and food components are key environmental factors that interact with the genome, transcriptome, proteome, metabolome and the microbiota, and this life-long interplay defines health and diseases state of the individual. Rheumatoid arthritis (RA) is a chronic autoimmune disease featured by a systemic immune-inflammatory response, in genetically susceptible individuals exposed to environmental triggers, including diet. In recent years increasing evidences suggested that nutritional factors and gut microbiome have a central role in RA risk and progression. The aim of this review is to summarize the main and most recent applications of 'omics' technologies in human nutrition and in RA research, examining the possible influences of some nutrients and nutritional patterns on RA pathogenesis, following a nutrigenomics approach. The opportunities and challenges of novel 'omics technologies' in the exploration of new avenues in RA and nutritional research to prevent and manage RA will be also discussed.

Peripheral blood mononuclear cells extracts VEGF protein levels and VEGF mRNA: Associations with inflammatory molecules in a healthy population

Background

Vascular endothelial growth factor (VEGF) is a signal protein, implicated in various physiological and pathophysiological processes together with other common inflammatory biomarkers. However, their associations have not yet been fully elucidated. In the present study, we investigated associations between VEGF and four specific VEGF mRNA isoforms with levels of 11 inflammation molecules, derived from peripheral blood mononuclear cells (PBMCs) extracts.

Methods

Healthy participants from the STANISLAS Family Study (n = 285) were included. Levels of VEGF (four mRNA isoforms and protein levels) and inflammatory molecules (IL-1α, IL-1β, IL-2, IL-4, IL-6, IL-8, IL-10, INF-γ, TNF-α, MCP-1, EGF) were measured in PBMCs extracts. Multiple regression analyses were performed, adjusted for age and gender.

Results

The analyses revealed significant associations between VEGF protein levels and levels of IL-4 (β = 0.028, P = 0.013), MCP-1 (β = 0.015, P<0.0001) and EGF (β = 0.017, P<0.0001). Furthermore, mRNA isoform VEGF₁₆₅ was associated with MCP-1 and IL-1α (P = 0.002 and P = 0.008, respectively); and mRNA isoform VEGF₁₈₉ was associated with IL-4 and IL-6 (P = 0.019 and P = 0.034, respectively).

Conclusions

To our knowledge, the present study represents the first investigation that successfully demonstrates links between VEGF protein levels and inflammatory molecules levels derived from PBMCs extracts and identifies associations between specific VEGF mRNA isoforms and inflammatory molecules.

Impact

These findings provide novel insights that may assist in the development of new tissue and mRNA isoform specific measurements of VEGF levels, which may positively contribute to predicting the risk of common complex diseases and response of currently used anti-VEGF agents, and developing of novel targeted therapies for VEGF-related pathophysiology.

Comparative analysis of differential gene expression tools for RNA sequencing time course data

RNA sequencing (RNA-seq) has become a standard procedure to investigate transcriptional changes between conditions and is routinely used in research and clinics. While standard differential expression (DE) analysis between two conditions has been extensively studied, and improved over the past decades, RNA-seq time course (TC) DE analysis algorithms are still in their early stages. In this study, we compare, for the first time, existing TC RNA-seq tools on an extensive simulation data set and validated the best performing tools on published data. Surprisingly, TC tools were outperformed by the classical pairwise comparison approach on short time series (<8 time points) in terms of overall performance and robustness to noise, mostly because of high number of false positives, with the exception of ImpulseDE2. Overlapping of candidate lists between tools improved this shortcoming, as the majority of false-positive, but not true-positive, candidates were unique for each method. On longer time series, pairwise approach was less efficient on the overall performance compared with splineTC and maSigPro, which did not identify any false-positive candidate.

Transcriptome analysis of peripheral blood from patients with rheumatoid arthritis: a systematic review

In the era of precision medicine, transcriptome analysis of whole gene expression is an essential technology. While DNA microarray has a limited dynamic range and a problem of background hybridization, RNA sequencing (RNA-seq) has a broader dynamic range and a lower background signal that increase the sensitivity and reproducibility. While transcriptome analyses in rheumatoid arthritis (RA) have generally focused on whole peripheral blood mononuclear cells (PBMC), analyses of detailed cell subsets have an increased need for understanding the pathophysiology of disease because the involvement of CD4⁺ T cells in the pathogenesis of RA has been established. Transcriptome analysis of detailed CD4⁺ T cell subsets or neutrophils shed new light on the pathophysiology of RA. There are several analyses about the effect of biological treatment. Many studies report the association between type I interferon signature gene expression and response to therapy.

Systematic Analysis of Differential Expression Profile in Rheumatoid Arthritis Chondrocytes Using Next-Generation Sequencing and Bioinformatics Approaches

Cartilage destruction in rheumatoid arthritis (RA) occurs primarily in the pannus-cartilage interface. The close contact of the synovium-cartilage interface implicates crosstalk between synovial fibroblasts and chondrocytes. The aim of this study is to explore the differentially expressed genes and novel microRNA regulations potentially implicated in the dysregulated cartilage homeostasis in joint destruction of RA. Total RNAs were extracted from human primary cultured normal and RA chondrocytes for RNA and small RNA expression profiling using next-generation sequencing. Using systematic bioinformatics analyses, we identified 463 differentially expressed genes in RA chondrocytes were enriched in biological functions related to altered cell cycle process, inflammatory response and hypoxic stimulation. Moreover, fibroblast growth factor 9 (FGF9), kynureninase (KYNU), and regulator of cell cycle (RGCC) were among the top dysregulated genes identified to be potentially affected in the RA joint microenvironment, having similar expression patterns observed in arrays of clinical RA synovial tissues from the Gene Expression Omnibus database. Additionally, among the 31 differentially expressed microRNAs and 10 candidate genes with potential microRNA-mRNA interactions in RA chondrocytes, the novel miR-140-3p-FGF9 interaction was validated in different microRNA prediction databases, and proposed to participate in the pathogenesis of joint destruction through dysregulated cell growth in RA. The findings provide new perspectives for target genes in the management of cartilage destruction in RA.

Integrative Approaches to Understanding the Pathogenic Role of Genetic Variation in Rheumatic Diseases

The use of high-throughput omics may help to understand the contribution of genetic variants to the pathogenesis of rheumatic diseases. We discuss the concept of missing heritability: that genetic variants do not explain the heritability of rheumatoid arthritis and related rheumatologic conditions. In addition to an overview of how integrative data analysis can lead to novel insights into mechanisms of rheumatic diseases, we describe statistical approaches to prioritizing genetic variants for future functional analyses. We illustrate how analyses of large datasets provide hope for improved approaches to the diagnosis, treatment, and prevention of rheumatic diseases.

Screening for key genes and transcription factors in ankylosing spondylitis by RNA-Seq

Ankylosing spondylitis (AS) is a chronic inflammatory arthritis and autoimmune disease, the etiology and pathogenesis of which remain largely unknown. In the present study, blood samples were harvested from patients with AS and from healthy volunteers as a normal control (NC) for RNA-sequencing. Differentially expressed genes (DEGs) in the AS group compared with the NC group were identified, and gene ontology (GO) term and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were subsequently performed. Protein-protein interaction (PPI) network and AS-specific transcriptional regulatory network construction was performed for the DEGs. A total of 503 DEGs, including 338 upregulated and 165 downregulated DEGs, were identified in patients with AS compared with the NC group. Three upregulated DEGs identified, interferon-induced protein with tetratricopeptide repeats (IFIT)1, IFIT3 and radical S-adenosyl methionine domain containing (RSAD)2, are interferon (IFN)-stimulated genes that serve a role in the IFN signaling pathway. The most significantly enriched GO term was response to other organisms. Osteoclast differentiation was a significantly enriched pathway for eight DEGs [High affinity immunoglobulin gamma Fc receptor (FCGR)1A, FCGR2B, four and a half LIM domains 2, integrin β3, signal transducer and activator of transcription 2 (STAT2), suppressor of cytokine signaling 3 (SOCS3), leukocyte immunoglobulin like receptor (LILR)A4 and LILRA6]. The six hub genes in the PPI network constructed were interferon-stimulated gene 15, heat shock protein β1, microtubule-associated proteins 1A/1B light chain 3A, IFIT1, IFIT3 and SOCS3. POU domain class 2 transcription factor 1 (1-Oct) and ecotropic virus integration site-1 (Evi-1) were identified as two important transcription factors (TFs) in AS according to the AS-specific transcriptional regulatory network constructed. In addition, IFIT1 and IFIT3 were identified as targets of 1-Oct. The results of the present study indicate that osteoclast differentiation, the IFN signaling pathway and genes associated with these two signaling pathways, particularly FCGR2B, STAT2, SOCS3, IFIT1 and IFIT3, may serve a role in AS. In addition, Evi-1 and 1-Oct may be two important TFs associated with AS. These results may provide a basis for elucidating the underlying mechanisms of and developing novel treatments for AS.

Transcriptome analysis of ankylosing spondylitis patients before and after TNF-α inhibitor therapy reveals the pathways affected

Tumor necrosis factor-α (TNF-α) inhibitors are highly effective in suppressing inflammation in ankylosing spondylitis (AS) patients, and operate by suppression of TFN-α and downstream immunological pathways. To determine the mechanisms of action of TNF-α inhibitors in AS patients, we used transcriptomic and bioinformatic approaches on peripheral blood mononuclear cells from AS patients pre and post treatment. We found 656 differentially expressed genes, including the genome-wide significant AS-associated genes, IL6R, NOTCH1, IL10, CXCR2 and TNFRSF1A. A distinctive gene expression profile was found between male and female patients, mainly because of sex chromosome-linked genes and interleukin 17 receptor C, potentially accounting for the differences in clinical manifestation and treatment response between the genders. In addition to immune and inflammation regulatory pathways, like intestinal immune network for IgA production, cytokine-cytokine receptor interaction, Ras signaling pathway, allograft rejection and hematopoietic cell lineage, KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway analyses revealed that infection-associated pathways (influenza A and toxoplasmosis) and metabolism-associated pathways were involved in response to TNF-α inhibitor treatment, providing insight into the mechanism of TNF-α inhibitors.

RNA sequencing to determine the contribution of kinase receptor transactivation to G protein coupled receptor signalling in vascular smooth muscle cells

G protein coupled receptor (GPCR) signalling covers three major mechanisms. GPCR agonist engagement allows for the G proteins to bind to the receptor leading to a classical downstream signalling cascade. The second mechanism is via the utilization of the β-arrestin signalling molecule and thirdly via transactivation dependent signalling. GPCRs can transactivate protein tyrosine kinase receptors (PTKR) to activate respective downstream signalling intermediates. In the past decade GPCR transactivation dependent signalling was expanded to show transactivation of serine/threonine kinase receptors (S/TKR). Kinase receptor transactivation enormously broadens the GPCR signalling paradigm. This work utilizes next generation RNA-sequencing to study the contribution of transactivation dependent signalling to total protease activated receptor (PAR)-1 signalling. Transactivation, assessed as gene expression, accounted for 50 percent of the total genes regulated by thrombin acting through PAR-1 in human coronary artery smooth muscle cells. GPCR transactivation of PTKRs is approximately equally important as the transactivation of the S/TKR with 209 and 177 genes regulated respectively, via either signalling pathway. This work shows that genome wide studies can provide powerful insights into GPCR mediated signalling pathways.

Diagnostic and prognostic tests in systemic lupus erythematosus

Systemic lupus erythematosus (SLE) is a chronic autoimmune inflammatory disease characterized by autoantibodies directed against numerous self-nuclear antigens. Because of the heterogeneous nature of lupus, it has been challenging to identify markers that are sensitive and specific enough for its diagnosis and monitoring. However, with the sequencing of the human genome, rapid development of high-throughput approaches has allowed for a better understanding of the etiopathogenesis of complex diseases, including SLE. Here we present a review of the latest advancements in biomarker discovery during the "omics" era, using these novel technologies, for assisting in the diagnosis and prognosis of patients with SLE.

Personalized medicine in rheumatology

In the era of the 21st century, rheumatoid arthritis (RA) is still poorly characterized. Rheumatoid arthritis is a common but heterogeneous disease, not only in the course and clinical symptoms, but also in the clinical response to treatment. Now it is known that early, correct diagnosis and starting treatment with disease-modifying drugs (DMARDs), of which methotrexate (MTX) remains the gold standard in the treatment of RA, is crucial in order to prevent joint destruction, functional disability and an unfavourable disease outcome. Early diagnosis of rheumatoid arthritis is significant in so much as the primary treatment can be started better. Pharmacogenetic and pharmacogenomic studies, which help determine the genetic profile of individual patients, may bring us closer to personalized medicine. Further studies on RA should allow for the identification of disease-specific genes at the stage when their tolerance by the organism is still preserved (before auto-aggression develops).

Limits of Peripheral Blood Mononuclear Cells for Gene Expression-Based Biomarkers in Juvenile Idiopathic Arthritis

Juvenile Idiopathic Arthritis (JIA) is one of the most common chronic disease conditions affecting children in the USA. As with many rheumatic diseases, there is growing interest in using genomic technologies to develop biomarkers for either diagnosis or to guide treatment ("personalized medicine"). Here, we explore the use of gene expression patterns in peripheral blood mononuclear cells (PBMC) as a first step approach to developing such biomarkers. Although PBMC carry many theoretical advantages for translational research, we have found that sample heterogeneity makes RNASeq on PBMC unsuitable as a first-step method for screening biomarker candidates in JIA. RNASeq studies of homogeneous cell populations are more likely to be useful and informative.