GENI: A web server to identify gene set enrichments in tumor samples.

The Cancer Genome Atlas (TCGA) and analogous projects have yielded invaluable tumor-associated genomic data. Despite several web-based platforms designed to enhance accessibility, certain analyses require prior bioinformatic expertise. To address this need, we developed Gene ENrichment Identifier (GENI, https://www.shaullab.com/geni), which is designed to promptly compute correlations for genes of interest against the entire transcriptome and rank them against well-established biological gene sets. Additionally, it generates comprehensive tables containing genes of interest and their corresponding correlation coefficients, presented in publication-quality graphs. Furthermore, GENI has the capability to analyze multiple genes simultaneously within a given gene set, elucidating their significance within a specific biological context. Overall, GENI's user-friendly interface simplifies the biological interpretation and analysis of cancer patient-associated data, advancing the understanding of cancer biology and accelerating scientific discoveries.

In vivo and in vitro genome editing to explore GNE functions.

GNE myopathy is an adult onset neuromuscular disorder characterized by slowly progressive distal and proximal muscle weakness, caused by missense recessive mutations in the GNE gene. Although the encoded bifunctional enzyme is well known as the limiting factor in the biosynthesis of sialic acid, no clear mechanisms have been recognized to account for the muscle atrophic pathology, and novel functions for GNE have been hypothesized. Two major issues impair studies on this protein. First, the expression of the GNE protein is minimal in human and mice muscles and there is no reliable antibody to follow up endogenous expression. Second, no reliable animal model is available for the disease and cellular models from GNE myopathy patients' muscle cells (expressing the mutated protein) are less informative than expected. In order to broaden our knowledge on GNE functions in muscle, we have taken advantage of the CRISPR/Cas9 method for genome editing to first, add a tag to the endogenous Gne gene in mouse, allowing the determination of the spatiotemporal expression of the protein in the organism, using well established and reliable antibodies against the specific tag. In addition we have generated a Gne knock out murine muscle cell lineage to identify the events resulting from the total lack of the protein. A thorough multi-omics analysis of both cellular systems including transcriptomics, proteomics, phosphoproteomics and ubiquitination, unraveled novel pathways for Gne, in particular its involvement in cell cycle control and in the DNA damage/repair pathways. The elucidation of fundamental mechanisms of Gne in normal muscle may contribute to the identification of the disrupted functions in GNE myopathy, thus, to the definition of novel biomarkers and possible therapeutic targets for this disease.

Generation and characterization of a novel gne Knockout Model in Zebrafish.

GNE Myopathy is a rare, recessively inherited neuromuscular worldwide disorder, caused by a spectrum of bi-allelic mutations in the human GNE gene. GNE encodes a bi-functional enzyme responsible for the rate-limiting step of sialic acid biosynthesis pathway. However, the process in which GNE mutations lead to the development of a muscle pathology is not clear yet. Cellular and mouse models for GNE Myopathy established to date have not been informative. Further, additional GNE functions in muscle have been hypothesized. In these studies, we aimed to investigate gne functions using zebrafish genetic and transgenic models, and characterized them using macroscopic, microscopic, and molecular approaches. We first established transgenic zebrafish lineages expressing the human GNE cDNA carrying the M743T mutation, driven by the zebrafish gne promoter. These fish developed entirely normally. Then, we generated a gne knocked-out (KO) fish using the CRISPR/Cas9 methodology. These fish died 8-10 days post-fertilization (dpf), but a phenotype appeared less than 24 h before death and included progressive body axis curving, deflation of the swim bladder and decreasing movement and heart rate. However, muscle histology uncovered severe defects, already at 5 dpf, with compromised fiber organization. Sialic acid supplementation did not rescue the larvae from this phenotype nor prolonged their lifespan. To have deeper insights into the potential functions of gne in zebrafish, RNA sequencing was performed at 3 time points (3, 5, and 7 dpf). Genotype clustering was progressive, with only 5 genes differentially expressed in gne KO compared to gne WT siblings at 3 dpf. Enrichment analyses of the primary processes affected by the lack of gne also at 5 and 7 dpf point to the involvement of cell cycle and DNA damage/repair processes in the gne KO zebrafish. Thus, we have established a gne KO zebrafish lineage and obtained new insights into gne functions. This is the only model where GNE can be related to clear muscle defects, thus the only animal model relevant to GNE Myopathy to date. Further elucidation of gne precise mechanism-of-action in these processes could be relevant to GNE Myopathy and allow the identification of novel therapeutic targets.

Upregulation of COX4-2 via HIF-1α in Mitochondrial COX4-1 Deficiency.

Cytochrome-c-oxidase (COX) subunit 4 (COX4) plays important roles in the function, assembly and regulation of COX (mitochondrial respiratory complex 4), the terminal electron acceptor of the oxidative phosphorylation (OXPHOS) system. The principal COX4 isoform, COX4-1, is expressed in all tissues, whereas COX4-2 is mainly expressed in the lungs, or under hypoxia and other stress conditions. We have previously described a patient with a COX4-1 defect with a relatively mild presentation compared to other primary COX deficiencies, and hypothesized that this could be the result of a compensatory upregulation of COX4-2. To this end, COX4-1 was downregulated by shRNAs in human foreskin fibroblasts (HFF) and compared to the patient's cells. COX4-1, COX4-2 and HIF-1α were detected by immunocytochemistry. The mRNA transcripts of both COX4 isoforms and HIF-1 target genes were quantified by RT-qPCR. COX activity and OXPHOS function were measured by enzymatic and oxygen consumption assays, respectively. Pathways were analyzed by CEL-Seq2 and by RT-qPCR. We demonstrated elevated COX4-2 levels in the COX4-1-deficient cells, with a concomitant HIF-1α stabilization, nuclear localization and upregulation of the hypoxia and glycolysis pathways. We suggest that COX4-2 and HIF-1α are upregulated also in normoxia as a compensatory mechanism in COX4-1 deficiency.

Slow Transcription of the 99a/let-7c/125b-2 Cluster Results in Differential MiRNA Expression and Promotes Melanoma Phenotypic Plasticity.

Almost half of the human microRNAs (miRNAs) are encoded in clusters. Although transcribed as a single unit, the levels of individual mature miRNAs often differ. The mechanisms underlying differential biogenesis of clustered miRNAs and the resulting physiological implications are mostly unknown. In this study, we report that the melanoma master transcription regulator MITF regulates the differential expression of the 99a/let-7c/125b-2 cluster by altering the distribution of RNA polymerase II along the cluster. We discovered that MITF interacts with TRIM28, a known inhibitor of RNA polymerase II transcription elongation, at the mIR-let-7c region, resulting in the pausing of RNA polymerase II activity and causing an elevation in mIR-let-7c expression; low levels of RNA polymerase II occupation over miR-99a and miR-125b-2 regions decreases their biogenesis. Furthermore, we showed that this differential expression affects the phenotypic state of melanoma cells. RNA-sequencing analysis of proliferative melanoma cells that express miR-99a and miR-125b mimics revealed a transcriptomic shift toward an invasive phenotype. Conversely, expression of a mIR-let-7c mimic in invasive melanoma cells induced a shift to a more proliferative state. We confirmed direct target genes of these miRNAs, including FGFR3, BAP1, Bcl2, TGFBR1, and CDKN1A. Our study demonstrates an MITF-governed biogenesis mechanism that results in differential expression of clustered 99a/let-7c/125b-2 miRNAs that control melanoma progression.

Chronic expression of p16INK4a in the epidermis induces Wnt-mediated hyperplasia and promotes tumor initiation.

p16INK4a (CDKN2A) is a central tumor suppressor, which induces cell-cycle arrest and senescence. Cells expressing p16INK4a accumulate in aging tissues and appear in premalignant lesions, yet their physiologic effects are poorly understood. We found that prolonged expression of transgenic p16INK4a in the mouse epidermis induces hyperplasia and dysplasia, involving high proliferation rates of keratinocytes not expressing the transgene. Continuous p16INK4a expression increases the number of epidermal papillomas formed after carcinogen treatment. Wnt-pathway ligands and targets are activated upon prolonged p16INK4a expression, and Wnt inhibition suppresses p16INK4a-induced hyperplasia. Senolytic treatment reduces p16INK4a-expressing cell numbers, and inhibits Wnt activation and hyperplasia. In human actinic keratosis, a precursor of squamous cell carcinoma, p16INK4a-expressing cells are found adjacent to dividing cells, consistent with paracrine interaction. These findings reveal that chronic p16INK4a expression is sufficient to induce hyperplasia through Wnt-mediated paracrine stimulation, and suggest that this tumor suppressor can promote early premalignant epidermal lesion formation.