PathwayKO: An integrated platform for deciphering the systems-level signaling pathways.

Systems characterization of immune landscapes in health, disease and clinical intervention cases is a priority in modern medicine. High-throughput transcriptomes accumulated from gene-knockout (KO) experiments are crucial for deciphering target KO signaling pathways that are impaired by KO genes at the systems-level. There is a demand for integrative platforms. This article describes the PathwayKO platform, which has integrated state-of-the-art methods of pathway enrichment analysis, statistics analysis, and visualizing analysis to conduct cutting-edge integrative pathway analysis in a pipeline fashion and decipher target KO signaling pathways at the systems-level. We focus on describing the methodology, principles and application features of PathwayKO. First, we demonstrate that the PathwayKO platform can be utilized to comprehensively analyze real-world mouse KO transcriptomes (GSE22873 and GSE24327), which reveal systemic mechanisms underlying the innate immune responses triggered by non-infectious extensive hepatectomy (2 hours after 85% liver resection surgery) and infectious CASP-model sepsis (12 hours after CASP-model surgery). Strikingly, our results indicate that both cases hit the same core set of 21 KO MyD88-associated signaling pathways, including the Toll-like receptor signaling pathway, the NFκB signaling pathway, the MAPK signaling pathway, and the PD-L1 expression and PD-1 checkpoint pathway in cancer, alongside the pathways of bacterial, viral and parasitic infections. These findings suggest common fundamental mechanisms between these immune responses and offer informative cues that warrant future experimental validation. Such mechanisms in mice may serve as models for humans and ultimately guide formulating the research paradigms and composite strategies to reduce the high mortality rates of patients in intensive care units who have undergone successful traumatic surgical treatments. Second, we demonstrate that the PathwayKO platform model-based assessments can effectively evaluate the performance difference of pathway analysis methods when benchmarked with a collection of proper transcriptomes. Together, such advances in methods for deciphering biological insights at the systems-level may benefit the fields of bioinformatics, systems immunology and beyond.

CASP-Model Sepsis Triggers Systemic Innate Immune Responses Revealed by the Systems-Level Signaling Pathways.

Colon ascendens stent peritonitis (CASP) surgery induces a leakage of intestinal contents which may cause polymicrobial sepsis related to post-operative failure of remote multi-organs (including kidney, liver, lung and heart) and possible death from systemic syndromes. Mechanisms underlying such phenomena remain unclear. This article aims to elucidate the mechanisms underlying the CASP-model sepsis by analyzing real-world GEO data (GSE24327_A, B and C) generated from mice spleen 12 hours after a CASP-surgery in septic MyD88-deficient and wildtype mice, compared with untreated wildtype mice. Firstly, we identify and characterize 21 KO MyD88-associated signaling pathways, on which true key regulators (including ligands, receptors, adaptors, transducers, transcriptional factors and cytokines) are marked, which were coordinately, significantly, and differentially expressed at the systems-level, thus providing massive potential biomarkers that warrant experimental validations in the future. Secondly, we observe the full range of polymicrobial (viral, bacterial, and parasitic) sepsis triggered by the CASP-surgery by comparing the coordinated up- or down-regulations of true regulators among the experimental treatments born by the three data under study. Finally, we discuss the observed phenomena of "systemic syndrome", "cytokine storm" and "KO MyD88 attenuation", as well as the proposed hypothesis of "spleen-mediated immune-cell infiltration". Together, our results provide novel insights into a better understanding of innate immune responses triggered by the CASP-model sepsis in both wildtype and MyD88-deficient mice at the systems-level in a broader vision. This may serve as a model for humans and ultimately guide formulating the research paradigms and composite strategies for the early diagnosis and prevention of sepsis.

CRISPR-Cas9 Based Bacteriophage Genome Editing.

Bacteriophages are the most abundant entities in the biosphere, and many genomes of rare and novel bacteriophages have been sequenced to date. However, bacteriophage functional genomics has been limited by a lack of effective research methods. Clustered regularly interspaced short palindromic repeat/CRISPR-associated gene (CRISPR-Cas) systems provide bacteriophages with a new mechanism for attacking host bacteria as well as new tools for study bacteriophage functional genomics. It has been reported that bacteriophages are not only the driving elements of the evolution of prokaryote CRISPR arrays but also the targets of CRISPR-Cas systems. In this study, a phage genome editing platform based on the heterologous CRISPR-Cas9 system was theoretically designed, and a Vibrio natriegens phage TT4P2 genome editing experiment was carried out in vivo in the host bacterium Vibrio natriegens TT4 to achieve phage gene deletion and replacement. The construction of this phage genome editing platform is expected to solve the problem of insufficient research on phage gene diversity, promote the development of phage synthetic biology and nanotechnology, and even accelerate the discovery of new molecular biology tools. IMPORTANCE Bacteriophages are the most numerous organisms on earth and are known for their diverse lifestyles. Since the discovery of bacteriophages, our knowledge of the wider biological world has undergone immense and unforeseen changes. A variety of V. natriegens phages have been detected, but few have been well characterized. CRISPR was first documented in Escherichia coli in 1987. It has been reported that the CRISPR-Cas system can target and cleave invaders, including bacteriophages, in a sequence-specific manner. Here, we show that the construction of a phage genome editing platform based on the heterologous CRISPR-Cas9 system can achieve V. natriegens phage TT4P2 gene editing and can also improve the efficiency and accuracy of phage TT4P2 gene editing.

MIMRDA: A Method Incorporating the miRNA and mRNA Expression Profiles for Predicting miRNA-Disease Associations to Identify Key miRNAs (microRNAs).

Identifying cancer-related miRNAs (or microRNAs) that precisely target mRNAs is important for diagnosis and treatment of cancer. Creating novel methods to identify candidate miRNAs becomes an imminent Frontier of researches in the field. One major obstacle lies in the integration of the state-of-the-art databases. Here, we introduce a novel method, MIMRDA, which incorporates the miRNA and mRNA expression profiles for predicting miRNA-disease associations to identify key miRNAs. As a proof-of-principle study, we use the MIMRDA method to analyze TCGA datasets of 20 types (BLCA, BRCA, CESE, CHOL, COAD, ESCA, HNSC, KICH, KIRC, KIRP, LIHC, LUAD, LUSC, PAAD, PRAD, READ, SKCM, STAD, THCA and UCEC) of cancer, which identified hundreds of top-ranked miRNAs. Some (as Category 1) of them are endorsed by public databases including TCGA, miRTarBase, miR2Disease, HMDD, MISIM, ncDR and mTD; others (as Category 2) are supported by literature evidences. miR-21 (representing Category 1) and miR-1258 (representing Category 2) display the excellent characteristics of biomarkers in multi-dimensional assessments focusing on the function similarity analysis, overall survival analysis, and anti-cancer drugs' sensitivity or resistance analysis. We compare the performance of the MIMRDA method over the Limma and SPIA packages, and estimate the accuracy of the MIMRDA method in classifying top-ranked miRNAs via the Random Forest simulation test. Our results indicate the superiority and effectiveness of the MIMRDA method, and recommend some top-ranked key miRNAs be potential biomarkers that warrant experimental validations.

Complementary Transcriptomic and Proteomic Analysis in the Substantia Nigra of Parkinson's Disease.

Parkinson's disease (PD) is a disease that involves brain damage and is associated with neuroinflammation, mitochondrial damage, and cell aging. However, the pathogenic mechanism of PD is still unknown. Sequencing data and proteomic data can describe the fluctuation of molecular abundance in diseases at the mRNA level and protein level, respectively. In order to explore new targets in the pathogenesis of PD, the study analyzed molecular changes from the database by combining transcriptomic and proteomic analysis. Differentially expressed genes and differentially abundant proteins were summarized and analyzed. Enrichment and cluster analysis emphasized the importance of neurotransmitter release, mitochondrial damage, and vesicle transport. The molecular network revealed a subnetwork of 9 molecules related to SCNA and TH and revealed hub gene with differential expression at both mRNA and protein levels. It found that ACHE and CADPS could be used as new targets in PD, emphasizing that impaired nerve signal transmission and vesicle transport affect the pathogenesis of PD. Our research emphasized that the joint analysis and verification of transcriptomics and proteomics were devoted to understanding the comprehensive views and mechanism of pathogenesis in PD.

Effects of Nrf-2 expression in trophoblast cells and vascular endothelial cells in preeclampsia.

The present study aimed to explore the role of kelch-like ECH-associated protein-1 (Keap1)/Nuclear factor erythroid 2-related factor 2 (Nrf-2) signaling pathway in regulating heme oxygenase-1 (HO-1) expression in adverse outcomes of preeclampsia (PE). Adult Wistar rats, HTR-8/SVneo and hESC cells were used for models in vitro and in vivo, respectively. Inhibition of Nrf-2 could slightly reduce the elevation of systolic blood pressure (SBP) and urinary protein in PE rats. The percentages of dead fetuses during pregnancy and within seven days of birth were decreased by Nrf-2 inhibitor. There was no significant effect on the pathology and HO-1 expression of Nrf-2 in placental tissue. Deficiency of Nrf-2 increased significantly the levels of chemokine 2 (CCL2), interleukin-1ß (IL-1ß), tumor necrosis factor-alpha (TNF-α), angiotensin II receptor type 1 (AT1R) and reactive oxygen species (ROS) in the embryonic tissues. Knockdown of Nrf-2 suppressed cell proliferation, improved cell apoptosis and invasion with an increase of ROS and HO-1, but the effect on cells apoptosis was greater. Activation of Nrf-2 pathway could reduce oxidative stress in PE rats and trophoblast cells induced by Ang II, and enhance the adverse outcome of PE via increasing HO-1. Nrf-2 silence reshaped blood vessels and achieved the effect of treating PE. Our results might provide theoretical guidance for the application of Nrf-2 in the treatment of PE.

Long non-coding RNA LINC01194 promotes the proliferation, migration and invasion of lung adenocarcinoma cells by targeting miR-641/SETD7 axis.

BACKGROUND: It is increasingly evidenced that long non-coding RNAs (lncRNAs) play an important role in various diseases. LncRNA LINC01194 acts as an oncogene in several cancer types. Nevertheless, the role of LINC01194 in lung adenocarcinoma (LUAD) has not yet been revealed. METHODS: qRT-PCR was used to detect the expression of LINC01194, miR-641 and SETD7 mRNA, while western blot was exploited to examine SETD7 protein level. Cell proliferation was detected by colony formation and EdU assays. Transwell assays detected cell migration and invasion. TUNEL assay and flow cytometry analysis were used to detect cell apoptosis. RIP, RNA pull down and luciferase reporter assays detected the binding among LINC01194, miR-641 and SETD7. RESULTS: LINC01194 was significantly upregulated in LUAD tissues and cell lines. Knockdown of LINC01194 resulted in decreased cell proliferation, migration and invasion, and increased apoptosis. Mechanistic experiments unveiled that LINC01194 augmented SETD7 expression in LUAD cells by competitively interacting with miR-641. Rescue experiments showed that miR-641 inhibition and SETD7 overexpression rescued the repressing impacts on LUAD cell proliferation, migration and invasion caused by LINC01194 knockdown. CONCLUSION: LINC01194 promotes the progression of LUAD by enhancing miR-641-targeted SETD7. The LINC01194/miR-641/SETD7 axis might provide new molecular targets for treating LUAD.

Pseudomonas sp. UW4 acdS gene promotes primordium initiation and fruiting body development of Agaricus bisporus.

To simplify industrial mushroom cultivation, we introduced a bacterial Pseudomonas sp. UW4 acdS gene, encoding 1-aminocyclopropane-1-carboxylic acid (ACC) deaminase (AcdS), into fungus Agaricus bisporus. Transformant A. bisporus-acdS14 cased with sterilized-vermiculite generated primordia 5 days sooner than wild-type strain, confirming the specific role of the AcdS enzyme. Being consistent with the AcdS enzyme activity increased by 84%, the mycelium growth rate was increased by 25%; but, the ACC and ethylene concentrations were reduced by 71% and 36%, respectively, in the A. bisporus-acdS14 transformant. And the bacterium P. sp. UW4 attachment on the mycelium of the A. bisporus-acdS14 transformant was drastically reduced. We conclude that the heterogeneously expressed bacterial acdS gene degrades ACC and reduces ethylene-synthesis, eliminating ethylene inhibition on the mycelium growth and primordium formation in A. bisporus. Our results provide new insights into the mechanism underlying casing soil bacterium, and help formulate a casing-less cultivation for the next-generation mushroom industry.

NT5E is associated with unfavorable prognosis and regulates cell proliferation and motility in gastric cancer.

Ecto-5'-nucleotidase (NT5E) is a glycosylphosphatidylinositol anchored cell surface protein, and has been suggested to be dysregulated in most types of human cancer including gastric cancer. The aim of the present study was to present more evidence about the clinical and prognostic value of Ecto-5'-nucleotidase in gastric cancer patients, and preliminarily explore the biological function of Ecto-5'-nucleotidase in gastric cancer cells. In our study, high Ecto-5'-nucleotidase expression was observed in gastric cancer tissues and cell lines, respectively, compared with normal gastric mucosa tissues cells. Meanwhile, TCGA database also indicated that Ecto-5'-nucleotidase expression levels were notably elevated in gastric cancer tissues compared with normal gastric mucosa tissues. Furthermore, high-expression of Ecto-5'-nucleotidase was obviously associated with advanced clinical stage, deep tumor invasion, lymph node metastasis and distant metastasis in gastric cancer patients. The survival analyses of TCGA database and our study consistent suggested high Ecto-5'-nucleotidase expression was negatively correlated with overall survival time in gastric cancer patients. The univariate and multivariate Cox proportional hazards regression model showed high Ecto-5'-nucleotidase expression was an independent poor prognostic factor for gastric cancer patients. Moreover, silencing of Ecto-5'-nucleotidase expression suppressed cell proliferation, migration and invasion in vitro in gastric cancer. In conclusion, Ecto-5'-nucleotidase is a credible prognostic biomarker, and serves as a potential therapeutic target in gastric cancer.

Dust mite-derived Der f 3 activates a pro-inflammatory program in airway epithelial cells via PAR-1 and PAR-2.

Protease activity of allergens has been suggested to be involved in the pathogenesis of allergic diseases. The major allergen Der f 3 from Dermatophagoides farinae harbors serine protease activity, but its immunopathogenesis remains unclear. This study aims to explore the effect of Der f 3 on the airway epithelial barrier and on the molecular pathways by which Der f 3 induces inflammation. RNA-seq was performed to identify differentially expressed genes in bronchial airway epithelial cells (AEC) between native Der f 3 and heat-inactivated (H) Der f 3, coupled with real-time PCR (RT-PCR) and ELISA for validation. Unlike other protease allergens such as that induce Th2-promoting alarmins (IL-25, IL-33, TSLP) in AECs, Der f 3 induced pro-inflammatory cytokines and chemokines including IL-6, IL-8 and GM-CSF, which are known to promote Th17 response. These pro-inflammatory mediators were induced by Der f 3 via the MAPK and NF-κB pathways as well as the store-operated calcium signaling. Gene silencing with small interfering RNA in A549 and BEAS-2B cells indicated that activation of AECs by Der f 3 was mainly dependent on protease-activated receptor 2 (PAR-2), while PAR-1 was also required for the full activation of AECs. Double knock-down of PAR-1 and PAR-2 largely impaired Der f 3-inducecd IL-8 production and subsequent signaling pathways. Our data suggest that Der f 3 induces pro-inflammatory mediators in human epithelial cell lines via the PARs-MAPK-NF-κB axis. Our results provide a molecular mechanism by which Der f 3 may trigger the Th17-skewed allergic response toward house dust mites.

Identification of DNA methylation signature to predict prognosis in gastric adenocarcinoma.

Gastric adenocarcinoma is an important death-related cancer. To find factors related to survival and prognosis, and thus improve recovery prospects, a powerful signature is needed. DNA methylation plays an important role in gastric adenocarcinoma processes and development, and here we report on the search for a significant DNA methylation gene to aid with the earlier diagnosis of gastric adenocarcinoma patients. A Cox proportional risk regression analysis and random survival forest algorithm were used to analyze gastric adenocarcinoma patients' DNA methylation data from The Cancer Genome Atlas, a public database. DNA methylation gene signature consisting of five genes (SERPINA3, AP000357.4, GZMA, AC004702.2, and GREB1L) were selected. As the most accurate predictor, the area under the curve in the training and test group were 0.72 and 0.61, respectively. The signature was able to sort patients into high- and low-risk groups with meaningful overall survival rates (median: 18.36 vs 72.23 months, log-rank test, P < 0.001) in the training group, which predictive ability was validated in a test data set (median: 25.56 vs 58.80 months, log-rank test, P < 0.016). A multivariate Cox regression analysis showed the significant DNA methylation was an independent prediction prognostic factor for gastric adenocarcinoma patients. Functional analysis suggests that these signature genes may be related to pathways and biological processes associated with tumorigenesis. The significant DNA methylation gene could be a novel prediction and prognostic biomarker that both aids in the treatment and predicts the overall survival likelihoods of gastric adenocarcinoma patients.

Exploration and optimization of mixed acid synergistic catalysis pretreatment for maximum C5 sugars.

The liquid hot water (LHW) pretreatment could be strengthened by acetic and lactic acids produced from the process. The synergistic effect of the mixed acid catalyst of lactic acid and acetic acid was investigated for the purpose of maximization of the overall C5 sugars yield. Individual acids (acetic and lactic acid) and mixed acid were used to strengthen the LHW pretreatment at different conditions. The results showed that the suitable conditions of mixed acid synergistic catalysis was at 180 °C for 60 min and 3 wt% mixed acid where the ratio of 40% (i.e. 0.40 in mass fraction of lactic acid in mixed acid). Response surface methodology (RSM) was applied to further optimize this process. The highest yield of C5 sugars of 93.83% according to theoretical predicted model, was close to the experiment value of 92.53% at 177 °C for 67 min and with the ratio of mixed acid of 40%.

GenomeLandscaper: Landscape analysis of genome-fingerprints maps assessing chromosome architecture.

Assessing correctness of an assembled chromosome architecture is a central challenge. We create a geometric analysis method (called GenomeLandscaper) to conduct landscape analysis of genome-fingerprints maps (GFM), trace large-scale repetitive regions, and assess their impacts on the global architectures of assembled chromosomes. We develop an alignment-free method for phylogenetics analysis. The human Y chromosomes (GRCh.chrY, HuRef.chrY and YH.chrY) are analysed as a proof-of-concept study. We construct a galaxy of genome-fingerprints maps (GGFM) for them, and a landscape compatibility among relatives is observed. But a long sharp straight line on the GGFM breaks such a landscape compatibility, distinguishing GRCh38p1.chrY (and throughout GRCh38p7.chrY) from GRCh37p13.chrY, HuRef.chrY and YH.chrY. We delete a 1.30-Mbp target segment to rescue the landscape compatibility, matching the antecedent GRCh37p13.chrY. We re-locate it into the modelled centromeric and pericentromeric region of GRCh38p10.chrY, matching a gap placeholder of GRCh37p13.chrY. We decompose it into sub-constituents (such as BACs, interspersed repeats, and tandem repeats) and trace their homologues by phylogenetics analysis. We elucidate that most examined tandem repeats are of reasonable quality, but the BAC-sized repeats, 173U1020C (176.46 Kbp) and 5U41068C (205.34 Kbp), are likely over-repeated. These results offer unique insights into the centromeric and pericentromeric regions of the human Y chromosomes.

Kinetic studies of the strengthening effect on liquid hot water pretreatments by organic acids.

The liquid hot water (LHW) pretreatments would be accelerated by the organic acids produced from the process. In the study, the organic acids included not only acetic acid but also lactic acid during LHW hydrolysis of reeds, at 180-220°C and for 15-135min. The lactic acid was presumably produced from xylose degradation in the pretreatment process. The different organic acids, such as acetic acid, lactic acid and acetic-lactic acids, were used to strengthen the LHW pretreatments for increasing xylose production. Moreover, the work presented kinetic models of xylose and hemicellulose at different conditions, considering the generation of lactic acid. The experimental and kinetic results both indicated that acetic-lactic acids had synergistic catalytic effect on the reaction, which could not only inhibit the degradation of xylose, but also promote the hydrolysis of hemicellulose. Besides, the highest concentration of xylose of 7.323g/L was obtained at 200°C, for 45min and with 1wt% acetic-lactic acids.

GenomeFingerprinter: the genome fingerprint and the universal genome fingerprint analysis for systematic comparative genomics.

BACKGROUND: No attention has been paid on comparing a set of genome sequences crossing genetic components and biological categories with far divergence over large size range. We define it as the systematic comparative genomics and aim to develop the methodology. RESULTS: First, we create a method, GenomeFingerprinter, to unambiguously produce a set of three-dimensional coordinates from a sequence, followed by one three-dimensional plot and six two-dimensional trajectory projections, to illustrate the genome fingerprint of a given genome sequence. Second, we develop a set of concepts and tools, and thereby establish a method called the universal genome fingerprint analysis (UGFA). Particularly, we define the total genetic component configuration (TGCC) (including chromosome, plasmid, and phage) for describing a strain as a systematic unit, the universal genome fingerprint map (UGFM) of TGCC for differentiating strains as a universal system, and the systematic comparative genomics (SCG) for comparing a set of genomes crossing genetic components and biological categories. Third, we construct a method of quantitative analysis to compare two genomes by using the outcome dataset of genome fingerprint analysis. Specifically, we define the geometric center and its geometric mean for a given genome fingerprint map, followed by the Euclidean distance, the differentiate rate, and the weighted differentiate rate to quantitatively describe the difference between two genomes of comparison. Moreover, we demonstrate the applications through case studies on various genome sequences, giving tremendous insights into the critical issues in microbial genomics and taxonomy. CONCLUSIONS: We have created a method, GenomeFingerprinter, for rapidly computing, geometrically visualizing, intuitively comparing a set of genomes at genome fingerprint level, and hence established a method called the universal genome fingerprint analysis, as well as developed a method of quantitative analysis of the outcome dataset. These have set up the methodology of systematic comparative genomics based on the genome fingerprint analysis.

Time-dependent fermentation control strategies for enhancing synthesis of marine bacteriocin 1701 using artificial neural network and genetic algorithm.

The artificial neural network (ANN) and genetic algorithm (GA) were combined to optimize the fermentation process for enhancing production of marine bacteriocin 1701 in a 5-L-stirred-tank. Fermentation time, pH value, dissolved oxygen level, temperature and turbidity were used to construct a "5-10-1" ANN topology to identify the nonlinear relationship between fermentation parameters and the antibiotic effects (shown as in inhibition diameters) of bacteriocin 1701. The predicted values by the trained ANN model were coincided with the observed ones (the coefficient of R(2) was greater than 0.95). As the fermentation time was brought in as one of the ANN input nodes, fermentation parameters could be optimized by stages through GA, and an optimal fermentation process control trajectory was created. The production of marine bacteriocin 1701 was significantly improved by 26% under the guidance of fermentation control trajectory that was optimized by using of combined ANN-GA method.

[Advances of development of phage display systems].

The phage display technology (PDT) was unique in genetic engineering and recombinant expression. The phage display systems (PDS) were platforms (kits) composed of genetic modified phages, helper phages, and host bacteria. This review concisely summarized the development of four types of PDS, based on M13, λ, T4, and T7 phages, in terms of phage molecular genetics and genetic (gene or genome) engineering. We addressed on the key components and their genetic (genomic) engineering for modifications, the technical features of different anchors, and the development progress and selection reference of those different kits.