Search

Arthritis is a common manifestation of systemic lupus erythematosus (SLE) yet understanding of the underlying pathogenic mechanisms remains incomplete. We, therefore, interrogated gene expression profiles of SLE synovium to gain insight into the nature of lupus arthritis (LA), using osteoarthritis (OA) and rheumatoid arthritis (RA) as comparators. Knee synovia from SLE, OA, and RA patients were analyzed for differentially expressed genes (DEGs) and also by Weighted Gene Co-expression Network Analysis (WGCNA) to identify modules of highly co-expressed genes. Genes upregulated and/or co-expressed in LA revealed numerous immune/inflammatory cells dominated by a myeloid phenotype, in which pathogenic macrophages, myeloid-lineage cells, and their secreted products perpetuate inflammation, whereas OA was characterized by fibroblasts and RA of lymphocytes. Genes governing trafficking of immune cells into the synovium by chemokines were identified, but not in situ generation of germinal centers (GCs). Gene Set Variation Analysis (GSVA) confirmed activation of specific immune cell types in LA. Numerous therapies were predicted to target LA, including TNF, NFκB, MAPK, and CDK inhibitors. Detailed gene expression analysis identified a unique pattern of cellular components and physiologic pathways operative in LA, as well as drugs potentially able to target this common manifestation of SLE.

Gene expression signatures can stratify patients with heterogeneous diseases, such as systemic lupus erythematosus (SLE), yet understanding the contributions of ancestral background to this heterogeneity is not well understood. We hypothesized that ancestry would significantly influence gene expression signatures and measured 34 gene modules in 1566 SLE patients of African ancestry (AA), European ancestry (EA), or Native American ancestry (NAA). Healthy subject ancestry-specific gene expression provided the transcriptomic background upon which the SLE patient signatures were built. Although standard therapy affected every gene signature and significantly increased myeloid cell signatures, logistic regression analysis determined that ancestral background significantly changed 23 of 34 gene signatures. Additionally, the strongest association to gene expression changes was found with autoantibodies, and this also had etiology in ancestry: the AA predisposition to have both RNP and dsDNA autoantibodies compared with EA predisposition to have only anti-dsDNA. A machine learning approach was used to determine a gene signature characteristic to distinguish AA SLE and was most influenced by genes characteristic of the perturbed B cell axis in AA SLE patients.

The integration of gene expression data to predict systemic lupus erythematosus (SLE) disease activity is a significant challenge because of the high degree of heterogeneity among patients and study cohorts, especially those collected on different microarray platforms. Here we deployed machine learning approaches to integrate gene expression data from three SLE data sets and used it to classify patients as having active or inactive disease as characterized by standard clinical composite outcome measures. Both raw whole blood gene expression data and informative gene modules generated by Weighted Gene Co-expression Network Analysis from purified leukocyte populations were employed with various classification algorithms. Classifiers were evaluated by 10-fold cross-validation across three combined data sets or by training and testing in independent data sets, the latter of which amplified the effects of technical variation. A random forest classifier achieved a peak classification accuracy of 83 percent under 10-fold cross-validation, but its performance could be severely affected by technical variation among data sets. The use of gene modules rather than raw gene expression was more robust, achieving classification accuracies of approximately 70 percent regardless of how the training and testing sets were formed. Fine-tuning the algorithms and parameter sets may generate sufficient accuracy to be informative as a standalone estimate of disease activity.

Systemic lupus erythematosus (SLE) is an autoimmune disease characterized by the presence of low-density granulocytes (LDGs) with a heightened capacity for spontaneous NETosis, but the contribution of LDGs to SLE pathogenesis remains unclear. To characterize LDGs in human SLE, gene expression profiles derived from isolated LDGs were characterized by weighted gene coexpression network analysis, and a 92-gene module was identified. The LDG gene signature was enriched in genes related to neutrophil degranulation and cell cycle regulation. This signature was assessed in gene expression datasets from two large-scale SLE clinical trials to study associations between LDG enrichment, SLE manifestations, and treatment regimens. LDG enrichment in the blood was associated with corticosteroid treatment as well as anti-dsDNA, low serum complement, renal manifestations, and vasculitis, but the latter two of these associations were dependent on concomitant corticosteroid treatment. In addition, LDG enrichment was associated with enrichment of gene signatures induced by type I IFN and TNF irrespective of corticosteroid treatment. Notably, LDG enrichment was not found in numerous tissues affected by SLE. Comparison with relevant reference datasets indicated that LDG enrichment is likely reflective of increased granulopoiesis in the bone marrow and not peripheral neutrophil activation. The results have uncovered important determinants of the appearance of LDGs in SLE and have emphasized the likely role of LDGs in specific aspects of lupus pathogenesis., (Copyright © 2019 by The American Association of Immunologists, Inc.)

A role for interferon (IFN) in systemic lupus erythematosus (SLE) pathogenesis is inferred from the prominent IFN gene signature (IGS), but the major IFN species and its relationship to disease activity are unknown. A bioinformatic approach employing individual IFN species gene signatures to interrogate SLE microarray datasets demonstrates a putative role for numerous IFN species, with prominent expression of IFNB1 and IFNW signatures. In contrast with other SLE-affected organs, the IGS is less prominent in lupus nephritis. SLE patients with active and inactive disease have readily detectable IGS and the IGS changes synchronously with a monocyte signature but not disease activity, and is significantly related to monocyte transcripts. Monocyte over-expression of three times as many IGS transcripts as T and B cells and IGS retention in monocytes, but not T and B cells from inactive SLE patients contribute to the lack of correlation between the IGS and SLE disease activity., Competing Interests: The authors declare no competing interests.

Systemic lupus erythematosus (SLE) is characterized by abnormalities in B cell and T cell function, but the role of disturbances in the activation status of macrophages (Mϕ) has not been well described in human patients. To address this, gene expression profiles from isolated lymphoid and myeloid populations were analyzed to identify differentially expressed (DE) genes between healthy controls and patients with either inactive or active SLE. While hundreds of DE genes were identified in B and T cells of active SLE patients, there were no DE genes found in B or T cells from patients with inactive SLE compared to healthy controls. In contrast, large numbers of DE genes were found in myeloid cells (MC) from both active and inactive SLE patients. Among the DE genes were several known to play roles in Mϕ activation and polarization, including the M1 genes STAT1 and SOCS3 and the M2 genes STAT3, STAT6, and CD163. M1-associated genes were far more frequent in data sets from active versus inactive SLE patients. To characterize the relationship between Mϕ activation and disease activity in greater detail, weighted gene co-expression network analysis (WGCNA) was used to identify modules of genes associated with clinical activity in SLE patients. Among these were disease activity-correlated modules containing activation signatures of predominantly M1-associated genes. No disease activity-correlated modules were enriched in M2-associated genes. Pathway and upstream regulator analysis of DE genes from both active and inactive SLE MC were cross-referenced with high-scoring hits from the drug discovery Library of Integrated Network-based Cellular Signatures (LINCS) to identify new strategies to treat both stages of SLE. A machine learning approach employing MC gene modules and a generalized linear model was able to predict the disease activity status in unrelated gene expression data sets. In summary, altered MC gene expression is characteristic of both active and inactive SLE. However, disease activity is associated with an alteration in the activation of MC, with a bias toward the M1 proinflammatory phenotype. These data suggest that while hyperactivity of B cells and T cells is associated with active SLE, MC potentially direct flare-ups and remission by altering their activation status toward the M1 state., Competing Interests: All authors are employed by AMPEL BioSolutions, LLC or by the RILITE Research Institute. This does not alter our adherence to PLoS One policies on data sharing and materials. The work presented in this study is not part of any patents or marketed products.

Background: Leishmania major infection induces robust interleukin-12 (IL12) production in human dendritic cells (hDC), ultimately resulting in Th1-mediated immunity and clinical resolution. The surface of Leishmania parasites is covered in a dense glycocalyx consisting of primarily lipophosphoglycan (LPG) and other phosphoglycan-containing molecules (PGs), making these glycoconjugates the likely pathogen-associated molecular patterns (PAMPS) responsible for IL12 induction., Methodology/principal Findings: Here we explored the role of parasite glycoconjugates on the hDC IL12 response by generating L. major Friedlin V1 mutants defective in LPG alone, (FV1 lpg1-), or generally deficient for all PGs, (FV1 lpg2-). Infection with metacyclic, infective stage, L. major or purified LPG induced high levels of IL12B subunit gene transcripts in hDCs, which was abrogated with FV1 lpg1- infections. In contrast, hDC infections with FV1 lpg2- displayed increased IL12B expression, suggesting other PG-related/LPG2 dependent molecules may act to dampen the immune response. Global transcriptional profiling comparing WT, FV1 lpg1-, FV1 lpg2- infections revealed that FV1 lpg1- mutants entered hDCs in a silent fashion as indicated by repression of gene expression. Transcription factor binding site analysis suggests that LPG recognition by hDCs induces IL-12 in a signaling cascade resulting in Nuclear Factor κ B (NFκB) and Interferon Regulatory Factor (IRF) mediated transcription., Conclusions/significance: These data suggest that L. major LPG is a major PAMP recognized by hDC to induce IL12-mediated protective immunity and that there is a complex interplay between PG-baring Leishmania surface glycoconjugates that result in modulation of host cellular IL12.

Leishmania are intracellular protozoa that influence host immune responses eliciting parasite species-specific pathologies. MicroRNAs (miRNAs) are short single-stranded ribonucleic acids that complement gene transcripts to block protein translation and have been shown to regulate immune system molecular mechanisms. Human monocyte-derived dendritic cells (DC) and macrophages (MP) were infected in vitro with Leishmania major or Leishmania donovani parasites. Small RNAs were isolated from total RNA and sequenced to identify mature miRNAs associated with leishmanial infections. Normalized sequence read count profiles revealed a global downregulation in miRNA expression among host cells following infection. Most identified miRNAs were expressed at higher levels in L. donovani-infected cells relative to L. major-infected cells. Pathway enrichments using in silico-predicted gene targets of differentially expressed miRNAs showed evidence of potentially universal MAP kinase signalling pathway effects. Whereas JAK-STAT and TGF-β signalling pathways were more highly enriched using targets of miRNAs upregulated in L. donovani-infected cells, these data provide evidence in support of a selective influence on host cell miRNA expression and regulation in response to differential Leishmania infections., (© 2014 John Wiley & Sons Ltd.)

Leishmania major-infected human dendritic cells (DCs) exhibit a marked induction of IL-12, ultimately promoting a robust Th1-mediated response associated with parasite killing and protective immunity. The host cell transcription machinery associated with the specific IL-12 induction observed during L. major infection remains to be thoroughly elucidated. In this study, we used Affymetrix GeneChip (Affymetrix) to globally assess the host cell genes and pathways associated with early L. major infection in human myeloid-derived DCs. Our data revealed 728 genes were significantly differentially expressed and molecular signaling pathway revealed that the type I IFN pathway was significantly enriched. Addition of a neutralizing type I IFN decoy receptor blocked the expression of IRF7 and IL-12p40 during DC infection, indicating the L. major-induced expression of IL-12p40 is dependent upon the type I IFN signaling pathway. In stark contrast, IL-12p40 expression is not elicited by L. donovani, the etiological agent of deadly visceral leishmaniasis. Therefore, we examined the gene expression profile for several IFN response genes in L. major versus L. donovani DC infections. Our data revealed that L. major, but not L. donovani, induces expression of IRF2, IRF7, and IFIT5, implicating the regulation of type I IFN-associated signaling pathways as mediating factors toward the production of IL-12., (Copyright © 2014 by The American Association of Immunologists, Inc.)

Phlebotomus papatasi sand flies are among the primary vectors of Leishmania major parasites from Morocco to the Indian subcontinent and from southern Europe to central and eastern Africa. Antibody-based immunity to sand fly salivary gland proteins in human populations remains a complex contextual problem that is not yet fully understood. We profiled the immunoreactivities of plasma antibodies to sand fly salivary gland sonicates (SGSs) from 229 human blood donors residing in different regions of sand fly endemicity throughout Jordan and Egypt as well as 69 US military personnel, who were differentially exposed to P. papatasi bites and L. major infections in Iraq. Compared with plasma from control region donors, antibodies were significantly immunoreactive to five salivary proteins (12, 26, 30, 38, and 44 kDa) among Jordanian and Egyptian donors, with immunoglobulin G4 being the dominant anti-SGS isotype. US personnel were significantly immunoreactive to only two salivary proteins (38 and 14 kDa). Using k-means clustering, donors were segregated into four clusters distinguished by unique immunoreactivity profiles to varying combinations of the significantly immunogenic salivary proteins. SGS-induced cellular proliferation was diminished among donors residing in sand fly-endemic regions. These data provide a clearer picture of human immune responses to sand fly vector salivary constituents.

Dendritic cells (DC) play a pivotal role in regulating immunity, establishing immunologically privileged tissue microenvironments and maintaining homoeostasis. It is becoming increasingly clear that one key mechanism that mediates many DC functions is production of the immunomodulatory enzyme indoleamine 2,3-dioxygenase (IDO). For pathogens that cause chronic infection, exploitation of host DCs is a solution to establish and persist within a host. Leishmania parasites cause a range of clinical manifestations, all involving chronic infection, and are proficient at avoiding immune responses. We demonstrate here that infection of human myeloid-derived DC with L. major and L. donovani induces IDO expression using a mechanism that involves autocrine or paracrine stimulation with a DC-secreted factor. Leishmania-induced IDO suppresses allogeneic and tetanus toxoid-specific lymphocyte proliferation, an inhibition that is reversed with the IDO inhibitor, 1-methyl tryptophan (1-MT). Furthermore, IDO expression by human DC does not require live Leishmania infection, as parasite lysates also up-regulate IDO mRNA production. Our data suggest that one mechanism Leishmania parasites utilize to circumvent immune clearance may be to promote the induction of IDO among host DC within the infection microenvironment., (© 2012 Blackwell Publishing Ltd.)

Fluorescent in-situ hybridization (FISH) using meiotic chromosome preparations and highly repetitive DNA from the southern cattle tick, Rhipicephalus microplus, was undertaken to investigate genome organization. Several classes of highly repetitive DNA elements were identified by screening a R. microplus bacterial artificial chromosome (BAC) library. A repeat unit of approximately 149 bp, RMR-1 was localized to the subtelomeric regions of R. microplus autosomes 1-6 and 8-10. A second repeat unit, RMR-2 was localized to the subtelomeric regions of all autosomes and the X chromosome. RMR-2 was composed of three distinct repeat populations, RMR-2a, RMR-2b and RMR-2c of 178, 177 and 216 bp in length, respectively. Localization of an rDNA probe identified a single nucleolar organizing region on one autosome. Using a combination of labeled probes, we developed a preliminary karyotype for R. microplus. We present evidence that R. microplus has holocentric chromosomes and explore the implications of these findings for tick chromosome biology and genomic research.

Ticks and mites (subphylum Chelicerata; subclass Acari) include important pests of animals and plants worldwide. The Ixodes scapularis (black-legged tick) genome sequencing project marks the beginning of the genomics era for the field of acarology. This project is the first to sequence the genome of a blood-feeding tick vector of human disease and a member of the subphylum Chelicerata. Genome projects for other species of Acari are forthcoming and their genome sequences will likely feature significantly in the future of tick research. Parasitologists interested in advancing the field of tick genomics research will be faced with specific challenges. The development of genetic tools and resources, and the size and repetitive nature of tick genomes are important considerations. Innovative approaches may be required to sequence, assemble, annotate and analyse tick genomes. Overcoming these challenges will enable scientists to investigate the genes and genome organisation of this important group of arthropods and may ultimately lead to new solutions for control of ticks and tick-borne diseases.

The suborder Ixodida includes many tick species of medical and veterinary importance, but little is known about the genomic characteristics of these ticks. We report the first study to determine genome size in two species of Argasidae (soft ticks) and seven species of Ixodidae (hard ticks) using flow cytometry analysis of fluorescent stained nuclei. Our results indicate a large haploid genome size (1C>1000 Mbp) for all Ixodida with a mean of 1281 Mbp (1.31+/-0.07 pg) for the Argasidae and 2671 Mbp (2.73+/-0.04 pg) for the Ixodidae. The haploid genome size of Ixodes scapularis was determined to be 2262 Mbp. We observed inter- and intra-familial variation in genome size as well as variation between strains of the same species. We explore the implications of these results for tick genome evolution and tick genomics research.