Interactive Clustered Heat Map Builder: An easy web-based tool for creating sophisticated clustered heat maps

Improved humanized mouse model of Staphylococcus aureus infection

Staphylococcus aureus is a significant cause of pulmonary infections, but existing mouse models fail to recapitulate human-specific responses. In this study, we developed a novel mouse model of S. aureus infection using humanized mice implanted with autologous fetal lung tissue. We show that these human lung implants support S. aureus survival and dissemination. Immunological profiling revealed extensive immune cell death after infection and an absence of chemokine induction. Transcriptomic profiling of the human lung implants revealed significant changes in gene expression, including NF-κB and JAK/STAT signaling. We identified upregulation of Cyp24a1, suggesting a role for vitamin D metabolism in host defense, but it had a mild effect on dissemination. Examination of the bacterial response to the host environment, found downregulation of virulence factors and metabolic genes, and upregulation of stress response pathways. The importance of the heat shock response in bacterial survival was shown as hrcA-deficient S. aureus exhibited reduced tissue colonization. These findings underscore the utility of this humanized lung model for studying S. aureus pathogenesis and bacterial adaptation to the human pulmonary environment.

Quality Evaluation of Qingwei Huanglian Pills Based on Fingerprint and Quantitative Analysis of Multi-Index Components Combined with Chemical Pattern Recognition Analysis

Qingwei Huanglian Pills (QHPs) is one of the most commonly used traditional Chinese medicine preparations for the treatment of mouth and tongue sores, but the existing quality evaluation standards have certain shortcomings and deficiencies. An effective and scientific quality evaluation method plays a vital role in medication safety. In this study, fingerprint and quantitative analysis of multi-index components combined with chemical pattern recognition analysis was used to comprehensively evaluate the quality of QHPs. The fingerprints of 15 batches of QHPs were generated and evaluated for similarity, with 10 characteristic peaks identified. Clustering hierarchical cluster analysis (HCA), principal component analysis (PCA) and orthogonal partial least squares discriminant analysis (OPLS-DA) were employed to cluster and rank the 15 batches, while simultaneously identifying the components responsible for differences between batches. The HPLC fingerprints of QHPs, along with the content determination of 10 components, were established. Twenty-eight common peaks were identified, and 10 components were specified. The similarity between the 15 batches of samples ranged from 0.983 to 0.999. Cluster analysis and comprehensive score ranking of 15 batches of samples were performed by HCA and PCA, respectively, and 13 chemical markers affecting batch differences were screened by OPLS-DA. The method established here can serve as a reference for the quality evaluation and product quality control of QHPs.

Topic modeling analysis of the Allen Human Brain Atlas

The human brain is a complex interconnected structure controlling all elementary and high-level cognitive tasks. It is composed of many regions that exhibit specific distributions of cell types and distinct patterns of functional connections. This complexity is rooted in differential transcription. The constituent cell types of different brain regions express distinctive combinations of genes as they develop and mature, ultimately shaping their functional state in adulthood. How precisely the genetic information of anatomical structures is connected to their underlying biological functions remains an open question in modern neuroscience. A major challenge is the identification of "universal patterns", which do not depend on the particular individual, but are instead basic structural properties shared by all brains. Despite the vast amount of gene expression data available at both the bulk and single-cell levels, this task remains challenging, mainly due to the lack of suitable data mining tools. In this paper, we propose an approach to address this issue based on a hierarchical version of Stochastic Block Modeling. Thanks to its specific choice of priors, the method is particularly effective in identifying these universal features. We use as a laboratory to test our algorithm a dataset obtained from six independent human brains from the Allen Human Brain Atlas. We show that the proposed method is indeed able to identify universal patterns much better than more traditional algorithms such as Latent Dirichlet Allocation or Weighted Correlation Network Analysis. The probabilistic association between genes and samples that we find well represents the known anatomical and functional brain organization. Moreover, leveraging the peculiar "fuzzy" structure of the gene sets obtained with our method, we identify examples of transcriptional and post-transcriptional pathways associated with specific brain regions, highlighting the potential of our approach.

Transcriptomic responses of lung mesenchymal cells during pneumonia

The role of mesenchymal cells during respiratory infection is not well defined, including whether, which, and how the different types of mesenchymal cells respond. We collected all mesenchymal cells from lung single-cell suspensions of mice that were naive (after receiving only saline vehicle), pneumonic (after intratracheal instillation of pneumococcus 24 hours previously), or resolved from infection (after nonlethal pneumococcal infections 6 weeks previously) and performed single-cell RNA sequencing. Cells clustered into 5 well-separated groups based on their transcriptomes: matrix fibroblasts, myofibroblasts, pericytes, smooth muscle cells, and mesothelial cells. Fibroblasts were the most abundant and could be further segregated into Pdgfra+Npnt+Ces1d+Col13a1+ alveolar fibroblasts and Cd9+Pi16+Sca1+Col14a1+ adventitial fibroblasts. The cells from naive and resolved groups overlapped in dimension reduction plots, suggesting the mesenchymal cells returned to baseline transcriptomes after resolution. During pneumonia, all mesenchymal cells responded with altered transcriptomes, revealing a core response that had been conserved across cell types as well as distinct mesenchymal cell type-specific responses. The different subsets of fibroblasts induced similar gene sets, but the alveolar fibroblasts responded more strongly than the adventitial fibroblasts. These data demonstrated diverse and specialized immune activities of lung mesenchymal cells during pneumonia.

Exploring the Interplay Between Cervicovaginal Microbiome, HPV Infection, and Cervical Intraepithelial Neoplasia in Taiwanese Women

This study examined the relationship between the vaginal microbiome, HPV infection, and cervical intraepithelial neoplasia (CIN) in 173 women. Subjects were grouped by HPV status and cervical lesion severity, ranging from HPV-negative to CIN Grade 2 or higher. Using VALENCIA classification, the study identified different community state types (CSTs) of vaginal microbiota, with CST IV subtypes (Staphylococcus dominated) showing high diversity and increased pathogenic bacteria. Results showed a higher microbial diversity in women with CIN compared to HPV-negative women, with specific microbial profiles like Enterococcus faecalis and Staphylococcus aureus correlating with HPV infection severity. LEfSe and STAMP analyses revealed that particular bacterial species, such as Fannyhessea vaginae and Sneathia vaginalis, could serve as biomarkers for HPV and CIN status. Furthermore, Streptococcus and Staphylococcus were associated with HPV genotypes α7 and α9, respectively. These findings highlight how the composition of vaginal microbiota may influence HPV persistence and CIN progression, emphasizing the need for further research into microbial dysbiosis as a potential target for cervical cancer prevention and management.

Targeting ERBB3 and AKT to overcome adaptive resistance in EML4-ALK-driven non-small cell lung cancer

The fusion event between EML4 and ALK drives a significant oncogenic activity in 5% of non-small cell lung cancer (NSCLC). Even though potent ALK-tyrosine kinase inhibitors (ALK-TKIs) are successfully used for the treatment of EML4-ALK-positive NSCLC patients, a subset of those patients eventually acquire resistance during their therapy. Here, we investigate the kinase responses in EML4-ALK V1 and V3-harbouring NSCLC cancer cells after acute inhibition with ALK TKI, lorlatinib (LOR). Using phosphopeptide chip array and upstream kinase prediction analysis, we identified a group of phosphorylated tyrosine peptides including ERBB and AKT proteins that are upregulated upon ALK-TKI treatment in EML4-ALK-positive NSCLC cell lines. Dual inhibition of ALK and ERBB receptors or AKT disrupts RAS/MAPK and AKT/PI3K signalling pathways, and enhances apoptosis in EML4-ALK + NSCLC cancer cells. Heregulin, an ERBB3 ligand, differentially modulates the sensitivity of EML4-ALK cell lines to ALK inhibitors. We found that EML4-ALK cells made resistant to LOR are sensitive to inhibition of ERBB and AKT. These findings emphasize the important roles of AKT and ERBB3 to regulate signalling after acute LOR treatment, identifying them as potential targets that may be beneficial to prevent adaptive resistance to EML4-ALK-targeted therapies in NSCLC.

PKMζ, a Brain-specific PKCζ Isoform, is Required for Glycolysis and Myofibroblastic Activation of Hepatic Stellate Cells

BACKGROUND & AIMS

Transforming growth factor (TGF)β1 induces plasma membrane (PM) accumulation of glucose transporter 1 (Glut1) required for glycolysis of hepatic stellate cells (HSCs) and HSC activation. This study aimed to understand how Glut1 is anchored/docked onto the PM of HSCs.

METHODS

HSC expression of protein kinase M zeta isoform (PKMζ) was detected by reverse transcription polymerase chain reaction (RT-PCR), Western blotting, and immunofluorescence. PKMζ level was manipulated by short hairpin RNA (shRNA) or overexpression; HSC activation was assessed by cell expression of activation markers; PM Glut1, glucose uptake, and glycolysis of HSCs were analyzed by biotinylation, 2-NBDG-based assay, and Seahorse Glycolysis Stress Test. Phospho-mutants of vasodilator-stimulated phosphorylated protein (VASP) were created by site-directed mutagenesis. TGFβ transcriptome was obtained by RNA sequencing. Single-cell RNA sequencing datasets and immunofluorescence were leveraged to analyze PKMζ expression in cancer-associated fibroblasts (CAFs) of colorectal liver metastases. Function of HSC PKMζ was determined by tumor/HSC co-implantation study.

RESULTS

Primary human and murine HSCs express PKMζ, but not full-length PKCζ. PKMζ knockdown suppresses, whereas PKMζ overexpression potentiates PM accumulation of Glut1, glycolysis, and HSC activation induced by TGFβ1. Mechanistically, PKMζ binds to and induces VASP phosphorylation at serines 157 and 239 facilitating anchoring/docking of Glut1 onto the PM of HSCs. PKMζ expression is increased in the CAFs of murine and patient colorectal liver metastases compared with quiescent HSCs. Targeting PKMζ suppresses transcriptome, CAF activation of HSCs, and colorectal tumor growth in mice.

CONCLUSIONS

Because HSCs are also a major contributor of liver fibrosis, our data highlight PKMζ and VASP as targets to inhibit metabolic reprogramming, HSC activation, liver fibrosis, and the pro-metastatic microenvironment of the liver.

Transcriptomic and Proteomic Insights into Host Immune Responses in Pediatric Severe Malarial Anemia: Dysregulation in HSP60-70-TLR2/4 Signaling and Altered Glutamine Metabolism

Severe malarial anemia (SMA, Hb < 6.0 g/dL) is a leading cause of childhood morbidity and mortality in holoendemic Plasmodium falciparum transmission zones. This study explored the entire expressed human transcriptome in whole blood from 66 Kenyan children with non-SMA (Hb ≥ 6.0 g/dL, n = 41) and SMA (n = 25), focusing on host immune response networks. RNA-seq analysis revealed 6862 differentially expressed genes, with equally distributed up-and down-regulated genes, indicating a complex host immune response. Deconvolution analyses uncovered leukocytic immune profiles indicative of a diminished antigenic response, reduced immune priming, and polarization toward cellular repair in SMA. Weighted gene co-expression network analysis revealed that immune-regulated processes are central molecular distinctions between non-SMA and SMA. A top dysregulated immune response signaling network in SMA was the HSP60-HSP70-TLR2/4 signaling pathway, indicating altered pathogen recognition, innate immune activation, stress responses, and antigen recognition. Validation with high-throughput gene expression from a separate cohort of Kenyan children (n = 50) with varying severities of malarial anemia (n = 38 non-SMA and n = 12 SMA) confirmed the RNA-seq findings. Proteomic analyses in 35 children with matched transcript and protein abundance (n = 19 non-SMA and n = 16 SMA) confirmed dysregulation in the HSP60-HSP70-TLR2/4 signaling pathway. Additionally, glutamine transporter and glutamine synthetase genes were differentially expressed, indicating altered glutamine metabolism in SMA. This comprehensive analysis underscores complex immune dysregulation and novel pathogenic features in SMA.

Loss of HD-PTP function results in lipodystrophy, defective cellular signaling and altered lipid homeostasis

His domain protein tyrosine phosphatase (HD-PTP; also known as PTPN23) facilitates function of the endosomal sorting complexes required for transport (ESCRTs) during multivesicular body (MVB) formation. To uncover its role in physiological homeostasis, embryonic lethality caused by a complete lack of HD-PTP was bypassed through generation of hypomorphic mice expressing reduced protein, resulting in animals that are viable into adulthood. These mice exhibited marked lipodystrophy and decreased receptor-mediated signaling within white adipose tissue (WAT), involving multiple prominent pathways including RAS/MAPK, phosphoinositide 3-kinase (PI3K)/AKT and receptor tyrosine kinases (RTKs), such as EGFR. EGFR signaling was dissected in vitro to assess the nature of defective signaling, revealing decreased trans-autophosphorylation and downstream effector activation, despite normal EGF binding. This corresponds to decreased plasma membrane cholesterol and increased lysosomal cholesterol, likely resulting from defective endosomal maturation necessary for cholesterol trafficking and homeostasis. The ESCRT components Vps4 and Hrs have previously been implicated in cholesterol homeostasis; thus, these findings expand knowledge on which ESCRT subunits are involved in cholesterol homeostasis and highlight a non-canonical role for HD-PTP in signal regulation and adipose tissue homeostasis.

Cell-type-specific effects of autism-associated 15q duplication syndrome in the human brain

Recurrent copy-number variation represents one of the most well-established genetic drivers in neurodevelopmental disorders, including autism spectrum disorder. Duplication of 15q11-q13 (dup15q) is a well-described neurodevelopmental syndrome that increases the risk of autism more than 40-fold. However, the effects of this duplication on gene expression and chromatin accessibility in specific cell types in the human brain remain unknown. To identify the cell-type-specific transcriptional and epigenetic effects of dup15q in the human frontal cortex, we conducted single-nucleus RNA sequencing and multi-omic sequencing on dup15q-affected individuals (n = 6) as well as individuals with non-dup15q autism (n = 7) and neurotypical control individuals (n = 7). Cell-type-specific differential expression analysis identified significantly regulated genes, critical biological pathways, and differentially accessible genomic regions. Although there was overall increased gene expression across the duplicated genomic region, cellular identity represented an important factor mediating gene-expression changes. As compared to other cell types, neuronal subtypes showed greater upregulation of gene expression across a critical region within the duplication. Genes that fell within the duplicated region and had high baseline expression in control individuals showed only modest changes in dup15q, regardless of cell type. Of note, dup15q and autism had largely distinct signatures of chromatin accessibility but shared the majority of transcriptional regulatory motifs, suggesting convergent biological pathways. However, the transcriptional binding-factor motifs implicated in each condition implicated distinct biological mechanisms: neuronal JUN and FOS networks in autism vs. an inflammatory transcriptional network in dup15q microglia. This work provides a cell-type-specific analysis of how dup15q changes gene expression and chromatin accessibility in the human brain, and it finds evidence of marked cell-type-specific effects of this genetic driver. These findings have implications for guiding therapeutic development in dup15q syndrome, as well as understanding the functional effects of copy-number variants more broadly in neurodevelopmental disorders.

Discovery and Verification of Sjögren's Syndrome Protein Biomarkers in Tears by Targeted LC-MRM

Sjögren's syndrome (SS) is an autoimmune rheumatic disorder characterized by exocrine gland dysfunction, mainly from the lacrimal and salivary glands. The disease causes severe aqueous dry eye syndrome (DED) and is associated with high rates of complications, including corneal ulceration, scaring, and perforation. Systemic complications may occur as well as a higher risk of developing lymphoma. Diagnosis of SS-DED is often delayed and difficult to establish. With the aim of discovering biomarkers to help discriminate SS-DED patients, a combination of untargeted and targeted LC-MS/MS analyses were performed on tear samples collected on Schirmer strips and subjected to tryptic digestion. Following the analysis of three cohorts and the development of two targeted LC-sMRM methods for the verification of putative biomarkers found in the first cohort of samples, 64 proteins could be linked to Sjögren's syndrome, in the hopes of helping to confirm diagnoses as well as potentially stratifying the severity of disease in these patients. Proteins that were increased in SS-DED showed activation of the immune system and alterations in homeostasis. Several proteases and protease inhibitors were found to be significantly changing in SS-DED, as well as a consistent decrease in specific proteins known to be secreted by the lacrimal gland.

Tumors Established in a Defective Immune Environment Reprogram the Oncogenic Signaling Pathways to Escalate Tumor Antigenicity

Tumors developed in immunocompromised hosts are more immunogenic. However, few studies have addressed the potential mechanisms underlying the high immunogenicity of tumors found in a suppressed immune system. Therefore, we aimed to elucidate the impacts of the immune system on tumor behaviors and immunogenicity sculpting. A murine colorectal adenocarcinoma cell line, CT26wt, was administrated into immunocompetent (BALB/c) and immunocompromised (NOD.SCID) mice, respectively. On day 11, the CT26 cells slowly progressed in the NOD.SCID mice compared to the BALB/c mice. We then performed liquid chromatography-tandem mass spectrometry (LC-MS/MS) and analyzed the differentially expressed proteins (DEPs). The DEPs participated in numerous oncogenic pathways, PI3K/AKT/mTOR cell signaling, and the silencing of several tumor suppressors, such as PTEN and RBL1, during tumorigenesis. On day 34, the CT26/SCID tumors inversely became malignant counterparts; then the CT26/SCID tumors were harvested and re-inoculated into immunocompetent mice (CT26/SCID-Re tumors) to determine the immunogenicity. The CT26/SCID-Re tumor growth rate significantly decreased. Furthermore, increased infiltrations of dendritic cells and tumor-infiltrating T lymphocytes were found in the CT26/SCID-Re tumors. These findings suggest that immunogenic tumors might express multiple tumor rejection antigens, unlike wild-type tumors, and attract more immune cells, therefore decreasing the growth rate. Collectively, our study first revealed that in immunodeficient hosts, tumor suppressors were silenced and oncogenic signaling pathways were changed during the initial phase of tumor development. With tumor progression, the tumor antigens were overexpressed, exhibiting elevated immunogenicity. This study offers a hint on the mechanisms of tumorigenesis and provides a niche for investigating the interaction between host immunity and cancer development.

Optimized bacterial community characterization through full-length 16S rRNA gene sequencing utilizing MinION nanopore technology

BACKGROUND

Accurate identification of bacterial communities is crucial for research applications, diagnostics, and clinical interventions. Although 16S ribosomal RNA (rRNA) gene sequencing is a widely employed technique for bacterial taxonomic classification, it often results in misclassified or unclassified bacterial taxa. This study sought to refine the full-length 16S rRNA gene sequencing protocol using the MinION sequencer, focusing on the V1-V9 regions. Our methodological enquiry examined several factors, including the number of PCR amplification cycles, choice of primers and Taq polymerase, and specific sequence databases and workflows employed. We used a microbial standard comprising eight bacterial strains (five gram-positive and three gram-negative) in known proportions as a validation control.

RESULTS

Based on the MinION protocol, we employed the microbial standard as the DNA template for the 16S rRNA gene amplicon sequencing procedure. Our analysis showed that an elevated number of PCR amplification cycles introduced PCR bias, and the selection of Taq polymerase and primer sets significantly affected the subsequent analysis. Bacterial identification at genus level demonstrated Pearson correlation coefficients ranging from 0.73 to 0.79 when assessed using BugSeq, Kraken-Silva and EPI2ME-16S workflows. Notably, the EPI2ME-16S workflow exhibited the highest Pearson correlation with the microbial standard, minimised misclassification, and increased alignment accuracy. At the species taxonomic level, the BugSeq workflow was superior, with a Pearson correlation coefficient of 0.92.

CONCLUSIONS

These findings emphasise the importance of careful selection of PCR settings and a well-structured analytical framework for 16S rRNA full-length gene sequencing. The results showed a robust correlation between the predicted and observed bacterial abundances at both the genus and species taxonomic levels, making these findings applicable across diverse research contexts and with clinical utility for reliable pathogen identification.

DendroX: multi-level multi-cluster selection in dendrograms

BACKGROUND

Cluster heatmaps are widely used in biology and other fields to uncover clustering patterns in data matrices. Most cluster heatmap packages provide utility functions to divide the dendrograms at a certain level to obtain clusters, but it is often difficult to locate the appropriate cut in the dendrogram to obtain the clusters seen in the heatmap or computed by a statistical method. Multiple cuts are required if the clusters locate at different levels in the dendrogram.

RESULTS

We developed DendroX, a web app that provides interactive visualization of a dendrogram where users can divide the dendrogram at any level and in any number of clusters and pass the labels of the identified clusters for functional analysis. Helper functions are provided to extract linkage matrices from cluster heatmap objects in R or Python to serve as input to the app. A graphic user interface was also developed to help prepare input files for DendroX from data matrices stored in delimited text files. The app is scalable and has been tested on dendrograms with tens of thousands of leaf nodes. As a case study, we clustered the gene expression signatures of 297 bioactive chemical compounds in the LINCS L1000 dataset and visualized them in DendroX. Seventeen biologically meaningful clusters were identified based on the structure of the dendrogram and the expression patterns in the heatmap. We found that one of the clusters consisting of mostly naturally occurring compounds is not previously reported and has its members sharing broad anticancer, anti-inflammatory and antioxidant activities.

CONCLUSIONS

DendroX solves the problem of matching visually and computationally determined clusters in a cluster heatmap and helps users navigate among different parts of a dendrogram. The identification of a cluster of naturally occurring compounds with shared bioactivities implicates a convergence of biological effects through divergent mechanisms.

Orthodontically induced changes to the genetic profile in periodontal ligament tissue and cytokine release in gingival crevicular fluid - A pilot investigation

Background/purpose

It has been known that genetic factors influence orthodontic tooth movement, however, scientific research on humans is lacking. Therefore, this study aimed to investigate dynamic changes to the genetic profile in human periodontal ligament (PDL) tissue and cytokine release in gingival crevicular fluid (GCF) during the first 28 days of orthodontic treatment.

Materials and methods

Fifteen teeth from three patients were recruited. Full-mouth fixed appliances with extraction of four premolars and one maxillary third molar was planned for orthodontic treatment. GCF collection and tooth extraction were performed following force application for 0, 1, 3, 7, and 28 days. GCF was analyzed using multiplex immunoassay for 27 cytokines. PDL tissue was collected after extraction and submitted for RNA exome-sequencing using Illumina sequencing platform. Further analysis of differentially expressed genes (DEGs), gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, and heatmaps were conducted.

Results

GCF cytokine levels varied among three patients; some patients exhibited a peak cytokine level on Day 0 whereas others did so on Days 1-3. In RNA exome sequencing data, GO and KEGG analyses showed that genes associated with sensory receptors were upregulated on Day 1, genes involved in bone remodeling were upregulated on Days 3 and 28, and genes related to osteoclast differentiation were upregulated on Day 7.

Conclusion

RNA sequencing data demonstrate that the specific types of genes are expressed at different time points, whereas the data on cytokine changes show a large variation in concentration levels and dynamic change patterns among the patients.

Protein Phosphorylation Orchestrates Acclimations of Arabidopsis Plants to Environmental pH

Environment pH (pHe) is a key parameter dictating a surfeit of conditions critical to plant survival and fitness. To elucidate the mechanisms that recalibrate cytoplasmic and apoplastic pH homeostasis, we conducted a comprehensive proteomic/phosphoproteomic inventory of plants subjected to transient exposure to acidic or alkaline pH, an approach that covered the majority of protein-coding genes of the reference plant Arabidopsis thaliana. Our survey revealed a large set-of so far undocumented pHe-dependent phospho-sites, indicative of extensive post-translational regulation of proteins involved in the acclimation to pHe. Changes in pHe altered both electrogenic H+ pumping via P-type ATPases and H+/anion co-transport processes, putatively leading to altered net trans-plasma membrane translocation of H+ ions. In pH 7.5 plants, the transport (but not the assimilation) of nitrogen via NRT2-type nitrate and AMT1-type ammonium transporters was induced, conceivably to increase the cytosolic H+ concentration. Exposure to both acidic and alkaline pH resulted in a marked repression of primary root elongation. No such cessation was observed in nrt2.1 mutants. Alkaline pH decreased the number of root hairs in the wild type but not in nrt2.1 plants, supporting a role of NRT2.1 in developmental signaling. Sequestration of iron into the vacuole via alterations in protein abundance of the vacuolar iron transporter VTL5 was inversely regulated in response to high and low pHe, presumptively in anticipation of associated changes in iron availability. A pH-dependent phospho-switch was also observed for the ABC transporter PDR7, suggesting changes in activity and, possibly, substrate specificity. Unexpectedly, the effect of pHe was not restricted to roots and provoked pronounced changes in the shoot proteome. In both roots and shoots, the plant-specific TPLATE complex components AtEH1 and AtEH2-essential for clathrin-mediated endocytosis-were differentially phosphorylated at multiple sites in response to pHe, indicating that the endocytic cargo protein trafficking is orchestrated by pHe.

Transcriptome Analysis of Sweet Cherry (Prunus avium L.) Cultivar 'Lapins' upon Infection of Pseudomonas syringae pv. syringae

Bacterial canker caused by Pseudomonas syringae pv. syringae (Pss) is responsible for substantial loss to the production of sweet cherry in Chile. To date, the molecular mechanisms of the Pss-sweet cherry interaction and the disease-related genes in the plant are poorly understood. In order to gain insight into these aspects, a transcriptomic analysis of the sweet cherry cultivar 'Lapins' for differentially expressed genes (DEGs) in response to Pss inoculation was conducted. Three Pss strains, A1M3, A1M197, and 11116_b1, were inoculated in young twigs, and RNA was extracted from tissue samples at the inoculation site and distal sections. RNA sequencing and transcriptomic expression analysis revealed that the three strains induced different patterns of responses in local and distal tissues. In the local tissues, A1M3 triggered a much more extensive response than the other two strains, enriching DEGs especially involved in photosynthesis. In the distal tissues, the three strains triggered a comparable extent of responses, among which 11116_b1 induced a group of DEGs involved in defense responses. Furthermore, tissues from various inoculations exhibited an enrichment of DEGs related to carbohydrate metabolism, terpene metabolism, and cell wall biogenesis. This study opened doors to future research on the Pss-sweet cherry interaction, immunity responses, and disease control.

Differential Protein and Glycan Packaging into Extracellular Vesicles in Response to 3D Gastric Cancer Cellular Organization

Alterations of the glycosylation machinery are common events in cancer, leading to the synthesis of aberrant glycan structures by tumor cells. Extracellular vesicles (EVs) play a modulatory role in cancer communication and progression, and interestingly, several tumor-associated glycans have already been identified in cancer EVs. Nevertheless, the impact of 3D tumor architecture in the selective packaging of cellular glycans into EVs has never been addressed. In this work, the capacity of gastric cancer cell lines with differential glycosylation is evaluated in producing and releasing EVs when cultured under conventional 2D monolayer or in 3D culture conditions. Furthermore, the proteomic content is identified and specific glycans are studied in the EVs produced by these cells, upon differential spatial organization. Here, it is observed that although the proteome of the analyzed EVs is mostly conserved, an EV differential packaging of specific proteins and glycans is found. In addition, protein-protein interaction and pathway analysis reveal individual signatures on the EVs released by 2D- and 3D-cultured cells, suggesting distinct biological functions. These protein signatures also show a correlation with clinical data. Overall, this data highlight the importance of tumor cellular architecture when assessing the cancer-EV cargo and its biological role.

Proteomic Profiling Reveals Specific Molecular Hallmarks of the Pig Claustrum

The present study, employing a comparative proteomic approach, analyzes the protein profile of pig claustrum (CLA), putamen (PU), and insula (IN). Pig brain is an interesting model whose key translational features are its similarities with cortical and subcortical structures of human brain. A greater difference in protein spot expression was observed in CLA vs PU as compared to CLA vs IN. The deregulated proteins identified in CLA resulted to be deeply implicated in neurodegenerative (i.e., sirtuin 2, protein disulfide-isomerase 3, transketolase) and psychiatric (i.e., copine 3 and myelin basic protein) disorders in humans. Metascape analysis of differentially expressed proteins in CLA vs PU comparison suggested activation of the α-synuclein pathway and L1 recycling pathway corroborating the involvement of these anatomical structures in neurodegenerative diseases. The expression of calcium/calmodulin-dependent protein kinase and dihydropyrimidinase like 2, which are linked to these pathways, was validated using western blot analysis. Moreover, the protein data set of CLA vs PU comparison was analyzed by Ingenuity Pathways Analysis to obtain a prediction of most significant canonical pathways, upstream regulators, human diseases, and biological functions. Interestingly, inhibition of presenilin 1 (PSEN1) upstream regulator and activation of endocannabinoid neuronal synapse pathway were observed. In conclusion, this is the first study presenting an extensive proteomic analysis of pig CLA in comparison with adjacent areas, IN and PUT. These results reinforce the common origin of CLA and IN and suggest an interesting involvement of CLA in endocannabinoid circuitry, neurodegenerative, and psychiatric disorders in humans.

Toxic effect and inability of L-homoserine to be a nitrogen source for growth of Escherichia coli resolved by a combination of in vivo evolution engineering and omics analyses

L-homoserine is a pivotal intermediate in the carbon and nitrogen metabolism of E. coli. However, this non-canonical amino acid cannot be used as a nitrogen source for growth. Furthermore, growth of this bacterium in a synthetic media is potently inhibited by L-homoserine. To understand this dual effect, an adapted laboratory evolution (ALE) was applied, which allowed the isolation of a strain able to grow with L-homoserine as the nitrogen source and was, at the same time, desensitized to growth inhibition by this amino acid. Sequencing of this evolved strain identified only four genomic modifications, including a 49 bp truncation starting from the stop codon of thrL. This mutation resulted in a modified thrL locus carrying a thrL* allele encoding a polypeptide 9 amino acids longer than the thrL encoded leader peptide. Remarkably, the replacement of thrL with thrL* in the original strain MG1655 alleviated L-homoserine inhibition to the same extent as strain 4E, but did not allow growth with this amino acid as a nitrogen source. The loss of L-homoserine toxic effect could be explained by the rapid conversion of L-homoserine into threonine via the thrL*-dependent transcriptional activation of the threonine operon thrABC. On the other hand, the growth of E. coli on a mineral medium with L-homoserine required an activation of the threonine degradation pathway II and glycine cleavage system, resulting in the release of ammonium ions that were likely recaptured by NAD(P)-dependent glutamate dehydrogenase. To infer about the direct molecular targets of L-homoserine toxicity, a transcriptomic analysis of wild-type MG1655 in the presence of 10 mM L-homoserine was performed, which notably identified a potent repression of locomotion-motility-chemotaxis process and of branched-chain amino acids synthesis. Since the magnitude of these effects was lower in a ΔthrL mutant, concomitant with a twofold lower sensitivity of this mutant to L-homoserine, it could be argued that growth inhibition by L-homoserine is due to the repression of these biological processes. In addition, L-homoserine induced a strong upregulation of genes in the sulfate reductive assimilation pathway, including those encoding its transport. How this non-canonical amino acid triggers these transcriptomic changes is discussed.

Transport of toxic metals in the bottom sediments and health risk assessment of Corbicula fluminea (Asiatic clam) collected from Laguna de Bay, Philippines

Laguna de Bay, the 3rd largest lake in Southeast Asia, is the most significant source of freshwater fish in the Philippines. With decades of unregulated discharge of industrial, domestic, and agricultural wastewaters into the lake, this study investigates the apportionment of heavy metals from the bottom sediments and its impact on the toxicity of Corbicula fluminea (Asiatic clam), a popular food item in the markets. The sediment samples from the western part of the lake contained higher Cd, Cu, Pb, and Zn and lower As and Cr concentrations compared to the eastern part. There were positive correlations for As, Cr, Cu, Pb, and Zn and negative correlations for Cd and Ni concentrations noted between sediments and C. fluminea. Human health risk associated with the consumption of C. fluminea collected from Laguna de Bay was attributed to the following heavy metals: Cu > As > Zn > Cd > Pb > Cr > Ni. Interestingly, the observed trend in toxicities of the shellfish was consistent with the transport phenomenon of heavy metals facilitated by the counterclockwise direction of the bottom current. This study strongly suggests the re-evaluation of the waste management plan in the industrial zones and policies regulating the sale of the shellfish harvest.

Xenograft cancer vaccines prepared from immunodeficient mice increase tumor antigen diversity and host T cell efficiency against colorectal cancers

Autologous cancer vaccines (ACVs) are a desirable approach for personalized medicine, but the efficiency of ACVs remains unsatisfactory due to their low immunogenicity. This study developed a platform that can enhance the immunogenicity of ACVs by transplanting the tumors into immunodeficient mice. The CT26 cell line was inoculated into severe combined immunodeficient mice (SCID) for vaccine preparation where escalates tumor development, subsequently diversifying the tumor antigenic topology. CT26/SCID cancer vaccines significantly inhibited tumor growth, increased the amount of tumor infiltrating lymphocytes, and triggered Th-1 predominant immune responses. Tumor antigenic profiles of CT26/SCID cells were further analyzed by liquid chromatography-tandem mass spectrometry. Compared to CT26 parental cells, a total of 428 differentially expressed proteins (DEPs) were detected. These DEPs revealed that CT26/SCID cells overexpressed several novel therapeutic targets, including KNG1, apoA-I and, β2-GPI, which can trigger cytotoxic T cells towards Th-1 predominant immune responses and directly suppress proliferation in tumors. CT26/SCID cancer vaccines can be easily manufactured, while traits of triggering stronger antigen-specific Th-1 immune activity against tumors, are retained. Results of this study provide an effective proof-of-concept of an ACV for personalized cancer immunotherapy.

HIV-1 Accessory Proteins Impart a Modest Interferon Response and Upregulate Cell Cycle-Related Genes in Macrophages

HIV-1 infection of myeloid cells is associated with the induction of an IFN response. How HIV-1 manipulates and subverts the IFN response is of key interest for the design of therapeutics to improve immune function and mitigate immune dysregulation in people living with HIV. HIV-1 accessory genes function to improve viral fitness by altering host pathways in ways that enable transmission to occur without interference from the immune response. We previously described changes in transcriptomes from HIV-1 infected and from IFN-stimulated macrophages and noted that transcription of IFN-regulated genes and genes related to cell cycle processes were upregulated during HIV-1 infection. In the present study, we sought to define the roles of individual viral accessory genes in upregulation of IFN-regulated and cell cycle-related genes using RNA sequencing. We observed that Vif induces a set of genes involved in mitotic processes and that these genes are potently downregulated upon stimulation with type-I and -II IFNs. Vpr also upregulated cell cycle-related genes and was largely responsible for inducing an attenuated IFN response. We note that the induced IFN response most closely resembled a type-III IFN response. Vpu and Nef-regulated smaller sets of genes whose transcriptomic signatures upon infection related to cytokine and chemokine processes. This work provides more insight regarding processes that are manipulated by HIV-1 accessory proteins at the transcriptional level.

Roles of Close Homologues SigB and SigD in Heat and High Light Acclimation of the Cyanobacterium Synechocystis sp. PCC 6803

Acclimation of cyanobacterium Synechocystis sp. PCC6803 to suboptimal conditions is largely dependent on adjustments of gene expression, which is highly controlled by the σ factor subunits of RNA polymerase (RNAP). The SigB and SigD σ factors are close homologues. Here we show that the sigB and sigD genes are both induced in high light and heat stresses. Comparison of transcriptomes of the control strain (CS), ΔsigB, ΔsigD, ΔsigBCE (containing SigD as the only functional group 2 σ factor), and ΔsigCDE (SigB as the only functional group 2 σ factor) strains in standard, high light, and high temperature conditions revealed that the SigB and SigD factors regulate different sets of genes and SigB and SigD regulons are highly dependent on stress conditions. The SigB regulon is bigger than the SigD regulon at high temperature, whereas, in high light, the SigD regulon is bigger than the SigB regulon. Furthermore, our results show that favoring the SigB or SigD factor by deleting other group 2 σ factors does not lead to superior acclimation to high light or high temperature, indicating that all group 2 σ factors play roles in the acclimation processes.

Make Interactive Complex Heatmaps in R

SUMMARY

Heatmap is a powerful visualization method on two-dimensional data to reveal patterns shared by subsets of rows and columns. In this work, we introduce a new R package InteractiveComplexHeatmap that brings interactivity to the widely used ComplexHeatmap package. InteractiveComplexHeatmap is designed with an easy-to-use interface where static complex heatmaps can be directly exported to an interactive Shiny web application only with one additional line of code. InteractiveComplexHeatmap also provides flexible functionalities for integrating interactive heatmap widgets to build more complex and customized Shiny web applications.

AVAILABILITY AND IMPLEMENTATION

The InteractiveComplexHeatmap package and documentations are freely available from the Bioconductor project: https://bioconductor.org/packages/InteractiveComplexHeatmap/. A complete and printer-friendly version of the documentation can also be found in Supplementary File S1.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Identification of conserved transcriptome features between humans and Drosophila in the aging brain utilizing machine learning on combined data from the NIH Sequence Read Archive

Aging is universal, yet characterizing the molecular changes that occur in aging which lead to an increased risk for neurological disease remains a challenging problem. Aging affects the prefrontal cortex (PFC), which governs executive function, learning, and memory. Previous sequencing studies have demonstrated that aging alters gene expression in the PFC, however the extent to which these changes are conserved across species and are meaningful in neurodegeneration is unknown. Identifying conserved, age-related genetic and morphological changes in the brain allows application of the wealth of tools available to study underlying mechanisms in model organisms such as Drosophila melanogaster. RNA sequencing data from human PFC and fly heads were analyzed to determine conserved transcriptome signatures of age. Our analysis revealed that expression of 50 conserved genes can accurately determine age in Drosophila (R2 = 0.85) and humans (R2 = 0.46). These transcriptome signatures were also able to classify Drosophila into three age groups with a mean accuracy of 88% and classify human samples with a mean accuracy of 69%. Overall, this work identifies 50 highly conserved aging-associated genetic changes in the brain that can be further studied in model organisms and demonstrates a novel approach to uncovering genetic changes conserved across species from multi-study public databases.

Synthesis, Characterization, and Specific Localization of Mitochondrial-Targeted Antioxidant Peptide SS31 Probe

The aim of this study is to investigate the targeting efficiency of FITC-SS31 to mitochondria in both normal and H2O2-induced oxidative damaged 661W cells, characterizing the properties of FITC-SS31 in the biological assays. The purity and molecular weight of FITC-SS31 were identified using HPLC and MS. MTT and LDH assays were used to evaluate the cytotoxicity and cell permeability. The binding ability of FITC-SS31 to cells was demonstrated by flow cytometry. The colocalization of FITC-SS31 and MitoTracker both in normal and oxidative cells was analyzed by a laser confocal microscope. We detected the DEGs between SS31+H2O2 and H2O2-alone-treated cells by RNA seq. GO and KEGG analyses were performed to predict the functional gene of SS31. The molecular weight of FITC-SS31 was 1142.2 with the 97.76% purity. The flow cytometry results showed that the MFI (mean fluorescence intensity) of FITC-SS31 in normal cells in the 4 h probe treatment group was higher than that in the 2 h and the 0 h group. The MFI in the 2 h probe treatment group was much higher than that in the 4 h and 0 h groups in damaged cells. The positive rate of 10 μM FITC-SS31 was higher than that of 1 μM and 5 μM. Fluorescein imaging analysis confirmed that FITC-SS31 was overlapped with MitoTracker. Through the analysis, DEGs were highly expressed in "localization, organelle, antioxidant activity, binding" functions and enriched in "AMPK signaling pathway, MAPK targets/nuclear events mediated by MAP kinase pathway and PI3K-Akt signaling pathway." It is speculated that SS31 exerts an antioxidant effect through one of these pathways. We hypothesized that SS31 could play a more efficient role in the pathological cells in the half-life period to avoid cell death due to oxidative damage. The functions of the DEGs in SS31+H2O2 and H2O2-alone samples are related to the localization and antioxidant activity of SS31. DEGs are mostly enriched in the AMPK signaling pathway, which needs further studies.