HMBOX1, homeobox transcription factor, negatively regulates interferon-γ production in natural killer cells

Inhibition of Hmbox1 Promotes Cardiomyocyte Survival and Glucose Metabolism Through Gck Activation in Ischemia/Reperfusion Injury

BACKGROUND

Exercise-induced physiological cardiac growth regulators may protect the heart from ischemia/reperfusion (I/R) injury. Homeobox-containing 1 (Hmbox1), a homeobox family member, has been identified as a putative transcriptional repressor and is downregulated in the exercised heart. However, its roles in exercise-induced physiological cardiac growth and its potential protective effects against cardiac I/R injury remain largely unexplored.

METHODS

We studied the function of Hmbox1 in exercise-induced physiological cardiac growth in mice after 4 weeks of swimming exercise. Hmbox1 expression was then evaluated in human heart samples from deceased patients with myocardial infarction and in the animal cardiac I/R injury model. Its role in cardiac I/R injury was examined in mice with adeno-associated virus 9 (AAV9) vector-mediated Hmbox1 knockdown and in those with cardiac myocyte-specific Hmbox1 ablation. We performed RNA sequencing, promoter prediction, and binding assays and identified glucokinase (Gck) as a downstream effector of Hmbox1. The effects of Hmbox1 together with Gck were examined in cardiomyocytes to evaluate their cell size, proliferation, apoptosis, mitochondrial respiration, and glycolysis. The function of upstream regulator of Hmbox1, ETS1, was investigated through ETS1 overexpression in cardiac I/R mice in vivo.

RESULTS

We demonstrated that Hmbox1 downregulation was required for exercise-induced physiological cardiac growth. Inhibition of Hmbox1 increased cardiomyocyte size in isolated neonatal rat cardiomyocytes and human embryonic stem cell-derived cardiomyocytes but did not affect cardiomyocyte proliferation. Under pathological conditions, Hmbox1 was upregulated in both human and animal postinfarct cardiac tissues. Furthermore, both cardiac myocyte-specific Hmbox1 knockout and AAV9-mediated Hmbox1 knockdown protected against cardiac I/R injury and heart failure. Therapeutic effects were observed when sh-Hmbox1 AAV9 was administered after I/R injury. Inhibition of Hmbox1 activated the Akt/mTOR/P70S6K pathway and transcriptionally upregulated Gck, leading to reduced apoptosis and improved mitochondrial respiration and glycolysis in cardiomyocytes. ETS1 functioned as an upstream negative regulator of Hmbox1 transcription, and its overexpression was protective against cardiac I/R injury.

CONCLUSIONS

Our studies unravel a new role for the transcriptional repressor Hmbox1 in exercise-induced physiological cardiac growth. They also highlight the therapeutic potential of targeting Hmbox1 to improve myocardial survival and glucose metabolism after I/R injury.

Uncovering structural variants in Creole cattle from Guadeloupe and their impact on environmental adaptation through whole genome sequencing

Structural variants play an important role in evolutionary processes. Besides, they constitute a large source of inter individual genetic variation that might represent a major factor in the aetiology of complex, multifactorial traits. Their importance in adaptation is becoming increasingly evident in literature. Yet, the characterization of the genomic landscape of structural variants in local breeds remains scarce to date. Herein, we investigate patterns and gene annotation of structural variants in the Creole cattle from Guadeloupe breed using whole genome sequences from 23 bulls representative of the population. In total, we detected 32821 ascertained SV defining 15258 regions, representing ~ 17% of the Creole cattle genome. Among these, 6639 regions have not been previously reported in the Database of Genomic Variants archive. Average number of structural variants detected per individual in the studied population is in the same order of magnitude of that observed in indicine populations and higher than that reported in taurine breeds. We observe an important within-individual variability where approximately half of the detected structural variants have low frequency (MAF < 0.25). Most of the detected structural variants (55%) occurred in intergenic regions. Genic structural variants overlapped with 7793 genes and the predicted effect of most of them is ranked as "modifier". Among the structural variants that were predicted to have a high functional impact on the protein, a 5.5 Kb in length, highly frequent deletion on chromosome 2, affects ALPI, a gene associated with the interaction between gut microbiota and host immune system. The 6639 newly identified structural variants regions include three deletions and three duplications shared by more than 80% of individuals that are significantly enriched for genes related to tRNA threonylcarbamoyladenosine metabolic process, important for temperature adaptation in thermophilic organisms, therefore suggesting a potential role in the thermotolerance of Creole cattle from Guadeloupe cattle to tropical climate. Overall, highly frequent structural variants that are specific to the Creole cattle population encompass olfactory receptor and immunity genes as well as genes involved in muscle tone, muscle development and contraction. Beyond mapping and characterizing structural variants in the Creole cattle from Guadeloupe breed, this study provides valuable information for a better understanding of the potential role of chromosomal rearrangements in adaptive traits in cattle.

Metabolomic analysis reveals the influence of HMBOX1 on RAW264.7 cells proliferation based on UPLC-MS/MS

Macrophages are important effector cells in tumor progression and immune regulation. Previously, we demonstrated that the transcription suppressor homeobox containing 1(HMBOX1) exhibits immunosuppressive activity in LPS-induced acute liver injury by impeding macrophage infiltration and activation. We also observed a lower proliferation in HMBOX1-overexpressed RAW264.7 cells. However, the specific mechanism was unclear. Here, a work was performed to characterize HMBOX1 function related to cell proliferation from a metabolomics standpoint by comparing the metabolic profiles of HMBOX1-overexpressed RAW264.7 cells to those of the controls. Firstly, we assessed HMBOX1 anti-proliferation activity in RAW264.7 cells with CCK8 assay and clone formation. Then, we performed metabolomic analyses by ultra-liquid chromatography coupled with mass spectrometry to explore the potential mechanisms. Our results indicated that HMBOX1 inhibited the macrophage growth curve and clone formation ability. Metabolomic analyses showed significant changes in HMBOX1-overexpressed RAW264.7 metabolites. A total of 1312 metabolites were detected, and 185 differential metabolites were identified based on the criterion of OPLS-DA VIP > 1 and p value < 0.05. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis indicated that the elevated HMBOX1 in RAW264.7 inhibited the pathways of amino acid and nucleotide metabolism. Glutamine concentrations decreased significantly in HMBOX1-overexpressed macrophages, and glutamine-related transporter SLC1A5 was also downregulated. Furthermore, SLC1A5 overexpression reversed HMBOX1 inhibition of macrophage proliferation. This study demonstrated the potential mechanism of the HMBOX1/SLC1A5 pathway in cell proliferation by regulating glutamine transportation. The results may help provide a new direction for therapeutic interventions in macrophage-related inflammatory diseases.

Single Nucleotide Polymorphisms of ALDH18A1 and MAT2A Genes and Their Genetic Associations with Milk Production Traits of Chinese Holstein Cows

Our preliminary work had suggested two genes, aldehyde dehydrogenase 18 family member A1 (ALDH18A1) and methionine adenosyltransferase 2A (MAT2A), related to amino acid synthesis and metabolism as candidates affecting milk traits by analyzing the liver transcriptome and proteome of dairy cows at different lactation stages. In this study, the single nucleotide polymorphisms (SNPs) of ALDH18A1 and MAT2A genes were identified and their genetic effects and underlying causative mechanisms on milk production traits in dairy cattle were analyzed, with the aim of providing effective genetic information for the molecular breeding of dairy cows. By resequencing the entire coding and partial flanking regions of ALDH18A1 and MAT2A, we found eight SNPs located in ALDH18A1 and two in MAT2A. Single-SNP association analysis showed that most of the 10 SNPs of these two genes were significantly associated with the milk yield traits, 305-day milk yield, fat yield, and protein yield in the first and second lactations (corrected p ≤ 0.0488). Using Haploview 4.2, we found that the seven SNPs of ALDH18A1 formed two haplotype blocks; subsequently, the haplotype-based association analysis showed that both haplotypes were significantly associated with 305-day milk yield, fat yield, and protein yield (corrected p ≤ 0.014). Furthermore, by Jaspar and Genomatix software, we found that 26:g.17130318 C>A and 11:g.49472723G>C, respectively, in the 5′ flanking region of ALDH18A1 and MAT2A genes changed the transcription factor binding sites (TFBSs), which might regulate the expression of corresponding genes to affect the phenotypes of milk production traits. Therefore, these two SNPs were considered as potential functional mutations, but they also require further verification. In summary, ALDH18A1 and MAT2A were proved to probably have genetic effects on milk production traits, and their valuable SNPs might be used as candidate genetic markers for dairy cattle’s genomic selection (GS).

Swine MicroRNAs ssc-miR-221-3p and ssc-miR-222 Restrict the Cross-Species Infection of Avian Influenza Virus

Avian influenza virus (AIV) can cross species barriers to infect humans and other mammals. However, these species-cross transmissions are most often dead-end infections due to host restriction. Current research about host restriction focuses mainly on the barriers of cell membrane, nuclear envelope, and host proteins; whether microRNAs (miRNAs) play a role in host restriction is largely unknown. In this study, we used porcine alveolar macrophage (PAM) cells as a model to elucidate the role of miRNAs in host range restriction. During AIV infection, 40 dysregulation expressed miRNAs were selected in PAM cells. Among them, two Sus scrofa (ssc; swine) miRNAs, ssc-miR-221-3p and ssc-miR-222, could inhibit the infection and replication of AIV in PAM cells by directly targeting viral genome and inducing cell apoptosis via inhibiting the expression of anti-apoptotic protein HMBOX1. Avian but not swine influenza virus caused upregulated expressions of ssc-miR-221-3p and ssc-miR-222 in PAM cells. We further found that NF-κB P65 was more effectively phosphorylated upon AIV infection and that P65 functioned as a transcription activator to regulate the AIV-induced expression of miR-221-3p/222 Importantly, we found that ssc-miR-221-3p and ssc-miR-222 could also be specifically upregulated upon AIV infection in newborn pig tracheal epithelial (NPTr) cells and also exerted anti-AIV function. In summary, our study indicated that miRNAs act as a host barrier during cross-species infection of influenza A virus.IMPORTANCE The host range of an influenza A virus is determined by species-specific interactions between virus and host cell factors. Host miRNAs can regulate influenza A virus replication; however, the role of miRNAs in host species specificity is unclear. Here, we show that the induced expression of ssc-miR-221-3p and ssc-miR-222 in swine cells is modulated by NF-κB P65 phosphorylation in response to AIV infection but not swine influenza virus infection. ssc-miR-221-3p and ssc-miR-222 exerted antiviral function via targeting viral RNAs and causing apoptosis by inhibiting the expression of HMBOX1 in host cells. These findings uncover miRNAs as a host range restriction factor that limits cross-species infection of influenza A virus.

WTAP promotes osteosarcoma tumorigenesis by repressing HMBOX1 expression in an m6A-dependent manner

N⁶-methyladenosine (m⁶A) regulators are involved in the progression of various cancers via regulating m⁶A modification. However, the potential role and mechanism of the m⁶A modification in osteosarcoma remains obscure. In this study, WTAP was found to be highly expressed in osteosarcoma tissue and it was an independent prognostic factor for overall survival in osteosarcoma. Functionally, WTAP, as an oncogene, was involved in the proliferation and metastasis of osteosarcoma in vitro and vivo. Mechanistically, M⁶A dot blot, RNA-seq and MeRIP-seq, MeRIP-qRT-PCR and luciferase reporter assays showed that HMBOX1 was identified as the target gene of WTAP, which regulated HMBOX1 stability depending on m⁶A modification at the 3′UTR of HMBOX1 mRNA. In addition, HMBOX1 expression was downregulated in osteosarcoma and was an independent prognostic factor for overall survival in osteosarcoma patients. Silenced HMBOX1 evidently attenuated shWTAP-mediated suppression on osteosarcoma growth and metastasis in vivo and vitro. Finally, WTAP/HMBOX1 regulated osteosarcoma growth and metastasis via PI3K/AKT pathway. In conclusion, this study demonstrated the critical role of the WTAP-mediated m⁶A modification in the progression of osteosarcoma, which could provide novel insights into osteosarcoma treatment.

Loss of HMBOX1 promotes LPS-induced apoptosis and inhibits LPS-induced autophagy of vascular endothelial cells in mouse

Our previous study revealed that Homeobox containing 1 (HMBOX1), essential for the survival of vascular endothelial cells (VECs), was involved in the progression of atherosclerosis. Knockdown of HMBOX1 promoted apoptosis and inhibited autophagy through regulating intracellular free zinc level in cultured VECs. In current study, in order to investigate the roles of HMBOX1 in vivo and in endothelium, we generated a knockout (KO) mouse for HMBOX1 by using transcription activator-like effector nucleases (TALENs) technology. Herein, we reported that the protein level of HMBOX1 was gradually increased during mouse development. The HMBOX1 KO mouse was viable and fertile. There existed no differences in apoptosis and autophagy of aortic endothelial cells between wild type and KO mouse. Whereas, loss of HMBOX1 promoted apoptosis and inhibited autophagy of aortic endothelial cells under lipopolysaccharide (LPS) stimulation in mouse. We also demonstrated that HMBOX1 deletion had no influence on the secretion of inflammatory cytokines TNF-α and IL-6. Moreover, overexpression or knockdown of HMBOX1 failed to regulate multiple pro-apoptotic genes expression in vitro. In conclusion, HMBOX1 participated in the functional maintenance of mouse aortic endothelial cells, the aortic endothelial cells of HMBOX1 KO mouse showed increased apoptosis and decreased autophagy with LPS treatment.

Prediction and analysis of weighted genes in isoflurane induced general anesthesia based on network analysis

Purpose: Isoflurane is still wildly used in the developing countries and isoflurane-induced general anesthesia gives rise to serious side effects. The aim of the present study was to investigate the molecular mechanism on isoflurane-induced general anesthesia.Materials and methods: The microarray data of GSE64617 dataset was downloaded from Gene Expression Omnibus (GEO) database. A total of 755 DEGs were identified using the limma package in the R programming language. Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes, and Genomes (KEGG) pathways enrichment were conducted for DEGs. A protein-protein interaction (PPI) network was constructed for DEGs and sensory perception related genes. A global miRNA-mRNA regulatory network was constructed to reveal the interactions in miRNA and mRNA in isoflurane treated samples. Degree was used to evaluate the importance of a gene in the PPI network and miRNA-mRNA regulatory network.Results and conclusions: HMBOX1, CSNK2A1, PNN, SRRM1, PRPF40A, APCNTRK1, MAPK1, hsa-miR-16-5p, hsa-miR-424-5p, hsa-miR-497-5p and hsa-miR-17-5p were selected as weighted genes. The expression changes were further vitrificated in the rat models by performing quantitative real-time PCR (qPCR) analysis. In conclusion, we find several weighted mRNAs and miRNAs involved in isoflurane induced general anesthesia through bioinformatics analysis.

Global long non-coding RNA expression in the rostral anterior cingulate cortex of depressed suicides

Long non-coding RNAs (lncRNAs) are an emerging class of regulatory RNA that may be implicated in psychiatric disorders. Here we performed RNA-sequencing in the rostral anterior cingulate cortex of 26 depressed suicides and 24 matched controls. We first performed differential lncRNA expression analysis, and then conducted Weighted Gene Co-expression Network Analysis (WGCNA) to identify co-expression modules associating with depression and suicide. We identified 23 differentially expressed lncRNAs (FDR < 0.1) as well as their differentially expressed overlapping and antisense protein-coding genes. Several of these overlapping or antisense genes were associated with interferon signaling, which is a component of the innate immune response. Using WGCNA, we identified modules of highly co-expressed genes associated with depression and suicide and found protein-coding genes highly connected to differentially expressed lncRNAs within these modules. These protein-coding genes were located distal to their associated lncRNAs and were found to be part of several GO terms enriched in the significant modules, which include: cytoskeleton organization, plasma membrane, cell adhesion, nucleus, DNA-binding, and regulation of dendrite development and morphology. Altogether, we report that lncRNAs are differentially expressed in the brains of depressed individuals who died by suicide and may represent regulators of important molecular functions and biological processes.

HMBOX1 in hepatocytes attenuates LPS/D-GalN-induced liver injury by inhibiting macrophage infiltration and activation

The HMBOX1 (Homeobox Containing 1) gene was first isolated from the human pancreatic cDNA libraries and is widely expressed in many tissues. Previously, we detected high expression of HMBOX1 in the liver, but its function was unclear. In this study, hepatocyte-specific HMBOX1 knockout mice (Hm^△hep mice) were generated and used to characterize the function of HMBOX1 in the LPS/D-GalN-induced acute liver failure model. HMBOX1-knockout exhibits exacerbated liver injury induced by LPS/D-GalN, accompanied with high levels of inflammatory cytokines both in the liver and in circulation. Further investigation demonstrated that HMBOX1 negatively regulates NF-κB signal transduction. Therefore, HMBOX1-knockout in hepatocytes promotes CCL2 expression through the activation of NF-κB signaling, which enhanced the infiltration of macrophages into the liver. In addition, the decrease of HMBOX1 in hepatocytes promotes the activation of macrophages, upregulating CD80 and MHCⅡ, as well as inflammatory factors TNF-α and IL-6. Importantly, overexpression of HMBOX1 rescued liver injury in Hm^△hep mice. These findings indicate that HMBOX1 in hepatocytes acts as a key immunosuppressive factor for inflammation and plays a critical protective role in LPS/D-GalN-induced liver injury.

Stochastic left-right neuronal asymmetry in Caenorhabditis elegans

Left-right asymmetry in the nervous system is observed across species. Defects in left-right cerebral asymmetry are linked to several neurological diseases, but the molecular mechanisms underlying brain asymmetry in vertebrates are still not very well understood. The Caenorhabditis elegans left and right amphid wing 'C' (AWC) olfactory neurons communicate through intercellular calcium signalling in a transient embryonic gap junction neural network to specify two asymmetric subtypes, AWC^OFF (default) and AWC^ON (induced), in a stochastic manner. Here, we highlight the molecular mechanisms that establish and maintain stochastic AWC asymmetry. As the components of the AWC asymmetry pathway are highly conserved, insights from the model organism C. elegans may provide a window onto how brain asymmetry develops in humans.This article is part of the themed issue 'Provocative questions in left-right asymmetry'.

Identification of high risk anaplastic gliomas by a diagnostic and prognostic signature derived from mRNA expression profiling

Anaplastic gliomas are characterized by variable clinical and genetic features, but there are few studies focusing on the substratification of anaplastic gliomas. To identify a more objective and applicable classification of anaplastic gliomas, we analyzed whole genome mRNA expression profiling of four independent datasets. Univariate Cox regression, linear risk score formula and receiver operating characteristic (ROC) curve were applied to derive a gene signature with best prognostic performance. The corresponding clinical and molecular information were further analyzed for interpretation of the different prognosis and the independence of the signature. Gene ontology (GO), Gene Set Variation Analysis (GSVA) and Gene Set Enrichment Analysis (GSEA) were performed for functional annotation of the differences. We found a three-gene signature, by applying which, the anaplastic gliomas could be divided into low risk and high risk groups. The two groups showed a high concordance with grade II and grade IV gliomas, respectively. The high risk group was more aggressive and complex. The three-gene signature showed diagnostic and prognostic value in anaplastic gliomas.

HMBOX1 interacts with MT2A to regulate autophagy and apoptosis in vascular endothelial cells

We previously found that Homeobox containing 1 (HMBOX1) was required for bone mesenchymal stem cell (BMSC) and mouse embryonic stem cell (ESC) differentiation into vascular endothelial cells (VECs). However, the function of HMBOX1 in VECs is still unknown. In this study, we found that HMBOX1 was abundantly expressed in the cytoplasm of human umbilical vascular endothelial cells (HUVECs). Knockdown of HMBOX1 induced apoptosis and inhibited autophagy. Overexpression of HMBOX1 inhibited apoptosis induced by fibroblast growth factor 2 deprivation and promoted autophagy. Metallothionein 2A (MT2A) was identified as an interaction protein with HMBOX1 by yeast two-hybrid assay, and confirmed by co-immunoprecipitation. Overexpression of HMBOX1 elevated intracellular free zinc level. Knockdown of MT2A inhibited this phenomenon. Moreover, N,N,N = ,N = -tetrakis (2-pyridylmethyl) ethylenediamine (TPEN), a zinc chelator, reversed the anti-apoptosis and pro-autophagy effects of HMBOX1. In conclusion, HMBOX1 regulated intracellular free zinc level by interacting with MT2A to inhibit apoptosis and promote autophagy in VECs.

Identification of a small molecule preventing BMSC senescence in vitro by improving intracellular homeostasis via ANXA7 and Hmbox1

ABO was discovered to be a novel anti-aging chemical in cultured BMSCs by improving intracellular homeostasis.

A chemical small molecule induces mouse embryonic stem cell differentiation into functional vascular endothelial cells via Hmbox1

Embryonic stem cells (ESCs) can differentiate to endothelial progenitor cells and vascular endothelial cells (VECs), but the mechanism is largely unknown. In this study, we synthesized 2 chiral compounds (R-ABO and S-ABO) and identified R-ABO as an effective inducer of ESC differentiation into VECs. Furthermore, we found that R-ABO induced ESC differentiation into VECs via homeobox containing 1 (Hmbox1) that acted upstream of fibroblast growth factor 2 (FGF-2). The data suggest that R-ABO is a novel tool for ESC differentiation into VECs, and Hmbox1 is a key regulator in this differentiation process. These findings provide information on a novel target and a new platform for further investigating the gene control of ESC differentiation to VECs.

miR-30c-1* promotes natural killer cell cytotoxicity against human hepatoma cells by targeting the transcription factor HMBOX1

Natural killer (NK) cells play a critical role in antitumor immunity, and the activation of NK cells is regulated by a series of NK cell receptors. Here, we show that crosslinking CD226, an important NK cell receptor, with the anti-CD226 mAb LeoA1 on NKL cells, regulated the expression of several microRNA and transmembrane tumor necrosis factor-α. Among them, miR-30c-1(*) was noticed because overexpression of miR-30c-1(*) triggered upregulation of transmembrane tumor necrosis factor-α expression and enhanced NK cell cytotoxicity against hepatoma cell lines SMMC-7721 and HepG2. Furthermore, we proved that the inhibitory transcription factor HMBOX1, which depressed the activation of NK cells, was the direct target gene of miR-30c-1(*). In conclusion, our results revealed a novel regulatory mechanism: miR-30c-1(*) promoted NK cell cytotoxicity against hepatoma cells by targeting HMBOX1.