Network-based approach to identify key candidate genes and pathways shared by thyroid cancer and chronic kidney disease

Identification of Biomarkers and Molecular Pathways Implicated in Smoking and COVID-19 Associated Lung Cancer Using Bioinformatics and Machine Learning Approaches

Lung cancer (LC) is a significant global health issue, with smoking as the most common cause. Recent epidemiological studies have suggested that individuals who smoke are more susceptible to COVID-19. In this study, we aimed to investigate the influence of smoking and COVID-19 on LC using bioinformatics and machine learning approaches. We compared the differentially expressed genes (DEGs) between LC, smoking, and COVID-19 datasets and identified 26 down-regulated and 37 up-regulated genes shared between LC and smoking, and 7 down-regulated and 6 up-regulated genes shared between LC and COVID-19. Integration of these datasets resulted in the identification of ten hub genes (SLC22A18, CHAC1, ROBO4, TEK, NOTCH4, CD24, CD34, SOX2, PITX2, and GMDS) from protein-protein interaction network analysis. The WGCNA R package was used to construct correlation network analyses for these shared genes, aiming to investigate the relationships among them. Furthermore, we also examined the correlation of these genes with patient outcomes through survival curve analyses. The gene ontology and pathway analyses were performed to find out the potential therapeutic targets for LC in smoking and COVID-19 patients. Moreover, machine learning algorithms were applied to the TCGA RNAseq data of LC to assess the performance of these common genes and ten hub genes, demonstrating high performances. The identified hub genes and molecular pathways can be utilized for the development of potential therapeutic targets for smoking and COVID-19-associated LC.

TidyGEO: preparing analysis-ready datasets from Gene Expression Omnibus

TidyGEO is a Web-based tool for downloading, tidying, and reformatting data series from Gene Expression Omnibus (GEO). As a freely accessible repository with data from over 6 million biological samples across more than 4000 organisms, GEO provides diverse opportunities for secondary research. Although scientists may find assay data relevant to a given research question, most analyses require sample-level annotations. In GEO, such annotations are stored alongside assay data in delimited, text-based files. However, the structure and semantics of the annotations vary widely from one series to another, and many annotations are not useful for analysis purposes. Thus, every GEO series must be tidied before it is analyzed. Manual approaches may be used, but these are error prone and take time away from other research tasks. Custom computer scripts can be written, but many scientists lack the computational expertise to create such scripts. To address these challenges, we created TidyGEO, which supports essential data-cleaning tasks for sample-level annotations, such as selecting informative columns, renaming columns, splitting or merging columns, standardizing data values, and filtering samples. Additionally, users can integrate annotations with assay data, restructure assay data, and generate code that enables others to reproduce these steps.

Small RNA sequencing evaluation of renal microRNA biomarkers in dogs with X-linked hereditary nephropathy

Dogs with X-linked hereditary nephropathy (XLHN) are an animal model for Alport syndrome in humans and progressive chronic kidney disease (CKD). Using mRNA sequencing (mRNA-seq), we have characterized the gene expression profile affecting the progression of XLHN; however, the microRNA (miRNA, miR) expression remains unknown. With small RNA-seq and quantitative RT-PCR (qRT-PCR), we used 3 small RNA-seq analysis tools (QIAGEN OmicSoft Studio, miRDeep2, and CPSS 2.0) to profile differentially expressed renal miRNAs, top-ranked miRNA target genes, and enriched biological processes and pathways in CKD progression. Twenty-three kidney biopsies were collected from 5 dogs with XLHN and 4 age-matched, unaffected littermates at 3 clinical time points (T1: onset of proteinuria, T2: onset of azotemia, and T3: advanced azotemia). We identified up to 23 differentially expressed miRNAs at each clinical time point. Five miRNAs (miR-21, miR-146b, miR-802, miR-142, miR-147) were consistently upregulated in affected dogs. We identified miR-186 and miR-26b as effective reference miRNAs for qRT-PCR. This study applied small RNA-seq to identify differentially expressed miRNAs that might regulate critical pathways contributing to CKD progression in dogs with XLHN.

Identification of hub genes in thyroid carcinoma to predict prognosis by integrated bioinformatics analysis

The aim of this study was to identify hub genes closely related to the pathogenesis and prognosis of thyroid carcinoma (THCA) by integrated bioinformatics analysis. In this study, through differential gene expression analysis, 1916 and 665 differentially expressed genes (DEGs) were obtained from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) database, and 7 and 11 co-expressed modules were identified from the TCGA-THCA and GSE153659 datasets, respectively, by weighted gene co-expression network analysis. We identified 162 overlapping genes between the DEGs and co-expression module genes as candidate hub genes. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses of the 162 overlapping DEGs identified significant functions and pathways of THCA, such as thyroid hormone generation and metabolic process. A protein-protein interaction (PPI) analysis detected the top 10 hub genes (QSOX1, WFS1, EVA1A, FSTL3, CHRDL1, FABP4, PRDM16, PPARGC1A, PPARG, COL23A1). Finally, survival analysis, clinical correlation analysis, and protein abundance validation confirmed that 3 of the 10 hub genes were associated with survival prognosis of patients with THCA, and 8 of them were associated with the clinical stages of THCA. In summary, we identified hub genes and key modules that were closely related to THCA, and validated these genes by survival analysis, clinical correlation analysis, and Human Protein Atlas image analysis. Our results provide important information that will help to elucidate the pathogenesis of THCA and identify novel candidate prognostic biomarkers and potential therapeutic targets.

Clinically Applicable Machine Learning Approaches to Identify Attributes of Chronic Kidney Disease (CKD) for Use in Low-Cost Diagnostic Screening

OBJECTIVE

Chronic kidney disease (CKD) is a major public health concern worldwide. High costs of late-stage diagnosis and insufficient testing facilities can contribute to high morbidity and mortality rates in CKD patients, particularly in less developed countries. Thus, early diagnosis aided by vital parameter analytics using affordable computer-aided diagnosis could not only reduce diagnosis costs but improve patient management and outcomes.

METHODS

In this study, we developed machine learning models using selective key pathological categories to identify clinical test attributes that will aid in accurate early diagnosis of CKD. Such an approach will save time and costs for diagnostic screening. We have also evaluated the performance of several classifiers with k-fold cross-validation on optimized datasets derived using these selected clinical test attributes.

RESULTS

Our results suggest that the optimized datasets with important attributes perform well in diagnosis of CKD using our proposed machine learning models. Furthermore, we evaluated clinical test attributes based on urine and blood tests along with clinical parameters that have low costs of acquisition. The predictive models with the optimized and pathologically categorized attributes set yielded high levels of CKD diagnosis accuracy with random forest (RF) classifier being the best performing.

CONCLUSIONS

Our machine learning approach has yielded effective predictive analytics for CKD screening which can be developed as a resource to facilitate improved CKD screening for enhanced and timely treatment plans.

Pathogenetic profiling of COVID-19 and SARS-like viruses

The novel coronavirus (2019-nCoV) has recently emerged, causing COVID-19 outbreaks and significant societal/global disruption. Importantly, COVID-19 infection resembles SARS-like complications. However, the lack of knowledge about the underlying genetic mechanisms of COVID-19 warrants the development of prospective control measures. In this study, we employed whole-genome alignment and digital DNA-DNA hybridization analyses to assess genomic linkage between 2019-nCoV and other coronaviruses. To understand the pathogenetic behavior of 2019-nCoV, we compared gene expression datasets of viral infections closest to 2019-nCoV with four COVID-19 clinical presentations followed by functional enrichment of shared dysregulated genes. Potential chemical antagonists were also identified using protein-chemical interaction analysis. Based on phylogram analysis, the 2019-nCoV was found genetically closest to SARS-CoVs. In addition, we identified 562 upregulated and 738 downregulated genes (adj. P ≤ 0.05) with SARS-CoV infection. Among the dysregulated genes, SARS-CoV shared ≤19 upregulated and ≤22 downregulated genes with each of different COVID-19 complications. Notably, upregulation of BCL6 and PFKFB3 genes was common to SARS-CoV, pneumonia and severe acute respiratory syndrome, while they shared CRIP2, NSG1 and TNFRSF21 genes in downregulation. Besides, 14 genes were common to different SARS-CoV comorbidities that might influence COVID-19 disease. We also observed similarities in pathways that can lead to COVID-19 and SARS-CoV diseases. Finally, protein-chemical interactions suggest cyclosporine, resveratrol and quercetin as promising drug candidates against COVID-19 as well as other SARS-like viral infections. The pathogenetic analyses, along with identified biomarkers, signaling pathways and chemical antagonists, could prove useful for novel drug development in the fight against the current global 2019-nCoV pandemic.

Network-based computational approach to identify genetic links between cardiomyopathy and its risk factors

Cardiomyopathy (CMP) is a group of myocardial diseases that progressively impair cardiac function. The mechanisms underlying CMP development are poorly understood, but lifestyle factors are clearly implicated as risk factors. This study aimed to identify molecular biomarkers involved in inflammatory CMP development and progression using a systems biology approach. The authors analysed microarray gene expression datasets from CMP and tissues affected by risk factors including smoking, ageing factors, high body fat, clinical depression status, insulin resistance, high dietary red meat intake, chronic alcohol consumption, obesity, high-calorie diet and high-fat diet. The authors identified differentially expressed genes (DEGs) from each dataset and compared those from CMP and risk factor datasets to identify common DEGs. Gene set enrichment analyses identified metabolic and signalling pathways, including MAPK, RAS signalling and cardiomyopathy pathways. Protein-protein interaction (PPI) network analysis identified protein subnetworks and ten hub proteins (CDK2, ATM, CDT1, NCOR2, HIST1H4A, HIST1H4B, HIST1H4C, HIST1H4D, HIST1H4E and HIST1H4L). Five transcription factors (FOXC1, GATA2, FOXL1, YY1, CREB1) and five miRNAs were also identified in CMP. Thus the authors' approach reveals candidate biomarkers that may enhance understanding of mechanisms underlying CMP and their link to risk factors. Such biomarkers may also be useful to develop new therapeutics for CMP.

Network-Based Genetic Profiling Reveals Cellular Pathway Differences Between Follicular Thyroid Carcinoma and Follicular Thyroid Adenoma

Molecular mechanisms underlying the pathogenesis and progression of malignant thyroid cancers, such as follicular thyroid carcinomas (FTCs), and how these differ from benign thyroid lesions, are poorly understood. In this study, we employed network-based integrative analyses of FTC and benign follicular thyroid adenoma (FTA) lesion transcriptomes to identify key genes and pathways that differ between them. We first analysed a microarray gene expression dataset (Gene Expression Omnibus GSE82208, n = 52) obtained from FTC and FTA tissues to identify differentially expressed genes (DEGs). Pathway analyses of these DEGs were then performed using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) resources to identify potentially important pathways, and protein-protein interactions (PPIs) were examined to identify pathway hub genes. Our data analysis identified 598 DEGs, 133 genes with higher and 465 genes with lower expression in FTCs. We identified four significant pathways (one carbon pool by folate, p53 signalling, progesterone-mediated oocyte maturation signalling, and cell cycle pathways) connected to DEGs with high FTC expression; eight pathways were connected to DEGs with lower relative FTC expression. Ten GO groups were significantly connected with FTC-high expression DEGs and 80 with low-FTC expression DEGs. PPI analysis then identified 12 potential hub genes based on degree and betweenness centrality; namely, TOP2A, JUN, EGFR, CDK1, FOS, CDKN3, EZH2, TYMS, PBK, CDH1, UBE2C, and CCNB2. Moreover, transcription factors (TFs) were identified that may underlie gene expression differences observed between FTC and FTA, including FOXC1, GATA2, YY1, FOXL1, E2F1, NFIC, SRF, TFAP2A, HINFP, and CREB1. We also identified microRNA (miRNAs) that may also affect transcript levels of DEGs; these included hsa-mir-335-5p, -26b-5p, -124-3p, -16-5p, -192-5p, -1-3p, -17-5p, -92a-3p, -215-5p, and -20a-5p. Thus, our study identified DEGs, molecular pathways, TFs, and miRNAs that reflect molecular mechanisms that differ between FTC and benign FTA. Given the general similarities of these lesions and common tissue origin, some of these differences may reflect malignant progression potential, and include useful candidate biomarkers for FTC and identifying factors important for FTC pathogenesis.

Machine Learning and Bioinformatics Models to Identify Pathways that Mediate Influences of Welding Fumes on Cancer Progression

Welding generates and releases fumes that are hazardous to human health. Welding fumes (WFs) are a complex mix of metallic oxides, fluorides and silicates that can cause or exacerbate health problems in exposed individuals. In particular, WF inhalation over an extended period carries an increased risk of cancer, but how WFs may influence cancer behaviour or growth is unclear. To address this issue we employed a quantitative analytical framework to identify the gene expression effects of WFs that may affect the subsequent behaviour of the cancers. We examined datasets of transcript analyses made using microarray studies of WF-exposed tissues and of cancers, including datasets from colorectal cancer (CC), prostate cancer (PC), lung cancer (LC) and gastric cancer (GC). We constructed gene-disease association networks, identified signaling and ontological pathways, clustered protein-protein interaction network using multilayer network topology, and analyzed survival function of the significant genes using Cox proportional hazards (Cox PH) model and product-limit (PL) estimator. We observed that WF exposure causes altered expression of many genes (36, 13, 25 and 17 respectively) whose expression are also altered in CC, PC, LC and GC. Gene-disease association networks, signaling and ontological pathways, protein-protein interaction network, and survival functions of the significant genes suggest ways that WFs may influence the progression of CC, PC, LC and GC. This quantitative analytical framework has identified potentially novel mechanisms by which tissue WF exposure may lead to gene expression changes in tissue gene expression that affect cancer behaviour and, thus, cancer progression, growth or establishment.