BIC: a database for the transcriptional landscape of bacteria in cancer.

Microbial communities are massively resident in the human body, yet dysbiosis has been reported to correlate with many diseases, including various cancers. Most studies focus on the gut microbiome, while the bacteria that participate in tumor microenvironments on site remain unclear. Previous studies have acquired the bacteria expression profiles from RNA-seq, whole genome sequencing, and whole exon sequencing in The Cancer Genome Atlas (TCGA). However, small-RNA sequencing data were rarely used. Using TCGA miRNA sequencing data, we evaluated bacterial abundance in 32 types of cancer. To uncover the bacteria involved in cancer, we applied an analytical process to align unmapped human reads to bacterial references and developed the BIC database for the transcriptional landscape of bacteria in cancer. BIC provides cancer-associated bacterial information, including the relative abundance of bacteria, bacterial diversity, associations with clinical relevance, the co-expression network of bacteria and human genes, and their associated biological functions. These results can complement previously published databases. Users can easily download the result plots and tables, or download the bacterial abundance matrix for further analyses. In summary, BIC can provide information on cancer microenvironments related to microbial communities. BIC is available at: http://bic.jhlab.tw/.

Macrophages as determinants and regulators of systemic sclerosis-related interstitial lung disease.

BACKGROUND: Interstitial lung disease (ILD) is the primary cause of mortality in systemic sclerosis (SSc), an autoimmune disease characterized by tissue fibrosis. SSc-related ILD (SSc-ILD) occurs more frequently in females aged 30-55 years, whereas idiopathic pulmonary fibrosis (IPF) is more prevalent in males aged 60-75 years. SSc-ILD occurs earlier than IPF and progresses rapidly. FCN1, FABP4, and SPP1 macrophages are involved in the pathogenesis of lung fibrosis; SPP1 macrophages demonstrate upregulated expression in both SSc-ILD and IPF. To identify the differences between SSc-ILD and IPF using single-cell analysis, clarify their distinct pathogeneses, and propose directions for prevention and treatment. METHODS: We performed single-cell RNA sequencing on NCBI Gene Expression Omnibus (GEO) databases GSE159354 and GSE212109, and analyzed lung tissue samples across healthy controls, IPF, and SSc-ILD. The primary measures were the filtered genes integrated with batch correction and annotated cell types for distinguishing patients with SSc-ILD from healthy controls. We proposed an SSc-ILD pathogenesis using cell-cell interaction inferences, and predicted transcription factors regulating target genes using SCENIC. Drug target prediction of the TF gene was performed using Drug Bank Online. RESULTS: A subset of macrophages activates the MAPK signaling pathway under oxidative stress. Owing to the lack of inhibitory feedback from ANNEXIN and the autoimmune characteristics, this leads to an earlier onset of lung fibrosis compared to IPF. During initial lung injury, fibroblasts begin to activate the IL6 pathway under the influence of SPP1 alveolar macrophages, but IL6 appears unrelated to other inflammatory and immune cells. This may explain why tocilizumab (an anti-IL6-receptor antibody) only preserves lung function in patients with early SSc-ILD. Finally, we identified BCLAF1 and NFE2L2 as influencers of MAPK activation in macrophages. Metformin downregulates NFE2L2 and could serve as a repurposed drug candidate. CONCLUSIONS: SPP1 alveolar macrophages play a role in the profibrotic activity of IPF and SSc-ILD. However, SSc-ILD is influenced by autoimmunity and oxidative stress, leading to the continuous activation of MAPK in macrophages. This may result in an earlier onset of lung fibrosis than in IPF. Such differences could serve as potential research directions for early prevention and treatment.

GenEpi: gene-based epistasis discovery using machine learning.

BACKGROUND: Genome-wide association studies (GWAS) provide a powerful means to identify associations between genetic variants and phenotypes. However, GWAS techniques for detecting epistasis, the interactions between genetic variants associated with phenotypes, are still limited. We believe that developing an efficient and effective GWAS method to detect epistasis will be a key for discovering sophisticated pathogenesis, which is especially important for complex diseases such as Alzheimer's disease (AD). RESULTS: In this regard, this study presents GenEpi, a computational package to uncover epistasis associated with phenotypes by the proposed machine learning approach. GenEpi identifies both within-gene and cross-gene epistasis through a two-stage modeling workflow. In both stages, GenEpi adopts two-element combinatorial encoding when producing features and constructs the prediction models by L1-regularized regression with stability selection. The simulated data showed that GenEpi outperforms other widely-used methods on detecting the ground-truth epistasis. As real data is concerned, this study uses AD as an example to reveal the capability of GenEpi in finding disease-related variants and variant interactions that show both biological meanings and predictive power. CONCLUSIONS: The results on simulation data and AD demonstrated that GenEpi has the ability to detect the epistasis associated with phenotypes effectively and efficiently. The released package can be generalized to largely facilitate the studies of many complex diseases in the near future.

Emergence of differentially regulated pathways associated with the development of regional specificity in chicken skin.

BACKGROUND: Regional specificity allows different skin regions to exhibit different characteristics, enabling complementary functions to make effective use of the integumentary surface. Chickens exhibit a high degree of regional specificity in the skin and can serve as a good model for when and how these regional differences begin to emerge. RESULTS: We used developing feather and scale regions in embryonic chickens as a model to gauge the differences in their molecular pathways. We employed cosine similarity analysis to identify the differentially regulated and co-regulated genes. We applied low cell techniques for expression validation and chromatin immunoprecipitation (ChIP)-based enhancer identification to overcome limited cell availabilities from embryonic chicken skin. We identified a specific set of genes demonstrating a high correlation as being differentially expressed during feather and scale development and maturation. Some members of the WNT, TGF-beta/BMP, and Notch family known to be involved in feathering skin differentiation were found to be differentially regulated. Interestingly, we also found genes along calcium channel pathways that are differentially regulated. From the analysis of differentially regulated pathways, we used calcium signaling pathways as an example for further verification. Some voltage-gated calcium channel subunits, particularly CACNA1D, are expressed spatio-temporally in the skin epithelium. These calcium signaling pathway members may be involved in developmental decisions, morphogenesis, or epithelial maturation. We further characterized enhancers associated with histone modifications, including H3K4me1, H3K27ac, and H3K27me3, near calcium channel-related genes and identified signature intensive hotspots that may be correlated with certain voltage-gated calcium channel genes. CONCLUSION: We demonstrated the applicability of cosine similarity analysis for identifying novel regulatory pathways that are differentially regulated during development. Our study concerning the effects of signaling pathways and histone signatures on enhancers suggests that voltage-gated calcium signaling may be involved in early skin development. This work lays the foundation for studying the roles of these gene pathways and their genomic regulation during the establishment of skin regional specificity.

Risk of dementia after anaesthesia and surgery.

BACKGROUND: The potential relationship between anaesthesia, surgery and onset of dementia remains elusive. AIMS: To determine whether the risk of dementia increases after surgery with anaesthesia, and to evaluate possible associations among age, mode of anaesthesia, type of surgery and risk of dementia. METHOD: The study cohort comprised patients aged 50 years and older who were anaesthetised for the first time since 1995 between 1 January 2004 and 31 December 2007, and a control group of randomly selected patients matched for age and gender. Patients were followed until 31 December 2010 to identify the emergence of dementia. RESULTS: Relative to the control group, patients who underwent anaesthesia and surgery exhibited an increased risk of dementia (hazard ratio = 1.99) and a reduced mean interval to dementia diagnosis. The risk of dementia increased in patients who received intravenous or intramuscular anaesthesia, regional anaesthesia and general anaesthesia. CONCLUSIONS: The results of our nationwide, population-based study suggest that patients who undergo anaesthesia and surgery may be at increased risk of dementia.

scDrug+: predicting drug-responses using single-cell transcriptomics and molecular structure.

Predicting drug responses based on individual transcriptomic profiles holds promise for refining prognosis and advancing precision medicine. Although many studies have endeavored to predict the responses of known drugs to novel transcriptomic profiles, research into predicting responses for newly discovered drugs remains sparse. In this study, we introduce scDrug+, a comprehensive pipeline that seamlessly integrates single-cell analysis with drug-response prediction. Importantly, scDrug+ is equipped to predict the response of new drugs by analyzing their molecular structures. The open-source tool is available as a Docker container, ensuring ease of deployment and reproducibility. It can be accessed at https://github.com/ailabstw/scDrugplus.

De novo motif discovery facilitates identification of interactions between transcription factors in Saccharomyces cerevisiae.

MOTIVATION: Gene regulation involves complicated mechanisms such as cooperativity between a set of transcription factors (TFs). Previous studies have used target genes shared by two TFs as a clue to infer TF-TF interactions. However, this task remains challenging because the target genes with low binding affinity are frequently omitted by experimental data, especially when a single strict threshold is employed. This article aims at improving the accuracy of inferring TF-TF interactions by incorporating motif discovery as a fundamental step when detecting overlapping targets of TFs based on ChIP-chip data. RESULTS: The proposed method, simTFBS, outperforms three naïve methods that adopt fixed thresholds when inferring TF-TF interactions based on ChIP-chip data. In addition, simTFBS is compared with two advanced methods and demonstrates its advantages in predicting TF-TF interactions. By comparing simTFBS with predictions based on the set of available annotated yeast TF binding motifs, we demonstrate that the good performance of simTFBS is indeed coming from the additional motifs found by the proposed procedures. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

MicroRNA-regulated protein-protein interaction networks and their functions in breast cancer.

MicroRNAs, which are small endogenous RNA regulators, have been associated with various types of cancer. Breast cancer is a major health threat for women worldwide. Many miRNAs were reported to be associated with the progression and carcinogenesis of breast cancer. In this study, we aimed to discover novel breast cancer-related miRNAs and to elucidate their functions. First, we identified confident miRNA-target pairs by combining data from miRNA target prediction databases and expression profiles of miRNA and mRNA. Then, miRNA-regulated protein interaction networks (PINs) were constructed with confident pairs and known interaction data in the human protein reference database (HPRD). Finally, the functions of miRNA-regulated PINs were elucidated by functional enrichment analysis. From the results, we identified some previously reported breast cancer-related miRNAs and functions of the PINs, e.g., miR-125b, miR-125a, miR-21, and miR-497. Some novel miRNAs without known association to breast cancer were also found, and the putative functions of their PINs were also elucidated. These include miR-139 and miR-383. Furthermore, we validated our results by receiver operating characteristic (ROC) curve analysis using our miRNA expression profile data, gene expression-based outcome for breast cancer online (GOBO) survival analysis, and a literature search. Our results may provide new insights for research in breast cancer-associated miRNAs.

Exploiting Machine Learning Technologies to Study the Compound Effects of Serum Creatinine and Electrolytes on the Risk of Acute Kidney Injury in Intensive Care Units.

Assessing the risk of acute kidney injury (AKI) has been a challenging issue for clinicians in intensive care units (ICUs). In recent years, a number of studies have been conducted to investigate the associations between several serum electrolytes and AKI. Nevertheless, the compound effects of serum creatinine, blood urea nitrogen (BUN), and clinically relevant serum electrolytes have yet to be comprehensively investigated. Accordingly, we initiated this study aiming to develop machine learning models that illustrate how these factors interact with each other. In particular, we focused on ICU patients without a prior history of AKI or AKI-related comorbidities. With this practice, we were able to examine the associations between the levels of serum electrolytes and renal function in a more controlled manner. Our analyses revealed that the levels of serum creatinine, chloride, and magnesium were the three major factors to be monitored for this group of patients. In summary, our results can provide valuable insights for developing early intervention and effective management strategies as well as crucial clues for future investigations of the pathophysiological mechanisms that are involved. In future studies, subgroup analyses based on different causes of AKI should be conducted to further enhance our understanding of AKI.

A Single-Cell Network-Based Drug Repositioning Strategy for Post-COVID-19 Pulmonary Fibrosis.

Post-COVID-19 pulmonary fibrosis (PCPF) is a long-term complication that appears in some COVID-19 survivors. However, there are currently limited options for treating PCPF patients. To address this problem, we investigated COVID-19 patients' transcriptome at single-cell resolution and combined biological network analyses to repurpose the drugs treating PCPF. We revealed a novel gene signature of PCPF. The signature is functionally associated with the viral infection and lung fibrosis. Further, the signature has good performance in diagnosing and assessing pulmonary fibrosis. Next, we applied a network-based drug repurposing method to explore novel treatments for PCPF. By quantifying the proximity between the drug targets and the signature in the interactome, we identified several potential candidates and provided a drug list ranked by their proximity. Taken together, we revealed a novel gene expression signature as a theragnostic biomarker for PCPF by integrating different computational approaches. Moreover, we showed that network-based proximity could be used as a framework to repurpose drugs for PCPF.

Stratification of lncRNA modulation networks in breast cancer.

BACKGROUND: Recently, non-coding RNAs are of growing interest, and more scientists attach importance to research on their functions. Long non-coding RNAs (lncRNAs) are defined as non-protein coding transcripts longer than 200 nucleotides. We already knew that lncRNAs are related to cancers and will be dysregulated in them. But most of their functions are still left to further study. A mechanism of RNA regulation, known as competing endogenous RNAs (ceRNAs), has been proposed to explain the complex relationships among mRNAs and lncRNAs by competing for binding with shared microRNAs (miRNAs). METHODS: We proposed an analysis framework to construct the association networks among lncRNA, mRNA, and miRNAs based on their expression patterns and decipher their network modules. RESULTS: We collected a large-scale gene expression dataset of 1,061 samples from breast invasive carcinoma (BRCA) patients, each consisted of the expression profiles of 4,359 lncRNAs, 16,517 mRNAs, and 534 miRNAs, and applied the proposed analysis approach to interrogate them. We have uncovered the underlying ceRNA modules and the key modulatory lncRNAs for different subtypes of breast cancer. CONCLUSIONS: We proposed a modulatory analysis to infer the ceRNA effects among mRNAs and lncRNAs and performed functional analysis to reveal the plausible mechanisms of lncRNA modulation in the four breast cancer subtypes. Our results might provide new directions for breast cancer therapeutics and the proposed method could be readily applied to other diseases.

Machine learning for emerging infectious disease field responses.

Emerging infectious diseases (EIDs), including the latest COVID-19 pandemic, have emerged and raised global public health crises in recent decades. Without existing protective immunity, an EID may spread rapidly and cause mass casualties in a very short time. Therefore, it is imperative to identify cases with risk of disease progression for the optimized allocation of medical resources in case medical facilities are overwhelmed with a flood of patients. This study has aimed to cope with this challenge from the aspect of preventive medicine by exploiting machine learning technologies. The study has been based on 83,227 hospital admissions with influenza-like illness and we analysed the risk effects of 19 comorbidities along with age and gender for severe illness or mortality risk. The experimental results revealed that the decision rules derived from the machine learning based prediction models can provide valuable guidelines for the healthcare policy makers to develop an effective vaccination strategy. Furthermore, in case the healthcare facilities are overwhelmed by patients with EID, which frequently occurred in the recent COVID-19 pandemic, the frontline physicians can incorporate the proposed prediction models to triage patients suffering minor symptoms without laboratory tests, which may become scarce during an EID disaster. In conclusion, our study has demonstrated an effective approach to exploit machine learning technologies to cope with the challenges faced during the outbreak of an EID.

ProteDNA: a sequence-based predictor of sequence-specific DNA-binding residues in transcription factors.

This article presents the design of a sequence-based predictor named ProteDNA for identifying the sequence-specific binding residues in a transcription factor (TF). Concerning protein-DNA interactions, there are two types of binding mechanisms involved, namely sequence-specific binding and nonspecific binding. Sequence-specific bindings occur between protein sidechains and nucleotide bases and correspond to sequence-specific recognition of genes. Therefore, sequence-specific bindings are essential for correct gene regulation. In this respect, ProteDNA is distinctive since it has been designed to identify sequence-specific binding residues. In order to accommodate users with different application needs, ProteDNA has been designed to operate under two modes, namely, the high-precision mode and the balanced mode. According to the experiments reported in this article, under the high-precision mode, ProteDNA has been able to deliver precision of 82.3%, specificity of 99.3%, sensitivity of 49.8% and accuracy of 96.5%. Meanwhile, under the balanced mode, ProteDNA has been able to deliver precision of 60.8%, specificity of 97.6%, sensitivity of 60.7% and accuracy of 95.4%. ProteDNA is available at the following websites: http://protedna.csbb.ntu.edu.tw/, http://protedna.csie.ntu.edu.tw/, http://bio222.esoe.ntu.edu.tw/ProteDNA/.

Predicting microRNA precursors with a generalized Gaussian components based density estimation algorithm.

BACKGROUND: MicroRNAs (miRNAs) are short non-coding RNA molecules, which play an important role in post-transcriptional regulation of gene expression. There have been many efforts to discover miRNA precursors (pre-miRNAs) over the years. Recently, ab initio approaches have attracted more attention because they do not depend on homology information and provide broader applications than comparative approaches. Kernel based classifiers such as support vector machine (SVM) are extensively adopted in these ab initio approaches due to the prediction performance they achieved. On the other hand, logic based classifiers such as decision tree, of which the constructed model is interpretable, have attracted less attention. RESULTS: This article reports the design of a predictor of pre-miRNAs with a novel kernel based classifier named the generalized Gaussian density estimator (G2DE) based classifier. The G2DE is a kernel based algorithm designed to provide interpretability by utilizing a few but representative kernels for constructing the classification model. The performance of the proposed predictor has been evaluated with 692 human pre-miRNAs and has been compared with two kernel based and two logic based classifiers. The experimental results show that the proposed predictor is capable of achieving prediction performance comparable to those delivered by the prevailing kernel based classification algorithms, while providing the user with an overall picture of the distribution of the data set. CONCLUSION: Software predictors that identify pre-miRNAs in genomic sequences have been exploited by biologists to facilitate molecular biology research in recent years. The G2DE employed in this study can deliver prediction accuracy comparable with the state-of-the-art kernel based machine learning algorithms. Furthermore, biologists can obtain valuable insights about the different characteristics of the sequences of pre-miRNAs with the models generated by the G2DE based predictor.

Risk assessment of dementia after hysterectomy: Analysis of 14-year data from the National Health Insurance Research Database in Taiwan.

BACKGROUND: Anesthesia and surgery may increase the risk of dementia in the elderly, but the higher prevalence of dementia in women and other evidence suggest that dementia risk increases in younger women undergoing hysterectomy. In this study, we assessed the risk of dementia after hysterectomy. METHODS: Hysterectomies registered in the National Health Insurance Research Database from 2000 to 2013 were evaluated using a retrospective generational research method. Multivariate Cox regression analysis was used to assess the effect of age at surgery, anesthesia method, and surgery type on the hazard ratio (HR) for the development of dementia. RESULTS: Among 280 308 patients who underwent hysterectomy, 4753 (1.7%) developed dementia. Age at surgery and anesthesia method were associated with the occurrence of dementia, independent of surgery type. Among patients 30-49 years of age, general anesthesia (GA) was associated with a higher risk of dementia than spinal anesthesia (SA). The HR for GA was 2.678 (95% confidence interval [CI] = 1.269-5.650) and the risk of dementia increased by 7.4% for every 1-year increase in age (HR = 1.074; 95% CI = 1.048-1.101). In patients >50 years of age, the HR for GA was 1.206 (95% CI = 1.057-1.376), and the risk of dementia increased by 13.0% for every 1-year increase in age (HR = 1.130; 95% CI = 1.126-1.134). CONCLUSION: The risk of dementia in women who underwent hysterectomy was significantly affected by older age at surgery, and the risk might not increase linearly with age, but show instead an S-curve with exponential increase at about 50 years of age. Although less significant, GA was associated with higher risk than SA, and the effect of the anesthesia method was greater in patients <50 years of age. In contrast, the surgical procedure used was not associated to the risk of dementia.

Comparing machine learning with case-control models to identify confirmed dengue cases.

In recent decades, the global incidence of dengue has increased. Affected countries have responded with more effective surveillance strategies to detect outbreaks early, monitor the trends, and implement prevention and control measures. We have applied newly developed machine learning approaches to identify laboratory-confirmed dengue cases from 4,894 emergency department patients with dengue-like illness (DLI) who received laboratory tests. Among them, 60.11% (2942 cases) were confirmed to have dengue. Using just four input variables [age, body temperature, white blood cells counts (WBCs) and platelets], not only the state-of-the-art deep neural network (DNN) prediction models but also the conventional decision tree (DT) and logistic regression (LR) models delivered performances with receiver operating characteristic (ROC) curves areas under curves (AUCs) of the ranging from 83.75% to 85.87% [for DT, DNN and LR: 84.60% ± 0.03%, 85.87% ± 0.54%, 83.75% ± 0.17%, respectively]. Subgroup analyses found all the models were very sensitive particularly in the pre-epidemic period. Pre-peak sensitivities (<35 weeks) were 92.6%, 92.9%, and 93.1% in DT, DNN, and LR respectively. Adjusted odds ratios examined with LR for low WBCs [≤ 3.2 (x103/µL)], fever (≥38°C), low platelet counts [< 100 (x103/µL)], and elderly (≥ 65 years) were 5.17 [95% confidence interval (CI): 3.96-6.76], 3.17 [95%CI: 2.74-3.66], 3.10 [95%CI: 2.44-3.94], and 1.77 [95%CI: 1.50-2.10], respectively. Our prediction models can readily be used in resource-poor countries where viral/serologic tests are inconvenient and can also be applied for real-time syndromic surveillance to monitor trends of dengue cases and even be integrated with mosquito/environment surveillance for early warning and immediate prevention/control measures. In other words, a local community hospital/clinic with an instrument of complete blood counts (including platelets) can provide a sentinel screening during outbreaks. In conclusion, the machine learning approach can facilitate medical and public health efforts to minimize the health threat of dengue epidemics. However, laboratory confirmation remains the primary goal of surveillance and outbreak investigation.

A sequence-based hybrid predictor for identifying conformationally ambivalent regions in proteins.

BACKGROUND: Proteins are dynamic macromolecules which may undergo conformational transitions upon changes in environment. As it has been observed in laboratories that protein flexibility is correlated to essential biological functions, scientists have been designing various types of predictors for identifying structurally flexible regions in proteins. In this respect, there are two major categories of predictors. One category of predictors attempts to identify conformationally flexible regions through analysis of protein tertiary structures. Another category of predictors works completely based on analysis of the polypeptide sequences. As the availability of protein tertiary structures is generally limited, the design of predictors that work completely based on sequence information is crucial for advances of molecular biology research. RESULTS: In this article, we propose a novel approach to design a sequence-based predictor for identifying conformationally ambivalent regions in proteins. The novelty in the design stems from incorporating two classifiers based on two distinctive supervised learning algorithms that provide complementary prediction powers. Experimental results show that the overall performance delivered by the hybrid predictor proposed in this article is superior to the performance delivered by the existing predictors. Furthermore, the case study presented in this article demonstrates that the proposed hybrid predictor is capable of providing the biologists with valuable clues about the functional sites in a protein chain. The proposed hybrid predictor provides the users with two optional modes, namely, the high-sensitivity mode and the high-specificity mode. The experimental results with an independent testing data set show that the proposed hybrid predictor is capable of delivering sensitivity of 0.710 and specificity of 0.608 under the high-sensitivity mode, while delivering sensitivity of 0.451 and specificity of 0.787 under the high-specificity mode. CONCLUSION: Though experimental results show that the hybrid approach designed to exploit the complementary prediction powers of distinctive supervised learning algorithms works more effectively than conventional approaches, there exists a large room for further improvement with respect to the achieved performance. In this respect, it is of interest to investigate the effects of exploiting additional physiochemical properties that are related to conformational ambivalence. Furthermore, it is of interest to investigate the effects of incorporating lately-developed machine learning approaches, e.g. the random forest design and the multi-stage design. As conformational transition plays a key role in carrying out several essential types of biological functions, the design of more advanced predictors for identifying conformationally ambivalent regions in proteins deserves our continuous attention.

DNA-binding residues and binding mode prediction with binding-mechanism concerned models.

BACKGROUND: Protein-DNA interactions are essential for fundamental biological activities including DNA transcription, replication, packaging, repair and rearrangement. Proteins interacting with DNA can be classified into two categories of binding mechanisms - sequence-specific and non-specific binding. Protein-DNA specific binding provides a mechanism to recognize correct nucleotide base pairs for sequence-specific identification. Protein-DNA non-specific binding shows sequence independent interaction for accelerated targeting by interacting with DNA backbone. Both sequence-specific and non-specific binding residues contribute to their roles for interaction. RESULTS: The proposed framework has two stage predictors: DNA-binding residues prediction and binding mode prediction. In the first stage - DNA-binding residues prediction, the predictor for DNA specific binding residues achieves 96.45% accuracy with 50.14% sensitivity, 99.31% specificity, 81.70% precision, and 62.15% F-measure. The predictor for DNA non-specific binding residues achieves 89.14% accuracy with 53.06% sensitivity, 95.25% specificity, 65.47% precision, and 58.62% F-measure. While combining prediction results of sequence-specific and non-specific binding residues with OR operation, the predictor achieves 89.26% accuracy with 56.86% sensitivity, 95.63% specificity, 71.92% precision, and 63.51% F-measure. In the second stage, protein-DNA binding mode prediction achieves 75.83% accuracy while using support vector machine with multi-class prediction. CONCLUSION: This article presents the design of a sequence based predictor aiming to identify sequence-specific and non-specific binding residues in a transcription factor with DNA binding-mechanism concerned. The protein-DNA binding mode prediction was introduced to help improve DNA-binding residues prediction. In addition, the results of this study will help with the design of binding-mechanism concerned predictors for other families of proteins interacting with DNA.

Effect of de novo transcriptome assembly on transcript quantification.

Correct quantification of transcript expression is essential to understand the functional elements in different physiological conditions. For the organisms without the reference transcriptome, de novo transcriptome assembly must be carried out prior to quantification. However, a large number of erroneous contigs produced by the assemblers might result in unreliable estimation. In this regard, this study investigates how assembly quality affects the performance of quantification based on de novo transcriptome assembly. We examined the over-extended and incomplete contigs, and demonstrated that assembly completeness has a strong impact on the estimation of contig abundance. Then we investigated the behavior of the quantifiers with respect to sequence ambiguity which might be originally presented in the transcriptome or accidentally produced by assemblers. The results suggested that the quantifiers often over-estimate the expression of family-collapse contigs and under-estimate the expression of duplicated contigs. For organisms without reference transcriptome, it remains challenging to detect the inaccurate estimation on family-collapse contigs. On the contrary, we observed that the situation of under-estimation on duplicated contigs can be warned through analyzing the read proportion of estimated abundance (RPEA) of contigs in the connected component inferenced by the quantifiers. In addition, we suggest that the estimated quantification results on the connected component level have better accuracy over sequence level quantification. The analytic results conducted in this study provides valuable insights for future development of transcriptome assembly and quantification.

Author Correction: Effect of de novo transcriptome assembly on transcript quantification.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.