Search

DEAD/H-box helicases are implicated in virtually every aspect of RNA metabolism, including transcription, pre-mRNA splicing, ribosomes biogenesis, nuclear export, translation initiation, RNA degradation, and mRNA editing. Most of these helicases are upregulated in various cancers and mutations in some of them are associated with several malignancies. Lately, synthetic lethality (SL) and synthetic dosage lethality (SDL) approaches, where genetic interactions of cancer-related genes are exploited as therapeutic targets, are emerging as a leading area of cancer research. Several DEAD/H-box helicases, including DDX3, DDX9 (Dbp9), DDX10 (Dbp4), DDX11 (ChlR1), and DDX41 (Sacy-1), have been subjected to SL analyses in humans and different model organisms. It remains to be explored whether SDL can be utilized to identity druggable targets in DEAD/H-box helicase overexpressing cancers. In this review, we analyze gene expression data of a subset of DEAD/H-box helicases in multiple cancer types and discuss how their SL/SDL interactions can be used for therapeutic purposes. We also summarize the latest developments in clinical applications, apart from discussing some of the challenges in drug discovery in the context of targeting DEAD/H-box helicases.

Synthetic antibodies have been engineered against a wide variety of antigens with desirable biophysical, biochemical, and pharmacological properties. Here, we describe the generation and characterization of synthetic antigen-binding fragments (Fabs) against Notch-1. Three single-framework synthetic Fab libraries, named S, F, and modified-F, were screened against the recombinant human Notch-1 extracellular domain using phage display. These libraries were built on a modified trastuzumab framework, containing two or four diversified complementarity-determining regions (CDRs) and different CDR diversity designs. In total, 12 Notch-1 Fabs were generated with 10 different CDRH3 lengths. These Fabs possessed a high affinity for Notch-1 (sub-nM to mid-nM KDapp values) and exhibited different binding profiles (mono-, bi-or tri-specific) toward Notch/Jagged receptors. Importantly, we showed that screening focused diversity libraries, implementing next-generation sequencing approaches, and fine-tuning the CDR length diversity provided improved binding solutions for Notch-1 recognition. These findings have implications for antibody library design and antibody phage display.

Individual variability in responses to vaccination can result in vaccinated subjects failing to develop a protective immune response. Vaccine non-responders can remain susceptible to infection and may compromise efforts to achieve herd immunity. Biomarkers of vaccine unresponsiveness could aid vaccine research and development as well as strategically improve vaccine administration programs. We previously vaccinated piglets (n = 117) against a commercial Mycoplasma hyopneumoniae vaccine (RespiSure-One) and observed in low vaccine responder piglets, as defined by serum IgG antibody titers, differential phosphorylation of peptides involved in pro-inflammatory cytokine signaling within peripheral blood mononuclear cells (PBMCs) prior to vaccination, elevated plasma interferon-gamma concentrations, and lower birth weight compared to high vaccine responder piglets. In the current study, we use kinome analysis to investigate signaling events within PBMCs collected from the same high and low vaccine responders at 2 and 6 days post-vaccination. Furthermore, we evaluate the use of inflammatory plasma cytokines, birthweight, and signaling events as biomarkers of vaccine unresponsiveness in a validation cohort of high and low vaccine responders. Differential phosphorylation events (FDR 0.5) and plasma IFNγ concentrations at the time of vaccination (p > 0.6) between high and low responders within the validation cohort. The results in this study suggest, at least within this study population, phosphorylation biomarkers are more robust predictors of vaccine responsiveness than other physiological markers.

Computational drug repositioning aims to identify potential applications of existing drugs for the treatment of diseases for which they were not designed. This approach can considerably accelerate the traditional drug discovery process by decreasing the required time and costs of drug development. Tensor decomposition enables us to integrate multiple drug- and disease-related data to boost the performance of prediction. In this study, a nonnegative tensor decomposition for drug repositioning, NTD-DR, is proposed. In order to capture the hidden information in drug-target, drug-disease, and target-disease networks, NTD-DR uses these pairwise associations to construct a three-dimensional tensor representing drug-target-disease triplet associations and integrates them with similarity information of drugs, targets, and disease to make a prediction. We compare NTD-DR with recent state-of-the-art methods in terms of the area under the receiver operating characteristic (ROC) curve (AUC) and the area under the precision and recall curve (AUPR) and find that our method outperforms competing methods. Moreover, case studies with five diseases also confirm the reliability of predictions made by NTD-DR. Our proposed method identifies more known associations among the top 50 predictions than other methods. In addition, novel associations identified by NTD-DR are validated by literature analyses.

Abstract Background Antimicrobial resistance (AMR) is a major threat to global public health because it makes standard treatments ineffective and contributes to the spread of infections. It is important to understand AMR’s biological mechanisms for the development of new drugs and more rapid and accurate clinical diagnostics. The increasing availability of whole-genome SNP (single nucleotide polymorphism) information, obtained from whole-genome sequence data, along with AMR profiles provides an opportunity to use feature selection in machine learning to find AMR-associated mutations. This work describes the use of a supervised feature selection approach using deep neural networks to detect AMR-associated genetic factors from whole-genome SNP data. Results The proposed method, DNP-AAP (deep neural pursuit – average activation potential), was tested on a Neisseria gonorrhoeae dataset with paired whole-genome sequence data and resistance profiles to five commonly used antibiotics including penicillin, tetracycline, azithromycin, ciprofloxacin, and cefixime. The results show that DNP-AAP can effectively identify known AMR-associated genes in N. gonorrhoeae, and also provide a list of candidate genomic features (SNPs) that might lead to the discovery of novel AMR determinants. Logistic regression classifiers were built with the identified SNPs and the prediction AUCs (area under the curve) for penicillin, tetracycline, azithromycin, ciprofloxacin, and cefixime were 0.974, 0.969, 0.949, 0.994, and 0.976, respectively. Conclusions DNP-AAP can effectively identify known AMR-associated genes in N. gonorrhoeae. It also provides a list of candidate genes and intergenic regions that might lead to novel AMR factor discovery. More generally, DNP-AAP can be applied to AMR analysis of any bacterial species with genomic variants and phenotype data. It can serve as a useful screening tool for microbiologists to generate genetic candidates for further lab experiments.

Abstract Background Synthetic lethal interactions (SLIs) that occur between gene pairs are exploited for cancer therapeutics. Studies in the model eukaryote yeast have identified ~ 550,000 negative genetic interactions that have been extensively studied, leading to characterization of novel pathways and gene functions. This resource can be used to predict SLIs that can be relevant to cancer therapeutics. Methods We used patient data to identify genes that are down-regulated in breast cancer. InParanoid orthology mapping was performed to identify yeast orthologs of the down-regulated genes and predict their corresponding SLIs in humans. The predicted network graphs were drawn with Cytoscape. CancerRXgene database was used to predict drug response. Results Harnessing the vast available knowledge of yeast genetics, we generated a Humanized Yeast Genetic Interaction Network (HYGIN) for 1009 human genes with 10,419 interactions. Through the addition of patient-data from The Cancer Genome Atlas (TCGA), we generated a breast cancer specific subnetwork. Specifically, by comparing 1009 genes in HYGIN to genes that were down-regulated in breast cancer, we identified 15 breast cancer genes with 130 potential SLIs. Interestingly, 32 of the 130 predicted SLIs occurred with FBXW7, a well-known tumor suppressor that functions as a substrate-recognition protein within a SKP/CUL1/F-Box ubiquitin ligase complex for proteasome degradation. Efforts to validate these SLIs using chemical genetic data predicted that patients with loss of FBXW7 may respond to treatment with drugs like Selumitinib or Cabozantinib. Conclusions This study provides a patient-data driven interpretation of yeast SLI data. HYGIN represents a novel strategy to uncover therapeutically relevant cancer drug targets and the yeast SLI data offers a major opportunity to mine these interactions.

Abstract Background Recent adva1nces in whole genome sequencing (WGS) based technologies have facilitated multi-step applications for predicting antimicrobial resistance (AMR) and investigating the molecular epidemiology of Neisseria gonorrhoeae. However, generating full scaffolds of N. gonorrhoeae genomes from short reads, and the assignment of molecular epidemiological information (NG-MLST, NG-MAST, and NG-STAR) to multiple assembled samples, is challenging due to required manual tasks such as annotating antimicrobial resistance determinants with standard nomenclature for a large number of genomes. Results We present Gen2Epi, a pipeline that assembles short reads into full scaffolds and automatically assigns molecular epidemiological and AMR information to the assembled genomes. Gen2Epi is a command-line tool integrating third-party software and tailored specifically for N. gonorrhoeae. For its evaluation, the Gen2Epi pipeline successfully assembled the WGS short reads from 1484 N. gonorrhoeae samples into full-length genomes for both chromosomes and plasmids and was able to assign in silico molecular determinant information to each dataset automatically. The assemblies were generated using raw as well as trimmed short reads. The median genome coverage of full-length scaffolds and “N” statistics (N50, NG50, and NGA50) were higher than, or comparable to, previously published results and the scaffolding process improved the quality of the draft genome assemblies. Molecular antimicrobial resistant (AMR) determinants identified by Gen2Epi reproduced information for the 1484 samples as previously reported, including NG-MLST, NG-MAST, and NG-STAR molecular sequence types. Conclusions Gen2Epi can be used to assemble short reads into full-length genomes and assign accurate molecular marker and AMR information automatically from NG-STAR, NG-MAST, and NG-MLST. Gen2Epi is publicly available under “CC BY-NC 2.0 CA” Creative Commons licensing as a VirtualBox image containing the constituent software components running on the LINUX operating system (CentOS 7). The image and associated documentation are available via anonymous FTP at ftp://www.cs.usask.ca/pub/combi or ftp://ftp.cs.usask.ca/pub/combi

Peptide microarrays consisting of defined phosphorylation target sites are an effective approach for high throughput analysis of cellular kinase (kinome) activity. Kinome peptide arrays are highly customizable and do not require species-specific reagents to measure kinase activity, making them amenable for kinome analysis in any species. Our group developed software, Platform for Integrated, Intelligent Kinome Analysis (PIIKA), to enable more effective extraction of meaningful biological information from kinome peptide array data. A subsequent version, PIIKA2, unveiled new statistical tools and data visualization options. Here we introduce PIIKA 2.5 to provide two essential quality control metrics and a new background correction technique to increase the accuracy and consistency of kinome results. The first metric alerts users to improper spot size and informs them of the need to perform manual resizing to enhance the quality of the raw intensity data. The second metric uses inter-array comparisons to identify outlier arrays that sometimes emerge as a consequence of technical issues. In addition, a new background correction method, background scaling, can sharply reduce spatial biases within a single array in comparison to background subtraction alone. Collectively, the modifications of PIIKA 2.5 enable identification and correction of technical issues inherent to the technology and better facilitate the extraction of meaningful biological information. We show that these metrics demonstrably enhance kinome analysis by identifying low quality data and reducing batch effects, and ultimately improve clustering of treatment groups and enhance reproducibility. The web-based and stand-alone versions of PIIKA 2.5 are freely accessible at via http://saphire.usask.ca.

Within human health research, the remarkable utility of kinase inhibitors as therapeutics has motivated efforts to understand biology at the level of global cellular kinase activity (the kinome). In contrast, the diminished potential for using kinase inhibitors in food animals has dampened efforts to translate this research approach to livestock species. This, in our opinion, was a lost opportunity for livestock researchers given the unique potential of kinome analysis to offer insight into complex biology. To remedy this situation, our lab developed user-friendly, cost-effective approaches for kinome analysis that can be readily incorporated into most research programs but with a specific priority to enable the technology to livestock researchers. These contributions include the development of custom software programs for the creation of species-specific kinome arrays as well as comprehensive deconvolution and analysis of kinome array data. Presented in this review are examples of the application of kinome analysis to highlight the utility of the technology to further our understanding of two key complex biological events of priority to the livestock industry: host immune responses to infectious diseases and animal stress responses. These advances and examples of application aim to provide both mechanisms and motivation for researchers, particularly livestock researchers, to incorporate kinome analysis into their research programs.

Abstract Objectives Hi-C is a proximity-based ligation reaction used to detect regions of the genome that are close in 3D space (or “interacting”). Typically, results from Hi-C experiments (contact maps) are visualized as heatmaps or Circos plots. While informative, these visualizations do not directly represent genomic structure and folding, making the interpretation of the underlying 3D genomic organization obscured. Our objective was to generate a graph-based contact map representation that leads to a more intuitive structural visualization. Results Normalized contact maps were converted into undirected graphs where each vertex represented a genomic region and each edge represented a detected (intra- and inter-chromosomal) or known (linear) interaction between two regions. Each edge was weighted by the inverse of the linear distance (Hi-C experimental resolution) or the interaction frequency from the contact map. Graphs were generated based on this representation scheme for contact maps from existing fission yeast datasets. Originally, these datasets were used to (1) identify specific principles influencing fission yeast genome organization and (2) uncover changes in fission yeast genome organization during the cell cycle. When compared to the equivalent heatmaps and/or Circos plots, the graph-based visualizations more intuitively depicted the changes in genome organization described in the original studies.

Microbiome data consists of operational taxonomic unit (OTU) counts characterized by zero-inflation, over-dispersion, and grouping structure among samples. Currently, statistical testing methods are commonly performed to identify OTUs that are associated with a phenotype. The limitations of statistical testing methods include that the validity of p-values/q-values depend sensitively on the correctness of models and that the statistical significance does not necessarily imply predictivity. Predictive analysis using methods such as LASSO is an alternative approach for identifying associated OTUs and for measuring the predictability of the phenotype variable with OTUs and other covariate variables. We investigate three strategies of performing predictive analysis: (1) LASSO: fitting a LASSO multinomial logistic regression model to all OTU counts with specific transformation; (2) screening+GLM: screening OTUs with q-values returned by fitting a GLMM to each OTU, then fitting a GLM model using a subset of selected OTUs; (3) screening+LASSO: fitting a LASSO to a subset of OTUs selected with GLMM. We have conducted empirical studies using three simulation datasets generated using Dirichlet-multinomial models and a real gut microbiome data related to Parkinson's disease to investigate the performance of the three strategies for predictive analysis. Our simulation studies show that the predictive performance of LASSO with appropriate variable transformation works remarkably well on zero-inflated data. Our results of real data analysis show that Parkinson's disease can be predicted based on selected OTUs after the binary transformation, age, and sex with high accuracy (Error Rate = 0.199, AUC = 0.872, AUPRC = 0.912). These results provide strong evidences of the relationship between Parkinson's disease and the gut microbiome.

Eleven primer pairs were developed for the identification of Neisseria gonorrhoeae. The sensitivity and specificity of these primers were evaluated by Real Time (RT)-PCR melt curve analyses with DNA from 145 N. gonorrhoeae isolates and 40 other Neisseria or non-Neisseria species. Three primer pairs were further evaluated in a hydrogel-based RT-PCR detection platform, using DNA extracted from 50 N. gonorrhoeae cultures. We observed 100% sensitivity and specificity in the hydrogel assay, confirming its potential as a point-of-care test (POCT) for N. gonorrhoeae diagnosis.

Protein kinase-mediated phosphorylation is among the most important post-translational modifications. However, few phosphorylation sites have been experimentally identified for most species, making it difficult to determine the degree to which phosphorylation sites are conserved. The goal of this study was to use computational methods to characterize the conservation of human phosphorylation sites in a wide variety of eukaryotes. Using experimentally-determined human sites as input, homologous phosphorylation sites were predicted in all 432 eukaryotes for which complete proteomes were available. For each pair of species, we calculated phosphorylation site conservation as the number of phosphorylation sites found in both species divided by the number found in at least one of the two species. A clustering of the species based on this conservation measure was concordant with phylogenies based on traditional genomic measures. For a subset of the 432 species, phosphorylation site conservation was compared to conservation of both protein kinases and proteins in general. Protein kinases exhibited the highest degree of conservation, while general proteins were less conserved and phosphorylation sites were least conserved. Although preliminary, these data tentatively suggest that variation in phosphorylation sites may play a larger role in explaining phenotypic differences among organisms than differences in the complements of protein kinases or general proteins.

DNA microarrays and RNA sequencing (RNA-seq) are major technologies for performing high-throughput analysis of transcript abundance. Recently, concerns have been raised regarding the concordance of data derived from the two techniques. Using cDNA libraries derived from normal human foreskin fibroblasts, we measured changes in transcript abundance as cells transitioned from proliferative growth to quiescence using both DNA microarrays and RNA-seq. The internal reproducibility of the RNA-seq data was greater than that of the microarray data. Correlations between the RNA-seq data and the individual microarrays were low, but correlations between the RNA-seq values and the geometric mean of the microarray values were moderate. The two technologies had good agreement when considering probes with the largest (both positive and negative) fold change (FC) values. An independent technique, quantitative reverse-transcription PCR (qRT-PCR), was used to measure the FC of 76 genes between proliferative and quiescent samples, and a higher correlation was observed between the qRT-PCR data and the RNA-seq data than between the qRT-PCR data and the microarray data.

Thrombin and hypoxia are important players in breast cancer progression. Breast cancers often develop drug resistance, but mechanisms linking thrombin and hypoxia to drug resistance remain unresolved. Our studies using Doxorubicin (DOX) resistant MCF7 breast cancer cells reveals a mechanism linking DOX exposure with hypoxic induction of DOX resistance. Global expression changes between parental and DOX resistant MCF7 cells were examined. Westerns, Northerns and immunocytochemistry were used to validate drug resistance and differentially expressed genes. A cluster of genes involved in the anticoagulation pathway, with Tissue Factor Pathway Inhibitor 1 (TFPI1) the top hit, was identified. Plasmids overexpressing TFPI1 were utilized, and 1% O2 was used to test the effects of hypoxia on drug resistance. Lastly, microarray datasets from patients with drug resistant breast tumors were interrogated for TFPI1 expression levels. TFPI1 protein levels were found elevated in 3 additional DOX resistant cells lines, from humans and rats, indicating evolutionarily conservation of the effect. Elevated TFPI1 in DOX resistant cells was active, as thrombin protein levels were coincidentally low. We observed elevated HIF1α protein in DOX resistant cells, and in cells with forced expression of TFPI1, suggesting TFPI1 induces HIF1α. TFPI1 also induced c-MYC, c-SRC, and HDAC2 protein, as well as DOX resistance in parental cells. Growth of cells in 1% O2 induced elevated HIF1α, BCRP and MDR-1 protein, and these cells were resistant to DOX. Our in vitro results were consistent with in vivo patient datasets, as tumors harboring increased BCRP and MDR-1 expression also had increased TFPI1 expression. Our observations are clinically relevant indicating that DOX treatment induces an anticoagulation cascade, leading to inhibition of thrombin and the expression of HIF1α. This in turn activates a pathway leading to drug resistance.

Kinome microarrays are comprised of peptides that act as phosphorylation targets for protein kinases. This platform is growing in popularity due to its ability to measure phosphorylation-mediated cellular signaling in a high-throughput manner. While software for analyzing data from DNA microarrays has also been used for kinome arrays, differences between the two technologies and associated biologies previously led us to develop Platform for Intelligent, Integrated Kinome Analysis (PIIKA), a software tool customized for the analysis of data from kinome arrays. Here, we report the development of PIIKA 2, a significantly improved version with new features and improvements in the areas of clustering, statistical analysis, and data visualization. Among other additions to the original PIIKA, PIIKA 2 now allows the user to: evaluate statistically how well groups of samples cluster together; identify sets of peptides that have consistent phosphorylation patterns among groups of samples; perform hierarchical clustering analysis with bootstrapping; view false negative probabilities and positive and negative predictive values for t-tests between pairs of samples; easily assess experimental reproducibility; and visualize the data using volcano plots, scatterplots, and interactive three-dimensional principal component analyses. Also new in PIIKA 2 is a web-based interface, which allows users unfamiliar with command-line tools to easily provide input and download the results. Collectively, the additions and improvements described here enhance both the breadth and depth of analyses available, simplify the user interface, and make the software an even more valuable tool for the analysis of kinome microarray data. Both the web-based and stand-alone versions of PIIKA 2 can be accessed via http://saphire.usask.ca.

Numerous aspects of the relationship between bacteria and human have been investigated. One aspect that has recently received attention is sequence overlap at the proteomic level. However, there has not yet been a study that comprehensively characterizes the level of sequence overlap between bacteria and human, especially as it relates to bacterial characteristics like pathogenicity, G-C content, and proteome size. In this study, we began by performing a general characterization of the range of bacteria-human similarity at the proteomic level, and identified characteristics of the most- and least-similar bacterial species. We then examined the relationship between proteomic similarity and numerous other variables. While pathogens and nonpathogens had comparable similarity to the human proteome, pathogens causing chronic infections were found to be more similar to the human proteome than those causing acute infections. Although no general correspondence between a bacterium's proteome size and its similarity to the human proteome was noted, no bacteria with small proteomes had high similarity to the human proteome. Finally, we discovered an interesting relationship between similarity and a bacterium's G-C content. While the relationship between bacteria and human has been studied from many angles, their proteomic similarity still needs to be examined in more detail. This paper sheds further light on this relationship, particularly with respect to immunity and pathogenicity.

Certain short peptides do not occur in humans and are rare or non-existent in the universal proteome. Antigens that contain rare amino acid sequences are in general highly immunogenic and may activate different arms of the immune system. We first generated a list of rare, semi-common, and common 5-mer peptides using bioinformatics tools to analyze the UniProtKB database. Experimental observations indicated that rare and semi-common 5-mers generated stronger cellular responses in comparison with common-occurring sequences. We hypothesized that the biological process responsible for this enhanced immunogenicity could be used to positively modulate immune responses with potential application for vaccine development. Initially, twelve rare 5-mers, 9-mers, and 13-mers were incorporated in frame at the end of an H5N1 hemagglutinin (HA) antigen and expressed from a DNA vaccine. The presence of some 5-mer peptides induced improved immune responses. Adding one 5-mer peptide exogenously also offered improved clinical outcome and/or survival against a lethal H5N1 or H1N1 influenza virus challenge in BALB/c mice and ferrets, respectively. Interestingly, enhanced anti-HBsAg antibody production by up to 25-fold in combination with a commercial Hepatitis B vaccine (Engerix-B, GSK) was also observed in BALB/c mice. Mechanistically, NK cell activation and dependency was observed with enhancing peptides ex vivo and in NK-depleted mice. Overall, the data suggest that rare or non-existent oligopeptides can be developed as immunomodulators and supports the further evaluation of some 5-mer peptides as potential vaccine adjuvants.