Palanikumar Ravindran

Expression profiling of human immune cell subsets identifies miRNA-mRNA regulatory relationships correlated with cell type specific expression

Blood consists of different cell populations with distinct functions and correspondingly, distinct gene expression profiles. In this study, global miRNA expression profiling was performed across a panel of nine human immune cell subsets (neutrophils, eosinophils, monocytes, B cells, NK cells, CD4 T cells, CD8 T cells, mDCs and pDCs) to identify cell-type specific miRNAs. mRNA expression profiling was performed on the same samples to determine if miRNAs specific to certain cell types down-regulated expression levels of their target genes. Six cell-type specific miRNAs (miR-143; neutrophil specific, miR-125; T cells and neutrophil specific, miR-500; monocyte and pDC specific, miR-150; lymphoid cell specific, miR-652 and miR-223; both myeloid cell specific) were negatively correlated with expression of their predicted target genes. These results were further validated using an independent cohort where similar immune cell subsets were isolated and profiled for both miRNA and mRNA expression. miRNAs which negatively correlated with target gene expression in both cohorts were identified as candidates for miRNA/mRNA regulatory pairs and were used to construct a cell-type specific regulatory network. miRNA/mRNA pairs formed two distinct clusters in the network corresponding to myeloid (nine miRNAs) and lymphoid lineages (two miRNAs). Several myeloid specific miRNAs targeted common genes including ABL2, EIF4A2, EPC1 and INO80D; these common targets were enriched for genes involved in the regulation of gene expression (p<9.0E-7). Those miRNA might therefore have significant further effect on gene expression by repressing the expression of genes involved in transcriptional regulation. The miRNA and mRNA expression profiles reported in this study form a comprehensive transcriptome database of various human blood cells and serve as a valuable resource for elucidating the role of miRNA mediated regulation in the establishment of immune cell identity.

Bleomycin induces molecular changes directly relevant to idiopathic pulmonary fibrosis: a model for "active" disease.

The preclinical model of bleomycin-induced lung fibrosis, used to investigate mechanisms related to idiopathic pulmonary fibrosis (IPF), has incorrectly predicted efficacy for several candidate compounds suggesting that it may be of limited value. As an attempt to improve the predictive nature of this model, integrative bioinformatic approaches were used to compare molecular alterations in the lungs of bleomycin-treated mice and patients with IPF. Using gene set enrichment analysis we show for the first time that genes differentially expressed during the fibrotic phase of the single challenge bleomycin model were significantly enriched in the expression profiles of IPF patients. The genes that contributed most to the enrichment were largely involved in mitosis, growth factor, and matrix signaling. Interestingly, these same mitotic processes were increased in the expression profiles of fibroblasts isolated from rapidly progressing, but not slowly progressing, IPF patients relative to control subjects. The data also indicated that TGFβ was not the sole mediator responsible for the changes observed in this model since the ALK-5 inhibitor SB525334 effectively attenuated some but not all of the fibrosis associated with this model. Although some would suggest that repetitive bleomycin injuries may more effectively model IPF-like changes, our data do not support this conclusion. Together, these data highlight that a single bleomycin instillation effectively replicates several of the specific pathogenic molecular changes associated with IPF, and may be best used as a model for patients with active disease.

Systemic Biomarkers of Neutrophilic Inflammation, Tissue Injury and Repair in COPD Patients with Differing Levels of Disease Severity

The identification and validation of biomarkers to support the assessment of novel therapeutics for COPD continues to be an important area of research. The aim of the current study was to identify systemic protein biomarkers correlated with measures of COPD severity, as well as specific protein signatures associated with comorbidities such as metabolic syndrome. 142 protein analytes were measured in serum of 140 patients with stable COPD, 15 smokers without COPD and 30 non-smoking controls. Seven analytes (sRAGE, EN-RAGE, NGAL, Fibrinogen, MPO, TGF-α and HB-EGF) showed significant differences between severe/very severe COPD, mild/moderate COPD, smoking and non-smoking control groups. Within the COPD subjects, univariate and multivariate analyses identified analytes significantly associated with FEV1, FEV1/FVC and DLCO. Most notably, a set of 5 analytes (HB-EGF, Fibrinogen, MCP-4, sRAGE and Sortilin) predicted 21% of the variability in DLCO values. To determine common functions/pathways, analytes were clustered in a correlation network by similarity of expression profile. While analytes related to neutrophil function (EN-RAGE, NGAL, MPO) grouped together to form a cluster associated with FEV1 related parameters, analytes related to the EGFR pathway (HB-EGF, TGF-α) formed another cluster associated with both DLCO and FEV1 related parameters. Associations of Fibrinogen with DLCO and MPO with FEV1/FVC were stronger in patients without metabolic syndrome (r  =  −0.52, p  = 0.005 and r  =  −0.61, p  = 0.023, respectively) compared to patients with coexisting metabolic syndrome (r  =  −0.25, p  = 0.47 and r  =  −0.15, p  = 0.96, respectively), and may be driving overall associations in the general cohort. In summary, our study has identified known and novel serum protein biomarkers and has demonstrated specific associations with COPD disease severity, FEV1, FEV1/FVC and DLCO. These data highlight systemic inflammatory pathways, neutrophil activation and epithelial tissue injury/repair processes as key pathways associated with COPD.

Transcriptomic analysis of the woodchuck model of chronic hepatitis B

The Eastern woodchuck (Marmota monax) is naturally infected with woodchuck hepatitis virus (WHV), a hepadnavirus closely related to the human hepatitis B virus (HBV). The woodchuck is used as an animal model for studying chronic hepatitis B (CHB) and HBV‐associated hepatocellular carcinoma (HCC) in humans, but the lack of sequence information has hitherto precluded functional genomics analysis. To address this major limitation of the model, we report here the sequencing, assembly, and annotation of the woodchuck transcriptome, together with the generation of custom woodchuck microarrays. Using this new platform, we characterized the transcriptional response to persistent WHV infection and WHV‐induced HCC. This revealed that chronic WHV infection, like HBV, is associated with (1) a limited intrahepatic type I interferon response; (2) intrahepatic induction of markers associated with T cell exhaustion; (3) elevated levels of suppressor of cytokine signaling 3 (SOCS3) in the liver; and (4) intrahepatic accumulation of neutrophils. Underscoring the translational value of the woodchuck model, this study also determined that WHV‐induced HCC shares molecular characteristics with a subtype of human HCC with poor prognosis. Conclusion: Our data establish the translational value of the woodchuck model and provide new insight into immune pathways which may play a role either in the persistence of HBV infection or the sequelae of CHB. (HEPATOLOGY 2012;56:820–830)

Determination of the Human Cardiomyocyte mRNA and miRNA Differentiation Network by Fine-Scale Profiling

To gain insight into the molecular regulation of human heart development, a detailed comparison of the mRNA and miRNA transcriptomes across differentiating human-induced pluripotent stem cell (hiPSC)-derived cardiomyocytes and biopsies from fetal, adult, and hypertensive human hearts was performed. Gene ontology analysis of the mRNA expression levels of the hiPSCs differentiating into cardiomyocytes revealed 3 distinct groups of genes: pluripotent specific, transitional cardiac specification, and mature cardiomyocyte specific. Hierarchical clustering of the mRNA data revealed that the transcriptome of hiPSC cardiomyocytes largely stabilizes 20 days after initiation of differentiation. Nevertheless, analysis of cells continuously cultured for 120 days indicated that the cardiomyocytes continued to mature toward a more adult-like gene expression pattern. Analysis of cardiomyocyte-specific miRNAs (miR-1, miR-133a/b, and miR-208a/b) revealed an miRNA pattern indicative of stem cell to cardiomyocyte specification. A biostatistitical approach integrated the miRNA and mRNA expression profiles revealing a cardiomyocyte differentiation miRNA network and identified putative mRNAs targeted by multiple miRNAs. Together, these data reveal the miRNA network in human heart development and support the notion that overlapping miRNA networks re-enforce transcriptional control during developmental specification.

Characterization of a novel CRAC inhibitor that potently blocks human T cell activation and effector functions

Store operated calcium entry (SOCE) downstream of T cell receptor (TCR) activation in T lymphocytes has been shown to be mediated mainly through the Calcium Release Activated Calcium (CRAC) channel. Here, we compared the effects of a novel, potent and selective CRAC current inhibitor, 2,6-Difluoro-N-{5-[4-methyl-1-(5-methyl-thiazol-2-yl)-1,2,5,6-tetrahydro-pyridin-3-yl]-pyrazin-2-yl}-benzamide (RO2959), on T cell effector functions with that of a previously reported CRAC channel inhibitor, YM-58483, and a calcineurin inhibitor Cyclosporin A (CsA). Using both electrophysiological and calcium-based fluorescence measurements, we showed that RO2959 is a potent SOCE inhibitor that blocked an IP3-dependent current in CRAC-expressing RBL-2H3 cells and CHO cells stably expressing human Orai1 and Stim1, as well as SOCE in human primary CD4(+) T cells triggered by either TCR stimulation or thapsigargin treatment. Furthermore, we demonstrated that RO2959 completely inhibited cytokine production as well as T cell proliferation mediated by TCR stimulation or MLR (mixed lymphocyte reaction). Lastly, we showed by gene expression array analysis that RO2959 potently blocked TCR triggered gene expression and T cell functional pathways similar to CsA and another calcineurin inhibitor FK506. Thus, both from a functional and transcriptional level, our data provide evidence that RO2959 is a novel and selective CRAC current inhibitor that potently inhibits human T cell functions.

Bayesian Analysis of Circular Data Using Wrapped Distributions

Circular data arise in a number of different areas such as geological, meteorological, biological and industrial sciences. Standard statistical techniques can not be used to model circular data due to the circular geometry of the sample space. One of the common methods to analyze circular data is known as the wrapping approach. This approach is based on a simple fact that a probability distribution on a circle can be obtained by wrapping a probability distribution defined on the real line. A large class of probability distributions that are flexible to account for different features of circular data can be obtained by the aforementioned approach. However, the likelihood-based inference for wrapped distributions can be very complicated and computationally intensive. The EM algorithm to compute the MLE is feasible, but is computationally unsatisfactory. A data augmentation method using slice sampling is proposed to overcome such computational difficulties. The proposed method turns out to be flexible and computationally efficient to fit a wide class of wrapped distributions. In addition, a new model selection criteria for circular data is developed. Results from an extensive simulation study are presented to validate the performance of the proposed estimation method and the model selection criteria. Application to a real data set is also presented and parameter estimates are compared to those that are available in the literature.

Identification of an intrahepatic transcriptional signature associated with self‐limiting infection in the woodchuck model of hepatitis B

The woodchuck model of hepatitis B virus (HBV) infection displays many characteristics of human infection and has particular value for characterizing the host immune responses during the development of chronic infection. Using the newly developed custom woodchuck microarray platform, we compared the intrahepatic transcriptional profiles of neonatal woodchucks with self‐limiting woodchuck hepatitis virus (WHV) infection to those woodchucks progressing to persistent WHV infection. This revealed that WHV does not induce significant intrahepatic gene expression changes during the early‐acute stage of infection (8 weeks), suggesting it is a stealth virus. At the mid‐acute phase of infection (14 weeks), resolution was associated with induction of a prominent cytotoxic T‐cell signature. Strikingly, this was accompanied by high‐level expression of PD‐1 and various other inhibitory T‐cell receptors, which likely act to minimize liver damage by cytotoxic T cells during viral clearance. In contrast to the expression of perforin and other cytotoxic effector genes, the interferon‐γ (IFN‐γ) signaling response in the mid‐acute phase was comparable to that in chronically infected adult animals. The absence of a strong IFN‐α/β transcriptional response indicated that type I IFN is not a critical mediator of self‐limiting infection. Nevertheless, a number of antiviral genes, including viperin, were differentially expressed during resolving infection, suggesting that a subset of IFN‐stimulated genes (ISG) may play a role in the control of WHV replication. Conclusion: We identified new immune pathways associated with the clearance of hepadnavirus infection revealing novel molecular targets with potential for the therapeutic treatment of chronic hepatitis B. (HEPATOLOGY 2013)

Bruton's tyrosine kinase regulates TLR9 but not TLR7 signaling in human plasmacytoid dendritic cells

Plasmacytoid dendritic cells (PDCs) represent a key cell type for both innate and adaptive immunity. PDCs express both TLR7 and TLR9 and the recognition of nucleic acids by these two receptors triggers the production of a large amount of type‐I IFN and the induction of PDC maturation into APCs. This unique feature of PDCs is at the basis of clinical development of both TLR7 and TLR9 agonists for infectious diseases, allergy, cancer, and asthma. However, TLR7 and TLR9 recognition of self‐nucleic acids is linked to many autoimmune diseases including lupus, and a better understanding of the signaling pathways of these two receptors in PDCs is thus important. We have identified Brutons tyrosine kinase (Btk) as an important player for TLR9 but not TLR7 signaling in human PDCs. Blocking Btk using a specific inhibitor leads to the reduction of all TLR9‐induced responses in PDCs, including cytokine production and expression of costimulatory molecules, while this has no impact on the TLR7 response. This identifies Btk as a key molecule in TLR9 signaling in PDCs and is the first demonstration that the TLR7 and TLR9 pathways can be dissociated in human PDCs.

Double-stranded RNA induces molecular and inflammatory signatures that are directly relevant to COPD

Polyinosinic:polycytidylic acid (poly I:C) is a synthetic analogue of double-stranded (ds)RNA, a molecular pattern associated with viral infections, that is used to exacerbate inflammation in lung injury models. Despite its frequent use, there are no detailed studies of the responses elicited by a single topical administration of poly I:C to the lungs of mice. Our data provides the first demonstration that the molecular responses in the airways induced by poly I:C correlate to those observed in the lungs of chronic obstructive pulmonary disease (COPD) patients. These expression data also revealed three distinct phases of response to poly I:C, consistent with the changing inflammatory cell infiltrate in the airways. Poly I:C induced increased numbers of neutrophils and natural killer cells in the airways, which were blocked by CXCR2 and CCR5 antagonists, respectively. Using gene set variation analysis on representative clinical data sets, gene sets defined by poly I:C–induced differentially expressed genes were enriched in the molecular profiles of COPD but not idiopathic pulmonary fibrosis patients. Collectively, these data represent a new approach for validating the clinical relevance of preclinical animal models and demonstrate that a dual CXCR2/CCR5 antagonist may be an effective treatment for COPD patients.