Stricturing Crohn's Disease Single-Cell RNA Sequencing Reveals Fibroblast Heterogeneity and Intercellular Interactions.

BACKGROUND & AIMS: Fibroblasts play a key role in stricture formation in Crohn's disease (CD) but understanding its pathogenesis requires a systems-level investigation to uncover new treatment targets. We studied full-thickness CD tissues to characterize fibroblast heterogeneity and function by generating the first single-cell RNA sequencing (scRNAseq) atlas of strictured bowel and providing proof of principle for therapeutic target validation. METHODS: We performed scRNAseq of 13 fresh full-thickness CD resections containing noninvolved, inflamed nonstrictured, and strictured segments as well as 7 normal non-CD bowel segments. Each segment was separated into mucosa/submucosa or muscularis propria and analyzed separately for a total of 99 tissue samples and 409,001 cells. We validated cadherin-11 (CDH11) as a potential therapeutic target by using whole tissues, isolated intestinal cells, NanoString nCounter, next-generation sequencing, proteomics, and animal models. RESULTS: Our integrated dataset revealed fibroblast heterogeneity in strictured CD with the majority of stricture-selective changes detected in the mucosa/submucosa, but not the muscle layer. Cell-cell interaction modeling revealed CXCL14+ as well as MMP/WNT5A+ fibroblasts displaying a central signaling role in CD strictures. CDH11, a fibroblast cell-cell adhesion molecule, was broadly expressed and up-regulated, and its profibrotic function was validated using NanoString nCounter, RNA sequencing, tissue target expression, in vitro gain- and loss-of-function experiments, proteomics, and knock-out and antibody-mediated CDH11 blockade in experimental colitis. CONCLUSIONS: A full-thickness bowel scRNAseq atlas revealed previously unrecognized fibroblast heterogeneity and interactions in CD strictures and CDH11 was validated as a potential therapeutic target. These results provide a new resource for a better understanding of CD stricture formation and open potential therapeutic developments. This work has been posted as a preprint on Biorxiv under doi: 10.1101/2023.04.03.534781.

Lactate dehydrogenase inhibition synergizes with IL-21 to promote CD8+ T cell stemness and antitumor immunity.

Interleukin (IL)-2 and IL-21 dichotomously shape CD8+ T cell differentiation. IL-2 drives terminal differentiation, generating cells that are poorly effective against tumors, whereas IL-21 promotes stem cell memory T cells (TSCM) and antitumor responses. Here we investigated the role of metabolic programming in the developmental differences induced by these cytokines. IL-2 promoted effector-like metabolism and aerobic glycolysis, robustly inducing lactate dehydrogenase (LDH) and lactate production, whereas IL-21 maintained a metabolically quiescent state dependent on oxidative phosphorylation. LDH inhibition rewired IL-2-induced effects, promoting pyruvate entry into the tricarboxylic acid cycle and inhibiting terminal effector and exhaustion programs, including mRNA expression of members of the NR4A family of nuclear receptors, as well as Prdm1 and Xbp1 While deletion of Ldha prevented development of cells with antitumor effector function, transient LDH inhibition enhanced the generation of memory cells capable of triggering robust antitumor responses after adoptive transfer. LDH inhibition did not significantly affect IL-21-induced metabolism but caused major transcriptomic changes, including the suppression of IL-21-induced exhaustion markers LAG3, PD1, 2B4, and TIM3. LDH inhibition combined with IL-21 increased the formation of TSCM cells, resulting in more profound antitumor responses and prolonged host survival. These findings indicate a pivotal role for LDH in modulating cytokine-mediated T cell differentiation and underscore the therapeutic potential of transiently inhibiting LDH during adoptive T cell-based immunotherapy, with an unanticipated cooperative antitumor effect of LDH inhibition and IL-21.

Immune Complex-Driven Generation of Human Macrophages with Anti-Inflammatory and Growth-Promoting Activity.

To maintain homeostasis, macrophages must be capable of assuming either an inflammatory or an anti-inflammatory phenotype. To better understand the latter, we stimulated human macrophages in vitro with TLR ligands in the presence of high-density immune complexes (IC). This combination of stimuli resulted in a broad suppression of inflammatory mediators and an upregulation of molecules involved in tissue remodeling and angiogenesis. Transcriptomic analysis of TLR stimulation in the presence of IC predicted the downstream activation of AKT and the inhibition of GSK3. Consequently, we pretreated LPS-stimulated human macrophages with small molecule inhibitors of GSK3 to partially phenocopy the regulatory effects of stimulation in the presence of IC. The upregulation of DC-STAMP and matrix metalloproteases was observed on these cells and may represent potential biomarkers for this regulatory activation state. To demonstrate the presence of these anti-inflammatory, growth-promoting macrophages in a human infectious disease, biopsies from patients with leprosy (Hanseniasis) were analyzed. The lepromatous form of this disease is characterized by hypergammaglobulinemia and defective cell-mediated immunity. Lesions in lepromatous leprosy contained macrophages with a regulatory phenotype expressing higher levels of DC-STAMP and lower levels of IL-12, relative to macrophages in tuberculoid leprosy lesions. Therefore, we propose that increased signaling by FcγR cross-linking on TLR-stimulated macrophages can paradoxically promote the resolution of inflammation and initiate processes critical to tissue growth and repair. It can also contribute to infectious disease progression.

Macrophages and the Recovery from Acute and Chronic Inflammation.

In recent years, researchers have devoted much attention to the diverse roles of macrophages and their contributions to tissue development, wound healing, and angiogenesis. What should not be lost in the discussions regarding the diverse biology of these cells is that when perturbed, macrophages are the primary contributors to potentially pathological inflammatory processes. Macrophages stand poised to rapidly produce large amounts of inflammatory cytokines in response to danger signals. The production of these cytokines can initiate a cascade of inflammatory mediator release that can lead to wholesale tissue destruction. The destructive inflammatory capability of macrophages is amplified by exposure to exogenous interferon-γ, which prolongs and heightens inflammatory responses. In simple terms, macrophages can thus be viewed as incendiary devices with hair triggers waiting to detonate. We have begun to ask questions about how these cells can be regulated to mitigate the collateral destruction associated with macrophage activation.

Exploring Choke Holds in Brazilian Jiujitsu Athletes: A Demographic Study.

Introduction Brazilian jiujitsu is a relatively new sport that has grown exponentially in popularity along with the growth of the Ultimate Fighting Championship (UFC). In jiujitsu, there are a variety of submissions with a choke hold being one of the most popular. There is a subset of athletes in jiujitsu who believes chokes are safe. However, there have been case reports of relatively young athletes suffering strokes secondary to internal carotid or vertebral artery dissections after being placed in choke holds. There have been manuscripts describing the injury profile in jiujitsu, but none mention stroke or dissections. This study evaluated how frequently chokes happen in jiujitsu and if athletes have ever experienced symptoms consistent with cervical artery dissection (CAD). Additionally, this study aimed to describe the training frequency and baseline demographics of jiujitsu athletes. Methods A survey was distributed throughout social media platforms which asked both quantitative and qualitative questions regarding athlete training. The survey consisted of 28 questions which collected largely baseline grappling information about the participants such as how long they trained, how often they spar, favorite submission, how frequently they are choked, etc. This data was then analyzed using odds ratio and one sample t-test to evaluate for statistical differences. Results A total of 521 participants were included in the analysis. The participants were mostly male (84.7%), trained for four years, four times per week; 99.8% (520) participated in sparring, with an average age of 37; and 55.7% (290) have experienced symptoms consistent with CAD. Descriptive statistics revealed that individuals who were 37 years of age or younger were more likely to experience symptoms consistent with CAD (odds ratio: 1.5337 (95% confidence interval (CI): 1.0827-2.1727). Athletes that were 37 years of age or younger have been training for fewer years (4.7 years vs 8.8 years) but train more days per week (4.03 times per week vs 3.76 time per week), drill for a longer amount of time (46.8 minutes per class vs 38.3 minutes per class), attend longer classes (81.12 minutes vs 72.3 minutes), and train for a longer period of time per week (338.5 minutes vs 274.6 minutes) than athletes over 37 years. All previously mentioned variables were analyzed using a one sample t-test and were significant at the α = 0.05 level. The lone qualitative question regarding the term "train brain" revealed that of those who experienced it, 84.1% (58) described it as a cognitive/physical impairing event. Conclusion Jiujitsu athletes train multiple times per week and are frequently exposed to choke holds. There is no literature to examine the long-term effects of these chokes on the athlete's cervical vasculature. Additional studies should be conducted to evaluate the effects of the repetitive stress placed on these vessels.

Peptide-directed preparation and X-ray structural study of Au nanoparticles on titanium surfaces.

We report the peptide-directed preparation and X-ray structural study of biofunctionalized Au nanoparticles (NPs) deposited on Ti surfaces. Au NPs were prepared by reduction of Au(3+) compound onto HCl-refreshed Ti in the presence of thiol-functionalized small peptides. A modified extended X-ray absorption fine-structure (EXAFS) technique, equipped with a rotating-stage and glancing-angle setup, was able to more sensitively detect the structure and bonding of Au NPs on Ti with low surface coverage. It was found that the use of the tripeptide glutathione (GSH) results in smaller NP size when compared to N-(2-mercapto-propionyl) glycine (MPG), a pseudodipeptide, over a wide range of Au/peptide molar ratios (20:1, 10:1, 5:1, and 2:1). By varying the ligand concentration, the Au NP structure in both systems can be controlled, generating nanocrystals, nanoclusters, and Au-thiolate polymer, which is unique for substrate-supported NP synthesis. This work presents a facile preparation of Au-peptide nanoparticles on biocompatible surfaces, and illustrates the high sensitivity of this modified EXAFS technique for structural studies of substrate-supported nanoparticles with low coverage.

Single-cell isolation from full-thickness human intestinal tissue resections for single-cell RNA sequencing.

Single-cell isolation techniques allow the investigation of physical and functional relationships between individual cells within a complex cell population. Here, we present a protocol for single-cell isolation from full-thickness intestinal tissue resections. We describe steps for pre-processing specimens, isolation of lamina propria and muscular layers, and red blood cell lysis. We then detail fixation of isolated cells and assessment of cell quality. The resulting cell suspension can be subjected to RNA sequencing on the 10× Chromium platform. For complete details on the use and execution of this protocol, please refer to Mukherjee et al.1.

Identification of a broadly fibrogenic macrophage subset induced by type 3 inflammation.

Macrophages are central orchestrators of the tissue response to injury, with distinct macrophage activation states playing key roles in fibrosis progression and resolution. Identifying key macrophage populations found in human fibrotic tissues could lead to new treatments for fibrosis. Here, we used human liver and lung single-cell RNA sequencing datasets to identify a subset of CD9+TREM2+ macrophages that express SPP1, GPNMB, FABP5, and CD63. In both human and murine hepatic and pulmonary fibrosis, these macrophages were enriched at the outside edges of scarring and adjacent to activated mesenchymal cells. Neutrophils expressing MMP9, which participates in the activation of TGF-ß1, and the type 3 cytokines GM-CSF and IL-17A coclustered with these macrophages. In vitro, GM-CSF, IL-17A, and TGF-ß1 drive the differentiation of human monocytes into macrophages expressing scar-associated markers. Such differentiated cells could degrade collagen IV but not collagen I and promote TGF-ß1-induced collagen I deposition by activated mesenchymal cells. In murine models blocking GM-CSF, IL-17A or TGF-ß1 reduced scar-associated macrophage expansion and hepatic or pulmonary fibrosis. Our work identifies a highly specific macrophage population to which we assign a profibrotic role across species and tissues. It further provides a strategy for unbiased discovery, triage, and preclinical validation of therapeutic targets based on this fibrogenic macrophage population.

Distinct use of super-enhancer elements controls cell type-specific CD25 transcription and function.

The IL-2 receptor α chain (IL-2Rα/CD25) is constitutively expressed on double-negative (DN2/DN3 thymocytes and regulatory T cells (Tregs) but induced by IL-2 on T and natural killer (NK) cells, with Il2ra expression regulated by a STAT5-dependent super-enhancer. We investigated CD25 regulation and function using a series of mice with deletions spanning STAT5-binding elements. Deleting the upstream super-enhancer region mainly affected constitutive CD25 expression on DN2/DN3 thymocytes and Tregs, with these mice developing autoimmune alopecia, whereas deleting an intronic region decreased IL-2-induced CD25 on peripheral T and NK cells. Thus, distinct super-enhancer elements preferentially control constitutive versus inducible expression in a cell type-specific manner. The mediator-1 coactivator colocalized with specific STAT5-binding sites. Moreover, both upstream and intronic regions had extensive chromatin interactions, and deletion of either region altered the super-enhancer structure in mature T cells. These results demonstrate differential functions for distinct super-enhancer elements, thereby indicating previously unknown ways to manipulate CD25 expression in a cell type-specific fashion.

Stricturing Crohn's disease single-cell RNA sequencing reveals fibroblast heterogeneity and intercellular interactions.

Background: Fibroblasts play a key role in stricture formation in Crohn's disease (CD) but understanding it's pathogenesis requires a systems-level investigation to uncover new treatment targets. We studied full thickness CD tissues to characterize fibroblast heterogeneity and function by generating the first single cell RNA sequencing (scRNAseq) atlas of strictured bowel and providing proof of principle for therapeutic target validation. Methods: We performed scRNAseq of 13 fresh full thickness CD resections containing non-involved, inflamed non-strictured, and strictured segments as well as 7 normal non-CD bowel segments. Each segment was separated into mucosa/submucosa or muscularis propria and analyzed separately for a total of 99 tissue samples and 409,001 cells. We validated cadherin-11 (CDH11) as a potential therapeutic target by using whole tissues, isolated intestinal cells, NanoString nCounter, next generation sequencing, proteomics and animal models. Results: Our integrated dataset revealed fibroblast heterogeneity in strictured CD with the majority of stricture-selective changes detected in the mucosa/submucosa, but not the muscle layer. Cell-cell interaction modeling revealed CXCL14+ as well as MMP/WNT5A+ fibroblasts displaying a central signaling role in CD strictures. CDH11, a fibroblast cell-cell adhesion molecule, was broadly expressed and upregulated, and its pro-fibrotic function was validated by NanoString nCounter, RNA sequencing, tissue target expression, in vitro gain- and loss-of-function experiments, proteomics, and two animal models of experimental colitis. Conclusion: A full-thickness bowel scRNAseq atlas revealed previously unrecognized fibroblast heterogeneity and interactions in CD strictures and CDH11 was validated as a potential therapeutic target. These results provide a new resource for a better understanding of CD stricture formation and opens potential therapeutic developments.

Standardization of epidermal growth factor receptor (EGFR) measurement by quantitative immunofluorescence and impact on antibody-based mutation detection in non-small cell lung cancer.

Challenges in measurement of epidermal growth factor receptor (EGFR) protein expression have led to conflicting data on its prognostic value and discontinuation of its use for prediction of response. Herein is described a quantitative standardized assay for EGFR and its use in a series of retrospective cohorts of patients with non-small cell lung cancer (NSCLC). The AQUA technology of quantitative immunofluorescence was used in conjunction with Western blot analysis to calculate the absolute concentration of EGFR in two independent NSCLC cohorts (170 from Yale New Haven Hospital and 335 from Sotiria and Patras Hospitals in Greece). EGFR and mutated EGFR were measured using D38B1 antibody and two mutation-specific antibodies. All patients positive or borderline for mutation-specific antibody were genotyped. A threshold for reproducible detection of EGFR was defined as 0.85 ng/µg total protein. EGFR expression demonstrated no prognostic value in either cohort. The mutation rate was 1.79% in the Yale cohort, and 1.52% in the Sotiria/Patras cohort, with no antibody detection-based false-positive cases. No mutations were detected for EGFR concentrations <1.46 ng/µg total protein. In summary, accurate measurement of EGFR still shows no prognostic value in NSCLC. In these two population-based cohorts, the antibody-based EGFR mutation rate was lower than has been frequently reported.

Biomolecule-coated metal nanoparticles on titanium.

Immobilizations of nanoparticles and biomolecules on biocompatible substrates such as titanium are two promising approaches to bringing new functionalities to Ti-based biomaterials. Herein, we used a variety of X-ray spectroscopic techniques to study and better understand metal-thiolate interactions in biofunctionalized metal nanoparticle systems supported on Ti substrates. Using a facile one-step procedure, a series of Au nanoparticle samples with varied biomolecule coatings ((2-mercatopropionyl)glycine (MPG) and bovine serum albumin (BSA)) and biomolecule concentrations are prepared. Ag and Pd systems are also studied to observe change with varying metal composition. The structure and properties of these biomolecule-coated nanoparticles are investigated with scanning electron microscopy (SEM) and element-specific X-ray techniques, including extended X-ray absorption fine structure (Au L(3)-edge), X-ray absorption near-edge structure (Au L(3), Ag L(3), Pd L(3), and S K-edge), and X-ray photoelectron spectroscopy (Au 4f, Ag 3d, Pd 3d, and S 2p core level). It was found that, by comparison of SEM and X-ray spectroscopy results, the coating of metal nanoparticles with varying model biomolecule systems can have a significant effect on both surface coverage and organization. This work offers a facile chemical method for bio- and nanofunctionalization of Ti substrates as well as provides a physical picture of the structure and bonding of biocoated metal nanoparticles, which may lead to useful applications in orthopedics and biomedicine.

IFNγ regulates NAD+ metabolism to promote the respiratory burst in human monocytes.

Interferon γ (IFNγ) is an essential and pleiotropic activator of human monocytes, but little is known about the changes in cellular metabolism required for IFNγ-induced activation. We sought to elucidate the mechanisms by which IFNγ reprograms monocyte metabolism to support its immunologic activities. We found that IFNγ increased oxygen consumption rates (OCR) in monocytes, indicative of reactive oxygen species generation by both mitochondria and nicotinamide adenine dinucleotide phosphate (NADPH) oxidase. Transcriptional profiling revealed that this oxidative phenotype was driven by IFNγ-induced reprogramming of NAD+ metabolism, which is dependent on nicotinamide phosphoribosyltransferase (NAMPT)-mediated NAD+ salvage to generate NADH and NADPH for oxidation by mitochondrial complex I and NADPH oxidase, respectively. Consistent with this pathway, monocytes from patients with gain-of-function mutations in STAT1 demonstrated higher-than-normal OCR, whereas chemical or genetic disruption of mitochondrial complex I (rotenone treatment or Leigh syndrome patient monocytes) or NADPH oxidase (diphenyleneiodonium treatment or chronic granulomatous disease [CGD] patient monocytes) reduced OCR. Interestingly, inhibition of NAMPT in healthy monocytes completely abrogated the IFNγ-induced oxygen consumption, comparable to levels observed in CGD monocytes. These data identify an IFNγ-induced, NAMPT-dependent, NAD+ salvage pathway that is critical for IFNγ activation of human monocytes.

Single cell sequencing identifies clonally expanded synovial CD4+ TPH cells expressing GPR56 in rheumatoid arthritis.

Rheumatoid arthritis (RA) is an autoimmune disease affecting synovial joints where different CD4+ T cell subsets may contribute to pathology. Here, we perform single cell sequencing on synovial CD4+ T cells from anti-citrullinated protein antibodies (ACPA)+ and ACPA- RA patients and identify two peripheral helper T cell (TPH) states and a cytotoxic CD4+ T cell subset. We show that the adhesion G-protein coupled receptor 56 (GPR56) delineates synovial CXCL13high TPH CD4+ T cells expressing LAG-3 and the tissue-resident memory receptors CXCR6 and CD69. In ACPA- SF, TPH cells display lower levels of GPR56 and LAG-3. Further, most expanded T cell clones in the joint are within CXCL13high TPH CD4+ T cells. Finally, RNA-velocity analyses suggest a common differentiation pathway between the two TPH clusters and effector CD4+ T cells. Our study provides comprehensive immunoprofiling of the synovial CD4+ T cell subsets in ACPA+ and ACPA- RA.

Humoral immunity in leishmaniasis - Prevention or promotion of parasite growth?

Leishmaniasis can present as a "spectrum" of clinical outcomes. There is evidence that these divergent clinical outcomes are attributable to genetic differences in the human host [1] as well the species of infecting parasite [2]. The spectrum of disease has largely been described by defining the polar opposites of T cell immune responses. In the mouse model, a TH1 immune response is associated with low numbers of Leishmania parasites in lesions, whereas a TH2 immune response has been associated with unrestricted parasite growth. In the present work, we revisit leishmaniasis and seek to better define the clinical spectrum as a function of divergent humoral immune responses. We describe examples in human, canine, and even some murine models of leishmaniasis that reveal a direct correlation between high anti-parasite antibody responses and unrestricted parasite growth. Therefore, we propose that the spectral nature of this disease may be due to quantitative and qualitative differences in the antibodies that are produced during disease. In human visceral leishmaniasis, a decrease in anti-parasite antibody levels may actually predict disease resolution. Thus, rather than defining this disease as a simple TH1/TH2 dichotomy, we propose that clinical leishmaniasis depends on the degree of humoral immunity, with high IgG predicting parasite persistence. These observations have obvious implications for vaccine development in leishmaniasis, and they may extend to other diseases caused by intracellular pathogens.

Host and parasite responses in human diffuse cutaneous leishmaniasis caused by L. amazonensis.

Diffuse cutaneous leishmaniasis (DCL) is a rare form of leishmaniasis where parasites grow uncontrolled in diffuse lesions across the skin. Meta-transcriptomic analysis of biopsies from DCL patients infected with Leishmania amazonensis demonstrated an infiltration of atypical B cells producing a surprising preponderance of the IgG4 isotype. DCL lesions contained minimal CD8+ T cell transcripts and no evidence of persistent TH2 responses. Whereas localized disease exhibited activated (so-called M1) macrophage presence, transcripts in DCL suggested a regulatory macrophage (R-Mϕ) phenotype with higher levels of ABCB5, DCSTAMP, SPP1, SLAMF9, PPARG, MMPs, and TM4SF19. The high levels of parasite transcripts in DCL and the remarkable uniformity among patients afforded a unique opportunity to study parasite gene expression in this disease. Patterns of parasite gene expression in DCL more closely resembled in vitro parasite growth in resting macrophages, in the absence of T cells. In contrast, parasite gene expression in LCL revealed 336 parasite genes that were differently upregulated, relative to DCL and in vitro macrophage growth, and these transcripts may represent transcripts that are produced by the parasite in response to host immune pressure.

Correction: Meta-transcriptome Profiling of the Human-Leishmania braziliensis Cutaneous Lesion.

[This corrects the article DOI: 10.1371/journal.pntd.0004992.].

Image-guided Coring for Large-scale Studies in Molecular Pathology.

Sampling of formalin-fixed paraffin-embedded (FFPE) tissue blocks is a critical initial step in molecular pathology. Image-guided coring (IGC) is a new method for using digital pathology images to guide tissue block coring for molecular analyses. The goal of our study is to evaluate the use of IGC for both tissue-based and nucleic acid-based projects in molecular pathology. First, we used IGC to construct a tissue microarray (TMA); second, we used IGC for FFPE block sampling followed by RNA extraction; and third, we assessed the correlation between nuclear counts quantitated from the IGC images and RNA yields. We used IGC to construct a TMA containing 198 normal and breast cancer cores. Histopathologic analysis showed high accuracy for obtaining tumor and normal breast tissue. Next, we used IGC to obtain normal and tumor breast samples before RNA extraction. We selected a random subset of tumor and normal samples to perform computational image analysis to quantify nuclear density, and we built regression models to estimate RNA yields from nuclear count, age of the block, and core diameter. Number of nuclei and core diameter were the strongest predictors of RNA yields in both normal and tumor tissue. IGC is an effective method for sampling FFPE tissue blocks for TMA construction and nucleic acid extraction. We identify significant associations between quantitative nuclear counts obtained from IGC images and RNA yields, suggesting that the integration of computational image analysis with IGC may be an effective approach for tumor sampling in large-scale molecular studies.

Meta-transcriptome Profiling of the Human-Leishmania braziliensis Cutaneous Lesion.

Host and parasite gene expression in skin biopsies from Leishmania braziliensis-infected patients were simultaneously analyzed using high throughput RNA-sequencing. Biopsies were taken from 8 patients with early cutaneous leishmaniasis and 17 patients with late cutaneous leishmaniasis. Although parasite DNA was found in all patient lesions at the time of biopsy, the patients could be stratified into two groups: one lacking detectable parasite transcripts (PTNeg) in lesions, and another in which parasite transcripts were readily detected (PTPos). These groups exhibited substantial differences in host responses to infection. PTPos biopsies contained an unexpected increase in B lymphocyte-specific and immunoglobulin transcripts in the lesions, and an upregulation of immune inhibitory molecules. Biopsies without detectable parasite transcripts showed decreased evidence for B cell activation, but increased expression of antimicrobial genes and genes encoding skin barrier functions. The composition and abundance of L. braziliensis transcripts in PTPos lesions were surprisingly conserved among all six patients, with minimal meaningful differences between lesions from patients with early and late cutaneous leishmaniasis. The most abundant parasite transcripts expressed in lesions were distinct from transcripts expressed in vitro in human macrophage cultures infected with L. amazonensis or L. major. Therefore in vitro gene expression in macrophage monolayers may not be a strong predictor of gene expression in lesions. Some of the most highly expressed in vivo transcripts encoded amastin-like proteins, hypothetical genes, putative parasite virulence factors, as well as histones and tubulin. In summary, RNA sequencing allowed us to simultaneously analyze human and L. braziliensis transcriptomes in lesions of infected patients, and identify unexpected differences in host immune responses which correlated with active transcription of parasite genes.

Extensive rewiring of epithelial-stromal co-expression networks in breast cancer.

BACKGROUND: Epithelial-stromal crosstalk plays a critical role in invasive breast cancer pathogenesis; however, little is known on a systems level about how epithelial-stromal interactions evolve during carcinogenesis. RESULTS: We develop a framework for building genome-wide epithelial-stromal co-expression networks composed of pairwise co-expression relationships between mRNA levels of genes expressed in the epithelium and stroma across a population of patients. We apply this method to laser capture micro-dissection expression profiling datasets in the setting of breast carcinogenesis. Our analysis shows that epithelial-stromal co-expression networks undergo extensive rewiring during carcinogenesis, with the emergence of distinct network hubs in normal breast, and estrogen receptor-positive and estrogen receptor-negative invasive breast cancer, and the emergence of distinct patterns of functional network enrichment. In contrast to normal breast, the strongest epithelial-stromal co-expression relationships in invasive breast cancer mostly represent self-loops, in which the same gene is co-expressed in epithelial and stromal regions. We validate this observation using an independent laser capture micro-dissection dataset and confirm that self-loop interactions are significantly increased in cancer by performing computational image analysis of epithelial and stromal protein expression using images from the Human Protein Atlas. CONCLUSIONS: Epithelial-stromal co-expression network analysis represents a new approach for systems-level analyses of spatially localized transcriptomic data. The analysis provides new biological insights into the rewiring of epithelial-stromal co-expression networks and the emergence of epithelial-stromal co-expression self-loops in breast cancer. The approach may facilitate the development of new diagnostics and therapeutics targeting epithelial-stromal interactions in cancer.