- Research article
- Open Access
Integrated analysis of mRNA, microRNA and protein in systemic lupus erythematosus-specific induced pluripotent stem cells from urine
BMC Genomics volume 17, Article number: 488 (2016)
In clinical practice, it is difficult to monitor the repeating relapse in patients who have been suffering from systemic lupus erythematosus (SLE). The underlying etiology remains largely unknown.
Aiming to understand the pathogenesis of SLE, a detailed study was conducted. Renal tubular cells–derived iPSCs were successfully obtained from the urine of SLE patients and healthy controls. With the purpose to identify simultaneous expression profiling of microRNA, mRNA and protein, Illumina HiSeq™ 2000 System and iTRAQ-coupled 2D LC-MS/MS analysis were utilized in systemic lupus erythematosus-specific induced pluripotent stem cells (SLE-iPSCs) and normal control-iPSCs (NC-iPSCs). The integration of multiple profiling datasets was realized since it could facilitate the identification of non-seed miRNA targets, as well as differentially expressed mRNAs and proteins.
For this study, profiling datasets of 1099 differentially expressed mRNAs, 223 differentially expressed microRNAs and 94 differentially expressed proteins were integrated. In order to investigate the influence of miRNA on the processes of regulating mRNAs and proteins’ levels, potential targets of differentially expressed mRNAs and proteins were predicted using miRanda, TargetScan and Pictar. Multiple profiling datasets were integrated to facilitate the identification of miRNA targets, as well as differentially expressed mRNAs and proteins. Through gene ontology (GO) analysis of differentially expressed mRNAs and proteins, biological processes that drive proliferation were identified, such as mRNA processing and translation. Western blot and Q-PCR confirmed AK4 protein and mRNA up-regulation. The findings also showed that TAGLN’s protein and mRNA level were down-regulated in SLE-iPSCs, both miR-371a-5p and let-7a-5p in SLE-iPSC were down-regulated and verified using Q-PCR. The up-regulation of AK4 involved in nucleotide biosynthesis suggested a general acceleration of anabolic metabolism induced by down-regulated miR-371a-5p, which might contribute to SLE.
Based on high throughput analysis, integrated miRNA, mRNA, and protein expression data were generated. Differentially expressed dates were also adopted in conjunction with in-silico tools to identify potential candidates for SLE-iPSCs. Representative miRNA, mRNA and proteins were verified. It was also expected that the knowledge gained from this study can be applied to assess the usefulness of pathogenesis and novel biomarker candidates of SLE, which may develop a new way for SLE diagnosis.
Systemic lupus erythematosus (SLE) is the prototype of complex autoimmune disease characterized by the production of autoantibodies which results in widespread immunologic abnormalities and immune complex formation . The patients can present variable manifestations and the nature courses are alternately remissions and relapses. Till now, although a lot of related researches have been undertaken , SLE patients have no effective cures, whose treatments are often based upon long-term broad-spectrum immune suppressive regimes in the current therapeutic management. It becomes a major public health problem.
The exact mechanism involved in SLE is needed to understand. a lot of work has been done in searching for biomarkers, from the aspects of DNA, mRNA and protein, expecting to illustrate the mechanism of SLE and find ideal biomarkers for diagnosis and as a precaution to SLE [3–5]. However, the underlying mechanism and pathogenesis of SLE are still far away from understanding. Novel methods should be looked into in this area.
As is known to all, the possibility of reprograming somatic cells to induced pluripotent stem cells (iPSCs) offers an opportunity to generate pluripotent patient-specific cell lines, which can be beneficial for studying pathogenesis of model human diseases . Also, these iPSCs lines are powerful tools for recapitulating disease conditions and thus better understanding the underlying mechanisms and pathogenesis of specific diseases . iPSCs can be derived from immune cells as equally as they can differentiate into specific immune cell types for modeling diseases or clinical immunotherapy . So far, generations of iPSCs from urine, fibroblasts and keratinocytes of disease patients have been reported [7, 8]. On this basis, renal tubular cells-derived iPSCs were successfully acquired from urine of SLE patients to study SLE pathogenesis .
Currently, extensive researches, including those involving system-wide genomic and transcriptomic approaches, have been conducted to characterize iPSCs . However, molecular mechanisms remain insufficiently understood. At the molecular level, iPSCs are more variable than embryonic stem cells (ESCs) . It is imperative that the compendium of differences described between iPSCs and ESCs be considered as further evidence for the fact that the reprogramming process requires a wide variety of molecular changes. Besides, cells can distinguish bona fide iPSCs from partial reprogrammed cells during an earlier stage . Therefore, systematic deciphering of genetic or epigenetic alternations would help to identify hotspots in iPSCs with consistent genetic background .
Briefly, renal tubular cells–derived iPSCs were successfully obtained from urine of SLE patients and healthy controls to understand the pathogenesis of SLE. For the integration of multiple profiling datasets, simultaneous expression profiling of microRNA, mRNA and protein in SLE-iPSCs and NC-iPSCs were identified with Illumina HiSeq™ 2000 System and iTRAQ-coupled 2D LC-MS/MS analysis.
mRNA, protein, miRNA expression profiles
In order to comprehensively reveal mRNA, microRNA and protein interactions, we obtained high quality, complete information and estimated the expression levels of mRNA, miRNA and protein between the SLE-iPSCs and control-iPSCs (Fig. 1).
For the mRNA expression profiling, P-value ≤ 0.05 and FDR ≤ 0.001 were set as threshold value, 4,254 genes were detected to have at least two-fold differences between SLE-iPSCs and control-iPSCs, The genes numbers of 2,856 and 1,398 respectively represent the higher and lower abundances of more than two fold compared with control-iPSCs. Among 4,254 genes, 1,099 differentially expressed ones built an interaction network, with 744 genes up-regulated and 355 down-regulated (Additional file 1: Table S1).
Compared with the control group, it was found in the microRNA libraries that 223 miRNAs were expressed at the level of significant difference, with 126 up-regulated miRNAs and 97 down-regulated ones in SLE-iPSCs (Additional file 2: Table S2).
In the protein expression profiling, the identification and quantification of 2,305 proteins in SLE-iPSCs and control-iPSCs were performed with iTRAQ technique. Given a confidence level of 95 %, 1.5-fold was set as the cut-off criterion for up-regulated and down-regulated proteins. Under these criteria, a total of 207 proteins were classified as differentially expressed ones between SLE-iPSCs and control-iPSCs, with 55 showing increased abundance and the remaining 152 decreased abundance (Additional file 3, Additional file 4: Table S3,4).
MicroRNA-target genes regulation network
With integrated mRNA and microRNA transcriptome, the expression change trend of miRNA target mRNAs were focused on. Then, a systematically investigation was conducted on the potential functional correlations among microRNA target mRNAs. Cytoscape software was used to construct the regulation network of differentially expressed mRNAs. Their microRNA-target mRNAs were predicted by using miRanda, TargetScan and Pictar. If the target mRNAs were successfully predicted at least with two kinds of software at the same time, they would be considered to be reliable. MicroRNA-target mRNAs regulation network was shown in Figs. 2, 3. It can be seen that microRNA and target genes were mutually cross-regulated.
A functional analysis was conducted based on Gene Ontology. Some gene ontologies of biological processes, molecular functions and cellular components were selected, which were enriched by the features of the transcriptomic and proteomic datasets with p value < 0.05. In the study, GO enrichment analysis was performed through functional annotation clustering of microRNA and target genes, respectively. Totally, 11 biological processes (BP), 4 molecular functions (MF) and 14 cellular components (CC) GO terms were enriched in target genes. 10 BP, 4 MF and 13 CC GO terms in microRNA are shown in Fig. 4 a–b, respectively. BP includes translation, mRNA processing, nucleocytoplasmic transportation, vesicle-mediated transportation, cell motility and so on. MF involves structural constituent of ribosome, translation factor activity, nucleic acid binding, rRNA binding, mRNA binding. CC contains nucleolus nucleus, ribosome, Golgi apparatus, chromosome, endoplasmic reticulum, cytoplasmic membrane-bounded vesicle, mitochondrion, cytoplasmic membrane-bounded vesicle and so on. According to the biological process, the microRNA and target genes are mostly involved in BP GO terms, structural constituent of ribosome of MF GO term and nucleus of CC GO terms.
Generally, microRNAs are classified as a class of small non-coding RNAs that bind complementary sequences in target mRNAs to specifically regulate gene expression through either mRNA degradation or translational inhibition . If a target gene is down-regulated, it suggests that effective activity of miRNA is enhanced under the treatment, while an up-regulation of a target gene indicates a decreased activity of the corresponding miRNAs. Consequently, a miRNA-mRNA interaction pair means anti-regulation of a miRNA and a corresponding mRNA [14, 15].
MicroRNA-target protein regulation network
In this study, a comparative proteome survey was performed on the SLE-iPSC and control-iPSC using iTRAQ technique. The identification and quantification of differentially expressed proteins were realized. In order to investigate the potential functional correlations among microRNA target proteins, an integrated systematic analysis was made on microRNA and protein data. In addition, Cytoscape software was applied to construct the regulation network of microRNA and target protein. Given a 95 % confidence level, the interaction network of 94 proteins in 207 differentially expressed proteins were classified (Fig. 5); 37 MicroRNAs which regulated 49 target proteins were predicted. MicroRNA-target proteins regulation network is shown in Fig. 5, which manifests that microRNA- target protein were mutually cross-regulated.
GO enrichment analyses of microRNA and target proteins were performed, respectively. Totally, 4 BP, 3 MF and 10 CC GO terms were enriched in target proteins. 4 BP, 1MF and 9 CC GO terms in microRNA are respectively shown in Fig. 6a–b. BP includes translation, nucleocytoplasmic transport, embryo development and mRNA processing; MF contains structural constituents of ribosome; CC involves nucleolusnucleus, ribosome, endoplasmic reticulum, mitochondrion, cytoskeleton and so on. According to biological process, the microRNA-target proteins are mostly involved in the translation of BP GO terms, structural constituent of ribosome of MF GO term and nucleus of CC GO terms.
mRNA-protein regulation network
For the purpose of identifying the correlations between mRNA and protein expressions, the expression change trend of transcriptomic and proteomic profiles were focused on. 1099 differentially expressed mRNAs and proteins were found to be correlated in the mRNA-protein library (Fig. 7).
A functional analysis was conducted through the KEGG pathway. Given a cut-off criterion of Q-value < 0.05, 1109 differentially-expressed mRNAs were involved in 235 pathways, and 99 differentially expressed proteins in 65 pathways (Additional file 5, Additional file 6: Table S5,6). According to the enrichment analysis, a few important pathways were significantly enriched in response to SLE-iPSCs. It was quite evident that mRNAs and proteins were both involved in ribosome, spliceosome, RNA transport, lysine degradation and so on.
Based on the central dogma, it was generally assumed that there was a direct correspondence between mRNA transcripts and generated protein expressions. In recent years, various studies have demonstrated the role of RNA-binding proteins (RBPs) . On the basis of mRNA-protein correlation, it was found that 11 proteins and 11 mRNAs were differentially expressed consistently (3 up-regulated and 8 down-regulated, Fig. 8). To further validate the expression data, 2 differentially expressed mRNAs and proteins were randomly chosen and verified with Q-PCR and western blotting (Figs. 9, 10).
The common KEGG pathway of differentially expressed proteins and mRNAs were analyzed using GenMAPP v2.1. KEGG enrichment analyses showed that proteins and mRNA significantly participated in ECM-receptor interaction, fatty acid metabolism, ribosome, spliceosome, RNA transport and so on(Additional file 5, Additional file 6: Table S5,6).
mRNA-protein-microRNA regulation network
In order to integrate profiling datasets of differentially expressed mRNAs, microRNAs and Proteins, the mRNA-protein-microRNA regulation network was constructed with 39 down-regulated target proteins, 42 down-regulated mRNAs and 21 up-regulated miRNAs (Fig. 11). The result indicated coherent miRNAs-target proteins and miRNAs-target mRNAs interaction pairs. This demonstrates that a single miRNA can target multiple proteins and mRNAs, while a single protein or mRNA can be targeted by multiple miRNAs, which can cooperatively repress a range of targets.
After Takahashi and Yamanaka demonstrate that murine fibroblasts can be reprogrammed into iPSCs in 2006 . The application of powerful technique to humans opened the possibility to reprogram fibroblasts and/or various specific immune cells isolated not only from healthy people but also from patients suffering of a genetic disease [6, 18]. What is more, the advances in stem cell technology of iPSCs offer the opportunity to generate and directly analyze human autoimmune disease in healthy and disease states . Despite the fact that human iPSC-derived cells remained controversial to the present day, the system has attested successful in vitro replication of the main cellular characteristics already known to be associated with the modeled disease . Araki et al.  verified that no pronounced differences in immunogenicity existed between the teratomas formed by seven iPSCs and five embryonic stem cell lines. Based on the research results, iPSCs can be derived from syngeneic autologous cells as equally as they can differentiate into specific immune cell types for modeling diseases or clinical immunotherapy .
Network modeling, analysis and subsequent target identification approach appear more promising in characterizing diseases . Biological networks comprise of interactions between genes and proteins that co-ordinate in the regulation of cellular processes . Computational and systems biological techniques can be employed to study the cellular level networks that mediate physiological system dynamics . In this study, an integrative methodology was employed to combine the data obtained through multiple expression profiling methods, genomic sequence and bioinformatics analysis of in-silico target miRNA, so as to investigate the impacts of miRNA, mRNA and protein expression on SLE-iPSCs.
Differentially expressed proteins, mRNAs and miRNAs were systematically investigated. Computational methods were adopted to predict miRNA interaction and prioritize potential direct targets. Finally, mRNA-proteomic mapping, miRNA-target mRNA network and miRNA-target protein interactions were achieved. Systematic comparison results demonstrated the common physiological changes in SLE-iPSCs and NC-iPSCs. Nucleus and nucleolus in cellular component, transmembrane transport in biological process and structural constituent of ribosome in molecular function were found to be significant clustered in GO terms. That is to say, systemic lupus erythematosus (SLE) is a systemic non-organ specific autoimmune disease associated with multiple autoantibodies targeting autoantigens from the nucleus . Anti-single and anti-double-stranded DNA antibodies and rheumatoid factor were marker proteins of SLE, which mainly derived from the cell nucleus. Our results showed that differentially expressed proteins, mRNAs and target miRNAs were mainly directed to nuclear constituents consistently. A subset of nuclear constituents is directed to proteins and RNAs of the ribosomes . While a part of auto-antibodies directed against some antigenic site composed of a portion of both ribosomal RNA and ribosomal protein . As a result, differentially expressed proteins, mRNAs and target miRNAs, which are clustered in structural constituent of ribosome, will become markers of SLE in the future.
Transport systems are essential to each living cell, which is vital for all life-endowing processes, such as communication, biosynthesis, reproduction, and cooperative behaviors . Given this, SLE is viewed as expressions of combined transport dysfunction syndromes . The results indicate that most differentially expressed proteins, mRNAs and target miRNAs are classified as transmembrane transport in biological process.
In the light of the central dogma, it was generally assumed that there was a direct correspondence between mRNA transcripts and generated protein expressions. By combining mRNA and protein expression data, it was found that 11 proteins and 11 RNAs were differentially expressed consistently in SLE-iPSCs and NC-iPSCs. Transgelin (TAGLN), also known as smooth muscle protein 22 (SM22), is a ubiquitous 22 kDa protein among smooth muscle tissues of normal adult vertebrates . TAGLN plays a role in cell differentiation by stabilizing the cytoskeleton through actin-binding . SM22/transgelin gene expression were down-regulated during a variety of cell types [30, 31], as well as in several human cancers, including lung, renal and breast cancer [30, 32, 33]. TAGLN has been proposed as a candidate tumor antigen . The study results showed that TAGLN was down-regulated in SLE-iPSCs. Western blot and RT-PCR confirmed the down-regulation of both TAGLN protein and mRNA. The former belongs to transcriptional down-regulation. The mechanisms controlling the regulations of other proteins require additional study in SLE-iPSCs.
Adenylate kinase (AK) is an ATP-ADP phosphotransferase that catalyzes the interconversion of adenine nucleotides, which are involved in maintaining the homeostasis of adenine nucleotide composition in various organisms . AKs are proved to be multiple isozymes, such as AK1, AK2, AK3 and AK4. These isozymes are characterized by significantly conserved sequences and their sequences are fairly comparable . The gene for human AK4 was identified in 1992 , which is expressed mainly in tissues rich in mitochondria like the brain, heart, kidney and liver . Abnormal expression of AK4 are proven in a variety of diseases, such as Parkinson's disease , lung tumors and so on. Liu RJ et al. show that AK4 protein levels are increased in cultured cells exposed to hypoxia and in an animal model of amyotrophic lateral sclerosis (ALS), a neurodegenerative disease . In our study, AK4 protein and mRNA were found to be up-regulated in SLE-iPSCs. Western blot and RT-PCR verified the up-regulation of both AK4 protein and mRNA. It can be inferred that AK4 is associated with the development of SLE. However, its mechanism needs in-depth studies with more detailed plans.
In biology process of cells, reproduction is controlled not only by protein-coding genes but also by non-coding regions, including loci that produce small RNAs . MiRNAs, a family of small non-coding RNAs (typically 19–23 nucleotides) , have critical functions over a wide range of biological and pathological processes , such as translational repression, deadenylation, mRNA cleavage [43, 44]. MiRNAs are synthesized from short hairpin precursors and that bind to complementary messenger RNAs (mRNAs) to stimulate their degradation or repress their translation . This complex circuitry of miRNA–mRNA interactions has been studied massively for over a decade in pursuit of identifying miRNA target genes and evaluating the contribution of miRNAs to the regulation of targeted genes . Similarly, several studies have identified miRNA-protein interactions. Trabucchi M et al.  show that RNA-binding protein promotes the biogenesis of a subset of microRNAs. Davis BN  report the smad proteins bind a conserved RNA sequence to promote microRNA maturation by Drosha. These studies indicate proteins that directly bind different subsets of miRNA precursors and enhance their processing by either Drosha, Dicer, or both. Thus, one miRNA can target different mRNAs or proteins, and one mRNA or protein can be under the regulation of multiple miRNAs through a positive-feedback loop , Dimitrios G et al. found that miR-let-7 in conjunction with argonaute protein was found to bind its own primary transcript to enhance its processing by Drosha . Our results indicated that miR-371a-5p, conjunct with AK4 protein, was involved in SLE.
The aim of presenting this report is to help predict miRNA-RNA targets through the integrated analysis of miRNA and mRNA expression profiling, as well as miRNA protein targets using miRNA and protein regulated network. The miRNA level of miR-371a-5p and let-7a-5p were down-regulated in SLE-iPSC and verified using Q-PCR. The up-regulation of AK4 involved in nucleotide biosynthesis may suggest a general acceleration of anabolic metabolism induced by down-regulated miR-371a-5p, which may contribute to SLE. Recent studies demonstrated that miR-371-5p was significantly elevated in gastric cancer (GC) patients and regulated hepatocellular carcinoma (HCC) cell proliferation, the miR-371-5 cluster has proved to be a novel prognostic factor and therapeutic target for tumors and disease [48, 49]. Let-7a is highly conserved among species in sequence and function , which operates as a tumor suppressor in human cells [51–53]. Serguienko A et al.  show that let-7a down-regulates key anabolic enzymes and increases both oxidative phosphorylation and glycolysis in triple-negative breast cancer and metastatic melanoma cell lines. Moreover, let-7a causes mitochondrial ROS production concomitant with the up-regulation of oxidative stress responsive genes. Loss or down-regulation of let-7 levels is associated with increased cancer aggressiveness and poor clinical outcome . This report also indicated that let-7a-5p was down-regulated, which may be involved in SLE development.
In summary, a range of production SLE-iPSCs has been analyzed to investigate the regulatory mechanism. In order to integrate the biological system, global miRNA, mRNA and proteomic expression profiling were employed in parallel. 1099 mRNAs, 223 microRNAs and 94 differentially expressed proteins were verified. Computational and systems biological techniques can be employed to study the regulated networks that mediate physiological system dynamics. To investigate the influence of miRNA on these processes, potential targets of differentially expressed mRNAs and Proteins were predicted using miRanda, TargetScan and Pictar. Gene ontology (GO) analysis of differentially expressed mRNAs and proteins found that identified biological processes drive proliferation including mRNA processing, translation and so on. Western blot and RT-PCR confirmed the up-regulation of AK4 protein and mRNA, as well as the down-regulation of TAGLN protein and mRNA levels in SLE-iPSCs. MiR-371a- 5p and let-7a-5p were down-regulated and verified using RT-PCR. The up-regulation of AK4 involved in nucleotide biosynthesis suggests a general acceleration of anabolic metabolism induced by down-regulated miR-371a-5p, which may contribute to SLE. It is believed that the integration of multiple profiling datasets will provide a useful tool for biologists who are investigating biological networks and seeking to identify direct regulators of mRNAs and proteins. It is also expected that the knowledge gained from this study can assess the usefulness of the pathogenesis and novel biomarker candidates of SLE, which may develop a new way for diagnosis of SLE.
All studies and other procedures were approved by the ethics committee of the Shenzhen People’s Hospital (Shenzhen, China) or the Guangzhou Institutes of Biomedicine and Health (Guangzhou, China). 4 patients (four women, mean age 39, range 30–46) were diagnosed as active SLE with SLEDAI > 8. Patients treated with immunosuppressant within 3 months were excluded. Equivalent subjects, age and sex matched, was recruited as healthy controls. All participating subjects were explained their participation rights. Written informed consent was obtained. Renal tubular cells from urine of participating subjects were reprogrammed to generate human iPSCs clone. One SLE-iPSC clone and control-iPSC clone were identified , and which morphology were identified. All of the samples were collected and immediately frozen in liquid nitrogen and stored at −80 °C.
Sequencing of mRNA and miRNA
Total RNA was extracted from SLE-iPSCs and Control-iPSCs using Trizol (Invitrogen, Carlsbad, USA) according to the manufacturer's instructions. The concentration and quality of total RNA were measured by the UV absorbance at 260 nm and 280 nm (A260/280) and checked by gel electrophoresis. The polyA mRNAs were selected by RNA Purification Beads (Illumina, SanDiego, CA). MicroRNA isolation was carried out from the total RNA using mirVana™ microRNA isolation kit (Ambion, Austin, TX) according to the manufacturer's instructions. Total mRNA and miRNA isolated from three independent cultures were, respectively, pooled for subsequent library construction and sequence analysis. The libraries of mRNA and miRNA were constructed by using the Illumina TruSeq RNA Sample Prep Kit v2-Set A and TruSeq Small RNA Sample Prep Kit Set A, respectively, sequenced by the Illumina HiSeq™ 2000 System. Adaptor sequences were subsequently trimmed to clean full-length reads and formatted into a non-redundant FASTQ format. The high-quality clean reads were mapped to the reference human genome using the SOAP (version 2.0) software. Sequences that perfectly matched the genome along their entire length were considered for next analyses. Differential expression of mRNA and miRNA between SLE-iPSCs and Control-iPSCs were calculated by relative expression analysis.
Protein extraction, identification and relative quantification
The process of protein extraction and iTRAQ sample labeling were performed as described by our previous research . In briefly, the iPSC cells were harvested and washed with ice-cold PBS, then lysed in ice-cold lysis buffer. Protein concentration was determined by Bradford Protein Assay Kit (Pierce). One hundred micrograms of protein from each corresponding group was digested with trypsin, and labeled according to the iTRAQ protocol (Applied Biosystems). The labeled digests were subjected to liquid chromatography coupled with tandem mass spectrometry (LC-MS/MS) analysis. Once dates had been acquired, all data files were processed by Mascot (version 2.3.02) and searched against the IPI_human v3.87 database for peptide identification and quantification, peptide grouping into proteins and protein ratio calculation. Only proteins with ≥ 1 peptides matched a ≥ 1 fold difference in abundance in both directions and a p-value < 0.05 were considered to be differentially expressed protein.
MicroRNA target predictions from differentially expressed mRNA and differentially expressed protein
MiRNA target prediction was performed using three software: PicTar: http://pictar.mdc-berlin.de/; miRanda v5: http://www.ebi.ac.uk/enright%20srv/microcosm/htdocs/targets/v5/; TargetScan 6.2: http://www.targetscan.org/vert/. When predicted at least by two software at the same time, the result was considered reliable. Candidate miRNA-target pairs were collected and further analyzed.
Data integration and analysis
The experiment, identified differentially expressed protein, mRNA and miRNA, was combined with microRNA target predictions to form integrated regulatory networks. Namely, those data used as the input for differential expression identification and miRNA target prediction. Finally, acquired data were analyzed with the KEGGSOAP software (http://www.bioconductor.org/packages/2.4/bioc/html/KEGG SOAP.html). The results of regulatory networks were handled in cytoscape software (http:// cytoscape. org).
The predicted miRNA target mRNA, protein and the differentially expressed mRNA, proteins were subjected to the gene ontology analysis using KEGGSOAP software. The target predictions and differentially expressed dates were mapped to the GO annotation dataset, and the enriched cellular component (CC), biological process (BP) and molecular function (MF) were extracted using the hypergeometric test. a P < 0.05 was considered to be significant.
The SLE-iPSCs and NC-iPSCs were lysed with ice-cold lysis buffer. Proteins were separated by SDS/PAGE, electrophoresed, and transferred onto a PVDF membrane. The membrane was incubated with primary antibodies to anti-TAGLN 1 : 500 (GeneTex) or AK4 1 : 2000 (abcam) at 4 °C overnight, followed by horseradish peroxidase (HRP)-conjugated goat anti-rabbit immunoglobulin G (IgG) 1 : 2000 (southern biotech) at room temperature for 2 h. The antibody-bound proteins were analyzed by using immobilon western chemilum HRP substrate (Millipore). The quantitative analysis of the immunoblot results were performed using Image J software.
Total RNA of SLE-iPSCs and NC-iPSCs was extracted with trizol (Invitrogen) and reversed to cDNA using iScript™ cDNA Synthesis Kit (Bio-Rad). The following specific primers were used: TAGLN, forward: 5’-ATGGCGTGATTCTGAGCAA-3’, reverse: 5’-ATCTGCTTGAAGACCATGG A-3’; AK4, forward: 5’-GCCCAGGCTAATCTATGAAG-3’, reverse: 5’-CAAGGAGCTCAAAAG CCTAT-3’; 18 s rRNA, forward: 5’-CCTGGATACCGCAGCTAGGA-3’, reverse: 5’-GCGGCGCAATACGAATGCC CC-3’; miR-371a-5p, forward: 5’-ACACTCCAGCTGGGTAGCTTATCAGACTGAT G-3’; let-7a-5p, forward: 5’-ACA CTCCAGCTGGGACTCAAACTGTGGGGGC-3’, miRNA, reverse: 5’-CTCAACT GGTGTCGTGGA-3’; U6, forward: 5’-CTCGCTTCGGCAGCACA-3’, reverse: 5’-AA CGCTTCACG AATTT GCGT-3’. Q-PCR was performed on a Real-time PCR system (Bio-Rad) using individual PCR tubes (Bio-Rad). The relative expression level was calculated from a relative standard curve obtained by using log dilutions of cDNA containing genes of TAGLN, AK4, miR-371a-5p, let-7a-5p. The relative expression was quantified using the “2-ΔΔCt” method . The average of three independent analyses for each gene was calculated.
Perl A. Pathogenic mechanisms in systemic lupus erythematosus. Autoimmunity. 2010;43(1):1–6.
de Macedo PA, Borba EF, Viana VDT, Leon EP, Testagrossa LD, Barros RT, Nascimento AP, Bonfa E. Antibodies to ribosomal P proteins in lupus nephritis: a surrogate marker for a better renal survival? Autoimmun Rev. 2011;10(3):126–30.
Herbst R, Liu Z, Jallal B, Yao YH. Biomarkers for systemic lupus erythematosus. Int J Rheum Dis. 2012;15(5):433–44.
Quan JX, Lakhanpal A, Reddy MM, Zaman S, Li QZ, German DC, Olsen NJ, Kodadek T, Karp DR. Discovery of biomarkers for systemic lupus erythematosus using a library of synthetic autoantigen surrogates. J Immunol Methods. 2014;402(1–2):23–34.
Abulaban KM, Brunner HI. Biomarkers for Childhood-Onset Systemic Lupus Erythematosus. Curr Rheumatol Rep. 2015;17:1.
Jiang ZP, Han YM, Cao XT. Induced pluripotent stem cell (iPSCs) and their application in immunotherapy. Cell Mol Immunol. 2014;11(1):17–24.
Guan X, Mack DL, Moreno CM, Strande JL, Mathieu J, Shi YG, Markert CD, Wang ZJ, Liu GH, Lawlor MW, et al. Dystrophin-deficient cardiomyocytes derived from human urine: New biologic reagents for drug discovery. Stem Cell Res. 2014;12(2):467–80.
Thatava T, Armstrong AS, De Lamo JG, Edukulla R, Khan YK, Sakuma T, Ohmine S, Sundsbak JL, Harris PC, Kudva YC, et al. Successful disease-specific induced pluripotent stem cell generation from patients with kidney transplantation. Stem Cell Res Therapy. 2011;2.
Chen YY, Luo RP, Xu Y, Cai XJ, Li WX, Tan KB, Huang JR, Dai Y. Generation of systemic lupus erythematosus-specific induced pluripotent stem cells from urine. Rheumatol Int. 2013;33(8):2127–34.
Yamana R, Iwasaki M, Wakabayashi M, Nakagawa M, Yamanaka S, Ishihama Y. Rapid and deep profiling of human induced pluripotent stem cell proteome by one-shot NanoLC-MS/MS analysis with meter-scale monolithic silica columns. J Proteome Res. 2013;12(1):214–21.
Plath K, Lowry WE. Progress in understanding reprogramming to the induced pluripotent state. Nat Rev Genet. 2011;12(4):253–65.
Chang G, Gao S, Hou XF, Xu ZJ, Liu YF, Kang L, Tao Y, Liu WQ, Huang B, Kou XC, et al. High-throughput sequencing reveals the disruption of methylation of imprinted gene in induced pluripotent stem cells. Cell Res. 2014;24(3):293–306.
Bartel DP. MicroRNAs: Genomics, biogenesis, mechanism, and function. Cell. 2004;116(2):281–97.
Wang YP, Li KB. Correlation of expression profiles between microRNAs and mRNA targets using NCI-60 data. BMC Genom. 2009;10:218.
Cheng C, Li LM. Inferring MicroRNA Activities by Combining Gene Expression with MicroRNA Target Prediction. PloS One. 2008;3:4.
Gerstberger S, Hafner M, Tuschl T. A census of human RNA-binding proteins. Nat Rev Genet. 2014;15(12):829–45.
Takahashi K, Yamanaka S. Induction of pluripotent stem cells from mouse embryonic and adult fibroblast cultures by defined factors. Cell. 2006;126(4):663–76.
Bui YK, Cordes KR, Hastie S, Srivastava D. Cardiac Disease Modeling of Familial Left Ventricular Noncompaction Cardiomyopathy Using Induced Pluripotent Stem Cell. Circulation. 2012;126:21.
Pessach IM, Ordovas-Montanes J, Zhang SY, Casanova JL, Giliani S, Gennery AR, Al-Herz W, Manos PD, Schlaeger TM, Park IH, et al. Induced pluripotent stem cells: A novel frontier in the study of human primary immunodeficiencies. J Allergy Clin Immun. 2011;127(6):1400–U1140.
Araki R, Uda M, Hoki Y, Sunayama M, Nakamura M, Ando S, Sugiura M, Ideno H, Shimada A, Nifuji A, et al. Negligible immunogenicity of terminally differentiated cells derived from induced pluripotent or embryonic stem cells. Nature. 2013;494(7435):100–4.
Mardinoglu A, Nielsen J. Systems medicine and metabolic modelling. J Internal Med. 2012;271(2):142–54.
Somvanshi PR, Venkatesh KV. A conceptual review on systems biology in health and diseases: from biological networks to modern therapeutics. Syst Synth Biol. 2014;8(1):99–116.
De Bandt M. Lessons for lupus from tumour necrosis factor blockade. Lupus. 2006;15(11):762–7.
Anderson CJ, Neas BR, Uchiumi T, Stafford HA. Autoantibodies to the 20-kDa ribosomal proteins: identification, characterization, and new aspects on prevalence in systemic Lupus erythematosus. Clin Immunol. 2001;98(2):249–57.
Sturgill BC, Carpenter RR. Antibody to Ribosomes in Systemic Lupus Erythematosus. Arthritis Rheum. 1965;8:213–8.
Yen MR, Choi J, Saier MH. Bioinformatic Analyses of Transmembrane Transport: Novel Software for Deducing Protein Phylogeny, Topology, and Evolution. J Mol Microb Biotech. 2009;17(4):163–76.
Rudin DO. The choroid plexus and system disease in mental illness. II. Systemic lupus erythematosus: a combined transport dysfunction model for schizophrenia. Biol Psychiat. 1981;16(4):373–97.
Lees-Miller JP, Heeley DH, Smillie LB, Kay CM. Isolation and characterization of an abundant and novel 22-kDa protein (SM22) from chicken gizzard smooth muscle. J Biol Chem. 1987;262(7):2988–93.
Camoretti-Mercado B, Forsythe SM, LeBeau MM, Espinosa 3rd R, Vieira JE, Halayko AJ, Willadsen S, Kurtz B, Ober C, Evans GA, et al. Expression and cytogenetic localization of the human SM22 gene (TAGLN). Genomics. 1998;49(3):452–7.
Shields JM, Rogers-Graham K, Der CJ. Loss of transgelin in breast and colon tumors and in RIE-1 cells by Ras deregulation of gene expression through Raf-independent pathways. J Biol Chem. 2002;277(12):9790–9.
Bregant E, Renzone G, Lonigro R, Passon N, Di Loreto C, Pandolfi M, Scaloni A, Tell G, Damante G. Down-regulation of SM22/transgelin gene expression during H9c2 cells differentiation. Mol Cell Biochem. 2009;327(1–2):145–52.
Li LS, Kim H, Rhee H, Kim SH, Shin DH, Chung KY, Park KS, Paik YK, Chang J, Kim H. Proteomic analysis distinguishes basaloid carcinoma as a distinct subtype of nonsmall cell lung carcinoma. Proteomics. 2004;4(11):3394–400.
Klade CS, Voss T, Krystek E, Ahorn H, Zatloukal K, Pummer K, Adolf GR. Identification of tumor antigens in renal cell carcinoma by serological proteome analysis. Proteomics. 2001;1(7):890–8.
Yoneda T, Sato M, Maeda M, Takagi H. Identification of a novel adenylate kinase system in the brain: cloning of the fourth adenylate kinase. Brain Res Mol Brain Res. 1998;62(2):187–95.
Schulz GE, Schiltz E, Tomasselli AG, Frank R, Brune M, Wittinghofer A, Schirmer RH. Structural relationships in the adenylate kinase family. Eur J Biochem. 1986;161(1):127–32.
Xu G, O'Connell P, Stevens J, White R. Characterization of human adenylate kinase 3 (AK3) cDNA and mapping of the AK3 pseudogene to an intron of the NF1 gene. Genomics. 1992;13(3):537–42.
Panayiotou C, Solaroli N, Karlsson A. The many isoforms of human adenylate kinases. Int J Biochem Cell B. 2014;49:75–83.
Garcia-Esparcia P, Hernandez-Ortega K, Ansoleaga B, Carmona M, Ferrer I. Purine metabolism gene deregulation in Parkinson's disease. Neuropathol Appl Neurobiol. 2015;41(7):926-40.
Greengard O, Head JF, Goldberg SL. Uridine kinase, adenylate kinase, and guanase in human lung tumors. Cancer Res. 1980;40(7):2295–9.
Liu RJ, Strom AL, Zhai JJ, Gal J, Bao SL, Gong WM, Zhu HN. Enzymatically inactive adenylate kinase 4 interacts with mitochondrial ADP/ATP translocase. Int J Biochem Cell B. 2009;41(6):1371–80.
Zhuang X, Li ZM, Lin HN, Gu L, Lin Q, Lu ZX, Tzeng CM. Integrated miRNA and mRNA expression profiling to identify mRNA targets of dysregulated miRNAs in non-obstructive azoospermia. Sci Rep-Uk. 2015;5:7922.
He ZP, Kokkinaki M, Pant D, Gallicano GI, Dym M. Small RNA molecules in the regulation of spermatogenesis. Reproduction. 2009;137(6):901–11.
Wilczynska A, Bushell M. The complexity of miRNA-mediated repression. Cell Death Differ. 2015;22(1):22–33.
Bak RO, Mikkelsen JG. miRNA sponges: soaking up miRNAs for regulation of gene expression. Wiley Interdisciplinary Rev RNA. 2014;5(3):317–33.
Trabucchi M, Briata P, Garcia-Mayoral M, Haase AD, Filipowicz W, Ramos A, Gherzi R, Rosenfeld MG. The RNA-binding protein KSRP promotes the biogenesis of a subset of microRNAs. Nature. 2009;459(7249):1010–U1144.
Davis BN, Hilyard AC, Nguyen PH, Lagna G, Hata A. Smad Proteins Bind a Conserved RNA Sequence to Promote MicroRNA Maturation by Drosha. Molecul Cell. 2010;39(3):373–84.
Zisoulis DG, Kai ZS, Chang RK, Pasquinelli AE. Autoregulation of microRNA biogenesis by let-7 and Argonaute. Nature. 2012;486(7404):541–U140.
He D, Miao HL, Xu YM, Xiong LH, Wang Y, Xiang HX, Zhang H, Zhang ZY. MiR-371-5p facilitates pancreatic cancer cell proliferation and decreases patient surviva. PloS One. 2014;9:11.
Liu RY, Diao CF, Zhang Y, Wu N, Wan HY, Nong XY, Liu M, Tang H. miR-371-5p down-regulates pre mRNA processing factor 4 homolog B (PRPF4B) and facilitates the G1/S transition in human hepatocellular carcinoma cells. Cancer Lett. 2013;335(2):351–60.
Pasquinelli AE, Reinhart BJ, Slack F, Martindale MQ, Kuroda MI, Maller B, Hayward DC, Ball EE, Degnan B, Muller P, et al. Conservation of the sequence and temporal expression of let-7 heterochronic regulatory RNA. Nature. 2000;408(6808):86–9.
Serguienko A, Grad I, Wennerstrom AB, Meza-Zepeda LA, Thiede B, Stratford EW, Myklebost O, Munthe E. Metabolic reprogramming of metastatic breast cancer and melanoma by let-7a microRNA. Oncotarget. 2015;6(4):2451–65.
Liu YC, Yin BD, Zhang CC, Zhou LB, Fan J. Hsa-let-7a functions as a tumor suppressor in renal cell carcinoma cell lines by targeting c-myc. Biochem Bioph Res Co. 2012;417(1):371–5.
Guled M, Lahti L, Lindholm PM, Salmenkivi K, Bagwan I, Nicholson AG, Knuutila S. CDKN2A, NF2, and JUN Are Dysregulated Among Other Genes by miRNAs in Malignant Mesothelioma-A miRNA Microarray Analysis. Gene Chromosome Canc. 2009;48(7):615–23.
Nair VS, Maeda LS, Ioannidis JPA. Clinical Outcome Prediction by MicroRNAs in Human Cancer: A Systematic Review. J Natl Cancer Inst. 2012;104(7):528–40.
Lan P, Li WF, Wen TN, Shiau JY, Wu YC, Lin WD, Schmidt W. iTRAQ Protein Profile Analysis of Arabidopsis Roots Reveals New Aspects Critical for Iron Homeostasis. Plant Physiol. 2011;155(2):821–34.
Schmittgen TD, Livak KJ. Analyzing real-time PCR data by the comparative C(T) method. Nat Protoc. 2008;3(6):1101–8.
We thank the patients and healthy volunteers who participated in this study.
Y D and DE T designed the experiments, analyzed the data and wrote the manuscript with input from all authors. DE T, YY C and JR H executed the experiments. WB C and WJ P performed the bioinformatics analysis. HY H contributed to the promoter analysis, QJ L and DE T led the analysis and interpretation of the data, assembled the figures and wrote the manuscript. All authors read and approved the final manuscript.
The authors declare they have no competing interests.
Table S1. Differentially expressed mRNA. (XLSX 271 kb)
Table S2. Differentially expressed microRNA. (XLSX 22 kb)
Table S3. Up-regulation of differentially expressed proteins. (XLSX 11 kb)
Table S4. Down-regulation of differentially expressed proteins. (XLSX 15 kb)
Table S5. KEGG analysis of differentially expressed mRNA. (XLSX 28 kb)
Table S6. KEGG analysis of differentially expressed protein. (XLSX 16 kb)
About this article
Cite this article
Tang, D., Chen, Y., He, H. et al. Integrated analysis of mRNA, microRNA and protein in systemic lupus erythematosus-specific induced pluripotent stem cells from urine. BMC Genomics 17, 488 (2016) doi:10.1186/s12864-016-2809-9
- Systemic lupus erythematosus (SLE)
- Induced pluripotent stem cells (iPSCs)