Search

Abstract Objective Systemic chemotherapy is the first-line therapeutic option for head and neck squamous cell carcinoma (HNSCC), but it often fails. This study aimed to develop an effective prognostic model for evaluating the therapeutic effects of systemic chemotherapy. Methods This study utilized CRISPR/cas9 whole gene loss-of-function library screening and data from The Cancer Genome Atlas (TCGA) HNSCC patients who have undergone systemic therapy to examine differentially expressed genes (DEGs). A lipid metabolism-related clustered polygenic model called the lipid metabolism related score (LMRS) model was established based on the identified functionally enriched DEGs. The prediction efficiency of the model for survival outcome, chemotherapy, and immunotherapy response was evaluated using HNSCC datasets, the GEO database and clinical samples. Results Screening results from the study demonstrated that genes those were differentially expressed were highly associated with lipid metabolism-related pathways, and patients receiving systemic therapy had significantly different prognoses based on lipid metabolism gene characteristics. The LMRS model, consisting of eight lipid metabolism-related genes, outperformed each lipid metabolism gene-based model in predicting outcome and drug response. Further validation of the LMRS model in HNSCCs confirmed its prognostic value. Conclusion In conclusion, the LMRS polygenic prognostic model is helpful to assess outcome and drug response for HNSCCs and could assist in the timely selection of the appropriate treatment for HNSCC patients. This study provides important insights for improving systemic chemotherapy and enhancing patient outcomes.

Cholesterol is the most abundant sterol molecule in mammalian cells, which not only constitutes the cell membrane but also plays essential roles in the synthesis of important hormones, synapse formation, and cell signal transduction. The effect of hypercholesterolemia on hearing has been studied extensively, and multiple studies have demonstrated that hypercholesterolemia is a risk factor for hearing loss. However, the impact of cholesterol homeostasis within auditory cells on peripheral auditory development and maintenance has not been evaluated in detail. Mutations in certain cholesterol metabolism-related genes, such as NPC1, SERAC1, DHCR7, and OSBPL2, as well as derivatives of cholesterol metabolism-related ototoxic drugs, such as β-cyclodextrin, can lead to disruptions of cholesterol homeostasis within auditory cells, resulting in hearing loss. This article aims to review the impact of cholesterol homeostasis within auditory cells on the peripheral auditory function from the following two perspectives: (1) changes in cholesterol homeostasis regulatory genes in various hearing loss models; (2) mechanisms underlying the effects of some drugs that have a therapeutic effect on hearing loss via regulating cholesterol homeostasis. This article aims to summarize and analyze the impact of disruption of cellular cholesterol homeostasis within auditory cells on hearing, in order to provide evidence regarding the underlying mechanisms.

Recent studies showed that lipid metabolism reprogramming contributes to tumorigenicity and malignancy by interfering energy production, membrane formation, and signal transduction in cancers. HNSCCs are highly reliant on aerobic glycolysis and glutamine metabolism. However, the mechanisms underlying lipid metabolism reprogramming in HNSCCs remains obscure. The present review summarizes and discusses the “vital” cellular signaling roles of the lipid metabolism reprogramming in HNSCCs. We also address the differences between HNSCCs regions caused by anatomical heterogeneity. We enumerate these recent findings into our current understanding of lipid metabolism reprogramming in HNSCCs and introduce the new and exciting therapeutic implications of targeting the lipid metabolism.

With the increase in life expectancy in the global population, aging societies have emerged in many countries, including China. As a common sensory defect in the elderly population, the prevalence of age-related hearing loss and its influence on society are increasing yearly. Metabolic syndrome is currently one of the main health problems in the world. Many studies have demonstrated that metabolic syndrome and its components are correlated with a variety of age-related diseases of the peripheral sensory system, including age-related hearing loss. Both age-related hearing loss and metabolic syndrome are high-prevalence chronic diseases, and many people suffer from both at the same time. In recent years, more and more studies have found that mitochondrial dysfunction occurs in both metabolic syndrome and age-related hearing loss. Therefore, to better understand the impact of metabolic syndrome on age-related hearing loss from the perspective of mitochondrial dysfunction, we reviewed the literature related to the relationship between age-related hearing loss and metabolic syndrome and their components to discern the possible role of mitochondria in both conditions.

The transcriptomic landscape of mice with primary auditory neurons degeneration (PAND) indicates key pathways in its pathogenesis, including complement cascades, immune responses, tumor necrosis factor (TNF) signaling pathway, and cytokine-cytokine receptor interaction. Toll-like receptors (TLRs) are important immune and inflammatory molecules that have been shown to disrupt the disease network of PAND. In a PAND model involving administration of kanamycin combined with furosemide to destroy cochlear hair cells, Tlr 2/4 double knockout (DKO) mice had auditory preservation advantages, which were mainly manifested at 4–16 kHz. DKO mice and wild type (WT) mice had completely damaged cochlear hair cells on the 30th day, but the density of spiral ganglion neurons (SGN) in the Rosenthal canal was significantly higher in the DKO group than in the WT group. The results of immunohistochemistry for p38 and p65 showed that the attenuation of SGN degeneration in DKO mice may not be mediated by canonical Tlr signaling pathways. The SGN transcriptome of DKO and WT mice indicated that there was an inverted gene set enrichment relationship between their different transcriptomes and the SGN degeneration transcriptome, which is consistent with the morphology results. Core module analysis suggested that DKO mice may modulate SGN degeneration by activating two clusters, and the involved molecules include EGF, STAT3, CALB2, LOX, SNAP25, CAV2, SDC4, MYL1, NCS1, PVALB, TPM4, and TMOD4.

Introduction: Atrial fibrillation (AF) is the most prevalent cardiac arrhythmia, and it significantly increases the risk of cardiovascular complications and morbidity, even with appropriate treatment. Tissue remodeling has been a significant topic, while its systematic transcriptional signature remains unclear in AF. Objectives: Our study aims to systematically investigate the molecular characteristics of AF at the cellular-level. Methods: We conducted single-nuclei RNA-sequencig (snRNA-seq) analysis using nuclei isolated from the left atrial appendage (LAA) of AF patients and sinus rhythm. Pathological staining was performed to validate the key findings of snRNA-seq. Results: A total of 30 cell subtypes were identified among 80, 592 nuclei. Within the LAA of AF, we observed a specific subtype of dedifferentiated cardiomyocytes (CMs) characterized by reduced expression of cardiac contractile proteins (TTN and TRDN) and heightened expression of extracellular-matrix related genes (COL1A2 and FBN1). Transcription factor prediction analysis revealed that gene expression patterns in dedifferentiated CMs were primarily regulated by CEBPG and GISLI. Additionally, we identified a distinct subtype of endothelial progenitor cells (EPCs) demonstrating elevated expression of PROM1 and KDR, a population decreased within the LAA of AF. Epicardial adipocytes disclosed a reduced release of the anti-inflammatory and anti-fibrotic factor PRG4, and an augmented secretion of VEGF signals targeting CMs. Additionally, we noted accumulation of M2-like macrophages and CD8+ T cells with high pro-inflammatory score in LAA of AF. Furthermore, the analysis of intercellular communication revealed specific pathways related to AF, such as inflammation, extracellular matrix, and vascular remodeling signals. Conclusions: This study has discovered the presence of dedifferentiated CMs, a decrease in endothelial progenitor cells, a shift in the secretion profile of adipocytes, and an amplified inflammatory response in AF. These findings could offer crucial insights for future research on AF and serve as valuable references for investigating novel therapeutic approaches for AF.

Abstract Background Osteoporotic vertebral compression fractures (OVCF) in the elderly increase refracture risk post-surgery, leading to higher mortality rates. Genome-wide association studies (GWAS) have identified susceptibility genes for osteoporosis, but the phenotypic variance explained by these genes has been limited, indicating the need to explore additional causal factors. Epigenetic modifications, such as DNA methylation, may influence osteoporosis and refracture risk. However, prospective cohorts for assessing epigenetic alterations in Chinese elderly patients are lacking. Here, we propose to conduct a prospective cohort study to investigate the causal network of DNA polymorphisms, DNA methylation, and environmental factors on the development of osteoporosis and the risk of refracture. Methods We will collect vertebral and peripheral blood from 500 elderly OVCF patients undergoing surgery, extract DNA, and generate whole genome genotype data and DNA methylation data. Observation indicators will be collected and combined with one-year follow-up data. A healthy control group will be selected from a natural population cohort. Epigenome-wide association studies (EWAS) of osteoporosis and bone mineral density will be conducted. Differential methylation analysis will compare candidate gene methylation patterns in patients with and without refracture. Multi-omics prediction models using genetic variants and DNA methylation sites will be built to predict OVCF risk. Discussion This study will be the first large-scale population-based study of osteoporosis and bone mineral density phenotypes based on genome-wide data, multi-time point methylation data, and phenotype data. By analyzing methylation changes related to osteoporosis and bone mineral density in OVCF patients, the study will explore the feasibility of DNA methylation in evaluating postoperative osteoporosis intervention effects. The findings may identify new molecular markers for effective anti-osteoporosis treatment and inform individualized prevention and treatment strategies. Trial registration chictr.org.cn ChiCTR2200065316, 02/11/2022.

Summary: Orofacial clefts (OFCs) are among the most common human congenital birth defects. Previous multiethnic studies have identified dozens of associated loci for both cleft lip with or without cleft palate (CL/P) and cleft palate alone (CP). Although several nearby genes have been highlighted, the “casual” variants are largely unknown. Here, we developed DeepFace, a convolutional neural network model, to assess the functional impact of variants by SNP activity difference (SAD) scores. The DeepFace model is trained with 204 epigenomic assays from crucial human embryonic craniofacial developmental stages of post-conception week (pcw) 4 to pcw 10. The Pearson correlation coefficient between the predicted and actual values for 12 epigenetic features achieved a median range of 0.50–0.83. Specifically, our model revealed that SNPs significantly associated with OFCs tended to exhibit higher SAD scores across various variant categories compared to less related groups, indicating a context-specific impact of OFC-related SNPs. Notably, we identified six SNPs with a significant linear relationship to SAD scores throughout developmental progression, suggesting that these SNPs could play a temporal regulatory role. Furthermore, our cell-type specificity analysis pinpointed the trophoblast cell as having the highest enrichment of risk signals associated with OFCs. Overall, DeepFace can harness distal regulatory signals from extensive epigenomic assays, offering new perspectives for prioritizing OFC variants using contextualized functional genomic features. We expect DeepFace to be instrumental in accessing and predicting the regulatory roles of variants associated with OFCs, and the model can be extended to study other complex diseases or traits.

Summary: Orofacial clefts (OFCs) are common congenital birth defects with various etiologies, including genetic variants. Online Mendelian Inheritance in Man (OMIM) annotated several hundred genes involving OFCs. Furthermore, several hundreds of de novo variants (DNVs) have been identified from individuals with OFCs. Some DNVs are related to known OFC genes or pathways, but there are still many DNVs whose relevance to OFC development is unknown. To explore novel gene functions and their cellular expression profiles, we focused on DNVs in genes that were not listed in OMIM. We collected 960 DNVs in 853 genes from published studies and curated these genes, based on the DNVs’ deleteriousness, into 230 and 23 genes related to cleft lip with or without cleft palate (CL/P) and cleft palate only (CPO), respectively. For comparison, we curated 178 CL/P and 277 CPO genes from OMIM. In CL/P, the pathways enriched in DNV and OMIM genes were significantly overlapped (p = 0.002). Single-cell RNA sequencing (scRNA-seq) analysis of mouse lip development revealed that both gene sets had abundant expression in the ectoderm (DNV genes: adjusted p = 0.032, OMIM genes: adjusted p

The incidence of bladder cancer and patient survival vary greatly among different populations, but the influence of the associated molecular features and evolutionary processes on its clinical treatment and prognostication remains unknown. Here, we analyze the genomic architectures of 505 bladder cancer patients from Asian/Black/White populations. We identify a previously unknown association between AHNAK mutations and activity of the APOBEC‐a mutational signature, the activity of which varied substantially across populations. All significantly mutated genes but only half of arm‐level somatic copy number alterations (SCNAs) are enriched with clonal events, indicating large‐scale SCNAs as rich sources of bladder cancer clonal diversities. The prevalence of TP53 and ATM clonal mutations as well as the associated burden of SCNAs is significantly higher in Whites/Blacks than in Asians. We identify a trans‐ancestry prognostic subtype of bladder cancer characterized by enrichment of non‐muscle‐invasive patients and muscle‐invasive patients with good prognosis, increased CREBBP/FGFR3/HRAS/NFE2L2 mutations, decreased intra‐tumor heterogeneity and genome instability, and an activated tumor microenvironment.

The dysregulation of the m6A epitranscriptomic networks have been reported to contribute to the development of gliomas. Here, the authors utilize induced pluripotent stem cell-derived astrocytes with a p53 mutation and demonstrate that mutant p53 upregulates the m6A reader YTHDF2, resulting in the initiation of gliomas.

With the rapid progress of sensor technology, underwater monitoring wireless sensor network (UMWSN) equipped with water resource monitoring capabilities have surfaced. UMWSN is currently emerging as a popular research and application area, with a wide range of applications in various fields, including but not limited to disaster warning, pollutant monitoring and hydrological data monitoring. However, due to the limited size and computing power of wireless sensors, UMWSN are exceedingly susceptible to network security attacks by contaminating entities. In order to tackle the problems of inadequate security, elevated network energy consumption, and high system latency in existing QoS routing research, this paper puts forward an innovative multi-objective routing protocol named multi-objective chaotic Boltzmann ant colony optimization routing protocol (MO-CBACORP). MO-CBACORP can minimize the energy consumption of sensor nodes in UMWSN routing, thereby effectively diminishing system latency and significantly enhancing routing security. During the routing optimization process, MO-CBACORP proposed an innovative design of the chaos operator and Boltzmann selection strategy, which effectively avoids the issue of being stuck in local optimal solutions and greatly enhances the convergence speed of the algorithm. Additionally, a novel experiment for underwater monitoring environments is presented that precisely replicates underwater monitoring conditions and attack scenarios. According to the experimental results, MO-CBACORP exhibits superior performance compared to the latest state-of-the-art research including THOA-MHR, ACO-MCMC, and MAP-ACO, with a minimum improvement of 7.25%, 9.31%, and 7.68% in terms of route security, network energy consumption, and delay, respectively. All in all, MO-CBACORP exhibits superior performance compared to the three latest research algorithms.

Abstract Background The rapid accumulation of single-cell RNA sequencing (scRNA-seq) data presents unique opportunities to decode the genetically mediated cell-type specificity in complex diseases. Here, we develop a new method, scGWAS, which effectively leverages scRNA-seq data to achieve two goals: (1) to infer the cell types in which the disease-associated genes manifest and (2) to construct cellular modules which imply disease-specific activation of different processes. Results scGWAS only utilizes the average gene expression for each cell type followed by virtual search processes to construct the null distributions of module scores, making it scalable to large scRNA-seq datasets. We demonstrated scGWAS in 40 genome-wide association studies (GWAS) datasets (average sample size N ≈ 154,000) using 18 scRNA-seq datasets from nine major human/mouse tissues (totaling 1.08 million cells) and identified 2533 trait and cell-type associations, each with significant modules for further investigation. The module genes were validated using disease or clinically annotated references from ClinVar, OMIM, and pLI variants. Conclusions We showed that the trait-cell type associations identified by scGWAS, while generally constrained to trait-tissue associations, could recapitulate many well-studied relationships and also reveal novel relationships, providing insights into the unsolved trait-tissue associations. Moreover, in each specific cell type, the associations with different traits were often mediated by different sets of risk genes, implying disease-specific activation of driving processes. In summary, scGWAS is a powerful tool for exploring the genetic basis of complex diseases at the cell type level using single-cell expression data.

Summary: Ischemia reperfusion injury (IRI), often related to surgical procedures, is one of the important causes of acute kidney injury (AKI). To decipher the dynamic process of AKI caused by IRI (with prolonged ischemia phase), we performed single-cell RNA sequencing (scRNA-seq) of clinically relevant IRI murine model with different ischemic intervals. We discovered that Slc5a2hi proximal tubular cells were susceptible to AKI and highly expressed neutral amino acid transporter gene Slc6a19, which was dramatically decreased over the time course. With the usage of mass spectrometry–based metabolomic analysis, we detected that the level of neutral amino acid isoleucine dropped off in AKI mouse plasma metabolites. And the reduction of plasma isoleucine was also verified in patients with cardiac surgery–associated acute kidney injury (CSA-AKI). The findings advanced the understanding of dynamic process of AKI and introduced reduction of isoleucine as a potential biomarker for CSA-AKI.

The SARS-CoV-2 virus rapidly spread worldwide, threatening public health. Since it emerged, the scientific community has been engaged in the development of effective therapeutics and vaccines. The subunit S1 in the spike protein of SARS-CoV-2 mediates the viral entry into the host and is therefore one of the major research targets. The S1 protein is extensively glycosylated, and there is compelling evidence that glycans protect the virus’ active site from the human defense system. Therefore, investigation of the S1 protein glycome alterations in the different virus variants will provide a view of the glycan evolution and its relationship with the virus pathogenesis. In this study, we explored the N-glycosylation expression of the S1 protein for eleven SARS-CoV-2 variants: five variants of concern (VOC), including alpha, beta, gamma, delta, and omicron, and six variants of interest (VOI), including epsilon, eta, iota, lambda, kappa, and mu. The results showed significant differences in the N-glycome abundance of all variants. The N-glycome of the VOC showed a large increase in the abundance of sialofucosylated glycans, with the greatest abundance in the omicron variant. In contrast, the results showed a large abundance of fucosylated glycans for most of the VOI. Two glycan compositions, GlcNAc4,Hex5,Fuc,NeuAc (4-5-1-1) and GlcNAc6,Hex8,Fuc,NeuAc (6-8-1-1), were the most abundant structures across all variants. We believe that our data will contribute to understanding the S1 protein’s structural differences between SARS-CoV-2 mutations.

Monitoring and identifying the freshness levels of meat holds significant importance in the field of food safety as it directly relates to human dietary safety. Traditional packaging methods for lamb meat quality assessment present issues such as cumbersome operations and irreversible damage. This research proposes a quality assessment method for modified atmosphere packaging lamb meat using near-infrared spectroscopy and multi-parameter fusion. Fresh lamb meat quality is taken as the research subject, comparing various physicochemical indicators and near-infrared spectroscopic information under different temperatures (4 °C and 10 °C) and different modified atmosphere packaging combinations. Through precision parameter comparison, rebound and TVB-N values are selected as the modeling parameters. Six spectral preprocessing methods (multi-scatter calibration, MSC; standard normal variate transformation, SNV; normalization; Savitzky–Golay smoothing, SG; Savitzky–Golay 1 derivative, SG-1st; and Savitzky–Golay 2 derivative, SG-2nd), and three feature wavelength selection methods (competitive adaptive reweighted sampling, CARS; successive projections algorithm, SPA; and uninformative variable elimination, UVE) are compared. Partial least squares (PLS) and support vector machine (SVM) are used to construct prediction models for chilled fresh lamb meat quality. The results show that when rebound is used as a parameter, the SG-2nd-SPA-PLSR model has the highest accuracy, with a determination coefficient R2p of 0.94 for the prediction set. When TVB-N is used as a parameter, the MSC-UVE-SVM model has the highest accuracy, with an R2p of 0.95 for the prediction set. In conclusion, the use of near-infrared spectroscopic analysis enables rapid and non-destructive prediction and evaluation of lamb meat freshness, including its textural characteristics and TVB-N content under different modified atmosphere packaging. This study provides a theoretical basis and technical support for further encapsulating the models into portable devices and developing portable near-infrared spectrometers to rapidly determine lamb meat freshness.

Drug response in cancer patients vary dramatically due to inter- and intra-tumor heterogeneity and transcriptome context plays a significant role in shaping the actual treatment outcome. Here, the authors develop a deep variational autoencoder model to compress gene signatures into latent vectors and accurately impute drug response.

Testicular germ cell tumours (TGCTs) are the most common cancer in young male adults (aged 15 to 40). Unlike most other cancer types, identification of molecular signatures in TGCT has rarely reported. In this study, we developed a novel integrative analysis framework to identify co-methylated and co-expressed genes [mRNAs and microRNAs (miRNAs)] modules in two TGCT subtypes: non-seminoma (NSE) and seminoma (SE). We first integrated DNA methylation and mRNA/miRNA expression data and then used a statistical method, CoMEx (Combined score of DNA Methylation and Expression), to assess differentially expressed and methylated (DEM) genes/miRNAs. Next, we identified co-methylation and co-expression modules by applying WGCNA (Weighted Gene Correlation Network Analysis) tool to these DEM genes/miRNAs. The module with the highest average Pearson’s Correlation Coefficient (PCC) after considering all pair-wise molecules (genes/miRNAs) included 91 molecules. By integrating both transcription factor and miRNA regulations, we constructed subtype-specific regulatory networks for NSE and SE. We identified four hub miRNAs (miR-182-5p, miR-520b, miR-520c-3p, and miR-7-5p), two hub TFs (MYC and SP1), and two genes (RECK and TERT) in the NSE-specific regulatory network, and two hub miRNAs (miR-182-5p and miR-338-3p), five hub TFs (ETS1, HIF1A, HNF1A, MYC, and SP1), and three hub genes (CDH1, CXCR4, and SNAI1) in the SE-specific regulatory network. miRNA (miR-182-5p) and two TFs (MYC and SP1) were common hubs of NSE and SE. We further examined pathways enriched in these subtype-specific networks. Our study provides a comprehensive view of the molecular signatures and co-regulation in two TGCT subtypes.

Amyotrophic lateral sclerosis (ALS) is a fatal progressive multisystem disorder with limited therapeutic options. Although genome-wide association studies (GWASs) have revealed multiple ALS susceptibility loci, the exact identities of causal variants, genes, cell types, tissues, and their functional roles in the development of ALS remain largely unknown. Here, we reported a comprehensive post-GWAS analysis of the recent large ALS GWAS (n = 80,610), including functional mapping and annotation (FUMA), transcriptome-wide association study (TWAS), colocalization (COLOC), and summary data-based Mendelian randomization analyses (SMR) in extensive multi-omics datasets. Gene property analysis highlighted inhibitory neuron 6, oligodendrocytes, and GABAergic neurons (Gad1/Gad2) as functional cell types of ALS and confirmed cerebellum and cerebellar hemisphere as functional tissues of ALS. Functional annotation detected the presence of multiple deleterious variants at three loci (9p21.2, 12q13.3, and 12q14.2) and highlighted a list of SNPs that are potentially functional. TWAS, COLOC, and SMR identified 43 genes at 24 loci, including 23 novel genes and 10 novel loci, showing significant evidence of causality. Integrating multiple lines of evidence, we further proposed that rs2453555 at 9p21.2 and rs229243 at 14q12 functionally contribute to the development of ALS by regulating the expression of C9orf72 in pituitary and SCFD1 in skeletal muscle, respectively. Together, these results advance our understanding of the biological etiology of ALS, feed into new therapies, and provide a guide for subsequent functional experiments.

Abstract Background Psychiatric disorders such as schizophrenia (SCZ), bipolar disorder (BIP), major depressive disorder (MDD), attention deficit-hyperactivity disorder (ADHD), and autism spectrum disorder (ASD) are often related to brain development. Both shared and unique biological and neurodevelopmental processes have been reported to be involved in these disorders. Methods In this work, we developed an integrative analysis framework to seek for the sensitive spatiotemporal point during brain development underlying each disorder. Specifically, we first identified spatiotemporal gene co-expression modules for four brain regions three developmental stages (prenatal, birth to 11 years old, and older than 13 years), totaling 12 spatiotemporal sites. By integrating GWAS summary statistics and the spatiotemporal co-expression modules, we characterized the risk genes and their co-expression partners for five disorders. Results We found that SCZ and BIP, ASD and ADHD tend to cluster with each other and keep a distance from other psychiatric disorders. At the gene level, we identified several genes that were shared among the most significant modules, such as CTNNB1 and LNX1, and a hub gene, ATF2, in multiple modules. Moreover, we pinpointed two spatiotemporal points in the prenatal stage with active expression activities and highlighted one postnatal point for BIP. Further functional analysis of the disorder-related module highlighted the apoptotic signaling pathway for ASD and the immune-related and cell-cell adhesion function for SCZ, respectively. Conclusion Our study demonstrated the dynamic changes of disorder-related genes at the network level, shedding light on the spatiotemporal regulation during brain development.

Long non-coding RNAs (lncRNAs) have recently been found to be important in gene regulation. lncRNA H19 has been reported to play an oncogenic role in many human cancers. Its specific regulatory role is still elusive. In this study, we developed a novel analytic approach by integrating the synergistic regulation among lncRNAs (e.g., H19), transcription factors (TFs), target genes, and microRNAs (miRNAs) and then applied it to the pan-cancer expression datasets from The Cancer Genome Atlas (TCGA). Using linear regression models, we identified 88 H19-TF-gene co-regulatory triplets, in which 93% of the TF-gene pairs were related to cancer, indicating that our approach was effective to identify disease-related lncRNA-TF-gene co-regulation mechanisms. lncRNAs can function as miRNA sponges. Our further experiments found that H19 might regulate SP1-TGFBR2 through let-7b and miR-200b, ETS1-TGFBR2 through miR-29a and miR-200b, and STAT3-KLF11 through miR-17 in breast cancer cell lines. Our work suggests that miRNA-mediated lncRNA-TF-gene co-regulation is complicated yet important in cancer.

Abstract Background Multiple sclerosis (MS) is a complex disease in which the immune system attacks the central nervous system. The molecular mechanisms contributing to the etiology of MS remain poorly understood. Genome-wide association studies (GWAS) of MS have identified a small number of genetic loci significant at the genome level, but they are mainly non-coding variants. Network-assisted analysis may help better interpret the functional roles of the variants with association signals and potential translational medicine application. The Dense Module Searching of GWAS tool (dmGWAS version 2.4) developed in our team is applied to 2 MS GWAS datasets (GeneMSA and IMSGC GWAS) using the human protein interactome as the reference network. A dual evaluation strategy is used to generate results with reproducibility. Results Approximately 7500 significant network modules were identified for each independent GWAS dataset, and 20 significant modules were identified from the dual evaluation. The top modules included GRB2, HDAC1, JAK2, MAPK1, and STAT3 as central genes. Top module genes were enriched with functional terms such as “regulation of glial cell differentiation” (adjusted p-value = 2.58 × 10− 3), “T-cell costimulation” (adjusted p-value = 2.11 × 10− 6) and “virus receptor activity” (adjusted p-value = 1.67 × 10− 3). Interestingly, top gene networks included several MS FDA approved drug target genes HDAC1, IL2RA, KEAP1, and RELA, Conclusions Our dmGWAS network analyses highlighted several genes (GRB2, HDAC1, IL2RA, JAK2, KEAP1, MAPK1, RELA and STAT3) in top modules that are promising to interpret GWAS signals and link to MS drug targets. The genes enriched with glial cell differentiation are important for understanding neurodegenerative processes in MS and for remyelination therapy investigation. Importantly, our identified genetic signals enriched in T cell costimulation and viral receptor activity supported the viral infection onset hypothesis for MS.

Abstract Background Cleft lip with or without cleft palate (CL/P) is one of the most common congenital human birth defects. A combination of genetic and epidemiology studies has contributed to a better knowledge of CL/P-associated candidate genes and environmental risk factors. However, the etiology of CL/P remains not fully understood. In this study, to identify new CL/P-associated genes, we conducted an integrative analysis using our in-house network tools, dmGWAS [dense module search for Genome-Wide Association Studies (GWAS)] and EW_dmGWAS (Edge-Weighted dmGWAS), in a combination with GWAS data, the human protein-protein interaction (PPI) network, and differential gene expression profiles. Results A total of 87 genes were consistently detected in both European and Asian ancestries in dmGWAS. There were 31.0% (27/87) showed nominal significance with CL/P (gene-based p

Salmonella enterica serovar Typhimurium (S. Typhimurium) is one of the leading causes of food-borne illnesses worldwide. To colonize the gastrointestinal tract, S. Typhimurium produces multiple virulence factors that facilitate cellular invasion. Chitinases have been recently emerging as virulence factors for various pathogenic bacterial species, and the S. Typhimurium genome contains two annotated chitinases: STM0018 (chiA) and STM0233. However, the role of these chitinases during S. Typhimurium pathogenesis is unknown. The putative chitinase STM0233 has not been studied previously, and only limited data exists on ChiA. Chitinases typically hydrolyze chitin polymers, which are absent in vertebrates. However, chiA expression was detected in infection models and purified ChiA cleaved carbohydrate subunits present on mammalian surface glycoproteins, indicating a role during pathogenesis. Here, we demonstrate that expression of chiA and STM0233 is upregulated in the mouse gut and that both chitinases facilitate epithelial cell adhesion and invasion. S. Typhimurium lacking both chitinases showed a 70% reduction in invasion of small intestinal epithelial cells in vitro. In a gastroenteritis mouse model, chitinase-deficient S. Typhimurium strains were also significantly attenuated in the invasion of small intestinal tissue. This reduced invasion resulted in significantly delayed S. Typhimurium dissemination to the spleen and the liver, but chitinases were not required for systemic survival. The invasion defect of the chitinase-deficient strain was rescued by the presence of wild-type S. Typhimurium, suggesting that chitinases are secreted. By analyzing N-linked glycans of small intestinal cells, we identified specific N-acetylglucosamine-containing glycans as potential extracellular targets of S. Typhimurium chitinases. This analysis also revealed a differential abundance of Lewis X/A-containing glycans that is likely a result of host cell modulation due to the detection of S. Typhimurium chitinases. Similar glycomic changes elicited by chitinase deficient strains indicate functional redundancy of the chitinases. Overall, our results demonstrate that S. Typhimurium chitinases contribute to intestinal adhesion and invasion through modulation of the host glycome.

Abstract Background Heart failure (HF) has been recognized as a global pandemic with a high rate of hospitalization, morbidity, and mortality. Although numerous advances have been made, its representative molecular signatures remain largely unknown, especially the role of genes in HF progression. The aim of the present prospective follow-up study was to reveal potential biomarkers associated with the progression of heart failure. Methods We generated multi-level transcriptomic data from a cohort of left ventricular heart tissue collected from 21 HF patients and 9 healthy donors. By using Masson staining to calculate the fibrosis percentage for each sample, we applied lasso regression model to identify the genes associated with fibrosis as well as progression. The genes were further validated by immunohistochemistry (IHC) staining in the same cohort and qRT-PCR using another independent cohort (20 HF and 9 healthy donors). Enzyme-linked immunosorbent assay (ELISA) was used to measure the plasma level in a validation cohort (139 HF patients) for predicting HF progression. Results Based on the multi-level transcriptomic data, we examined differentially expressed genes [mRNAs, microRNAs, and long non-coding RNAs (lncRNAs)] in the study cohort. The follow-up functional annotation and regulatory network analyses revealed their potential roles in regulating extracellular matrix. We further identified several genes that were associated with fibrosis. By using the survival time before transplantation, COL1A1 was identified as a potential biomarker for HF progression and its upregulation was confirmed by both IHC and qRT-PCR. Furthermore, COL1A1 content ≥ 256.5 ng/ml in plasma was found to be associated with poor survival within 1 year of heart transplantation from heart failure [hazard ratio (HR) 7.4, 95% confidence interval (CI) 3.5 to 15.8, Log-rank p value

Abstract Background Cleft lip (CL), one of the most common congenital birth defects, shows considerable geographic and ethnic variation, with contribution of both genetic and environmental factors. Mouse genetic studies have identified several CL-associated genes. However, it remains elusive how these CL-associated genes are regulated and involved in CL. Environmental factors may regulate these genes at the post-transcriptional level through the regulation of non-coding microRNAs (miRNAs). In this study, we sought to identify miRNAs associated with CL in mice. Results Through a systematic literature review and a Mouse Genome Informatics (MGI) database search, we identified 55 genes that were associated with CL in mice. Subsequent bioinformatic analysis of these genes predicted that a total of 33 miRNAs target multiple CL-associated genes, with 20 CL-associated genes being potentially regulated by multiple miRNAs. To experimentally validate miRNA function in cell proliferation, we conducted cell proliferation/viability assays for the selected five candidate miRNAs (miR-124-3p, let-7a-5p, let-7b-5p, let-7c-5p, and let-7d-5p). Overexpression of miR-124-3p, but not of the others, inhibited cell proliferation through suppression of CL-associated genes in cultured mouse embryonic lip mesenchymal cells (MELM cells) isolated from the developing mouse lip region. By contrast, miR-124-3p knockdown had no effect on MELM cell proliferation. This miRNA-gene regulatory mechanism was mostly conserved in O9–1 cells, an established cranial neural crest cell line. Expression of miR-124-3p was low in the maxillary processes at E10.5, when lip mesenchymal cells proliferate, whereas it was greatly increased at later developmental stages, suggesting that miR-124-3p expression is suppressed during the proliferation phase in normal palate development. Conclusions Our findings indicate that upregulated miR-124-3p inhibits cell proliferation in cultured lip cells through suppression of CL-associated genes. These results will have a significant impact, not only on our knowledge about lip morphogenesis, but also on the development of clinical approaches for the diagnosis and prevention of CL.

Abstract Background The etiology of cleft lip with or without cleft palate (CL/P), a common congenital birth defect, is complex and involves the contribution of genetic and environmental factors. Although many candidate genes have been identified, the regulation and interaction of these genes in CL/P remain unclear. In addition, the contribution of microRNAs (miRNAs), non-coding RNAs that regulate the expression of multiple genes, to the etiology of CL/P is largely unknown. Methods To identify the signatures of causative biological pathways for human CL/P, we conducted a systematic literature review for human CL/P candidate genes and subsequent bioinformatics analyses. Functional enrichment analyses of the candidate CL/P genes were conducted using the pathway databases GO and KEGG. The miRNA-mediated post-transcriptional regulation of the CL/P candidate genes was analyzed with miRanda, PITA, and TargetScan, and miRTarbase. Genotype-phenotype association analysis was conducted using GWAS. The functional significance of the candidate miRNAs was evaluated experimentally in cell proliferation and target gene regulation assays in human lip fibroblasts. Results Through an extensive search of the main biomedical databases, we mined 177 genes with mutations or association/linkage reported in individuals with CL/P, and considered them as candidate genes for human CL/P. The genotype-phenotype association study revealed that mutations in 12 genes (ABCA4, ADAM3A, FOXE1, IRF6, MSX2, MTHFR, NTN1, PAX7, TP63, TPM1, VAX1, and WNT9B) were significantly associated with CL/P. In addition, our bioinformatics analysis predicted 16 microRNAs (miRNAs) to be post-transcriptional regulators of CL/P genes. To validate the bioinformatics results, the top six candidate miRNAs (miR-124-3p, miR-369-3p, miR-374a-5p, miR-374b-5p, miR-497-5p, and miR-655-3p) were evaluated by cell proliferation/survival assays and miRNA-gene regulation assays in cultured human lip fibroblasts. We found that miR-497-5p and miR-655-3p significantly suppressed cell proliferation in these cells. Furthermore, the expression of the predicted miRNA-target genes was significantly downregulated by either miR-497-5p or miR-655-3p mimic. Conclusion Expression of miR-497-5p and miR-655-3p suppresses cell proliferation through the regulation of human CL/P-candidate genes. This study provides insights into the role of miRNAs in the etiology of CL/P and suggests possible strategies for the diagnosis of CL/P.

Clear cell renal cell carcinoma (ccRCC) is a urogenital cancer with a well-defined genetic landscape. Here, the authors analyse the clonal architecture of ccRCC patients from three populations, and find prognostic subtypes linked to immune infiltrates and clonal architecture.

Abstract Background Genome-wide association studies (GWAS) have been successful in identifying disease-associated genetic variants. Recently, an increasing number of GWAS summary statistics have been made available to the research community, providing extensive repositories for studies of human complex diseases. In particular, cross-trait associations at the genetic level can be beneficial from large-scale GWAS summary statistics by using genetic variants that are associated with multiple traits. However, direct assessment of cross-trait associations using susceptibility loci has been challenging due to the complex genetic architectures in most diseases, calling for advantageous methods that could integrate functional interpretation and imply biological mechanisms. Results We developed an analytical framework for systematic integration of cross-trait associations. It incorporates two different approaches to detect enriched pathways and requires only summary statistics. We demonstrated the framework using 25 traits belonging to four phenotype groups. Our results revealed an average of 54 significantly associated pathways (ranged between 18 and 175) per trait. We further proved that pathway-based analysis provided increased power to estimate cross-trait associations compared to gene-level analysis. Based on Fisher’s Exact Test (FET), we identified a total of 24 (53) pairs of trait-trait association at adjusted p FET

Abstract Background Existing functional description of genes are categorical, discrete, and mostly through manual process. In this work, we explore the idea of gene embedding, distributed representation of genes, in the spirit of word embedding. Results From a pure data-driven fashion, we trained a 200-dimension vector representation of all human genes, using gene co-expression patterns in 984 data sets from the GEO databases. These vectors capture functional relatedness of genes in terms of recovering known pathways - the average inner product (similarity) of genes within a pathway is 1.52X greater than that of random genes. Using t-SNE, we produced a gene co-expression map that shows local concentrations of tissue specific genes. We also illustrated the usefulness of the embedded gene vectors, laden with rich information on gene co-expression patterns, in tasks such as gene-gene interaction prediction. Conclusions We proposed a machine learning method that utilizes transcriptome-wide gene co-expression to generate a distributed representation of genes. We further demonstrated the utility of our distribution by predicting gene-gene interaction based solely on gene names. The distributed representation of genes could be useful for more bioinformatics applications.

Abstract Approximately 15% of human cancers are estimated to be attributed to viruses. Virus sequences can be integrated into the host genome, leading to genomic instability and carcinogenesis. Here, a new deep convolutional neural network (CNN) model is developed with attention architecture, namely DeepVISP, for accurately predicting oncogenic virus integration sites (VISs) in the human genome. Using the curated benchmark integration data of three viruses, hepatitis B virus (HBV), human herpesvirus (HPV), and Epstein‐Barr virus (EBV), DeepVISP achieves high accuracy and robust performance for all three viruses through automatically learning informative features and essential genomic positions only from the DNA sequences. In comparison, DeepVISP outperforms conventional machine learning methods by 8.43–34.33% measured by area under curve (AUC) value enhancement in three viruses. Moreover, DeepVISP can decode cis‐regulatory factors that are potentially involved in virus integration and tumorigenesis, such as HOXB7, IKZF1, and LHX6. These findings are supported by multiple lines of evidence in literature. The clustering analysis of the informative motifs reveales that the representative k‐mers in clusters could help guide virus recognition of the host genes. A user‐friendly web server is developed for predicting putative oncogenic VISs in the human genome using DeepVISP.

Background: Genome-wide association studies have successfully identified variants associated with multiple conditions. However, generalizing discoveries across diverse populations remains challenging due to large variations in genetic composition. Methods that perform gene expression imputation have attempted to address the transferability of gene discoveries across populations, but with limited success. Methods: Here, we introduce a pipeline that combines gene expression imputation with gene module discovery, including a dense gene module search and a gene set variation analysis, to address the transferability issue. Our method feeds association probabilities of imputed gene expression with a selected phenotype into tissue-specific gene-module discovery over protein interaction networks to create higher-level gene modules. Results: We demonstrate our method’s utility in three case-control studies of Alzheimer’s disease (AD) for three different race/ethnic populations (Whites, African descent and Hispanics). We discovered 182 AD-associated genes from gene modules shared between these populations, highlighting new gene modules associated with AD. Conclusions: Our innovative framework has the potential to identify robust discoveries across populations based on gene modules, as demonstrated in AD.

With steady increase of new COVID-19 cases around the world, especially in the United States, health care resources in areas with the disease outbreak are quickly exhausted by overwhelming numbers of COVID-19 patients. Therefore, strategies that can effectively and quickly predict the disease progression and stratify patients for appropriate health care arrangements are urgently needed. We explored the features and evolutionary difference of viral gene expression in the SARS-CoV-2 infected cells from the bronchoalveolar lavage fluids of patients with moderate and severe COVID-19 using both single cell and bulk tissue transcriptome data. We found SARS-CoV-2 sequences were detectable in 8 types of immune related cells, including macrophages, T cells, and NK cells. We first reported that the SARS-CoV-2 ORF10 gene was differentially expressed in the severe vs. moderate samples. Specifically, ORF10 was abundantly expressed in infected cells of severe cases, while it was barely detectable in the infected cells of moderate cases. Consequently, the expression ratio of ORF10 to nucleocapsid (N) was significantly higher in severe than moderate cases (p = 0.0062). Moreover, we found transcription regulatory sequences (TRSs) of the viral leader sequence-independent fusions with a 5’ joint point at position 1073 of SARS-CoV-2 genome were detected mainly in the patients with death outcome, suggesting its potential indication of clinical outcome. Finally, we identified the motifs in TRS of the viral leader sequence-dependent fusion events of SARS-CoV-2 and compared with that in SARS-CoV, suggesting its evolutionary trajectory. These results implicated potential roles and predictive features of viral transcripts in the pathogenesis of COVID-19 moderate and severe patients. Such features and evolutionary patterns require more data to validate in future.

Abstract Background Cleft lip (CL) is one of the most common congenital birth defects with complex etiology. While genome-wide association studies (GWAS) have made significant advances in our understanding of mutations and their related genes with potential involvement in the etiology of CL, it remains unknown how these genes are functionally regulated and interact with each other in lip development. Currently, identifying the disease-causing genes in human CL is urgently needed. So far, the causative CL genes have been largely undiscovered, making it challenging to design experiments to validate the functional influence of the mutations identified from large genomic studies such as CL GWAS. Results Transcription factors (TFs) and microRNAs (miRNAs) are two important regulators in cellular system. In this study, we aimed to investigate the genetic interactions among TFs, miRNAs and the CL genes curated from the previous studies. We constructed miRNA-TF co-regulatory networks, from which the critical regulators as putative drivers in CL were examined. Based on the constructed networks, we identified ten critical hub genes with prior evidence in CL. Furthermore, the analysis of partitioned regulatory modules highlighted a number of biological processes involved in the pathology of CL, including a novel pathway “Signaling pathway regulating pluripotency of stem cells”. Our subnetwork analysis pinpointed two candidate miRNAs, hsa-mir-27b and hsa-mir-497, activating the Wnt pathway that was associated with CL. Our results were supported by an independent gene expression dataset in CL. Conclusions This study represents the first regulatory network analysis of CL genes. Our work presents a global view of the CL regulatory network and a novel approach on investigating critical miRNAs, TFs and genes via combinatory regulatory networks in craniofacial development. The top genes and miRNAs will be important candidates for future experimental validation of their functions in CL.

Abstract Background There are two main types of lung cancer: small cell lung cancer (SCLC) and non-small cell lung cancer (NSCLC). NSCLC has many subtypes, but the two most common are lung adenocarcinoma (LUAD) and lung squamous cell carcinoma (LUSC). These subtypes are mainly classified by physiological and pathological characteristics, although there is increasing evidence of genetic and molecular differences as well. Although some work has been done at the somatic level to explore the genetic and biological differences among subtypes, little work has been done that interrogates these differences at the germline level to characterize the unique and shared susceptibility genes for each subtype. Methods We used single-nucleotide polymorphisms (SNPs) from a genome-wide association study (GWAS) of European samples to interrogate the similarity of the subtypes at the SNP, gene, pathway, and regulatory levels. We expanded these genotyped SNPs to include all SNPs in linkage disequilibrium (LD) using data from the 1000 Genomes Project. We mapped these SNPs to several lung tissue expression quantitative trait loci (eQTL) and enhancer datasets to identify regulatory SNPs and their target genes. We used these genes to perform a biological pathway analysis for each subtype. Results We identified 8295, 8734, and 8361 SNPs with moderate association signals for LUAD, LUSC, and SCLC, respectively. Those SNPs had p < 1 × 10− 3 in the original GWAS or were within LD (r 2 > 0.8, Europeans) to the genotyped SNPs. We identified 215, 320, and 172 disease-associated genes for LUAD, LUSC, and SCLC, respectively. Only five genes (CHRNA5, IDH3A, PSMA4, RP11-650 L12.2, and TBC1D2B) overlapped all subtypes. Furthermore, we observed only two pathways from the Kyoto Encyclopedia of Genes and Genomes shared by all subtypes. At the regulatory level, only three eQTL target genes and two enhancer target genes overlapped between all subtypes. Conclusions Our results suggest that the three lung cancer subtypes do not share much genetic signal at the SNP, gene, pathway, or regulatory level, which differs from the common subtype classification based upon histology. However, three (CHRNA5, IDH3A, and PSMA4) of the five genes shared between the subtypes are well-known lung cancer genes that may act as general lung cancer genes regardless of subtype.

Abstract Background Genome–phenome studies have identified thousands of variants that are statistically associated with disease or traits; however, their functional roles are largely unclear. A comprehensive investigation of regulatory mechanisms and the gene regulatory networks between phenome-wide association study (PheWAS) and genome-wide association study (GWAS) is needed to identify novel regulatory variants contributing to risk for human diseases. Methods In this study, we developed an integrative functional genomics framework that maps 215,107 significant single nucleotide polymorphism (SNP) traits generated from the PheWAS Catalog and 28,870 genome-wide significant SNP traits collected from the GWAS Catalog into a global human genome regulatory map via incorporating various functional annotation data, including transcription factor (TF)-based motifs, promoters, enhancers, and expression quantitative trait loci (eQTLs) generated from four major functional genomics databases: FANTOM5, ENCODE, NIH Roadmap, and Genotype-Tissue Expression (GTEx). In addition, we performed a tissue-specific regulatory circuit analysis through the integration of the identified regulatory variants and tissue-specific gene expression profiles in 7051 samples across 32 tissues from GTEx. Results We found that the disease-associated loci in both the PheWAS and GWAS Catalogs were significantly enriched with functional SNPs. The integration of functional annotations significantly improved the power of detecting novel associations in PheWAS, through which we found a number of functional associations with strong regulatory evidence in the PheWAS Catalog. Finally, we constructed tissue-specific regulatory circuits for several complex traits: mental diseases, autoimmune diseases, and cancer, via exploring tissue-specific TF-promoter/enhancer-target gene interaction networks. We uncovered several promising tissue-specific regulatory TFs or genes for Alzheimer’s disease (e.g. ZIC1 and STX1B) and asthma (e.g. CSF3 and IL1RL1). Conclusions This study offers powerful tools for exploring the functional consequences of variants generated from genome–phenome association studies in terms of their mechanisms on affecting multiple complex diseases and traits.

Abstract Background In 2016, it is estimated that there will be 62,700 new cases of kidney cancer in the United States, and 14,240 patients will die from the disease. Because the incidence of kidney renal clear cell carcinoma (KIRC), the most common type of kidney cancer, is expected to continue to increase in the US, there is an urgent need to find effective diagnostic biomarkers for KIRC that could help earlier detection of and customized treatment strategies for the disease. Accordingly, in this study we systematically investigated KIRC’s prognostic biomarkers for survival using the reverse phase protein array (RPPA) data and the high throughput sequencing data from The Cancer Genome Atlas (TCGA). Results With comprehensive data available in TCGA, we systematically screened protein expression based survival biomarkers in 10 major cancer types, among which KIRC presented many protein prognostic biomarkers of survival time. This is in agreement with a previous report that expression level changes (mRNAs, microRNA and protein) may have a better performance for prognosis of KIRC. In this study, we also identified 52 prognostic genes for KIRC, many of which are involved in cell-cycle and cancer signaling, as well as 15 tumor-stage-specific prognostic biomarkers. Notably, we found fewer prognostic biomarkers for early-stage than for late-stage KIRC. Four biomarkers (the RPPA protein IDs: FASN, ACC1, Cyclin_B1 and Rad51) were found to be prognostic for survival based on both protein and mRNA expression data. Conclusions Through pan-cancer screening, we found that many protein biomarkers were prognostic for patients’ survival in KIRC. Stage-specific survival biomarkers in KIRC were also identified. Our study indicated that these protein biomarkers might have potential clinical value in terms of predicting survival in KIRC patients and developing individualized treatment strategies. Importantly, we found many biomarkers in KIRC at both the mRNA expression level and the protein expression level. These biomarkers shared a significant overlap, indicating that they were technically replicable.

Abstract Background HIV-associated neurocognitive disorder (HAND) remains common, despite antiretroviral therapy (ART). HIV dysregulates iron metabolism, but cerebrospinal fluid (CSF) levels of iron and iron-transport proteins in HIV-infected (HIV+) persons are largely unknown. The objectives of this study were to characterize CSF iron-related biomarkers in HIV+ adults and explore their relationships to known predictors of HAND. Methods We quantified total iron, transferrin and heavy-chain (H)-ferritin by immunoassay in CSF sampled by lumbar puncture in 403 HIV+ participants in a multi-center, observational study and evaluated biomarker associations with demographic and HIV-related correlates of HAND [e.g., age, sex, self-reported race/ethnicity, ART, and detectable plasma virus and CSF viral load (VL)] by multivariable regression. In a subset (N = 110) with existing CSF: serum albumin (QAlb) measurements, QAlb and comorbidity severity were also included as covariates to account for variability in the blood–CSF-barrier. Results Among 403 individuals (median age 43 years, 19% women, 56% non-Whites, median nadir CD4+ T cell count 180 cells/µL, 46% with undetectable plasma virus), men had 25% higher CSF transferrin (median 18.1 vs. 14.5 µg/mL), and 71% higher H-ferritin (median 2.9 vs. 1.7 ng/mL) than women (both p-values ≤0.01). CSF iron was 41% higher in self-reported Hispanics and 27% higher in (non-Hispanic) Whites than in (non-Hispanic) Blacks (median 5.2 and 4.7 µg/dL in Hispanics and Whites, respectively, vs. 3.7 µg/dL in Blacks, both p ≤ 0.01); these findings persisted after adjustment for age, sex, and HIV-specific factors. Median H-ferritin was 25% higher (p 50 years) than in younger persons (age ≤ 35 years; both p

Crohn’s Disease (CD) is one of the predominant forms of inflammatory bowel disease (IBD). A combination of genetic and non-genetic risk factors have been reported to contribute to the development of CD. Many high-throughput omics studies have been conducted to identify disease associated risk variants that might contribute to CD, such as genome-wide association studies (GWAS) and next generation sequencing studies. A pressing need remains to prioritize and characterize candidate genes that underlie the etiology of CD. In this study, we collected a comprehensive multi-dimensional data from GWAS, gene expression, and methylation studies and generated transcriptome-wide association study (TWAS) data to further interpret the GWAS association results. We applied our previously developed method called mega-analysis of Odds Ratio (MegaOR) to prioritize CD candidate genes (CDgenes). As a result, we identified consensus sets of CDgenes (62–235 genes) based on the evidence matrix. We demonstrated that these CDgenes were significantly more frequently interact with each other than randomly expected. Functional annotation of these genes highlighted critical immune-related processes such as immune response, MHC class II receptor activity, and immunological disorders. In particular, the constitutive photomorphogenesis 9 (COP9) signalosome related genes were found to be significantly enriched in CDgenes, implying a potential role of COP9 signalosome involved in the pathogenesis of CD. Finally, we found some of the CDgenes shared biological functions with known drug targets of CD, such as the regulation of inflammatory response and the leukocyte adhesion to vascular endothelial cell. In summary, we identified highly confident CDgenes from multi-dimensional evidence, providing insights for the understanding of CD etiology.

Recently, accumulating evidence has supported that circular RNA (circRNA) plays important roles in tumorigenesis by regulating gene expression at transcriptional and post-transcriptional levels. Expression of circRNAs can be epigenetically silenced by DNA methylation; however, the underlying regulatory mechanisms of circRNAs by DNA methylation remains largely unknown. We explored this regulation in hepatocellular carcinoma (HCC) using genome-wide DNA methylation and RNA sequencing data of the primary tumor and matched adjacent normal tissues from 20 HCC patients. Our pipeline identified 1012 upregulated and 747 downregulated circRNAs (collectively referred to as differentially expressed circRNAs, or DE circRNAs) from HCC RNA-seq data. Among them, 329 DE circRNAs covered differentially methylated sites (adjusted p-value < 0.05, |ΔM| > 0.5) in circRNAs’ interior and/or flanking regions. Interestingly, the corresponding parental genes of 46 upregulated and 31 downregulated circRNAs did not show significant expression change in the HCC tumor versus normal samples. Importantly, 34 of the 77 DE circRNAs (44.2%) had significant correlation with DNA methylation change in HCC (Spearman’s rank-order correlation, p-value < 0.05), suggesting that aberrant DNA methylation might regulate circular RNA expression in HCC. Our study revealed genome-wide differential circRNA expression in HCC. The significant correlation with DNA methylation change suggested that epigenetic regulation might act on both mRNA and circRNA expression. The specific regulation in HCC and general view in other cancer or disease requires further investigation.