IL-17 signaling in primary sclerosing cholangitis patient-derived organoids.

BACKGROUND: The pathogenesis of primary sclerosing cholangitis (PSC) is unclear, although studies implicate IL-17A as an inflammatory mediator in this disease. However, a direct assessment of IL-17 signaling in PSC cholangiocytes is lacking. In this study, we aimed to investigate and characterize the response of PSC extrahepatic cholangiocyte organoids (ECO) to IL-17A stimulation. METHODS: Cholangiocytes obtained from patients with PSC and without PSC by endoscopic retrograde cholangiography were cultured as ECO. The ECO were treated with vehicle or IL-17A and assessed by transcriptomics, secretome analysis, and genome sequencing. RESULTS: Unsupervised clustering of all integrated single-cell RNA sequencing data identified 8 cholangiocyte clusters that did not differ between PSC and non-PSC ECO. However, PSC ECO cells demonstrated a robust response to IL-17 treatment, as noted by an increased number of differentially expressed genes by transcriptomics and more abundant chemokine and cytokine expression and secretion. After rigorous filtering, genome sequencing identified candidate somatic variants shared among PSC ECO from unrelated individuals. However, no candidate rare variants in genes regulating the IL-17 pathway were identified, but rare variants regulating the MAPK signaling pathway were present in all PSC ECO. CONCLUSIONS: PSC and non-PSC patient-derived ECO respond differently to IL-17 stimulation, implicating this pathway in the pathogenesis of PSC.

Tumor necrosis factor-related apoptosis-inducing ligand (TRAIL) deletion in myeloid cells augments cholestatic liver injury.

Ductular reactive (DR) cells exacerbate cholestatic liver injury and fibrosis. Herein, we posit that tumor necrosis factor-related apoptosis-inducing ligand (TRAIL) emanates from recruited macrophages and restrains DR cell expansion, thereby limiting cholestatic liver injury. Wild type (WT), Trailfl/fl and myeloid-specific Trail deleted (TrailΔmye) C57BL/6 mice were exposed to DDC diet-induced cholestatic liver injury, which induced hepatomegaly and liver injury as compared to control diet-fed mice. However, parameters of liver injury, fibrosis, and inflammation were all increased in the TrailΔmye mice as compared to the WT and Trailfl/fl mice. High dimensional mass cytometry indicated that cholestasis resulted in increased hepatic recruitment of subsets of macrophages and neutrophils in the TrailΔmye mice. Spatial transcriptomics analysis revealed that the PanCK+ cholangiocytes from TrailΔmye mice had increased expression of the known myeloid attractants S100a8, Cxcl5, Cx3cl1, and Cxcl1. Additionally, in situ hybridization of Cxcl1, a potent neutrophil chemoattractant, demonstrated an increased expression in CK19+ cholangiocytes of TrailΔmye mice. Collectively, these data suggest that TRAIL from myeloid cells, particularly macrophages, restrains a subset of DR cells (i.e., Cxcl1 positive cells), limiting liver inflammation and fibrosis. Reprogramming macrophages to express TRAIL may be salutary in cholestasis.

Interspecies regulatory landscapes and elements revealed by novel joint systematic integration of human and mouse blood cell epigenomes.

Knowledge of locations and activities of cis -regulatory elements (CREs) is needed to decipher basic mechanisms of gene regulation and to understand the impact of genetic variants on complex traits. Previous studies identified candidate CREs (cCREs) using epigenetic features in one species, making comparisons difficult between species. In contrast, we conducted an interspecies study defining epigenetic states and identifying cCREs in blood cell types to generate regulatory maps that are comparable between species, using integrative modeling of eight epigenetic features jointly in human and mouse in our V al i dated S ystematic I ntegrati on (VISION) Project. The resulting catalogs of cCREs are useful resources for further studies of gene regulation in blood cells, indicated by high overlap with known functional elements and strong enrichment for human genetic variants associated with blood cell phenotypes. The contribution of each epigenetic state in cCREs to gene regulation, inferred from a multivariate regression, was used to estimate epigenetic state Regulatory Potential (esRP) scores for each cCRE in each cell type, which were used to categorize dynamic changes in cCREs. Groups of cCREs displaying similar patterns of regulatory activity in human and mouse cell types, obtained by joint clustering on esRP scores, harbored distinctive transcription factor binding motifs that were similar between species. An interspecies comparison of cCREs revealed both conserved and species-specific patterns of epigenetic evolution. Finally, we showed that comparisons of the epigenetic landscape between species can reveal elements with similar roles in regulation, even in the absence of genomic sequence alignment.

IL-17 Signaling in Primary Sclerosing Cholangitis Patient-Derived Organoids.

The pathogenesis of primary sclerosing cholangitis (PSC) is unclear, although studies implicate IL-17A as an inflammatory mediator in this disease. However, a direct assessment of IL-17 signaling in PSC cholangiocytes is lacking. In this study we aimed to investigate the response of PSC extrahepatic cholangiocyte organoids (ECO) to IL-17A stimulation. Cholangiocytes obtained from PSC and non-PSC patients by endoscopic retrograde cholangiography (ERC) were cultured as ECO. The ECO were treated with vehicle or IL-17A and assessed by transcriptomics, secretome analysis, and genome sequencing (GS). Unsupervised clustering of all integrated scRNA-seq data identified 8 cholangiocyte clusters which did not differ between PSC and non-PSC ECO. However, PSC ECO cells demonstrated a robust response to IL-17 treatment, noted by an increased number of differentially expressed genes (DEG) by transcriptomics, and more abundant chemokine and cytokine expression and secretion. After rigorous filtering, GS identified candidate somatic variants shared among PSC ECO from unrelated individuals. However, no candidate rare variants in genes regulating the IL-17 pathway were identified, but rare variants regulating the MAPK signaling pathway were present in all PSC ECO. In conclusion, PSC and non-PSC patient derived ECO respond differently to IL-17 stimulation implicating this pathway in the pathogenesis of PSC.

Telomere dysfunction promotes cholangiocyte senescence and biliary fibrosis in primary sclerosing cholangitis.

Cellular senescence and biliary fibrosis are prototypical features of obliterative cholangiopathies, such as primary sclerosing cholangitis (PSC). Telomere dysfunction can lead to senescence either through telomere erosion or damaged telomeres. Our goal was to investigate a mechanistic relationship between telomere damage and biliary fibrosis in PSC. Telomere attrition was observed in the bile ducts of patients with PSC along with a reduction in telomerase reverse transcriptase (TERT) expression, compared with that in normal livers. Similarly, liver tissue from mouse models of biliary fibrosis showed telomere attrition with increased damage at telomeres measured as telomere-associated foci (TAF). Cellular models of senescence induction increased the TAF in cholangiocytes. This coincided with decreased TERT expression and increased senescence, which was rescued by modulating TERT levels. Epigenetic analysis revealed increased acquisition of repressive histone methylation at the TERT promoter, which correlated with decreased TERT transcription. Cholangiocyte-selective deletion of TERT in mice exacerbated fibrosis, whereas androgen therapy toward telomerase rescued liver fibrosis and liver function in a genetic mouse model of PSC. Our results demonstrate a mechanistic role for telomere dysfunction in cellular senescence and fibrosis that characterize PSC. This suggests that PSC may be, in part, a telomere biology disorder, and identifies TERT as a potential therapeutic target.

Identification of a ΔNp63-Dependent Basal-Like A Subtype-Specific Transcribed Enhancer Program (B-STEP) in Aggressive Pancreatic Ductal Adenocarcinoma.

A major hurdle to the application of precision oncology in pancreatic cancer is the lack of molecular stratification approaches and targeted therapy for defined molecular subtypes. In this work, we sought to gain further insight and identify molecular and epigenetic signatures of the Basal-like A pancreatic ductal adenocarcinoma (PDAC) subgroup that can be applied to clinical samples for patient stratification and/or therapy monitoring. We generated and integrated global gene expression and epigenome mapping data from patient-derived xenograft models to identify subtype-specific enhancer regions that were validated in patient-derived samples. In addition, complementary nascent transcription and chromatin topology (HiChIP) analyses revealed a Basal-like A subtype-specific transcribed enhancer program in PDAC characterized by enhancer RNA (eRNA) production that is associated with more frequent chromatin interactions and subtype-specific gene activation. Importantly, we successfully confirmed the validity of eRNA detection as a possible histologic approach for PDAC patient stratification by performing RNA-ISH analyses for subtype-specific eRNAs on pathologic tissue samples. Thus, this study provides proof-of-concept that subtype-specific epigenetic changes relevant for PDAC progression can be detected at a single-cell level in complex, heterogeneous, primary tumor material. IMPLICATIONS: Subtype-specific enhancer activity analysis via detection of eRNAs on a single-cell level in patient material can be used as a potential tool for treatment stratification.

Cell Type-Specific Induction of Inflammation-Associated Genes in Crohn's Disease and Colorectal Cancer.

Based on the rapid increase in incidence of inflammatory bowel disease (IBD), the identification of susceptibility genes and cell populations contributing to this condition is essential. Previous studies suggested multiple genes associated with the susceptibility of IBD; however, due to the analysis of whole-tissue samples, the contribution of individual cell populations remains widely unresolved. Single-cell RNA sequencing (scRNA-seq) provides the opportunity to identify underlying cellular populations. We determined the enrichment of Crohn's disease (CD)-induced genes in a publicly available Crohn's disease scRNA-seq dataset and detected the strongest induction of these genes in innate lymphoid cells (ILC1), highly activated T cells and dendritic cells, pericytes and activated fibroblasts, as well as epithelial cells. Notably, these genes were highly enriched in IBD-associated neoplasia, as well as sporadic colorectal cancer (CRC). Indeed, the same six cell populations displayed an upregulation of CD-induced genes in a CRC scRNA-seq dataset. Finally, after integrating and harmonizing the CD and CRC scRNA-seq data, we demonstrated that these six cell types display a gradual increase in gene expression levels from a healthy state to an inflammatory and tumorous state. Together, we identified cell populations that specifically upregulate CD-induced genes in CD and CRC patients and could, therefore, contribute to inflammation-associated tumor development.

BMI1 maintains the Treg epigenomic landscape to prevent inflammatory bowel disease.

FOXP3+ Tregs are expanded within the inflamed intestine of human Crohn's disease, yet FOXP3-mediated gene repression within these cells is lost. The polycomb repressive complexes play a role in FOXP3 target gene regulation, but deeper mechanistic insight is incomplete. We have now specifically identified the polycomb-repressive complex 1 (PRC1) family member, BMI1 in the regulation of a proinflammatory enhancer network in both human and murine Tregs. Using human Tregs and lamina propria T cells, we inferred PRC1 to regulate Crohn's associated gene networks through assays of chromatin accessibility. Conditional deletion of BMI1 in murine FOXP3+ cells led to systemic inflammation. BMI1-deficient Tregs beared a TH1/TH17-like phenotype as assessed by assays of genome wide transcription, chromatin accessibility and proteomic techniques. Finally, BMI1 mutant FOXP3+ cells did not suppress colitis in the adoptive transfer model of human inflammatory bowel disease. We propose that BMI1 plays an important role in enforcing Treg identity in vitro and in vivo. Loss of Treg identity via genetic or transient BMI1 depletion perturbs the epigenome and converts Tregs into Th1/Th17-like proinflammatory cells, a transition relevant to human Crohn's disease associated CD4+ T cells.

Rhinovirus Reduces the Severity of Subsequent Respiratory Viral Infections by Interferon-Dependent and -Independent Mechanisms.

Coinfection by heterologous viruses in the respiratory tract is common and can alter disease severity compared to infection by individual virus strains. We previously found that inoculation of mice with rhinovirus (RV) 2 days before inoculation with a lethal dose of influenza A virus [A/Puerto Rico/8/34 (H1N1) (PR8)] provides complete protection against mortality. Here, we extended that finding to a second lethal respiratory virus, pneumonia virus of mice (PVM), and analyzed potential mechanisms of RV-induced protection. RV completely prevented mortality and weight loss associated with PVM infection. Major changes in host gene expression upon PVM infection were delayed compared to PR8. RV induced earlier recruitment of inflammatory cells, which were reduced at later times in RV-inoculated mice. Findings common to both virus pairs included the upregulated expression of mucin-associated genes and dampening of inflammation-related genes in mice that were inoculated with RV before lethal virus infection. However, type I interferon (IFN) signaling was required for RV-mediated protection against PR8 but not PVM. IFN signaling had minor effects on PR8 replication and contributed to controlling neutrophilic inflammation and hemorrhagic lung pathology in RV/PR8-infected mice. These findings, combined with differences in virus replication levels and disease severity, suggest that the suppression of inflammation in RV/PVM-infected mice may be due to early, IFN-independent suppression of viral replication, while that in RV/PR8-infected mice may be due to IFN-dependent modulation of immune responses. Thus, a mild upper respiratory viral infection can reduce the severity of a subsequent severe viral infection in the lungs through virus-dependent mechanisms. IMPORTANCE Respiratory viruses from diverse families cocirculate in human populations and are frequently detected within the same host. Although clinical studies suggest that infection by multiple different respiratory viruses may alter disease severity, animal models in which we can control the doses, timing, and strains of coinfecting viruses are critical to understanding how coinfection affects disease severity. Here, we compared gene expression and immune cell recruitment between two pairs of viruses (RV/PR8 and RV/PVM) inoculated sequentially in mice, both of which result in reduced severity compared to lethal infection by PR8 or PVM alone. Reduced disease severity was associated with suppression of inflammatory responses in the lungs. However, differences in disease kinetics and host and viral gene expression suggest that protection by coinfection with RV may be due to distinct molecular mechanisms. Indeed, we found that antiviral cytokine signaling was required for RV-mediated protection against lethal infection by PR8 but not PVM.

Effects of sheared chromatin length on ChIP-seq quality and sensitivity.

Chromatin immunoprecipitation followed by massively parallel, high throughput sequencing (ChIP-seq) is the method of choice for genome-wide identification of DNA segments bound by specific transcription factors or in chromatin with particular histone modifications. However, the quality of ChIP-seq datasets varies widely, with a substantial fraction being of intermediate to poor quality. Thus, it is important to discern and control the factors that contribute to variation in ChIP-seq. In this study, we focused on sonication, a user-controlled variable, to produce sheared chromatin. We systematically varied the amount of shearing of fixed chromatin from a mouse erythroid cell line, carefully measuring the distribution of resultant fragment lengths prior to ChIP-seq. This systematic study was complemented with a retrospective analysis of additional experiments. We found that the level of sonication had a pronounced impact on the quality of ChIP-seq signals. Over-sonication consistently reduced quality, while the impact of under-sonication differed among transcription factors, with no impact on sites bound by CTCF but frequently leading to the loss of sites occupied by TAL1 or bound by POL2. The bound sites not observed in low-quality datasets were inferred to be a mix of both direct and indirect binding. We leveraged these findings to produce a set of CTCF ChIP-seq datasets in rare, primary hematopoietic progenitor cells. Our observation that the amount of chromatin sonication is a key variable in success of ChIP-seq experiments indicates that monitoring the level of sonication can improve ChIP-seq quality and reproducibility and facilitate ChIP-seq in rare cell types.

STIM1 Mediates Calcium-Dependent Epigenetic Reprogramming in Pancreatic Cancer.

Pancreatic ductal adenocarcinoma (PDAC) displays a dismal prognosis due to late diagnosis and high chemoresistance incidence. For advanced disease stages or patients with comorbidities, treatment options are limited to gemcitabine alone or in combination with other drugs. While gemcitabine resistance has been widely attributed to the levels of one of its targets, RRM1, the molecular consequences of gemcitabine resistance in PDAC remain largely elusive. Here we sought to identify genomic, epigenomic, and transcriptomic events associated with gemcitabine resistance in PDAC and their potential clinical relevance. We found that gemcitabine-resistant cells displayed a coamplification of the adjacent RRM1 and STIM1 genes. Interestingly, RRM1, but not STIM1, was required for gemcitabine resistance, while high STIM1 levels caused an increase in cytosolic calcium concentration. Higher STIM1-dependent calcium influx led to an impaired endoplasmic reticulum stress response and a heightened nuclear factor of activated T-cell activity. Importantly, these findings were confirmed in patient and patient-derived xenograft samples. Taken together, our study uncovers previously unknown biologically relevant molecular properties of gemcitabine-resistant tumors, revealing an undescribed function of STIM1 as a rheostat directing the effects of calcium signaling and controlling epigenetic cell fate determination. It further reveals the potential benefit of targeting STIM1-controlled calcium signaling and its downstream effectors in PDAC. SIGNIFICANCE: Gemcitabine-resistant and some naïve tumors coamplify RRM1 and STIM1, which elicit gemcitabine resistance and induce a calcium signaling shift, promoting ER stress resistance and activation of NFAT signaling.

Initial responses of the trap-crop, Solanum sisymbriifolium, to Globodera pallida invasions.

Many researchers today are looking for mechanisms underlying plant defenses against nematodes by identifying differentially expressed genes in domesticated hosts. In order to provide a different perspective, we analyzed the root transcriptome of an undomesticated non-host species, Solanum sisymbriifolium Lamark (SSI) before and after Globodera pallida infection. Utilizing RNAseq analyses, we identified changes in the expression of 277 transcripts. Many of these genes were not annotated; however, the annotated set included peroxidases, reactive oxygen species-producing proteins, and regulators of cell death. Importantly, 60% of the nematode-responsive genes did not respond to physical damage to root tissues, or to exogenous treatments with either salicylic acid or methyl jasmonate. Based on this, we speculate that the majority of changes in SSI gene expression were promoted by either nematode effectors, pathogen-associated molecular patterns (PAMPs), or by exposure to untested endogenous signaling molecules such as ethylene, or by exposure to multiple stimuli. This study incorporates our findings into a model that accounts for part of this plant's unusual resistance to nematodes.

An integrative view of the regulatory and transcriptional landscapes in mouse hematopoiesis.

Thousands of epigenomic data sets have been generated in the past decade, but it is difficult for researchers to effectively use all the data relevant to their projects. Systematic integrative analysis can help meet this need, and the VISION project was established for validated systematic integration of epigenomic data in hematopoiesis. Here, we systematically integrated extensive data recording epigenetic features and transcriptomes from many sources, including individual laboratories and consortia, to produce a comprehensive view of the regulatory landscape of differentiating hematopoietic cell types in mouse. By using IDEAS as our integrative and discriminative epigenome annotation system, we identified and assigned epigenetic states simultaneously along chromosomes and across cell types, precisely and comprehensively. Combining nuclease accessibility and epigenetic states produced a set of more than 200,000 candidate cis-regulatory elements (cCREs) that efficiently capture enhancers and promoters. The transitions in epigenetic states of these cCREs across cell types provided insights into mechanisms of regulation, including decreases in numbers of active cCREs during differentiation of most lineages, transitions from poised to active or inactive states, and shifts in nuclease accessibility of CTCF-bound elements. Regression modeling of epigenetic states at cCREs and gene expression produced a versatile resource to improve selection of cCREs potentially regulating target genes. These resources are available from our VISION website to aid research in genomics and hematopoiesis.

Assessment of an Organ-Specific de Novo Transcriptome of the Nematode Trap-Crop, Solanum sisymbriifolium.

Solanum sisymbriifolium, also known as "Litchi Tomato" or "Sticky Nightshade," is an undomesticated and poorly researched plant related to potato and tomato. Unlike the latter species, S. sisymbriifolium induces eggs of the cyst nematode, Globodera pallida, to hatch and migrate into its roots, but then arrests further nematode maturation. In order to provide researchers with a partial blueprint of its genetic make-up so that the mechanism of this response might be identified, we used single molecule real time (SMRT) sequencing to compile a high quality de novo transcriptome of 41,189 unigenes drawn from individually sequenced bud, root, stem, and leaf RNA populations. Functional annotation and BUSCO analysis showed that this transcriptome was surprisingly complete, even though it represented genes expressed at a single time point. By sequencing the 4 organ libraries separately, we found we could get a reliable snapshot of transcript distributions in each organ. A divergent site analysis of the merged transcriptome indicated that this species might have undergone a recent genome duplication and re-diploidization. Further analysis indicated that the plant then retained a disproportionate number of genes associated with photosynthesis and amino acid metabolism in comparison to genes with characteristics of R-proteins or involved in secondary metabolism. The former processes may have given S. sisymbriifolium a bigger competitive advantage than the latter did.