Identifying novel regulatory effects for clinically relevant genes through the study of the Greek population.

BACKGROUND: Expression quantitative trait loci (eQTL) studies provide insights into regulatory mechanisms underlying disease risk. Expanding studies of gene regulation to underexplored populations and to medically relevant tissues offers potential to reveal yet unknown regulatory variants and to better understand disease mechanisms. Here, we performed eQTL mapping in subcutaneous (S) and visceral (V) adipose tissue from 106 Greek individuals (Greek Metabolic study, GM) and compared our findings to those from the Genotype-Tissue Expression (GTEx) resource. RESULTS: We identified 1,930 and 1,515 eGenes in S and V respectively, over 13% of which are not observed in GTEx adipose tissue, and that do not arise due to different ancestry. We report additional context-specific regulatory effects in genes of clinical interest (e.g. oncogene ST7) and in genes regulating responses to environmental stimuli (e.g. MIR21, SNX33). We suggest that a fraction of the reported differences across populations is due to environmental effects on gene expression, driving context-specific eQTLs, and suggest that environmental effects can determine the penetrance of disease variants thus shaping disease risk. We report that over half of GM eQTLs colocalize with GWAS SNPs and of these colocalizations 41% are not detected in GTEx. We also highlight the clinical relevance of S adipose tissue by revealing that inflammatory processes are upregulated in individuals with obesity, not only in V, but also in S tissue. CONCLUSIONS: By focusing on an understudied population, our results provide further candidate genes for investigation regarding their role in adipose tissue biology and their contribution to disease risk and pathogenesis.

In Vitro Fermentation of Pleurotus eryngii Mushrooms by Human Fecal Microbiota: Metataxonomic Analysis and Metabolomic Profiling of Fermentation Products.

Edible mushrooms contain biologically active compounds with antioxidant, antimicrobial, immunomodulatory and anticancer properties. The link between their anticancer and immunomodulatory properties with their possible prebiotic activity on gut micro-organisms has been the subject of intense research over the last decade. Lyophilized Pleurotus eryngii (PE) mushrooms, selected due to their strong lactogenic effect and anti-genotoxic, immunomodulatory properties, underwent in vitro static batch fermentation for 24 h by fecal microbiota from eight elderly apparently healthy volunteers (>65 years old). The fermentation-induced changes in fecal microbiota communities were examined using Next Generation Sequencing of the hypervariable regions of the 16S rRNA gene. Primary processing and analysis were conducted using the Ion Reporter Suite. Changes in the global metabolic profile were assessed by 1H NMR spectroscopy, and metabolites were assigned by 2D NMR spectroscopy and the MetaboMiner platform. PLS-DA analysis of both metataxonomic and metabolomic data showed a significant cluster separation of PE fermented samples relative to controls. DEseq2 analysis showed that the abundance of families such as Lactobacillaceae and Bifidobacteriaceae were increased in PE samples. Accordingly, in metabolomics, more than twenty metabolites including SCFAs, essential amino acids, and neurotransmitters discriminate PE samples from the respective controls, further validating the metataxonomic findings.

Interactive Analysis, Exploration, and Visualization of RNA-Seq Data with SeqCVIBE.

The rise of modern gene expression profiling techniques, such as RNA-Seq, has generated a wealth of high-quality datasets spanning all fields of current biological research. The large data sets and the continually expanding applications for which they can be mined, such as the investigation of alternative splicing and others, have created novel challenges for data management, exploration, analysis, and visualization. Although a large variety of RNA-Seq data analysis software packages has emerged, both open-source and commercial, most fail to simultaneously address the above challenges, while they lack obvious functionalities, such as estimating RNA abundance over non-annotated genomic regions of interest in real time. We have developed SeqCVIBE, an R Shiny web application for the interactive exploration, analysis, visualization, and genome browsing of large RNA-Seq datasets. SeqCVIBE allows for multiple on-the-fly visualizations and calculations, such as differential expression analysis, averaging genomic signals over specific regions of the genome, and calculating RNA abundances over custom, potentially non-annotated regions, such as novel long non-coding RNAs. In addition, SeqCVIBE comprises a database for pre-analyzed data, where users can navigate and explore results, as well as perform a variety of basic on-the-fly analyses and export the outcomes. Finally, we demonstrate the value of SeqCVIBE in the elucidation of the interplay of a novel lincRNA, WiNTRLINC1, and Wnt signaling in colon cancer.

The transcription factor Xrp1 orchestrates both reduced translation and cell competition upon defective ribosome assembly or function.

Ribosomal Protein (Rp) gene haploinsufficiency affects translation rate, can lead to protein aggregation, and causes cell elimination by competition with wild type cells in mosaic tissues. We find that the modest changes in ribosomal subunit levels observed were insufficient for these effects, which all depended on the AT-hook, bZip domain protein Xrp1. Xrp1 reduced global translation through PERK-dependent phosphorylation of eIF2α. eIF2α phosphorylation was itself sufficient to enable cell competition of otherwise wild type cells, but through Xrp1 expression, not as the downstream effector of Xrp1. Unexpectedly, many other defects reducing ribosome biogenesis or function (depletion of TAF1B, eIF2, eIF4G, eIF6, eEF2, eEF1α1, or eIF5A), also increased eIF2α phosphorylation and enabled cell competition. This was also through the Xrp1 expression that was induced in these depletions. In the absence of Xrp1, translation differences between cells were not themselves sufficient to trigger cell competition. Xrp1 is shown here to be a sequence-specific transcription factor that regulates transposable elements as well as single-copy genes. Thus, Xrp1 is the master regulator that triggers multiple consequences of ribosomal stresses and is the key instigator of cell competition.

Increased Autotaxin Levels in Severe COVID-19, Correlating with IL-6 Levels, Endothelial Dysfunction Biomarkers, and Impaired Functions of Dendritic Cells.

Autotaxin (ATX; ENPP2) is a secreted lysophospholipase D catalyzing the extracellular production of lysophosphatidic acid (LPA), a pleiotropic signaling phospholipid. Genetic and pharmacologic studies have previously established a pathologic role for ATX and LPA signaling in pulmonary injury, inflammation, and fibrosis. Here, increased ENPP2 mRNA levels were detected in immune cells from nasopharyngeal swab samples of COVID-19 patients, and increased ATX serum levels were found in severe COVID-19 patients. ATX serum levels correlated with the corresponding increased serum levels of IL-6 and endothelial damage biomarkers, suggesting an interplay of the ATX/LPA axis with hyperinflammation and the associated vascular dysfunction in COVID-19. Accordingly, dexamethasone (Dex) treatment of mechanically ventilated patients reduced ATX levels, as shown in two independent cohorts, indicating that the therapeutic benefits of Dex include the suppression of ATX. Moreover, large scale analysis of multiple single cell RNA sequencing datasets revealed the expression landscape of ENPP2 in COVID-19 and further suggested a role for ATX in the homeostasis of dendritic cells, which exhibit both numerical and functional deficits in COVID-19. Therefore, ATX has likely a multifunctional role in COVID-19 pathogenesis, suggesting that its pharmacological targeting might represent an additional therapeutic option, both during and after hospitalization.

Application of human liver organoids as a patient-derived primary model for HBV infection and related hepatocellular carcinoma.

The molecular events that drive hepatitis B virus (HBV)-mediated transformation and tumorigenesis have remained largely unclear, due to the absence of a relevant primary model system. Here we propose the use of human liver organoids as a platform for modeling HBV infection and related tumorigenesis. We first describe a primary ex vivo HBV-infection model derived from healthy donor liver organoids after challenge with recombinant virus or HBV-infected patient serum. HBV-infected organoids produced covalently closed circular DNA (cccDNA) and HBV early antigen (HBeAg), expressed intracellular HBV RNA and proteins, and produced infectious HBV. This ex vivo HBV-infected primary differentiated hepatocyte organoid platform was amenable to drug screening for both anti-HBV activity and drug-induced toxicity. We also studied HBV replication in transgenically modified organoids; liver organoids exogenously overexpressing the HBV receptor sodium taurocholate co-transporting polypeptide (NTCP) after lentiviral transduction were not more susceptible to HBV, suggesting the necessity for additional host factors for efficient infection. We also generated transgenic organoids harboring integrated HBV, representing a long-term culture system also suitable for viral production and the study of HBV transcription. Finally, we generated HBV-infected patient-derived liver organoids from non-tumor cirrhotic tissue of explants from liver transplant patients. Interestingly, transcriptomic analysis of patient-derived liver organoids indicated the presence of an aberrant early cancer gene signature, which clustered with the hepatocellular carcinoma (HCC) cohort on The Cancer Genome Atlas Liver Hepatocellular Carcinoma dataset and away from healthy liver tissue, and may provide invaluable novel biomarkers for the development of HCC and surveillance in HBV-infected patients.

The transcription factor BCL-6 controls early development of innate-like T cells.

Innate T cells, including invariant natural killer T (iNKT) and mucosal-associated innate T (MAIT) cells, are a heterogeneous T lymphocyte population with effector properties preprogrammed during their thymic differentiation. How this program is initiated is currently unclear. Here, we show that the transcription factor BCL-6 was transiently expressed in iNKT cells upon exit from positive selection and was required for their proper development beyond stage 0. Notably, development of MAIT cells was also impaired in the absence of Bcl6. BCL-6-deficient iNKT cells had reduced expression of genes that were associated with the innate T cell lineage, including Zbtb16, which encodes PLZF, and PLZF-targeted genes. BCL-6 contributed to a chromatin accessibility landscape that was permissive for the expression of development-related genes and inhibitory for genes associated with naive T cell programs. Our results revealed new functions for BCL-6 and illuminated how this transcription factor controls early iNKT cell development.

Refining a steroidogenic model: an analysis of RNA-seq datasets from insect prothoracic glands.

BACKGROUND: The prothoracic gland (PG), the principal steroidogenic organ of insects, has been proposed as a model for steroid hormone biosynthesis and regulation. RESULTS: To validate the robustness of the model, we present an analysis of accumulated transcriptomic data from PGs of two model species, Drosophila melanogaster and Bombyx mori. We identify that the common core components of the model in both species are encoded by nine genes. Five of these are Halloween genes whose expression differs substantially between the PGs of these species. CONCLUSIONS: We conclude that the PGs can be a model for steroid hormone synthesis and regulation within the context of mitochondrial cholesterol transport and steroid biosynthesis but beyond these core mechanisms, gene expression in insect PGs is too diverse to fit in a context-specific model and should be analysed within a species-specific framework.

Smyd3-associated regulatory pathways in cancer.

SMYD3 is a member of the SET and MYND-domain family of methyl-transferases, the increased expression of which correlates with poor prognosis in various types of cancer. In liver and colon tumors, SMYD3 is localized in the nucleus, where it interacts with RNA Pol II and H3K4me3 and functions as a selective transcriptional amplifier of oncogenes and genes that control cell proliferation and metastatic spread. Smyd3 expression has a high discriminative power for the characterization of liver tumors and positively correlates with poor prognosis. In lung and pancreatic cancer, SMYD3 acts in the cytoplasm, potentiating oncogenic Ras/ERK signaling through the methylation of the MAP3K2 kinase and the subsequent release from its inhibitor. A clinico-pathological analysis of lung cancer patients uncovers prognostic significance of SMYD3 only for first progression survival. However, stratification of patients according to their smoking history significantly expands the prognostic value of SMYD3 to overall survival and other features, suggesting that smoking-related effects saturate the clinical analysis and mask the function of SMYD3 as an oncogenic potentiator.

A Positive Regulatory Loop between a Wnt-Regulated Non-coding RNA and ASCL2 Controls Intestinal Stem Cell Fate.

The canonical Wnt pathway plays a central role in stem cell maintenance, differentiation, and proliferation in the intestinal epithelium. Constitutive, aberrant activity of the TCF4/ß-catenin transcriptional complex is the primary transforming factor in colorectal cancer. We identify a nuclear long non-coding RNA, termed WiNTRLINC1, as a direct target of TCF4/ß-catenin in colorectal cancer cells. WiNTRLINC1 positively regulates the expression of its genomic neighbor ASCL2, a transcription factor that controls intestinal stem cell fate. WiNTRLINC1 interacts with TCF4/ß-catenin to mediate the juxtaposition of its promoter with the regulatory regions of ASCL2. ASCL2, in turn, regulates WiNTRLINC1 transcriptionally, closing a feedforward regulatory loop that controls stem cell-related gene expression. This regulatory circuitry is highly amplified in colorectal cancer and correlates with increased metastatic potential and decreased patient survival. Our results uncover the interplay between non-coding RNA-mediated regulation and Wnt signaling and point to the diagnostic and therapeutic potential of WiNTRLINC1.

Smyd3 Is a Transcriptional Potentiator of Multiple Cancer-Promoting Genes and Required for Liver and Colon Cancer Development.

Smyd3 is a protein methyltransferase implicated in cancer development. Here we show that Smyd3 expression in mice is required for chemically induced liver and colon cancer formation. In these organs Smyd3 functions in the nucleus, stimulating the transcription of several key regulators involved in cell proliferation, epithelial-mesenchymal transition, the JAK/Stat3 oncogenic pathway, as well as the Myc and Ctnnb1 oncogenes. Smyd3 interacts with H3K4Me3-modified histone tails, which facilitates its recruitment to the core promoter regions of most active genes. Smyd3 binding density on target genes positively correlates with increased RNA polymerase-II density and transcriptional outputs. Despite its widespread distribution, the transcription-potentiating function of Smyd3 is restricted to a particular set of genes, whose expression is induced specifically during carcinogenesis.

Spontaneous development of hepatocellular carcinoma with cancer stem cell properties in PR-SET7-deficient livers.

PR-SET7-mediated histone 4 lysine 20 methylation has been implicated in mitotic condensation, DNA damage response and replication licensing. Here, we show that PR-SET7 function in the liver is pivotal for maintaining genome integrity. Hepatocyte-specific deletion of PR-SET7 in mouse embryos resulted in G2 phase arrest followed by massive cell death and defect in liver organogenesis. Inactivation at postnatal stages caused cell duplication-dependent hepatocyte necrosis, accompanied by inflammation, fibrosis and compensatory growth induction of neighboring hepatocytes and resident ductal progenitor cells. Prolonged necrotic regenerative cycles coupled with oncogenic STAT3 activation led to the spontaneous development of hepatic tumors composed of cells with cancer stem cell characteristics. These include a capacity to self-renew in culture or in xenografts and the ability to differentiate to phenotypically distinct hepatic cells. Hepatocellular carcinoma in PR-SET7-deficient mice displays a cancer stem cell gene signature specified by the co-expression of ductal progenitor markers and oncofetal genes.

Coactivation of GR and NFKB alters the repertoire of their binding sites and target genes.

Glucocorticoid receptor (GR) exerts anti-inflammatory action in part by antagonizing proinflammatory transcription factors such as the nuclear factor kappa-b (NFKB). Here, we assess the crosstalk of activated GR and RELA (p65, major NFKB component) by global identification of their binding sites and target genes. We show that coactivation of GR and p65 alters the repertoire of regulated genes and results in their association with novel sites in a mutually dependent manner. These novel sites predominantly cluster with p65 target genes that are antagonized by activated GR and vice versa. Our data show that coactivation of GR and NFKB alters signaling pathways that are regulated by each factor separately and provide insight into the networks underlying the GR and NFKB crosstalk.

A transcriptomic computational analysis of mastic oil-treated Lewis lung carcinomas reveals molecular mechanisms targeting tumor cell growth and survival.

BACKGROUND: Mastic oil from Pistacia lentiscus variation chia, a blend of bioactive terpenes with recognized medicinal properties, has been recently shown to exert anti-tumor growth activity through inhibition of cancer cell proliferation, survival, angiogenesis and inflammatory response. However, no studies have addressed its mechanisms of action at genome-wide gene expression level. METHODS: To investigate molecular mechanisms triggered by mastic oil, Lewis Lung Carcinoma cells were treated with mastic oil or DMSO and RNA was collected at five distinct time points (3-48 h). Microarray expression profiling was performed using Illumina mouse-6 v1 beadchips, followed by computational analysis. For a number of selected genes, RT-PCR validation was performed in LLC cells as well as in three human cancer cell lines of different origin (A549, HCT116, K562). PTEN specific inhibition by a bisperovanadium compound was applied to validate its contribution to mastic oil-mediated anti-tumor growth effects. RESULTS: In this work we demonstrated that exposure of Lewis lung carcinomas to mastic oil caused a time-dependent alteration in the expression of 925 genes. GO analysis associated expression profiles with several biological processes and functions. Among them, modifications on cell cycle/proliferation, survival and NF-kappaB cascade in conjunction with concomitant regulation of genes encoding for PTEN, E2F7, HMOX1 (up-regulation) and NOD1 (down-regulation) indicated some important mechanistic links underlying the anti-proliferative, pro-apoptotic and anti-inflammatory effects of mastic oil. The expression profiles of Hmox1, Pten and E2f7 genes were similarly altered by mastic oil in the majority of test cancer cell lines. Inhibition of PTEN partially reversed mastic oil effects on tumor cell growth, indicating a multi-target mechanism of action. Finally, k-means clustering, organized the significant gene list in eight clusters demonstrating a similar expression profile. Promoter analysis in a representative cluster revealed shared putative cis-elements suggesting a common regulatory transcription mechanism. CONCLUSIONS: Present results provide novel evidence on the molecular basis of tumor growth inhibition mediated by mastic oil and set a rational basis for application of genomics and bioinformatic methodologies in the screening of natural compounds with potential cancer chemopreventive activities.