cfSNV: a software tool for the sensitive detection of somatic mutations from cell-free DNA.

Cell-free DNA (cfDNA) in blood, viewed as a surrogate for tumor biopsy, has many clinical applications, including diagnosing cancer, guiding cancer treatment and monitoring treatment response. All these applications depend on an indispensable, yet underdeveloped task: detecting somatic mutations from cfDNA. The task is challenging because of the low tumor fraction in cfDNA. Recently, we developed the computational method cfSNV, the first method that comprehensively considers the properties of cfDNA for the sensitive detection of mutations from cfDNA. cfSNV vastly outperformed the conventional methods that were developed primarily for calling mutations from solid tumor tissues. cfSNV can accurately detect mutations in cfDNA even with medium-coverage (e.g., ≥200×) sequencing, which makes whole-exome sequencing (WES) of cfDNA a viable option for various clinical utilities. Here, we present a user-friendly cfSNV package that exhibits fast computation and convenient user options. We also built a Docker image of it, which is designed to enable researchers and clinicians with a limited computational background to easily carry out analyses on both high-performance computing platforms and local computers. Mutation calling from a standard preprocessed WES dataset (~250× and ~70 million base pair target size) can be carried out in 3 h on a server with eight virtual CPUs and 32 GB of random access memory.

Cost-effective methylome sequencing of cell-free DNA for accurately detecting and locating cancer.

Early cancer detection by cell-free DNA faces multiple challenges: low fraction of tumor cell-free DNA, molecular heterogeneity of cancer, and sample sizes that are not sufficient to reflect diverse patient populations. Here, we develop a cancer detection approach to address these challenges. It consists of an assay, cfMethyl-Seq, for cost-effective sequencing of the cell-free DNA methylome (with > 12-fold enrichment over whole genome bisulfite sequencing in CpG islands), and a computational method to extract methylation information and diagnose patients. Applying our approach to 408 colon, liver, lung, and stomach cancer patients and controls, at 97.9% specificity we achieve 80.7% and 74.5% sensitivity in detecting all-stage and early-stage cancer, and 89.1% and 85.0% accuracy for locating tissue-of-origin of all-stage and early-stage cancer, respectively. Our approach cost-effectively retains methylome profiles of cancer abnormalities, allowing us to learn new features and expand to other cancer types as training cohorts grow.

Downregulation of Methionine Cycle Genes MAT1A and GNMT Enriches Protein-Associated Translation Process and Worsens Hepatocellular Carcinoma Prognosis.

The major biological methyl donor, S-adenosylmethionine (adoMet) synthesis occurs mainly in the liver. Methionine adenosyltransferase 1A (MAT1A) and glycine N-methyltransferase (GNMT) are two key enzymes involved in the functional implications of that variation. We collected 42 RNA-seq data from paired hepatocellular carcinoma (HCC) and its adjacent normal liver tissue from the Cancer Genome Atlas (TCGA). There was no mutation found in MAT1A or GNMT RNA in the 42 HCC patients. The 11,799 genes were annotated in the RNA-Seq data, and their expression levels were used to investigate the phenotypes of low MAT1A and low GNMT by Gene Set Enrichment Analysis (GSEA). The REACTOME_TRANSLATION gene set was enriched and visualized in a heatmap along with corresponding differences in gene expression between low MAT1A versus high MAT1A and low GNMT versus high GNMT. We identified 43 genes of the REACTOME_TRANSLATION gene set that are powerful prognosis factors in HCC. The significantly predicted genes were referred into eukaryotic translation initiation (EIF3B, EIF3K), eukaryotic translation elongation (EEF1D), and ribosomal proteins (RPs). Cell models expressing various MAT1A and GNMT proved that simultaneous restoring the expression of MAT1A and GNMT decreased cell proliferation, invasion, as well as the REACTOME_TRANSLATION gene EEF1D, consistent with a better prognosis in human HCC. We demonstrated new findings that downregulation or defect in MAT1A and GNMT genes can enrich the protein-associated translation process that may account for poor HCC prognosis. This is the first study demonstrated that MAT1A and GNMT, the 2 key enzymes involved in methionine cycle, could attenuate the function of ribosome translation. We propose a potential novel mechanism by which the diminished GNMT and MAT1A expression may confer poor prognosis for HCC.

Genomic characterization of a rare Carica papaya X chromosome mutant reveals a candidate monodehydroascorbate reductase 4 gene involved in all-hermaphrodite phenomenon.

Sex form is one of the most important characteristics in papaya cultivation in which hermaphrodite is the preferable form. Self-pollination of H*-TSS No.7, an inbred line derived from a rare X chromosome mutant SR*, produced all-hermaphrodite progeny. The recessive lethal allele controlling the all-hermaphrodite phenomenon was proposed to be the recessive Germination suppressor (gs) locus. This study employed next-generation sequencing technology and genome comparison to identify the candidate Gs gene. One specific gene, monodehydroascorbate reductase 4 (MDAR4) harboring a unique polymorphic 3 bp deletion in H*-TSS No.7 was identified. The function of MDAR4 is known to be involved in the hydrogen peroxide (H2O2) scavenging pathway and is associated with seed germination. Furthermore, MDAR4 showed higher expression in the imbibed seeds than that in the dry seeds indicating its potential role in the seed germination. Perhaps this is the very first report providing the evidences that MDAR4 is the candidate of Gs locus in H*-TSS No.7. In addition, Gs allele-specific markers were developed which would be facilitated for breeding all-hermaphrodite lines.

MTHFR Knockdown Assists Cell Defense against Folate Depletion Induced Chromosome Segregation and Uracil Misincorporation in DNA.

Folate depletion causes chromosomal instability by increasing DNA strand breakage, uracil misincorporation, and defective repair. Folate mediated one-carbon metabolism has been suggested to play a key role in the carcinogenesis and progression of hepatocellular carcinoma (HCC) through influencing DNA integrity. Methylenetetrahydrofolate reductase (MTHFR) is the enzyme catalyzing the irreversible conversion of 5,10-methylenetetrahydrofolate to 5-methyltetrahydrofolate that can control folate cofactor distributions and modulate the partitioning of intracellular one-carbon moieties. The association between MTHFR polymorphisms and HCC risk is inconsistent and remains controversial in populational studies. We aimed to establish an in vitro cell model of liver origin to elucidate the interactions between MTHFR function, folate status, and chromosome stability. In the present study, we (1) examined MTHFR expression in HCC patients; (2) established cell models of liver origin with stabilized inhibition of MTHFR using small hairpin RNA delivered by a lentiviral vector, and (3) investigated the impacts of reduced MTHFR and folate status on cell cycle, methyl group homeostasis, nucleotide biosynthesis, and DNA stability, all of which are pathways involved in DNA integrity and repair and are critical in human tumorigenesis. By analyzing the TCGA/GTEx datasets available within GEPIA2, we discovered that HCC cancer patients with higher MTHFR had a worse survival rate. The shRNA of MTHFR (shMTHFR) resulted in decreased MTHFR gene expression, MTHFR protein, and enzymatic activity in human hepatoma cell HepG2. shMTHFR tended to decrease intracellular S-adenosylmethionine (SAM) contents but folate depletion similarly decreased SAM in wildtype (WT), negative control (Neg), and shMTHFR cells, indicating that in cells of liver origin, shMTHFR does not exacerbate the methyl group supply in folate depletion. shMTHFR caused cell accumulations in the G2/M, and cell population in the G2/M was inversely correlated with MTHFR gene level (r = -0.81, p < 0.0001), MTHFR protein expression (r = -0.8; p = 0.01), and MTHFR enzyme activity (r = -0.842; p = 0.005). Folate depletion resulted in G2/M cell cycle arrest in WT and Neg but not in shMTHFR cells, indicating that shMTHFR does not exacerbate folate depletion-induced G2/M cell cycle arrest. In addition, shMTHFR promoted the expression and translocation of nuclei thymidine synthetic enzyme complex SHMT1/DHFR/TYMS and assisted folate-dependent de novo nucleotide biosynthesis under folate restriction. Finally, shMTHFR promoted nuclear MLH1/p53 expression under folate deficiency and further reduced micronuclei formation and DNA uracil misincorporation under folate deficiency. In conclusion, shMTHFR in HepG2 induces cell cycle arrest in G2/M that may promote nucleotide supply and assist cell defense against folate depletion-induced chromosome segregation and uracil misincorporation in the DNA. This study provided insight into the significant impact of MTHFR function on chromosome stability of hepatic tissues. Data from the present study may shed light on the potential regulatory mechanism by which MTHFR modulates the risk for hepatic malignancies.

Pharmacophore-based virtual screening for the identification of the novel Src inhibitor SJG-136 against lung cancer cell growth and motility.

Aberrant elevated Src activity is related to lung cancer growth and metastasis. Therefore, the development of potent small molecule inhibitors to target Src kinase is a potential therapeutic strategy for lung cancer. This study aimed to develop a computational model for the in silico screening of Src inhibitors and then assess the suppressive effect of candidate compounds on cellular functions. A 3D-quantitative structure-activity relationship (QSAR) pharmacophore model consisting of two hydrogen bond acceptors and two hydrophobic regions was constructed by using 28 structurally diverse compounds with IC50 values spanning four orders of magnitude. A National Cancer Institute (NCI) compound dataset was employed for virtual screening by applying the pharmacophore model and molecular docking. Candidate compounds were chosen from the top 20% of scored hits. Among these compounds, the suppressive effects of 30 compounds available in the NCI on Src phosphorylation were validated by using an enzyme-linked immunosorbent assay. Among these compounds, SJG-136, a pyrrolobenzodiazepine dimer, showed a significant inhibitory effect against Src activity in a dose-dependent manner. Further investigations showed that SJG-136 can inhibit lung cancer cell proliferation, clonogenicity, invasion and migration in vitro and tumour growth in vivo. Furthermore, SJG-136 also had an inhibitory effect on Src-related signaling pathways, including the FAK, paxillin, p130Cas, PI3K, AKT, and MEK pathways. In conclusion, we have established a pharmacophore-based virtual screening approach to identify novel Src inhibitors that can inhibit lung cancer cell growth and motility through suppressing Src-related pathways. These findings may contribute to the development of targeted drugs for lung cancer treatment, such as lead compounds.

Spatial Distribution of Heavy Metals and Associated Risks in Sediment of the Urban River Flowing into the Pearl River Estuary, China.

Twenty-eight transects in six urban rivers located in the Qianhai area of Baoan District, Shenzhen, China, were selected to study the contamination of heavy metals in its sediment, reveal the spatial distribution of heavy metals in these rivers, and assess the ecological risk of heavy metals in the urban river sediment. The result showed that the concentration of heavy metals displayed considerable variation in the study area. The results showed that southern airport drainage river had exposure to high concentrations of Cu, Nanchang Chung exposure to high concentrations of Zn, Cr, and Pb, Gongle Chung exposure to high concentrations of Ni, and Xixiang River exposure to high concentrations of Hg. The results of correlation analysis and cluster analysis indicated that Cu had significant correlations with Zn and Cr respectively (P < 0.05). There was significant correlation, and the same cluster existed among Zn, Cr, and Pb (P < 0.01), whereas Hg was not significantly correlated with other metals. The potential risk of the six rivers could be sorted as Gushu Chung < Gongle Chung < Tiegang reservoir flood discharge river < Xixiang River (a), Nanchang Chung < and Southern Airport drainage river.

Metabolic Pathways Enhancement Confers Poor Prognosis in p53 Exon Mutant Hepatocellular Carcinoma.

RNA-Sequencing (RNA-Seq), the most commonly used sequencing application tool, is not only a method for measuring gene expression but also an excellent media to detect important structural variants such as single nucleotide variants (SNVs), insertion/deletion (Indels), or fusion transcripts. The Cancer Genome Atlas (TCGA) contains genomic data from a variety of cancer types and also provides the raw data generated by TCGA consortium. p53 is among the top 10 somatic mutations associated with hepatocellular carcinoma (HCC). The aim of the present study was to analyze concordant different gene profiles and the priori defined set of genes based on p53 mutation status in HCC using RNA-Seq data. In the study, expression profile of 11 799 genes on 42 paired tumor and adjacent normal tissues was collected, processed, and further stratified by the mutated versus normal p53 expression. Furthermore, we used a knowledge-based approach Gene Set Enrichment Analysis (GSEA) to compare between normal and p53 mutation gene expression profiles. The statistical significance (nominal P value) of the enrichment score (ES) genes was calculated. The ranked gene list that reflects differential expression between p53 wild-type and mutant genotypes was then mapped to metabolic process by KEGG, an encyclopedia of genes and genomes to assign functional meanings. These approaches enable us to identify pathways and potential target gene/pathways that are highly expressed in p53 mutated HCC. Our analysis revealed 2 genes, the hexokinase 2 (HK2) and Enolase 1 (ENO1), were conspicuous of red pixel in the heatmap. To further explore the role of these genes in HCC, the overall survival plots by Kaplan-Meier method were performed for HK2 and ENO1 that revealed high HK2 and ENO1 expression in patients with HCC have poor prognosis. These results suggested that these glycolysis genes are associated with mutated-p53 in HCC that may contribute to poor prognosis. In this proof-of-concept study, we proposed an approach for identifying novel potential therapeutic targets in human HCC with mutated p53. These approaches can take advantage of the massive next-generation sequencing (NGS) data generated worldwide and make more out of it by exploring new potential therapeutic targets.

An Association of Gut Microbiota with Different Phenotypes in Chinese Patients with Rheumatoid Arthritis.

We aimed to investigate the association of gut microbiota with disease activity, inflammatory parameters, and auto-antibodies profile in rheumatoid arthritis (RA). A total of 138 RA patients and 21 healthy controls (HC) were enrolled. Fecal samples were collected for bacterial DNA extraction and 16S ribosome (r)RNA sequencing, followed by analyses of gut microbiota composition. Serum levels of tumor necrosis factor (TNF)-α, interleukin (IL)-6, and IL-17A were determined by using ELISA. Our results indicated that RA patients had lower diversity index, which reflects both evenness and richness of gut microbiota, compared to HC. The alpha-diversity was lower in anti-citrullinated peptide antibodies (ACPA)-positive patients than in HC. The phylum Verrucomicrobiae and genus Akkermansia were more abundant in patients compared to HC. There was increased relative abundance of Enterobacteriaceae as well as Klebsiella, and less abundance of Bifidobacterium in patients with high levels of TNF-α or IL-17A compared to those who had low levels of these cytokines. In addition, ACPA-positive patients had higher proportions of Blautia, Akkermansia, and Clostridiales than ACPA-negative patients. Gut dysbiosis in RA patients was presented as different microbial composition and its association with inflammatory parameters as well as ACPA seropositivity. These findings support the involvement of gut microbiota in RA pathogenesis.

CMEP: a database for circulating microRNA expression profiling.

MOTIVATION: In recent years, several experimental studies have revealed that the microRNAs (miRNAs) in serum, plasma, exosome and whole blood are dysregulated in various types of diseases, indicating that the circulating miRNAs may serve as potential noninvasive biomarkers for disease diagnosis and prognosis. However, no database has been constructed to integrate the large-scale circulating miRNA profiles, explore the functional pathways involved and predict the potential biomarkers using feature selection between the disease conditions. Although there have been several studies attempting to generate a circulating miRNA database, they have not yet integrated the large-scale circulating miRNA profiles or provided the biomarker-selection function using machine learning methods. RESULTS: To fill this gap, we constructed the Circulating MicroRNA Expression Profiling (CMEP) database for integrating, analyzing and visualizing the large-scale expression profiles of phenotype-specific circulating miRNAs. The CMEP database contains massive datasets that were manually curated from NCBI GEO and the exRNA Atlas, including 66 datasets, 228 subsets and 10 419 samples. The CMEP provides the differential expression circulating miRNAs analysis and the KEGG functional pathway enrichment analysis. Furthermore, to provide the function of noninvasive biomarker discovery, we implemented several feature-selection methods, including ridge regression, lasso regression, support vector machine and random forests. Finally, we implemented a user-friendly web interface to improve the user experience and to visualize the data and results of CMEP. AVAILABILITY AND IMPLEMENTATION: CMEP is accessible at http://syslab5.nchu.edu.tw/CMEP.

Plant stress RNA-seq Nexus: a stress-specific transcriptome database in plant cells.

BACKGROUND: Abiotic and biotic stresses severely affect the growth and reproduction of plants and crops. Determining the critical molecular mechanisms and cellular processes in response to stresses will provide biological insight for addressing both climate change and food crises. RNA sequencing (RNA-Seq) is a revolutionary tool that has been used extensively in plant stress research. However, no existing large-scale RNA-Seq database has been designed to provide information on the stress-specific differentially expressed transcripts that occur across diverse plant species and various stresses. RESULTS: We have constructed a comprehensive database, the plant stress RNA-Seq nexus (PSRN), which includes 12 plant species, 26 plant-stress RNA-Seq datasets, and 937 samples. All samples are assigned to 133 stress-specific subsets, which are constructed into 254 subset pairs, a comparison between selected two subsets, for stress-specific differentially expressed transcript identification. CONCLUSIONS: PSRN is an open resource for intuitive data exploration, providing expression profiles of coding-transcript/lncRNA and identifying which transcripts are differentially expressed between different stress-specific subsets, in order to support researchers generating new biological insights and hypotheses in molecular breeding or evolution. PSRN is freely available at http://syslab5.nchu.edu.tw/PSRN .

CancerDetector: ultrasensitive and non-invasive cancer detection at the resolution of individual reads using cell-free DNA methylation sequencing data.

The detection of tumor-derived cell-free DNA in plasma is one of the most promising directions in cancer diagnosis. The major challenge in such an approach is how to identify the tiny amount of tumor DNAs out of total cell-free DNAs in blood. Here we propose an ultrasensitive cancer detection method, termed 'CancerDetector', using the DNA methylation profiles of cell-free DNAs. The key of our method is to probabilistically model the joint methylation states of multiple adjacent CpG sites on an individual sequencing read, in order to exploit the pervasive nature of DNA methylation for signal amplification. Therefore, CancerDetector can sensitively identify a trace amount of tumor cfDNAs in plasma, at the level of individual reads. We evaluated CancerDetector on the simulated data, and showed a high concordance of the predicted and true tumor fraction. Testing CancerDetector on real plasma data demonstrated its high sensitivity and specificity in detecting tumor cfDNAs. In addition, the predicted tumor fraction showed great consistency with tumor size and survival outcome. Note that all of those testing were performed on sequencing data at low to medium coverage (1× to 10×). Therefore, CancerDetector holds the great potential to detect cancer early and cost-effectively.

The therapeutic potential of CETP inhibitors: a patent review.

INTRODUCTION: Epidemiological studies have identified that high levels of low-density lipoprotein-cholesterol (LDL-C) and low levels of high-density lipoprotein-cholesterol (HDL-C) are two independent causes of cardiovascular disease (CVD). Statins, niacin and fibrate are used for the treatment of CVD. However, some defects are shown in the treatment process. Thus, there is a demand for better treatment strategies that confer preferable efficacy with fewer side effects. Cholesteryl ester transfer protein (CETP) promotes the movement of CEs from HDL to LDL and VLDL in exchange for triglycerides (TGs). AREAS COVERED: In this review, we reviewed the development and therapeutic applications of CETP inhibitors. A comprehensive review of the patents and pharmaceutical applications between 2009 and 2017 has been highlighted. EXPERT OPINION: Recently, CETP inhibitors have attracted considerable interest in atherosclerosis-related disease. There are four drugs (torcetrapib, anacetrapib, evacetrapib and dalcetrapib) that have been clinically evaluated in phase III clinical trials and showed promising results in raising HDL-C levels, but there were suboptimal performances in reducing the risk of cardiovascular events with all the compounds. The correlation between plasma HDL-C levels and CVD incidence needs further verification. The timeline is still long for CETP inhibitors to emerge from the treatment of CVD.

Biclustering of transcriptome sequencing data reveals human tissue-specific circular RNAs.

BACKGROUND: Emerging evidence has been experimentally confirmed the tissue-specific expression of circRNAs (circRNAs). Global identification of human tissue-specific circRNAs is crucial for the functionality study, which facilitates the discovery of circRNAs for potential diagnostic biomarkers. RESULTS: In this study, circRNA back-splicing junctions were identified from 465 publicly available transcriptome sequencing samples. The number of reads aligned to these identified junctions was normalized with the read length and sequence depth for each sample. We generated 66 models representing enriched circRNAs among human tissue transcriptome through biclustering algorithm. The result provides thousands of newly identified human tissue-specific circRNAs. CONCLUSIONS: This result suggests that expression of circRNAs is not prompted by random splicing error but serving molecular functional roles. We also identified circRNAs enriched within circulating system, which, along with identified tissue-specific circRNAs, can serve as potential diagnostic biomarkers.

Design, Synthesis, and Biological Evaluation of N,N-Disubstituted-4-Arylthiazole-2-Methylamine Derivatives as Cholesteryl Ester Transfer Inhibitors.

Cholesteryl ester transfer protein (CETP) has been identified as a potential target for cardiovascular disease (CVD) for its important role in the reverse cholesteryl transfer (RCT) process. In our previous work, compound 5 was discovered as a moderate CETP inhibitor. The replacement of the amide linker by heterocyclic aromatics and then a series of N,N-substituted-4-arylthiazole-2-methylamine derivatives were designed by utilizing a conformational restriction strategy. Thirty-six compounds were synthesized and evaluated for their CETP inhibitory activities. Structure-activity relationship studies indicate that electron donor groups substituted ring A, and electron-withdrawing groups at the 4-position of ring B were critical for potency. Among these compounds, compound 30 exhibited excellent CETP inhibitory activity (IC50 = 0.79 ± 0.02 µM) in vitro and showed an acceptable metabolic stability.

Design, Synthesis and Biological Evaluation of N,N-Substituted Amine Derivatives as Cholesteryl Ester Transfer Protein Inhibitors.

N,N-Substituted amine derivatives were designed by utilizing a bioisosterism strategy. Consequently, twenty-two compounds were synthesized and evaluated for their inhibitory activity against CETP. Structure-activity relationship (SAR) studies indicate that hydrophilic groups at the 2-position of the tetrazole and 3,5-bistrifluoromethyl groups on the benzene ring provide important contributions to the potency. Among these compounds, compound 17 exhibited excellent CETP inhibitory activity (IC50 = 0.38 ± 0.08 µM) in vitro. Furthermore, compound 17 was selected for an in vitro metabolic stability study.

Design, synthesis and evaluation of aromatic heterocyclic derivatives as potent antifungal agents.

To further enhance the anti-Aspergillus efficacy of our previously discovered antifungal lead compounds (1), a series of aromatic heterocyclic derivatives were designed, synthesized and evaluated for in vitro antifungal activity. Many of the target compounds showed good inhibitory activity against Candida albicans and Cryptococcus neoformans. In particular, the isoxazole nuclei were more suited for improving the activity against Aspergillus spp. Among these compounds, 2-F substituted analogues 23g and 23h displayed the most remarkable in vitro activity against Candida spp., C. neoformans, A. fumigatus and fluconazole-resistant C.alb. strains, which is superior or comparable to the activity of the reference drugs fluconazole and voriconazole. Notably, the compounds 23g and 23h exhibited low inhibition profiles for various isoforms of human cytochrome P450 and excellent blood plasma stability.

Improved n-butanol production via co-expression of membrane-targeted tilapia metallothionein and the clostridial metabolic pathway in Escherichia coli.

BACKGROUND: N-Butanol has favorable characteristics for use as either an alternative fuel or platform chemical. Bio-based n-butanol production using microbes is an emerging technology that requires further development. Although bio-industrial microbes such as Escherichia coli have been engineered to produce n-butanol, reactive oxygen species (ROS)-mediated toxicity may limit productivity. Previously, we show that outer-membrane-targeted tilapia metallothionein (OmpC-TMT) is more effective as an ROS scavenger than human and mouse metallothioneins to reduce oxidative stress in the host cell. RESULTS: The host strain (BUT1-DE) containing the clostridial n-butanol pathway displayed a decreased growth rate and limited n-butanol productivity, likely due to ROS accumulation. The clostridial n-butanol pathway was co-engineered with inducible OmpC-TMT in E. coli (BUT3-DE) for simultaneous ROS removal, and its effect on n-butanol productivity was examined. The ROS scavenging ability of cells overexpressing OmpC-TMT was examined and showed an approximately twofold increase in capacity. The modified strain improved n-butanol productivity to 320 mg/L, whereas the control strain produced only 95.1 mg/L. Transcriptomic analysis revealed three major KEGG pathways that were significantly differentially expressed in the BUT3-DE strain compared with their expression in the BUT1-DE strain, including genes involved in oxidative phosphorylation, fructose and mannose metabolism and glycolysis/gluconeogenesis. CONCLUSIONS: These results indicate that OmpC-TMT can increase n-butanol production by scavenging ROS. The transcriptomic analysis suggested that n-butanol causes quinone malfunction, resulting in oxidative-phosphorylation-related nuo operon downregulation, which would diminish the ability to convert NADH to NAD+ and generate proton motive force. However, fructose and mannose metabolism-related genes (fucA, srlE and srlA) were upregulated, and glycolysis/gluconeogenesis-related genes (pfkB, pgm) were downregulated, which further assisted in regulating NADH/NAD+ redox and preventing additional ATP depletion. These results indicated that more NADH and ATP were required in the n-butanol synthetic pathway. Our study demonstrates a potential approach to increase the robustness of microorganisms and the production of toxic chemicals through the ability to reduce oxidative stress.