BRSK2 is regulated by ER stress in protein level and involved in ER stress-induced apoptosis

Mycoplasma glycine cleavage system key subunit GcvH is an apoptosis inhibitor targeting host endoplasmic reticulum

Mycoplasmas are minimal but notorious bacteria that infect humans and animals. These genome-reduced organisms have evolved strategies to overcome host apoptotic defense and establish persistent infection. Here, using Mycoplasma bovis as a model, we demonstrate that mycoplasma glycine cleavage system (GCS) H protein (GcvH) targets the endoplasmic reticulum (ER) to hijack host apoptosis facilitating bacterial infection. Mechanically, GcvH interacts with the ER-resident kinase Brsk2 and stabilizes it by blocking its autophagic degradation. Brsk2 subsequently disturbs unfolded protein response (UPR) signaling, thereby inhibiting the key apoptotic molecule CHOP expression and ER-mediated intrinsic apoptotic pathway. CHOP mediates a cross-talk between ER- and mitochondria-mediated intrinsic apoptosis. The GcvH N-terminal amino acid 31-35 region is necessary for GcvH interaction with Brsk2, as well as for GcvH to exert anti-apoptotic and potentially pro-infective functions. Notably, targeting Brsk2 to dampen apoptosis may be a conserved strategy for GCS-containing mycoplasmas. Our study reveals a novel role for the conserved metabolic route protein GcvH in Mycoplasma species. It also sheds light on how genome-reduced bacteria exploit a limited number of genomic proteins to resist host cell apoptosis thereby facilitating pathogenesis.

Genome-wide DNA methylation analysis related to ALS patient progression and survival

BACKGROUND

Epigenetics contributes to the pathogenesis of amyotrophic lateral sclerosis (ALS). We aimed to characterize the DNA methylation profiles associated with clinical heterogeneity in disease progression and survival among patients.

METHODS

We included a cohort of 41 patients with sporadic ALS, with a median follow-up of 86.9 months, and 27 rigorously matched healthy controls. Blood-based genome-wide DNA methylation analysis was conducted.

RESULTS

A total of 948 progression rate-associated differentially methylated positions, 298 progression rate-associated differentially methylated regions (R-DMRs), 590 survival time-associated DMPs, and 197 survival time-associated DMRs (S-DMRs) were identified, using complementary grouping strategies. Enrichment analysis of differentially methylated genes highlighted the involvement of synapses and axons in ALS progression and survival. Clinical analysis revealed a positive correlation between the average methylation levels of the R-DMR in PRDM8 and disease progression rate (r = 0.479, p = 0.002). Conversely, there was an inverse correlation between the average methylation levels of the R-DMR in ANKRD33 and disease progression rate (r = - 0.476, p = 0.002). In addition, patients with higher methylation levels within the S-DMR of ZNF696 experienced longer survival (p = 0.016), while those with elevated methylation levels in the S-DMR of RAI1 had shorter survival (p = 0.006).

CONCLUSION

DNA methylation holds promise as a potential biomarker for tracking disease progression and predicting survival outcome and also offers targets for precision medicine.

Effects of different iodine levels on the DNA methylation of intrinsic apoptosis-associated genes and analysis of gene-environment interactions in patients with autoimmune thyroiditis

Iodine is an essential nutrient that may change the occurrence of autoimmune thyroiditis (AIT). Apoptosis and DNA methylation participate in the pathogenesis and destructive mechanism of AIT. We detected the methylation and the expression of mRNA of intrinsic apoptosis-associated genes (YWHAG, ING4, BRSK2 and GJA1) to identify the potential interactions between the levels of methylation in these genes and different levels of iodine. 176 adult patients with AIT in Shandong Province, China, were included. The MethylTargetTM assay was used to verify the levels of methylation. We used PCR to detect the mRNA levels of the candidate genes. Interactions between methylation levels of the candidate genes and iodine levels were evaluated with multiplicative and addictive interaction models and GMDR. In the AIT group, YWHAG_1 and six CpG sites and BRSK2_1 and eight CpG sites were hypermethylated, whereas ING4_1 and one CpG site were hypomethylated. A negative correlation was found between methylation levels of YWHAG and mRNA expression. The combination of iodine fortification, YWHAG_1 hypermethylation and BRSK2_1 hypermethylation was significantly associated with elevated AIT risk. A four-locus model (YWHAG_1 × ING4_1 × BRSK2_1 × iodine level) was found to be the best model of the gene-environment interactions. We identified abnormal changes in the methylation status of YWHAG, ING4 and BRSK2 in patients with AIT in different iodine levels. Iodine fortification not only affected the methylation levels of YWHAG and BRSK2 but also interacted with the methylation levels of these genes and may ultimately increase the risk of AIT.

Correlation of endoplasmic reticulum stress patterns with the immune microenvironment in hepatocellular carcinoma: a prognostic signature analysis

Backgrounds

The extended duration of endoplasmic reticulum stress (ERS) can impact the progression of hepatocellular carcinoma (HCC) and the efficacy of immunotherapies by interacting with immune cells that have infiltrated the tumor microenvironment (TME).

Methods and results

The study utilized a training cohort of 364 HCC patients with complete information from The Cancer Genome Atlas Program (TCGA) database, and a validation cohort of 231 HCC patients from the International Cancer Genome Consortium (ICGC) database. The genes related to ERS exhibiting a strong correlation with overall survival (OS) were identified using univariate Cox regression analysis. A 13-gene predictive signature was then produced through the least absolute shrinkage and selection operator (LASSO) regression approach. The data revealed that the ERS-associated gene signature effectively stratified patients into high- or low-risk groups regarding OS in both the training and validation cohorts (P < 0.0001 and P = 0.00029, respectively). Using the multivariate method, it is still an independent prognostic factor in both the training and validation cohorts (P < 0.001 and P = 0.008, respectively). Moreover, several metabolic pathways were identified to be enriched among the 13 genes in the predictive signature. When the ERS-associated gene signature was combined with the tumor-node-metastasis (TNM) stage, the ERS nomogram performed better than either the gene signature or the TNM stage alone (C-index values: 0.731, 0.729, and 0.573, respectively). Further analysis revealed that patients in the high-risk group exhibited increased infiltration of immune cells. Additionally, GP6 was downregulated in HCC tissues among these signature genes (P < 0.05), which was related to poor OS.

Conclusions

The data suggest that this novel ERS-associated gene signature could contribute to personalized cancer management for HCC. Moreover, targeting GP6 inhibition might be a potential method for HCC therapy.

Identification and Characterization of Novel Mutations in Chronic Kidney Disease (CKD) and Autosomal Dominant Polycystic Kidney Disease (ADPKD) in Saudi Subjects by Whole-Exome Sequencing

Background: Autosomal dominant polycystic kidney disease (ADPKD) is a condition usually caused by a single gene mutation and manifested by both renal and extrarenal features, eventually leading to end-stage renal disease (ESRD) by the median age of 60 years worldwide. Approximately 89% of ADPKD patients had either PKD1 or PKD2 gene mutations. The majority (85%) of the mutations are in the PKD1 gene, especially in the context of family history. Objectives: This study investigated the genetic basis and the undiscovered genes that are involved in ADPKD development among the Saudi population. Materials and Methods: In this study, 11 patients with chronic kidney disease were enrolled. The diagnosis of ADPKD was based on history and diagnostic images: CT images include enlargement of renal outlines, renal echogenicity, and presence of multiple renal cysts with dilated collecting ducts, loss of corticomedullary differentiation, and changes in GFR and serum creatinine levels. Next-generation whole-exome sequencing was conducted using the Ion Torrent PGM platform. Results: Of the 11 Saudi patients diagnosed with chronic kidney disease (CKD) and ADPKD, the most common heterozygote nonsynonymous variant in the PKD1 gene was exon15: (c.4264G > A). Two missense mutations were identified with a PKD1 (c.1758A > C and c.9774T > G), and one patient had a PKD2 mutation (c.1445T > G). Three detected variants were novel, identified at PKD1 (c.1758A > C), PKD2L2 (c.1364A > T), and TSC2 (deletion of a’a at the 3’UTR, R1680C) genes. Other variants in PKD1L1 (c.3813_381 4delinsTG) and PKD1L2 (c.404C > T) were also detected. The median age of end-stage renal disease for ADPK patients in Saudi Arabia was 30 years. Conclusion: This study reported a common variant in the PKD1 gene in Saudi patients with typical ADPKD. We also reported (to our knowledge) for the first time two novel missense variants in PKD1 and PKD2L2 genes and one indel mutation at the 3’UTR of the TSC2 gene. This study establishes that the reported mutations in the affected genes resulted in ADPKD development in the Saudi population by a median age of 30. Nevertheless, future protein−protein interaction studies to investigate the influence of these mutations on PKD1 and PKD2 functions are required. Furthermore, large-scale population-based studies to verify these findings are recommended.

Gain-of-function genetic screen of the kinome reveals BRSK2 as an inhibitor of the NRF2 transcription factor

Nuclear factor erythroid 2-related factor 2 (NFE2L2, also known as NRF2) is a transcription factor and master regulator of cellular antioxidant response. Aberrantly high NRF2-dependent transcription is recurrent in human cancer, but conversely NRF2 activity diminishes with age and in neurodegenerative and metabolic disorders. Although NRF2-activating drugs are clinically beneficial, NRF2 inhibitors do not yet exist. Here, we describe use of a gain-of-function genetic screen of the kinome to identify new druggable regulators of NRF2 signaling. We found that the under-studied protein kinase brain-specific kinase 2 (BRSK2) and the related BRSK1 kinases suppress NRF2-dependent transcription and NRF2 protein levels in an activity-dependent manner. Integrated phosphoproteomics and RNAseq studies revealed that BRSK2 drives 5'-AMP-activated protein kinase α2 (AMPK) signaling and suppresses the mTOR pathway. As a result, BRSK2 kinase activation suppresses ribosome-RNA complexes, global protein synthesis and NRF2 protein levels. Collectively, our data illuminate the BRSK2 and BRSK1 kinases, in part by functionally connecting them to NRF2 signaling and mTOR. This signaling axis might prove useful for therapeutically targeting NRF2 in human disease.This article has an associated First Person interview with the first author of the paper.

Increased Expression of Lipid Metabolism Genes in Early Stages of Wooden Breast Links Myopathy of Broilers to Metabolic Syndrome in Humans

Wooden breast is a muscle disorder affecting modern commercial broiler chickens that causes a palpably firm pectoralis major muscle and severe reduction in meat quality. Most studies have focused on advanced stages of wooden breast apparent at market age, resulting in limited insights into the etiology and early pathogenesis of the myopathy. Therefore, the objective of this study was to identify early molecular signals in the wooden breast transcriptional cascade by performing gene expression analysis on the pectoralis major muscle of two-week-old birds that may later exhibit the wooden breast phenotype by market age at 7 weeks. Biopsy samples of the left pectoralis major muscle were collected from 101 birds at 14 days of age. Birds were subsequently raised to 7 weeks of age to allow sample selection based on the wooden breast phenotype at market age. RNA-sequencing was performed on 5 unaffected and 8 affected female chicken samples, selected based on wooden breast scores (0 to 4) assigned at necropsy where affected birds had scores of 2 or 3 (mildly or moderately affected) while unaffected birds had scores of 0 (no apparent gross lesions). Differential expression analysis identified 60 genes found to be significant at an FDR-adjusted p-value of 0.05. Of these, 26 were previously demonstrated to exhibit altered expression or genetic polymorphisms related to glucose tolerance or diabetes mellitus in mammals. Additionally, 9 genes have functions directly related to lipid metabolism and 11 genes are associated with adiposity traits such as intramuscular fat and body mass index. This study suggests that wooden breast disease is first and foremost a metabolic disorder characterized primarily by ectopic lipid accumulation in the pectoralis major.

Selection of growth-related genes and dominant genotypes in transgenic Yellow River carp Cyprinus carpio L

Transgenic Yellow River carp is characterized by rapid growth rate and high feed-conversion efficiency and exhibits a great application prospect. However, there is still a significant separation of growth traits in the transgenic Yellow River carp family; as such, growth-related genotypes must be screened for molecular marker-assisted selection. In this study, 23 growth-related candidate genes containing 48 SNP markers were screened through bulked segregant analysis (BSA) among transgenic Yellow River carp family members showing significant separation of growth traits. Then, two growth-related genes (Nos. 17 and 14 genes) were identified through combined genome-wide association study (GWAS) of candidate genes and validation of the full-sibling family approach. Nos. 17 and 14 genes encode BR serine/threonine-protein kinase 2 (BRSK2) and eukaryotic translation-initiation factor 2-alpha kinase 3 (Eif2ak3), respectively. The average body weight of three subgroups carrying the genotypes 17GG, 17GG + 14CC, and 17GG + 14TT of these two genes increased by 27.96, 38.28, and 33.72%, respectively, compared with the controls. The proportion of individuals with body weight > 500 g in these subgroups increased by 19.22, 26.82, and 30.92%, respectively. The results showed that appropriate genotype carriers can be selected from the progeny population through BSA sequencing combined with simplified GWAS analysis. Hence, basic population for breeding can be constructed and transgenic Yellow River carp strains with stable production performance and uniform phenotypic properties can be bred.

Identification of Upstream Kinases by Fluorescence Complementation Mass Spectrometry

Protein kinases and their substrates comprise extensive signaling networks that regulate many diverse cellular functions. However, methods and techniques to systematically identify kinases directly responsible for specific phosphorylation events have remained elusive. Here we describe a novel proteomic strategy termed fluorescence complementation mass spectrometry (FCMS) to identify kinase-substrate pairs in high throughput. The FCMS strategy employs a specific substrate and a kinase library, both of which are fused with fluorescence complemented protein fragments. Transient and weak kinase-substrate interactions in living cells are stabilized by the association of fluorescence protein fragments. These kinase-substrate pairs are then isolated with high specificity and are identified and quantified by LC-MS. FCMS was applied to the identification of both known and novel kinases of the transcription factor, cAMP response element-binding protein (CREB). Novel CREB kinases were validated by in vitro kinase assays, and the phosphorylation sites were unambiguously located. These results uncovered possible new roles for CREB in multiple important signaling pathways and demonstrated the great potential of this new proteomic strategy.

Mitomycin C induces apoptosis in human epidural scar fibroblasts after surgical decompression for spinal cord injury

Numerous studies have shown that topical application of mitomycin C after surgical decompression effectively reduces scar adhesion. However, the underlying mechanisms remain unclear. In this study, we investigated the effect of mitomycin C on the proliferation and apoptosis of human epidural scar fibroblasts. Human epidural scar fibroblasts were treated with various concentrations of mitomycin C (1, 5, 10, 20, 40 μg/mL) for 12, 24 and 48 hours. Mitomycin C suppressed the growth of these cells in a dose- and time-dependent manner. Mitomycin C upregulated the expression levels of Fas, DR4, DR5, cleaved caspase-8/9, Bax, Bim and cleaved caspase-3 proteins, and it downregulated Bcl-2 and Bcl-xL expression. In addition, inhibitors of caspase-8 and caspase-9 (Z-IETD-FMK and Z-LEHD-FMK, respectively) did not fully inhibit mitomycin C-induced apoptosis. Furthermore, mitomycin C induced endoplasmic reticulum stress by increasing the expression of glucose-regulated protein 78, CAAT/enhancer-binding protein homologous protein (CHOP) and caspase-4 in a dose-dependent manner. Salubrinal significantly inhibited the mitomycin C-induced cell viability loss and apoptosis, and these effects were accompanied by a reduction in CHOP expression. Our results support the hypothesis that mitomycin C induces human epidural scar fibroblast apoptosis, at least in part, via the endoplasmic reticulum stress pathway.

Epigenetic Variability across Human Populations: A Focus on DNA Methylation Profiles of the KRTCAP3, MAD1L1 and BRSK2 Genes

Natural epigenetic diversity has been suggested as a key mechanism in microevolutionary processes due to its capability to create phenotypic variability within individuals and populations. It constitutes an important reservoir of variation potentially useful for rapid adaptation in response to environmental stimuli. The analysis of population epigenetic structure represents a possible tool to study human adaptation and to identify external factors that are able to naturally shape human DNA methylation variability. The aim of this study is to investigate the dynamics that create epigenetic diversity between and within different human groups. To this end, we first used publicly available epigenome-wide data to explore population-specific DNA methylation changes that occur at macro-geographic scales. Results from this analysis suggest that nutrients, UVA exposure and pathogens load might represent the main environmental factors able to shape DNA methylation profiles. Then, we evaluated DNA methylation of candidate genes (KRTCAP3, MAD1L1, and BRSK2), emerged from the previous analysis, in individuals belonging to different populations from Morocco, Nigeria, Philippines, China, and Italy, but living in the same Italian city. DNA methylation of the BRSK2 gene is significantly different between Moroccans and Nigerians (pairwise t-test: CpG 6 P-value = 5.2*10 (-) (3); CpG 9 P-value = 2.6*10 (-) (3); CpG 10 P-value = 3.1*10 (-) (3); CpG 11 P-value = 2.8*10 (-) (3)). Comprehensively, these results suggest that DNA methylation diversity is a source of variability in human groups at macro and microgeographical scales and that population demographic and adaptive histories, as well as the individual ancestry, actually influence DNA methylation profiles.

Metabolomic Quantitative Trait Loci (mQTL) Mapping Implicates the Ubiquitin Proteasome System in Cardiovascular Disease Pathogenesis

Levels of certain circulating short-chain dicarboxylacylcarnitine (SCDA), long-chain dicarboxylacylcarnitine (LCDA) and medium chain acylcarnitine (MCA) metabolites are heritable and predict cardiovascular disease (CVD) events. Little is known about the biological pathways that influence levels of most of these metabolites. Here, we analyzed genetics, epigenetics, and transcriptomics with metabolomics in samples from a large CVD cohort to identify novel genetic markers for CVD and to better understand the role of metabolites in CVD pathogenesis. Using genomewide association in the CATHGEN cohort (N = 1490), we observed associations of several metabolites with genetic loci. Our strongest findings were for SCDA metabolite levels with variants in genes that regulate components of endoplasmic reticulum (ER) stress (USP3, HERC1, STIM1, SEL1L, FBXO25, SUGT1) These findings were validated in a second cohort of CATHGEN subjects (N = 2022, combined p = 8.4x10^-6–2.3x10^-10). Importantly, variants in these genes independently predicted CVD events. Association of genomewide methylation profiles with SCDA metabolites identified two ER stress genes as differentially methylated (BRSK2 and HOOK2). Expression quantitative trait loci (eQTL) pathway analyses driven by gene variants and SCDA metabolites corroborated perturbations in ER stress and highlighted the ubiquitin proteasome system (UPS) arm. Moreover, culture of human kidney cells in the presence of levels of fatty acids found in individuals with cardiometabolic disease, induced accumulation of SCDA metabolites in parallel with increases in the ER stress marker BiP. Thus, our integrative strategy implicates the UPS arm of the ER stress pathway in CVD pathogenesis, and identifies novel genetic loci associated with CVD event risk.

Cardiovascular disease is a strongly heritable trait. Despite application of the latest genomic technologies, the genetic architecture of disease risk remains poorly defined, and mechanisms underlying this susceptibility are incompletely understood. In this study, we performed genome-wide mapping of heart disease-related metabolites measured in the blood as the genetic traits of interest (instead of the disease itself), in a large cohort of 3512 patients at risk of heart disease from the CATHGEN study. Our goal was to discover new cardiovascular disease genes and thereby mechanisms of disease pathogenesis by understanding the genes that regulate levels of these metabolites. These analyses identified novel genetic variants associated with metabolite levels and with cardiovascular disease itself. Importantly, by utilizing an unbiased systems-based approach integrating genetics, gene expression, epigenetics and metabolomics, we uncovered a novel pathway of heart disease pathogenesis, that of endoplasmic reticulum (ER) stress, represented by elevated levels of circulating short-chain dicarboxylacylcarnitine (SCDA) metabolites.

An information theoretic method to identify combinations of genomic alterations that promote glioblastoma

Tumors are the result of accumulated genomic alterations that cooperate synergistically to produce uncontrollable cell growth. Although identifying recurrent alterations among large collections of tumors provides a way to pinpoint genes that endow a selective advantage in oncogenesis and progression, it fails to address the genetic interactions behind this selection process. A non-random pattern of co-mutated genes is evidence for selective forces acting on tumor cells that harbor combinations of these genetic alterations. Although existing methods have successfully identified mutually exclusive gene sets, no current method can systematically discover more general genetic relationships. We develop Genomic Alteration Modules using Total Correlation (GAMToC), an information theoretic framework that integrates copy number and mutation data to identify gene modules with any non-random pattern of joint alteration. Additionally, we present the Seed-GAMToC procedure, which uncovers the mutational context of any putative cancer gene. The software is publicly available. Applied to glioblastoma multiforme samples, GAMToC results show distinct subsets of co-occurring mutations, suggesting distinct mutational routes to cancer and providing new insight into mutations associated with proneural, proneural/G-CIMP, and classical types of the disease. The results recapitulate known relationships such as mutual exclusive mutations, place these alterations in the context of other mutations, and find more complex relationships such as conditional mutual exclusivity.

BRSK2 is a valosin-containing protein (VCP)-interacting protein that affects VCP functioning in endoplasmic reticulum-associated degradation

Endoplasmic reticulum-associated protein degradation (ERAD) removes improperly-folded proteins from the ER membrane into the cytosol where they undergo proteasomal degradation. Valosin-containing protein (VCP)/p97 mediates in the extraction of ERAD substrates from the ER. BRSK2 (also known as SAD-A), a serine/threonine kinase of the AMP-activated protein kinase family affected VCP/p97 activity in ERAD. In addition, BRSK2 interacted with VCP/p97 via three of the four functional domains of VCP/p97. Immunofluorescence demonstrated that BRSK2 and VCP/p97 were co-localized and also that knockdown of endogenous BRSK2 induced increased levels of CD3δ, a substrate in ERAD for VCP/p97. Thus, BRSK2 might affect the activity of VCP/p97 in ERAD.