Genomic analysis of host – Peste des petits ruminants vaccine viral transcriptome uncovers transcription factors modulating immune regulatory pathways
- Siddappa Manjunath†1,
- Gandham Ravi Kumar†1, 2Email author,
- Bishnu Prasad Mishra1,
- Bina Mishra1,
- Aditya Prasad Sahoo1,
- Chaitanya G Joshi3,
- Ashok K Tiwari1,
- Kaushal Kishore Rajak4 and
- Sarath Chandra Janga2, 5, 6Email author
© Manjunath et al.; licensee BioMed Central. 2015
Received: 8 September 2014
Accepted: 16 January 2015
Published: 24 February 2015
Peste des petits ruminants (PPR), is an acute transboundary viral disease of economic importance, affecting goats and sheep. Mass vaccination programs around the world resulted in the decline of PPR outbreaks. Sungri 96 is a live attenuated vaccine, widely used in Northern India against PPR. This vaccine virus, isolated from goat works efficiently both in sheep and goat. Global gene expression changes under PPR vaccine virus infection are not yet well defined. Therefore, in this study we investigated the host-vaccine virus interactions by infecting the peripheral blood mononuclear cells isolated from goat with PPRV (Sungri 96 vaccine virus), to quantify the global changes in the transcriptomic signature by RNA-sequencing. Viral genome of Sungri 96 vaccine virus was assembled from the PPRV infected transcriptome confirming the infection and demonstrating the feasibility of building a complete non-host genome from the blood transcriptome. Comparison of infected transcriptome with control transcriptome revealed 985 differentially expressed genes. Functional analysis showed enrichment of immune regulatory pathways under PPRV infection. Key genes involved in immune system regulation, spliceosomal and apoptotic pathways were identified to be dysregulated. Network analysis revealed that the protein - protein interaction network among differentially expressed genes is significantly disrupted in infected state. Several genes encoding TFs that govern immune regulatory pathways were identified to co-regulate the differentially expressed genes. These data provide insights into the host - PPRV vaccine virus interactome for the first time. Our findings suggested dysregulation of immune regulatory pathways and genes encoding Transcription Factors (TFs) that govern these pathways in response to viral infection.
Peste des petits ruminants (PPR), an acute viral disease of goats and sheep, characterized by fever, erosive stomatitis, conjunctivitis, gastroenteritis and pneumonia is caused by Peste des petits ruminants virus (PPRV). PPRV, a linear single stranded, negative sense RNA virus, belongs to the genus Morbilivirus, family Paramyxoviridae . The PPRV genome of 15948 bp encodes six structural proteins, nucleocapsid (N), phosphoprotein (P), matrix (M), fusion (F), hemagglutinin (H) and large (L) protein in 3′ to 5′ direction (3′-N-P-M-F-H-L-5′); and two non-structural proteins (C/V) due to RNA editing of the phosphoprotein gene. Morbidity and mortality in this highly contagious disease is as high as 100% and 90%, respectively [1,2]. PPR was first reported in the Ivory Coast in 1942 and is now widespread across sub-Saharan Africa, Arabian Peninsula, Middle East and the Indian sub-continent. Given its economic relevance and severity, the disease is classified as a World Organisation for Animal Health (OIE) listed disease. Four different lineages of this virus (I–IV) have been defined all over the world based on molecular epidemiology of “N”, “F” and “H” gene sequences of the virus . Viruses from Asia are predominantly classified into lineage IV and hence refereed as Asian lineage. However, few recent reports indicate that Asian lineage is now present in Africa as well .
In India, this transboundary contagious disease is endemic, causing severe economic losses to small ruminant production. PPR outbreaks have increased between the year 1996 and 2005 . The economic loss was estimated to be around 1800 million Indian rupees (US$ 39 million) every year in India . The number of outbreaks after 2005 has declined due to the mass vaccination efforts carried out from 2004 under the Assistance to States for Control of Animal Diseases programme, funded by the government of India. Three attenuated live cell culture vaccines, Sungri 96, Arasur 87 and Coimbatore 97 are commercially available in India. Sungri 96 (isolate of goat origin) was developed by Indian Veterinary Research Institute, Arasur 87 (isolate of sheep origin) and Coimbatore 97 (isolate of goat origin) were developed by Tamil Nadu Veterinary and Animal Sciences University. All these vaccines are considered to be potent and safe .
Sungri 96 vaccine, widely used in Northern India, was developed by attenuating the Sungri isolate of PPRV up to 60th passage in Vero cells . This virus, isolated from goat does not have a restriction in inducing immunity in either of the species, sheep or goat. Single dose of PPR live attenuated vaccine (Sungri 96) contains ~103 TCID50 of Vero cell-attenuated PPRV and elicits robust immune response for up to 78 months , indicating that a single vaccination is sufficient to provide long-term immunity in sheep and goats. Cell-mediated immunity is suggested to play a major role in eliciting this response, which warrants further investigation. PPRV like other morbilliviruses is cell associated and reaches the other organs and the tissues by piggybacking on Peripheral blood mononuclear cells (PBMCs) . PBMCs consist of immune cells T- and B-lymphocytes, as well as dendritic cells, which express SLAM, a (co) receptor for morbilliviruses. PBMCs are widely used as a model for studying the immunosuppressive nature, immune response and molecular pathways triggered following morbillivirus infection [9,10]. SLAM constitutively expressed on these cells facilitates the entry of Morbilliviruses. Basal level of SLAM expression in PBMCs has been shown to have high correlation with their ability to replicate PPRV . Apart from SLAM, morbilliviruses use alternate receptors for entry. In addition to the SLAM, which acts as cellular receptor for PPR virus on immune cells, Ovine nectin-4 was identified as a novel epithelial receptor for PPR virus, which highlights on pathogenesis factors like tissue distribution and genetic variation of morbillivirus receptors . CD46 (a membrane co-factor protein) acts as a cellular receptor for vaccine and laboratory passaged strains of measles virus. Understanding the molecular mechanisms of PPR vaccine virus (Sungri-96) induced protection, using PBMCs as a model will shed light on immune genes and cellular pathways involved in the protective mechanisms in hosts. This understanding of the host virus interactome would pave way for developing novel tools for better clinical management of PPR.
Recently, global transcriptome analysis was used to understand the molecular events of host vaccine interactome in Measles virus, Rinderpest virus [9,13]. However, no reports on PPRV-host interactome are available in the literature till date. Therefore, in the present study, as PPRV is lymphotropic we infected the peripheral blood mononuclear cells (PBMCs) isolated from goat with Sungri 96 vaccine virus and investigated the host-virus crosstalk vis-à-vis changes in global gene expression signatures at cellular level.
Materials and methods
Ethics statement and animals
The experimental procedures in the present study were approved by Institute Animal Ethics Committee (I.A.E.C No.F.1.53/2012-13-J.D.). Goat kids (5 months old) used for blood collection were housed in appropriate containment facilities with feed and water adlibitum.
Virus and PBMC culture
PPRV strain-Sungri/96, a vaccine virus adapted to Vero cells (African Green Monkey Kidney cells), available at the Division of Biological Products, IVRI, Izatnagar and maintained at the passage level 60–62 was used in the present study. Blood was collected from goats (n = 5) in heparin coated vacutainer vials. The whole blood was layered on 3 mL of histopaque-1077 and was centrifuged at 2200 rpm for 40 min. The interphase layer rich in peripheral blood mononuclear cells (PBMCs) was transferred into a separate tube and was washed three times with RPMI-1640 medium at 1800 rpm for 10 min. The final pellet was re-suspended in a complete medium containing RPMI-1640 and 10% fetal calf serum (FCS). The live cells were counted in haemocytometer and 1 × 106 cells per mL per well were seeded into two six well plates for culture in 5% CO2 incubator at 37 °C.
PBMCs infection with PPRV and cytopathic effect
One plate served as a control (uninfected cells) and the other was used for infection. Sungri-96 vaccine virus, was grown in Vero cells and purified by banding on sucrose gradient (ultracentrifuged). The purified virus was titrated and tested for its infectivity in Vero cells and subsequently used in the experiment to infect PBMCs. Goat PBMC’s were infected with purified PPR virus (Sungri/96) at MOI of 1.0 in RPMI medium and incubated at room temperature for 1 h of adsorption. The cells were then washed, resuspended in complete RPMI medium and incubated at 37 °C in 5% CO2 incubator for 120 h. Cells were observed for morphological changes like clumping and ballooning at 24 h, 72 h and 120 h post infection (pi). As mRNA for N protein, is most abundantly transcribed and is required for viral replication and infection, at each of these time points, PCR and real-time PCR (qRT-PCR) for N gene expression was done to evaluate viral infection. Total RNA from cells was extracted using RNeasy mini Kit (Qiagen). The concentration and quality of the extracted RNA samples was determined using Bioanalyzer 2100 (Agilent Technologies Inc). The cDNA was synthesized using Revert Aid First Strand cDNA Kit (Thermo Scientific-#K1622) following manufacturers instructions. Fold change in qRT-PCR was calculated taking 24 h time point as calibrator. It was observed that the expression of N gene increased significantly from 24 h pi to 120 h pi. The viral infection in PBMCs at 120 h pi, was further confirmed by identifying PPRV specific proteins by western blot using polyclonal serum raised in goats against the Sungri/96 vaccine strain. Briefly, the cells after 120 h pi were harvested and lysed in lysis buffer (10 mM Tris/HCl, pH 8.0, 150 mM NaCl, 1% TritonX-100 and 1 mM DTT supplemented with protease and phosphatase inhibitors (2 mM sodium orthovanadate, 100 nM okadaic acid, 1 mM NaF, 1 mM β-glycerophosphate and cocktail (Sigma)). The protein samples were separated on 12% SDS-PAGE and electroblotted onto nitrocellulose membrane. After blocking overnight in 2% BSA the membrane was incubated with diluted polyclonal antiserum raised in goats (1: 100) for 1 h. The membrane was then washed three times with PBST buffer and then incubated in anti-goat IgG conjugated with HRP at a dilution of 1:10 000 for 1 h at 37 °C. The membrane was then washed thrice with PBST for 10 min and bands were visualized following incubation with diaminobenzidine tetrahydrochloride (DAB system, GeNei). PPRV specific bands confirmed viral infection in PBMCs after 120 h pi.
Library preparation and RNA sequencing
Infected cells were harvested at 120 h pi along with the mock control cells. The harvested cells were pelleted at 2000 rpm for 10 min and stored in RNA later at −80 °C for RNA isolation. Total RNA was isolated from infected and mock-control cells using RNeasy Kit (Qiagen) according to the manufacturer’s instructions. RNA was quantified and checked for RNA integrity number (RIN) using Agilent 2100 Bioanalyser. As the RIN number for both the infected and control total RNA was ≥ 8.0, cDNA library was prepared using Ion Total RNA-Seq Kit according to the standard protocol and purified by magnetic bead module. The average fragment size of the library was observed to be 377 bp on the bioanalyzer. The beads were amplified by em-PCR after calculating the bead to fragment ratio. The amplified beads were enriched in Ion one Touch ES and the enriched template positive Ion sphere particles were subjected to sequencing using Ion 318 chip. The RNA sequencing data generated in FASTQ format was used for further analysis.
Transcriptome quantification and functional analysis of differentially expressed genes
Viral genome assembly
As PPRV is an RNA virus we could assemble the whole genome of PPRV (Sungri/96) from the infected PBMC transcriptome and not from the control transcriptome. The quality processed reads were aligned with the available PPRV reference genome from NCBI using bowtie 2.0 . The aligned reads were assembled using MacVector 12.7.5 (@ 2012 MacVector Inc).
Protein-Protein interaction network among the differentially expressed genes
The Biological General Repository for Interaction Datasets (BioGRID) is a curated biological database of protein-protein and genetic interactions created in 2003. It provides a comprehensive resource of protein–protein and genetic interactions for all major model organism species. BioGRID currently holds 347 966 interactions (170 162 genetic, 177 804 protein) curated from both high-throughput data sets and individual focused studies derived from over 23 000 publications in the primary literature . In this repository, protein-protein interactions in human are well defined and those in Bos taurus are very few. Since protein interactions have been shown to be well-conserved across species , g:Orth in g:Profiler web server was employed to produce orthology (functionally equivalent genes) predictions between species to facilitate functional and interaction annotation transfer across species . To construct the protein-protein interaction network with the differentially expressed genes in the present study, Bos taurus orthologs in human were queried using g:Orth. Customized perl scripts were used to extract interactions involving the differentially expressed genes. The complete interaction network visualized in Cytoscape 3.0.2  was found to have an average degree of 4.4 with 6756 nodes and 30 023 edges (Additional file 2). Degree or connectivity of a node/gene in a network is the number of connections it has to other nodes/genes. The network involving only the 985 differentially expressed genes was found to comprise of 905 nodes and 5340 edges (Additional file 3). The network properties of each of these networks including clustering coefficient and graph density were calculated using igraph package in R . Out of the 985 differentially expressed genes, 105 genes with 3 fold change in expression (positive or negative) and a degree of at least 5 were filtered to generate a sub-network using the original protein - protein interaction network. This set of 105 genes was designated as differentially expressed highly connected (DEHC) gene network.
Extracting 5 kb upstream of DEHC genes
Ensembl is a comprehensive database addressing the challenges of decoding a eukaryotic genome from the set of functional elements it represents to providing access to the vast sea of data . Ensemble BioMart is a hub for data retrieval across the taxonomic space. We extracted 5 kb upstream of 105 DEHC genes from Ensemble BioMart using the Bos taurus ensemble IDs. The sequences containing a string of undefined bases (N’s) due to partially complete contigs were eliminated from further analysis using a customized perl script.
DNA motif discovery using MEME and prediction of transcription factors binding to upstream regions of DEHC genes
MEME (Multiple EM for Motif Elicitation) Suite is a comprehensive collection of tools most widely used for discovery of new transcription factor and protein domain binding sites . The motif discovery algorithm MEME, finds ungapped motifs within DNA or protein sequences. TOMTOM  tool in the suite allows comparing the discovered motif to a database of motifs (JASPAR, JOLMA, etc.) to find matches to the established Position Weight Matrices (PWMs) of transcription factors. To predict TFs that bind to the upstream of the 105 DEHC genes identified in the present study, initially over-represented conserved motifs among these genes were identified using MEME followed by TOMTOM using a locally installed version of MEME suite. The TFs identified were mapped onto their orthologs in Bos taurus. The presence of these TF binding sites across the DEHC genes was depicted by a heatmap using cluster and Java tree view.
Validation by Quantitative Real time PCR (qRT-PCR)
Quantitative Real time PCR was carried out on the same biological material that was used in RNA-Seq experiment. RNA was extracted from the harvested cells using RNeasy mini kit (Qiagen) and was quantified using nanodrop spectrophotometer (Thermo Scientific). cDNA was synthesized using Revert Aid First Strand cDNA synthesis kit according to the manufacturer’s instructions and qRT-PCR was performed using Applied Biosystems 7500 Fast system using 2X SYBR Green Master mix (USB, Sigma). Some DEHC genes were validated using GAPDH as an endogenous control. A panel of six housekeeping genes was tested for their stable expression across the infected and control samples at different time points in the experiment. RefFinder was used to generate stability values for all genes. GAPDH, was found to be the most stably expressed gene on analysis (Data not shown).
Genes and their primer sequence for validation
Validation by qRT-PCR in goats vaccinated with Sungri/96 PPR vaccine
Animal experimentation in the present study was conducted under the approval mentioned above. Goats (n = 6, 5–6 months old), three served as control and three were used for vaccination with Sungri/96 vaccine virus. The goats were previously screened and were found to be negative for PPRV antibodies using c-ELISA. PBMCs were collected from both the groups (control and vaccinated) on 5th day post-vaccination and differentially expressed genes validated in the in vitro study, were further evaluated for their expression in the in vivo experiment by qRT-PCR. The vaccinated animals were tested for presence of PPRV antibodies using c-ELISA and were found to be positive based on percentage inhibition (PI) values.
PCR, Real time PCR, Western blot analysis and Viral genome assembly confirm viral infection in PBMCs
Functional analysis of differentially expressed genes reveals enrichment for immune regulatory pathways
The enriched biological pathways in infected PBMCs were identified using DAVID. Among the 985 differentially expressed genes (DEGs), DAVID provided functional annotation for 896 and 921 genes, associated to Bos taurus and Homo sapiens, respectively. A total of 268 genes were mapped to 14 statistically significant categories (p < 0.05; Figure 3B, Additional file 6) among the documented canonical pathways found in KEGG. Ribosome, proteasome and T-cell signaling pathways were the top three annotated with 36, 14 and 22 genes, respectively. There were also a statistically significant number of mapped genes representing antigen processing and presentation, spliceosome, chemokine and JaK-STAT signaling pathways. When mapped to the signaling pathways in BioCarta, 128 differentially expressed genes representing 18 statistically significant categories with most of them related to immune system regulation were detected (p < 0.05; Figure 3B, Additional file 6). Out of the 985 differentially expressed, a total of 576 genes could also be mapped to 19 significant (p < 0.05) categories in the Reactome database. Among these categories, gene expression, metabolism of proteins and 3′-UTR-mediated translational regulation were found to be over-represented with 90, 69 and 48 genes, respectively (Figure 3B, Additional file 6).
Dense cross-talk in the protein - protein interaction network of differentially expressed genes is disrupted in infected state
The global interaction network of the 985 differentially expressed genes (GN) and the network involving interactions only between the 985 differentially expressed genes (DEN) were both constructed and analyzed (see Materials and methods). Both the graph density and the average clustering coefficient for the DEN was found to be much higher than that observed for GN (Graph Density: 0.013 vs 0.001; and Clustering Coefficient: 0.066 vs 0.026). Among the 985 differentially expressed genes 407 genes were upregulated and 578 genes were downregulated. The average connectivity was found to be 15.43 and 9.25 for the up- and downregulated genes, respectively. Though, the degree distribution was significantly different (p = 2.983e-05; Wilcox test), functional enrichment of the up and downregulated genes, independently, revealed significant enrichment of similar cellular, metabolic and immune regulatory processes.
Genes in the DEHCnet involved in cell profileration, apoptosis and immune regulatory pathways
Plays key role in cellular growth and differentation. Suppressor of Apoptosis
Promotes cell cycle progression and inhibits induction of apoptosis or senescence.
Role in cell cycle progression.
Interacts with Ki-67 in proliferating cells and has a role in mitosis
Regulator of cell proliferation. Marker for breast cancer and acute lymphoblastic leukemia.
Role in apoptosis. Reacts with CDK11p46 and sensitizes cells for Fas ligand-induced apoptosis
Regulates cell cycle progression and apoptosis.
Represses function of Tim3 expressed on T cells during Hepatitis C virus infection. Regulates Apoptosis.
Regulates signaling cascades that suppress growth or induces apoptosis.
Induces antiviral response (anti-viral signaling) by inducing IFN-alpha.
Involved in T cell development. Mutation in this gene leads to Immunodeficiency (SCID)
Signaling molecule Involved in synaptic adhesion. Processed to form ß-amyloid peptides
Chaperon protein that regulates DNA damage and repair
In addition we found several major classes of regulators at transcriptional, signaling and post-transcriptional levels dysregulated in the DEHC network. For instance, 3 transcriptional factors viz. CAT, IRF3 and YY1 as well as important signaling proteins, TRRAP – a phosphoinositide 3-kinase, CARD11- a membrane-associated guanylate kinase, PRPF4B – a CDK like kinase and CSNK2A1- casein kinase II were found to be significantly altered. Also, our DEHCNet analysis uncovered 16 RNA binding proteins (RBPs) viz. RPS8, RPS10, RPS18, SRSF1, RPL18A, RNPS1, DEK, EIF3D, RPL4, TRA2B, MKI67IP, RBM10, HUWE1, LBR, PRPF4B and RRP12 that are involved in post transcriptional regulation of gene expression .
Dysregulated TFs contribute to rewiring the expression landscape of the infected PBMCs
Motifs upstream of 105 DEHC genes and predicted TFs that bind to these motifs
NFKB1, RELA, MOT3, NFATC2, dl_2, dl_1, REL, Stat3, E2F7_DBD*, NKX2-8_DBD*, NFATC1_full_3*, NKX2-8_full*, ZNF306_full*, ZNF410_DBD*, DPRX_DBD_2*
EBF1, REST, ADR1, YPR022C, REST, Smad3_secondary, RUNX2_DBD_2*, RUNX3_DBD_3*, Vdr_DBD*, EBF1_full*, VDR_full*
Hoxc12_3480.1, ADR1, usp, Hoxc10_2779.2, Sp4_secondary, Rfx3_secondary, FOXK1_DBD*, RARG_DBD_3*, HOXC12_DBD_2*, HOXC11_full*, HOXD12_DBD_2*
Zfp691_secondary, NRG1,ESR1, Smad3_secondary, Pax5, Zbtb7b_secondary, Myf6_secondary, Zic3_secondary, TFAP2A_DBD_3*, TFAP2C_full_2*, KLF14_DBD*, KLF13_full*, Tcfap2a_DBD_3*, SP8_DBD*
Sp4_secondary, Egr1_secondary, FOXK1_DBD*, PAX5_DBD*, PAX2_DBD*, PAX1_DBD*, Foxk1_DBD*
RFX1, Hic1_primary, Duxl_1286.2, FEV, Titf1_1722.2, IRF2, Irf3_primary, Gamyb, Nkx2-4_3074.1, PPARG, IRF3_full*, Hic1_DBD*, Hic1_DBD_2*, NKX2-8_DBD*, NKX2-8_full*
ELF5, Rfxdc2_secondary, NFATC2, REI1, FEV, YRM1, SPI1, STAT1, EBF1, Elf3_primary, IRF3_full*, FOXB1_DBD*, ETV6_full*
Zfp187_secondary, Ddit3::Cebpa, SOK2, EBF1, Hnf4a_primary, Zic1_secondary, PUT3,Zic2_secondary, Esrra_secondary, Irf4_primary, ZNF524_full_2*, ZIC4_DBD*, ZIC1_full*, RORA_DBD*, Zic3_DBD*
NF-kappaB, EBF1, SPI1, PUT3, GCR2, Bapx1_2343.1, Nkx3-1_primary, Lag1, GCR1, REL, EBF1_full*, Tcfap2a_DBD_3*, ZBTB7B_full*, MAFK_DBD_2*
SPI1, dl_1, Ddit3::Cebpa, RELA, EBF1, NF-kappaB,dl_2,REL, NFKB1,PUT3, FOXO3_full_3*, DPRX_DBD_2*, RHOXF1_DBD*, RHOXF1_full*, RHOXF1_full_2*
Validation of RNA sequencing data by real-time RT-PCR
Validation of RNA sequencing data and fold change comparison
RNA Sequencing fold change
Real-time RT-PCR fold change
In this study our approach to dissect the host-virus interactome broadly comprised of the following steps as shown in the flowchart (Figure 1) and discussed in Materials and methods: 1) Confirming the viral infection 2) RNA-sequencing for comparative transcriptomic analysis in infected and uninfected goat PBMCs 3) Functional analysis to identify the processes enriched in differentially expressed genes and uncovering the associated protein interaction network and 4) Discovery of the transcription factors contributing to the dysregulation of genes.
Viral infection in the PBMCs infected with PPRV at all time points (24 h, 72 h and 120 h pi) was confirmed by amplification and expression of N gene, which is most abundantly transcribed and required for PPR viral replication and infection . The cytopathic changes in the PBMCs infected with PPRV corroborated the view that PBMCs were infected by the virus. Similar cytopathic effect was reported from 4th day pi, in Sungri vaccine virus strain infected Vero cells . PPRV infection in PBMCs was further confirmed by western blot analysis of PPRV specific proteins using antibody raised against purified virus as documented earlier . The assembly of the PPR viral genome from the RNA-seq reads of the infected PBMCs ratified the viral infection. Our ability to build the complete viral genome based on the transcriptome of the infected PBMCs strengthens the notion that the transcriptomes of blood as well as other body fluids could be used as surrogates to know the extent of infection by various pathogens. To our knowledge this is the first global analysis demonstrating the feasibility of building a complete non-host genome from a blood transcriptome.
Analysis of transcriptome data generated from 120 h pi of PBMCs infected with PPRV revealed 985 differentially expressed genes. On assigning the gene ontology terms using g- profiler 117 genes exhibited a significant enrichment to the immune system processes. The enrichment of T-cell signaling , chemokine , antigen processing and presentation , Jak-STAT , IL-7  pathways on annotation by DAVID, indicated the involvement of immune response regulatory pathways in PBMCs infected with PPRV. In this study, we found that PPRV infection induced, interferon regulatory factors - IRF4 and IRF5, interferon induced tetricopeptide - IFIT3 (ISG60) and tripartite motif protein - TRIM56. It has been reported that Interferon regulatory factors and TRIM56 are activated in response to invading viruses (ssRNA or dsRNA viruses) and lead to increase in transcription of interferon stimulated genes (ISGs) by binding to the interferon stimulated responsive elements (ISREs). TRIM56, a virus inducible E3 ubiquitin ligase, was found to restrict bovine viral diarrhea virus replication by inducing ISGs that are transcriptionally regulated by IRFs . Thus, TRIM proteins together with IRFs play an important role in innate immunity to virus infections by stimulating interferon stimulated genes . This suggests that TRIM56, IRF4/5 may act in restricting PPRV replication by inducing ISGs (IFIT3/ISG60) thus augmenting the innate immune response. However, to identify the role of transcriptome signatures in understanding cross protection with different PPRV strains, independent studies have to be conducted for each of the strains under consideration to identify common signatures across different strains. These common signatures may then help in understanding cross protection within different PPRV strains as studied in other viral strains .
Protein interaction network (interactome) analysis provides an effective way to understand the interrelationships between genes . The graph density and the average clustering coefficient for the network involving interactions only between the 985 differentially expressed genes (DEN) was found to be much higher than that observed for the global interaction network of the 985 differentially expressed genes (GN) suggesting that the network of interactions between differently expressed genes is much more intricately connected and is likely contributing to the observed phenotype due to the disruption in the interactome. High clustering coefficient of the DEN compared to the GN also suggests that the former network is more modular. The average connectivity among the 407 genes that were upregulated and 578 genes that were downregulated among the differentially expressed genes was significantly different. However, similar functional enrichment of the up and downregulated genes, independently, suggested no evidence of functional differences in the gene pool contributing to these differences. Nevertheless, it is possible to speculate from this data that up-regulated genes due to their already high number of interactions might be increasing their interactome while down-regulated genes might be losing their functionality thereby contributing to the overall disruption of the protein-protein interactions in the infected state.
The DEHC network represented genes with a fold change ≥ 3 and degree ≥ 5 among the differentially expressed genes. Genes involved in cell proliferation such as EP400, NUB1, MK1671P and TPD52L2 were found to be downregulated while genes such as RPS8, PP1CA, BAG6 and PPP5C involved in progression of apoptosis were upregulated (Table 2). VIM, Viementin-involved in maintaining cell shape and integrity was also found to be downregulated. These findings were in line with the ability of PPRV virus to induce apoptosis in goat PBMCs . IFIT3 that induces IFN-γ antiviral response, CD3D that is involved in T-cell development and GNAI3, critical for normal B cell function were all upregulated reasserting the involvement of the immune response regulatory pathways in PPRV infected PBMCs.
APP, which is a precursor of β-amyloid protein - known for its pro-inflammatory activity, was downregulated. Our analysis identified 147 and 25 interacting proteins for APP among the 985 differentially expressed genes and 105 DEHC genes, respectively (Additional file 8). This observation supports the notion that some viruses suppress the expression of proinflammatory cytokines thereby potentially impeding the antiviral host response to infection . Chaperone protein, VCP that regulates DNA damage and repair was also upregulated. It was found to interact with 16 proteins in the DEHCNet and 77 proteins among the 985 differentially expressed genes. It was recently established that VCP is a host factor required for viral RNA replication of polio virus and that it provides a novel link between cellular protein secretion and viral RNA replication  suggesting its probable role in PPRV infected PBMCs. Further, the presence of 16 dysregulated RNA Binding Proteins (RBPs) reaffirmed the enrichment of spliceosome pathway among the KEGG pathways.
Transcription factors (TFs) are key regulators of gene expression and 12 TFs were identified among the 985 differentially expressed genes. TFs, IRF3 , IRF4 , MAF , NFKB1 , PRDM1 , SP1 , STAT1  and YY1 are involved in the immune regulatory pathways strongly supporting the observation that immune system associated genes are differentially expressed in PPRV infected PBMCs.
An analysis of the overlap of the 41 TFs identified from footprinting the upstream regions above, among the 985 differentially expressed genes, indicated a strong enrichment for TFs in this set (p = 1.34E-05, Hypergeometric test). As the activity of TFs is governed by a number of factors including their protein levels, sub-cellular localization and post-translational modifications, the RNA levels may not always explain their activity. Our analysis and the observations presented here nevertheless provide a strong support that the 12 identified TFs could be significantly contributing to rewiring the host transcriptomic landscape under PPRV infection. Some of the 41 identified TFs likely controlling the expression of the differentially expressed genes may not always be differentially expressed because of our thresholds and/or they are farther upstream in the complex hierarchical transcriptional network  and hence their effect may not be as direct.
The quantity of data generated from RNA sequencing is large and therefore it is important to validate differential expression by real-time RT-PCR. Real time results in our study were found to be consistent with RNA sequencing analysis, but, we observed fold variations between two methods, which could be due to intrinsic differences between two techniques. Further, expression of in vitro validated genes (IFIT3, RB1CC1, VIM, CD44, CXCR4 and YY1) when evaluated for their expression in vaccinated goats (in vivo) showed similar expression pattern for IFIT3, RB1CC1 and VIM. However, CD44, CXCR4 and YY1 all of which were downregulated in in vitro were found be upregulated in vivo. This shows there could be a considerable difference in expression pattern of differentially expressed genes seen in vitro from in vivo experiments. Thus, it is necessary to analyze the expression of molecules which are predicted as anti-viral or pro-viral in, in vitro studies for their similar functions in vivo as in in vivo experiments the regulation of these molecules might be different due to a multitude of processes and interactions.
In this study, we investigated the host-virus interactome underlying PPRV vaccine virus infection. RNA-seq data analysis followed by functional analysis revealed enrichment of immune system regulatory pathways in PBMCs infected with PPRV. IFIT3 and TRIM56 can be suggested as the important early innate immune molecules following PPRV infection. The protein-protein interaction network of differentially expressed and highly connected genes revealed key dysregulated genes involved in immune system regulatory pathways, spliceosome pathways and apoptotic pathways. Genes encoding TFs that co-regulate the differentially expressed and highly connected genes were identified among the differentially expressed genes. Most of these TFs were found to govern immune regulatory pathways. Our results show that one can build the complete viral genome based on the transcriptome of the infected PBMCs indicating that the transcriptomes of blood as well as other body fluids could be used as markers to the know the extent of infection by various pathogens. This observation extends on recent studies on the impact of non-host genomic signals in the transcriptomes of humans and other hosts in uncovering the host-pathogen interplay .
Peste des petits ruminants
Peste des petits ruminants virus
Peripheral blood mononuclear cells
Tissue culture infective dose
Genomic mapping and alignment program
Fragments per Kilobase per Millions of reads
Database for s
Differentially expressed and highly connected genes
Multiple em and motif elicitation
Position weight matrices
This study was supported in part by Department of Biotechnology (BT/PR7729/AAQ/1/542/2013), Government of India. The authors would like to acknowledge Dr RP Singh for providing Sungri-96 vaccine virus for the study.
- Kerur N, Jhala MK, Joshi CG (2008) Genetic characterization of Indian peste des petits ruminants virus (PPRV) by sequencing and phylogenetic analysis of fusion protein and nucleoprotein gene segments. Res Vet Sci 85:176–183View ArticlePubMedGoogle Scholar
- Abu Elzein EM, Hassanien MM, Al-Afaleq AI, Abd Elhadi MA, Housawi FM (1990) Isolation of peste des petits ruminants from goats in Saudi Arabia. Vet Rec 127:309–310PubMedGoogle Scholar
- Senthil Kumar K, Babu A, Sundarapandian G, Roy P, Thangavelu A, Siva Kumar K, Arumugam R, Chandran ND, Muniraju M, Mahapatra M, Banyard AC, Manohar BM, Parida S (2014) Molecular characterisation of lineage IV peste des petits ruminants virus using multi gene sequence data. Vet Microbiol 174:39–49View ArticlePubMedGoogle Scholar
- Muniraju M, Mahapatra M, Ayelet G, Babu A, Olivier G, Munir M, Libeau G, Batten C, Banyard AC, Parida S. Emergence of lineage IV peste des petits ruminants virus in Ethiopia: complete genome sequence of an Ethiopian isolate 2010. Transbound Emerg Dis, in press
- Singh RP (2011) Control strategies for peste des petits ruminants in small ruminants of India. Rev Sci Tech 30:879–887PubMedGoogle Scholar
- Venkataraman R, Bandyopadhyay SK, Oberoi MS (2005) Present status and strategies for the control of transboundary and other economically important animal diseases in India: a review. Indian J Anim Sci 75:456–464Google Scholar
- Saravanan P, Sen A, Balamurugan V, Rajak KK, Bhanuprakash V, Palaniswami KS, Nachimuthu K, Thangavelu A, Dhinakarraj G, Hegde R, Singh RK (2010) Comparative efficacy of peste des petits ruminants (PPR) vaccines. Biologicals 38:479–485View ArticlePubMedGoogle Scholar
- Kumar N, Maherchandani S, Kashyap SK, Singh SV, Sharma S, Chaubey KK, Ly H (2014) Peste des petits ruminants virus infection of small ruminants: a comprehensive review. Viruses 6:2287–2327View ArticlePubMed CentralPubMedGoogle Scholar
- Bolt G, Berg K, Blixenkrone-Moller M (2002) Measles virus induced modulation of host-cell gene expression. J Gen Virol 83:1157–1165PubMedGoogle Scholar
- Iwasa T, Suga S, Qi L, Komada Y (2010) Apoptosis of human peripheral blood mononuclear cells by wild-type measles virus infection is induced by interaction of haemagglutinin and cellular receptor, SLAM via caspase-dependent pathway. Microbiol Immunol 54:405–416PubMedGoogle Scholar
- Pawar RM, Dhinakar Raj G, Balachandran C (2008) Relationship between the level of signaling lymphocyte activation molecule and replication of peste-des-petits-ruminants virus in peripheral blood mononuclear cells of host animals. Acta Virol 52:231–236PubMedGoogle Scholar
- Birch J, Juleff N, Heaton MP, Kalbfleisch T, Kijas J, Bailey D (2013) Characterization of ovine Nectin-4, a novel peste des petits ruminants virus receptor. J Virol 87:4756–4761View ArticlePubMed CentralPubMedGoogle Scholar
- Nanda SK, Baron J, Royall E, Robinson L, Falciani F, Baron MD (2009) In fection of bovine dendritic cells by rinderpest or measles virus induces different changes in host transcription. Virology 395:223–231View ArticlePubMedGoogle Scholar
- Schmeider R, Edwards R (2011) Quality control and preprocessing of metagenomic datasets. Bioinformatics 27:863–864View ArticleGoogle Scholar
- Liu H, Wang T, Wang J, Quan F, Zhang Y (2013) Characterization of liaoning cashmere goat transcriptome: sequencing, de novo assembly, functional annotation and comparative analysis. PLoS One 8:e77062View ArticlePubMed CentralPubMedGoogle Scholar
- Xu T, Guo X, Wang H, Du X, Gao X, Liu D (2013) De Novo transcriptome assembly and differential gene expression profiling of three Capra hircus skin types during anagen of the hair growth cycle. Int J Genomics 2013:269191View ArticlePubMed CentralPubMedGoogle Scholar
- Fontanesi L, Martelli PL, Beretti F, Riggio V, Dall’Olio S, Colombo M, Casadio R, Russo V, Portolano B (2010) An initial comparative map of copy number variations in the goat (Capra hircus) genome. BMC Genomics 11:639View ArticlePubMed CentralPubMedGoogle Scholar
- Wu TD, Watanabe CK (2005) GMAP: a genomic mapping and alignment program for mRNA and EST sequences. Bioinformatics 21:1859–1875View ArticlePubMedGoogle Scholar
- Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, Pimentel H, Salzberg SL, Rinn JL, Pachter L (2012) Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat Protoc 7:562–578View ArticlePubMed CentralPubMedGoogle Scholar
- Reimand J, Arak T, Vilo J (2011) g:Profiler–a web server for functional interpretation of gene lists (2011 update). Nucleic Acids Res 39:W307–W315View ArticlePubMed CentralPubMedGoogle Scholar
- da Huang W, Sherman BT, Lempicki RA (2009) Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc 4:44–57View ArticleGoogle Scholar
- Langmead B, Salzberg SL (2012) Fast gapped-read alignment with Bowtie 2. Nat Methods 9:357–359View ArticlePubMed CentralPubMedGoogle Scholar
- Stark C, Breitkreutz BJ, Chatr-Aryamontri A, Boucher L, Oughtred R, Livstone MS, Nixon J, Van Auken K, Wang X, Shi X, Reguly T, Rust JM, Winter A, Dolinski K, Tyers M (2011) The BioGRID interaction database: 2011 update. Nucleic Acids Res 39:D698–D704View ArticlePubMed CentralPubMedGoogle Scholar
- Yu H, Luscombe NM, Lu HX, Zhu X, Xia Y, Han JD, Bertin N, Chung S, Vidal M, Gerstein M (2004) Annotation transfer between genomes: protein-protein interologs and protein-DNA regulogs. Genome Res 14:1107–1118View ArticlePubMed CentralPubMedGoogle Scholar
- Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T (2003) Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res 13:2498–2504View ArticlePubMed CentralPubMedGoogle Scholar
- Csardi G, Nepusz T (2006) The igraph software package for complex network research. Inter Journal, Complex Systems 2006:1695Google Scholar
- Birney E, Andrews TD, Bevan P, Caccamo M, Chen Y, Clarke L, Coates G, Cuff J, Curwen V, Cutts T, Down T, Eyras E, Fernandez-Suarez XM, Gane P, Gibbins B, Gilbert J, Hammond M, Hotz HR, Iyer V, Jekosch K, Kahari A, Kasprzyk A, Keefe D, Keenan S, Lehvaslaiho H, McVicker G, Melsopp C, Meidl P, Mongin E, Pettett R, et al. (2004) An overview of Ensembl. Genome Res 14:925–928View ArticlePubMed CentralPubMedGoogle Scholar
- Bailey TL, Williams N, Misleh C, Li WW (2006) MEME: discovering and analyzing DNA and protein sequence motifs. Nucleic Acids Res 34:W369–W373View ArticlePubMed CentralPubMedGoogle Scholar
- Gupta S, Stamatoyannopoulos JA, Bailey TL, Noble WS (2007) Quantifying similarity between motifs. Genome Biol 8:R24View ArticlePubMed CentralPubMedGoogle Scholar
- Schmittgen TD, Livak KJ (2008) Analyzing real-time PCR data by the comparative C(T) method. Nat Protoc 3:1101–1108View ArticlePubMedGoogle Scholar
- Siddappa M, Gandham RK, Sarsani V, Mishra BP, Mishra B, Joshi CG, Sahoo AP, Tiwari AK, Janga SC (2014) Whole-Genome sequence of Sungri/96 vaccine starin of Peste des petits ruminants virus. Genome Announc 2:e0056–14View ArticleGoogle Scholar
- Chatterjee A, Chatterjee U, Ghosh MK (2013) Activation of protein kinase CK2 attenuates FOXO3a functioning in a PML-dependent manner: implications in human prostate cancer. Cell Death Dis 4:e543View ArticlePubMed CentralPubMedGoogle Scholar
- Mattera L, Courilleau C, Legube G, Ueda T, Fukunaga R, Chevillard-Briet M, Canitrot Y, Escaffit F, Trouche D (2010) The E1A-associated p400 protein modulates cell fate decisions by the regulation of ROS homeostasis. PLoS Genet 6:e1000983View ArticlePubMed CentralPubMedGoogle Scholar
- Tanji K, Tanaka T, Mori F, Kito K, Takahashi H, Wakabayashi K, Kamitani T (2006) NUB1 suppresses the formation of Lewy body-like inclusions by proteasomal degradation of synphilin-1. Am J Pathol 169:553–565View ArticlePubMed CentralPubMedGoogle Scholar
- Magnusson C, Norrby M, Libelius R, Tagerud S (2003) The nifk gene is widely expressed in mouse tissues and is up-regulated in denervated hind limb muscle. Cell Biol Int 27:469–475View ArticlePubMedGoogle Scholar
- Nourse CR, Mattei MG, Gunning P, Byrne JA (1998) Cloning of a third member of the D52 gene family indicates alternative coding sequence usage in D52-like transcripts. Biochim Biophys Acta 1443:155–168View ArticlePubMedGoogle Scholar
- Hao Y, Kong X, Ruan Y, Gan H, Chen H, Zhang C, Ren S, Gu J (2011) CDK11p46 and RPS8 associate with each other and suppress translation in a synergistic manner. Biochem Biophys Res Commun 407:169–174View ArticlePubMedGoogle Scholar
- Hsu LC, Huang X, Seasholtz S, Potter DM, Gollin SM (2006) Gene amplification and overexpression of protein phosphatase 1alpha in oral squamous cell carcinoma cell lines. Oncogene 25:5517–5526View ArticlePubMedGoogle Scholar
- Rangachari M, Zhu C, Sakuishi K, Xiao S, Karman J, Chen A, Angin M, Wakeham A, Greenfield EA, Sobel RA, Okada H, McKinnon PJ, Mak TW, Addo MM, Anderson AC, Kuchroo VK (2012) Bat3 promotes T cell responses and autoimmunity by repressing Tim-3-mediated cell death and exhaustion. Nat Med 18:1394–1400View ArticlePubMed CentralPubMedGoogle Scholar
- Golden T, Aragon IV, Rutland B, Tucker JA, Shevde LA, Samant RS, Zhou G, Amable L, Skarra D, Honkanen RE (2008) Elevated levels of Ser/Thr protein phosphatase 5 (PP5) in human breast cancer. Biochim Biophys Acta 1782:259–270View ArticlePubMed CentralPubMedGoogle Scholar
- Schmeisser H, Mejido J, Balinsky CA, Morrow AN, Clark CR, Zhao T, Zoon KC (2010) Identification of alpha interferon-induced genes associated with antiviral activity in Daudi cells and characterization of IFIT3 as a novel antiviral gene. J Virol 84:10671–10680View ArticlePubMed CentralPubMedGoogle Scholar
- Gil J, Busto EM, Garcillan B, Chean C, Garcia-Rodriguez MC, Diaz-Alderete A, Navarro J, Reine J, Mencia A, Gurbindo D, Beléndez C, Gordillo I, Duchniewicz M, Höhne K, García-Sánchez F, Fernández-Cruz E, López-Granados E, Schamel WW, Moreno-Pelayo MA, Recio MJ, Regueiro JR (2011) A leaky mutation in CD3D differentially affects alphabeta and gammadelta T cells and leads to a Talphabeta-Tgammadelta+B+NK+ human SCID. J Clin Invest 121:3872–3876View ArticlePubMed CentralPubMedGoogle Scholar
- Zhang H, Ma Q, Zhang Y, Xu H (2012) Proteolytic processing of Alzheimer’s β-amyloid precursor protein. J Neurochem 120:9–21View ArticlePubMed CentralPubMedGoogle Scholar
- Jiang N, Shen Y, Fei X, Sheng K, Sun P, Qiu Y, Larner J, Cao L, Kong X, Mi J (2013) Valosin-containing protein regulates the proteasome-mediated degradation of DNA-PKcs in glioma cells. Cell Death Dis 4:e647View ArticlePubMed CentralPubMedGoogle Scholar
- Glisovic T, Bachorik JL, Yong J, Dreyfuss G (2008) RNA-binding proteins and post-transcriptional gene regulation. FEBS Lett 582:1977–1986View ArticlePubMed CentralPubMedGoogle Scholar
- Wiegand MA, Bossow S, Schlecht S, Neubert WJ (2007) De novo synthesis of N and P proteins as a key step in Sendai virus gene expression. J Virol 81:13835–13844View ArticlePubMed CentralPubMedGoogle Scholar
- Hegde R, Gomes AR, Byregowda SM, Santhosh AK, Renukaprasad C (2009) Cytopathic effect of PPR vaccine virus strains in Vero cells. Vet World 2:93–94Google Scholar
- Yunus M, Shaila MS (2012) Establishment of an in vitro transcription system for Peste des petits ruminant virus. Virol J 9:302View ArticlePubMed CentralPubMedGoogle Scholar
- Wilkinson B, Downey JS, Rudd CE (2005) T-cell signalling and immune system disorders. Expert Rev Mol Med 7:1–29View ArticlePubMedGoogle Scholar
- Esche C, Stellato C, Beck LA (2005) Chemokines: key players in innate and adaptive immunity. J Invest Dermatol 125:615–628View ArticlePubMedGoogle Scholar
- Myers CD (1991) Role of B cell antigen processing and presentation in the humoral immune response. FASEB J 5:2547–2553PubMedGoogle Scholar
- Shuai K, Liu B (2003) Regulation of JAK-STAT signalling in the immune system. Nat Rev Immunol 3:900–911View ArticlePubMedGoogle Scholar
- Katzman SD, Hoyer KK, Dooms H, Gratz IK, Rosenblum MD, Paw JS, Isakson SH, Abbas AK (2011) Opposing functions of IL-2 and IL-7 in the regulation of immune responses. Cytokine 56:116–121View ArticlePubMed CentralPubMedGoogle Scholar
- Wang J, Liu B, Wang N, Lee YM, Liu C, Li K (2011) TRIM56 is a virus- and interferon-inducible E3 ubiquitin ligase that restricts pestivi- rus infection. J Virol 85:3733–3745View ArticlePubMed CentralPubMedGoogle Scholar
- Shen Y, Li NL, Wang J, Liu B, Lester S, Li K (2012) TRIM56 is an essential component of TLR3 antiviral signalling pathway. J Biol Chem 287:36404–36413View ArticlePubMed CentralPubMedGoogle Scholar
- Morrison J, Josset L, Tchitchek N, Chang J, Belser JA, Swayne DE, Pantin-Jackwood MJ, Tumpey TM, Katze MG (2014) H7N9 and other pathogenic avian influenza viruses elicit a three-pronged transcriptomic signature that is reminiscent of 1918 influenza virus and is associated with lethal outcome in mice. J Virol 88:10556–10568View ArticlePubMedGoogle Scholar
- Wachi S, Yoneda K, Wu R (2005) Interactome-transcriptome analysis reveals the high centrality of genes differentially expressed in lung cancer tissues. Bioinformatics 21:4205–4208View ArticlePubMedGoogle Scholar
- Mondal B, Sreenivasa BP, Dhar P, Singh RP, Bandyopadhyay SK (2001) Apoptosis induced by peste des petits ruminants virus in goat peripheral blood mononuclear cells. Virus Res 73:113–119View ArticlePubMedGoogle Scholar
- Mogensen TH, Melchjorsen J, Malmgaard L, Casola A, Paludan SR (2004) Suppression of proinflammatory cytokine expression by herpes simplex virus type 1. J Virol 78:5883–5890View ArticlePubMed CentralPubMedGoogle Scholar
- Arita M, Wakita T, Shimizu H (2012) Valosin-containing protein (VCP/p97) is required for poliovirus replication and is involved in cellular protein secretion pathway in poliovirus infection. J Virol 86:5541–5553View ArticlePubMed CentralPubMedGoogle Scholar
- Inoue K, Tsukiyama-Kohara K, Matsuda C, Yoneyama M, Fujita T, Kuge S, Yoshiba M, Kohara M (2012) Impairment of interferon regulatory factor-3 activation by hepatitis C virus core protein basic amino acid region 1. Biochem Biophys Res Commun 428:494–499View ArticlePubMedGoogle Scholar
- Forero A, Moore PS, Sarkar SN (2013) Role of IRF4 in IFN-stimulated gene induction and maintenance of Kaposi sarcoma-associated herpesvirus latency in primary effusion lymphoma cells. J Immunol 191:1476–1485View ArticlePubMed CentralPubMedGoogle Scholar
- Sato K, Miyoshi F, Yokota K, Araki Y, Asanuma Y, Akiyama Y, Yoh K, Takahashi S, Aburatani H, Mimura T (2011) Marked induction of c-Maf protein during Th17 cell differentiation and its implication in memory Th cell development. J Biol Chem 286:14963–14971View ArticlePubMed CentralPubMedGoogle Scholar
- Oh H, Ghosh S (2013) NF-kappaB: roles and regulation in different CD4(+) T-cell subsets. Immunol Rev 252:41–51View ArticlePubMed CentralPubMedGoogle Scholar
- Lin MH, Chou FC, Yeh LT, Fu SH, Chiou HY, Lin KI, Chang DM, Sytwu HK (2013) B lymphocyte-induced maturation protein 1 (BLIMP-1) attenuates autoimmune diabetes in NOD mice by suppressing Th1 and Th17 cells. Diabetologia 56:136–146View ArticlePubMedGoogle Scholar
- Villa C, Ridolfi E, Fenoglio C, Ghezzi L, Vimercati R, Clerici F, Marcone A, Gallone S, Serpente M, Cantoni C, Bonsi R, Cioffi S, Cappa S, Franceschi M, Rainero I, Mariani C, Scarpini E, Galimberti D (2013) Expression of the transcription factor Sp1 and its regulatory hsa-miR-29b in peripheral blood mononuclear cells from patients with Alzheimer’s disease. J Alzheimers Dis 35:487–494PubMedGoogle Scholar
- Chen NX, Kiattisunthorn K, O’Neill KD, Chen X, Moorthi RN, Gattone VH, 2nd, Allen MR, Moe SM (2013) Decreased microRNA is involved in the vascular remodeling abnormalities in chronic kidney disease (CKD). PLoS One 8:e64558View ArticlePubMed CentralPubMedGoogle Scholar
- Gerstein MB, Kundaje A, Hariharan M, Landt SG, Yan KK, Cheng C, Mu XJ, Khurana E, Rozowsky J, Alexander R, Min R, Alves P, Abyzov A, Addleman N, Bhardwaj N, Boyle AP, Cayting P, Charos A, Chen DZ, Cheng Y, Clarke D, Eastman C, Euskirchen G, Frietze S, Fu Y, Gertz J, Grubert F, Harmanci A, Jain P, Kasowski M, et al. (2012) Architecture of the human regulatory network derived from ENCODE data. Nature 489:91–100View ArticlePubMed CentralPubMedGoogle Scholar
- Wang K, Li H, Yuan Y, Etheridge A, Zhou Y, Huang D, Wilmes P, Galas D (2012) The complex exogenous RNA spectra in human plasma: an interface with human gut biota? PLoS One 7:e51009View ArticlePubMed CentralPubMedGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.