Set Of Human Genes Research Articles

Enriched G4 forming repeats in the human genome are associated with robust well-coordinated transcription and reduced cancer transcriptome variation

Non-B DNA G-quadruplex (G4) structures with guanine (G) runs of 2-4 repeats can trigger opposing experimental transcriptional impacts. Here, we employed bioinformatic algorithms to comprehensively assess correlations of steady-state RNA transcript levels with all putative G4 sequence (pG4) locations genome-wide in three mammalian genomes and in normal and tumor human tissues. The human pG4-containing gene set displays higher expression levels than the set without pG4, supporting and extending some prior observations. pG4 enrichment at transcription start sites (TSS) in human, but not chimpanzee and mouse genomes, suggests possible positive selection pressure for pG4 at human TSS, potentially driving genome rewiring and gene expression divergence between human and chimpanzee. Comprehensive bioinformatic analyses revealed lower pG4-containing gene set variability in humans and among different pG4 genes in tumors. As G4 stabilizers are under therapeutic consideration for cancer and pathogens, such distinctions between human normal and tumor G4s along with other species merit attention. Furthermore, in germline and cancer sequences, the most mutagenic pG4 mapped to regions promoting alternative DNA structures. Overall findings establish high pG4 at TSS as a human genome attribute statistically associated with robust well-coordinated transcription and reduced cancer transcriptome variation with implications for biology, model organisms, and medicine.

Genome-Wide Transcriptional Roles of KSHV Viral Interferon Regulatory Factors in Oral Epithelial Cells.

The viral interferon regulatory factors (vIRFs) of KSHV are known to dysregulate cell signaling pathways to promote viral oncogenesis and to block antiviral immune responses to facilitate infection. However, it remains unknown to what extent each vIRF plays a role in gene regulation. To address this, we performed a comparative analysis of the protein structures and gene regulation of the four vIRFs. Our structure prediction analysis revealed that despite their low amino acid sequence similarity, vIRFs exhibit high structural homology in both their DNA-binding domain (DBD) and IRF association domain. However, despite this shared structural homology, we demonstrate that each vIRF regulates a distinct set of KSHV gene promoters and human genes in epithelial cells. We also found that the DBD of vIRF1 is essential in regulating the expression of its target genes. We propose that the structurally similar vIRFs evolved to possess specialized transcriptional functions to regulate specific genes.

Repositioning of ezetimibe for the treatment of idiopathic pulmonary fibrosis.

We previously identified ezetimibe, an inhibitor of Niemann-Pick C1-like intracellular cholesterol transporter 1 and European Medicines Agency-approved lipid-lowering agent, as a potent autophagy activator. However, its efficacy against pulmonary fibrosis has not yet been evaluated. This study aimed to determine whether ezetimibe has therapeutic potential against idiopathic pulmonary fibrosis. Primary lung fibroblasts isolated from both humans and mice were employed for mechanistic in vitro experiments. mRNA sequencing of human lung fibroblasts and gene set enrichment analysis were performed to explore the therapeutic mechanism of ezetimibe. A bleomycin-induced pulmonary fibrosis mouse model was used to examine in vivo efficacy of the drug. Tandem fluorescent-tagged microtubule-associated protein 1 light chain 3 transgenic mice were used to measure autophagic flux. Finally, the medical records of patients with idiopathic pulmonary fibrosis from three different hospitals were reviewed retrospectively, and analyses on survival and lung function were conducted to determine the benefits of ezetimibe. Ezetimibe inhibited myofibroblast differentiation by restoring the mechanistic target of rapamycin complex 1-autophagy axis with fine control of intracellular cholesterol distribution. Serum response factor, a potential autophagic substrate, was identified as a primary downstream effector in this process. Similarly, ezetimibe ameliorated bleomycin-induced pulmonary fibrosis in mice by inhibiting mechanistic target of rapamycin complex 1 activity and increasing autophagic flux, as observed in mouse lung samples. Patients with idiopathic pulmonary fibrosis who regularly used ezetimibe showed decreased rates of all-cause mortality and lung function decline. Our study presents ezetimibe as a potential novel therapeutic for idiopathic pulmonary fibrosis.

The T2T-CHM13 reference assembly uncovers essential WASH1 and GPRIN2 paralogues.

The recently published T2T-CHM13 reference assembly completed the annotation of the final 8% of the human genome. It introduced 1956 genes, close to 100 of which are predicted to be coding because they have a protein coding parent gene. Here, we confirm the coding status and functional relevance of two of these genes, paralogues of WASHC1 and GPRIN2. We find that LOC124908094, one of four novel subtelomeric WASH1 genes uncovered in the new assembly, produces the WASH1 protein that forms part of the vital actin-regulatory WASH complex. Its coding status is supported by abundant proteomics, conservation, and cDNA evidence. It was previously assumed that gene WASHC1 produced the functional WASH1 protein, but new evidence shows that WASHC1 is a human-derived duplication and likely to be one of 12 WASH1 pseudogenes in the human gene set. We also find that the T2T-CHM13 assembly has added a functionally important copy of GPRIN2 to the human gene set. We demonstrate that uniquely mapping peptides from proteomics databases support the novel LOC124900631 rather than the GRCh38 assembly GPRIN2 gene. These new additions to the set of human coding genes underlines the importance of the new T2T-CHM13 assembly. None.

Probable human origin of the SARS-CoV-2 polybasic furin cleavage motif

BackgroundThe key evolutionary step leading to the pandemic virus was the acquisition of the PRRA furin cleavage motif at the spike glycoprotein S1/S2 junction by a progenitor of SARS-CoV-2. Two of its features draw attention: (i) it is absent in other known lineage B beta-coronaviruses, including the newly discovered coronaviruses in bats from Laos and Vietnam, which are the closest known relatives of the covid virus; and, (ii) it introduced the pair of arginine codons (CGG-CGG), whose usage is extremely rare in coronaviruses. With an occurrence rate of only 3%, the arginine CGG codon is considered a minority in SARS CoV-2. On the other hand, Laos and Vietnam bat coronaviruses contain receptor-binding domains that are almost identical to that of SARS-CoV-2 and can therefore infect human cells despite the absence of the furin cleavage motif.ResultsBased on these data, the aim of this work is to provide a detailed sequence analysis between the SARS-CoV-2 S gene insert encoding PRRA and the human mRNA transcripts. The result showed a 100% match to several mRNA transcripts. The set of human genes whose mRNAs match this S gene insert are ubiquitous and highly expressed, e.g., the ATPase F1 (ATP5F1) and the ubiquitin specific peptidase 21 (USP21) genes; or specific genes of target organs or tissues of the SARS-CoV-2 infection (e.g., MEMO1, SALL3, TRIM17, CWC15, CCDC187, FAM71E2, GAB4, PRDM13). Results suggest that a recombination between the genome of a SARS-CoV-2 progenitor and human mRNA transcripts could be the origin of the S gene 12-nucleotide insert encoding the S protein PRRA motif.ConclusionsThe hypothesis of probable human origin of the SARS-CoV-2 polybasic furin cleavage motif is supported by: (i) the nature of human genes whose mRNA sequence 100% match the S gene insert; (ii) the synonymous base substitution in the arginine codons (CGG-CGG); and (iii) further spike glycoprotein PRRA-like insertions suggesting that the acquisition of PRRA may not have been a single recombination event.

Abstract 132: Proteomic Profiles Of Human Arterioles Isolated From Fresh Adipose Tissue Or Following Overnight Storage

Resistance arteries, or arterioles, are key determinants of the total peripheral vascular resistance, which, in turn, is a key determinant of arterial blood pressure (BP). However, the amount of protein available from one isolated human arteriole may be less than 5 μg, making proteomic analysis challenging. In addition, obtaining human arterioles requires manual dissection of unfrozen clinical specimens. This limits its feasibility, especially for powerful multi-center clinical studies in which clinical specimens need to be shipped overnight to a research lab for arteriole isolation. We performed a study to address low input, test overnight tissue storage, and develop a reference human arteriolar proteomic profile. We found that, in tandem mass tag proteomic analysis, the use of a booster channel consisting of endothelial and vascular smooth muscle cells (1:5 ratio) increased the number of proteins detected in a human arteriole segment with FDR < 0.01 from 1,051 to more than 3,000. We collected adipose tissues from three human subjects and isolated two arterioles from each fresh aliquot of the tissue or after 24h of cold storage of unfrozen aliquots in MACS tissue solution. The correlation coefficient of proteomic profile was similar (p=0.6) between replicate arterioles isolated freshly, following the cold storage, or before and after the cold storage. We built a human arteriolar proteomic profile consisting of 3,836 proteins based on the analysis of 12 arteriole samples from the three subjects. The average arteriolar protein had a mean copy number of 2.76х10 6 per cell using the histone proteomic ruler, which is based on the fact that the mass of DNA per cell is approximately equal to the protein mass of histones. Transgelin was the most abundant protein detected. We curated a set of BP-relevant human genes, which encode 1,945 proteins. Of these BP-relevant proteins, 476 (12.5%) were detected in the arteriolar proteome, which was a significant overrepresentation (p<0.05, Chi squared test). These findings demonstrate that proteomic analysis is feasible with arterioles isolated from human adipose tissue following cold overnight storage and provide a reference human arteriolar proteome profile highly valuable for studies of arteriole-related traits.

Semaglutide Has Beneficial Effects on Non-Alcoholic Steatohepatitis in Ldlr-/-.Leiden Mice.

Semaglutide, a glucagon-like peptide-1 receptor agonist, is an antidiabetic medication that has recently been approved for the treatment of obesity as well. Semaglutide is postulated to be a promising candidate for the treatment of non-alcoholic steatohepatitis (NASH). Here, Ldlr-/-.Leiden mice received a fast-food diet (FFD) for 25 weeks, followed by another 12 weeks on FFD with daily subcutaneous injections of semaglutide or vehicle (control). Plasma parameters were evaluated, livers and hearts were examined, and hepatic transcriptome analysis was performed. In the liver, semaglutide significantly reduced macrovesicular steatosis (-74%, p < 0.001) and inflammation (-73%, p < 0.001) and completely abolished microvesicular steatosis (-100%, p < 0.001). Histological and biochemical assessment of hepatic fibrosis showed no significant effects of semaglutide. However, digital pathology revealed significant improvements in the degree of collagen fiber reticulation (-12%, p < 0.001). Semaglutide did not affect atherosclerosis relative to controls. Additionally, we compared the transcriptome profile of FFD-fed Ldlr-/-.Leiden mice with a human gene set that differentiates human NASH patients with severe fibrosis from those with mild fibrosis. In FFD-fed Ldlr-/-.Leiden control mice, this gene set was upregulated as well, while semaglutide predominantly reversed this gene expression. Using a translational model with advanced NASH, we demonstrated that semaglutide is a promising candidate with particular potential for the treatment of hepatic steatosis and inflammation, while for the reversal of advanced fibrosis, combinations with other NASH agents may be necessary.

Human gene age dating reveals an early and rapid evolutionary construction of the adaptive immune system.

T cells are a type of white blood cell that play a critical role in the immune response against foreign pathogens through a process called T Cell Adaptive Immunity (TCAI). However, the evolution of the genes and nucleotide sequences involved in TCAI is not well understood. To investigate this, we performed comparative studies of gene annotations and genome assemblies of 28 vertebrate species and identified sets of human genes that are involved in TCAI, carcinogenesis, and ageing. We found that these gene sets share interaction pathways which may have contributed to the evolution of longevity in the vertebrate lineage leading to humans. Our human gene age dating analyses revealed that there was rapid origination of genes with TCAI-related functions prior to the Cretaceous eutherian radiation and these new genes mainly encode negative regulators. We identified no new TCAI-related genes after the divergence of placental mammals, but we did detect an extensive number of amino acid substitutions under strong positive selection in recently evolved human immunity genes suggesting they are co-evolving with adaptive immunity. More specifically, we observed that antigen processing and presentation and checkpoint genes are significantly enriched among new genes evolving under positive selection. These observations reveal an evolutionary process of T Cell Adaptive Immunity that were associated with rapid gene duplication in the early stages of vertebrates and subsequent sequence changes in TCAI-related genes. These processes together suggest an early genetic construction of the vertebrate immune system and subsequent molecular adaptation to diverse antigens.

Abstract 708: Murine Atherosclerosis Data Integration: The Sum Is Greater Than The Parts

Translation of murine atherosclerosis studies to human disease has limited reproducibility. Our prior work showed that aggregated data from one mouse atherosclerosis model identified etiologic pathways and networks associated with human atherosclerotic cardiovascular disease (ASCVD). Therefore, we extracted data from 991 ApoE knockout (KO) manuscripts and 544 LDLR KO manuscripts to identify 405 unique genes associated with significant changes in plaque size in mice. The impact of each gene on plaque size was analyzed using Ingenuity Pathway Analysis to contextualize the gene data. The most significantly activated pathways for genes that increased plaque size were neuroinflammation (p=3.02 E-55) and hepatic fibrosis signaling (p=2.25 E-52), while the LXR/RXR pathway was the most significantly downregulated (p=1.32 E-37). We then performed association analyses from Vanderbilt University Medical Center’s BioVU using genetically predicted gene expression of the individual genes from these datasets. Neither the aggregated genes identified in each individual mouse model nor the genes directly studied in the murine models that mapped to the top pathways were significantly enriched for association with ASCVD in humans. However, when homologous human gene sets derived from the top identified murine pathways were considered as a whole, they were highly significantly enriched for association with ASCVD (Fig 1). Further, this method identified human genes significantly associated with ASCVD that were not directly studied in mice. For example, in the neuroinflammation pathway, the genes studied in mice were not enriched for ASCVD, yet 61/117 genes (52%) in the global pathway were associated with ASCVD (p=0.00017), 43 of which were not directly studied in mice. Thus, integration of murine data identifies novel pathways significantly associated with human ASCVD thereby providing mechanistic insight into human disease only possible through contextualization.

Abstract 6568: Extended support for model organisms: The next frontier for the molecular signatures database

Abstract In 2005, we introduced the Gene Set Enrichment Analysis (GSEA) method to enable the identification and estimation of significance of activated biological pathways and processes in molecular data. Serving a community of over 300,000 registered users, and contributing scientific analysis to more than 30,000 publications, GSEA has become ubiquitous in gene expression analysis and particularly in cancer research. The power of GSEA relies, in a large part, on the Molecular Signatures Database (MSigDB) that provides tens of thousands of expertly curated and annotated gene sets representing cellular processes, canonical pathways, response signatures, etc., derived from prior studies, public data, and curated pathway databases. Here we describe an important recent addition to the GSEA/MSigDB resource. Historically, the GSEA-MSigDB resource has focused specifically on human data, offering guidance and support for analysis of model organism data through limited provision of orthology mapping files. In recognition of the importance of the mouse as a model organism for cancer research, we recently significantly expanded our support for mouse data in two ways. (1) We introduced Mouse MSigDB with a total of about 15,000 sets that are provided in the native mouse gene symbol space and are derived from published mouse datasets and mouse-focused resources. This allows GSEA analysis of mouse datasets without the need for orthology conversion. (2) For investigators who wish to use MSigDB’s human gene sets to analyze mouse data, or vice versa, we also implemented a new procedure for orthology mapping gene identifiers in an input dataset to match the identifiers of the gene sets. The new mapping considers the particular requirements of gene set enrichment analysis, and incorporates high confidence gene ortholog data from the Alliance of Genome Resources and ortholog-match refinement utilizing Ensembl’s ortholog datasets. Together, these new additions offer a substantial step forward for making mouse a first-class citizen of the GSEA-MSigDB ecosystem. The new Mouse MSigDB, and the full GSEA-MSigDB resource, is available at gsea-msigdb.org. Citation Format: Anthony S. Castanza, Jill M. Recla, David Eby, Alexander T. Wenzel, Helga Thorvaldsdottir, Carol J. Bult, Jill P. Mesirov. Extended support for model organisms: The next frontier for the molecular signatures database [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 6568.

Reprogramming the Circadian Dynamics of Epileptic Genes in Mouse Temporal Lobe Epilepsy

Temporal lobe epilepsy (TLE) is a common and severe epilepsy displaying rhythmicity in humans and animals. However, how the circadian clock contributes to TLE remains elusive. A recent circadian analysis of the ventral hippocampal transcriptome of pilocarpine-induced TLE mice revealed as many as 1650 rhythmically expressed transcripts. Here, a comparison of the mouse ventral hippocampal transcriptome with the human epilepsy-related gene set identified 315 possible mouse epilepsy-related genes. Rhythmicity analysis classified them into arrhythmicity, loss-of-rhythmicity, gain-of-rhythmicity, and rhythmicity-maintaining groups. KEGG and GO analyses of these mouse epilepsy genes suggest their involvement in circadian entrainment. In TLE mice, Htr1d, Drd2, and Chrna3 lose rhythmicity, but P2rx7 gains rhythmicity; the up-regulation of Htr1d and Drd2 and down-regulation of Chrna3 inhibit adenylate cyclase (AC), and up-regulation of Htr1d, Drd2, and P2rx7 activates protein kinase C (PKC). Together, these results suggest that epilepsy can disrupt the circadian dynamics of the epileptic genes, shed light on possible TLE pathogenesis, and provide potential targets for TLE diagnosis and chronotherapy.

RNA sequencing least shrew (Cryptotis parva) brainstem and gut transcripts following administration of a selective substance P neurokinin NK1 receptor agonist and antagonist expands genomics resources for emesis research.

The least shrew is among the subset of animals that are capable of vomiting and therefore serves as a valuable research model for investigating the biochemistry, molecular biology, pharmacology, and genomics of emesis. Both nausea and vomiting are associated with a variety of illnesses (bacterial/viral infections, bulimia, exposure to toxins, gall bladder disease), conditions (pregnancy, motion sickness, emotional stress, overeating) and reactions to drugs (chemotherapeutics, opiates). The severe discomfort and intense fear associated with the stressful symptoms of nausea and emesis are the major reason for patient non-compliance when being treated with cancer chemotherapeutics. Increased understanding of the physiology, pharmacology and pathophysiology underlying vomiting and nausea can accelerate progress for developing new antiemetics. As a major animal model for emesis, expanding genomic knowledge associated with emesis in the least shrew will further enhance the laboratory utility of this model. A key question is which genes mediate emesis, and are they expressed in response to emetics/antiemetics. To elucidate the mediators of emesis, in particular emetic receptors, their downstream signaling pathways, as well as the shared emetic signals, we carried out an RNA sequencing study focused on the central and peripheral emetic loci, the brainstem and gut. Thus, we sequenced RNA extracted from brainstem and gut tissues from different groups of least shrews treated with either a neurokinin NK1 receptor selective emetic agonist, GR73632 (5mg/kg, i.p.), its corresponding selective antagonist netupitant (5mg/kg, i.p.), a combination of these two agents, versus their corresponding vehicle-pretreated controls and drug naïve animals. The resulting sequences were processed using a de novo transcriptome assembly and used it to identify orthologs within human, dog, mouse, and ferret gene sets. We compared the least shrew to human and a veterinary species (dog) that may be treated with vomit-inducing chemotherapeutics, and the ferret, another well-established model organism for emesis research. The mouse was included because it does not vomit. In total, we identified a final set of 16,720 least shrew orthologs. We employed comparative genomics analyses as well as gene ontology enrichment, KEGG pathway enrichment and phenotype enrichment to better understand the molecular biology of genes implicated in vomiting.

Identification of TAC1 Associated with Alzheimer's Disease Using a Robust Rank Aggregation Approach.

Alzheimer's disease (AD) brings heavy burden to society and family. There is an urgent need to find effective methods for disease diagnosis and treatment. The robust rank aggregation (RRA) approach that could aggregate the resulting gene lists has been widely utilized in genomic data analysis. To identify hub genes using RRA approach in AD. Seven microarray datasets in frontal cortex from GEO database were used to identify differential expressed genes (DEGs) in AD patients using RRA approach. STRING was performed to explore the protein-to-protein interaction (PPI). Gene Ontology enrichment and Kyoto Encyclopedia of Genes and Genomes pathway analyses were utilized for enrichment analysis. Human Gene Connectome and Gene Set Enrichment Analysis were used for functional annotation. Finally, the expression levels of hub genes were validated in the cortex of 5xFAD mice by quantitative real-time polymerase chain reaction. After RRA analysis, 473 DEGs (216 upregulated and 257 downregulated) were identified in AD samples. PPI showed that DEGs had a total of 416 nodes and 2750 edges. These genes were divided into 17 clusters, each of which contains at least three genes. After functional annotation and enrichment analysis, TAC1 is identified as the hub gene and may be related to synaptic function and inflammation. In addition, Tac1 was found downregulated in cortices of 5xFAD mice. In the current study, TAC1 is identified as a key gene in the frontal cortex of AD, providing insight into the possible pathogenesis and potential therapeutic targets for this disease.

Exploiting Multi-Omics Profiling and Systems Biology to Investigate Functions of TOMM34.

Although modern biology is now in the post-genomic era with vastly increased access to high-quality data, the set of human genes with a known function remains far from complete. This is especially true for hundreds of mitochondria-associated genes, which are under-characterized and lack clear functional annotation. However, with the advent of multi-omics profiling methods coupled with systems biology algorithms, the cellular role of many such genes can be elucidated. Here, we report genes and pathways associated with TOMM34, Translocase of Outer Mitochondrial Membrane, which plays role in the mitochondrial protein import as a part of cytosolic complex together with Hsp70/Hsp90 and is upregulated in various cancers. We identified genes, proteins, and metabolites altered in TOMM34-/- HepG2 cells. To our knowledge, this is the first attempt to study the functional capacity of TOMM34 using a multi-omics strategy. We demonstrate that TOMM34 affects various processes including oxidative phosphorylation, citric acid cycle, metabolism of purine, and several amino acids. Besides the analysis of already known pathways, we utilized de novo network enrichment algorithm to extract novel perturbed subnetworks, thus obtaining evidence that TOMM34 potentially plays role in several other cellular processes, including NOTCH-, MAPK-, and STAT3-signaling. Collectively, our findings provide new insights into TOMM34's cellular functions.

Comprehensive Molecular Analyses of Notch Pathway-Related Genes to Predict Prognosis and Immunotherapy Response in Patients with Gastric Cancer.

Gastric cancer (GC) is a highly molecular heterogeneous tumor with unfavorable outcomes. The Notch signaling pathway is an important regulator of immune cell differentiation and has been associated with autoimmune disorders, the development of tumors, and immunomodulation caused by tumors. In this study, by developing a gene signature based on genes relevant to the Notch pathway, we could improve our ability to predict the outcome of patients with GC. From the TCGA database, RNA sequencing data of GC tumors and associated normal tissues were obtained. Microarray data were collected from GEO datasets. The Molecular Signature Database (MSigDB) was accessed in order to retrieve sets of human Notch pathway-related genes (NPRGs). The LASSO analysis performed on the TCGA cohort was used to generate a multigene signature based on prognostic NPRGs. In order to validate the gene signature, the GEO cohort was utilized. Using the CIBERSORT method, we were able to determine the amounts of immune cell infiltration in the GC. In this study, a total of 21 differentially expressed NPRGs were obtained between GC specimens and nontumor specimens. The construction of a prognostic prediction model for patients with GC involved the identification and selection of three different NPRGs. According to the appropriate cutoff value, the patients with GC were divided into two groups: those with a low risk and those with a high risk. The time-dependent ROC curves demonstrated that the new model had satisfactory performance when it came to prognostic prediction. Multivariate assays confirmed that the risk score was an independent marker that may be used to predict the outcome of GC. In addition, the generated nomogram demonstrated a high level of predictive usefulness. Moreover, the scores of immunological infiltration of the majority of immune cells were distinctly different between the two groups, and the low-risk group responded to immunotherapy in a significantly greater degree. According to the results of a functional enrichment study of candidate genes, there are multiple pathways and processes associated with cancer. Taken together, a new gene model associated with the Notch pathway may be utilized for the purpose of predicting the prognosis of GC. One potential method of treatment for GC is to focus on NPRGs.

A pan-cancer analysis of the oncogenic role of ribonucleotide reductase subunit M2 in human tumors.

Recent studies have identified ribonucleotide reductase subunit M2 (RRM2) as a putative promoter of tumors. However, no systematic analysis of its carcinogenicity has been conducted. The potential functions of RRM2 in various tumor types were investigated using data from the Genotype-Tissue Expression (GTEx), the Clinical Proteomic Tumor Analysis Consortium (CPTAC), the Cancer Genome Atlas (TCGA), the Human Protein Atlas (HPA), cBioPortal, GEPIA, String, and Gene Set Enrichment Analysis (GSEA). We analyzed the difference in mRNA and protein expression, pathological stage, survival, mutation, tumor microenvironment (TME), and immune cell infiltration in relation to RRM2. Meanwhile, using TCGA and the Tumor Immune Estimation Resource 2 (TIMER 2), the associations between RRM2 expression, immune infiltration, and immune-related genes were assessed. Additionally, CCK-8, Edu and RT-PCR assays were used to validate that RRM2 acts as an oncogene in liver cancer cells and its association with HBx. A cohort of liver hepatocellular carcinoma (LIHC) patients (n=154) from Huashan Hospital was analyzed for the expression of RRM2 and the association between RRM2 and immune infiltration. Using the GTEx and TCGA databases, we discovered that 28 tumors expressed RRM2 at significantly higher levels than the corresponding normal tissues. Increased RRM2 expression may be predictive of a poor overall survival (OS) in patients with seven different cancers. GO, KEGG, and GSEA analyses revealed that the biological process of RRM2 was associated with the regulation of carcinogenic processes and immune pathways in a variety of tumor types. The expression of RRM2 was highly correlated with maker genes involved in immune activation and immunosuppression, immune checkpoints, DNA mismatch repair system (MMR), and the infiltration levels of Tregs and macrophages (TAMs), suggesting that the carcinogenic effect of RRM2 may be achieved by regulating immune related genes. Moreover, as demonstrated by CCK-8 and Edu assays, RRM2 was an oncogene in liver cancer cells. We confirmed for the first time that RRM2 was significantly upregulated by HBx, suggesting that RRM2 may be a key regulator of LIHC induced by HBV. IHC analysis validated the upregulated expression of RRM2 protein and its correlation with immune infiltration makers in a LIHC patient cohort. RRM2 may be a valuable molecular biomarker for predicting prognosis and immunotherapeutic efficacy in pan-cancer, particularly in LIHC.

Origins and Evolution of Human Tandem Duplicated Exon Substitution Events.

The mutually exclusive splicing of tandem duplicated exons produces protein isoforms that are identical save for a homologous region that allows for the fine tuning of protein function. Tandem duplicated exon substitution events are rare, yet highly important alternative splicing events. Most events are ancient, their isoforms are highly expressed, and they have significantly more pathogenic mutations than other splice events. Here, we analyzed the physicochemical properties and functional roles of the homologous polypeptide regions produced by the 236 tandem duplicated exon substitutions annotated in the human gene set. We find that the most important structural and functional residues in these homologous regions are maintained, and that most changes are conservative rather than drastic. Three quarters of the isoforms produced from tandem duplicated exon substitution events are tissue-specific, particularly in nervous and cardiac tissues, and tandem duplicated exon substitution events are enriched in functional terms related to structures in the brain and skeletal muscle. We find considerable evidence for the convergent evolution of tandem duplicated exon substitution events in vertebrates, arthropods, and nematodes. Twelve human gene families have orthologues with tandem duplicated exon substitution events in both Drosophila melanogaster and Caenorhabditis elegans. Six of these gene families are ion transporters, suggesting that tandem exon duplication in genes that control the flow of ions into the cell has an adaptive benefit. The ancient origins, the strong indications of tissue-specific functions, and the evidence of convergent evolution suggest that these events may have played important roles in the evolution of animal tissues and organs.

Role of liquid fructose/sucrose in regulating the hepatic transcriptome in a high-fat Western diet model of NAFLD

Nonalcoholic fatty liver disease (NAFLD), which ranges from simple steatosis to nonalcoholic steatohepatitis (NASH), is the most common chronic liver disease. Yet, the molecular mechanisms for the progression of steatosis to NASH remain largely undiscovered. Thus, there is a need for identifying specific gene and pathway changes that drive the progression of NAFLD. This study uses high-fat Western diet (HFWD) together with liquid sugar [fructose and sucrose (F/S)] feeding for 12 weeks in mice to induce obesity and examine hepatic transcriptomic changes that occur in NAFLD progression. The combination of a HFWD+F/S in the drinking water exacerbated HFWD-induced obesity, hyperinsulinemia, hyperglycemia, hepatic steatosis, inflammation, and human and murine fibrosis gene set enrichment that is consistent with progression to NASH. RNAseq analysis revealed differentially expressed genes (DEGs) associated with HFWD and HFWD+F/S dietary treatments compared to Chow-fed mice. However, liquid sugar consumption resulted in a unique set of hepatic DEGs in HFWD+F/S-fed mice, which were enriched in the complement and coagulation cascades using network and biological analysis. Cluster analysis identified Orosomucoid (ORM) as a HFWD+F/S upregulated complement and coagulation cascades gene that was also upregulated in hepatocytes treated with TNFα or free fatty acids in combination with hypoxia. ORM expression was found to correlate with NAFLD parameters in obese mice. Taken together, this study examined key genes, biological processes, and pathway changes in the liver of HFWD+F/S mice in an effort to provide insight into the molecular basis for which the addition of liquid sugar promotes the progression of NAFLD.

Toxicity by descent: A comparative approach for chemical hazard assessment.

Toxicology is traditionally divided between human and eco-toxicology. In the shared pursuit of environmental health, this separation does not account for discoveries made in the comparative studies of animal genomes. Here, we provide evidence on the feasibility of understanding the health impact of chemicals on all animals, including ecological keystone species and humans, based on a significant number of conserved genes and their functional associations to health-related outcomes across much of animal diversity. We test four conditions to understand the value of comparative genomics data to inform mechanism-based human and environmental hazard assessment: (1) genes that are most fundamental for health evolved early during animal evolution; (2) the molecular functions of pathways are better conserved among distantly related species than the individual genes that are members of these pathways; (3) the most conserved pathways among animals are those that cause adverse health outcomes when disrupted; (4) gene sets that serve as molecular signatures of biological processes or disease-states are largely enriched by evolutionarily conserved genes across the animal phylogeny. The concept of homology is applied in a comparative analysis of gene families and pathways among invertebrate and vertebrate species compared with humans. Results show that over 70% of gene families associated with disease are shared among the greatest variety of animal species through evolution. Pathway conservation between invertebrates and humans is based on the degree of conservation within vertebrates and the number of interacting genes within the human network. Human gene sets that already serve as biomarkers are enriched by evolutionarily conserved genes across the animal phylogeny. By implementing a comparative method for chemical hazard assessment, human and eco-toxicology converge towards a more holistic and mechanistic understanding of toxicity disrupting biological processes that are important for health and shared among animals (including humans).

OS10.6.A What is the initial cell in the subventricular zone for human glioblastoma genesis?

Abstract Background We all have a fundamental question about why glioblastoma (GBM) develops. In order to find the answer to this fundamental question, if you find out what the first cell is, you will get closer to the answer. GBM arises from the subventricular zone (SVZ). GBM is one of the most devastating tumour of human brain as the most optimal treatment barely prolongs the survival, and it does not cure the disease. As the majority of GBM tissues show copy number variations (CNV) of co-altered chromosomal 7 gain and 10 loss, we hypothesized the origin cell (Oc) of SVZ may be traced back with these markers. The cellular identity of the Oc is still unknown and it is different from the tumour-derived progenitor-like cells. We aimed to define these cells from the SVZ that have a potential to get activated into GBM. Material and Methods We compared bulk RNA sequencing (RNAseq) data of IDHwt GBM tumor tissue (n=122), tumor free SVZ from GBM patients (n=40), tumor-free control SVZ of non-glial tumor (n=9). Pared single nucleus RNAseq (snRNAseq) or single cell RNAseq (scRNAseq) samples of tumor free SVZ (n=11) and GBM tumor (n=8), were done to see cell specific CNVs. We developed genetically engineered mouse models for GBM genesis introducing three driver mutations (TP53, PTEN, and EGFRviii) into SVZ to isolate mouse Oc (mOc) and mouse cancer cells (mCc). The biological characteristics of separated mOc and mCc were compared. Bulk RNAseq and scRNAseq were performed on these cells (mOc, mCc), and their cellular state was compared with the human gene set. Results In this work, we found two types of the Oc in the RNA sequencing of 60 human tumour free-SVZ samples. Furthermore, single-cell level analysis revealed that two Oc types in SVZ harbor ongoing patterns of CNV co-alterations from Oc1 to Oc2, and finally to GBM. The Oc1 type cells contained the CNV signature of Oc2 ancestor with neural progenitor cell (NPC) signature. Oc2 type cells expressed a high level of EGFR than other cells with astrocyte-like cell signature. Both of these cells expressed oligodendrocyte progenitor cell (OPC)-like signatures in the SVZ. We validated the human-based findings by using the P53/PTEN/EGFR-mutant mouse model with EGFR/tdTomato overexpression and P53/PTEN knockout in the SVZ cells. As a result, non-tumourigenic and highly motile Oc-like cell-states are found in the mouse models, supporting the firework-like migration pattern from the SVZ. Conclusion Our results demonstrate how members of Oc preoccupy the SVZ, known as the stem cell niche and give rise to the tumour. We anticipate that a new therapy may emerge by targeting the Oc in the SVZ.