GeneMANIA prediction server 2013 update

Exploring gene regulatory interaction networks and predicting therapeutic molecules for hypopharyngeal cancer and EGFR-mutated lung adenocarcinoma

Hypopharyngeal cancer gene regulatory networks and therapeutic molecules are a disease that is associated with EGFR-mutated lung adenocarcinoma. Here we utilized a bioinformatics approach to identify genetic commonalities between these two diseases. To this end, we examined microarray datasets from GEO (Gene Expression Omnibus) to identify differentially expressed genes, common genes, and hub genes between the selected two diseases. Our analyses identified potential therapeutic molecules for the selected diseases based on 10 hub genes with the highest interactions according to the degree topology method and the maximum clique centrality (MCC). These therapeutic molecules may have the potential for simultaneous treatment of these diseases.

In silico functional, structural and pathogenicity analysis of missense single nucleotide polymorphisms in human MCM6 gene

Single nucleotide polymorphisms (SNPs) are one of the most common determinants and potential biomarkers of human disease pathogenesis. SNPs could alter amino acid residues, leading to the loss of structural and functional integrity of the encoded protein. In humans, members of the minichromosome maintenance (MCM) family play a vital role in cell proliferation and have a significant impact on tumorigenesis. Among the MCM members, the molecular mechanism of how missense SNPs of minichromosome maintenance complex component 6 (MCM6) contribute to DNA replication and tumor pathogenesis is underexplored and needs to be elucidated. Hence, a series of sequence and structure-based computational tools were utilized to determine how mutations affect the corresponding MCM6 protein. From the dbSNP database, among 15,009 SNPs in the MCM6 gene, 642 missense SNPs (4.28%), 291 synonymous SNPs (1.94%), and 12,500 intron SNPs (83.28%) were observed. Out of the 642 missense SNPs, 33 were found to be deleterious during the SIFT analysis. Among these, 11 missense SNPs (I123S, R207C, R222C, L449F, V456M, D463G, H556Y, R602H, R633W, R658C, and P815T) were found as deleterious, probably damaging, affective and disease-associated. Then, I123S, R207C, R222C, V456M, D463G, R602H, R633W, and R658C missense SNPs were found to be highly harmful. Six missense SNPs (I123S, R207C, V456M, D463G, R602H, and R633W) had the potential to destabilize the corresponding protein as predicted by DynaMut2. Interestingly, five high-risk mutations (I123S, V456M, D463G, R602H, and R633W) were distributed in two domains (PF00493 and PF14551). During molecular dynamics simulations analysis, consistent fluctuation in RMSD and RMSF values, high Rg and hydrogen bonds in mutant proteins compared to wild-type revealed that these mutations might alter the protein structure and stability of the corresponding protein. Hence, the results from the analyses guide the exploration of the mechanism by which these missense SNPs of the MCM6 gene alter the structural integrity and functional properties of the protein, which could guide the identification of ways to minimize the harmful effects of these mutations in humans.

Mapping protein-protein interactions by mass spectrometry

Protein-protein interactions (PPIs) are essential for numerous biological activities, including signal transduction, transcription control, and metabolism. They play a pivotal role in the organization and function of the proteome, and their perturbation is associated with various diseases, such as cancer, neurodegeneration, and infectious diseases. Recent advances in mass spectrometry (MS)-based protein interactomics have significantly expanded our understanding of the PPIs in cells, with techniques that continue to improve in terms of sensitivity, and specificity providing new opportunities for the study of PPIs in diverse biological systems. These techniques differ depending on the type of interaction being studied, with each approach having its set of advantages, disadvantages, and applicability. This review highlights recent advances in enrichment methodologies for interactomes before MS analysis and compares their unique features and specifications. It emphasizes prospects for further improvement and their potential applications in advancing our knowledge of PPIs in various biological contexts.

Genetic Modifiers of Sickle Cell Anemia Phenotype in a Cohort of Angolan Children

The aim of this study was to identify genetic markers in the HBB Cluster; HBS1L-MYB intergenic region; and BCL11A, KLF1, FOX3, and ZBTB7A genes associated with the heterogeneous phenotypes of Sickle Cell Anemia (SCA) using next-generation sequencing, as well as to assess their influence and prevalence in an Angolan population. Hematological, biochemical, and clinical data were considered to determine patients' severity phenotypes. Samples from 192 patients were sequenced, and 5,019,378 variants of high quality were registered. A catalog of candidate modifier genes that clustered in pathophysiological pathways important for SCA was generated, and candidate genes associated with increasing vaso-occlusive crises (VOC) and with lower fetal hemoglobin (HbF) were identified. These data support the polygenic view of the genetic architecture of SCA phenotypic variability. Two single nucleotide polymorphisms in the intronic region of 2q16.1, harboring the BCL11A gene, are genome-wide and significantly associated with decreasing HbF. A set of variants was identified to nominally be associated with increasing VOC and are potential genetic modifiers harboring phenotypic variation among patients. To the best of our knowledge, this is the first investigation of clinical variation in SCA in Angola using a well-customized and targeted sequencing approach.

Skeletal muscle of young females under resistance exercise exhibits a unique innate immune cell infiltration profile compared to males and elderly individuals

Muscle damage resulting from physical activities such as exercise triggers an immune response crucial for tissue repair and recovery. This study investigates the immune cell profiles in muscle biopsies of individuals engaged in resistance exercise (RE) and explores the impact of age and sex on the immune response following exercise-induced muscle damage. Microarray datasets from muscle biopsies of young and old subjects were analyzed, focusing on the gene expression patterns associated with immune cell activation. Genes were compared with immune cell signatures to reveal the cellular landscape during exercise. Results show that the most significant modulated gene after RE was Folliculin Interacting Protein 2 (FNIP2) a crucial regulator in cellular homeostasis. Moreover, the transcriptome was stratified based on the expression of FNIP2 and the 203 genes common to the groups obtained based on sex and age. Gene ontology analysis highlighted the FLCN-FNIP1-FNIP2 complex, which exerts as a negative feedback loop to Pi3k-Akt-mTORC1 pathway. Furthermore, we highlighted that the young females exhibit a distinct innate immune cell activation signature compared to males after a RE session. Specifically, young females demonstrate a notable overlap with dendritic cells (DCs), M1 macrophages, M2 macrophages, and neutrophils, while young males overlap with M1 macrophages, M2 macrophages, and motor neurons. Interestingly, in elderly subjects, both sexes display M1 macrophage activation signatures. Comparison of young and elderly signatures reveals an increased M1 macrophage percentage in young subjects. Additionally, common genes were identified in both sexes across different age groups, elucidating biological functions related to cell remodeling and immune activation. This study underscores the intricate interplay between sex, age, and the immune response in muscle tissue following RE, offering potential directions for future research. Nevertheless, there is a need for further studies to delve deeper and confirm the dynamics of immune cells in response to exercise-induced muscle damage.

Evolutionary conservation of putative suicidality-related risk genes that produce diminished motivation corrected by clozapine, lithium and antidepressants

Background

Genome wide association studies (GWAS) and candidate gene analyses have identified genetic variants and genes that may increase the risk for suicidal thoughts and behaviors (STBs). Important unresolved issues surround these tentative risk variants such as the characteristics of the associated genes and how they might elicit STBs.

Methods

Putative suicidality-related risk genes (PSRGs) were identified by comprehensive literature search and were characterized with respect to evolutionary conservation, participation in gene interaction networks and associated phenotypes. Evolutionary conservation was established with database searches and BLASTP queries, whereas gene-gene interactions were ascertained with GeneMANIA. We then examined whether mutations in risk-gene counterparts in C. elegans produced a diminished motivation phenotype previously connected to suicide risk factors.

Results and conclusions

From the analysis, 105 risk-gene candidates were identified and found to be: 1) highly conserved during evolution, 2) enriched for essential genes, 3) involved in significant gene-gene interactions, and 4) associated with psychiatric disorders, metabolic disturbances and asthma/allergy. Evaluation of 17 mutant strains with loss-of-function/deletion mutations in PSRG orthologs revealed that 11 mutants showed significant evidence of diminished motivation that manifested as immobility in a foraging assay. Immobility was corrected in some or all of the mutants with clozapine, lithium and tricyclic antidepressant drugs. In addition, 5-HT2 receptor and muscarinic receptor antagonists restored goal-directed behavior in most or all of the mutants. These studies increase confidence in the validity of the PSRGs and provide initial clues about possible mechanisms that mediate STBs.

Identification of the shared genes in type 2 diabetes mellitus and osteoarthritis and the role of quercetin

This study investigated the underlying comorbidity mechanism between type 2 diabetes mellitus (T2DM) and osteoarthritis (OA), while also assessing the therapeutic potential of quercetin for early intervention and treatment of these two diseases. The shared genes were obtained through GEO2R, limma and weighted gene co-expression network analysis (WGCNA), and validated using clinical databases and the area under the curves (ROC). Functional enrichment analysis was conducted to elucidate the underlying mechanisms of comorbidity between T2DM and OA. The infiltration of immune cells was analysed using the CIBERSORT algorithm in conjunction with ESTIMATE algorithm. Subsequently, transcriptional regulation analysis, potential chemical prediction, gene-disease association, relationships between the shared genes and ferroptosis as well as immunity-related genes were investigated along with molecular docking. We identified the 12 shared genes (EPHA3, RASIP1, PENK, LRRC17, CEBPB, EFEMP2, UBAP1, PPP1R15A, SPEN, MAFF, GADD45B and KLF4) across the four datasets. Our predictions suggested that targeting these shared genes could potentially serve as therapeutic interventions for both T2DM and OA. Specifically, they are involved in key signalling pathways such as p53, IL-17, NF-kB and MAPK signalling pathways. Furthermore, the regulation of ferroptosis and immunity appears to be interconnected in both diseases. Notably, in this context quercetin emerges as a promising drug candidate for treating T2DM and OA by specifically targeting the shared genes. We conducted a bioinformatics analysis to identify potential therapeutic targets, mechanisms and drugs for T2DM and OA, thereby offering novel insights into molecular therapy for these two diseases.

Identification of basement membrane-related signatures for estimating prognosis, immune infiltration landscape and drug candidates in pancreatic adenocarcinoma

Background: Pancreatic adenocarcinoma (PAAD) is a frequent digestive system cancer, which has high mortality and bad outcome. However, the role of basement membrane (BM)-related gene in assessing patient's outcome, microenvironment and treatment response remain unclear. Methods: Basement membrane (BM)-associated genes were detected by univariate and least absolute shrinkage and selection operator (LASSO) Cox regression analyses using data from the TCGA databases. A risk score system was constructed to distinguish patients in the high- and low-risk groups. Prognostic gene distribution in various immune cell forms was explored through scRNA-seq. Immune cell infiltration was assessed using CIBERSORT and ESTIMATE. The IC50 of common chemotherapeutic drugs and useful molecule compounds were evaluated. The mRNA and protein expression of important signatures were validated utilizing GEPIA and HPA databases. Results: Compared to low risk PAAD patients, PAAD patients with high risk showed a significant much worse overall survival (OS). Risk score of BM-associated genes could estimate patient outcome well, and areas under the curve (AUC) of receiver operating characteristic (ROC) survival curve were 0.76, 0.85, and 0.85 at 1-, 3-, and 5-year. Clinical impact curve (CIC) curve demonstrated the clinical importance of risk score. scRNA-seq revealed that BM-related genes were mainly distributed in malignant cells. Significant variations existed in the immune microenvironment, immune checkpoint expression and chemotherapy response between the studied groups. Furthermore, the mRNA expression levels of FAM83A, LY6D, MET, MUC16, MYEOV, and WNT7A were elevated in PAAD tissues, while the protein expression patterns of LY6D, MET, MUC16, and WNT7A were higher in tumor sample. RO-90-7501, Scriptaid, TG-101348, XMD-892, and XMD-1150 may be valuable small molecule drugs to treat PAAD. Conclusions: In conclusion, we develop a novel BM-related gene signature provide new insights and targets for the diagnosis, outcome estimation, candidate drugs and therapy management of PAAD patients.

Involvement of toxic metals and PCBs mixture in the thyroid and male reproductive toxicity: In silico toxicogenomic data mining

Toxicological research is mostly limited to considering the effects of a single substance, even though the real exposure of people is reflected in their daily exposure to many different chemical substances in low-doses. This in silico toxicogenomic study aims to provide evidence for the selected environmental (organo)metals (lead, cadmium, methyl mercury) + polychlorinated biphenyls mixture involvement in the possible alteration of thyroid, and male reproductive system function, and furthermore to predict the possible toxic mechanisms of the environmental cocktail. The Comparative Toxicogenomic Database, GeneMANIA online software, and ToppGene Suite portal were used as the main tools for toxicogenomic data mining and gene ontology analysis. The results show that 35 annotated common genes between selected chemicals and endocrine system diseases can interact on the co-expression level. Our study highlighted the disruption of the cytokines, the cell's response to oxidative stress, and the influence of the transcription factors as the potential core of toxicological mechanisms of the discussed mixture's effects. The connected toxicological effects of the tested mixture were abnormal sperm cells, a disrupted level of testosterone, and thyroid hormones. The core mechanisms of these effects were inflammation, oxidative stress, disruption of androgen receptor signaling, and the alteration of the FOXO3-Keap-1/NRF2-HMOX1-NQO1 pathway signaling most likely controlled by the co-expression of overlapped genes among used chemicals. This in silico research can be used as a potential core for the determination of biomarkers that can be monitored in future further in vitro and in vivo experiments.

Evolutionary dynamics of genome size and content during the adaptive radiation of Heliconiini butterflies

Heliconius butterflies, a speciose genus of Müllerian mimics, represent a classic example of an adaptive radiation that includes a range of derived dietary, life history, physiological and neural traits. However, key lineages within the genus, and across the broader Heliconiini tribe, lack genomic resources, limiting our understanding of how adaptive and neutral processes shaped genome evolution during their radiation. Here, we generate highly contiguous genome assemblies for nine Heliconiini, 29 additional reference-assembled genomes, and improve 10 existing assemblies. Altogether, we provide a dataset of annotated genomes for a total of 63 species, including 58 species within the Heliconiini tribe. We use this extensive dataset to generate a robust and dated heliconiine phylogeny, describe major patterns of introgression, explore the evolution of genome architecture, and the genomic basis of key innovations in this enigmatic group, including an assessment of the evolution of putative regulatory regions at the Heliconius stem. Our work illustrates how the increased resolution provided by such dense genomic sampling improves our power to generate and test gene-phenotype hypotheses, and precisely characterize how genomes evolve.

Testing sulforaphane as a strategy against toxic chemicals of public health concern by toxicogenomic data analysis: Friend or foe at the gene level - Colorectal carcinoma case study

Toxic metals (cadmium (Cd), lead (Pb), mercury (Hg) and arsenic (As)) and plastificators (bis (2 - ethylhexyl) phthalate (DEHP), dibutyl phthalate (DBP)) and bisphenol A (BPA)) have been suggested to aid in colorectal carcinoma (CRC) advancement. Sulforaphane (SFN), isothiocyanate from cruciferous vegetables, diminishes chemical carcinogenesis susceptibility, but has been shown to act as a friend or a foe depending on various factors. By conducting the mechanistic toxicogenomic data mining approach, this research aimed to determine if SFN can alleviate toxic-metal and/or phthalate/BPA mixture-induced CRC at the gene level. Comparative Toxicogenomics Database, ToppGene Suite portal, Cytoscape software, InteractiVenn and Gene Expression Omnibus (GEO) database (GEO2R tool) was used. Among the mutual genes for all the investigated substances, SFN had a protective impact only through PTGS2. Other proposed protective SFN-targets included ABCA1, ALDH2, BMP2, DPYD, MYC, SLCO2A1, and SOD2, only in the case of phthalates/BPA exposure. The only additional gene relevant for SFN protection against the toxic metal mixture-induced CRC was ABCB1. Additionally, the majority of the top 15 molecular pathways extracted for SFN impact on phthalate and BPA mixture-linked CRC development were directly linked with cancer development, which was not the case with the toxic metal mixture. The current research has indicated that SFN is a more effective chemoprotective agent against CRC induced by phthalates/BPA mixture than by toxic-metal mixture. It has also presented the value of computational methods as a simple tool for directing further research, selecting appropriate biomarkers and exploring the mechanisms of toxicity.

In silico assessment of mixture toxicity mechanisms involved in the pathogenesis of thyroid diseases: the combination of toxic metal(oid)s and decabrominated diphenyl ether

The current study aimed to assess the connection between the mixture of lead (Pb), cadmium (Cd), arsenic (As), methylmercury (MeHg) and decabrominated diphenyl ether (decaBDE) and thyroid function, by using in silico toxicogenomic data-mining approach. To obtain the linkage between investigated toxic mixture and thyroid diseases (TDs), the Comparative Toxicogenomics Database (CTD) was used, while gene ontology (GO) enrichment analysis was performed by ToppGeneSuite portal. The analysis has shown 10 genes connected to all chemicals present in the mixture and TDs (CAT, GSR, IFNG, IL1B, IL4, IL6, MAPK1, SOD2, TGFB1, TNF), most of which were in co-expression (45.68%), or belonged to the same pathway (30.47%). Top 5 biological processes and molecular functions affected by the investigated mixture emphasized the role of two common mechanisms - oxidative stress and inflammation. Cytokines and inflammatory response was listed as the main molecular pathway that may be triggered by simultaneous exposure to toxic metal(oid)s and decaBDE and connected to TDs. The direct relations between Pb/decaBDE and redox status impairment in thyroid tissue was confirmed by our chemical-phenotype interaction analysis, while the strongest linkage between Pb, As and decaBDE and thyroid disorders was found. The obtained results provide better understanding of molecular mechanisms involved in the thyrotoxicity of the investigated mixture, and can be used to direct further research.

A new view on functions of the lysine demalonylase activity of SIRT5

AIMS

The specificity of the lysine demalonylation substrates of the pharmaceutically attractive tumor promoter/suppressor SIRT5 is not comprehensively clarified. The present study re-analyses publicly available data and highlights potentially pharmaceutically interesting outcomes by the use of bioinformatics.

MATERIALS AND METHODS

The interaction networks of SIRT5 malonylome from the wild-type and ob/ob (obese pre-diabetic type) mice were subjected to the pathway enrichment and gene function prediction analysis using GeneMania (3.5.2) application run under Cytoscape (3.9.1) environment.

KEY FINDINGS

The analysis in the wild-type mice revealed the involvement of SIRT5 malonylome in Eukaryotic translation elongation (ETE; the nodes EF1A1, EEF2, EEF1D, and EEF1G), Amino acid and derivative metabolism (AADM), and Selenoamino acid metabolism (SAM). The tumor promoter/suppressor activity of SIRT5 is mediated through the tumor promoter substrates included in AADM (GLUD1, SHMT1, ACAT1), and the tumor suppressor substrates involved in AADM and SAM (ALDH9A1, BHMT, GNMT). Selen stimulates the expression of SIRT5 and other sirtuins. SIRT5 in turn regulates the selenocysteine synthesis, which creates a regulatory loop. The analysis of SIRT5 malonylome in pre-diabetic ob/ob mice identifies the mTORC1 pathway as a mechanism, which facilitates SIRT5 functions. The comparison of the outcomes of SIRT5 malonylome, succinylome, and glutarylome analysis disclosed several differences.

SIGNIFICANCE

The analysis showed additional aspects of SIRT5 malonylome functions besides the control of glucose metabolism. It defined several unique substrates and pathways, and it showed differences compared to other enzymatic activities of SIRT5, which could be used for pharmaceutical benefits.

From Classical to Modern Computational Approaches to Identify Key Genetic Regulatory Components in Plant Biology

The selection of plant genotypes with improved productivity and tolerance to environmental constraints has always been a major concern in plant breeding. Classical approaches based on the generation of variability and selection of better phenotypes from large variant collections have improved their efficacy and processivity due to the implementation of molecular biology techniques, particularly genomics, Next Generation Sequencing and other omics such as proteomics and metabolomics. In this regard, the identification of interesting variants before they develop the phenotype trait of interest with molecular markers has advanced the breeding process of new varieties. Moreover, the correlation of phenotype or biochemical traits with gene expression or protein abundance has boosted the identification of potential new regulators of the traits of interest, using a relatively low number of variants. These important breakthrough technologies, built on top of classical approaches, will be improved in the future by including the spatial variable, allowing the identification of gene(s) involved in key processes at the tissue and cell levels.

Prioritization of New Candidate Genes for Rare Genetic Diseases by a Disease-Aware Evaluation of Heterogeneous Molecular Networks

Screening for pathogenic variants in the diagnosis of rare genetic diseases can now be performed on all genes thanks to the application of whole exome and genome sequencing (WES, WGS). Yet the repertoire of gene-disease associations is not complete. Several computer-based algorithms and databases integrate distinct gene-gene functional networks to accelerate the discovery of gene-disease associations. We hypothesize that the ability of every type of information to extract relevant insights is disease-dependent. We compiled 33 functional networks classified into 13 knowledge categories (KCs) and observed large variability in their ability to recover genes associated with 91 genetic diseases, as measured using efficiency and exclusivity. We developed GLOWgenes, a network-based algorithm that applies random walk with restart to evaluate KCs' ability to recover genes from a given list associated with a phenotype and modulates the prediction of new candidates accordingly. Comparison with other integration strategies and tools shows that our disease-aware approach can boost the discovery of new gene-disease associations, especially for the less obvious ones. KC contribution also varies if obtained using recently discovered genes. Applied to 15 unsolved WES, GLOWgenes proposed three new genes to be involved in the phenotypes of patients with syndromic inherited retinal dystrophies.

Comparative analysis of differentially expressed miRNAs in leaves of three sugarcanes (Saacharum officinarum L.) cultivars during salinity stress

BACKGROUND

Sugarcane is an important industrial plant cultivated mostlyin the arid and semiarid regions. Due to climate change and anthropogenic activities, the sugarcane fieldsare prone to be damagedas a result of salt deposition. The consequence of such phenomena is turning to become a major thread in sugarcane cultivation. To address this issue, the identification of salinity tolerant cultivars would be a suitable strategy to minimize yield loss in the area. It is well known thatthe expression of abiotic stress-responsive genes including noncoding microRNAs (miRNAs) and their codingtargetscould lead to enhancement of stress tolerance in crops. Therefore, the expression study of those noncoding and coding genes under stress conditions is an appropriate approach to screen the tolerant cultivars. In addition, the examination of the expression of miRNA's target genes could provide deeper insight into the molecular stress mechanism and facilitate the identification of tolerant cultivars.

METHODS AND RESULTS

We aimedto assess the expression of nine candidate miRNAsand their corresponding targeted genes among the studied sugarcane cultivars under high salinity conditions, leading to the identification of the salt-tolerant cultivar. To achieve our goal, a two-factorial experiment with three sugarcane cultivars (CP-48, CP-57, CP-69) and two salinity levels (0 and 8 ds/m) was conducted. The result indicated significant differences in expression with in miRNAs and also their target genes. The highest reduction of miRNAs expression occurred in miR160 while the lowest oneappeared in miR1432. The data also indicated that the higher and the lowest expression of targeted genes occurred in miR160 and miR393 respectively. Among studied cultivars, the CP-57 showed poor performance while CP-69 expresses a superior tolerance to salt stress.

CONCLUSIONS

Taken together, these results suggested that the monitoring of microRNA expressioncould provide a new approach forthe screening of well-adapted cultivars under salt conditions. Such an approach would be the appropriate solutionto combat plant stress inhigh salinity region/soil. Our result indicated that the miR160 generates sugarcane tolerant to salt stress, can be potentially be used as a biomarker to salt stress.

In silico and in vitro analysis of microRNAs with therapeutic potential in atherosclerosis

Atherosclerosis is a chronic inflammatory disease in which aberrant lipid metabolism plays a key role. MicroRNAs (miRNAs), micro-coordinators of gene expression, have been recently proposed as novel clinical biomarkers and potential therapeutic tools for a broad spectrum of diseases. This study aimed to identify miRNAs with therapeutic potential in atherosclerosis. Bioinformatic databases, including experimentally validated and computational prediction tools as well as a novel combination method, were used to identify miRNAs that are able to simultaneously inhibit key genes related to the pathogenesis of atherosclerosis. Further validation of genes and miRNAs was conducted using the STRING online tool, KEGG pathway analysis and DIANA-miRPath. The inhibitory effects of the identified miRNAs in HepG2 and Huh7 cells were verified by real-time PCR. The MTT assay was utilized to evaluate cell cytotoxicity effects of miRNAs. Atherosclerotic drug-targeted genes were selected as key genes. Strong interactions between genes were confirmed using STRING. These genes were shown to be integral to critical pathological processes involved in atherosclerosis. A novel combined method of validated and predicted tools for the identification of effective miRNAs was defined as the combination score (C-Score). Bioinformatic analysis showed that hsa-miR-124-3p and hsa-miR-16-5p possessed the best C-Score (0.68 and 0.62, respectively). KEGG and DIANA-miRPath analysis showed that selected genes and identified miRNAs were involved in atherosclerosis-related pathways. Compared with the controls in both HepG2 and Huh7 cell lines, miR-124 significantly reduced the expression of CETP, PCSK9, MTTP, and APOB, and miR-16 significantly reduced the expression of APOCIII, CETP, HMGCR, PCSK9, MTTP, and APOB, respectively. The cytotoxicity assay showed that miR-124 reduced cell viability, especially after 72 h; however, miR-16 did not show any significant cytotoxicity in either cell line. Our findings indicate that hsa-miR-124 and miR-16 have potential for use as therapeutic candidates in the treatment of atherosclerosis.

Gene Co-Expression Network Modular Analysis Reveals Altered Immune Mechanisms in HIV-HAND

Although the introduction of HAART has completely changed the natural course of HIV infection, the number of chronic forms of HIV-associated neurocognitive disorder (HAND) has risen. It is estimated that up to half of subjects undergoing HAART therapy exhibit mild cognitive impairments. In the current study, we apply the gene co-expression network modular analysis, a well-established system biology approach, to the gene expression profiles of cases from the National NeuroAIDS Tissue Consortium (NNTC). We observed a negative enrichment for genes associated with the control of immune responses and putatively regulated by the transcription factors IRF8 and SPI1 and by both type I and II interferons. Our study provides evidence of altered immune responses, which are likely associated with the occurrence of HAND in the absence of HIV encephalitis (HIVE).

The genetics of occupational asthma development among workers exposed to diisocyanates: A systematic literature review with meta-analysis

Diisocyanates are widely used compounds that pose a safety concern for workers in occupations within the spray-paint, spray-foam insulation, and furniture varnish industries. Epidemiological studies show that only a subset of workers exposed to diisocyanates develop diisocyanate-induced occupational asthma (diisocyanate asthma, DA), indicating that genetic susceptibility may play a role. The purpose of this systematic literature review was to compile and meta-analyze the reported data on genetic susceptibility markers for DA. Three databases (Embase, Pubmed, and Scopus) were searched and 169 non-duplicate publications were identified, of which 22 relevant occupational studies were included in this review. Researchers reported prevalence odds ratios (PORs) for 943 comparisons in 82 different genes/serotypes. Protein network functions for the DA-associated genes from this review include: antigen processing, lymphocyte activation, cytokine production regulation, and response to oxidative stress. Meta-analysis of comparisons between workers with DA and controls was conducted for 23 genetic markers within: CTNNA3, GSTM1, GSTP1, GSTT1, HLA-C, HLA-DQB1, HLA-DR1, HLA-DR3, HLA-DR4, HLA-DR7, and HLA-DR8. These genes code for proteins that are involved in cell-cell adhesions (CTNNA3), glutathione conjugation for xenobiotic metabolism (GST gene family), and immune system response (HLA gene family). The most compelling pooled PORs were for two studies on CTNNA3 (increased DA risk: rs10762058 GG, rs7088181 GG, rs4378283 TT; PORs 4.38–4.97) and three studies on HLA-DR1 (decreased DA risk, POR 0.24). Bioinformatics of the predicted protein pathways for DA shows overlap with biomarker-associated pathways in workers before development of asthma, suggesting overlap in toxicokinetic and toxicodynamic pathways of diisocyanates. The control groups were also compared against each other and differences were negligible. Suggestions for improving future research are also presented. Of the highest importance, the literature was found to be profoundly publication-biased, in which researchers need to report the data for all studied markers regardless of the statistical significance level. We demonstrate the utility of evaluating the overlap in predicted protein pathway functions for identifying more consistency across the reported literature including for asthma research, biomarker research, and in vitro studies. This will serve as an important resource for researchers to use when generating new hypothesis-driven research about diisocyanate toxicology.

Metformin regulates multiple signaling pathways within castration-resistant human prostate cancer cells

Background

The biguanide metformin has been shown to not only reduce circulating glucose levels but also suppress in vitro and in vivo growth of prostate cancer. However, the mechanisms underlying the anti-tumor effects of metformin in advanced prostate cancers are not fully understood. The goal of the present study was to define the signaling pathways regulated by metformin in androgen-receptor (AR) positive, castration-resistant prostate cancers.

Methods

Our group used RNA sequencing (RNA-seq) to examine genes regulated by metformin within the C4–2 human prostate cancer cell line. Western blot analysis and quantitative RT-PCR were used to confirm alterations in gene expression and further explore regulation of protein expression by metformin.

Results

Data from the RNA-seq analysis revealed that metformin alters the expression of genes products involved in metabolic pathways, the spliceosome, RNA transport, and protein processing within the endoplasmic reticulum. Gene products involved in ErbB, insulin, mTOR, TGF-β, MAPK, and Wnt signaling pathways are also regulated by metformin. A subset of metformin-regulated gene products were genes known to be direct transcriptional targets of p53 or AR. Western blot analyses and quantitative RT-PCR indicated these alterations in gene expression are due in part to metformin-induced reductions in AR mRNA and protein levels.

Conclusions

Together, our results suggest metformin regulates multiple pathways linked to tumor growth and progression within advanced prostate cancer cells.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-022-10115-3.

A sex-stratified analysis of neuroimmune gene expression signatures in Alzheimer's disease brains

Alzheimer's disease (AD) is the most common form of progressively disabling dementia. The chitinases CHI3L1 and CHI3L2 have long been known as biomarkers for microglial and astrocytic activation in neurodegeneration. Here, we collected microarray datasets from the National Center for Biotechnology Information (NCBI) brain samples of non-demented controls (NDC) (n = 460), and of deceased patients with AD (n = 697). The AD patients were stratified according to sex. Comparing the high CHI3L1 and CHI3L2 expression group (75^th percentile), and low CHI3L1 and CHI3L2 expression group (25^th percentile), we obtained eight signatures according to the sex of patients and performed a genomic deconvolution analysis using neuroimmune signatures (NIS) belonging to twelve cell populations. Expression analysis revealed significantly higher CHI3L1 and CHI3L2 expression in AD compared with NDC, and positive correlations of these genes with GFAP and TMEM119. Furthermore, deconvolution analysis revealed that CHI3L1 and CHI3L2 high expression was associated with inflammatory signatures in both sexes. Neuronal activation profiles were significantly activated in AD patients with low CHI3L1 and CHI3L2 expression levels. Furthermore, gene ontology analysis of common genes regulated by the two chitinases unveiled immune response as a main biological process. Finally, microglia NIS significantly correlated with CHI3L2 expression levels and were more than 98% similar to microglia NIS determined by CHI3L1. According to our results, high levels of CHI3L1 and CHI3L2 in the brains of AD patients are associated with inflammatory transcriptomic signatures. The high correlation between CHI3L1 and CHI3L2 suggests strong co-regulation.

Comprehensive analysis identifies as a critical prognostic prediction gene in breast cancer

BACKGROUND

Aurora kinases (AURKs) family plays a vital role not only in cell division but also in tumorigenesis. However, there are still rare systematic analyses of the diverse expression patterns and prognostic value of the AURKs family in breast cancer (BC). Systematic bioinformatics analysis was conducted to explore the biological role, prognostic value, and immunologic function of AURKs family in BC.

METHODS

The expression, prognostic value, and clinical functions of AURKs family in BC were evaluated with several bioinformatics web portals: ONCOMINE Gene Expression Profiling Interactive Analysis, Kaplan-Meier plotter, cBioPortal, Metascape, GeneMANIA, and LinkedOmics; and the result was verified using human tissues.

RESULTS

The expression of AURKA and AURKB were upregulated in BC in subgroup analyses based on tumor stage (all P   <  0.05). BC patients with high AURKA and AURKB expression had a worse overall survival, relapse-free survival, and distant metastasis-free survival (all P   <  0.05). Verification experiment revealed that AURKA and AURKB were upregulated in BC ( P  < 0.05). AURKA and AURKB were specifically associated with several tumor-associated kinases (polo-like kinase 1 and cyclin-dependent kinase 1), miRNAs (miR-507 and miR-381), and E2F transcription factor 1. Moreover, AURKA and AURKB were correlated with immune cell infiltration. Functional enrichment analysis revealed that AURKA and AURKB were involved in the cell cycle signaling pathway, platinum drug resistance signaling pathway, ErbB signaling pathway, Hippo signaling pathway, and nucleotide-binding and oligomerization domain-like receptor signaling pathway.

CONCLUSIONS

Aurora kinases AURKA and AURKB could be employed as novel prognostic biomarkers or promising therapeutic targets for BC.

GJA1/CX43 High Expression Levels in the Cervical Spinal Cord of ALS Patients Correlate to Microglia-Mediated Neuroinflammatory Profile

Amyotrophic lateral sclerosis (ALS) is a progressive neurodegenerative disorder affecting motoneurons (MNs) with a fatal outcome. The typical degeneration of cortico-spinal, spinal, and bulbar MNs, observed in post-mortem biopsies, is associated with the activation of neuroimmune cells. GJA1, a member of the connexins (Cxs) gene family, encodes for connexin 43 (Cx43), a core gap junctions (GJs)- and hemichannels (HCs)-forming protein, involved in cell death, proliferation, and differentiation. Recently, Cx43 expression was found to play a role in ALS pathogenesis. Here, we used microarray and RNA-seq datasets from the NCBI of the spinal cord of control (NDC) and ALS patients, which were stratified according to the GJA1 gene expression. Genes that positively or negatively correlated to GJA1 expression were used to perform a genomic deconvolution analysis (GDA) using neuroimmune signatures. Expression analysis revealed a significantly higher GJA1 expression in the MNs of ALS patients as compared to NDC. Gene deconvolution analysis revealed that positively correlated genes were associated with microglia activation, whereas negatively correlated genes were associated with neuronal activation profiles. Moreover, gene ontology analysis, performed on genes characterizing either microglia or neuronal signature, indicated immune activation or neurogenesis as main biological processes. Finally, using a synthetic analysis of drugs able to revert the GJA1 transcriptomic signatures, we found a specific drug profile for ALS patients with high GJA1 expression levels, composed of amlodipine, sertraline, and prednisolone. In conclusion, our exploratory study suggests GJA1 as a new neuro-immunological gene correlated to microglial cellular profile in the spinal cord of ALS patients. Further studies are warranted to confirm these results and to evaluate the therapeutic potential of drugs able to revert typical GJA1/CX43 signature in ALS patients

Activation of CCL21-GPR174/CCR7 on cardiac fibroblasts underlies myocardial ischemia/reperfusion injury

Background: The mechanisms underlying myocardial ischemia/reperfusion (I/R) injury are not fully understood. This study aims to explore key candidate genes and potential therapeutic targets for treatment of myocardial I/R injury.

Methods: The transcriptional profiles of ventricular myocardium during cardiac arrest, ischemia, and reperfusion were obtained from the Gene Expression Omnibus database. Based on the transcriptional data of GSE6381, functional pathway and process enrichment analyses, protein–protein interaction network, and gene set enrichment analyses were conducted. In the animal experiments, we established the myocardial I/R injury model in mice. We validated the mRNA and protein expression of the key genes using the qPCR and western blots. We further assessed the expression and localization of CCL21 and its receptors using immunofluorescence staining experiments.

Results: The microarray analyses identified five key genes (CCL21, XCR1, CXCL13, EDN1, and CASR). Myocardial I/R process in mice resulted in significant myocardial infraction, histological damage, and myocardial apoptosis. The results of qPCR and western blots showed that the expression of CCL21 and CXCL13 were increased following myocardial I/R injury in mice. Furthermore, the immunofluorescence staining results revealed that the expression of GPR174/CCR7 (CCL21 receptors), but not CXCR5 (CXCL13 receptor), was elevated following myocardial I/R injury. Moreover, the activated CCL21-GPR174/CCR7 signaling was located on the cardiac fibroblasts of the myocardium with I/R injury.

Conclusion: This study revealed several key factors underlying myocardial I/R injury. Of these, the activation of CCL21-GPR174/CCR7 signaling on cardiac fibroblasts was highlighted, which provides potential therapeutic targets for cardioprotection.

Nine quick tips for pathway enrichment analysis

Pathway enrichment analysis (PEA) is a computational biology method that identifies biological functions that are overrepresented in a group of genes more than would be expected by chance and ranks these functions by relevance. The relative abundance of genes pertinent to specific pathways is measured through statistical methods, and associated functional pathways are retrieved from online bioinformatics databases. In the last decade, along with the spread of the internet, higher availability of computational resources made PEA software tools easy to access and to use for bioinformatics practitioners worldwide. Although it became easier to use these tools, it also became easier to make mistakes that could generate inflated or misleading results, especially for beginners and inexperienced computational biologists. With this article, we propose nine quick tips to avoid common mistakes and to out a complete, sound, thorough PEA, which can produce relevant and robust results. We describe our nine guidelines in a simple way, so that they can be understood and used by anyone, including students and beginners. Some tips explain what to do before starting a PEA, others are suggestions of how to correctly generate meaningful results, and some final guidelines indicate some useful steps to properly interpret PEA results. Our nine tips can help users perform better pathway enrichment analyses and eventually contribute to a better understanding of current biology.

Characterization of methylation patterns associated with lifestyle factors and vitamin D supplementation in a healthy elderly cohort from Southwest Sweden

Numerous studies have shown that lifestyle factors, such as regular physical activity and vitamin D intake, may remarkably improve overall health and mental wellbeing. This is especially important in older adults whose vitamin D deficiency occurs with a high prevalence. This study aimed to examine the influence of lifestyle and vitamin D on global DNA methylation patterns in an elderly cohort in Southwest of Sweden. We also sought to examine the methylation levels of specific genes involved in vitamin D's molecular and metabolic activated pathways. We performed a genome wide methylation analysis, using Illumina Infinium DNA Methylation EPIC 850kBeadChip array, on 277 healthy individuals from Southwest Sweden at the age of 70–95. The study participants also answered queries on lifestyle, vitamin intake, heart medication, and estimated health. Vitamin D intake did not in general affect methylation patterns, which is in concert with other studies. However, when comparing the group of individuals taking vitamin supplements, including vitamin D, with those not taking supplements, a difference in methylation in the solute carrier family 25 (SCL25A24) gene was found. This confirms a previous finding, where changes in expression of SLC25A24 were associated with vitamin D treatment in human monocytes. The combination of vitamin D intake and high physical activity increased methylation of genes linked to regulation of vitamin D receptor pathway, the Wnt pathway and general cancer processes. To our knowledge, this is the first study detecting epigenetic markers associated with the combined effects of vitamin D supplementation and high physical activity. These results deserve to be further investigated in an extended, interventional study cohort, where also the levels of 25(OH)D₃ can be monitored.

Approaches in Gene Coexpression Analysis in Eukaryotes

Simple Summary

Genes whose expression levels rise and fall similarly in a large set of samples, may be considered coexpressed. Gene coexpression analysis refers to the en masse discovery of coexpressed genes from a large variety of transcriptomic experiments. The type of biological networks that studies gene coexpression, known as Gene Coexpression Networks, consist of an undirected graph depicting genes and their coexpression relationships. Coexpressed genes are clustered in smaller subnetworks, the predominant biological roles of which can be determined through enrichment analysis. By studying well-annotated gene partners, the attribution of new roles to genes of unknown function or assumption for participation in common metabolic pathways can be achieved, through a guilt-by-association approach. In this review, we present key issues in gene coexpression analysis, as well as the most popular tools that perform it.

Abstract

Gene coexpression analysis constitutes a widely used practice for gene partner identification and gene function prediction, consisting of many intricate procedures. The analysis begins with the collection of primary transcriptomic data and their preprocessing, continues with the calculation of the similarity between genes based on their expression values in the selected sample dataset and results in the construction and visualisation of a gene coexpression network (GCN) and its evaluation using biological term enrichment analysis. As gene coexpression analysis has been studied extensively, we present most parts of the methodology in a clear manner and the reasoning behind the selection of some of the techniques. In this review, we offer a comprehensive and comprehensible account of the steps required for performing a complete gene coexpression analysis in eukaryotic organisms. We comment on the use of RNA-Seq vs. microarrays, as well as the best practices for GCN construction. Furthermore, we recount the most popular webtools and standalone applications performing gene coexpression analysis, with details on their methods, features and outputs.

How Variation in Risk Allele Output and Gene Interactions Shape the Genetic Architecture of Schizophrenia

Schizophrenia is a highly heritable polygenic psychiatric disorder. Characterization of its genetic architecture may lead to a better understanding of the overall burden of risk variants and how they determine susceptibility to disease. A major goal of this project is to develop a modeling approach to compare and quantify the relative effects of single nucleotide polymorphisms (SNPs), copy number variants (CNVs) and other factors. We derived a mathematical model for the various genetic contributions based on the probability of expressing a combination of risk variants at a frequency that matched disease prevalence. The model included estimated risk variant allele outputs (VAOs) adjusted for population allele frequency. We hypothesized that schizophrenia risk genes would be more interactive than random genes and we confirmed this relationship. Gene–gene interactions may cause network ripple effects that spread and amplify small individual effects of risk variants. The modeling revealed that the number of risk alleles required to achieve the threshold for susceptibility will be determined by the average functional locus output (FLO) associated with a risk allele, the risk allele frequency (RAF), the number of protective variants present and the extent of gene interactions within and between risk loci. The model can account for the quantitative impact of protective variants as well as CNVs on disease susceptibility. The fact that non-affected individuals must carry a non-trivial burden of risk alleles suggests that genetic susceptibility will inevitably reach the threshold for schizophrenia at a recurring frequency in the population.

Differential gene expression profile of multinodular goiter

Introduction

The goiter, a neglected heterogeneous molecular disease, remains a major indication for thyroidectomies in its endemic regions.

Objectives

This study analyzed differential gene expression in surgical specimens diagnosed with multi nodular and compared the data to that of thyroid tissue without multinodular goiter from patients undergoing thyroidectomy in Manaus-AM, Brazil using RNA-seq technology.

Methodology

The transcriptome information of the surgical specimen fragments with and without multinodular goiter was accessed by Illumina HiSeq 2000 New Generation Sequencing (NGS) using the RNA-seq NEBNext^® Ultra^™ RNA Library Prep Kit for Illumina^®—#E7530L protocol and differential gene expression analysis.

Results

Differences were found between the gene expression profiles of the diseased tissues and those of the healthy control tissues; at least 70 genes were differentially expressed. The HOTS gene was expressed only in multinodular goiter tissues (p < 0.05).

Conclusion

These results demonstrate that the gene expression profile of multinodular goiter is pro-tumoral and that HOTS can play a central role in multinodular goiter development.

Integrated analysis of 14 lymphoma datasets revealed high expression of CXCL14 promotes cell migration in mantle cell lymphoma

Lymphoma is accompanied by the impairment of multiple immune functions. Cytokines play an important role in a variety of immune-related functions and affect the tumor microenvironment. However, the exact regulatory mechanisms between them remain unclear. This study aimed to explore the cytokines expression and function in Hodgkin's lymphoma (HL), diffuse large B-cell lymphoma (DLBCL), and mantle cell lymphoma (MCL). We performed a transcriptome integration analysis of 14 lymphoma datasets including 240 Hodgkin's lymphoma, 891 diffuse large B-cell lymphoma, 216 mantle cell lymphoma, and 64 health samples. The results showed that multiple immune functions and signal pathway damage were shared by all three types of lymphoma, and these functions were related to cytokines. Furthermore, through co-expression network and functional interaction network analysis, we identified CXCL14 as a key regulator and it affects cell chemotaxis and migration functions. The functional experiment showed that CXCL14 knockdown inhibited cell migration in MCL cell lines. This study suggested that high expression of CXCL14 may aggravate MCL via promoting cell migration. Our findings provide novel insights into the biology of this disease and would be helpful for the pathogenesis study and drug discovery of lymphomas.

NcRNA-Mediated High Expression of HMMR as a Prognostic Biomarker Correlated With Cell Proliferation and Cell Migration in Lung Adenocarcinoma

Background

Hyaluronan-mediated motility receptor (HMMR) plays a pivotal role in cell proliferation in various cancers, including lung cancer. However, its function and biological mechanism in lung adenocarcinoma (LUAD) remain unclear.

Methods

Data on HMMR expression from several public databases were extensively analyzed, including the prognosis of HMMR in the Gene Expression Profiling Interactive Analysis (GEPIA) database. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were analyzed using DAVID and gene set enrichment analysis (GSEA) software. The correlation between HMMR expression and immune cell infiltration was analyzed in the Tumor Immune Estimation Resource (TIMER) database, and the gene and protein networks were examined using the GeneMANIA and STRING databases. Experimentally, the expression of HMMR in LUAD and lung cancer cell lines was determined using immunohistochemistry and quantitative RT-PCR assays. Besides, the function of HMMR on cancer cell proliferation and migration was examined using cell growth curve and colony formation, Transwell, and wound healing assays.

Results

In this study, we found that HMMR was elevated in LUAD and that its high expression was associated with poor clinicopathological features and adverse outcomes in LUAD patients. Furthermore, our results demonstrated that the expression of HMMR was positively correlated with immune cell infiltration and immune modulation. Interestingly, diverse immune cell infiltration affects the prognosis of LUAD. In the functional assay, depletion of HMMR significantly repressed the cancer cell growth and migration of LUAD. Mechanically, we found that that the DNA methylation/TMPO-AS1/let-7b-5p axis mediated the high expression of HMMR in LUAD. Depletion of TMPO-AS1 and overexpression of let-7b-5p could result in the decreased expression of HMMR in LUAD cells. Furthermore, we found that TMPO-AS1 was positively correlated with HMMR, yet negatively correlated with let-7b-5p expression in LUAD.

Conclusions

Our findings elucidated that the DNA methylation/TMPO-AS1/let-7b-5p axis mediated the high expression of HMMR, which may be considered as a biomarker to predict prognosis in LUAD.

Finding New Ways How to Control BACE1

Recently, all applications of BACE1 inhibitors failed as therapeutical targets for Alzheimer´s disease (AD) due to severe side effects. Therefore, alternative ways for treatment development are a hot research topic. The present analysis investigates BACE1 protein-protein interaction networks and attempts to solve the absence of complete knowledge about pathways involving BACE1. A bioinformatics analysis matched the functions of the non-substrate interaction network with Voltage-gated potassium channels, which also appear as top priority protein nodes. Targeting BACE1 interactions with PS1 and GGA-s, blocking of BACE1 access to APP by BRI3 and RTN-s, activation of Wnt signaling and upregulation of β-catenin, and brain delivery of the extracellular domain of p75^NTR, are the main alternatives to the use of BACE 1 inhibitors highlighted by the analysis. The pathway enrichment analysis also emphasized substrates and substrate candidates with essential biological functions, which cleavage must remain controlled. They include ephrin receptors, ROBO1, ROBO2, CNTN-s, CASPR-s, CD147, CypB, TTR, APLP1/APLP2, NRXN-s, and PTPR-s. The analysis of the interaction subnetwork of BACE1 functionally related to inflammation identified a connection to three cardiomyopathies, which supports the hypothesis of the common molecular mechanisms with AD. A lot of potential shows the regulation of BACE1 activity through post-translational modifications. The interaction network of BACE1 and its phosphorylation enzyme CSNK1D functionally match the Circadian clock, p53, and Hedgehog signaling pathways. The regulation of BACE1 glycosylation could be achieved through N-acetylglucosamine transferases, α-(1→6)-fucosyltransferase, β-galactoside α-(2→6)-sialyltransferases, galactosyltransferases, and mannosidases suggested by the interaction network analysis of BACE1-MGAT3. The present analysis proposes possibilities for the alternative control of AD pathology.

PathwAX II: Network-based pathway analysis with interactive visualization of network crosstalk

Motivation

Pathway annotation tools are indispensable for the interpretation of a wide range of experiments in life sciences. Network-based algorithms have recently been developed which are more sensitive than traditional overlap-based algorithms, but there is still a lack of good online tools for network-based pathway analysis.

Results

We present PathwAX II—a pathway analysis web tool based on network crosstalk analysis using the BinoX algorithm. It offers several new features compared with the first version, including interactive graphical network visualization of the crosstalk between a query gene set and an enriched pathway, and the addition of Reactome pathways.

Availability and implementation

PathwAX II is available at http://pathwax.sbc.su.se.

Supplementary information

Supplementary data are available at Bioinformatics online.

A Novel Role of Arrhythmia-Related Gene KCNQ1 Revealed by Multi-Omic Analysis: Theragnostic Value and Potential Mechanisms in Lung Adenocarcinoma

The early diagnosis, prognostic prediction, and personalized therapy of lung adenocarcinoma (LUAD) remains a challenging issue. KCNQ1 (potassium voltage-gated channel subfamily Q Member 1) is implicated in long QT syndrome (LQTS) and cardiac arrhythmia, while its significance in LUAD remains unclear. In this study, we aimed to explore the significance of KCNQ1 in terms of clinical value, tumor immunity, underlying mechanisms, and a precision medicine approach by means of multi-omics analysis. The association of KCNQ1 with LUAD was first explored. Both altered variants and high expression of KCNQ1 in a TCGA-LUAD cohort indicated a favorable outcome. KCNQ1 levels had a negative correlation with tumor proliferation index Ki67 levels. siRNA-knockdown of KCNQ1 promoted the migration ability of lung cancer cells. KCNQ1 levels were decreased in LUAD tissue compared to normal tissue. A receiver operating characteristic (ROC) curve indicated good diagnostic efficiency of KCNQ1. High KCNQ1 is associated with an immunoactive profile of immune infiltration and immunomodulators and is involved in the inhibition of the cell cycle and DNA replication. Lapatinib was identified as a potent drug for LUAD in the context of low KCNQ1. This study unveiled the significance of KCNQ1 in diagnosis and prognosis and provided a corresponding precision medicine strategy for LUAD.

FOXM1/lncRNA TYMSOS/miR-214-3p–Mediated High Expression of NCAPG Correlates With Poor Prognosis and Cell Proliferation in Non–Small Cell Lung Carcinoma

Lung cancer is the most common cancer with high mortality. Increasing evidence has demonstrated that nonstructural maintenance of chromosomes condensin I complex subunit G (NCAPG) plays a crucial role in the progression of human cancers. However, the biological function and underlying mechanism of NCAPG in non–small cell lung cancer (NSCLC) are still unclear. Here, we utilized diverse public databases to analyze the expression of NCAPG in pan-cancer. We found that NCAPG was highly expressed in various human cancers, especially in NSCLC. NCAPG expression was significantly positively correlated with poor clinical-pathological features, poor prognosis, tumor mutational burden, DNA microsatellite instability, and immune cell infiltration in NSCLC. In addition, our results showed that depletion of NCAPG significantly inhibited NSCLC cell proliferation, migration, and self-renewal abilities, yet these could be reversed by adding microRNA (miRNA)-214-3p. Knockdown of long noncoding RNA (lncRNA) thymidylate synthetase opposite strand (TYMSOS) also inhibits the NSCLC cell proliferation, migration, and self-renewal abilities. In summary, our findings demonstrated that the crucial roles of the FOXM1/lncRNA-TYMSOS/miRNA-214-3p/NCAPG axis in NSCLC may shed light on how NCAPG may act as a therapeutic target for NSCLC.

Focus on Molecular Functions of Anti-Aging Deacetylase SIRT3

SIRT3 is a protein lysine deacetylase with a prominent role in the maintenance of mitochondrial integrity, which is a vulnerable target in many diseases. Intriguingly, cellular aging is reversible just by SIRT3 overexpression, which raises many questions about the role of SIRT3 in the molecular anti-aging mechanisms. Therefore, functions of SIRT3 were analyzed through the interaction network of 407 substrates collected by data mining. Results of the pathway enrichment and gene function prediction confirmed functions in the primary metabolism and mitochondrial ATP production. However, it also suggested involvement in thermogenesis, brain-related neurodegenerative diseases Alzheimer's (AD), Parkinson's, Huntington's disease (HD), and non-alcoholic fatty liver disease. The protein node prioritization analysis identified subunits of the complex I of the mitochondrial respiratory chain (MRC) as the nodes with the main regulatory effect within the entire interaction network. Additional high-ranked nodes were succinate dehydrogenase subunit B (SDHB), complex II, and ATP5F1, complex V of MRC. The analysis supports existence of the NADH/NAD+ driven regulatory feedback loop between SIRT3, complex I (MRC), and acetyl-CoA synthetases, and existence of the nuclear substrates of SIRT3. Unexplored functions of SIRT3 substrates such as LMNA and LMNB; HIF-1a, p53, DNA-PK, and PARK7 are highlighted for further scientific advances. SIRT3 acts as a repressor of BACE1 through the SIRT3-LKB1-AMPK-CREB-PGC1A-PPARG-BACE1 (SIRT3-BACE1), which functions are fitted the best by the Circadian Clock pathway. It forms a new working hypothesis as the therapeutical target for AD treatment. Other important pathways linked to SIRT3 activity are highlighted for therapeutical interventions.

Systematic analysis of the potential off-target activities of osimertinib by computational target fishing

Osimertinib is a third-generation epidermal growth factor receptor (EGFR) tyrosine kinase inhibitor used to treat non-small cell lung cancer. However, its off-targets are obscure, and systematic analysis of off-target activities remains to be performed. Here, we identified the off-targets of osimertinib using PharmMapper and DRAR-CPI and analyzed the intersected targets using the GeneMANIA and DAVID servers. A drug-target-pathway network was constructed to visualize the associations. The results showed that osimertinib is associated with 31 off-targets, 40 Kyoto Encyclopedia of Genes and Genomes pathways, and 9 diseases. Network analysis revealed that the targets were involved in cancer and other physiological processes. In addition to EGFR, molecular docking analysis showed that seven proteins, namely Janus kinase 3, peroxisome proliferator-activated receptor alpha, renin, mitogen-activated protein kinases, lymphocyte-specific protein tyrosine kinase, cell division protein kinase 2 and proto-oncogene tyrosine-protein kinase Src, could also be potential targets of osimertinib. In conclusion, osimertinib is predicted to target multiple proteins and pathways, resulting in the formation of an action network via which it exerts systematic pharmacological effects.

Multi-Omics Reveal the Immunological Role and the Theragnostic Value of miR-216a/GDF15 Axis in Human Colon Adenocarcinoma

Colon adenocarcinoma (COAD) is the most common type of gastrointestinal cancer and is still the third leading cause of cancer-related mortality worldwide. Accurate screening tools for early diagnosis and prediction of prognosis and precision treatment strategies are urgently required to accommodate the unmet medical needs of COAD management. We herein aimed to explore the significance of the microRNA (miR)-216a/growth differentiation factor 15 (GDF15) axis in terms of clinical value, tumor immunity, and potential mechanisms in COAD by using multi-omic analysis. The gene expression levels of miR-216a and GDF15 showed an increase in the COAD group compared to those of the normal group. The expression of miR-216a presented a negative correlation with GDF15 in COAD tumor tissue. The use of an in vitro luciferase reporter assay and bioinformatic prediction revealed that miR-216a-3p acted toward translational inhibition on GDF15 by targeting its 3′untranslated region (UTR) site. High miR-216a expression was associated with decreased overall survival (OS), while the high expression of GDF15 was associated with increased OS. Enriched type 1 T-helper (Th1), enriched regulatory T (Treg), enriched eosinophils, and decreased nature killer T-cells (NKTs) in COAD tumor tissue may play counteracting factors on the tumor-regulatory effects of miR-216a and GDF15. In addition, high GDF15 expression had associations with suppressed immunoinhibitory genes and negative correlations with the infiltration of macrophages and endothelial cells. The enrichment analysis revealed that GDF15 and its co-expression network may be implicated in mitochondrial organization, apoptosis signaling, and endoplasmic reticulum (ER) stress response. The Genomics of Drug Sensitivity in Cancer (GDSC) and Cancer Therapeutics Response Portal (CTRP) analysis identified that Gemcitabine acted as a precision treatment for COAD when GDF15 expression was low. This study supports the miR-216a/GDF15 axis as a diagnostic/prognostic panel for COAD, identifies Th1, Treg, eosinophils, and NKTs as counteracting factors, indicates potential relationships underlying immunomodulation, mitochondrial organization, apoptotic signaling, and ER stress and unveil Gemcitabine as a potential drug for the development of treatment strategy when combined with targeting GDF15.

Up-regulation of Arl4a gene expression by broccoli aqueous extract is associated with improved spermatogenesis in mouse testes

INTRODUCTION

Broccoli (Brassica oleracea) is well known for its properties as an anticancer, antioxidant, and scavenger of free radicals. However, its benefits in enhancing spermatogenesis have not been well established.

OBJECTIVE

To study broccoli aqueous extract effects on sperm factors and the expression of genes Catsper1, Catsper2, Arl4a, Sox5, and Sox9 in sperm factors in mice.

MATERIAL AND METHODS

Male mice were divided randomly into six groups: (1) Control; (2) cadmium (3 mg/kg of mouse body weight); (3) orally treated with 200 μl broccoli aqueous extract (1 g ml-1); (4) orally treated with 400 μl of broccoli aqueous extract; (5) orally treated with 200 broccoli aqueous extract plus cadmium, and (6) orally treated with 400 μl of broccoli aqueous extract plus cadmium. We analyzed the sperms factors and Catsper1, Catsper2, Arl4a, Sox5, and Sox9 gene expression.

RESULTS

An obvious improvement in sperm count and a slight enhancement in sperm motility were observed in mice treated with broccoli extract alone or with cadmium. Sperm viability was reduced by broccoli extract except for the 200 μl dose with cadmium, which significantly increased it. Interestingly, Arl4a gene expression increased in the 400 μl broccoli-treated group. Likewise, the Arl4a mRNA level in mice treated with cadmium and 200 μl of broccoli extract was higher than in the cadmium-treated mice. Furthermore, broccoli extract enhanced the mRNA level of Catsper2 and Sox5 genes in mice treated with 200 μl and 400 μl broccoli extract plus cadmium compared with the group treated solely with cadmium.

CONCLUSION

The higher sperm count in broccoli-treated mice opens the way for the development of pharmaceutical products for infertile men.

Single-Cell RNA Sequencing Defines the Regulation of Spermatogenesis by Sertoli-Cell Androgen Signaling

Androgen receptor (AR) signaling is essential for maintaining spermatogenesis and male fertility. However, the molecular mechanisms by which AR acts between male germ cells and somatic cells during spermatogenesis have not begun to be revealed until recently. With the advances obtained from the use of transgenic mice lacking AR in Sertoli cells (SCARKO) and single-cell transcriptomic sequencing (scRNA-seq), the cell specific targets of AR action as well as the genes and signaling pathways that are regulated by AR are being identified. In this study, we collected scRNA-seq data from wild-type (WT) and SCARKO mice testes at p20 and identified four somatic cell populations and two male germ cell populations. Further analysis identified that the distribution of Sertoli cells was completely different and uncovered the cellular heterogeneity and transcriptional changes between WT and SCARKO Sertoli cells. In addition, several differentially expressed genes (DEGs) in SCARKO Sertoli cells, many of which have been previously implicated in cell cycle, apoptosis and male infertility, have also been identified. Together, our research explores a novel perspective on the changes in the transcription level of various cell types between WT and SCARKO mice testes, providing new insights for the investigations of the molecular and cellular processes regulated by AR signaling in Sertoli cells.

Protein Receptors Evolved from Homologous Cohesion Modules That Self-Associated and Are Encoded by Interactive Networked Genes

Previously, it was proposed that protein receptors evolved from self-binding peptides that were encoded by self-interacting gene segments (inverted repeats) widely dispersed in the genome. In addition, self-association of the peptides was thought to be mediated by regions of amino acid sequence similarity. To extend these ideas, special features of receptors have been explored, such as their degree of homology to other proteins, and the arrangement of their genes for clues about their evolutionary origins and dynamics in the genome. As predicted, BLASTP searches for homologous proteins detected a greater number of unique hits for queries with receptor sequences than for sequences of randomly-selected, non-receptor proteins. This suggested that the building blocks (cohesion modules) for receptors were duplicated, dispersed, and maintained in the genome, due to structure/function relationships discussed here. Furthermore, the genes coding for a representative panel of receptors participated in a larger number of gene-gene interactions than for randomly-selected genes. This could conceivably reflect a greater evolutionary conservation of the receptor genes, with their more extensive integration into networks along with inherent properties of the genes themselves. In support of the latter possibility, some receptor genes were located in active areas of adaptive gene relocation/amalgamation to form functional blocks of related genes. It is suggested that adaptive relocation might allow for their joint regulation by common promoters and enhancers, and affect local chromatin structural domains to facilitate or repress gene expression. Speculation is included about the nature of the coordinated communication between receptors and the genes that encode them.

Systematic Genome-Wide Profiles Reveal Alternative Splicing Landscape and Implications of Splicing Regulator DExD-Box Helicase 21 in Aggressive Progression of Adrenocortical Carcinoma

Alternative splicing (AS) in the tumor biological process has provided a novel perspective on carcinogenesis. However, the clinical significance of individual AS patterns of adrenocortical carcinoma (ACC) has been underestimated, and in-depth investigations are lacking. We selected 76 ACC samples from the Cancer Genome Atlas (TCGA) SpliceSeq and SpliceAid2 databases, and 39 ACC samples from Fudan University Shanghai Cancer Center (FUSCC). Prognosis-related AS events (PASEs) and survival analysis were evaluated based on prediction models constructed by machine-learning algorithm. In total, 23,984 AS events and 3,614 PASEs were detected in the patients with ACC. The predicted risk score of each patient suggested that eight PASEs groups were significantly correlated with the clinical outcomes of these patients (p < 0.001). Prognostic models produced AUC values of 0.907 in all PASEs' groups. Eight splicing factors (SFs), including BAG2, CXorf56, DExD-Box Helicase 21 (DDX21), HSPB1, MBNL3, MSI1, RBMXL2, and SEC31B, were identified in regulatory networks of ACC. DDX21 was identified and validated as a novel clinical promoter and therapeutic target in 115 patients with ACC from TCGA and FUSCC cohorts. In conclusion, the strict standards used in this study ensured the systematic discovery of profiles of AS events using genome-wide cohorts. Our findings contribute to a comprehensive understanding of the landscape and underlying mechanism of AS, providing valuable insights into the potential usages of DDX21 for predicting prognosis for patients with ACC.

Supplementary Information

The online version contains supplementary material available at 10.1007/s43657-021-00026-x.

Prioritization of human well-being spectrum related GWAS-SNVs using ENCODE-based web-tools predict interplay between PSMC3, ITIH4, and SERPINC1 genes in modulating well-being

Several traits related to positive and negative affect show a high genetic as well as phenotypic correlation with well-being in humans, and are therefore collectively termed as "Well-being spectrum". Genome-Wide Association studies (GWA studies) on "well-being measurement" have led to identification of several genomic variants (Single Nucleotide Variants - SNVs), but very little has been explained with respect to their functionality and mode of alteration of well-being. Utilizing a pool of 1258 GWA studies based SNVs on "well-being measurement", we prioritized the SNVs and tried to annotate well-being related functionality through several bioinformatic tools to predict whether a protein sequence variation affects protein function, as well as experimentally validated datasets available in ENCODE based web-tools namely rSNPBase, RegulomeDB, Haploreg, along with GTEx Portal and STRING based protein interaction networks. Prioritization yielded three key SNVs; rs3781627-A, rs13072536-T and 5877-C potentially regulating three genes, PSMC3, ITIH4 and SERPINC1, respectively. Interestingly, the genes showed well clustered protein-protein interaction (maximum combined confidence score >0.4) with other well-being candidate genes, namely TNF and CRP genes suggesting their important role in modulation of well-being. PSMC3 and ITIH4 genes are also involved in driving acute phase responses signifying a probable cross-talk between well-being and psychoneuroimmunological system. To best of our knowledge this study is the first of its kind where the well-being associated GWA studies-SNVs were prioritized and functionally annotated, majorly based on functional data available in public domain, which revealed PSMC3, ITIH4 and SERPINC1 genes as probable candidates in regulation of well-being spectrum.

Multi-omics data integration and network-based analysis drives a multiplex drug repurposing approach to a shortlist of candidate drugs against COVID-19

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic is undeniably the most severe global health emergency since the 1918 Influenza outbreak. Depending on its evolutionary trajectory, the virus is expected to establish itself as an endemic infectious respiratory disease exhibiting seasonal flare-ups. Therefore, despite the unprecedented rally to reach a vaccine that can offer widespread immunization, it is equally important to reach effective prevention and treatment regimens for coronavirus disease 2019 (COVID-19). Contributing to this effort, we have curated and analyzed multi-source and multi-omics publicly available data from patients, cell lines and databases in order to fuel a multiplex computational drug repurposing approach. We devised a network-based integration of multi-omic data to prioritize the most important genes related to COVID-19 and subsequently re-rank the identified candidate drugs. Our approach resulted in a highly informed integrated drug shortlist by combining structural diversity filtering along with experts' curation and drug-target mapping on the depicted molecular pathways. In addition to the recently proposed drugs that are already generating promising results such as dexamethasone and remdesivir, our list includes inhibitors of Src tyrosine kinase (bosutinib, dasatinib, cytarabine and saracatinib), which appear to be involved in multiple COVID-19 pathophysiological mechanisms. In addition, we highlight specific immunomodulators and anti-inflammatory drugs like dactolisib and methotrexate and inhibitors of histone deacetylase like hydroquinone and vorinostat with potential beneficial effects in their mechanisms of action. Overall, this multiplex drug repurposing approach, developed and utilized herein specifically for SARS-CoV-2, can offer a rapid mapping and drug prioritization against any pathogen-related disease.

Targeting cathepsins: A potential link between COVID-19 and associated neurological manifestations

Many studies have shown that the lysosomal cathepsins, especially cathepsins B/L (CTSB/L) are required for SARS-CoV-2 entry into host cells. Lysosomal proteases, cathepsins are indispensable for normal health and are involved in several brain disorders occurring at different development age periods. On the other hand, it has been well known that COVID-19 infection is largely associated with several neurological disorders. Taken together these findings and given the high levels of expression of CTSB/L in the brain, we here proposed a reasonable hypothesis about the involvement of CTSB/L in the neurological manifestations linked to COVID-19. Pharmacological inhibitions of the CTSB/L could be a potential therapeutic target to block the virus entry as well as to mitigate the brain disorders. To this end, we utilized the network-based drug repurposing analyses to identify the possible drugs that can target CTSB/L. This study identifies the molecules like cyclosporine, phenytoin, and paclitaxel as potential drugs with binding ability to the CTSB/L. Further, we have performed molecular docking and all-atom molecular dynamics (MD) simulations to investigate the stability of CTSL-drug complexes. The results showed strong and stable binding of drugs with CTSL.

Multi-Omics Reveals the Immunological Role and Prognostic Potential of Mitochondrial Ubiquitin Ligase MARCH5 in Human Breast Cancer

E3 ubiquitin-linked enzyme MARCH5, also known as membrane-associated circular finger 5, is an enzyme encoded by the human MARCH5 gene. The main objective of this study was to visualize the prognosis of MARCH5 in breast cancer and to determine the relationship between MARCH5 expression and tumor immunity. MARCH5 expression was significantly higher in several cancers, including breast cancer (BRCA), compared with corresponding normal tissues. Not only was high MARCH5 expression associated with poorer overall survival, but also MARCH5 expression was positively correlated with the number of tumor-infiltrating immune cells in BRCA malignant tissues. Furthermore, MARCH5 expression showed a strong correlation with various immune markers of BRCA, suggesting its role in regulating tumor immunity. MARCH5 is a useful prognostic biomarker in several cancers, and its expression is highly correlated with tumor immune cell infiltration, and increased MARCH5 expression may serve as a new biomarker for BRCA diagnosis and prognosis.

Ist2 recruits the lipid transporters Osh6/7 to ER-PM contacts to maintain phospholipid metabolism

ER-plasma membrane (PM) contacts are proposed to be held together by distinct families of tethering proteins, which in yeast include the VAP homologues Scs2/22, the extended-synaptotagmin homologues Tcb1/2/3, and the TMEM16 homologue Ist2. It is unclear whether these tethers act redundantly or whether individual tethers have specific functions at contacts. Here, we show that Ist2 directly recruits the phosphatidylserine (PS) transport proteins and ORP family members Osh6 and Osh7 to ER-PM contacts through a binding site located in Ist2's disordered C-terminal tethering region. This interaction is required for phosphatidylethanolamine (PE) production by the PS decarboxylase Psd2, whereby PS transported from the ER to the PM by Osh6/7 is endocytosed to the site of Psd2 in endosomes/Golgi/vacuoles. This role for Ist2 and Osh6/7 in nonvesicular PS transport is specific, as other tethers/transport proteins do not compensate. Thus, we identify a molecular link between the ORP and TMEM16 families and a role for endocytosis of PS in PE synthesis.

Full-length isoform transcriptome of the developing human brain provides further insights into autism

Alternative splicing plays an important role in brain development, but its global contribution to human neurodevelopmental diseases (NDDs) requires further investigation. Here we examine the relationships between splicing isoform expression in the brain and de novo loss-of-function mutations from individuals with NDDs. We analyze the full-length isoform transcriptome of the developing human brain and observe differentially expressed isoforms and isoform co-expression modules undetectable by gene-level analyses. These isoforms are enriched in loss-of-function mutations and microexons, are co-expressed with a unique set of partners, and have higher prenatal expression. We experimentally test the effect of splice-site mutations and demonstrate exon skipping in five NDD risk genes, including SCN2A, DYRK1A, and BTRC. Our results suggest that the splice site mutation in BTRC reduces translational efficiency, likely affecting Wnt signaling through impaired degradation of β-catenin. We propose that functional effects of mutations should be investigated at the isoform- rather than gene-level resolution.

Integrated analysis of the expression, involved functions, and regulatory network of RUNX3 in melanoma

BACKGROUND

As a tumor suppressor or oncogenic gene, abnormal expression of RUNX family transcription factor 3 (RUNX3) has been reported in various cancers. <p> Introduction: This study aimed to investigate the role of RUNX3 in melanoma. <p> Methods: The expression level of RUNX3 in melanoma tissues was analyzed by immunohistochemistry and the Oncomine database. Based on microarray datasets GSE3189 and GSE7553, differentially expressed genes (DEGs) in melanoma samples were screened, followed by functional enrichment analysis. Gene Set Enrichment Analysis (GSEA) was performed for RUNX3. DEGs that co-expressed with RUNX3 were analyzed, and the transcription factors (TFs) of RUNX3 and its co-expressed genes were predicted. The protein-protein interactions (PPIs) for RUNX3 were analyzed utilizing the GeneMANIA database. MicroRNAs (miRNAs) that could target RUNX3 expression, were predicted. <p> Results: RUNX3 expression was significantly up-regulated in melanoma tissues. GSEA showed that RUNX3 expression was positively correlated with melanogenesis and melanoma pathways. Eleven DEGs showed significant co-expression with RUNX3 in melanoma, for example, TLE4 was negatively co-expressed with RUNX3. RUNX3 was identified as a TF that regulated the expression of both itself and its co-expressed genes. PPI analysis showed that 20 protein-encoding genes interacted with RUNX3, among which 9 genes were differentially expressed in melanoma, such as CBFB and SMAD3. These genes were significantly enriched in transcriptional regulation by RUNX3, RUNX3 regulates BCL2L11 (BIM) transcription, regulation of I-kappaB kinase/NF-kappaB signaling, and signaling by NOTCH. A total of 31 miRNAs could target RUNX3, such as miR-326, miR-330-5p, and miR-373-3p. <p> Conclusion: RUNX3 expression was up-regulated in melanoma and was implicated in the development of melanoma.