Single-cell transcriptomics of the immune system in ME/CFS at baseline and following symptom provocation.

Myalgic encephalomyelitis/chronic fatigue syndrome (ME/CFS) is a serious and poorly understood disease. To understand immune dysregulation in ME/CFS, we use single-cell RNA sequencing (scRNA-seq) to examine immune cells in patient and control cohorts. Postexertional malaise (PEM), an exacerbation of symptoms following strenuous exercise, is a characteristic symptom of ME/CFS. To detect changes coincident with PEM, we applied scRNA-seq on the same cohorts following exercise. At baseline, ME/CFS patients display classical monocyte dysregulation suggestive of inappropriate differentiation and migration to tissue. We identify both diseased and more normal monocytes within patients, and the fraction of diseased cells correlates with disease severity. Comparing the transcriptome at baseline and postexercise challenge, we discover patterns indicative of improper platelet activation in patients, with minimal changes elsewhere in the immune system. Taken together, these data identify immunological defects present at baseline in patients and an additional layer of dysregulation in platelets.

Chemical, Molecular, and Single-nucleus Analysis Reveal Chondroitin Sulfate Proteoglycan Aberrancy in Fibrolamellar Carcinoma.

Fibrolamellar carcinoma (FLC) is an aggressive liver cancer with no effective therapeutic options. The extracellular environment of FLC tumors is poorly characterized and may contribute to cancer growth and/or metastasis. To bridge this knowledge gap, we assessed pathways relevant to proteoglycans, a major component of the extracellular matrix. We first analyzed gene expression data from FLC and nonmalignant liver tissue (n = 27) to identify changes in glycosaminoglycan (GAG) biosynthesis pathways and found that genes associated with production of chondroitin sulfate, but not other GAGs, are significantly increased by 8-fold. We then implemented a novel LC/MS-MS based method to quantify the abundance of different types of GAGs in patient tumors (n = 16) and found that chondroitin sulfate is significantly more abundant in FLC tumors by 6-fold. Upon further analysis of GAG-associated proteins, we found that versican (VCAN) expression is significantly upregulated at the mRNA and protein levels, the latter of which was validated by IHC. Finally, we performed single-cell assay for transposase-accessible chromatin sequencing on FLC tumors (n = 3), which revealed for the first time the different cell types in FLC tumors and also showed that VCAN is likely produced not only from FLC tumor epithelial cells but also activated stellate cells. Our results reveal a pathologic aberrancy in chondroitin (but not heparan) sulfate proteoglycans in FLC and highlight a potential role for activated stellate cells. Significance: This study leverages a multi-disciplinary approach, including state-of-the-art chemical analyses and cutting-edge single-cell genomic technologies, to identify for the first time a marked chondroitin sulfate aberrancy in FLC that could open novel therapeutic avenues in the future.

Multiple stages of evolutionary change in anthrax toxin receptor expression in humans.

The advent of animal husbandry and hunting increased human exposure to zoonotic pathogens. To understand how a zoonotic disease may have influenced human evolution, we study changes in human expression of anthrax toxin receptor 2 (ANTXR2), which encodes a cell surface protein necessary for Bacillus anthracis virulence toxins to cause anthrax disease. In immune cells, ANTXR2 is 8-fold down-regulated in all available human samples compared to non-human primates, indicating regulatory changes early in the evolution of modern humans. We also observe multiple genetic signatures consistent with recent positive selection driving a European-specific decrease in ANTXR2 expression in multiple tissues affected by anthrax toxins. Our observations fit a model in which humans adapted to anthrax disease following early ecological changes associated with hunting and scavenging, as well as a second period of adaptation after the rise of modern agriculture.

Parallel evolution of ancient, pleiotropic enhancers underlies butterfly wing pattern mimicry.

Color pattern mimicry in Heliconius butterflies is a classic case study of complex trait adaptation via selection on a few large effect genes. Association studies have linked color pattern variation to a handful of noncoding regions, yet the presumptive cis-regulatory elements (CREs) that control color patterning remain unknown. Here we combine chromatin assays, DNA sequence associations, and genome editing to functionally characterize 5 cis-regulatory elements of the color pattern gene optix We were surprised to find that the cis-regulatory architecture of optix is characterized by pleiotropy and regulatory fragility, where deletion of individual cis-regulatory elements has broad effects on both color pattern and wing vein development. Remarkably, we found orthologous cis-regulatory elements associate with wing pattern convergence of distantly related comimics, suggesting that parallel coevolution of ancestral elements facilitated pattern mimicry. Our results support a model of color pattern evolution in Heliconius where changes to ancient, multifunctional cis-regulatory elements underlie adaptive radiation.

Chromatin conformation remains stable upon extensive transcriptional changes driven by heat shock.

Heat shock (HS) initiates rapid, extensive, and evolutionarily conserved changes in transcription that are accompanied by chromatin decondensation and nucleosome loss at HS loci. Here we have employed in situ Hi-C to determine how heat stress affects long-range chromatin conformation in human and Drosophila cells. We found that compartments and topologically associating domains (TADs) remain unchanged by an acute HS. Knockdown of Heat Shock Factor 1 (HSF1), the master transcriptional regulator of the HS response, identified HSF1-dependent genes and revealed that up-regulation is often mediated by distal HSF1 bound enhancers. HSF1-dependent genes were usually found in the same TAD as the nearest HSF1 binding site. Although most interactions between HSF1 binding sites and target promoters were established in the nonheat shock (NHS) condition, a subset increased contact frequency following HS. Integrating information about HSF1 binding strength, RNA polymerase abundance at the HSF1 bound sites (putative enhancers), and contact frequency with a target promoter accurately predicted which up-regulated genes were direct targets of HSF1 during HS. Our results suggest that the chromatin conformation necessary for a robust HS response is preestablished in NHS cells of diverse metazoan species.