The ENCODE Uniform Analysis Pipelines.

The Encyclopedia of DNA elements (ENCODE) project is a collaborative effort to create a comprehensive catalog of functional elements in the human genome. The current database comprises more than 19000 functional genomics experiments across more than 1000 cell lines and tissues using a wide array of experimental techniques to study the chromatin structure, regulatory and transcriptional landscape of the Homo sapiens and Mus musculus genomes. All experimental data, metadata, and associated computational analyses created by the ENCODE consortium are submitted to the Data Coordination Center (DCC) for validation, tracking, storage, and distribution to community resources and the scientific community. The ENCODE project has engineered and distributed uniform processing pipelines in order to promote data provenance and reproducibility as well as allow interoperability between genomic resources and other consortia. All data files, reference genome versions, software versions, and parameters used by the pipelines are captured and available via the ENCODE Portal. The pipeline code, developed using Docker and Workflow Description Language (WDL; https://openwdl.org/) is publicly available in GitHub, with images available on Dockerhub (https://hub.docker.com), enabling access to a diverse range of biomedical researchers. ENCODE pipelines maintained and used by the DCC can be installed to run on personal computers, local HPC clusters, or in cloud computing environments via Cromwell. Access to the pipelines and data via the cloud allows small labs the ability to use the data or software without access to institutional compute clusters. Standardization of the computational methodologies for analysis and quality control leads to comparable results from different ENCODE collections - a prerequisite for successful integrative analyses.

The ENCODE Uniform Analysis Pipelines.

The Encyclopedia of DNA elements (ENCODE) project is a collaborative effort to create a comprehensive catalog of functional elements in the human genome. The current database comprises more than 19000 functional genomics experiments across more than 1000 cell lines and tissues using a wide array of experimental techniques to study the chromatin structure, regulatory and transcriptional landscape of the Homo sapiens and Mus musculus genomes. All experimental data, metadata, and associated computational analyses created by the ENCODE consortium are submitted to the Data Coordination Center (DCC) for validation, tracking, storage, and distribution to community resources and the scientific community. The ENCODE project has engineered and distributed uniform processing pipelines in order to promote data provenance and reproducibility as well as allow interoperability between genomic resources and other consortia. All data files, reference genome versions, software versions, and parameters used by the pipelines are captured and available via the ENCODE Portal. The pipeline code, developed using Docker and Workflow Description Language (WDL; https://openwdl.org/) is publicly available in GitHub, with images available on Dockerhub (https://hub.docker.com), enabling access to a diverse range of biomedical researchers. ENCODE pipelines maintained and used by the DCC can be installed to run on personal computers, local HPC clusters, or in cloud computing environments via Cromwell. Access to the pipelines and data via the cloud allows small labs the ability to use the data or software without access to institutional compute clusters. Standardization of the computational methodologies for analysis and quality control leads to comparable results from different ENCODE collections - a prerequisite for successful integrative analyses.

Tumor loss-of-function mutations in STK11/LKB1 induce cachexia.

Cancer cachexia (CC), a wasting syndrome of muscle and adipose tissue resulting in weight loss, is observed in 50% of patients with solid tumors. Management of CC is limited by the absence of biomarkers and knowledge of molecules that drive its phenotype. To identify such molecules, we injected 54 human non-small cell lung cancer (NSCLC) lines into immunodeficient mice, 17 of which produced an unambiguous phenotype of cachexia or non-cachexia. Whole-exome sequencing revealed that 8 of 10 cachexia lines, but none of the non-cachexia lines, possessed mutations in serine/threonine kinase 11 (STK11/LKB1), a regulator of nutrient sensor AMPK. Silencing of STK11/LKB1 in human NSCLC and murine colorectal carcinoma lines conferred a cachexia phenotype after cell transplantation into immunodeficient (human NSCLC) and immunocompetent (murine colorectal carcinoma) models. This host wasting was associated with an alteration in the immune cell repertoire of the tumor microenvironments that led to increases in local mRNA expression and serum levels of CC-associated cytokines. Mutational analysis of circulating tumor DNA from patients with NSCLC identified 89% concordance between STK11/LKB1 mutations and weight loss at cancer diagnosis. The current data provide evidence that tumor STK11/LKB1 loss of function is a driver of CC, simultaneously serving as a genetic biomarker for this wasting syndrome.

Novel long non-coding RNAs associated with inflammation and macrophage activation in human.

Inflammation plays a central role in immune response and macrophage activation. Emerging studies demonstrate that along with proteins and genomic factors, noncoding RNA are potentially involved in regulation of immune response and inflammation. Our recent study demonstrated that lncRNA HOTAIR plays key roles in cytokine expression and inflammation in macrophages. The primary goal of this study is to discover novel lncRNAs that are crucial players in inflammation, macrophage activation, and immune response in humans. Towards this, we have stimulated THP1-derived macrophages (THP1-MΦ) with lipopolysaccharides (LPS) and performed the whole transcriptome RNA-seq analysis. Based on this analysis, we discovered that along with well-known marker for inflammation (such as cytokines), a series of long noncoding RNAs (lncRNAs) expression were highly induced upon LPS-stimulation of macrophages, suggesting their potential roles in inflammation and macrophage activation. We termed these family of lncRNAs as Long-noncoding Inflammation Associated RNA (LinfRNA). Dose and time dependent analysis demonstrated that many human LinfRNA (hLinfRNAs) expressions follow similar patterns as cytokine expressions. Inhibition of NF-κB suppressed the expression of most hLinfRNAs suggesting their potential regulation via NF-κB activation during inflammation and macrophage activation. Antisense-mediated knockdown of hLinfRNA1 suppressed the LPS-induced expression of cytokines and pro-inflammatory genes such as IL6, IL1ß, and TNFα expression, suggesting potential functionality of the hLinfRNAs in cytokine regulation and inflammation. Overall, we discovered a series of novel hLinfRNAs that are potential regulators of inflammation and macrophage activation and may be linked to inflammatory and metabolic diseases.

Enhancer and super-enhancer landscape in polycystic kidney disease.

Widespread aberrant gene expression is a pathological hallmark of polycystic kidney disease (PKD). Numerous pathogenic signaling cascades, including c-Myc, Fos, and Jun, are transactivated. However, the underlying epigenetic regulators are poorly defined. Here we show that H3K27ac, an acetylated modification of DNA packing protein histone H3 that marks active enhancers, is elevated in mouse and human samples of autosomal dominant PKD. Using comparative H3K27ac ChIP-Seq analysis, we mapped over 16000 active intronic and intergenic enhancer elements in Pkd1-mutant mouse kidneys. We found that the cystic kidney epigenetic landscape resembles that of a developing kidney, and over 90% of upregulated genes in Pkd1-mutant kidneys are co-housed with activated enhancers in the same topologically associated domains. Furthermore, we identified an evolutionarily conserved enhancer cluster downstream of the c-Myc gene and super-enhancers flanking both Jun and Fos loci in mouse and human models of autosomal dominant PKD. Deleting these regulatory elements reduced c-Myc, Jun, or Fos abundance and suppressed proliferation and 3D cyst growth of Pkd1-mutant cells. Finally, inhibiting glycolysis and glutaminolysis or activating Ppara in Pkd1-mutant cells lowerd global H3K27ac levels and its abundance on c-Myc enhancers. Thus, our work suggests that epigenetic rewiring mediates the transcriptomic dysregulation in PKD, and the regulatory elements can be targeted to slow cyst growth.

Targeting ESR1 mutation-induced transcriptional addiction in breast cancer with BET inhibition.

Acquired mutations in the ligand-binding domain (LBD) of the gene encoding estrogen receptor α (ESR1) are common mechanisms of endocrine therapy resistance in patients with metastatic ER+ breast cancer. The ESR1 Y537S mutation, in particular, is associated with development of resistance to most endocrine therapies used to treat breast cancer. Employing a high-throughput screen of nearly 1,200 Federal Drug Administration-approved (FDA-approved) drugs, we show that OTX015, a bromodomain and extraterminal domain (BET) inhibitor, is one of the top suppressors of ESR1 mutant cell growth. OTX015 was more efficacious than fulvestrant, a selective ER degrader, in inhibiting ESR1 mutant xenograft growth. When combined with abemaciclib, a CDK4/6 inhibitor, OTX015 induced more potent tumor regression than current standard-of-care treatment of abemaciclib + fulvestrant. OTX015 has preferential activity against Y537S mutant breast cancer cells and blocks their clonal selection in competition studies with WT cells. Thus, BET inhibition has the potential to both prevent and overcome ESR1 mutant-induced endocrine therapy resistance in breast cancer.

Epigenetic Repression of STING by MYC Promotes Immune Evasion and Resistance to Immune Checkpoint Inhibitors in Triple-Negative Breast Cancer.

The MYC oncogene is frequently amplified in triple-negative breast cancer (TNBC). Here, we show that MYC suppression induces immune-related hallmark gene set expression and tumor-infiltrating T cells in MYC-hyperactivated TNBCs. Mechanistically, MYC repressed stimulator of interferon genes (STING) expression via direct binding to the STING1 enhancer region, resulting in downregulation of the T-cell chemokines CCL5, CXCL10, and CXCL11. In primary and metastatic TNBC cohorts, tumors with high MYC expression or activity exhibited low STING expression. Using a CRISPR-mediated enhancer perturbation approach, we demonstrated that MYC-driven immune evasion is mediated by STING repression. STING repression induced resistance to PD-L1 blockade in mouse models of TNBC. Finally, a small-molecule inhibitor of MYC combined with PD-L1 blockade elicited a durable response in immune-cold TNBC with high MYC expression, suggesting a strategy to restore PD-L1 inhibitor sensitivity in MYC-overexpressing TNBC.

Decoding Human Genome Regulatory Features That Influence HIV-1 Proviral Expression and Fate Through an Integrated Genomics Approach.

Fundamental principles of HIV-1 integration into the human genome have been revealed in the past 2 decades. However, the impact of the integration site on proviral transcription and expression remains poorly understood. Solving this problem requires the analysis of multiple genomic datasets for thousands of proviral integration sites. Here, we generated and combined large-scale datasets, including epigenetics, transcriptome, and 3-dimensional genome architecture to interrogate the chromatin states, transcription activity, and nuclear sub-compartments around HIV-1 integrations in Jurkat CD4+ T cells to decipher human genome regulatory features shaping the transcription of proviral classes based on their position and orientation in the genome. Through a Hidden Markov Model and ranked informative values prior to a machine learning logistic regression model, we defined nuclear sub-compartments and chromatin states contributing to genomic architecture, transcriptional activity, and nucleosome density of regions neighboring the integration site, as additive features influencing HIV-1 expression. Our integrated genomics approach also allows for a robust experimental design, in which HIV-1 can be genetically introduced into precise genomic locations with known regulatory features to assess the relationship of integration positions to viral transcription and fate.

SeqWho: reliable, rapid determination of sequence file identity using k-mer frequencies in Random Forest classifiers.

MOTIVATION: With the vast improvements in sequencing technologies and increased number of protocols, sequencing is being used to answer complex biological problems. Subsequently, analysis pipelines have become more time consuming and complicated, usually requiring highly extensive prevalidation steps. Here, we present SeqWho, a program designed to assess heuristically the quality of sequencing files and reliably classify the organism and protocol type by using Random Forest classifiers trained on biases native in k-mer frequencies and repeat sequence identities. RESULTS: Using one of our primary models, we show that our method accurately and rapidly classifies human and mouse sequences from nine different sequencing libraries by species, library and both together, 98.32%, 97.86% and 96.38% of the time, respectively. Ultimately, we demonstrate that SeqWho is a powerful method for reliably validating the quality and identity of the sequencing files used in any pipeline. AVAILABILITY AND IMPLEMENTATION: https://github.com/DaehwanKimLab/seqwho. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Glucose starvation induces a switch in the histone acetylome for activation of gluconeogenic and fat metabolism genes.

Acetyl-CoA is a key intermediate situated at the intersection of many metabolic pathways. The reliance of histone acetylation on acetyl-CoA enables the coordination of gene expression with metabolic state. Abundant acetyl-CoA has been linked to the activation of genes involved in cell growth or tumorigenesis through histone acetylation. However, the role of histone acetylation in transcription under low levels of acetyl-CoA remains poorly understood. Here, we use a yeast starvation model to observe the dramatic alteration in the global occupancy of histone acetylation following carbon starvation; the location of histone acetylation marks shifts from growth-promoting genes to gluconeogenic and fat metabolism genes. This reallocation is mediated by both the histone deacetylase Rpd3p and the acetyltransferase Gcn5p, a component of the SAGA transcriptional coactivator. Our findings reveal an unexpected switch in the specificity of histone acetylation to promote pathways that generate acetyl-CoA for oxidation when acetyl-CoA is limiting.

Germline and sporadic mTOR pathway mutations in low-grade oncocytic tumor of the kidney.

Low-grade oncocytic tumor (LOT) of the kidney is a recently described entity with poorly understood pathogenesis. Using next-generation sequencing (NGS) and complementary approaches, we provide insight into its biology. We describe 22 LOT corresponding to 7 patients presenting with a median age of 75 years (range 63-86 years) and male to female ratio 2:5. All 22 tumors demonstrated prototypical microscopic features. Tumors were well-circumscribed and solid. They were composed of sheets of tumor cells in compact nests. Tumor cells had eosinophilic cytoplasm, round to oval nuclei (without nuclear membrane irregularities), focal subtle perinuclear halos, and occasional binucleation. Sharply delineated edematous stromal islands were often observed. Tumor cells were positive for PAX8, negative for CD117, and exhibited diffuse and strong cytokeratin-7 expression. Six patients presented with pT1 tumors. At a median follow-up of 29 months, four patients were alive without recurrence (three patients had died from unrelated causes). All tumors were originally classified as chromophobe renal cell carcinoma, eosinophilic variant (chRCC-eo). While none of the patients presented with known syndromic features, one patient with multiple bilateral LOTs was subsequently found to have a likely pathogenic germline TSC1 mutation. Somatic, likely activating, mutations in MTOR and RHEB were identified in all other evaluable LOTs. As assessed by phospho-S6 and phospho-4E-BP1, mTOR complex 1 (mTORC1) was activated across all cases but to different extent. MTOR mutant LOT exhibited lower levels of mTORC1 activation, possibly related to mTORC1 dimerization and the preservation of a wild-type MTOR copy (retained chromosome 1). Supporting its distinction from related entities, gene expression analyses showed that LOT clustered separately from classic chRCC, chRCC-eo, and RO. In summary, converging mTORC1 pathway mutations, mTORC1 complex activation, and a distinctive gene expression signature along with characteristic phenotypic features support LOT designation as a distinct entity with both syndromic and non-syndromic cases associated with an indolent course.

Improving Renal Tumor Biopsy Prognostication With BAP1 Analyses.

CONTEXT.­: Active surveillance of small renal masses highlights the need for accurate prognostication of biopsies. OBJECTIVE.­: To comprehensively evaluate the accuracy of biopsies in assessing known prognostic parameters including histologic subtype by comparison with subsequent nephrectomy samples. DESIGN.­: We retrospectively identified patients at University of Texas Southwestern Medical Center, Dallas, Texas, who had a biopsy for a renal mass between 2004-2018. Biopsy samples were evaluated for known prognostic factors such as tumor grade, necrosis, sarcomatoid/rhabdoid change, and BRCA1-associated protein-1 (BAP1) status, which we previously showed is an independent prognostic factor for clear cell renal cell carcinoma. Accuracy was determined by comparison with subsequent analyses of nephrectomy specimens. Statistical analyses were performed to assess biopsy accuracy and correlation with tumor size and pathologic stage. RESULTS.­: From 805 biopsies with a diagnosis of renal neoplasm, 178 had subsequent resection of the biopsied tumor. Concordance rate for histologic subtype was 96.9% (κ [w], 0.90; 95% CI, 0.82-0.99) and excellent for small renal masses (98.8%; κ [w], 0.97; 95% CI, 0.90-1). Amongst the prognostic variables evaluated, BAP1 immunohistochemistry in clear cell renal cell carcinoma had the highest agreement (94.8%; κ [w], 0.83; 95% CI, 0.66-0.99). The presence of 1 or more aggressive features (grade 3-4, tumor necrosis, BAP1 loss, sarcomatoid/rhabdoid change) in a biopsy significantly correlated with pT stage (P = .004). CONCLUSIONS.­: Biopsy analyses showed high accuracy for subtyping renal tumors, but it underestimated several poor prognostic features. Addition of BAP1 for clear cell renal cell carcinoma may increase prognostic accuracy. If validated, routine incorporation of BAP1 immunohistochemistry in clear cell renal cell carcinoma biopsies may refine prognosis and aid in the selection of patients for active surveillance.

A multi-omic single-cell landscape of human gynecologic malignancies.

Deconvolution of regulatory mechanisms that drive transcriptional programs in cancer cells is key to understanding tumor biology. Herein, we present matched transcriptome (scRNA-seq) and chromatin accessibility (scATAC-seq) profiles at single-cell resolution from human ovarian and endometrial tumors processed immediately following surgical resection. This dataset reveals the complex cellular heterogeneity of these tumors and enabled us to quantitatively link variation in chromatin accessibility to gene expression. We show that malignant cells acquire previously unannotated regulatory elements to drive hallmark cancer pathways. Moreover, malignant cells from within the same patients show substantial variation in chromatin accessibility linked to transcriptional output, highlighting the importance of intratumoral heterogeneity. Finally, we infer the malignant cell type-specific activity of transcription factors. By defining the regulatory logic of cancer cells, this work reveals an important reliance on oncogenic regulatory elements and highlights the ability of matched scRNA-seq/scATAC-seq to uncover clinically relevant mechanisms of tumorigenesis in gynecologic cancers.

A scalable open-source MATLAB toolbox for reconstruction and analysis of multispectral optoacoustic tomography data.

Multispectral photoacoustic tomography enables the resolution of spectral components of a tissue or sample at high spatiotemporal resolution. With the availability of commercial instruments, the acquisition of data using this modality has become consistent and standardized. However, the analysis of such data is often hampered by opaque processing algorithms, which are challenging to verify and validate from a user perspective. Furthermore, such tools are inflexible, often locking users into a restricted set of processing motifs, which may not be able to accommodate the demands of diverse experiments. To address these needs, we have developed a Reconstruction, Analysis, and Filtering Toolbox to support the analysis of photoacoustic imaging data. The toolbox includes several algorithms to improve the overall quantification of photoacoustic imaging, including non-negative constraints and multispectral filters. We demonstrate various use cases, including dynamic imaging challenges and quantification of drug effect, and describe the ability of the toolbox to be parallelized on a high performance computing cluster.

PARP-1 Regulates Estrogen-Dependent Gene Expression in Estrogen Receptor α-Positive Breast Cancer Cells.

Poly(ADP-ribose) polymerase-1 (PARP-1) has gained considerable attention as a target for therapeutic inhibitors in breast cancers. Previously we showed that PARP-1 localizes to active gene promoters to regulate histone methylation and RNA polymerase II activity (Pol II), altering the expression of various tumor-related genes. Here we report a role for PARP-1 in estrogen-dependent transcription in estrogen receptor alpha (ERα)-positive (ER+) breast cancers. Global nuclear run-on and sequencing analyses functionally linked PARP-1 to the direct control of estrogen-regulated gene expression in ER+ MCF-7 breast cancer cells by promoting transcriptional elongation by Pol II. Furthermore, chromatin immunoprecipitation sequencing analyses revealed that PARP-1 regulates the estrogen-dependent binding of ERα and FoxA1 to a subset of genomic ERα binding sites, promoting active enhancer formation. Moreover, we found that the expression levels of the PARP-1- and estrogen-coregulated gene set are enriched in the luminal subtype of breast cancer, and high PARP-1 expression in ER+ cases correlates with poor survival. Finally, treatment with a PARP inhibitor or a transcriptional elongation inhibitor attenuated estrogen-dependent growth of multiple ER+ breast cancer cell lines. Taken together, our results show that PARP-1 regulates critical molecular pathways that control the estrogen-dependent gene expression program underlying the proliferation of ER+ breast cancer cells. IMPLICATIONS: PARP-1 regulates the estrogen-dependent genomic binding of ERα and FoxA1 to regulate critical gene expression programs by RNA Pol II that underlie the proliferation of ER+ breast cancers, providing a potential therapeutic opportunity for PARP inhibitors in estrogen-responsive breast cancers.

Single-cell analysis of mouse and human prostate reveals novel fibroblasts with specialized distribution and microenvironment interactions.

Stromal-epithelial interactions are critical to the morphogenesis, differentiation, and homeostasis of the prostate, but the molecular identity and anatomy of discrete stromal cell types is poorly understood. Using single-cell RNA sequencing, we identified and validated the in situ localization of three smooth muscle subtypes (prostate smooth muscle, pericytes, and vascular smooth muscle) and two novel fibroblast subtypes in human prostate. Peri-epithelial fibroblasts (APOD+) wrap around epithelial structures, whereas interstitial fibroblasts (C7+) are interspersed in extracellular matrix. In contrast, the mouse displayed three fibroblast subtypes with distinct proximal-distal and lobe-specific distribution patterns. Statistical analysis of mouse and human fibroblasts showed transcriptional correlation between mouse prostate (C3+) and urethral (Lgr5+) fibroblasts and the human interstitial fibroblast subtype. Both urethral fibroblasts (Lgr5+) and ductal fibroblasts (Wnt2+) in the mouse contribute to a proximal Wnt/Tgfb signaling niche that is absent in human prostate. Instead, human peri-epithelial fibroblasts express secreted WNT inhibitors SFRPs and DKK1, which could serve as a buffer against stromal WNT ligands by creating a localized signaling niche around individual prostate glands. We also identified proximal-distal fibroblast density differences in human prostate that could amplify stromal signaling around proximal prostate ducts. In human benign prostatic hyperplasia, fibroblast subtypes upregulate critical immunoregulatory pathways and show distinct distributions in stromal and glandular phenotypes. A detailed taxonomy of leukocytes in benign prostatic hyperplasia reveals an influx of myeloid dendritic cells, T cells and B cells, resembling a mucosal inflammatory disorder. A receptor-ligand interaction analysis of all cell types revealed a central role for fibroblasts in growth factor, morphogen, and chemokine signaling to endothelia, epithelia, and leukocytes. These data are foundational to the development of new therapeutic targets in benign prostatic hyperplasia. © 2021 The Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

IRF1 governs the differential interferon-stimulated gene responses in human monocytes and macrophages by regulating chromatin accessibility.

Myeloid lineage cells use TLRs to recognize and respond to diverse microbial ligands. Although unique transcription factors dictate the outcome of specific TLR signaling, whether lineage-specific differences exist to further modulate the quality of TLR-induced inflammation remains unclear. Comprehensive analysis of global gene transcription in human monocytes, monocyte-derived macrophages, and monocyte-derived dendritic cells stimulated with various TLR ligands identifies multiple lineage-specific, TLR-responsive gene programs. Monocytes are hyperresponsive to TLR7/8 stimulation that correlates with the higher expression of the receptors. While macrophages and monocytes express similar levels of TLR4, macrophages, but not monocytes, upregulate interferon-stimulated genes (ISGs) in response to TLR4 stimulation. We find that TLR4 signaling in macrophages uniquely engages transcription factor IRF1, which facilitates the opening of ISG loci for transcription. This study provides a critical mechanistic basis for lineage-specific TLR responses and uncovers IRF1 as a master regulator for the ISG transcriptional program in human macrophages.

Perivascular mesenchymal cells control adipose-tissue macrophage accrual in obesity.

Chronic low-grade white adipose tissue (WAT) inflammation is a hallmark of metabolic syndrome in obesity. Here, we demonstrate that a subpopulation of mouse WAT perivascular (PDGFRß+) cells, termed fibro-inflammatory progenitors (FIPs), activate proinflammatory signalling cascades shortly after the onset of high-fat diet feeding and regulate proinflammatory macrophage accumulation in WAT in a TLR4-dependent manner. FIPs activation in obesity is mediated by the downregulation of zinc-finger protein 423 (ZFP423), identified here as a transcriptional corepressor of NF-κB. ZFP423 suppresses the DNA-binding capacity of the p65 subunit of NF-κB by inducing a p300-to-NuRD coregulator switch. Doxycycline-inducible expression of Zfp423 in PDGFRß+ cells suppresses inflammatory signalling in FIPs and attenuates metabolic inflammation of visceral WAT in obesity. Inducible inactivation of Zfp423 in PDGFRß+ cells increases FIP activity, exacerbates adipose macrophage accrual and promotes WAT dysfunction. These studies implicate perivascular mesenchymal cells as important regulators of chronic adipose-tissue inflammation in obesity and identify ZFP423 as a transcriptional break on NF-κB signalling.

Acute interstitial nephritis, a potential predictor of response to immune checkpoint inhibitors in renal cell carcinoma.

Immune checkpoint inhibitors (ICIs) such as nivolumab and ipilimumab have improved outcomes in metastatic renal cell carcinoma (mRCC) patients, but they are also associated with immune-related adverse events (irAEs). As observed in melanoma, we hypothesized that patients experiencing an autoimmune reaction directed against the tissue of origin may be more likely to benefit from ICI. Specifically, we asked whether patients with immune-related acute interstitial nephritis (irAIN) exhibited improved outcomes. Using Kidney Cancer Explorer (KCE), a data portal and i2b2-based central database for clinical, pathological and experimental genetic data, we systematically identified all patients with mRCC at UT Southwestern Medical Center (UTSW) from 2014-2018 who received at least one dose of ICI. More recent cases were identified through a provider query. We extracted creatinine (Cr) values at baseline and over the entirety of each patient ICI treatment course using KCE. Patients with ≥ 1.5-fold Cr increase over baseline were investigated. The likelihood of irAIN was determined based on the work-up (biopsy, if available), or by clinical criteria (timing of kidney injury, exclusion of other etiologies, treatment with immunosuppressants and response). We identified 177 mRCC patients who received at least one dose of ICI, 36 of whom had ≥ 1.5-fold increase in Cr over baseline while on treatment. Of those, two had biopsy-proven irAIN and one was clinically diagnosed, resulting in an incidence of 1.7%. One additional biopsy-proven case past 2018 was identified through a provider query, for a total of four patients. Two received combination nivolumab and ipilimumab in the first line, whereas the remaining received nivolumab after first line therapy. irAIN onset ranged from 1.5 to 12 months. All four patients stopped ICI with recovery of renal function, at least partially, three after receiving systemic steroids. Notably, all four patients had a deep response. In conclusion, irAIN is a rare event, but it may portend a higher likelihood of response. One possible explanation is antigenic overlap between normal renal tubular cells and tumor cells.

Total Functional Score of Enhancer Elements Identifies Lineage-Specific Enhancers That Drive Differentiation of Pancreatic Cells.

The differentiation of embryonic stem cells into various lineages is highly dependent on the chromatin state of the genome and patterns of gene expression. To identify lineage-specific enhancers driving the differentiation of progenitors into pancreatic cells, we used a previously described computational framework called Total Functional Score of Enhancer Elements (TFSEE), which integrates multiple genomic assays that probe both transcriptional and epigenomic states. First, we evaluated and compared TFSEE as an enhancer-calling algorithm with enhancers called using GRO-seq-defined enhancer transcripts (method 1) versus enhancers called using histone modification ChIP-seq data (method 2). Second, we used TFSEE to define the enhancer landscape and identify transcription factors (TFs) that maintain the multipotency of a subpopulation of endodermal stem cells during differentiation into pancreatic lineages. Collectively, our results demonstrate that TFSEE is a robust enhancer-calling algorithm that can be used to perform multilayer genomic data integration to uncover cell type-specific TFs that control lineage-specific enhancers.