Auto-phylo v2 and auto-phylo-pipeliner: building advanced, flexible, and reusable pipelines for phylogenetic inferences, estimation of variability levels and identification of positively selected amino acid sites.

The vast amount of genome sequence data that is available, and that is predicted to drastically increase in the near future, can only be efficiently dealt with by building automated pipelines. Indeed, the Earth Biogenome Project will produce high-quality reference genome sequences for all 1.8 million named living eukaryote species, providing unprecedented insight into the evolution of genes and gene families, and thus on biological issues. Here, new modules for gene annotation, further BLAST search algorithms, further multiple sequence alignment methods, the adding of reference sequences, further tree rooting methods, the estimation of rates of synonymous and nonsynonymous substitutions, and the identification of positively selected amino acid sites, have been added to auto-phylo (version 2), a recently developed software to address biological problems using phylogenetic inferences. Additionally, we present auto-phylo-pipeliner, a graphical user interface application that further facilitates the creation and running of auto-phylo pipelines. Inferences on S-RNase specificity, are critical for both cross-based breeding and for the establishment of pollination requirements. Therefore, as a test case, we develop an auto-phylo pipeline to identify amino acid sites under positive selection, that are, in principle, those determining S-RNase specificity, starting from both non-annotated Prunus genomes and sequences available in public databases.

Multiple Lines of Evidence Support 199 SARS-CoV-2 Positively Selected Amino Acid Sites.

SARS-CoV-2 amino acid variants that contribute to an increased transmissibility or to host immune system escape are likely to increase in frequency due to positive selection and may be identified using different methods, such as codeML, FEL, FUBAR, and MEME. Nevertheless, when using different methods, the results do not always agree. The sampling scheme used in different studies may partially explain the differences that are found, but there is also the possibility that some of the identified positively selected amino acid sites are false positives. This is especially important in the context of very large-scale projects where hundreds of analyses have been performed for the same protein-coding gene. To account for these issues, in this work, we have identified positively selected amino acid sites in SARS-CoV-2 and 15 other coronavirus species, using both codeML and FUBAR, and compared the location of such sites in the different species. Moreover, we also compared our results to those that are available in the COV2Var database and the frequency of the 10 most frequent variants and predicted protein location to identify those sites that are supported by multiple lines of evidence. Amino acid changes observed at these sites should always be of concern. The information reported for SARS-CoV-2 can also be used to identify variants of concern in other coronaviruses.

Proteomic analysis of chromophobe renal cell carcinoma and benign renal oncocytoma biopsies reveals shared metabolic dysregulation.

BACKGROUND: This study investigates the proteomic landscapes of chromophobe renal cell carcinoma (chRCC) and renal oncocytomas (RO), two subtypes of renal cell carcinoma that together account for approximately 10% of all renal tumors. Despite their histological similarities and shared origins, chRCC is a malignant tumor necessitating aggressive intervention, while RO, a benign growth, is often subject to overtreatment due to difficulties in accurate differentiation. METHODS: We conducted a label-free quantitative proteomic analysis on solid biopsies of chRCC (n = 5), RO (n = 5), and normal adjacent tissue (NAT, n = 5). The quantitative analysis was carried out by comparing protein abundances between tumor and NAT specimens. Our analysis identified a total of 1610 proteins across all samples, with 1379 (85.7%) of these proteins quantified in at least seven out of ten LC‒MS/MS runs for one renal tissue type (chRCC, RO, or NAT). RESULTS: Our findings revealed significant similarities in the dysregulation of key metabolic pathways, including carbohydrate, lipid, and amino acid metabolism, in both chRCC and RO. Compared to NAT, both chRCC and RO showed a marked downregulation in gluconeogenesis proteins, but a significant upregulation of proteins integral to the citrate cycle. Interestingly, we observed a distinct divergence in the oxidative phosphorylation pathway, with RO showing a significant increase in the number and degree of alterations in proteins, surpassing that observed in chRCC. CONCLUSIONS: This study underscores the value of integrating high-resolution mass spectrometry protein quantification to effectively characterize and differentiate the proteomic landscapes of solid tumor biopsies diagnosed as chRCC and RO. The insights gained from this research offer valuable information for enhancing our understanding of these conditions and may aid in the development of improved diagnostic and therapeutic strategies.

MyBrain-Seq: A Pipeline for MiRNA-Seq Data Analysis in Neuropsychiatric Disorders.

High-throughput sequencing of small RNA molecules such as microRNAs (miRNAs) has become a widely used approach for studying gene expression and regulation. However, analyzing miRNA-Seq data can be challenging because it requires multiple steps, from quality control and preprocessing to differential expression and pathway-enrichment analyses, with many tools and databases available for each step. Furthermore, reproducibility of the analysis pipeline is crucial to ensure that the results are accurate and reliable. Here, we present myBrain-Seq, a comprehensive and reproducible pipeline for analyzing miRNA-Seq data that incorporates miRNA-specific solutions at each step of the analysis. The pipeline was designed to be flexible and user-friendly, allowing researchers with different levels of expertise to perform the analysis in a standardized and reproducible manner, using the most common and widely used tools for each step. In this work, we describe the implementation of myBrain-Seq and demonstrate its capacity to consistently and reproducibly identify differentially expressed miRNAs and enriched pathways by applying it to a real case study in which we compared schizophrenia patients who responded to medication with treatment-resistant schizophrenia patients to obtain a 16-miRNA treatment-resistant schizophrenia profile.

PanDrugs2: prioritizing cancer therapies using integrated individual multi-omics data.

Genomics studies routinely confront researchers with long lists of tumor alterations detected in patients. Such lists are difficult to interpret since only a minority of the alterations are relevant biomarkers for diagnosis and for designing therapeutic strategies. PanDrugs is a methodology that facilitates the interpretation of tumor molecular alterations and guides the selection of personalized treatments. To do so, PanDrugs scores gene actionability and drug feasibility to provide a prioritized evidence-based list of drugs. Here, we introduce PanDrugs2, a major upgrade of PanDrugs that, in addition to somatic variant analysis, supports a new integrated multi-omics analysis which simultaneously combines somatic and germline variants, copy number variation and gene expression data. Moreover, PanDrugs2 now considers cancer genetic dependencies to extend tumor vulnerabilities providing therapeutic options for untargetable genes. Importantly, a novel intuitive report to support clinical decision-making is generated. PanDrugs database has been updated, integrating 23 primary sources that support >74K drug-gene associations obtained from 4642 genes and 14 659 unique compounds. The database has also been reimplemented to allow semi-automatic updates to facilitate maintenance and release of future versions. PanDrugs2 does not require login and is freely available at https://www.pandrugs.org/.

Negative Samples for Improving Object Detection-A Case Study in AI-Assisted Colonoscopy for Polyp Detection.

Deep learning object-detection models are being successfully applied to develop computer-aided diagnosis systems for aiding polyp detection during colonoscopies. Here, we evidence the need to include negative samples for both (i) reducing false positives during the polyp-finding phase, by including images with artifacts that may confuse the detection models (e.g., medical instruments, water jets, feces, blood, excessive proximity of the camera to the colon wall, blurred images, etc.) that are usually not included in model development datasets, and (ii) correctly estimating a more realistic performance of the models. By retraining our previously developed YOLOv3-based detection model with a dataset that includes 15% of additional not-polyp images with a variety of artifacts, we were able to generally improve its F1 performance in our internal test datasets (from an average F1 of 0.869 to 0.893), which now include such type of images, as well as in four public datasets that include not-polyp images (from an average F1 of 0.695 to 0.722).

On the identification of potential novel therapeutic targets for spinocerebellar ataxia type 1 (SCA1) neurodegenerative disease using EvoPPI3.

EvoPPI (http://evoppi.i3s.up.pt), a meta-database for protein-protein interactions (PPI), has been upgraded (EvoPPI3) to accept new types of data, namely, PPI from patients, cell lines, and animal models, as well as data from gene modifier experiments, for nine neurodegenerative polyglutamine (polyQ) diseases caused by an abnormal expansion of the polyQ tract. The integration of the different types of data allows users to easily compare them, as here shown for Ataxin-1, the polyQ protein involved in spinocerebellar ataxia type 1 (SCA1) disease. Using all available datasets and the data here obtained for Drosophila melanogaster wt and exp Ataxin-1 mutants (also available at EvoPPI3), we show that, in humans, the Ataxin-1 network is much larger than previously thought (380 interactors), with at least 909 interactors. The functional profiling of the newly identified interactors is similar to the ones already reported in the main PPI databases. 16 out of 909 interactors are putative novel SCA1 therapeutic targets, and all but one are already being studied in the context of this disease. The 16 proteins are mainly involved in binding and catalytic activity (mainly kinase activity), functional features already thought to be important in the SCA1 disease.

MiRNA Differences Related to Treatment-Resistant Schizophrenia.

Schizophrenia (SZ) is a serious mental disorder that is typically treated with antipsychotic medication. Treatment-resistant schizophrenia (TRS) is the condition where symptoms remain after pharmacological intervention, resulting in long-lasting functional and social impairments. As the identification and treatment of a TRS patient requires previous failed treatments, early mechanisms of detection are needed in order to quicken the access to effective therapy, as well as improve treatment adherence. In this study, we aim to find a microRNA (miRNA) signature for TRS, as well as to shed some light on the molecular pathways potentially involved in this severe condition. To do this, we compared the blood miRNAs of schizophrenia patients that respond to medication and TRS patients, thus obtaining a 16-miRNA TRS profile. Then, we assessed the ability of this signature to separate responders and TRS patients using hierarchical clustering, observing that most of them are grouped correctly (~70% accuracy). We also conducted a network, pathway analysis, and bibliography search to spot molecular pathways potentially altered in TRS. We found that the response to stress seems to be a key factor in TRS and that proteins p53, SIRT1, MDM2, and TRIM28 could be the potential mediators of such responses. Finally, we suggest a molecular pathway potentially regulated by the miRNAs of the TRS profile.

Pathway-guided monitoring of the disease course in bladder cancer with longitudinal urine proteomics.

BACKGROUND: Monitoring bladder cancer over time requires invasive and costly procedures. Less invasive approaches are required using readily available biological samples such as urine. In this study, we demonstrate a method for longitudinal analysis of the urine proteome to monitor the disease course in patients with bladder cancer. METHODS: We compared the urine proteomes of patients who experienced recurrence and/or progression (n = 13) with those who did not (n = 17). We identified differentially expressed proteins within various pathways related to the hallmarks of cancer. The variation of such pathways during the disease course was determined using our differential personal pathway index (dPPi) calculation, which could indicate disease progression and the need for medical intervention. RESULTS: Seven hallmark pathways are used to develop the dPPi. We demonstrate that we can successfully longitudinally monitor the disease course in bladder cancer patients through a combination of urine proteomic analysis and the dPPi calculation, over a period of 62 months. CONCLUSIONS: Using the information contained in the patient's urinary proteome, the dPPi reflects the individual's course of bladder cancer, and helps to optimise the use of more invasive procedures such as cystoscopy.

Bladder cancer must be closely monitored for progression, but this requires expensive and invasive procedures such as cystoscopy. Less invasive procedures using readily available samples such as urine are needed. Here, we present an approach that measures the levels of various proteins in the urine. We compare protein levels at different points during the disease course in patients with bladder cancer, and show this helps to flag disease recurrence and the need for medical intervention. Our approach could help clinicians to determine which patients require more invasive testing and treatment.

Performance of Convolutional Neural Networks for Polyp Localization on Public Colonoscopy Image Datasets.

Colorectal cancer is one of the most frequent malignancies. Colonoscopy is the de facto standard for precancerous lesion detection in the colon, i.e., polyps, during screening studies or after facultative recommendation. In recent years, artificial intelligence, and especially deep learning techniques such as convolutional neural networks, have been applied to polyp detection and localization in order to develop real-time CADe systems. However, the performance of machine learning models is very sensitive to changes in the nature of the testing instances, especially when trying to reproduce results for totally different datasets to those used for model development, i.e., inter-dataset testing. Here, we report the results of testing of our previously published polyp detection model using ten public colonoscopy image datasets and analyze them in the context of the results of other 20 state-of-the-art publications using the same datasets. The F1-score of our recently published model was 0.88 when evaluated on a private test partition, i.e., intra-dataset testing, but it decayed, on average, by 13.65% when tested on ten public datasets. In the published research, the average intra-dataset F1-score is 0.91, and we observed that it also decays in the inter-dataset setting to an average F1-score of 0.83.

SEDA: A Desktop Tool Suite for FASTA Files Processing.

SEDA (SEquence DAtaset builder) is a multiplatform desktop application for the manipulation of FASTA files containing DNA or protein sequences. The convenient graphical user interface gives access to a collection of simple (filtering, sorting, or file reformatting, among others) and advanced (BLAST searching, protein domain annotation, gene annotation, and sequence alignment) utilities not present in similar applications, which eases the work of life science researchers working with DNA and/or protein sequences, especially those who have no programming skills. This paper presents general guidelines on how to build efficient data handling protocols using SEDA, as well as practical examples on how to prepare high-quality datasets for single gene phylogenetic studies, the characterization of protein families, or phylogenomic studies. The user-friendliness of SEDA also relies on two important features: (i) the availability of easy-to-install distributable versions and installers of SEDA, including a Docker image for Linux, and (ii) the facility with which users can manage large datasets. SEDA is open-source, with GNU General Public License v3.0 license, and publicly available at GitHub (https://github.com/sing-group/seda). SEDA installers and documentation are available at https://www.sing-group.org/seda/.

Compi: a framework for portable and reproducible pipelines.

Compi is an application framework to develop end-user, pipeline-based applications with a primary emphasis on: (i) user interface generation, by automatically generating a command-line interface based on the pipeline specific parameter definitions; (ii) application packaging, with compi-dk, which is a version-control-friendly tool to package the pipeline application and its dependencies into a Docker image; and (iii) application distribution provided through a public repository of Compi pipelines, named Compi Hub, which allows users to discover, browse and reuse them easily. By addressing these three aspects, Compi goes beyond traditional workflow engines, having been specially designed for researchers who want to take advantage of common workflow engine features (such as automatic job scheduling or logging, among others) while keeping the simplicity and readability of shell scripts without the need to learn a new programming language. Here we discuss the design of various pipelines developed with Compi to describe its main functionalities, as well as to highlight the similarities and differences with similar tools that are available. An open-source distribution under the Apache 2.0 License is available from GitHub (available at https://github.com/sing-group/compi). Documentation and installers are available from https://www.sing-group.org/compi. A specific repository for Compi pipelines is available from Compi Hub (available at https://www.sing-group.org/compihub.

Application of miRNA-seq in neuropsychiatry: A methodological perspective.

MiRNAs are emerging as key molecules to study neuropsychiatric diseases. However, despite the large number of methodologies and software for miRNA-seq analyses, there is little supporting literature for researchers in this area. This review focuses on evaluating how miRNA-seq has been used to study neuropsychiatric diseases to date, analyzing both the main findings discovered and the bioinformatics workflows and tools used from a methodological perspective. The objective of this review is two-fold: first, to evaluate current miRNA-seq procedures used in neuropsychiatry; and second, to offer comprehensive information that can serve as a guide to new researchers in bioinformatics. After conducting a systematic search (from 2016 to June 30, 2020) of articles using miRNA-seq in neuropsychiatry, we have seen that it has already been used for different types of studies in three main categories: diagnosis, prognosis, and mechanism. We carefully analyzed the bioinformatics workflows of each study, observing a high degree of variability with respect to the tools and methods used and several methodological complexities that are identified and discussed in this review.

On the Identification of Clinically Relevant Bacterial Amino Acid Changes at the Whole Genome Level Using Auto-PSS-Genome.

The identification of clinically relevant bacterial amino acid changes can be performed using different methods aimed at the identification of genes showing positively selected amino acid sites (PSS). Nevertheless, such analyses are time consuming, and the frequency of genes showing evidence for PSS can be low. Therefore, the development of a pipeline that allows the quick and efficient identification of the set of genes that show PSS is of interest. Here, we present Auto-PSS-Genome, a Compi-based pipeline distributed as a Docker image, that automates the process of identifying genes that show PSS using three different methods, namely codeML, FUBAR, and omegaMap. Auto-PSS-Genome accepts as input a set of FASTA files, one per genome, containing all coding sequences, thus minimizing the work needed to conduct positively selected sites analyses. The Auto-PSS-Genome pipeline identifies orthologous gene sets and corrects for multiple possible problems in input FASTA files that may prevent the automated identification of genes showing PSS. A FASTA file containing all coding sequences can also be given as an external global reference, thus easing the comparison of results across species, when gene names are different. In this work, we use Auto-PSS-Genome to analyse Mycobacterium leprae (that causes leprosy), and the closely related species M. haemophilum, that mainly causes ulcerating skin infections and arthritis in persons who are severely immunocompromised, and in children causes cervical and perihilar lymphadenitis. The genes identified in these two species as showing PSS may be those that are partially responsible for virulence and resistance to drugs.

DREIMT: a drug repositioning database and prioritization tool for immunomodulation.

MOTIVATION: Drug immunomodulation modifies the response of the immune system and can be therapeutically exploited in pathologies such as cancer and autoimmune diseases. RESULTS: DREIMT is a new hypothesis-generation web tool, which performs drug prioritization analysis for immunomodulation. DREIMT provides significant immunomodulatory drugs targeting up to 70 immune cells subtypes through a curated database that integrates 4960 drug profiles and ∼2600 immune gene expression signatures. The tool also suggests potential immunomodulatory drugs targeting user-supplied gene expression signatures. Final output includes drug-signature association scores, FDRs and downloadable plots and results tables. AVAILABILITYAND IMPLEMENTATION: http://www.dreimt.org. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

An RNA-seq Based Machine Learning Approach Identifies Latent Tuberculosis Patients With an Active Tuberculosis Profile.

A better understanding of the response against Tuberculosis (TB) infection is required to accurately identify the individuals with an active or a latent TB infection (LTBI) and also those LTBI patients at higher risk of developing active TB. In this work, we have used the information obtained from studying the gene expression profile of active TB patients and their infected -LTBI- or uninfected -NoTBI- contacts, recruited in Spain and Mozambique, to build a class-prediction model that identifies individuals with a TB infection profile. Following this approach, we have identified several genes and metabolic pathways that provide important information of the immune mechanisms triggered against TB infection. As a novelty of our work, a combination of this class-prediction model and the direct measurement of different immunological parameters, was used to identify a subset of LTBI contacts (called TB-like) whose transcriptional and immunological profiles are suggestive of infection with a higher probability of developing active TB. Validation of this novel approach to identifying LTBI individuals with the highest risk of active TB disease merits further longitudinal studies on larger cohorts in TB endemic areas.

Metatax: Metataxonomics with a Compi-Based Pipeline for Precision Medicine.

The human body immune system, metabolism and homeostasis are affected by microbes. Dysbiosis occurs when the homeostatic equilibrium is disrupted due to an alteration in the normal microbiota of the intestine. Dysbiosis can cause cancer, and also affect a patient's ability to respond to treatment. Metataxonomics seeks to identify the bacteria present in a biological sample, based on the sequencing of the 16S rRNA genetic marker. Precision medicine attempts to find relationships between the microbiota and the risk of acquiring cancer, and design new therapies targeting bacteria. Flexible and portable bioinformatic pipelines are necessary to be able to bring metataxonomics to the clinical field, which allow groups of biological samples to be classified according to their diversity in the microbiota. With this aim we implemented Metatax, a new pipeline to analyze biological samples based on 16S rRNA gene sequencing. The results obtained with our pipeline should complement those obtained by sequencing a patient's DNA and RNA, in addition to clinical data, to improve knowledge of the possible reasons for a disease or a worse response to treatment.

Ultrasonic-assisted extraction and digestion of proteins from solid biopsies followed by peptide sequential extraction hyphenated to MALDI-based profiling holds the promise of distinguishing renal oncocytoma from chromophobe renal cell carcinoma.

A novel analytical approach is proposed to discriminate between solid biopsies of chromophobe renal cell carcinoma (chRCC) and renal oncocytoma (RO). The method comprises the following steps: (i) ultrasonic extraction of proteins from solid biopsies, (ii) protein depletion with acetonitrile, (iii) ultrasonic assisted in-solution digestion using magnetic nanoparticle with immobilized trypsin, (iv) C18 tip-based preconcentration of peptides, (v) sequential extraction of the peptides with ACN, (vi) MALDI-snapshot of the extracts and (vii) investigation of the extract containing the most discriminating features using high resolution mass spectrometry. With this approach we have been able to differentially cluster renal oncocytoma and chromophobe renal cell carcinoma and identified 18 proteins specific to chromophobe and seven unique to renal oncocytoma. Chromophobes express proteins associated with ATP function (ATP5I & 5E; VATE1 & G2; ADT2), glycolysis (PGK1) and neuromedin whilst oncocytomas express ATP5H, ATPA, DEPD7 and TRIPB thyroid receptor interacting protein.

ATXN1 N-terminal region explains the binding differences of wild-type and expanded forms.

BACKGROUND: Wild-type (wt) polyglutamine (polyQ) regions are implicated in stabilization of protein-protein interactions (PPI). Pathological polyQ expansion, such as that in human Ataxin-1 (ATXN1), that causes spinocerebellar ataxia type 1 (SCA1), results in abnormal PPI. For ATXN1 a larger number of interactors has been reported for the expanded (82Q) than the wt (29Q) protein. METHODS: To understand how the expanded polyQ affects PPI, protein structures were predicted for wt and expanded ATXN1, as well as, for 71 ATXN1 interactors. Then, the binding surfaces of wt and expanded ATXN1 with the reported interactors were inferred. RESULTS: Our data supports that the polyQ expansion alters the ATXN1 conformation and that it enhances the strength of interaction with ATXN1 partners. For both ATXN1 variants, the number of residues at the predicted binding interface are greater after the polyQ, mainly due to the AXH domain. Moreover, the difference in the interaction strength of the ATXN1 variants was due to an increase in the number of interactions at the N-terminal region, before the polyQ, for the expanded form. CONCLUSIONS: There are three regions at the AXH domain that are essential for ATXN1 PPI. The N-terminal region is responsible for the strength of the PPI with the ATXN1 variants. How the predicted motifs in this region affect PPI is discussed, in the context of ATXN1 post-transcriptional modifications.

Predicting Specificities Under the Non-self Gametophytic Self-Incompatibility Recognition Model.

Non-self gametophytic self-incompatibility (GSI) recognition system is characterized by the presence of multiple F-box genes tandemly located in the S-locus, that regulate pollen specificity. This reproductive barrier is present in Solanaceae, Plantaginacea and Maleae (Rosaceae), but only in Petunia functional assays have been performed to get insight on how this recognition mechanism works. In this system, each of the encoded S-pollen proteins (called SLFs in Solanaceae and Plantaginaceae /SFBBs in Maleae) recognizes and interacts with a sub-set of non-self S-pistil proteins, called S-RNases, mediating their ubiquitination and degradation. In Petunia there are 17 SLF genes per S-haplotype, making impossible to determine experimentally each SLF specificity. Moreover, domain -swapping experiments are unlikely to be performed in large scale to determine S-pollen and S-pistil specificities. Phylogenetic analyses of the Petunia SLFs and those from two Solanum genomes, suggest that diversification of SLFs predate the two genera separation. Here we first identify putative SLF genes from nine Solanum and 10 Nicotiana genomes to determine how many gene lineages are present in the three genera, and the rate of origin of new SLF gene lineages. The use of multiple genomes per genera precludes the effect of incompleteness of the genome at the S-locus. The similar number of gene lineages in the three genera implies a comparable effective population size for these species, and number of specificities. The rate of origin of new specificities is one per 10 million years. Moreover, here we determine the amino acids positions under positive selection, those involved in SLF specificity recognition, using 10 Petunia S-haplotypes with more than 11 SLF genes. These 16 amino acid positions account for the differences of self-incompatible (SI) behavior described in the literature. When SLF and S-RNase proteins are divided according to the SI behavior, and the positively selected amino acids classified according to hydrophobicity, charge, polarity and size, we identified fixed differences between SI groups. According to the in silico 3D structure of the two proteins these amino acid positions interact. Therefore, this methodology can be used to infer SLF/S-RNase specificity recognition.