CSI NGS Portal: An Online Platform for Automated NGS Data Analysis and Sharing

Methylation Mesa define functional regulatory elements for targeted gene activation

DNA methylation and mRNA expression correlations are often presented with inconsistent evidence supporting causal regulation. We hypothesized that causal regulatory methylation elements would exhibit heightened demethylation sensitivity. To investigate, we analyzed 20 whole-genomic bisulfite sequenced samples before and after demethylation and identified narrow-width (45-294 bp) elements within a short plateau, termed Methylation Mesa (MM). The Mesa signature was conserved across species and was independent of CpG islands. Mesa also demonstrate high concordance with primed and active histone marks. To assess causality, we developed CRISPR-DiR, a highly precise targeted demethylation technology. Targeted demethylation of a Mesa triggers locus and distal chromatin rewiring events that initiate mRNA expression significantly greater than promoter-CpG island targeting. Thus, Mesa are self-sustaining epigenetic regulatory elements that maintain long-term gene activation through focused demethylation only within the Mesa core, resulting in subsequent histone modifications and chromatin rewiring events that interact with distal elements also marked as Mesas.

A Review of Web-Based Metagenomics Platforms for Analysing Next-Generation Sequence Data

Metagenomics has now evolved as a promising technology for understanding the microbial population in the environment. By metagenomics, a number of extreme and complex environment has been explored for their microbial population. Using this technology, researchers have brought out novel genes and their potential characteristics, which have robust applications in food, pharmaceutical, scientific research, and other biotechnological fields. A sequencing platform can provide a sequence of microbial populations in any given environment. The sequence needs to be analysed computationally to derive meaningful information. It is presumed that only bioinformaticians with extensive computational skills can process the sequencing data till the downstream end. However, numerous open-source software and online servers are available to analyse the metagenomic data developed for a biologist with less computational skills. This review is focused on bioinformatics tools such as Galaxy, CSI-NGS portal, ANASTASIA and SHAMAN, EBI- metagenomics, IDseq, and MG-RAST for analysing metagenomic data.

FEZ1 participates in human embryonic brain development by modulating neuronal progenitor subpopulation specification and migrations

Mutations in the human fasciculation and elongation protein zeta 1 (FEZ1) gene are found in schizophrenia and Jacobsen syndrome patients. Here, using human cerebral organoids (hCOs), we show that FEZ1 expression is turned on early during brain development and is detectable in both neuroprogenitor subtypes and immature neurons. FEZ1 deletion disrupts expression of neuronal and synaptic development genes. Using single-cell RNA sequencing, we detected abnormal expansion of homeodomain-only protein homeobox (HOPX)- outer radial glia (oRG), concurrent with a reduction of HOPX+ oRG, in FEZ1-null hCOs. HOPX- oRGs show higher cell mobility as compared to HOPX+ oRGs. Ectopic localization of neuroprogenitors to the outer layer is seen in FEZ1-null hCOs. Anomalous encroachment of TBR2+ intermediate progenitors into CTIP2+ deep layer neurons further indicated abnormalities in cortical layer formation these hCOs. Collectively, our findings highlight the involvement of FEZ1 in early cortical brain development and how it contributes to neurodevelopmental disorders.

JAK2 as a surface marker for enrichment of human pluripotent stem cells-derived ventricular cardiomyocytes

BACKGROUND

Human pluripotent stem cell (hPSC)-derived cardiomyocytes (CMs) hold great promise for cardiac disease modelling, drug discovery and regenerative medicine. Despite the advancement in various differentiation protocols, the heterogeneity of the generated population composed of diverse cardiac subtypes poses a significant challenge to their practical applications. Mixed populations of cardiac subtypes can compromise disease modelling and drug discovery, while transplanting them may lead to undesired arrhythmias as they may not integrate and synchronize with the host tissue's contractility. It is therefore crucial to identify cell surface markers that could enable high purity of ventricular CMs for subsequent applications.

METHODS

By exploiting the fact that immature CMs expressing myosin light chain 2A (MLC2A) will gradually express myosin light chain 2 V (MLC2V) protein as they mature towards ventricular fate, we isolated signal regulatory protein alpha (SIRPA)-positive CMs expressing intracellular MLC2A or MLC2V using MARIS (method for analysing RNA following intracellular sorting). Subsequently, RNA sequencing analysis was performed to examine the gene expression profile of MLC2A + and MLC2V + sorted CMs. We identified genes that were significantly up-regulated in MLC2V + samples to be potential surface marker candidates for ventricular specification. To validate these surface markers, we performed immunostaining and western blot analysis to measure MLC2A and MLC2V protein expressions in SIRPA + CMs that were either positive or negative for the putative surface markers, JAK2 (Janus kinase 2) or CD200. We then characterized the electrophysiological properties of surface marker-sorted CMs, using fluo-4 AM, a green-fluorescent calcium indicator, to measure the cellular calcium transient at the single cell level. For functional validation, we investigated the response of the surface marker-sorted CMs to vernakalant, an atrial-selective anti-arrhythmic agent.

RESULTS

In this study, while JAK2 and CD200 were identified as potential surface markers for the purification of ventricular-like CMs, the SIRPA+/JAK2+ population showed a higher percentage of MLC2V-expressing cells (~ 90%) compared to SIRPA+/CD200+ population (~ 75%). SIRPA+/JAK2+ sorted CMs exhibited ventricular-like electrophysiological properties, including slower beating rate, slower calcium depolarization and longer calcium repolarization duration. Importantly, vernakalant had limited to no significant effect on the calcium repolarization duration of SIRPA+/JAK2+ population, indicating their enrichment for ventricular-like CMs.

CONCLUSION

Our study lays the groundwork for the identification of cardiac subtype surface markers that allow purification of cardiomyocyte sub-populations. Our findings suggest that JAK2 can be employed as a cell surface marker for enrichment of hPSC-derived ventricular-like CMs.

RUNX3 inactivates oncogenic MYC through disruption of MYC/MAX complex and subsequent recruitment of GSK3β-FBXW7 cascade

MYC is one of the most commonly dysregulated proto-oncogenes in cancer. MYC promotes cancer initiation and maintenance by regulating multiple biological processes, such as proliferation and stem cell function. Here, we show that developmental regulator RUNX3 targets MYC protein for rapid degradation through the glycogen synthase kinase-3 beta-F-box/WD repeat-containing protein 7 (GSK3β-FBXW7) proteolytic pathway. The evolutionarily conserved Runt domain of RUNX3 interacts directly with the basic helix-loop-helix leucine zipper of MYC, resulting in the disruption of MYC/MAX and MYC/MIZ-1 interactions, enhanced GSK3β-mediated phosphorylation of MYC protein at threonine-58 and its subsequent degradation via the ubiquitin-proteasomal pathway. We therefore uncover a previously unknown mode of MYC destabilization by RUNX3 and provide an explanation as to why RUNX3 inhibits early-stage cancer development in gastrointestinal and lung mouse cancer models.

Resources and tools for rare disease variant interpretation

Collectively, rare genetic disorders affect a substantial portion of the world's population. In most cases, those affected face difficulties in receiving a clinical diagnosis and genetic characterization. The understanding of the molecular mechanisms of these diseases and the development of therapeutic treatments for patients are also challenging. However, the application of recent advancements in genome sequencing/analysis technologies and computer-aided tools for predicting phenotype-genotype associations can bring significant benefits to this field. In this review, we highlight the most relevant online resources and computational tools for genome interpretation that can enhance the diagnosis, clinical management, and development of treatments for rare disorders. Our focus is on resources for interpreting single nucleotide variants. Additionally, we present use cases for interpreting genetic variants in clinical settings and review the limitations of these results and prediction tools. Finally, we have compiled a curated set of core resources and tools for analyzing rare disease genomes. Such resources and tools can be utilized to develop standardized protocols that will enhance the accuracy and effectiveness of rare disease diagnosis.

Patterns of Oncogene Coexpression at Single-Cell Resolution Influence Survival in Lymphoma

Cancers often overexpress multiple clinically relevant oncogenes, but it is not known if combinations of oncogenes in cellular subpopulations within a cancer influence clinical outcomes. Using quantitative multispectral imaging of the prognostically relevant oncogenes MYC, BCL2, and BCL6 in diffuse large B-cell lymphoma (DLBCL), we show that the percentage of cells with a unique combination MYC+BCL2+BCL6- (M+2+6-) consistently predicts survival across four independent cohorts (n = 449), an effect not observed with other combinations including M+2+6+. We show that the M+2+6- percentage can be mathematically derived from quantitative measurements of the individual oncogenes and correlates with survival in IHC (n = 316) and gene expression (n = 2,521) datasets. Comparative bulk/single-cell transcriptomic analyses of DLBCL samples and MYC/BCL2/BCL6-transformed primary B cells identify molecular features, including cyclin D2 and PI3K/AKT as candidate regulators of M+2+6- unfavorable biology. Similar analyses evaluating oncogenic combinations at single-cell resolution in other cancers may facilitate an understanding of cancer evolution and therapy resistance.

SIGNIFICANCE

Using single-cell-resolved multiplexed imaging, we show that selected subpopulations of cells expressing specific combinations of oncogenes influence clinical outcomes in lymphoma. We describe a probabilistic metric for the estimation of cellular oncogenic coexpression from IHC or bulk transcriptomes, with possible implications for prognostication and therapeutic target discovery in cancer. This article is highlighted in the In This Issue feature, p. 1027.

Systematic benchmarking of statistical methods to assess differential expression of circular RNAs

Circular RNAs (circRNAs) are covalently closed transcripts involved in critical regulatory axes, cancer pathways and disease mechanisms. CircRNA expression measured with RNA-seq has particular characteristics that might hamper the performance of standard biostatistical differential expression assessment methods (DEMs). We compared 38 DEM pipelines configured to fit circRNA expression data's statistical properties, including bulk RNA-seq, single-cell RNA-seq (scRNA-seq) and metagenomics DEMs. The DEMs performed poorly on data sets of typical size. Widely used DEMs, such as DESeq2, edgeR and Limma-Voom, gave scarce results, unreliable predictions or even contravened the expected behaviour with some parameter configurations. Limma-Voom achieved the most consistent performance throughout different benchmark data sets and, as well as SAMseq, reasonably balanced false discovery rate (FDR) and recall rate. Interestingly, a few scRNA-seq DEMs obtained results comparable with the best-performing bulk RNA-seq tools. Almost all DEMs' performance improved when increasing the number of replicates. CircRNA expression studies require careful design, choice of DEM and DEM configuration. This analysis can guide scientists in selecting the appropriate tools to investigate circRNA differential expression with RNA-seq experiments.

Paracrine Senescence of Mesenchymal Stromal Cells Involves Inflammatory Cytokines and the NF-κB Pathway

It has been known that senescence-associated secretory phenotype (SASP) triggers senescence of the surrounding normal cells. However, SASP signaling regarding mesenchymal stromal cell aging remains to be fully elucidated. Therefore, the present study aimed to clarify the molecular mechanism of late (passage) MSC-induced paracrine SASP-mediated senescence of early (passage) MSCs during ex vivo expansion. Here, we conducted an extensive characterization of senescence features in bone-marrow (BM)-derived MSCs from healthy human donors. Late MSCs displayed an enlarged senescent-like morphology, induced SASP-related proinflammatory cytokines (IL-1α and IL-8), and reduced clonogenic capacity and osteogenic differentiation when compared to early MSCs. Of note, paracrine effects of SASP-related IL-1α and IL-8 from late MSCs induced cellular senescence of early MSCs via an NF-κB-dependent manner. Moreover, cellular senescence of early MSCs was promoted by the synergistic action of IL-1α and IL-8. However, inhibition of NF-κB by shRNA transfection or using inhibitors in early MSCs blocked early MSCs cellular senescence caused by paracrine SASP of late MSCs. In conclusion, these findings reveal that late MSCs display features of senescence and that, during ex vivo expansion, SASP-related proinflammatory cytokines contribute to activate a cellular senescence program in early MSCs that may ultimately impair their functionality.

Multilayered control of splicing regulatory networks by DAP3 leads to widespread alternative splicing changes in cancer

The dynamic regulation of alternative splicing requires coordinated participation of multiple RNA binding proteins (RBPs). Aberrant splicing caused by dysregulation of splicing regulatory RBPs is implicated in numerous cancers. Here, we reveal a frequently overexpressed cancer-associated protein, DAP3, as a splicing regulatory RBP in cancer. Mechanistically, DAP3 coordinates splicing regulatory networks, not only via mediating the formation of ribonucleoprotein complexes to induce substrate-specific splicing changes, but also via modulating splicing of numerous splicing factors to cause indirect effect on splicing. A pan-cancer analysis of alternative splicing across 33 TCGA cancer types identified DAP3-modulated mis-splicing events in multiple cancers, and some of which predict poor prognosis. Functional investigation of non-productive splicing of WSB1 provides evidence for establishing a causal relationship between DAP3-modulated mis-splicing and tumorigenesis. Together, our work provides critical mechanistic insights into the splicing regulatory roles of DAP3 in cancer development.

RNA binding proteins (RBPs) can participate in regulatory networks to control alternative splicing. Here the authors show that DAP3 functions as an RBP splicing modulator via two mechanisms, and that its overexpression leads to mis-splicing events in cancers.

MYC overexpression leads to increased chromatin interactions at superenhancers and MYC binding sites

The MYC oncogene encodes for the MYC protein and is frequently dysregulated across multiple cancer cell types, making it an attractive target for cancer therapy. MYC overexpression leads to MYC binding at active enhancers, resulting in a global transcriptional amplification of active genes. Because super-enhancers are frequently dysregulated in cancer, we hypothesized that MYC preferentially invades into super-enhancers and alters the cancer genome organization. To that end, we performed ChIP-seq, RNA-seq, circular chromosome conformation capture (4C-seq), and Spike-in Quantitative Hi-C (SIQHiC) on the U2OS osteosarcoma cell line with tetracycline-inducible MYC. MYC overexpression in U2OS cells modulated histone acetylation and increased MYC binding at super-enhancers. SIQHiC analysis revealed increased global chromatin contact frequency, particularly at chromatin interactions connecting MYC binding sites at promoters and enhancers. Immunofluorescence staining showed that MYC molecules formed punctate foci at these transcriptionally active domains after MYC overexpression. These results demonstrate the accumulation of overexpressed MYC at promoter–enhancer hubs and suggest that MYC invades into enhancers through spatial proximity. At the same time, the increased protein–protein interactions may strengthen these chromatin interactions to increase chromatin contact frequency. CTCF siRNA knockdown in MYC-overexpressed U2OS cells demonstrated that removal of architectural proteins can disperse MYC and abrogate the increase in chromatin contacts. By elucidating the chromatin landscape of MYC-driven cancers, we can potentially target MYC-associated chromatin interactions for cancer therapy.

The emerging roles of NGS in clinical oncology and personalized medicine

Next-generation sequencing (NGS) has been increasingly popular in genomics studies over the last decade, as new sequencing technology has been created and improved. Recently, NGS started to be used in clinical oncology to improve cancer therapy through diverse modalities ranging from finding novel and rare cancer mutations, discovering cancer mutation carriers to reaching specific therapeutic approaches known as personalized medicine (PM). PM has the potential to minimize medical expenses by shifting the current traditional medical approach of treating cancer and other diseases to an individualized preventive and predictive approach. Currently, NGS can speed up in the early diagnosis of diseases and discover pharmacogenetic markers that help in personalizing therapies. Despite the tremendous growth in our understanding of genetics, NGS holds the added advantage of providing more comprehensive picture of cancer landscape and uncovering cancer development pathways. In this review, we provided a complete overview of potential NGS applications in scientific and clinical oncology, with a particular emphasis on pharmacogenomics in the direction of precision medicine treatment options.

ChIP-AP: an integrated analysis pipeline for unbiased ChIP-seq analysis

Chromatin immunoprecipitation coupled with sequencing (ChIP-seq) is a technique used to identify protein–DNA interaction sites through antibody pull-down, sequencing and analysis; with enrichment ‘peak’ calling being the most critical analytical step. Benchmarking studies have consistently shown that peak callers have distinct selectivity and specificity characteristics that are not additive and seldom completely overlap in many scenarios, even after parameter optimization. We therefore developed ChIP-AP, an integrated ChIP-seq analysis pipeline utilizing four independent peak callers, which seamlessly processes raw sequencing files to final result. This approach enables (1) better gauging of peak confidence through detection by multiple algorithms, and (2) more thoroughly surveys the binding landscape by capturing peaks not detected by individual callers. Final analysis results are then integrated into a single output table, enabling users to explore their data by applying selectivity and sensitivity thresholds that best address their biological questions, without needing any additional reprocessing. ChIP-AP therefore presents investigators with a more comprehensive coverage of the binding landscape without requiring additional wet-lab observations.

Targeting an Inducible SALL4-Mediated Cancer Vulnerability with Sequential Therapy

Oncofetal protein SALL4 is critical for cancer cell survival. Targeting SALL4, however, is only applicable in a fraction of cancer patients who are positive for this gene. To overcome this limitation, we propose to induce a cancer vulnerability by engineering a partial dependency upon SALL4. Following exogenous expression of SALL4, SALL4-negative cancer cells became partially dependent on SALL4. Treatment of SALL4-negative cells with the FDA-approved hypomethylating agent 5-aza-2'-deoxycytidine (DAC) resulted in transient upregulation of SALL4. DAC pretreatment sensitized SALL4-negative cancer cells to entinostat, which negatively affected SALL4 expression through a microRNA, miRNA-205, both in culture and in vivo. Moreover, SALL4 was essential for the efficiency of sequential treatment of DAC and entinostat. Overall, this proof-of-concept study provides a framework whereby the targeting pathways such as SALL4-centered therapy can be expanded, sensitizing cancer cells to treatment by transient target induction and engineering a dependency. SIGNIFICANCE: These findings provide a therapeutic approach for patients harboring no suitable target by induction of a SALL4-mediated vulnerability.

Chronic Benzene Exposure Aggravates Pressure Overload-Induced Cardiac Dysfunction

Benzene is a ubiquitous environmental pollutant abundant in household products, petrochemicals and cigarette smoke. Benzene is a well-known carcinogen in humans and experimental animals; however, little is known about the cardiovascular toxicity of benzene. Recent population-based studies indicate that benzene exposure is associated with an increased risk for heart failure. Nonetheless, it is unclear whether benzene exposure is sufficient to induce and/or exacerbate heart failure. We examined the effects of benzene (50 ppm, 6 h/day, 5 days/week, 6 weeks) or HEPA-filtered air exposure on transverse aortic constriction (TAC)-induced pressure overload in male C57BL/6J mice. Our data show that benzene exposure had no effect on cardiac function in the Sham group; however, it significantly compromised cardiac function as depicted by a significant decrease in fractional shortening and ejection fraction, as compared with TAC/Air-exposed mice. RNA-seq analysis of the cardiac tissue from the TAC/benzene-exposed mice showed a significant increase in several genes associated with adhesion molecules, cell-cell adhesion, inflammation, and stress response. In particular, neutrophils were implicated in our unbiased analyses. Indeed, immunofluorescence studies showed that TAC/benzene exposure promotes infiltration of CD11b+/S100A8+/myeloperoxidase+-positive neutrophils in the hearts by 3-fold. In vitro, the benzene metabolites, hydroquinone and catechol, induced the expression of P-selectin in cardiac microvascular endothelial cells by 5-fold and increased the adhesion of neutrophils to these endothelial cells by 1.5-2.0-fold. Benzene metabolite-induced adhesion of neutrophils to the endothelial cells was attenuated by anti-P-selectin antibody. Together, these data suggest that benzene exacerbates heart failure by promoting endothelial activation and neutrophil recruitment.

Quantitative imaging of RAD51 expression as a marker of platinum resistance in ovarian cancer

Early relapse after platinum chemotherapy in epithelial ovarian cancer (EOC) portends poor survival. A-priori identification of platinum resistance is therefore crucial to improve on standard first-line carboplatin-paclitaxel treatment. The DNA repair pathway homologous recombination (HR) repairs platinum-induced damage, and the HR recombinase RAD51 is overexpressed in cancer. We therefore designed a REMARK-compliant study of pre-treatment RAD51 expression in EOC, using fluorescent quantitative immunohistochemistry (qIHC) to overcome challenges in quantitation of protein expression in situ. In a discovery cohort (n = 284), RAD51-High tumours had shorter progression-free and overall survival compared to RAD51-Low cases in univariate and multivariate analyses. The association of RAD51 with relapse/survival was validated in a carboplatin monotherapy SCOTROC4 clinical trial cohort (n = 264) and was predominantly noted in HR-proficient cancers (Myriad HRDscore < 42). Interestingly, overexpression of RAD51 modified expression of immune-regulatory pathways in vitro, while RAD51-High tumours showed exclusion of cytotoxic T cells in situ. Our findings highlight RAD51 expression as a determinant of platinum resistance and suggest possible roles for therapy to overcome immune exclusion in RAD51-High EOC. The qIHC approach is generalizable to other proteins with a continuum instead of discrete/bimodal expression.

Suppression of adenosine-to-inosine (A-to-I) RNA editome by death associated protein 3 (DAP3) promotes cancer progression

RNA editing introduces nucleotide changes in RNA sequences. Recent studies have reported that aberrant A-to-I RNA editing profiles are implicated in cancers. Albeit changes in expression and activity of ADAR genes are thought to have been responsible for the dysregulated RNA editome in diseases, they are not always correlated, indicating the involvement of secondary regulators. Here, we uncover DAP3 as a potent repressor of editing and a strong oncogene in cancer. DAP3 mainly interacts with the deaminase domain of ADAR2 and represses editing via disrupting association of ADAR2 with its target transcripts. PDZD7, an exemplary DAP3-repressed editing target, undergoes a protein recoding editing at stop codon [Stop →Trp (W)]. Because of editing suppression by DAP3, the unedited PDZD7WT, which is more tumorigenic than edited PDZD7Stop518W, is accumulated in tumors. In sum, cancer cells may acquire malignant properties for their survival advantage through suppressing RNA editome by DAP3.