Voluntary running exercise modifies astrocytic population and features in the peri-infarct cortex

Rehabilitative exercise following a brain stroke has beneficial effects on the morphological plasticity of neurons. Particularly, voluntary running exercise after focal cerebral ischemia promotes functional recovery and ameliorates ischemia-induced dendritic spine loss in the peri-infarct motor cortex layer 5. Moreover, neuronal morphology is affected by changes in the perineuronal environment. Glial cells, whose phenotypes may be altered by exercise, are known to play a pivotal role in the formation of this perineuronal environment. Herein, we investigated the effects of voluntary running exercise on glial cells after middle cerebral artery occlusion. Voluntary running exercise increased the population of glial fibrillary acidic protein-positive astrocytes born between post-operative days (POD) 0 and 3 on POD15 in the peri-infarct cortex. After exercise, transcriptomic analysis of post-ischemic astrocytes revealed 10 upregulated and 70 downregulated genes. Furthermore, gene ontology analysis showed that the 70 downregulated genes were significantly associated with neuronal morphology. In addition, exercise reduced the number of astrocytes expressing lipocalin 2, a regulator of dendritic spine density, on POD15. Our results suggest that exercise modifies the composition of astrocytic population and their phenotype.

Peripheral Blood Gene Expression Profile of Infants with Atopic Dermatitis

To enhance the understanding of molecular mechanisms and mine previously unidentified biomarkers of pediatric atopic dermatitis, PBMC gene expression profiles were generated by RNA sequencing in infants with atopic dermatitis and age-matched controls. A total of 178 significantly differentially expressed genes (DEGs) (115 upregulations and 63 downregulations) were seen, compared with those in healthy controls. The DEGs identified included IL1β, TNF, TREM1, IL18R1, and IL18RAP. DEGs were validated by real-time RT- qPCR in a larger number of samples from PBMCs of infants with atopic dermatitis aged <12 months. Using the DAVID (Database for Annotation, Visualization and Integrated Discovery) database, functional and pathway enrichment analyses of DEGs were performed. Gene ontology enrichment analysis showed that DEGs were associated with immune responses, inflammatory responses, regulation of immune responses, and platelet activation. Pathway analysis indicated that DEGs were enriched in cytokine‒cytokine receptor interaction, immunoregulatory interactions between lymphoid and nonlymphoid cells, hematopoietic cell lineage, phosphoinositide 3-kinase‒protein kinase B signaling pathway, NK cell‒mediated cytotoxicity, and platelet activation. Furthermore, the protein‒protein interaction network was predicted using the STRING (Search Tool for the Retrieval of Interacting Genes/Proteins) database and visualized with Cytoscape software. Finally, on the basis of the protein‒protein interaction network, 18 hub genes were selected, and two significant modules were obtained. In conclusion, this study sheds light on the molecular mechanisms of pediatric atopic dermatitis and may provide diagnostic biomarkers and therapeutic targets.

Fibrosis resolution in the mouse liver: Role of Mmp12 and potential role of calpain 1/2

Although most work has focused on resolution of collagen ECM, fibrosis resolution involves changes to several ECM proteins. The purpose of the current study was twofold: 1) to examine the role of MMP12 and elastin; and 2) to investigate the changes in degraded proteins in plasma (i.e., the "degradome") in a preclinical model of fibrosis resolution. Fibrosis was induced by 4 weeks carbon tetrachloride (CCl4) exposure, and recovery was monitored for an additional 4 weeks. Some mice were treated with daily MMP12 inhibitor (MMP408) during the resolution phase. Liver injury and fibrosis was monitored by clinical chemistry, histology and gene expression. The release of degraded ECM peptides in the plasma was analyzed using by 1D-LC-MS/MS, coupled with PEAKS Studio (v10) peptide identification. Hepatic fibrosis and liver injury rapidly resolved in this mouse model. However, some collagen fibrils were still present 28d after cessation of CCl4. Despite this persistent collagen presence, expression of canonical markers of fibrosis were also normalized. The inhibition of MMP12 dramatically delayed fibrosis resolution under these conditions. LC-MS/MS analysis identified that several proteins were being degraded even at late stages of fibrosis resolution. Calpains 1/2 were identified as potential new proteases involved in fibrosis resolution. CONCLUSION. The results of this study indicate that remodeling of the liver during recovery from fibrosis is a complex and highly coordinated process that extends well beyond the degradation of the collagenous scar. These results also indicate that analysis of the plasma degradome may yield new insight into the mechanisms of fibrosis recovery, and by extension, new "theragnostic" targets. Lastly, a novel potential role for calpain activation in the degradation and turnover of proteins was identified.

Effects of galactosyltransferase on EPS biosynthesis and freeze-drying resistance of Lactobacillus acidophilus NCFM

Galactosyltransferase (GalT) is an important enzyme in synthesizing exopolysaccharide (EPS), the major polymer of biofilms protecting cells from severe conditions. However, the contribution to, and regulatory mechanism of GalT, in stressor resistance are still unclear. Herein, we successfully overexpressed GalT in Lactobacillus acidophilus NCFM by genetic engineering. The GalT activity and freeze-drying survival rate of the recombinant strain were significantly enhanced. The EPS yield also increased by 17.8%, indicating a positive relationship between freeze-drying resistance and EPS. RNA-Seq revealed that GalT could regulate the flux of the membrane transport system, pivotal sugar-related metabolic pathways, and promote quorum sensing to facilitate EPS biosynthesis, which enhanced freeze-drying resistance. The findings concretely prove that the mechanism of GalT regulating EPS biosynthesis plays an important role in protecting lactic acid bacteria from freeze-drying stress.

Unravelling the neuroprotective mechanisms of carotenes in differentiated human neural cells: Biochemical and proteomic approaches

•

Total mixed carotenes (TMC) protect differentiated human neural cells against 6-hydroxydopamine-induced toxicity.

•

TMC elevated the antioxidant enzymes activities and suppressed generation of reactive oxygen species.

•

TMC augmented the dopamine and tyrosine hydroxylase levels.

•

TMC exerted differential protein expression in human neural cells.

Carotenoids, fat-soluble pigments found ubiquitously in plants and fruits, have been reported to exert significant neuroprotective effects against free radicals. However, the neuroprotective effects of total mixed carotenes complex (TMC) derived from virgin crude palm oil have not been studied extensively. Therefore, the present study was designed to establish the neuroprotective role of TMC on differentiated human neural cells against 6-hydroxydopamine (6-OHDA)-induced cytotoxicity. The human neural cells were differentiated using retinoic acid for six days. Then, the differentiated neural cells were pre-treated for 24 hr with TMC before exposure to 6-OHDA. TMC pre-treated neurons showed significant alleviation of 6-OHDA-induced cytotoxicity as evidenced by enhanced activity of the superoxide dismutase (SOD) and catalase (CAT) enzymes. Furthermore, TMC elevated the levels of intra-neuronal dopamine and tyrosine hydroxylase (TH) in differentiated neural cells. The 6-OHDA induced overexpression of α-synuclein was significantly hindered in neural cells pre-treated with TMC. In proteomic analysis, TMC altered the expression of ribosomal proteins, α/β isotypes of tubulins, protein disulphide isomerases (PDI) and heat shock proteins (HSP) in differentiated human neural cells. The natural palm phytonutrient TMC is a potent antioxidant with significant neuroprotective effects against free radical-induced oxidative stress.

MSPJ: Discovering potential biomarkers in small gene expression datasets via ensemble learning

In transcriptomics, differentially expressed genes (DEGs) provide fine-grained phenotypic resolution for comparisons between groups and insights into molecular mechanisms underlying the pathogenesis of complex diseases or phenotypes. The robust detection of DEGs from large datasets is well-established. However, owing to various limitations (e.g., the low availability of samples for some diseases or limited research funding), small sample size is frequently used in experiments. Therefore, methods to screen reliable and stable features are urgently needed for analyses with limited sample size. In this study, MSPJ, a new machine learning approach for identifying DEGs was proposed to mitigate the reduced power and improve the stability of DEG identification in small gene expression datasets. This ensemble learning-based method consists of three algorithms: an improved multiple random sampling with meta-analysis, SVM-RFE (support vector machines-recursive feature elimination), and permutation test. MSPJ was compared with ten classical methods by 94 simulated datasets and large-scale benchmarking with 165 real datasets. The results showed that, among these methods MSPJ had the best performance in most small gene expression datasets, especially those with sample size below 30. In summary, the MSPJ method enables effective feature selection for robust DEG identification in small transcriptome datasets and is expected to expand research on the molecular mechanisms underlying complex diseases or phenotypes.

Comprehensive analysis of gene expression profiles to identify differential prognostic factors of primary and metastatic breast cancer

Breast cancer accounts for nearly half of all cancer-related deaths in women worldwide. However, the molecular mechanisms that lead to tumour development and progression remain poorly understood and there is a need to identify candidate genes associated with primary and metastatic breast cancer progression and prognosis. In this study, candidate genes associated with prognosis of primary and metastatic breast cancer were explored through a novel bioinformatics approach. Primary and metastatic breast cancer tissues and adjacent normal breast tissues were evaluated to identify biomarkers characteristic of primary and metastatic breast cancer. The Cancer Genome Atlas-breast invasive carcinoma (TCGA-BRCA) dataset (ID: HS-01619) was downloaded using the mRNASeq platform. Genevestigator 8.3.2 was used to analyse TCGA-BRCA gene expression profiles between the sample groups and identify the differentially-expressed genes (DEGs) in each group. For each group, Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses were used to determine the function of DEGs. Networks of protein-protein interactions were constructed to identify the top hub genes with the highest degree of interaction. Additionally, the top hub genes were validated based on overall survival and immunohistochemistry using The Human Protein Atlas. Of the top 20 hub genes identified, four (KRT14, KIT, RAD51, and TTK) were considered as prognostic risk factors based on overall survival. KRT14 and KIT expression levels were upregulated while those of RAD51 and TTK were downregulated in patients with breast cancer. The four proposed candidate hub genes might aid in further understanding the molecular changes that distinguish primary breast tumours from metastatic tumours as well as help in developing novel therapeutics. Furthermore, they may serve as effective prognostic risk markers based on the strong correlation between their expression and patient overall survival.

Research on the mechanism of berberine in the treatment of COVID-19 pneumonia pulmonary fibrosis using network pharmacology and molecular docking

Purpose Pulmonary fibrosis caused by COVID-19 pneumonia is a serious complication of COVID-19 infection, there is a lack of effective treatment methods clinically. This article explored the mechanism of action of berberine in the treatment of COVID-19 (Corona Virus Disease 2019, COVID-19) pneumonia pulmonary fibrosis with the help of the network pharmacology and molecular docking. Methods We predicted the role of berberine protein targets with the Pharmmapper database and the 3D structure of berberine in the Pubchem database. And GeneCards database was used in order to search disease target genes and screen common target genes. Then we used STRING web to construct PPI interaction network of common target protein. The common target genes were analyzed by GO and KEGG by DAVID database. The disease-core target gene-drug network was established and molecular docking was used for prediction. We also analyzed the binding free energy and simulates molecular dynamics of complexes. Results Berberine had 250 gene targets, COVID-19 pneumonia pulmonary fibrosis had 191 gene targets, the intersection of which was 23 in common gene targets. Molecular docking showed that berberine was associated with CCl2, IL-6, STAT3 and TNF-α. GO and KEGG analysis reveals that berberine mainly plays a vital role by the signaling pathways of influenza, inflammation and immune response. Conclusion Berberine acts on TNF-α, STAT3, IL-6, CCL2 and other targets to inhibit inflammation and the activation of fibrocytes to achieve the purpose of treating COVID-19 pneumonia pulmonary fibrosis.

Discovery of tumor immune infiltration-related snoRNAs for predicting tumor immune microenvironment status and prognosis in lung adenocarcinoma

Lung adenocarcinoma (LUAD) has a high mortality rate and is difficult to diagnose and treat in its early stage. Previous studies have demonstrated that small nucleolar RNAs (snoRNAs) play a critical role in tumor immune infiltration and the development of a variety of solid tumors. However, there have been no studies on the correlation between tumor-infiltrating immune-related snoRNAs (TIISRs) and LUAD. In this study, we filtered six immune-related snoRNAs based on the tissue specificity index (TSI) and expression profile of all snoRNAs between all LUAD cell lines from the Cancer Cell Line Encyclopedia and 21 types of immune cells from the Gene Expression Omnibus database. Further, we performed real-time quantitative polymerase chain reaction (RT-qPCR) to validate the expression status of these snoRNAs on peripheral blood mononuclear cells (PBMCs) and lung cancer cell lines. Next, we developed a TIISR signature based on the expression profiles of snoRNAs from 479 LUAD patients filtered by the random survival forest algorithm. We then analyzed the value of this TIISR signature (TIISR risk score) for assessing tumor immune infiltration, immune checkpoint inhibitor (ICI) treatment response, and the prognosis of LUAD between groups with high and low TIISR risk score. Further, we found that the TIISR risk score groups showed significant differences in biological characteristics and that the risk score could be used to assess the level of tumor immune cell infiltration, thereby predicting prognosis and responsiveness to immunotherapy in LUAD patients.

Proteomics and metabolic phenotyping define principal roles for the aryl hydrocarbon receptor in mouse liver

Dioxin-like molecules have been associated with endocrine disruption and liver disease. To better understand aryl hydrocarbon receptor (AHR) biology, metabolic phenotyping and liver proteomics were performed in mice following ligand-activation or whole-body genetic ablation of this receptor. Male wild type (WT) and Ahr^–/– mice (Taconic) were fed a control diet and exposed to 3,3′,4,4′,5-pentachlorobiphenyl (PCB126) (61 nmol/kg by gavage) or vehicle for two weeks. PCB126 increased expression of canonical AHR targets (Cyp1a1 and Cyp1a2) in WT but not Ahr^–/–. Knockouts had increased adiposity with decreased glucose tolerance; smaller livers with increased steatosis and perilipin-2; and paradoxically decreased blood lipids. PCB126 was associated with increased hepatic triglycerides in Ahr^–/–. The liver proteome was impacted more so by Ahr^–/– genotype than ligand-activation, but top gene ontology (GO) processes were similar. The PCB126-associated liver proteome was Ahr-dependent. Ahr principally regulated liver metabolism (e.g., lipids, xenobiotics, organic acids) and bioenergetics, but it also impacted liver endocrine response (e.g., the insulin receptor) and function, including the production of steroids, hepatokines, and pheromone binding proteins. These effects could have been indirectly mediated by interacting transcription factors or microRNAs. The biologic roles of the AHR and its ligands warrant more research in liver metabolic health and disease.

Abnormal Flow Conditions Promote Endocardial Fibroelastosis Via Endothelial-to-Mesenchymal Transition, Which Is Responsive to Losartan Treatment

•

EFE is a congenital cardiac pathology contributing to increased morbidity and mortality. The pathologic triggers of EFE remain to be characterized.

•

To determine whether abnormal flow promotes EFE development, we used in vivo neonatal rodent surgical models and an in vitro model using human primary endocardial cells

•

We established novel surgical model with flow profiles seen in patients that develop EFE. Static and turbulent flow conditions promoted EFE development in neonatal rodent hearts.

•

Losartan treatment is shown to significantly ameliorate EFE progression and decreases mRNA and protein expression of EndoMT markers in neonatal rodent hearts.

•

RNAseq analysis of human endocardial cells subjected to different flow conditions show that normal flow suppresses gene expression critical for mesenchymal differentiation and Notch signaling.

Endocardial fibroelastosis (EFE) is defined by fibrotic tissue on the endocardium and forms partly through aberrant endothelial-to-mesenchymal transition. However, the pathologic triggers are still unknown. In this study, we showed that abnormal flow induces EFE partly through endothelial-to-mesenchymal transition in a rodent model, and that losartan can abrogate EFE development. Furthermore, we translated our findings to human endocardial endothelial cells, and showed that laminar flow promotes the suppression of genes associated with mesenchymal differentiation. These findings emphasize the role of flow in promoting EFE in endocardial endothelial cells and provide a novel potential therapy to treat this highly morbid condition.

3D disorganization and rearrangement of genome provide insights into pathogenesis of NAFLD by integrated Hi-C, Nanopore, and RNA sequencing

The three-dimensional (3D) conformation of chromatin is integral to the precise regulation of gene expression. The 3D genome and genomic variations in non-alcoholic fatty liver disease (NAFLD) are largely unknown, despite their key roles in cellular function and physiological processes. High-throughput chromosome conformation capture (Hi-C), Nanopore sequencing, and RNA-sequencing (RNA-seq) assays were performed on the liver of normal and NAFLD mice. A high-resolution 3D chromatin interaction map was generated to examine different 3D genome hierarchies including A/B compartments, topologically associated domains (TADs), and chromatin loops by Hi-C, and whole genome sequencing identifying structural variations (SVs) and copy number variations (CNVs) by Nanopore sequencing. We identified variations in thousands of regions across the genome with respect to 3D chromatin organization and genomic rearrangements, between normal and NAFLD mice, and revealed gene dysregulation frequently accompanied by these variations. Candidate target genes were identified in NAFLD, impacted by genetic rearrangements and spatial organization disruption. Our data provide a high-resolution 3D genome interaction resource for NAFLD investigations, revealed the relationship among genetic rearrangements, spatial organization disruption, and gene regulation, and identified candidate genes associated with these variations implicated in the pathogenesis of NAFLD. The newly findings offer insights into novel mechanisms of NAFLD pathogenesis and can provide a new conceptual framework for NAFLD therapy.

The molecular feature of macrophages in tumor immune microenvironment of glioma patients

Background

Gliomas are one of the most common types of primary tumors in central nervous system. Previous studies have found that macrophages actively participate in tumor growth.

Methods

Weighted gene co-expression network analysis was used to identify meaningful macrophage-related gene genes for clustering. Pamr, SVM, and neural network were applied for validating clustering results. Somatic mutation and methylation were used for defining the features of identified clusters. Differentially expressed genes (DEGs) between the stratified groups after performing elastic regression and principal component analyses were used for the construction of MScores. The expression of macrophage-specific genes were evaluated in tumor microenvironment based on single cell sequencing analysis. A total of 2365 samples from 15 glioma datasets and 5842 pan-cancer samples were used for external validation of MScore.

Results

Macrophages were identified to be negatively associated with the survival of glioma patients. Twenty-six macrophage-specific DEGs obtained by elastic regression and PCA were highly expressed in macrophages at single-cell level. The prognostic value of MScores in glioma was validated by the active proinflammatory and metabolic profile of infiltrating microenvironment and response to immunotherapies of samples with this signature. MScores managed to stratify patient survival probabilities in 15 external glioma datasets and pan-cancer datasets, which predicted worse survival outcome. Sequencing data and immunohistochemistry of Xiangya glioma cohort confirmed the prognostic value of MScores. A prognostic model based on MScores demonstrated high accuracy rate.

Conclusion

Our findings strongly support a modulatory role of macrophages, especially M2 macrophages in glioma progression and warrants further experimental studies.

Molecular Signature of HFpEF: Systems Biology in a Cardiac-Centric Large Animal Model

In this study the authors used systems biology to define progressive changes in metabolism and transcription in a large animal model of heart failure with preserved ejection fraction (HFpEF). Transcriptomic analysis of cardiac tissue, 1-month post-banding, revealed loss of electron transport chain components, and this was supported by changes in metabolism and mitochondrial function, altogether signifying alterations in oxidative metabolism. Established HFpEF, 4 months post-banding, resulted in changes in intermediary metabolism with normalized mitochondrial function. Mitochondrial dysfunction and energetic deficiencies were noted in skeletal muscle at early and late phases of disease, suggesting cardiac-derived signaling contributes to peripheral tissue maladaptation in HFpEF. Collectively, these results provide insights into the cellular biology underlying HFpEF progression.

Identification of subtypes correlated with tumor immunity and immunotherapy in cutaneous melanoma

Because immune checkpoint inhibitors (ICIs) are effective for a subset of melanoma patients, identification of melanoma subtypes responsive to ICIs is crucial. We performed clustering analyses to identify immune subtypes of melanoma based on the enrichment levels of 28 immune cells using transcriptome datasets for six melanoma cohorts, including four cohorts not treated with ICIs and two cohorts treated with ICIs. We identified three immune subtypes (Im-H, Im-M, and Im-L), reproducible in these cohorts. Im-H displayed strong immune signatures, low stemness and proliferation potential, genomic stability, high immunotherapy response rate, and favorable prognosis. Im-L showed weak immune signatures, high stemness and proliferation potential, genomic instability, low immunotherapy response rate, and unfavorable prognosis. The pathways highly enriched in Im-H included immune, MAPK, apoptosis, calcium, VEGF, cell adhesion molecules, focal adhesion, gap junction, and PPAR. The pathways highly enriched in Im-L included Hippo, cell cycle, and ErbB. Copy number alterations correlated inversely with immune signatures in melanoma, while tumor mutation burden showed no significant correlation. The molecular features correlated with favorable immunotherapy response included immune-promoting signatures and pathways of PPAR, MAPK, VEGF, calcium, and glycolysis/gluconeogenesis. Our data recapture the immunological heterogeneity in melanoma and provide clinical implications for the immunotherapy of melanoma.

Analysis of differentially expressed genes responsible for the suppressive effect of anisomycin on cell proliferation of DLD-1 cells

Anisomycin is used as a chemical compound that possesses c-Jun N-terminal kinase (JNK)-activating effects. Recently, the potent anti-tumor effects of anisomycin have received much attention. In addition to its JNK-activating effects, anisomycin has been reported to affect gene expression in osteosarcoma, leukemia, hepatocellular carcinoma, ovarian cancer and other cancers. We previously demonstrated that anisomycin induced the degradation of transcription factor GATA-6 in DLD-1 cells (a colorectal cancer cell line) and inhibited their proliferation. However, the details of the gene network involved in the process remain unclear. In this study, we conducted an RNA-seq analysis of differentially expressed genes (DEGs) in anisomycin-treated DLD-1 cells to identify the molecular process of growth-suppressive genes. We found that LAMB3, which regulates cell adhesion and migration, and NFKB2 were down-regulated by anisomycin. In addition, the mRNA expression of several tumor suppressor genes (ATF3, ERRFI1, KLF6, and AKAP12) was transiently enhanced at 3 h after anisomycin treatment. These results suggest that anisomycin blocks a PI3K/Akt-signaling cascade to lead to the suppression of cell growth.

Single-cell RNA sequencing of human liver reveals hepatic stellate cell heterogeneity

Background & Aims

The multiple vital functions of the human liver are performed by highly specialised parenchymal and non-parenchymal cells organised in complex collaborative sinusoidal units. Although crucial for homeostasis, the cellular make-up of the human liver remains to be fully elucidated. Here, single-cell RNA-sequencing was used to unravel the heterogeneity of human liver cells, in particular of hepatocytes (HEPs) and hepatic stellate cells (HSCs).

Method

The transcriptome of ~25,000 freshly isolated human liver cells was profiled using droplet-based RNA-sequencing. Recently published data sets and RNA in situ hybridisation were integrated to validate and locate newly identified cell populations.

Results

In total, 22 cell populations were annotated that reflected the heterogeneity of human parenchymal and non-parenchymal liver cells. More than 20,000 HEPs were ordered along the portocentral axis to confirm known, and reveal previously undescribed, zonated liver functions. The existence of 2 subpopulations of human HSCs with unique gene expression signatures and distinct intralobular localisation was revealed (i.e. portal and central vein-concentrated GPC3⁺ HSCs and perisinusoidally located DBH⁺ HSCs). In particular, these data suggest that, although both subpopulations collaborate in the production and organisation of extracellular matrix, GPC3⁺ HSCs specifically express genes involved in the metabolism of glycosaminoglycans, whereas DBH⁺ HSCs display a gene signature that is reminiscent of antigen-presenting cells.

Conclusions

This study highlights metabolic zonation as a key determinant of HEP transcriptomic heterogeneity and, for the first time, outlines the existence of heterogeneous HSC subpopulations in the human liver. These findings call for further research on the functional implications of liver cell heterogeneity in health and disease.

Lay summary

This study resolves the cellular landscape of the human liver in an unbiased manner and at high resolution to provide new insights into human liver cell biology. The results highlight the physiological heterogeneity of human hepatic stellate cells.

•

A cell atlas from the near-native transcriptome of >25,000 human liver cells is presented.

•

Hepatocytes were ordered along the portocentral axis to reveal previously undescribed gene expression patterns and zonated liver functions.

•

Two subpopulations of human hepatic stellate cells (HSCs) are reported, characterised by different spatial distribution in the native tissue.

•

Characteristic gene signatures of HSC subpopulations are suggestive of far-reaching functional differences.

Integrative approaches for analysis of mRNA and microRNA high-throughput data

•

Review on tools and databases linking miRNA and its mRNA targetome.

•

Databases show little overlap in miRNA targetome predictions suggesting strong contextual effects.

•

Deconvolution and deep learning approaches are promising new approaches to improve miRNA targetome predictions.

Advanced sequencing technologies such as RNASeq provide the means for production of massive amounts of data, including transcriptome-wide expression levels of coding RNAs (mRNAs) and non-coding RNAs such as miRNAs, lncRNAs, piRNAs and many other RNA species. In silico analysis of datasets, representing only one RNA species is well established and a variety of tools and pipelines are available. However, attaining a more systematic view of how different players come together to regulate the expression of a gene or a group of genes requires a more intricate approach to data analysis. To fully understand complex transcriptional networks, datasets representing different RNA species need to be integrated. In this review, we will focus on miRNAs as key post-transcriptional regulators summarizing current computational approaches for miRNA:target gene prediction as well as new data-driven methods to tackle the problem of comprehensively and accurately dissecting miRNome-targetome interactions.

Proteomic and mitochondrial adaptations to early-life stress are distinct in juveniles and adults

Exposure to early-life stress (ELS) increases risk for poor mental and physical health outcomes that emerge at different stages across the lifespan. Yet, how age interacts with ELS to impact the expression of specific phenotypes remains largely unknown. An established limited-bedding paradigm was used to induce ELS in mouse pups over the early postnatal period. Initial analyses focused on the hippocampus, based on documented sensitivity to ELS in humans and various animal models, and the large body of data reporting anatomical and physiological outcomes in this structure using this ELS paradigm. An unbiased discovery proteomics approach revealed distinct adaptations in the non-nuclear hippocampal proteome in male versus female offspring at two distinct developmental stages: juvenile and adult. Gene ontology and KEGG pathway analyses revealed significant enrichment in proteins associated with mitochondria and the oxidative phosphorylation (OXPHOS) pathway in response to ELS in female hippocampus only. To determine whether the protein adaptations to ELS reflected altered function, mitochondrial respiration (driven through complexes II-IV) and complex I activity were measured in isolated hippocampal mitochondria using a Seahorse X96 Flux analyzer and immunocapture ELISA, respectively. ELS had no effect on basal respiration in either sex at either age. In contrast, ELS increased OXPHOS capacity in juvenile males and females, and reduced OXPHOS capacity in adult females but not adult males. A similar pattern of ELS-induced changes was observed for complex I activity. These data suggest that initial adaptations in juvenile hippocampus due to ELS were not sustained in adults. Mitochondrial adaptations to ELS were also exhibited peripherally by liver. Overall, the temporal distinctions in mitochondrial responses to ELS show that ELS-generated adaptations and outcomes are complex over the lifespan. This may contribute to differences in the timing of appearance of mental and physical disturbances, as well as potential sex differences that influence only select outcomes.

The SARS-CoV-2 host cell receptor ACE2 correlates positively with immunotherapy response and is a potential protective factor for cancer progression. Comput Struct Biotechnol J. 2020;18:2438-2444. [PMID: 32905022 PMCID: PMC7462778 DOI: 10.1016/j.csbj.2020.08.024]

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has infected more than 29 million people and has caused more than 900,000 deaths worldwide as of September 14, 2020. The SARS-CoV-2 human cell receptor ACE2 has recently received extensive attention for its role in SARS-CoV-2 infection. Many studies have also explored the association between ACE2 and cancer. However, a systemic investigation into associations between ACE2 and oncogenic pathways, tumor progression, and clinical outcomes in pan-cancer remains lacking. Using cancer genomics datasets from the Cancer Genome Atlas (TCGA) program, we performed computational analyses of associations between ACE2 expression and antitumor immunity, immunotherapy response, oncogenic pathways, tumor progression phenotypes, and clinical outcomes in 13 cancer cohorts. We found that ACE2 upregulation was associated with increased antitumor immune signatures and PD-L1 expression, and favorable anti-PD-1/PD-L1/CTLA-4 immunotherapy response. ACE2 expression levels inversely correlated with the activity of cell cycle, mismatch repair, TGF-β, Wnt, VEGF, and Notch signaling pathways. Moreover, ACE2 expression levels had significant inverse correlations with tumor proliferation, stemness, and epithelial-mesenchymal transition. ACE2 upregulation was associated with favorable survival in pan-cancer and in multiple individual cancer types. These results suggest that ACE2 is a potential protective factor for cancer progression. Our data may provide potential clinical implications for treating cancer patients infected with SARS-CoV-2.

Alterations in extracellular matrix composition during aging and photoaging of the skin

Human skin is composed of the cell-rich epidermis, the extracellular matrix (ECM) rich dermis, and the hypodermis. Within the dermis, a dense network of ECM proteins provides structural support to the skin and regulates a wide variety of signaling pathways which govern cell proliferation and other critical processes. Both intrinsic aging, which occurs steadily over time, and extrinsic aging (photoaging), which occurs as a result of external insults such as solar radiation, cause alterations to the dermal ECM. In this study, we utilized both quantitative and global proteomics, alongside single harmonic generation (SHG) and two-photon autofluorescence (TPAF) imaging, to assess changes in dermal composition during intrinsic and extrinsic aging. We find that both intrinsic and extrinsic aging result in significant decreases in ECM-supporting proteoglycans and structural ECM integrity, evidenced by decreasing collagen abundance and increasing fibril fragmentation. Intrinsic aging also produces changes distinct from those produced by photoaging, including reductions in elastic fiber and crosslinking enzyme abundance. In contrast, photoaging is primarily defined by increases in elastic fiber-associated protein and pro-inflammatory proteases. Changes associated with photoaging are evident even in young (mid 20s) sun-exposed forearm skin, indicating that proteomic evidence of photoaging is present decades prior to clinical signs of photoaging. GO term enrichment revealed that both intrinsic aging and photoaging share common features of chronic inflammation. The proteomic data has been deposited to the ProteomeXchange Consortium via the PRIDE partner repository with the data set identifier PXD015982.

•

Intrinsic aging and photoaging both decrease ECM-supporting proteoglycans and structural ECM.

•

Intrinsic aging produces reductions in elastic fiber and crosslinking enzyme abundance.

•

Photoaging results in increases in pro-inflammatory proteases and elastic fiber abundance.

•

Intrinsic aging and photoaging share common features associated with chronic inflammation.

•

Proteomic changes associated with photoaging are evident decades prior to clinical aging signs.

Predicting gene essentiality in Caenorhabditis elegans by feature engineering and machine-learning

Defining genes that are essential for life has major implications for understanding critical biological processes and mechanisms. Although essential genes have been identified and characterised experimentally using functional genomic tools, it is challenging to predict with confidence such genes from molecular and phenomic data sets using computational methods. Using extensive data sets available for the model organism Caenorhabditis elegans, we constructed here a machine-learning (ML)-based workflow for the prediction of essential genes on a genome-wide scale. We identified strong predictors for such genes and showed that trained ML models consistently achieve highly-accurate classifications. Complementary analyses revealed an association between essential genes and chromosomal location. Our findings reveal that essential genes in C. elegans tend to be located in or near the centre of autosomal chromosomes; are positively correlated with low single nucleotide polymorphim (SNP) densities and epigenetic markers in promoter regions; are involved in protein and nucleotide processing; are transcribed in most cells; are enriched in reproductive tissues or are targets for small RNAs bound to the argonaut CSR-1. Based on these results, we hypothesise an interplay between epigenetic markers and small RNA pathways in the germline, with transcription-based memory; this hypothesis warrants testing. From a technical perspective, further work is needed to evaluate whether the present ML-based approach will be applicable to other metazoans (including Drosophila melanogaster) for which comprehensive data sets (i.e. genomic, transcriptomic, proteomic, variomic, epigenetic and phenomic) are available.

Targeting autophagy-related protein kinases for potential therapeutic purpose

Autophagy, defined as a scavenging process of protein aggregates and damaged organelles mediated by lysosomes, plays a significant role in the quality control of macromolecules and organelles. Since protein kinases are integral to the autophagy process, it is critically important to understand the role of kinases in autophagic regulation. At present, intervention of autophagic processes by small-molecule modulators targeting specific kinases has becoming a reasonable and prevalent strategy for treating several varieties of human disease, especially cancer. In this review, we describe the role of some autophagy-related kinase targets and kinase-mediated phosphorylation mechanisms in autophagy regulation. We also summarize the small-molecule kinase inhibitors/activators of these targets, highlighting the opportunities of these new therapeutic agents.

PreMSIm: An R package for predicting microsatellite instability from the expression profiling of a gene panel in cancer

Microsatellite instability (MSI) is a genomic property of the cancers with defective DNA mismatch repair and is a useful marker for cancer diagnosis and treatment in diverse cancer types. In particular, MSI has been associated with the active immune checkpoint blockade therapy response in cancer. Most of computational methods for predicting MSI are based on DNA sequencing data and a few are based on mRNA expression data. Using the RNA-Seq pan-cancer datasets for three cancer cohorts (colon, gastric, and endometrial cancers) from The Cancer Genome Atlas (TCGA) program, we developed an algorithm (PreMSIm) for predicting MSI from the expression profiling of a 15-gene panel in cancer. We demonstrated that PreMSIm had high prediction performance in predicting MSI in most cases using both RNA-Seq and microarray gene expression datasets. Moreover, PreMSIm displayed superior or comparable performance versus other DNA or mRNA-based methods. We conclude that PreMSIm has the potential to provide an alternative approach for identifying MSI in cancer.

Data from differentially expressed proteins in platelet components associated with adverse transfusion reactions

The presented dataset was used for the study focused on the search for differentially expressed proteins in blood platelet components (PCs) associated with adverse transfusion reactions (ATRs). Pellets of ATR platelet components and their controls were subjected to high-throughput proteomics analysis using a Q Exactive high-resolution tandem mass spectrometer. The data reported here constitutes an extension of “Differential protein expression of blood platelet components associated with adverse transfusion reactions” article Aloui et al., 2018. The reported data herein have been deposited into the ProteomeXchange Consortium via the PRIDE partner repository with the dataset identifiers PXD003510 for the pooled platelet components (PPCs) and PXD008886 for the apheresis platelet components (SDA-PCs) associated with ATRs.

Data showing atherosclerosis-associated differentially methylated regions are often at enhancers

Atherosclerosis involves phenotypic modulation and transdifferentiation of vascular smooth muscle cells (SMCs). Data are given in tabular or figure format that illustrate genome-wide DNA methylation alterations in atherosclerotic vs. control aorta (athero DMRs). Data based upon publicly available chromatin state profiles are also shown for normal aorta, monocyte, and skeletal muscle tissue-specific DMRs and for aorta-specific chromatin features (enhancer chromatin, promoter chromatin, repressed chromatin, actively transcribed chromatin). Athero hypomethylated and hypermethylated DMRs as well as epigenetic and transcription profiles are described for the following genes: ACTA2, MYH10, MYH11 (SMC-associated genes); SMAD3 (a signaling gene for SMCs and other cell types); CD79B and SH3BP2 (leukocyte-associated genes); and TBX20 and genes in the HOXA, HOXB, HOXC, and HOXD clusters (T-box and homeobox developmental genes). The data reveal strong correlations between athero hypermethylated DMRs and regions of enhancer chromatin in aorta, which are discussed in the linked research article “Atherosclerosis-associated differentially methylated regions can reflect the disease phenotype and are often at enhancers” (M. Lacey et al., 2019).

Hypocrellin A-based photodynamic action induces apoptosis in A549 cells through ROS-mediated mitochondrial signaling pathway

Over recent decades, many studies have reported that hypocrellin A (HA) can eliminate cancer cells with proper irradiation in several cancer cell lines. However, the precise molecular mechanism underlying its anticancer effect has not been fully defined. HA-mediated cytotoxicity and apoptosis in human lung adenocarcinoma A549 cells were evaluated after photodynamic therapy (PDT). A temporal quantitative proteomics approach by isobaric tag for relative and absolute quantitation (iTRAQ) 2D liquid chromatography with tandem mass spectrometric (LC–MS/MS) was introduced to help clarify molecular cytotoxic mechanisms and identify candidate targets of HA-induced apoptotic cell death. Specific caspase inhibitors were used to further elucidate the molecular pathway underlying apoptosis in PDT-treated A549 cells. Finally, down-stream apoptosis-related protein was evaluated. Apoptosis induced by HA was associated with cell shrinkage, externalization of cell membrane phosphatidylserine, DNA fragmentation, and mitochondrial disruption, which were preceded by increased intracellular reactive oxygen species (ROS) generations. Further studies showed that PDT treatment with 0.08 µmol/L HA resulted in mitochondrial disruption, pronounced release of cytochrome c, and activation of caspase-3, -9, and -7. Together, HA may be a possible therapeutic agent directed toward mitochondria and a promising photodynamic anticancer candidate for further evaluation.

Variations in Copepod Proteome and Respiration Rate in Association with Diel Vertical Migration and Circadian Cycle

The diel vertical migration of zooplankton is a process during which individuals spend the night in surface waters and retreat to depth during the daytime, with substantial implications for carbon transport and the ecology of midwater ecosystems. The physiological consequences of this daily pattern have, however, been poorly studied beyond investigations of speed and the energetic cost of swimming. Many other processes are likely influenced, such as fuel use, energetic trade-offs, underlying diel (circadian) rhythms, and antioxidant responses. Using a new reference transcriptome, proteomic analyses were applied to compare the physiological state of a migratory copepod, Pleuromamma xiphias, immediately after arriving to the surface at night and six hours later. Oxygen consumption was monitored semi-continuously to explore underlying cyclical patterns in metabolic rate under dark-dark conditions. The proteomic analysis suggests a distinct shift in physiology that reflects migratory exertion and changes in metabolism. These proteomic analyses are supported by the respiration experiments, which show an underlying cycle in metabolic rate, with a peak at dawn. This project generates molecular tools (transcriptome and proteome) that will allow for more detailed understanding of the underlying physiological processes that influence and are influenced by diel vertical migration. Further, these studies suggest that P. xiphias is a tractable model for continuing investigations of circadian and diel vertical migration influences on plankton physiology. Previous studies did not account for this cyclic pattern of respiration and may therefore have unrepresented respiratory carbon fluxes from copepods by about 24%.

siRNA Library Screening Identifies a Druggable Immune-Signature Driving Esophageal Adenocarcinoma Cell Growth

BACKGROUND & AIMS

Effective therapeutic approaches are urgently required to tackle the alarmingly poor survival outcomes in esophageal adenocarcinoma (EAC) patients. EAC originates from within the intestinal-type metaplasia, Barrett's esophagus, a condition arising on a background of gastroesophageal reflux disease and associated inflammation.

METHODS

This study used a druggable genome small interfering RNA (siRNA) screening library of 6022 siRNAs in conjunction with bioinformatics platforms, genomic studies of EAC tissues, somatic variation data of EAC from The Cancer Genome Atlas data of EAC, and pathologic and functional studies to define novel EAC-associated, and targetable, immune factors.

RESULTS

By using a druggable genome library we defined genes that sustain EAC cell growth, which included an unexpected immunologic signature. Integrating Cancer Genome Atlas data with druggable siRNA targets showed a striking concordance and an EAC-specific gene amplification event associated with 7 druggable targets co-encoded at Chr6p21.1. Over-representation of immune pathway-associated genes supporting EAC cell growth included leukemia inhibitory factor, complement component 1, q subcomponent A chain (C1QA), and triggering receptor expressed on myeloid cells 2 (TREM2), which were validated further as targets sharing downstream signaling pathways through genomic and pathologic studies. Finally, targeting the triggering receptor expressed on myeloid cells 2-, C1q-, and leukemia inhibitory factor-activated signaling pathways (TYROBP-spleen tyrosine kinase and JAK-STAT3) with spleen tyrosine kinase and Janus-activated kinase inhibitor fostamatinib R788 triggered EAC cell death, growth arrest, and reduced tumor burden in NOD scid gamma mice.

CONCLUSIONS

These data highlight a subset of genes co-identified through siRNA targeting and genomic studies of expression and somatic variation, specifically highlighting the contribution that immune-related factors play in support of EAC development and suggesting their suitability as targets in the treatment of EAC.

Strategy to Identify and Test Putative Light-Sensitive Non-Opsin G-Protein-Coupled Receptors: A Case Study

The rise of high-throughput RNA sequencing (RNA-seq) and de novo transcriptome assembly has had a transformative impact on how we identify and study genes in the phototransduction cascade of non-model organisms. But the advantage provided by the nearly automated annotation of RNA-seq transcriptomes may at the same time hinder the possibility for gene discovery and the discovery of new gene functions. For example, standard functional annotation based on domain homology to known protein families can only confirm group membership, not identify the emergence of new biochemical function. In this study, we show the importance of developing a strategy that circumvents the limitations of semiautomated annotation and apply this workflow to photosensitivity as a means to discover non-opsin photoreceptors. We hypothesize that non-opsin G-protein-coupled receptor (GPCR) proteins may have chromophore-binding lysines in locations that differ from opsin. Here, we provide the first case study describing non-opsin light-sensitive GPCRs based on tissue-specific RNA-seq data of the common bay scallop Argopecten irradians (Lamarck, 1819). Using a combination of sequence analysis and three-dimensional protein modeling, we identified two candidate proteins. We tested their photochemical properties and provide evidence showing that these two proteins incorporate 11-cis and/or all-trans retinal and react to light photochemically. Based on this case study, we demonstrate that there is potential for the discovery of new light-sensitive GPCRs, and we have developed a workflow that starts from RNA-seq assemblies to the discovery of new non-opsin, GPCR-based photopigments.

Identification of Genes for Synthesis of the Blue Pigment, Biliverdin IXα, in the Blue Coral Heliopora coerulea

Heliopora coerulea is the only species in the subclass Octocorallia that has a crystalline aragonite skeleton. The skeleton has been reported to contain the blue pigment, biliverdin IXα, which is formed by heme oxygenase (HO) during heme decomposition. There is little information regarding gene expression in H. coerulea; therefore, the biosynthesis pathway for biliverdin IXα is poorly understood. To identify the genes related to heme synthesis and degradation, metatranscripts of H. coerulea and its symbiont Symbiodinium spp. were sequenced and separated from the host- and symbiont-derived sequences. From the metatranscriptome analyses, all genes for heme synthesis and three HOs were isolated from the host and symbiont. From our phylogenetic and amino acid analysis, we noted that one of the HO isoforms in the host coral was predicted to possess HO activity. However, biliverdin reductase, which reduces biliverdin to bilirubin, was not identified in the present study. Similarly, biliverdin reductase was not identified in the transcripts of the red coral Corallium rubrum, a species that also belongs to Octocorallia. However, genes related to heme synthesis and HO were found in C. rubrum. We speculate that Heliopora coerulea can produce biliverdin and accumulate it in the skeleton, while red corals and other Octocorallia species cannot. Further information from molecular studies of H. coerulea will provide insights into the synthesis of biliverdin IXα, the blue pigment in the hard crystalline aragonite skeleton, and will be fundamental to future ecological and physiological studies.

Transcriptional Regulation by ATOH1 and its Target SPDEF in the Intestine

BACKGROUND & AIMS

The transcription factor atonal homolog 1 (ATOH1) controls the fate of intestinal progenitors downstream of the Notch signaling pathway. Intestinal progenitors that escape Notch activation express high levels of ATOH1 and commit to a secretory lineage fate, implicating ATOH1 as a gatekeeper for differentiation of intestinal epithelial cells. Although some transcription factors downstream of ATOH1, such as SPDEF, have been identified to specify differentiation and maturation of specific cell types, the bona fide transcriptional targets of ATOH1 still largely are unknown. Here, we aimed to identify ATOH1 targets and to identify transcription factors that are likely to co-regulate gene expression with ATOH1.

METHODS

We used a combination of chromatin immunoprecipitation and messenger RNA-based high-throughput sequencing (ChIP-seq and RNA-seq), together with cell sorting and transgenic mice, to identify direct targets of ATOH1, and establish the epistatic relationship between ATOH1 and SPDEF.

RESULTS

By using unbiased genome-wide approaches, we identified more than 700 genes as ATOH1 transcriptional targets in adult small intestine and colon. Ontology analysis indicated that ATOH1 directly regulates genes involved in specification and function of secretory cells. De novo motif analysis of ATOH1 targets identified SPDEF as a putative transcriptional co-regulator of ATOH1. Functional epistasis experiments in transgenic mice show that SPDEF amplifies ATOH1-dependent transcription but cannot independently initiate transcription of ATOH1 target genes.

CONCLUSIONS

This study unveils the direct targets of ATOH1 in the adult intestines and illuminates the transcriptional events that initiate the specification and function of intestinal secretory lineages.

Comparative analysis of transcriptomes in aerial stems and roots of Ephedra sinica based on high-throughput mRNA sequencing

Ephedra plants are taxonomically classified as gymnosperms, and are medicinally important as the botanical origin of crude drugs and as bioresources that contain pharmacologically active chemicals. Here we show a comparative analysis of the transcriptomes of aerial stems and roots of Ephedra sinica based on high-throughput mRNA sequencing by RNA-Seq. De novo assembly of short cDNA sequence reads generated 23,358, 13,373, and 28,579 contigs longer than 200 bases from aerial stems, roots, or both aerial stems and roots, respectively. The presumed functions encoded by these contig sequences were annotated by BLAST (blastx). Subsequently, these contigs were classified based on gene ontology slims, Enzyme Commission numbers, and the InterPro database. Furthermore, comparative gene expression analysis was performed between aerial stems and roots. These transcriptome analyses revealed differences and similarities between the transcriptomes of aerial stems and roots in E. sinica. Deep transcriptome sequencing of Ephedra should open the door to molecular biological studies based on the entire transcriptome, tissue- or organ-specific transcriptomes, or targeted genes of interest.

The satiating hormone amylin enhances neurogenesis in the area postrema of adult rats

OBJECTIVE

Adult neurogenesis in the subgranular zone and subventricular zone is generally accepted, but its existence in other brain areas is still controversial. Circumventricular organs, such as the area postrema (AP) have recently been described as potential neurogenic niches in the adult brain. The AP is the major site of action of the satiating hormone amylin. Amylin has been shown to promote the formation of neuronal projections originating from the AP in neonatal rodents but the role of amylin in adult neurogenesis remains unknown.

METHODS

To test this, we first performed an RNA-sequencing of the AP of adult rats acutely injected with either amylin (20 μg/kg), amylin plus the amylin receptor antagonist AC187 (500 μg/kg) or vehicle. Second, animals were subcutaneously equipped with minipumps releasing either amylin (50 μg/kg/day) or vehicle for 3 weeks to assess cell proliferation and differentiation with the 5'-bromo-2-deoxyuridine (BrdU) technique.

RESULTS

Acute amylin injections affected genes involved in pathways and processes that control adult neurogenesis. Amylin consistently upregulated NeuroD1 transcript and protein in the adult AP, and this effect was blocked by the co-administration of AC187. Further, chronic amylin treatment increased the number of newly proliferated AP-cells and significantly promoted their differentiation into neurons rather than astrocytes.

CONCLUSION

Our findings revealed a novel role of the satiating hormone amylin in promoting neurogenesis in the AP of adult rats.

Transcriptomic Analysis Identifies the Effect of Beta-Blocking Agents on a Molecular Pathway of Contraction in the Heart and Predicts Response to Therapy

Over the last decades, beta-blockers have been a key component of heart failure therapy. However, currently there is no method to identify patients who will benefit from beta-blocking therapy versus those who will be unresponsive or worsen. Furthermore, there is an unmet need to better understand molecular mechanisms through which heart failure therapies, such as beta-blockers, improve cardiac function, in order to design novel targeted therapies. Solving these issues is an important step towards personalized medicine. Here, we present a comprehensive transcriptomic analysis of molecular pathways that are affected by beta-blocking agents and a transcriptomic biomarker to predict therapy response.

A plastic SQSTM1/p62-dependent autophagic reserve maintains proteostasis and determines proteasome inhibitor susceptibility in multiple myeloma cells

Multiple myeloma (MM) is the paradigmatic proteasome inhibitor (PI) responsive cancer, but many patients fail to respond. An attractive target to enhance sensitivity is (macro)autophagy, recently found essential to bone marrow plasma cells, the normal counterpart of MM. Here, integrating proteomics with hypothesis-driven strategies, we identified the autophagic cargo receptor and adapter protein, SQSTM1/p62 as an essential component of an autophagic reserve that not only synergizes with the proteasome to maintain proteostasis, but also mediates a plastic adaptive response to PIs, and faithfully reports on inherent PI sensitivity. Lentiviral engineering revealed that SQSTM1 is essential for MM cell survival and affords specific PI protection. Under basal conditions, SQSTM1-dependent autophagy alleviates the degradative burden on the proteasome by constitutively disposing of substantial amounts of ubiquitinated proteins. Indeed, its inhibition or stimulation greatly sensitized to, or protected from, PI-induced protein aggregation and cell death. Moreover, under proteasome stress, myeloma cells selectively enhanced SQSTM1 de novo expression and reset its vast endogenous interactome, diverting SQSTM1 from signaling partners to maximize its association with ubiquitinated proteins. Saturation of such autophagic reserve, as indicated by intracellular accumulation of undigested SQSTM1-positive aggregates, specifically discriminated patient-derived myelomas inherently susceptible to PIs from primarily resistant ones. These aggregates correlated with accumulation of the endoplasmic reticulum, which comparative proteomics identified as the main cell compartment targeted by autophagy in MM. Altogether, the data integrate autophagy into our previously established proteasome load-versus-capacity model, and reveal SQSTM1 aggregation as a faithful marker of defective proteostasis, defining a novel prognostic and therapeutic framework for MM.

TXNDC17 promotes paclitaxel resistance via inducing autophagy in ovarian cancer

Paclitaxel is recommended as a first-line chemotherapeutic agent against ovarian cancer, but drug resistance becomes a major limitation of its success clinically. The key molecule or mechanism associated with paclitaxel resistance in ovarian cancer still remains unclear. Here, we showed that TXNDC17 screened from 356 differentially expressed proteins by LC-MS/MS label-free quantitative proteomics was more highly expressed in paclitaxel-resistant ovarian cancer cells and tissues, and the high expression of TXNDC17 was associated with poorer prognostic factors and exhibited shortened survival in 157 ovarian cancer patients. Moreover, paclitaxel exposure induced upregulation of TXNDC17 and BECN1 expression, increase of autophagosome formation, and autophagic flux that conferred cytoprotection for ovarian cancer cells from paclitaxel. TXNDC17 inhibition by siRNA or enforced overexpression by a pcDNA3.1(+)-TXNDC17 plasmid correspondingly decreased or increased the autophagy response and paclitaxel resistance. Additionally, the downregulation of BECN1 by siRNA attenuated the activation of autophagy and cytoprotection from paclitaxel induced by TXNDC17 overexpression in ovarian cancer cells. Thus, our findings suggest that TXNDC17, through participation of BECN1, induces autophagy and consequently results in paclitaxel resistance in ovarian cancer. TXNDC17 may be a potential predictor or target in ovarian cancer therapeutics.

Type I interferon related genes are common genes on the early stage after vaccination by meta-analysis of microarray data

The objective of this study was to find common immune mechanism across different kinds of vaccines. A meta-analysis of microarray datasets was performed using publicly available microarray Gene Expression Omnibus (GEO) and Array Express data sets of vaccination records. Seven studies (out of 35) were selected for this meta-analysis. A total of 447 chips (145 pre-vaccination and 302 post-vaccination) were included. Significance analysis of microarrays (SAM) program was used for screening differentially expressed genes (DEGs). Functional pathway enrichment for the DEGs was conducted in DAVID Gene Ontology (GO) database. Twenty DEGs were identified, of which 10 up-regulated genes involved immune response. Six of which were type I interferon (IFN) related genes, including LY6E, MX1, OAS3, IFI44L, IFI6 and IFITM3. Ten down-regulated genes mainly mediated negative regulation of cell proliferation and cell motion. Results of a subgroup analysis showed that although the kinds of genes varied widely between days 3 and 7 post vaccination, the pathways between them are basically the same, such as immune response and response to viruses, etc. For an independent verification of these 6 type I IFN related genes, peripheral blood mononuclear cells (PBMCs) were collected at baseline and day 3 after the vaccination from 8 Enterovirus 71(EV71) vaccinees and were assayed by RT-PCR. Results showed that the 6 DEGs were also upregulated in EV71 vaccinees. In summary, meta-analysis methods were used to explore the immune mechanism of vaccines and results indicated that the type I IFN related genes and corresponding pathways were common in early immune responses for different kinds of vaccines.

Prediction of the anti-inflammatory mechanisms of curcumin by module-based protein interaction network analysis

Curcumin, the medically active component from Curcuma longa (Turmeric), is widely used to treat inflammatory diseases. Protein interaction network (PIN) analysis was used to predict its mechanisms of molecular action. Targets of curcumin were obtained based on ChEMBL and STITCH databases. Protein–protein interactions (PPIs) were extracted from the String database. The PIN of curcumin was constructed by Cytoscape and the function modules identified by gene ontology (GO) enrichment analysis based on molecular complex detection (MCODE). A PIN of curcumin with 482 nodes and 1688 interactions was constructed, which has scale-free, small world and modular properties. Based on analysis of these function modules, the mechanism of curcumin is proposed. Two modules were found to be intimately associated with inflammation. With function modules analysis, the anti-inflammatory effects of curcumin were related to SMAD, ERG and mediation by the TLR family. TLR9 may be a potential target of curcumin to treat inflammation.

Dynamics of re-constitution of the human nuclear proteome after cell division is regulated by NLS-adjacent phosphorylation

Phosphorylation by the cyclin-dependent kinase 1 (Cdk1) adjacent to nuclear localization signals (NLSs) is an important mechanism of regulation of nucleocytoplasmic transport. However, no systematic survey has yet been performed in human cells to analyze this regulatory process, and the corresponding cell-cycle dynamics have not yet been investigated. Here, we focused on the human proteome and found that numerous proteins, previously not identified in this context, are associated with Cdk1-dependent phosphorylation sites adjacent to their NLSs. Interestingly, these proteins are involved in key regulatory events of DNA repair, epigenetics, or RNA editing and splicing. This finding indicates that cell-cycle dependent events of genome editing and gene expression profiling may be controlled by nucleocytoplasmic trafficking. For in-depth investigations, we selected a number of these proteins and analyzed how point mutations, expected to modify the phosphorylation ability of the NLS segments, perturb nucleocytoplasmic localization. In each case, we found that mutations mimicking hyper-phosphorylation abolish nuclear import processes. To understand the mechanism underlying these phenomena, we performed a video microscopy-based kinetic analysis to obtain information on cell-cycle dynamics on a model protein, dUTPase. We show that the NLS-adjacent phosphorylation by Cdk1 of human dUTPase, an enzyme essential for genomic integrity, results in dynamic cell cycle-dependent distribution of the protein. Non-phosphorylatable mutants have drastically altered protein re-import characteristics into the nucleus during the G1 phase. Our results suggest a dynamic Cdk1-driven mechanism of regulation of the nuclear proteome composition during the cell cycle.

Identification of an epigenetic signature of early mouse liver regeneration that is disrupted by Zn-HDAC inhibition

Liver regeneration has been well studied with hope of discovering strategies to improve liver disease outcomes. Nevertheless, the signals that initiate such regeneration remain incompletely defined, and translation of mechanism-based pro-regenerative interventions into new treatments for hepatic diseases has not yet been achieved. We previously reported the isoform-specific regulation and essential function of zinc-dependent histone deacetylases (Zn-HDACs) during mouse liver regeneration. Those data suggest that epigenetically regulated anti-proliferative genes are deacetylated and transcriptionally suppressed by Zn-HDAC activity or that pro-regenerative factors are acetylated and induced by such activity in response to partial hepatectomy (PH). To investigate these possibilities, we conducted genome-wide interrogation of the liver histone acetylome during early PH-induced liver regeneration in mice using acetyL-histone chromatin immunoprecipitation and next generation DNA sequencing. We also compared the findings of that study to those seen during the impaired regenerative response that occurs with Zn-HDAC inhibition. The results reveal an epigenetic signature of early liver regeneration that includes both hyperacetylation of pro-regenerative factors and deacetylation of anti-proliferative and pro-apoptotic genes. Our data also show that administration of an anti-regenerative regimen of the Zn-HDAC inhibitor suberoylanilide hydroxamic acid (SAHA) not only disrupts gene-specific pro-regenerative changes in liver histone deacetylation but also reverses PH-induced effects on histone hyperacetylation. Taken together, these studies offer new insight into and suggest novel hypotheses about the epigenetic mechanisms that regulate liver regeneration.

UHRF1 promotes proliferation of gastric cancer via mediating tumor suppressor gene hypermethylation

Epigenetic changes play significant roles in cancer development. UHRF1, an epigenetic regulator, has been shown to be overexpressed and to coordinate tumor suppressor gene (TSG) silencing in several cancers. In a previous study, we found that UHRF1 promoted gastric cancer (GC) invasion and metastasis. However, the role and underlying mechanism of UHRF1 in GC carcinogenesis remain largely unknown. In the present study, we investigated UHRF1 expression and function in GC proliferation and explored its downstream regulatory mechanism. The results demonstrated that UHRF1 overexpression was an independent and significant predictor of GC prognosis. Downregulation of UHRF1 suppressed GC proliferation and growth in vitro and in vivo, and UHRF1 upregulation showed opposite effects. Furthermore, downregulation of UHRF1 reactivated 7 TSGs, including CDX2, CDKN2A, RUNX3, FOXO4, PPARG, BRCA1 and PML, via promoter demethylation. These results provide insight into the GC proliferation process, and suggest that targeting UHRF1 represents a new therapeutic approach to block GC development.

Comparative gene expression profiles in pancreatic islets associated with agouti yellow mutation and PACAP overexpression in mice

In diabetes mellitus, pituitary adenylate cyclase-activating polypeptide (PACAP) has insulinotropic and glucose-lowering properties. We previously demonstrated that transgenic mice overexpressing PACAP in pancreatic β-cells (PACAP-Tg) show attenuated pancreatic islet hyperplasia and hyperinsulinemia in type 2 diabetic models. To explore the underlying mechanisms, here we crossed PACAP-Tg mice with lethal yellow agouti (KKAy) diabetic mice, and performed gene chip analysis of laser capture microdissected pancreatic islets from four F₁ offspring genotypes (wild-type, PACAP-Tg, KKAy, and PACAP-Tg:KKAy). We identified 1371 probes with >16-fold differences between at least one pair of genotypes, and classified the probes into five clusters with characteristic expression patterns. Gene ontology enrichment analysis showed that genes involved in the terms ribosome and intracellular organelles such as ribonucleoprotein complex, mitochondrion, and chromosome organization were significantly enriched in clusters characterized by up-regulated genes in PACAP-Tg:KKAy mice compared with KKAy mice. These results may provide insight into the mechanisms of diabetes that accompany islet hyperplasia and amelioration by PACAP.

•

PACAP overexpressed in KKAy diabetic mice is known to exert antidiabetic effects.

•

We performed gene chip analysis of pancreatic islets in these mice.

•

Gene ontology analysis was performed for genes classified into five clusters.

•

Genes involved in the terms ribosome, mitochondrion, and chromosome were enriched.

•

These pathways may be involved in the mechanism by which PACAP ameliorates diabetes.

PDL1 expression is an independent prognostic factor in localized GIST

Gastrointestinal stromal tumors (GIST) are the most frequently occurring digestive sarcomas. The prognosis of localized GIST is heterogeneous, notably for patients with an Armed Forces Institute of Pathology (AFIP) intermediate or high risk of relapse. Despite imatinib effectiveness, it is crucial to develop therapies able to overcome the resistance mechanisms. The immune system represents an attractive prognostic and therapeutic target. The Programmed cell Death 1 (PD1)/programmed cell death ligand 1 (PDL1) pathway is a key inhibitor of the immune response; recently, anti-PD1 and anti-PDL1 drugs showed very promising results in patients with solid tumors. However, PDL1 expression has never been studied in GIST. Our objective was to analyze PDL1 expression in a large series of clinical samples. We analyzed mRNA expression data of 139 operated imatinib-untreated localized GIST profiled using DNA microarrays and searched for correlations with histoclinical features including postoperative metastatic relapse. PDL1 expression was heterogeneous across tumors and was higher in AFIP low-risk than in high-risk samples, and in samples without than with metastatic relapse. PDL1 expression was associated with immunity-related parameters such as T-cell-specific and CD8⁺ T-cell-specific gene expression signatures and probabilities of activation of interferon α (IFNα), IFNγ, and tumor necrosis factor α (TNFα) pathways, suggesting positive correlation with a cytotoxic T-cell response. In multivariate analysis, the PDL1-low group was associated with a higher metastatic risk independently of the AFIP classification and the KIT mutational status. In conclusion, PDL1 expression refines the prediction of metastatic relapse in localized GIST and might improve our ability to better tailor adjuvant imatinib. In the metastatic setting, PDL1 expression might guide the use of PDL1 inhibitors, alone or associated with tyrosine kinase inhibitors.

Novel DNA methylation profiles associated with key gene regulation and transcription pathways in blood and placenta of growth-restricted neonates

Fetal growth is determined by the feto-placental genome interacting with the maternal in utero environment. Failure of this interplay leads to poor placental development and fetal growth restriction (FGR), which is associated with future metabolic disease. We investigated whether whole genome methylation differences existed in umbilical cord blood and placenta, between gestational-matched, FGR, and appropriately grown (AGA) neonates. Using the Infinium HumanMethylation450 BeadChip®, we found that DNA from umbilical cord blood of FGR born at term (n = 19) had 839 differentially methylated positions (DMPs) that reached genome-wide significance compared with AGA (n = 18). Using gestational age as a continuous variable, we identified 76,249 DMPs in cord blood (adj. P < 0.05) of which 121 DMPs were common to the 839 DMPs and were still evident when comparing 12 FGR with 12 AGA [39.9 ± 1.2 vs. 40.0 ± 1.0 weeks (mean ± SD), respectively]. A total of 53 DMPs had a β methylation difference >10% and 25 genes were co-methylated more than twice within 1000 base pairs. Gene Ontology (GO) analysis of DMPs supported their involvement in gene regulation and transcription pathways related to organ development and metabolic function. A similar profile of DMPs was found across different cell types in the cord blood. At term, no DMPs between FGR and AGA placentae reached genome-wide significance, validated with an external dataset. GO analysis of 284 pre-term, placental DMPs associated with autophagy, oxidative stress and hormonal responses. Growth restricted neonates have distinct DNA methylation profiles in pre-term placenta and in cord blood at birth, which may predispose to future adult disease.

Improved abiotic stress tolerance of bermudagrass by exogenous small molecules

As a widely used warm-season turfgrass in landscapes and golf courses, bermudagrass encounters multiple abiotic stresses during the growth and development. Physiology analysis indicated that abiotic stresses induced the accumulation of ROS and decline of photosynthesis, resulting in increased cell damage and inhibited growth. Proteomic and metabolomic approaches showed that antioxidant enzymes and osmoprotectant contents (sugar, sucrose, dehydrin, proline) were extensively changed under abiotic stress conditions. Exogenous application of small molecules, such as ABA, NO, CaCl2, H2S, polyamine and melatonin, could effectively alleviate damages caused by multiple abiotic stresses, including drought, salt, heat and cold. Based on high through-put RNA seq analysis, genes involved in ROS, transcription factors, hormones, and carbohydrate metabolisms were largely enriched. The data indicated that small molecules induced the accumulation of osmoprotectants and antioxidants, kept cell membrane integrity, increased photosynthesis and kept ion homeostasis, which protected bermudagrass from damages caused by abiotic stresses.

Insights in the etiopathology of galactosyltransferase II (GalT-II) deficiency from transcriptome-wide expression profiling of skin fibroblasts of two sisters with compound heterozygosity for two novel B3GALT6 mutations

Mutations in B3GALT6, encoding the galactosyltransferase II (GalT-II) involved in the synthesis of the glycosaminoglycan (GAG) linkage region of proteoglycans (PGs), have recently been associated with a spectrum of connective tissue disorders, including spondyloepimetaphyseal dysplasia with joint laxity type 1 (SEMDJL1) and Ehlers–Danlos-like syndrome. Here, we report on two sisters compound heterozygous for two novel B3GALT6 mutations that presented with severe short stature and progressive kyphoscoliosis, joint hypermobility and laxity, hyperextensible skin, platyspondyly, short ilia, and elbow malalignment. Microarray-based transcriptome analysis revealed the differential expression of several genes encoding extracellular matrix (ECM) structural components, including COMP, SPP1, COL5A1, and COL15A1, enzymes involved in GAG synthesis and in ECM remodeling, such as CSGALNACT1, CHPF, LOXL3, and STEAP4, signaling transduction molecules of the TGFβ/BMP pathway, i.e., GDF6, GDF15, and BMPER, and transcription factors of the HOX and LIM families implicated in skeletal and limb development. Immunofluorescence analyses confirmed the down-regulated expression of some of these genes, in particular of the cartilage oligomeric matrix protein and osteopontin, encoded by COMP and SPP1, respectively, and showed the predominant reduction and disassembly of the heparan sulfate specific GAGs, as well as of the PG perlecan and type III and V collagens. The key role of GalT-II in GAG synthesis and the crucial biological functions of PGs are consistent with the perturbation of many physiological functions that are critical for the correct architecture and homeostasis of various connective tissues, including skin, bone, cartilage, tendons, and ligaments, and generates the wide phenotypic spectrum of GalT-II-deficient patients.

•

Clinical features/molecular characterization of two patients with spondyloepimetaphyseal dysplasia with joint laxity type 1

•

Identification of two novel B3GALT6 mutations

•

First report of transcriptome-wide gene expression profiling on GalT-II-deficient fibroblasts

•

Immunofluorescence studies of several ECM structural components in GalT-II-deficient cells

•

Enlargement of the knowledge on the GalT-II deficiency’s molecular pathogenesis

Influence of the Prader-Willi syndrome imprinting center on the DNA methylation landscape in the mouse brain

Reduced representation bisulfite sequencing (RRBS) was used to analyze DNA methylation patterns across the mouse brain genome in mice carrying a deletion of the Prader-Willi syndrome imprinting center (PWS-IC) on either the maternally- or paternally-inherited chromosome. Within the ~3.7 Mb imprinted Angelman/Prader-Willi syndrome (AS/PWS) domain, 254 CpG sites were interrogated for changes in methylation due to PWS-IC deletion. Paternally-inherited deletion of the PWS-IC increased methylation levels ~2-fold at each CpG site (compared to wild-type controls) at differentially methylated regions (DMRs) associated with 5' CpG island promoters of paternally-expressed genes; these methylation changes extended, to a variable degree, into the adjacent CpG island shores. Maternal PWS-IC deletion yielded little or no changes in methylation at these DMRs, and methylation of CpG sites outside of promoter DMRs also was unchanged upon maternal or paternal PWS-IC deletion. Using stringent ascertainment criteria, ~750,000 additional CpG sites were also interrogated across the entire mouse genome. This analysis identified 26 loci outside of the imprinted AS/PWS domain showing altered DNA methylation levels of ≥25% upon PWS-IC deletion. Curiously, altered methylation at 9 of these loci was a consequence of maternal PWS-IC deletion (maternal PWS-IC deletion by itself is not known to be associated with a phenotype in either humans or mice), and 10 of these loci exhibited the same changes in methylation irrespective of the parental origin of the PWS-IC deletion. These results suggest that the PWS-IC may affect DNA methylation at these loci by directly interacting with them, or may affect methylation at these loci through indirect downstream effects due to PWS-IC deletion. They further suggest the PWS-IC may have a previously uncharacterized function outside of the imprinted AS/PWS domain.

EDI3 links choline metabolism to integrin expression, cell adhesion and spreading

Endometrial carcinoma differential 3 (EDI3) was the first member of the glycerophosphodiesterase (GDE) protein family shown to be associated with cancer. Our initial work demonstrated that endometrial and ovarian cancer patients with primary tumors overexpressing EDI3 had a higher risk of developing metastasis and decreased survival. Further analysis indicated that EDI3 cleaves glycerophosphocholine to choline and glycerol-3-phosphate, increases the levels of active PKC, and enhances the migratory activity of tumor cells. Despite these initial findings, EDI3 remained mainly uncharacterized. Therefore, to obtain an overview of processes in which EDI3 may be involved, gene array analysis was performed using MCF-7 breast cancer cells after EDI3 knockdown compared with a non-targeting control siRNA. Several biological motifs were altered, including an enrichment of genes involved in integrin-mediated signaling. More specifically, silencing of EDI3 in MCF-7 and OVCAR-3 cells was associated with reduced expression of the key receptor subunit integrin β1, leading to decreased cell attachment and spreading accompanied by delayed formation of cell protrusions. To confirm these results, we stably overexpressed EDI3 in MCF-7 cells which led to elevated integrin β1 expression associated with enhanced cell attachment and spreading - two processes critical for metastasis. In conclusion, our data provide further insight into the role of EDI3 during cancer progression.

Integrated mRNA and microRNA transcriptome sequencing characterizes sequence variants and mRNA-microRNA regulatory network in nasopharyngeal carcinoma model systems

Nasopharyngeal carcinoma (NPC) is a prevalent malignancy in Southeast Asia among the Chinese population. Aberrant regulation of transcripts has been implicated in many types of cancers including NPC. Herein, we characterized mRNA and miRNA transcriptomes by RNA sequencing (RNASeq) of NPC model systems. Matched total mRNA and small RNA of undifferentiated Epstein–Barr virus (EBV)-positive NPC xenograft X666 and its derived cell line C666, well-differentiated NPC cell line HK1, and the immortalized nasopharyngeal epithelial cell line NP460 were sequenced by Solexa technology. We found 2812 genes and 149 miRNAs (human and EBV) to be differentially expressed in NP460, HK1, C666 and X666 with RNASeq; 533 miRNA–mRNA target pairs were inversely regulated in the three NPC cell lines compared to NP460. Integrated mRNA/miRNA expression profiling and pathway analysis show extracellular matrix organization, Beta-1 integrin cell surface interactions, and the PI3K/AKT, EGFR, ErbB, and Wnt pathways were potentially deregulated in NPC. Real-time quantitative PCR was performed on selected mRNA/miRNAs in order to validate their expression. Transcript sequence variants such as short insertions and deletions (INDEL), single nucleotide variant (SNV), and isomiRs were characterized in the NPC model systems. A novel TP53 transcript variant was identified in NP460, HK1, and C666. Detection of three previously reported novel EBV-encoded BART miRNAs and their isomiRs were also observed. Meta-analysis of a model system to a clinical system aids the choice of different cell lines in NPC studies. This comprehensive characterization of mRNA and miRNA transcriptomes in NPC cell lines and the xenograft provides insights on miRNA regulation of mRNA and valuable resources on transcript variation and regulation in NPC, which are potentially useful for mechanistic and preclinical studies.

•

Using RNASeq we characterized the mRNA and miRNA transcriptomes in NPC and NP models.

•

2812 Genes and 149 miRNAs (human and EBV) were differentially expressed in NPC vs NP models.

•

533 miRNA–mRNA target pairs were inversely regulated in HK1, C666, and X666 vs NP460.

•

ECM, β1 integrin, PI3K/AKT, EGFR, ErbB, and Wnt pathways appeared to be deregulated in NPC.

•

A novel TP53 mutation was identified in NP460, HK1, and C666.