Word-based GWAS harnesses the rich potential of genomic data for E. coli quinolone resistance.

Quinolone resistance presents a growing global health threat. We employed word-based GWAS to explore genomic data, aiming to enhance our understanding of this phenomenon. Unlike traditional variant-based GWAS analyses, this approach simultaneously captures multiple genomic factors, including single and interacting resistance mutations and genes. Analyzing a dataset of 92 genomic E. coli samples from a wastewater treatment plant in Dresden, we identified 54 DNA unitigs significantly associated with quinolone resistance. Remarkably, our analysis not only validated known mutations in gyrA and parC genes and the results of our variant-based GWAS but also revealed new (mutated) genes such as mdfA, the AcrEF-TolC multidrug efflux system, ptrB, and hisI, implicated in antibiotic resistance. Furthermore, our study identified joint mutations in 14 genes including the known gyrA gene, providing insights into potential synergistic effects contributing to quinolone resistance. These findings showcase the exceptional capabilities of word-based GWAS in unraveling the intricate genomic foundations of quinolone resistance.

Netrank: network-based approach for biomarker discovery.

BACKGROUND: Integrating multi-omics data is fast becoming a powerful approach for predicting disease progression and treatment outcomes. In light of that, we introduce a modified version of the NetRank algorithm, a network-based algorithm for biomarker discovery that incorporates the protein associations, co-expressions, and functions with its phenotypic association to differentiate different types of cancer. NetRank is introduced here as a robust feature selection method for biomarker selection in cancer prediction. We assess the robustness and suitability of the RNA gene expression data through scanning genomic data for 19 cancer types with more than 3000 patients from The Cancer Genome Atlas (TCGA). RESULTS: The results of evaluating different cancer type profiles from the TCGA data demonstrate the strength of our approach to identifying interpretable biomarker signatures for cancer outcome prediction. NetRank's biomarkers segregate most cancer types with an area under the curve (AUC) above 90% using compact signatures. CONCLUSION: In this paper we provide a fast and efficient implementation of NetRank, with a case study from The Cancer Genome Atlas, to assess the performance. We incorporated complete functionality for pre and post-processing for RNA-seq gene expression data with functions for building protein-protein interaction networks. The source code of NetRank is freely available (at github.com/Alfatlawi/Omics-NetRank) with an installable R library. We also deliver a comprehensive practical user manual with examples and data attached to this paper.

NetRank Recovers Known Cancer Hallmark Genes as Universal Biomarker Signature for Cancer Outcome Prediction.

Gene expression can serve as a powerful predictor for disease progression and other phenotypes. Consequently, microarrays, which capture gene expression genome-wide, have been used widely over the past two decades to derive biomarker signatures for tasks such as cancer grading, prognosticating the formation of metastases, survival, and others. Each of these signatures was selected and optimized for a very specific phenotype, tissue type, and experimental set-up. While all of these differences may naturally contribute to very heterogeneous and different biomarker signatures, all cancers share characteristics regardless of particular cell types or tissue as summarized in the hallmarks of cancer. These commonalities could give rise to biomarker signatures, which perform well across different phenotypes, cell and tissue types. Here, we explore this possibility by employing a network-based approach for pan-cancer biomarker discovery. We implement a random surfer model, which integrates interaction, expression, and phenotypic information to rank genes by their suitability for outcome prediction. To evaluate our approach, we assembled 105 high-quality microarray datasets sampled from around 13,000 patients and covering 13 cancer types. We applied our approach (NetRank) to each dataset and aggregated individual signatures into one compact signature of 50 genes. This signature stands out for two reasons. First, in contrast to other signatures of the 105 datasets, it is performant across nearly all cancer types and phenotypes. Second, It is interpretable, as the majority of genes are linked to the hallmarks of cancer in general and proliferation specifically. Many of the identified genes are cancer drivers with a known mutation burden linked to cancer. Overall, our work demonstrates the power of network-based approaches to compose robust, compact, and universal biomarker signatures for cancer outcome prediction.

A genome-wide scan of wastewater E. coli for genes under positive selection: focusing on mechanisms of antibiotic resistance.

Antibiotic resistance is a global health threat and consequently, there is a need to understand the mechanisms driving its emergence. Here, we hypothesize that genes and mutations under positive selection may contribute to antibiotic resistance. We explored wastewater E. coli, whose genomes are highly diverse. We subjected 92 genomes to a statistical analysis for positively selected genes. We obtained 75 genes under positive selection and explored their potential for antibiotic resistance. We found that eight genes have functions relating to antibiotic resistance, such as biofilm formation, membrane permeability, and bacterial persistence. Finally, we correlated the presence/absence of non-synonymous mutations in positively selected sites of the genes with a function in resistance against 20 most prescribed antibiotics. We identified mutations associated with antibiotic resistance in two genes: the porin ompC and the bacterial persistence gene hipA. These mutations are located at the surface of the proteins and may hence have a direct effect on structure and function. For hipA, we hypothesize that the mutations influence its interaction with hipB and that they enhance the capacity for dormancy as a strategy to evade antibiotics. Overall, genomic data and positive selection analyses uncover novel insights into mechanisms driving antibiotic resistance.

Mutations in bdcA and valS Correlate with Quinolone Resistance in Wastewater Escherichia coli.

Single mutations can confer resistance to antibiotics. Identifying such mutations can help to develop and improve drugs. Here, we systematically screen for candidate quinolone resistance-conferring mutations. We sequenced highly diverse wastewater E. coli and performed a genome-wide association study (GWAS) to determine associations between over 200,000 mutations and quinolone resistance phenotypes. We uncovered 13 statistically significant mutations including 1 located at the active site of the biofilm dispersal gene bdcA and 6 silent mutations in the aminoacyl-tRNA synthetase valS. The study also recovered the known mutations in the topoisomerases gyrase (gyrA) and topoisomerase IV (parC). In summary, we demonstrate that GWAS effectively and comprehensively identifies resistance mutations without a priori knowledge of targets and mode of action. The results suggest that mutations in the bdcA and valS genes, which are involved in biofilm dispersal and translation, may lead to novel resistance mechanisms.

Deep Learning Improves Pancreatic Cancer Diagnosis Using RNA-Based Variants.

For optimal pancreatic cancer treatment, early and accurate diagnosis is vital. Blood-derived biomarkers and genetic predispositions can contribute to early diagnosis, but they often have limited accuracy or applicability. Here, we seek to exploit the synergy between them by combining the biomarker CA19-9 with RNA-based variants. We use deep sequencing and deep learning to improve differentiating pancreatic cancer and chronic pancreatitis. We obtained samples of nucleated cells found in peripheral blood from 268 patients suffering from resectable, non-resectable pancreatic cancer, and chronic pancreatitis. We sequenced RNA with high coverage and obtained millions of variants. The high-quality variants served as input together with CA19-9 values to deep learning models. Our model achieved an area under the curve (AUC) of 96% in differentiating resectable cancer from pancreatitis using a test cohort. Moreover, we identified variants to estimate survival in resectable cancer. We show that the blood transcriptome harbours variants, which can substantially improve noninvasive clinical diagnosis.

Integrating evolutionary game theory into an agent-based model of ductal carcinoma in situ: Role of gap junctions in cancer progression.

BACKGROUND AND OBJECTIVE: There are many cells with various phenotypic behaviors in cancer interacting with each other. For example, an apoptotic cell may induce apoptosis in adjacent cells. A living cell can also protect cells from undergoing apoptosis and necrosis. These survival and death signals are propagated through interaction pathways between adjacent cells called gap junctions. The function of these signals depends on the cellular context of the cell receiving them. For instance, a receiver cell experiencing a low level of oxygen may interpret a received survival signal as an apoptosis signal. In this study, we examine the effect of these signals on tumor growth. METHODS: We make an evolutionary game theory component in order to model the signal propagation through gap junctions. The game payoffs are defined as a function of cellular context. Then, the game theory component is integrated into an agent-based model of tumor growth. After that, the integrated model is applied to ductal carcinoma in situ, a type of early stage breast cancer. Different scenarios are explored to observe the impact of the gap junction communication and parameters of the game theory component on cancer progression. We compare these scenarios by using the Wilcoxon signed-rank test. RESULTS: The Wilcoxon signed-rank test succeeds in proving a significant difference between the tumor growth of the model before and after considering the gap junction communication. The Wilcoxon signed-rank test also proves that the tumor growth significantly depends on the oxygen threshold of turning survival signals into apoptosis. CONCLUSIONS: In this study, the gap junction communication is modeled by using evolutionary game theory to illustrate its role at early stage cancers such as ductal carcinoma in situ. This work indicates that the gap junction communication and the oxygen threshold of turning survival signals into apoptosis can notably affect cancer progression.