Transcription factors (TFs) are crucial in modulating gene expression and sculpting cellular and organismal phenotypes. The identification of TF-target gene interactions is pivotal for comprehending molecular pathways and disease etiologies but has been hindered by the demanding nature of traditional experimental approaches. This paper introduces a novel web application, utilizing the R programming language, which predicts TF-target gene relationships and vice versa. Our application integrates the predictive power of various bioinformatic tools, leveraging their combined strengths to provide robust predictions. It merges databases for enhanced precision, incorporates gene expression correlation for accuracy, and employs pan-tissue correlation analysis for context-specific insights. The application also enables the integration of user data with established resources to analyze TF-target gene networks. Despite its current limitation to human data, it provides a platform for exploring gene regulatory mechanisms comprehensively. This integrated, systematic approach offers researchers an invaluable tool for dissecting the complexities of gene regulation, with the potential for future expansions to include a broader range of species.

Understanding gene expression changes in chicks after vaccination against Newcastle Disease (ND) can reveal vaccine biomarkers. There is limited data on chicks’ early immune response after ND vaccination. Two trials aimed into this knowledge gap. In experiment one, forty-two 13 days-old specific pathogen free (SPF) chicks were used. Harderian glands (Hg) and tracheas (Tc) from five birds per group were sampled at 12, 24, and 48 hours post-vaccination (hpv) to evaluate the gene transcription levels by RNA sequencing (RNA-seq) and RT-qPCR. Results of RNA-seq were compared by glmFTest, while results of RT-qPCR were compared by t-test. With RNA-seq, significant up-regulation of interferon related genes along with JAK-STAT signaling pathways regulation was observed in the Hg at 24- hpv. None of the DEGs identified by RNA-seq were positive for RT-qPCR. Experiment 2 used one-hundred and twelve SPF and commercial chickens at 1 day-old and 14 days-old. Only the commercial birds had maternal antibodies for NDV. By RNA-seq, twenty core DEGs associated with innate immunity and viral genome replication inhibition were identified. Genes previously unlinked to NDV response, such as USP41, were identified. This research present genes with potential as immunity biomarkers for vaccines, yet further investigation is needed to correlate the core gene expression with viral shedding post-vaccination.

Utilizing bioinformatics tools, this study expands our understanding of secondary metabolism in Botrytis cinerea, identifying novel genes within polyketide synthase (PKS), non-ribosomal peptide synthetase (NRPS), sesquiterpene cyclase (STC), diterpene cyclase (DTC), and dimethylal-lyltryptophan synthase (DMATS) families. These findings enrich the genetic framework associ-ated with B. cinerea's pathogenicity and ecological adaptation, offering insights into uncharted metabolic pathways. Significantly, the discovery of previously unannotated genes provides new molecular targets for developing targeted antifungal strategies, promising to enhance crop pro-tection and advance our understanding of fungal biochemistry. This research not only broadens the scope of known secondary metabolites but also opens avenues for future exploration into B. cinerea's biosynthetic capabilities, potentially leading to novel antifungal compounds. Our work underscores the importance of integrating bioinformatics and genomics for fungal research, paving the way for sustainable agricultural practices by pinpointing precise molecular interven-tions against B. cinerea. This study sets a foundation for further investigations into the fungus's secondary metabolism, with implications for biotechnology and crop disease management.

The study of rare diseases is important not only for the individuals affected but also for the advancement of medical knowledge and a deeper understanding of human biology and genetics. The wide repertoire of structural information now available from reliable and accurate prediction methods provides the opportunity to investigate the molecular origins of most of the rare diseases reviewed in the Orpha.net database. Thus, it has been possible to analyze the topology of the missense mutations found in the 2,535 proteins involved in Mendelian rare diseases (MRD), which form the database for our structural bioinformatics study. The amino acid substitutions responsible for MRD show different mutation site distributions at different three-dimensional protein depths. We then predicted the depth-dependent effects of mutations for the 20,248 pathogenic variants that are present in our database. The results of this structural bioinformatics investigation are relevant, as they provide additional clues to mitigate the damage caused by MRD.

DNA Topoisomerase IIa (Top2A) is a nuclear enzyme that is a cancer drug target, and there is interest in identifying novel sites on the enzyme to inhibit cancer cells more selectively and to reduce off-target toxicity. The C-terminal domain (CTD) is one potential target, but it is an intrinsically disordered domain, which prevents structural analysis. Therefore, we set out to analyze the sequence of Top2A from 105 species using bioinformatic analysis including the PSICalc algorithm, entropy analysis, and other approaches. Our results demonstrate that large (10th-order) interdependent clusters are found including non-proximal positions across the major domains of Top2A. Further, CTD-specific clusters of 3rd, 4th, and 5th-order were identified that included positions that had been previously analyzed via mutation and biochemical assays. Some of these clusters coincided with positions that when mutated either increased or decreased relaxation activity. Finally, sites of low entropy (i.e., high information content) were identified and mapped as key positions in the CTD. Included in the low entropy sites are phosphorylation sites and charged positions. Together, these results help to build a clearer picture of the critical positions in the CTD and provide potential sites/regions for further analysis.

Recent advancements in molecular diagnostics have highlighted the potential of circulating cell-free DNA (ccfDNA) as a biomarker in various diseases. This study aims to explore the role of ccfDNA in the early diagnosis, prognosis, and therapeutic monitoring of diseases, with a focus on oncology. We conducted a comprehensive review and analysis of current literature, combined with experimental approaches involving advanced genomic sequencing and bioinformatics techniques. Our research systematically characterizes the molecular composition of ccfDNA, including its origins, size distribution, and methylation patterns, across different disease states. Furthermore, we investigate the dynamics of ccfDNA release and clearance in the bloodstream, elucidating mechanisms that could influence its diagnostic accuracy. In the context of cancer, we analyze the correlation between tumor-derived ccfDNA and disease progression, response to therapy, and overall survival. Our findings indicate that ccfDNA holds promise as a non-invasive biomarker for early disease detection, real-time monitoring of treatment efficacy, and potentially, predicting relapse. This study contributes to the growing field of liquid biopsies and opens new avenues for personalized medicine, emphasizing the need for standardized protocols in ccfDNA analysis to maximize its clinical utility. The implications of our findings extend beyond oncology, suggesting potential applications of ccfDNA in other pathological conditions.

This review article offers an in-depth look at DNA, its discovery, structure, functions, and the significant role it plays in modern biomedical research. It begins by revisiting the landmark discovery of the DNA double helix by scientists Watson and Crick, a structure resembling a twisted ladder that forms the blueprint of life. The review explains the basic processes of DNA replication (copying DNA), transcription (converting DNA to RNA), and translation (turning RNA into proteins), which are fundamental to understanding how genetic information is passed on and used by living organisms. The article then explores gene regulation and epigenetics, which involve how genes are turned on or off and how environmental factors can affect gene expression without changing the DNA sequence itself. It examines how DNA interacts with various proteins and modifiers, which play crucial roles in these processes. A major focus is on the advancements in DNA sequencing technologies. From the initial Sanger sequencing method, which was labor-intensive and time-consuming, to the highly efficient next-generation sequencing (NGS) and the latest single-molecule sequencing techniques, the review covers how these technologies have revolutionized genomics and personalized medicine. It also highlights the importance of computational genomics and bioinformatics, fields that deal with analyzing and interpreting the vast amounts of data generated by DNA sequencing. The review discusses DNA's crucial role in understanding human diseases. It covers DNA-based research in cancer, rare genetic disorders, infectious diseases, and pharmacogenomics (how genes affect a person's response to drugs). It emphasizes studies that have led to new ways of diagnosing and treating diseases. The article looks at the future of DNA research. It talks about emerging technologies like CRISPR-Cas9, a revolutionary gene-editing tool that allows scientists to modify DNA with high precision. Synthetic biology, an interdisciplinary field that combines biology and engineering for designing and constructing new biological parts and systems, is also discussed. Furthermore, DNA nanotechnology, which involves designing DNA molecules for technological uses, is explored. Finally, the review addresses ethical issues in DNA research, such as privacy concerns and debates over genome editing.

This research ventures into the critical public health challenge of childhood obesity by exploring the dynamic interplay between psychological well-being and Body Mass Index (BMI) throughout various developmental stages of childhood. It delves into how emotional regulation, attachment dynamics, and social relationships correlate with obesity from early childhood to adolescence. Highlighting key findings, such as the negative correlation between psychological resilience and higher BMI in young children, the impact of social relationships on obesity risk during pre-adolescence, and the link between adaptive emotional strategies and higher BMI in adolescents, this study brings to the fore the nuanced relationship between psychological factors and obesity. Utilizing the interdisciplinary process of bioinformatics, this research synergizes psychometric and biomedical data to unearth psychological markers critical for crafting targeted, age-appropriate interventions. This study advocates for a holistic healthcare approach, emphasizing the integration of psychological support within obesity prevention and management strategies, thereby underscoring the indispensable role of psychological factors in the fight against childhood obesity. The application of bioinformatics methods to analyze complex datasets demonstrates how collaboration across medical specialties can enrich our understanding and response to childhood obesity, contributing significantly to the development of comprehensive, bioinformatics-enhanced healthcare solutions.

Bioinformatics is a rapidly developing field enabling scientific experiments through computer models and simulations. Considering the vast databases of biological data available, it is ex-tremely important to develop efficient methods and algorithms for their processing. Sequence comparison is the best method for studying the evolutionary interaction between genes. It is based on alignment – the process of arranging two or more sequences to achieve the maximum level of identity and degree of similarity. The paper presents a new version of the algorithm for pairwise DNA sequences alignment, based on a new method called CAT, where a dependency with a previous match and the closes neighbor are taken in consideration to increase uniqueness of the CAT profile and to reduce possible collisions, i.e. two or more sequence having same CAT profiles. This makes proposed algorithm suitable for finding exact match of a concrete DNA se-quence in a big set of DNA data faster. The generation of CAT profiles is made once before data has been uploaded to the database, allowing the profiles to be used as metadata for the sequenc-es. It consists of an algorithm to calculate a CAT profile against the selected reference sequences and an algorithm to compare two sequences based on the calculated CAT profiles. Improve-ments in generation of the CAT profiles, are detailed described in the paper. Block scheme, pseudo code tables and figures were updated according to the proposed new version and ex-perimental results. Experiments have been carried out with the new version of the method and different datasets to align DNA sequences based on the CAT method. New experimental results in terms of collisions, speed, and efficiency of the proposed solutions are presented. Experiments related to the performance comparison with Needleman-Wunsch were re-executed with the new version of the algorithm to confirm that we have the same performance. An analysis of the per-formance of the proposed CAT based algorithm against Knuth–Morris–Pratt algorithm, which has a complexity of O(n) and is one of the most widely used for searching biological data, was performed. The impact of prior matching dependencies on uniqueness for generated CAT pro-files is investigated. The analysis of the experimental results obtained by sequence alignment shows a small deviation of the proposed algorithm based on the CAT method, which can be ig-nored if this deviation is acceptable at the expense of performance. The time efficiency of the CAT algorithm remains constant, regardless of the length of the sequences. Therefore, the ad-vantage of the proposed method is its fast processing in the alignment of large sequences, for which the execution of the exact algorithms takes a long time. The example code realization of the CAT Method, under the protection of the GNU General Public License v3.0, can be accessed on GitHub at: https://github.com/HristoS/CATSequenceAnalysis.

Metabolites are at the end of the gene-transcript-protein-metabolism cascade. As such, metabo-lomics is the omic approach that offers the most direct correlation with phenotype. This allows that, where genomics, transcriptomics and proteomics fail to explain a trait, metabolomics might give an answer. Complex phenotypes, which are determined by the influence of multiple small effect alleles, are an example of these situations. Consequently, the interest in metabolomics has increased ex-ponentially in the last years. As a newer discipline, the metabolomics bioinformatics analysis pipe-lines are not as standardized as in the other omic approaches. In this review we synthesized the different steps that need to be carried out to obtain biological insight from the annotated metabolite abundance raw data. These steps were grouped in three different modules: preprocessing, statistical analyses and metabolic pathway enrichment. We included within each one of them the different state-of-the art procedures and tools that can be used depending on the characteristics of the study, providing details about the characteristics of each method has, as well as the issues the reader might encounter. Finally, we introduce genome scale metabolic modeling as a tool for obtaining pseu-do-metabolomic data in situations where its acquisition is difficult, being possible to analyze the resulting data with the modules of the described workflow.

Background: Advancement of next-generation sequencing (NGS) technologies provides opportunities for large-scale Pharmacogenetic (PGx) studies and pre-emptive PGx testing to cover a wide range of genotypes present in diverse populations. However, NGS-based PGx testing is limited by the lack of comprehensive computational tools to support genetic data analysis and clinical decisions. Methods: Bioinformatics utilities specialized for human genomics and the latest cloud-based technologies are used for developing a bioinformatics pipeline for analyzing the genomic sequence data and reporting PGx genotypes. A database was created and integrated in the pipeline for filtering the actionable PGx variants and clinical interpretations. Strict quality verification procedures were conducted on variant calls with whole genome sequencing (WGS) dataset of the 1000 Genomes Project (G1K). The accuracy of PGx allele identification was validated using the whole genome sequencing dataset of the Pharmacogenetics Reference Materials from the Centers for Disease Control and Prevention (CDC). Results: The newly created bioinformatics pipeline, Pgxtools, can analyze genomic sequence data, identify actionable variants in 13 PGx relevant genes, and generate reports annotated with specific interpretations and recommendations based on clinical practice guidelines. Verified with two independent methods, we have found that Pgxtools consistently identifies variants more accurately than the results in the G1K dataset on GRCh37 and GRCh38. Conclusions: Pgxtools provides an integrated workflow for large-scale genomic data analysis and PGx clinical decision support. Implemented with cloud-native technologies, it is highly portable in a wide variety of environments from a single laptop to high performance computing clusters and cloud platforms for different production scales and requirements.

Metal-catalyzed lipid oxidation is a major factor in food waste, as it reduces shelf life. Addressing this issue, our study investigates the potential of hydrolysates derived from potato protein, a by-product of potato starch production, as metal chelating antioxidants. Through sequential enzymatic hydrolysis using Alcalase or Trypsin combined with Flavourzyme, we produced various hydrolysates, which were then fractionated via ultrafiltration. Using a combination of peptidomics and bioinformatics, we predicted the presence of metal-chelating and free radical scavenging peptides across all hydrolysate fractions, with a trend indicating a higher content of antioxidant peptides in lower molecular weight fractions. To validate these predictions, we utilized Surface Plasmon Resonance (SPR) and a 9-day emulsion storage experiment. While SPR demonstrated potential in identifying antioxidant activity, it faced challenges in differentiating between hydrolysate fractions due to significant standard errors. In the storage experiment, all hydrolysates showed lipid oxidation inhibition, though not as effectively as EDTA. Remarkably, one fraction (AF13) was not significantly different (p < 0.05) from EDTA in suppressing hexanal formation. These results highlight SPR and peptidomics/bioinformatics as promising yet limited methods for antioxidant screening. Importantly, this study reveals the potential of potato protein hydrolysates as antioxidants in food products, warranting further research.

: F-box proteins are a large gene family in plants, and play crucial roles in plant growth, development, and stress response. However, a comprehensive investigation of F-box family genes in cultivated peanut (Arachis hypogaea L.), has not been undertaken. In this study, we used a bioinformatics approach to identify and analyze F-box protein family genes in the cultivated peanut genome, and examined their expression profiles. Herein, a total of 95 F-box protein family members on 18 chromosomes were identified in cultivated peanut, which were sequentially designed as AhFBX1-AhFBX95. Phylogenetic analysis classified these AhFBXs into four different groups (Group I - IV), and eight subfamilies based on the structural domains of the C-terminal. Structural analysis revealed the presence of ten unique motifs within these genes. Additionally, the investigation of non-synonymous and synonymous substitution rates (Ka/Ks), and synteny in the F-box proteins, provided valuable insights into the evolutionary and domestication processes of AhFBX genes in peanut. Furthermore, in silico expression profiling of 22 different tissues demonstrated that most genes were predominantly expressed in leaves, pistils, and the aerial gynophore tip, while some members showed tissue-specific expression at the apex of both vegetative and reproductive shoots. Subsequently, we validated the expression of 9 differentially expressed genes (DEGs) using transcriptomeic data of lateral branch development and qRT-PCR analysis. These DEGs may be involved in regulating the development of lateral branches and subsequently affect the architecture of peanut plants.

Microbiomes, comprised of diverse microbial species and viruses, play pivotal roles in human health, environmental processes, and biotechnological applications and interact with themselves, their environment, and hosts via metabolites and signaling molecules. Our understanding of microbiomes is still limited and hampered by their complexity. A concept improving this understanding is systems biology, which focuses on the holistic description of biological systems utilizing experimental and computational methods. An important set of such experimental methods are metaomics methods which analyze microbiomes and output lists of molecular features. These lists of data are integrated, interpreted, and compiled into computational microbiome models, to predict, optimize, and control microbiome behavior. There exists a gap in understanding between microbiologists and modelers/bioinformaticians, stemming from a lack of interdisciplinary knowledge. This knowledge gap hinders the establishment of computational models in microbiome analysis. This review aims to bridge this gap and is tailored for microbiologists, researchers new to microbiome modeling, and bioinformaticians. To achieve this goal, it provides an interdisciplinary overview of microbiome modeling, starting with fundamental knowledge of microbiomes, metagenomics methods, and modeling formalisms. Furthermore, the review explains model building, examples of microbiome model applications for prediction, optimization, and control. It concludes with guidelines, software, and repositories for modeling. Each section provides entry-level information, serving as a valuable resource for comprehending and navigating the complex landscape of microbiome research and modeling.

DNA-binding transcription factors are broadly characterized as proteins that bind to specific sequences within genomic DNA and modulate the expression of downstream genes. This study focused on KojR—a transcription factor involved in the metabolism of kojic acid, which is an organic acid synthesized in Aspergillus oryzae and is known for its tyrosinase-inhibitory properties. However, the regulatory mechanism underlying KojR-mediated kojic acid synthesis in remains unclear. Hence, we aimed to obtain a comprehensive identification of KojR-associated genes using genomic Systematic Evolution of Ligands by Exponential enrichment (gSELEX-Seq) and RNA-Seq. During the genome-wide exploration of KojR-binding sites via gSELEX-Seq and identification of KojR-dependent differentially expressed genes (DEGs) using RNA-Seq, we confirmed that KojR directly regulates kojT, as was previously reported. We also observed that kojA expression, which may be controlled by KojR, was significantly reduced in a ΔkojR strain. Notably, no binding of KojR to the kojA promoter region was detected. Furthermore, certain KojR-dependent DEGs identified in the present study were associated with enzymes implicated in the carbon metabolic pathway of A. oryzae. This strongly indicates that KojR plays a central role in carbon metabolism in A. oryzae.

Many challenges in proteomics result from the high-throughput nature of the experiments. The paper first presents pre-analytical problems, which still occur, although the call for standardisation in omics has been ongoing for many years. This article also discusses aspects, which affect bioinformatic analysis, based on three sets of reference data measured with different orbitrap instruments. Despite continuous advances in mass spectrometer technology as well as analysis software, data set-wise quality control is still necessary and decoy-based estimation, although challenged by modern instruments, should be utilized. We draw attention to the fact that numerous young researchers perceive proteomics as a mature, readily applicable technology. However, it is important to emphasize that the maximum potential of the technology can only be realized by an educated handling of its limitations.

Аnalysis of small RNA sequencing data across a range of biofluids is a significant research area, given the diversity of RNA types that holds potential diagnostic, prognostic, and predictive value. The intricate task of segregating the complex mixture of small RNAs from both human and other species, including bacteria, fungi, and viruses, poses one of the most formidable challenges in the analysis of small RNA sequencing data, currently lacking satisfactory solutions. This study introduces sRNAflow, a user-friendly bioinformatic tool with a web interface designed for the analysis of small RNAs obtained from biological fluids. Tailored to the unique requirements of such samples, the proposed pipeline addresses various challenges, including filtering potential RNAs from reagents and environment, classifying small RNA types, managing small RNA annotation overlap, conducting differential expression assays, analysing isomiRs, and presenting an approach to identify the sources of small RNAs within samples. sRNAflow also encompasses an alternative alignment-free analysis of RNA-seq data, featuring clustering and initial RNA source identification using BLAST. This comprehensive approach facilitates meaningful comparisons of results between different analytical methods. The source code can be accessed at https://github.com/zajakin/sRNAflow under the GPL3 licence.

Universal Stress Protein (USP) plays an important regulatory role in the response to abiotic stress. Most of the research related to USP so far has been conducted on plant models such as Arabidopsis, rice, and cotton. The potato is one of the four major food crops in the world. Therefore, it is of great significance to study the expression pattern of the potato StUSP gene family under abiotic stress conditions. In this study, a total of 108 USP genes were identified from the genome of the Atlantic potato, divided into four subgroups. And based on the genetic structure, the physical and chemical properties of proteins and other aspects of their biological characteristics are comprehensively analyzed. Collinear analysis showed that the homologous genes of StUSPs and four other representative species (Arabidopsis, rice, coyote tobacco, and tomato) were highly conserved. The cis-regulatory elements of the StUSPs promoter are involved in plant hormones, environmental stress, mechanical damage, and light response. RNA-seq analysis showed that under different abiotic stresses, there are differences in the expression patterns of members of each subgroup. WGCNA co-expression analysis of the central gene showed that the differential co-expression gene is mainly involved in the plant-pathogen response process, plant hormone signal transduction, and the biosynthesis process of secondary metabolites. Through qRT-PCR analysis, it was confirmed that StUSP13, StUSP14, StUSP15, and StUSP41 may be important candidate genes involved in the response to adversity stress in potatoes. The results of this study provide the basis for further research on the functional analysis of StUSPs in the potato response to adversity stress.

Proteins in the crowded environment of human cells have often been studied regarding nonspecific interactions, misfolding, and aggregation, which may cause cellular malfunction and disease. Specifically, proteins with high abundance are more susceptible to these issues due to the law of mass action. Therefore, the surfaces of highly abundant cytoplasmic (HAC) proteins directly exposed to the environment can exhibit specific physicochemical, structural, and geometrical characteristics that reduce nonspecific interactions and adapt to the environment. However, the quantitative relationships between the overall surface descriptors still need clarification. Here, we used machine learning to identify HAC proteins using hydrophobicity, charge, roughness, secondary structures, and B-factor from the protein surfaces and quantify the contribution of each descriptor. First, several supervised learning algorithms were compared to solve binary classification problems for the surfaces of HAC and extracellular proteins. Then, logistic regression was adopted for the feature importance analysis of descriptors considering model performance (80.2% accuracy and 87.6% AUC) and interpretability. The HAC proteins showed positive correlations with negatively and positively charged areas but negative correlations with hydrophobicity, B-factor, beta structure proportion, roughness, and disordered regions proportion. Finally, the details of each descriptor could be explained concerning adaptative surface strategies of HAC proteins to regulate nonspecific interactions, protein folding, flexibility, stability, and adsorption. This study presented a novel approach using various surface descriptors to identify HAC proteins and provided quantitative design rules for the surfaces well-suited to human cellular crowded environments.

The historical and cultural use of food as a remedy for diseases is evident in societies, such as using Hibiscus sabdariffa to address conditions like hypertension and high blood glucose. The natural biocompounds in this plant, including Delphinidin-3-Sambubioside (DS3), Quercetin (QRC), and Hibiscus Acid (HA), have been linked to various health benefits. Despite individual attention, molecular targets for these compounds remain unclear. In this study, in-silico analysis employed bioinformatic tools; including Swiss Target Prediction, ShinnyGo 0.77, KEGG, and Stringdb, to identify molecular targets, pathways, and hub genes. A PubMed literature search complemented the results. DS3 demonstrated potential modifications in genes related to nitrogen and glucose metabolism, inflammation, angiogenesis, and cell proliferation, particularly affecting the PI3K-AKT pathway. QRC displayed interconnected targets across multiple pathways, with overlap with DS3 and a focus on cancer-related pathways. HA showed distinct targets, particularly associated with nervous system-related pathways. These findings highlight the need for targeted research on the molecular effects of DS3, QRC, and HA, providing valuable insights into potential therapeutic pathways.

The receptor tyrosine kinase (EphB4), which produces erythropoietin, was identified as a molecular target for cancer. We identified six peptides as potential inhibitors of EphB4 using an integrated approach combining online bioinformatics tools. According to PeptideRanker, CSM-peptides, SCMTHP, and AcPEP analysis, all peptides are effective against cancer, with the exception of LTAPGQATLPT and VEQNLVPGLK. Furthermore, Toxin-Pred predicts that none of the peptides are toxic, and that only VQLPGWRVFP and TMKGPPDTLQ are potentially allergenic using AllerTop server. PepSite2, HPEPDOCK, and MDockPeP servers were utilized to perform molecular docking, and PacDock was used to examine molecular interactions. The docking score (kcal/mol) range of six promising anticancer peptides against EphB4 was recorded as follows: VQLPGWRVFP (-240.72), FTKDDEWSCFPF (-201.74), SYLPPLSAEVTAK (-177.52), LTAPGQATLPT (-162.23), VEQNLVPGLK (-157.30), TMKGPPDTLQ (-149.67). The main bindings on the allosteric site to VQLPGWRVFP were with R706, S815, R819 and P843 on EphB4. Collectively, our results demonstrated that M. oleifera peptides ameliorate cancer by inhibiting EphB4, and they are promising candidates for the prevention and treatment of cancer.

Introduction: The distribution of typical bacterial redox enzymes such as pyruvate-ferredoxin oxidoreductase (PFOR) in protozoa remains interestingly puzzling. Previous studies have demonstrated diverse cellular localizations of PFOR in amitochondriate anaerobic protozoa such as Trichomonas vaginalis. PFOR is of particular pharmacological importance because it catalyzes the reductive bio-activation of nitro-based prodrugs to cytotoxic radical metabolites. Metronidazole was developed primarily as an antiprotozoal agent against infections caused by T. vaginalis. However, its antimicrobial spectrum was subsequently expanded to cover anaerobic bacterial infections. It has been shown that mutations in the genes encoding PFOR result in inherent resistance of T. vaginalis and anaerobic bacteria to metronidazole. Aims: To analyze the phylogenetic distribution of PFOR in selected protozoa and bacteria using proteins encoded by housekeeping genes as controls. Using comparative bioinformatics to test the hypothesis that PFOR was most likely acquired through horizontal gene transfer from bacteria, proteome-wide analysis and gene expression analysis to identify other genes that were putatively acquired through horizontal gene transfer from bacteria to protozoa. Methods: Sequence similarity queries were performed using the proteins of interest against the NCBI non-redundant protein sequence database with BLASTP version 2.12.0+ from the NCBI website. Multiple sequence alignments were performed with MEGA version X software using the Muscle algorithm. These were then exported in Mega format to construct phylogenetic trees using the neighbor-joining algorithm and the Poisson substitution model. Complete proteomes of representative protozoan and bacterial species were downloaded in fasta format from the ensemblgenomes database. HMM profile libraries were used to screen the proteomes with hmmscan of the HMMer package. Whole proteome BLAST searches were performed between protozoan species of interest against closely related protozoa. The results were compared to the proteome BLAST against the intestinal bacteria Desulfovibrio vulgaris. Gene enrichment analysis was performed between the exclusively present proteins in the protozoa and D. vulgaris. Only results with p-values <0.05 were considered for analysis. Results and conclusions: A plausible explanation for the restricted occurrence of PFOR in protozoa is based on the hypothesis that bacteria serve as potential sources of genes that enhance optimal adaptation of protozoa in hostile environments. The expanded cladograms of Entamoeba and Cryptosporidium with their closely related genera substantiated this hypothesis. The exclusively expressed proteins obtained from E. histolytica and the putative bacterial gene donor, D. vulgaris, showed an over-representation of eleven genes involved in small molecule metabolism, generation of precursor metabolites, and carbohydrate metabolism. If these results can be reproduced in other PFOR-possessing protozoa, it would provide a more validated evidence to support the horizontal transfer of pfor from bacteria. Subsequent syntenic analyses of the significantly enriched genes would be required to provide further information regarding the positional relatedness of these genes at the chromosomal level. Since metronidazole is an established and well tolerated drug for treating infections caused by PFOR-possessing pathogens, it can be considered as a potential drug candidate for the treatment of infections caused by Cryptosporidium, Spironucleus and Blastocystis.

In the modern era of genomic research, the scientific community is witnessing an explosive growth in the volume of published findings. While this abundance of data offers invaluable insights, it also places a pressing responsibility on genetic professionals and researchers to stay informed about the latest findings and their clinical significance. Genomic variant interpretation is currently facing a challenge in identifying the most up-to-date and relevant scientific papers, while also extracting meaningful information to accelerate the process from clinical assessment to reporting. Computer-aided literature search and summarization can play a pivotal role in this context. By synthesizing complex genomic findings into concise, interpretable summaries, this approach facilitates the translation of extensive genomic datasets into clinically relevant insights. To bridge this gap, we present VarChat (varchat.engenome.com), an innovative tool based on generative AI, developed to find and summarize the fragmented scientific literature associated with genomic variants into brief yet informative texts.VarChat provides users with a concise description of specific genetic variants, detailing their impact on related proteins and possible effects on human health. Additionally, VarChat offers direct links to related scientific trustable sources, and encourages deeper research. Availability: VarChat is freely available at varchat.engenome.com

Neisseria meningitidis (N. meningitidis) serogroup B (MenB) is the leading cause of invasive meningococcal disease worldwide. The pathogen has a wide range of virulence factors, which are potential vaccine components. Studying the genetic variability of antigens within a population, especially their long-term persistence, is necessary to develop new vaccines and predict the effectiveness of existing ones. The multicomponent 4CMenB vaccine (Bexsero), used since 2014, contains three major genome-derived recombinant proteins: factor H-binding protein (fHbp), Neisserial Heparin-Binding Antigen (NHBA) and Neisserial adhesin A (NadA). Here, we assessed the prevalence and sequence variations of these vaccine antigens in a panel of 5667 meningococcal isolates collected worldwide over the past 10 years and deposited in the PubMLST database. Using multiple amino acid sequence alignments and machine learning methods, we estimated the potential strain coverage by fHbp and NHBA vaccine variants (51 and about 30%, respectively); the NadA antigen sequence was found in only 17.8% of MenB genomes analyzed, but cross-reactive variants were present in less than 1% of isolates. Based on our findings, we proposed various strategies to improve the 4CMenB vaccine and broaden coverage of N. meningitidis strains.

Oxidative stress is increasingly implicated in male secondary infertility, where the seminal plasma altered proteome composition plays a critical role. This study aims to explore proteomic changes in the seminal plasma of fertile individuals with oxidative stress and secondary infertility patients, hypothesizing oxidative stress as a key factor in secondary infertility. Pooled semen samples from both groups underwent quantitative proteomics analysis using advanced mass spectrometry, with subsequent bioinformatic analysis using tools like DAVID, STRING, and IPA for identifying differentially expressed proteins (DEPs). The study identified 377 DEPs in the secondary infertile group and 523 DEPs in healthy controls, revealing seven mutual pathways including acute phase response signaling and free radical scavenging. Significant alterations were noted in proteins like C3 and SERPINA3, associated with sperm function and fertility. The findings suggest that oxidative stress-induced proteomic alterations in seminal plasma are linked to secondary infertility, highlighting potential biomarkers and therapeutic targets. This research underscores the importance of understanding the interplay between oxidative stress and seminal plasma proteins, potentially guiding personalized treatment strategies to improve fertility outcomes in affected individuals.

Transposons are mobile DNA sequences that contribute large fractions of many plant genomes. They provide exclusive resources for tracking gene and genome evolution and for developing molecular tools for basic and applied research. Despite extensive efforts, it is still challenging to accurately annotate transposons, especially for the beginners as transposon prediction requires necessary expertise in both transposon biology and bioinformatics. Moreover, the complexity of plant genomes and the dynamic evolution of transposons also bring difficulties for genome-wide transposon discovery. This review summarizes the three major strategies for transposon detections including repeat-based, structure-based, and homology-based annotation, and introduces the transposon superfamilies identified in plants thus far, and some related bioinformatics resources for detecting plant transposons. Furthermore, it describes the transposon classification and explains why the terms of ‘autonomous’ and ‘non-autonomous’ cannot be used to classify the superfamilies of transposons. Last, this review also discusses how to identify misannotated transposons and improve the quality of transposon database. This review provides helpful information about plant transposons and a beginner’s guide on annotating these repetitive sequences.

Cytochrome P450 oxidoreductase (POR) is an essential redox partner for steroid and drug-metabolizing cytochromes P450 located in the endoplasmic reticulum. Mutations in POR lead to metabolic disorders, including congenital adrenal hyperplasia, and affect the metabolism of steroids, drugs, and xenobiotics. In this study, we examined approximately 450 missense variants of the POR gene listed in the Genome Aggregation Database (gnomAD) using eleven different in silico prediction tools. We found that 64 novel variants were consistently predicted to be disease-causing by most tools. To validate our findings, we conducted population analysis and selected two variations in POR for further investigation. The human POR wild type, along with the R268W and L577P variants, were expressed in bacteria and subjected to enzyme kinetic assays using model substrate. We also examined the activities of several cytochrome P450 proteins in the presence of POR (WT or variants) by combining P450 and reductase proteins in liposomes. We observed a decrease in enzymatic activities (ranging from 35% to 85%) of key drug metabolizing enzymes, supported by POR variants R288W and L577P compared to WT-POR. These results validate our approach of curating a vast amount of data from genome projects and provide an updated and reliable reference for diagnosing POR deficiency.

Trehalose is a non-reducing disaccharide and participates in physiological activities such as or-gan formation, energy metabolism, and stress resistance of insects. In the present study, phylo-genetic analysis divided 21 Bmtret1 orthologs into three clades. These genes are equally distrib-uted on the nine chromosomes. Cis-elements in the promoter regions of Bmtret1s indicated the possible function of Bmtret1s in response to hormones and environmental stimulus. The qPCR analysis showed the significantly different expression levels of Bmtret1s in different tissues and organs, indicating possible functional divergence. In addition, most Bmtret1s showed disturbed expression levels in response to BmNPV stresses. Our results provide a foundation for further functional dissection of Tret1s in Bombyx mori and implicate them as potential regulators for an-tiviral responses.

The increased cultivation of Cannabis sativa L. in North America, represented by high Δ9-tetrahydrocannabinol-containing (high-THC) cannabis genotypes and low THC-containing hemp genotypes, has been impacted by an increasing number of plant pathogens. These include fungi which destroy roots, stems and leaves, in some cases causing a build-up of populations and mycotoxins in the inflorescences that can negatively impact quality. Viroids and viruses have also increased in prevalence and severity and can reduce plant growth and product quality. Rapid diagnosis of the occurrence and spread of these pathogens is critical. Techniques in the area of molecular diagnostics have been applied to study these pathogens in both cannabis and hemp. These include polymerase chain reaction (PCR)-based technologies, including RT-PCR, multiplex RT-PCR, RT-qPCR, and ddPCR, as well as whole genome sequencing (NGS) and bioinformatics. In this study, examples of how these technologies have enhanced the rapidity and sensitivity of pathogen diagnosis on cannabis and hemp will be illustrated. These molecular tools have also enabled studies on the diversity and origins of specific pathogens, specifically viruses and viroids, and these will be illustrated. Comparative studies on the genomics and metabolomics of healthy and diseased plants are urgently needed to provide insight into their impact on quality and composition of cannabis and hemp-derived products. Management of these pathogens will require monitoring of their spread and survival using the appropriate technologies to allow accurate detection, followed by appropriate implementation of disease control measures.

The integration of bioinformatics in drug discovery has revolutionized the field of pharmaceutical research. The use of computational tools and techniques has enabled researchers to analyze vast amounts of data, identify potential drug targets, and design new drugs with greater precision and efficiency. The integration of bioinformatics has also facilitated the development of personalized medicine, where drugs can be tailored to individual patients based on their genetic makeup. However, there are still challenges that need to be addressed, such as the need for more accurate predictive models and the ethical considerations surrounding the use of patient data. Overall, the integration of bioinformatics in drug discovery holds great promise for improving human health and advancing our understanding of disease mechanisms. In this mini-review, we discuss how Bioinformatics plays a crucial role in each step of drug discovery by providing tools and techniques to analyze large amounts of data generated from various sources, as well as the challenges and the opportunities offered by bioinformatics.

CagY is the largest and most complex protein from Helicobacter pylori's type IV secretion system (T4SS) and may participate in the modulation of gastric tissue inflammation. A three-dimensional structure has been reported for only two segments of the protein. To build a more complete model, particularly the region that spans between the outer membrane (OM) and the inner membrane (IM), we employed different approaches, including homology modeling, ab initio, and deep learning techniques. For the long-middle repeat region (MRR), modeling was performed using deep learning techniques and Molecular Dynamics. The modeled segments were assembled into a chain of 1595 aa, and a 14-chain CagY multimer structure was composed by structural alignment. The final multimer structure correlated with previously published struc-tures and allows to show how the multimer may form the T4SS channel through which CagA and other molecules are translocated to gastric epithelial cells. The model further confirmed that MRR, the most polymorphic and complex region of CagY, presents numerous cysteine residues forming disulfide bonds that stabilize the protein and suggest this domain probably functions as a contractile region that may play an essential role in the modulatory activity of CagY on tissue inflammation.

Gabapentin (GBP), a GABA analogue, is primarily used as an anticonvulsant to treat partial seizures and neuropathic pain. While a majority of the side effects are associated with the nervous system, emerging evidence suggests a high risk for heart diseases in the patients taking GBP. In the present study, we used a preclinical model of rats to investigate (1) the acute cardiovascular responses to GBP (bolus i.v. injection, 50 mg/kg) and (2) the effects of chronic GBP treatment (i.p. 100 mg/kg/day x 7 days) on cardiovascular function and the myocardial proteome. Under isoflurane-anesthesia, rat blood pressure (BP), heart rate (HR), and left ventricular (LV) hemodynamics were measured using Millar pressure transducers. The LV myocardium and brain cortex were analyzed by proteomics, bioinformatics, and western blot to explore the molecular mechanisms underlying GBP-induced cardiac dysfunction. In experiment (1), we found that i.v. GBP significantly decreased BP, HR, maximal LV pressure, and maximal and minimal dP/dt, whereas it increased IRP-AdP/dt, Tau, systolic, diastolic, and cycle durations (*p < 0.05 and **p < 0.01 vs baseline; n = 4/group). In experiment (2), we found that chronic GBP treatment resulted in hypotension, bradycardia, and LV systolic dysfunction, with no change in plasma norepinephrine. In the myocardium, we identified 109 differentially expressed proteins involved in calcium pathways, cholesterol metabolism, and galactose metabolism. Particularly, we found that calmodulin, a key protein of intracellular calcium signaling, was significantly upregulated by GBP in the heart but not in the brain. In summary, we found that acute and chronic GBP treatments suppressed cardiovascular function in rats, which is attributed to abnormal calcium signaling in cardiomyocytes. These data reveal a novel side effect of GBP independent of the nervous system, providing important translational evidence to suggest that GBP can evoke adverse cardiovascular events by depression of myocardial function.

Background: Myasthenia Gravis (MG) is a rare autoimmune disease presenting with au-to-antibodies that affect the neuromuscular junction. Aside from symptomatic treatment options, novel therapeutics include monoclonal antibodies (mAbs). IMGT®, the international ImMunoGe-neTics information system® (https://www.imgt.org), extends the characterization of therapeutic antibodies with a systematic description of their mechanisms of action (MOA) and makes them available through its database for mAbs and fusion proteins, IMGT/mAb-DB. Methods: Using the available literature data combined with the amino acid sequence analyses from mAbs managed in IMGT/2Dstructure-DB, the IMGT® protein database, biocuration allowed to define in a standardized way descriptions of MOA of mAbs that target molecules towards MG treatment. Results: New therapeutic targets include FcRn and molecules such as CD38, CD40, CD19, MS4A1 and interleukin-6 receptor. A standardized graphical representation of the MOA of selected mAbs was created and integrated within IMGT/mAb-DB. The main mechanisms involved in these mAbs are either blocking or neutralizing. Therapies directed to B cell depletion and plasma cells have a blocking MOA with an immunosuppressant effect along with Fc-effector function (MS4A1, CD38) or FcγRIIb engager effect (CD19). Monoclonal antibodies targeting the complement also have blocking MOA with a complement inhibitor effect and treatments targeting T cells have a blocking MOA with an immunosuppressant effect (CD40) and Fc-effector function (IL6R). On the other hand, FcRn antagonists present a neutralizing MOA with an FcRn inhibitor effect. Conclusion: The MOA of each new mAb needs to be considered in association with the immuno-pathogenesis of each of the subtypes of MG in order to integrate the new mAbs as a viable and safe option in the therapy decision process. In IMGT/mAb-DB, mAbs for MG are characterized by their sequence, domains, chains and their MOA is described.

The role of desmoglein-3 (DSG3) in oncogenesis is unclear. This study aimed to uncover molecular mechanisms through comparative transcriptome analysis in oral cancer cells, defining potential key genes and associated biological processes related to DSG3 expression. Four mRNA libraries of oral squamous carcinoma H413 cell lines were sequenced and 599 candidate genes exhibited differential expression between DSG3-overexpressing and matched control lines, with 12 genes highly significantly differentially expressed, including 9 upregulated and 3 downregulated. Genes with known implications in cancer, such as MMP-13, KRT84, OLFM4, GJA1, AMOT, and ADAMTS1, were strongly linked to DSG3 overexpression. Gene ontology analysis indicated that the DSG3-associated candidate gene products participated in crucial cellular processes such as junction assembly, focal adhesion, extracellular matrix formation, intermediate filament organization, and keratinocyte differentiation. Validation of RNA-Seq was performed through qRT-PCR, Western blotting, and immunofluorescence analyses. Furthermore, transmission electron microscopy meticulously examined desmosome morphology, revealing slightly immature desmosome structure in DSG3-overexpressing cells compared to controls. No changes in desmosome frequency and diameter were observed between the two conditions. This study underscores intricate and multifaceted alterations associated with DSG3 in oral squamous carcinoma cells, implying a potential oncogenic role of this gene in biological processes that enable cell communication, motility and survival.

The term new genomic techniques (NGTs), is an umbrella term to describe a variety of techniques that can alter the genetic material of an organism and that have emerged or have been developed since 2001, when the existing genetically modified organism (GMO) legislation was adopted. The analytical framework to detect GMOs in Europe is a established single harmonized procedure which is mandatory for authorization of GM food and feed, thus generating reliable, transparent and effective labelling scheme for GMO products. However, NGT products challenge the implementation and enforcement of the current regulatory system in the EU, relating in particular to the detection of NGT products that contain no foreign genetic material. Consequently, current detection methods will fail to meet minimum performance requirements. Although existing detection methods may be able to detect and quantify even small alterations in the genome, this does not necessarily confirm the distinction between products resulting from NGTs subject to the GMO legislation and other products. Therefore, this study provides a stepwise approach for the in silico prediction of PCR systems specificity by testing a bioinformatics pipeline for amplicon and primer set searches in current genomic databases. We have tested our methodology in two mutant genotypes produced by CRISPR/Cas9 in Arabidopsis thaliana. Overall, our results demonstrated that the single PCR system developed for the identification of a nucleotide insertion in the grf1-3 genotype has sufficient hits in the databases that would not allow discrimination of this mutanted event. Differently, the second mutated genotype grf8-61 which contains a -3bp deletion, found no matches in the sequence variant database. Our approach is the first step in decision making on PCR system and modification of them for empirical testings.

With an ever increasing amount of research data available, it becomes constantly more important to possess data literacy skills to benefit from this valuable resource. An integrative course was developed to teach students the fundamentals of data literacy through an engaging genome sequencing project. Each cohort of students performed planning of the experiment, DNA extraction, nanopore sequencing, genome sequence assembly, prediction of genes in the assembled sequence, and assignment of functional annotation terms to predicted genes. Students learned how to communicate science through writing a protocol in the form of a scientific paper, providing comments during a peer-review process, and presenting their findings as part of an international symposium. Many students enjoyed the opportunity to own a project and to work towards a meaningful objective.

In order to enhance understanding of RNA-binding proteins in rice, a comprehensive investigation was conducted on the RRM1 gene family of rice, encompassing genome-wide identification and exploration of its role in rice blast resistance. Physical and chemical properties of the OsRRM1 gene family in rice was analyzed, including conserved domain, motif, location information, gene structure, phylogenetic tree, collinearity analysis, cis-acting elements, GO, and KEGG. Furthermore, the expression patterns of the OsRRM1 gene were examined at different time intervals following rice blast treatment. Furthermore, the alterations in expression patterns of selected OsRRM1 genes were assessed using quantitative real-time PCR(qRT-PCR). A total of 212 members of the OsRRM1 gene family were identified, which were dispersed across 12 chromosomes. Many of these genes exhibit multiple exons and introns, all of which encompass the conserved RRM1 domain and share analogous motifs. This observation suggests a high degree of conservation within the encoded sequence domain of these genes. Phylogenetic analysis revealed the existence of five subfamilies within the OsRRM1 gene family. Furthermore, the investigation of the promoter region identified homeopathic elements that are involved in nucleic acid binding and interaction with multiple transcription factors. By employing GO and KEGG analysis, four RRM1 genes were tentatively identified as crucial contributors to plant immunity, while the RRM1 gene family was also found to have a significant involvement in the complex of alternative splicing. Additionally, gene expression analysis indicated that the majority of OsRRM1 genes exhibited constitutive expression. The results of the qRT-PCR analysis revealed distinct temporal changes in the expression pattern of the OsRRM1 gene following rice blast treatment. These findings contribute to the existing knowledge of the OsRRM1 gene family, establish a foundation for further investigation into the role of the OsRRM1 gene in response to rice blast infection, and hold theoretical significance for future studies on the functionality of the OsRRM1 gene.

Precision medicine approaches often relies on complex and integrative analysis of multiple biomarkers from “omics” data to generate insights that can help either diagnostics, prognostics or therapeutical decisions. Such insights are often made using Machine learning (ML) models that make sample classification for a particular phenotype (yes/no). Building such models is a challenge and time-consuming, requiring advanced coding skills and mathematical modelling expertise. Artificial intelligence (AI) is a methodological solution that has the potential to facilitate, optimize and scale model development. In this work, we developed an AI-based, user-friendly and code-free platform (https://digitalphenomics.com) that fully automates the development of predictive models from quantitative “omics” data. Here, we show the application of this tool with the development of cancer survival prognostics models using real-life data from breast, lung and renal cancer transcriptomes. We report and compare their sensitivity, specificity, accuracy and Receiver Operating Characteristic (ROC) curve Area Under the Curve (AUC). Further, we report the associated sets of genes (biomarkers) and their expression pattern that are predictive of cancer survival. Moreover, we made our models available as online tools to generate prognostic predictions based on the gene expression of the biomarkers. In conclusion, we demonstrated that our tool is a robust user-friendly solution to develop bespoke predictive tools from “omics” data which facilitate precision medicine introduction to the point-of-care.

Gastric cancer (GC) is a highly heterogeneous, complex disease and the fifth most common cancer worldwide (about one million cases and 784 000 deaths worldwide in 2018). CG is lately diagnosed and guarantees a poor prognosis for GC (the 5-year survival rate is less than 20%, but in early detection can reach 90%). This study evaluated the transcriptional profile in tumor gastric samples to find genes highly expressed during tumor establishment and use the related proteins as targets to find new anticancer molecules. Data was collected at Gene Expression Omnibus (GEO) bank to obtain 3 dataset matrices that analyze gastric tumor tissue versus normal gastric tissue, performed microarray using GPL570 platform, and from different sources. The genes found in silico analysis were confirmed in several lines of GC cells by RT-qPCR. The protein data bank was used to find the correspondent proteins. Then a structural-based virtual screening was done to find molecules, and docking analysis was done using the DockThor server. Our results of transcriptomic analysis, together with RT-qPCR, confirm the high expression of the genes AJUBA, FBXL13, CCDC69, CD80 and NOLC1 in GC lines. Based on that, the correspondent proteins were used in SBVS analysis. Five molecules, one each target, MCULE-2386589557-0-6, MCULE-7343047040-0-1, MCULE-5230409338-0-3, MCULE-9178344200-0-1 and MCULE -5881513100-0-29. All molecules have favorable pharmacokinetic, pharmacodynamic and toxicological properties. Molecular docking analysis revealed that the molecules interact with proteins in critical sites for their activity. Using a virtual screening approach, a molecular docking study was performed for proteins encoded by genes that play important roles in cellular functions for carcinogenesis. Combining a systematic collection of public microarray data with a comparative meta-profiling, RT-qPCR, SBVS, and molecular docking analysis provided a suitable approach to find genes involved in GC and work the correspondent protein to strive for new molecules with anticancer properties.

In order to provide a new direction for the diagnosis and treatment of gastric cancer, the differentially expressed blood-based miRNAs in gastric carcinoma were obtained and analyzed from the datasets GSE113486, GSE112264 and GSE113740 by the SangerBox, the miRNA - target genes prediction, GO and KEGG functional enrichment analyses with DAVID 6.8, and the STRING database with Cytoscape, respectively. The correlation between the key miRNAs expression and the prognosis and diagnostic value of gastric cancer patients were detected by survival prognosis and ROC curve. QRT-PCR was used to detect the possible key miRNAs expression in human blood species. The dataset GSE108307 which shows the miRNAs affected by H. pylori infection was analyzed by GEO2R. The expression level of CagA and the clinically used gastric cancer markers were detected in human blood species for contrast. Five key blood-based miRNAs were determined, including hsa-miR-124-3p, hsa-miR-125a-3p, hsa-miR-29b-3p, hsa-miR-4276, and hsa-miR-575.The detection in human blood species as well as the cross analysis with H. pylori infection and clinically used markers expression suggested the importance and efficacy of the specific miRNAs in the early diagnosis and monitoring of gastric carcinoma. The study identified potential novel blood miRNA biomarkers of gastric cancer by bioinformatics and human blood species detection, which may give a new possibility of the important biomarkers for the diagnosis and prognosis of gastric cancer.

Objective: To verify SPRR2A is a hub gene for endometrial cancer (EC). Methods: Bioinformatics analysis was performed on the high-throughput sequencing dataset of 587 EC cases in TCGA. RT-qPCR verified the differential expression of the SPRR2A gene in EC cell lines. Results: Enrichment of differentially expressed genes in comparison with normal tissues has been found in several biological process pathways, including keratinization and epidermal growth. The signalling pathways of differentially expressed genes were mainly focused on the receptor-ligand pathway and passive transmembrane transport pathway. High SPRR2A expression was a poor prognostic indicator for EC, according to a one-way Cox analysis of TCGA EC data. Compared with normal tissues, EC tissues had considerably elevated SPRR2A expression according to the differential expressed genetic analysis of TCGA-UCEC data (P &lt; 0.05). RT-qPCR was used to experimentally confirm that the Ishikawa EC cell line expressed SPRR2A mRNA at a greater level than the human endometrial epithelial cell line and that this difference was statistically obvious (P &lt; 0.01). Conclusion: SPRR2A has been identified as hub gene for EC and potentially become a new marker for early diagnosis, precise treatment and prognosis predictor of EC.

Background： Sepsis-induced acute lung injury (ALI) is characterized by disruption of the epithelial barrier and activation of alveolar macrophages (AMs), which leads to uncontrolled pulmonary inflammation. However, effective treatments for ALI are unavailable. This study aimed to discover potential diagnostic molecular biomarkers based on bioinformatics, which will benefit the diagnosis and treatment of sepsis-induced ALI. Methods: GSE10474 was analyzed for differentially expressed genes (DEGs) in sepsis patients with ALI (sepsis + ALI) compared with sepsis patients without ALI. Functional enrichment analysis and protein-protein interaction (PPI) network were performed. on the DEGs via R package “clusterProfiler”, and visualized via Cytoscape. Prediction analysis of microarrays (PAM) was performed to identify diagnostic biomarkers and the diagnostic ability of diagnostic biomarkers was accessed via receiver operating characteristic (ROC) curves. Moreover, interactions among diagnostic biomarkers were analyzed via GeneMANIA. We also analyzed the function of diagnostic biomarkers and predicted their corresponding drugs via Cytoscape plugin BiNGO and web tool DGIdb. At last, we analyzed the transcriptional regulation of the diagnostic biomarkers via the web tool miRNet. Results: 71 genes were found to be differentially expressed in the sepsis + ALI group, mainly involved in immune-related biological processes and pathways. STRING database indicated 31 DEGs have protein-protein interactions. In addition, the PAM identified 6 diagnostic biomarkers, including HIST1H4H, CDKN1A, HMOX1, NQO2, RHOB, and TREM1, from these 31 DEGs. Conclusion: In conclusion, through bioinformatics analyses, we identified 6 potential diagnostic biomarkers and targets for sepsis-induced ALI.

The field of mitochondrial genomics has advanced rapidly and revolutionized disciplines from molecular anthropology, population genetics, to medical/oncogenetics. However, mtDNA next-generation sequencing (NGS) analysis for matrilineal haplotyping and phylogeographic inference remains hindered by the lack of a consolidated, mitogenome database and efficient bioinformatics pipeline. To address this, we developed a customized human mitogenome database (hMITO DB) embedded in a CLC Genomics workflow for read mapping, variant analysis, haplotyping, and geo-mapping. The database was constructed from 4,286 mitogenomes. The macro-haplogroup (A to Z) distribution and representative phylogenetic tree were found consistent with published literature. The hMITO DB automated workflow was tested using mtDNA-NGS sequences derived from Pap smears and cervical cancer cell lines. The auto-generated read mapping, variants track, and table of haplotypes and geo-origins were completed in 15 min for 47 samples. The mtDNA workflow proved to be a rapid, efficient and accurate means of sequence analysis for translational mitogenomics.

It has been previously shown that aldehyde dehydrogenase (ALDH) family member ALDH1A1 has a significant association with acute myeloid leukemia (AML) patient risk group classification, and that AML cells lacking ALDH1A1 expression can be readily killed by chemotherapy. In the past, however, a redundancy between the activities of subgroup members of the ALDH family has hampered conclusive evidence for addressing the role of specific ALDH genes. We here describe the bioinformatics evaluation of all nineteen member genes of the ALDH family as prospective actionable targets for the development of methods aimed to improve AML treatment. We implicate ALDH1A1 in the development of recurrent AML, and we show that from the nineteen members of the ALDH family, ALDH1A1 and ALDH2 have the strongest association with AML patient risk group classification. Furthermore, we discover that the sum of the expression values for RNA from the genes ALDH1A1 and ALDH2 has a stronger association with AML patient risk group classification and survival than either one gene alone. In conclusion, we identify ALDH1A1 and ALDH2 as prospective actionable targets for the treatment of AML in high risk patients. Substances that inhibit both enzymatic activities can constitute potentially effective pharmaceutics.

Oral squamous cell carcinoma (OSCC) is one of the most prevalent human malignancies and a global health concern with a poor prognosis despite therapeutic advances, highlighting the need for a better understanding of its molecular background. The genomic landscape of OSCC is well-established and recent research has focused on miRNAs, which regulate gene expression and may be useful as non-invasive biomarkers. A plethora of findings regarding miRNA expression have been generated, posing challenges for their interpretation and identification of disease-specific molecules. Hence, we opted to identify the most important miRNA molecules by bridging genetics and epigenetics, focusing on the key genes implicated in OSCC development. Based on published reports, we have developed a custom panel of 15 major oncogenes and a second panel of 5 major tumor suppressor genes. Following a miRNA/target interaction analysis and a comprehensive study of the literature, we selected the miRNA molecules, which target the majority of each gene panel and are reported to be downregulated or upregulated in OSCC, respectively. As a result, miR-34a-5p, miR-155-5p, miR-124-3p, miR-1-3p and miR-16-5p appeared to be the most OSCC-specific. Their expressional patterns, their verified targets and the signaling pathways affected by their dysregulation in OSCC, are thoroughly discussed in this review.

The advancement of mass spectrometry technologies has revolutionised plant metabolomics research by enabling the acquisition of raw metabolomics data. However, the identification, analysis, and visualisation of these data require specialised tools. Existing solutions lack a dedicated plant-specific metabolite database and pose usability challenges. To address these limitations, we developed PlantMetSuite, a web-based tool for comprehensive metabolomics analysis and visualisation. PlantMetSuite encompasses interactive bioinformatics tools and databases specifically tailored for plant metabolomics data, facilitating upstream-to-downstream analysis in metabolomics and supporting integrative multi-omics investigations. PlantMetSuite can be accessed directly through a user's browser without the need for installation or programming skills. The tool is freely available at https://plantmetsuite.verygenome.com/ and will undergo regular updates and expansions to incorporate additional libraries and newly published metabolomics analysis methods. The tool's significance lies in empowering researchers with an accessible and customisable platform for unlocking plant metabolomics insights.

In the 21st century, the effects of HIV-associated neurocognitive disorders (HAND) have been significantly reduced in individuals due to the development of antiretroviral therapies (ARTs). However, the growing epidemic of polysubstance use (PSU) has led to concern for the effects of PSU on HIV-seropositive individuals. To effectively treat individuals affected by HAND, it is critical to understand the biological mechanisms affected by PSU including identification of novel markers. To fill this important knowledge gap, we used an in vivo HIV-Transgenic (HIV-Tg) animal model, to investigate the effects of the combined use of chronic methamphetamine (METH) and oxycodone (oxy). RNA-Seq analysis on the striatum - a brain region that is primarily targeted by both HIV and drugs of abuse - identified key differentially expressed markers post-METH and oxy exposure. Furthermore, ClueGO analysis and Ingenuity Pathway Analysis (IPA) revealed crucial molecular and biological functions associated with ATP-activated adenosine receptor, neuropeptide hormone activity, and the oxytocin signaling pathway to be altered between the different treatment groups. The current study further reveals the harmful effects of chronic PSU and HIV infection that can subsequently impact neurological outcomes in polysubstance users with HAND.

Hypertrophic Cardiomyopathy (HCM) is a common inherited disorder that can lead to heart failure and sudden cardiac death, characterized at the histological level by focal areas of myocyte disarray, hypertrophy and fibrosis, and few disease-targeted therapies exist. To identify, focal, spatially restricted alterations in transcriptional pathways and reveal novel therapeutic targets, we performed a spatial transcriptomic analysis of areas of focal myocyte disarray compared to areas of normal tissue, using a commercially available platform (GeoMx, nanoString). We analyzed surgical myectomy tissue from four patients with HCM and control interventricular septum tissue from two unused organ donor hearts that were free of cardiovascular disease. Histological sections were reviewed by an expert pathologist and 72 focal areas with varying degrees of myocyte disarray (normal, mild, moderate, severe) were chosen for analysis. Areas of interest were interrogated with the Human Cancer Transcriptome Atlas designed to profile 1800 transcripts. Differential expression analysis revealed significant changes in gene expression between HCM and Control tissue, and functional enrichment analysis indicated these genes were primarily involved in interferon production and mitochondrial energetics. Within HCM tissue, differentially expressed genes between areas of mild and moderate disarray were enriched for genes related to mitochondrial energetics (moderate disarray) and response to oxygen/cytokine levels (mild disarray). The comparison between areas of moderate and severe disarray were enriched for genes related to the c-Jun N-terminal kinase (JNK) cascade in severe disarray. Analysis of ligand-receptor pair gene expression revealed that HCM tissue exhibited downregulation of platelet-derived growth factor (PDGF), NOTCH, junctional adhesion molecule, and CD46 signaling, while showing upregulation of fibronectin, CD99, cadherin, and amyloid precursor protein signaling. A deconvolution analysis utilizing the matched single nuclei RNA-sequencing (snRNA-seq) data to determine cell type composition in areas of interest revealed significant differences in fibroblast and vascular cell composition in areas of severe disarray when compared to normal areas in HCM samples. Cell composition in normal areas from control tissue was also divergent from normal areas in HCM samples, which was consistent with the differential expression results. Overall, our data identify novel and potential disease-modifying targets for therapy in HCM.

Spermatogonial stem cells (SSCs) are the only primitive spermatogonial cells in males that can naturally transmit genetic information to their offspring and replicate throughout their lives. Phospholipase D family member 6 (PLD6) has recently been found to be a surface marker for SSCs in mice and boars; however, it has not been validated in cattle. The results of RT-PCR and qRT-PCR found that the relative expression of PLD6 gene in the testicular tissues of 2-year-old Simmental calves was significantly higher than that of 6-month-old calves. Immunofluorescent staining further verified the expression of PLD6 protein in bovine spermatogenic cells like germ cell marker VASA. Based on multiple bioinformatic databases, PLD6 is a conservative protein which has high homology with mouse Q5SWZ9 protein. It is closely involved in the normal functioning of the reproductive system. Molecular dynamics simulation analyzed the binding of PLD6 as a phospholipase to cardiolipin (CL), and PLD6-CL complex showed high stability. The protein interaction network analysis showed that there is a significant relationship between PLD6 and piRNA binding protein. PLD6 acts as an endonuclease and participates in piRNA production. In addition, PLD6 in bovine and mouse testes has a similar expression pattern with the spermatogonium-related genes VASA and PIWIL2. In conclusion, these analyses imply that PLD6 has a relatively high expression in bovine testes and could be used as a biomarker for spermatogenic cells including SSCs.

Inflammatory bowel disease (IBD), including Crohn’s Disease (CD) and Ulcerative Colitis (UC) are chronic multifactorial disorders which affect the gastrointestinal tract with variable extent. Despite extensive research, their etiology and exact pathogenesis are still unknown. Cell-free DNAs (cfDNAs) are defined as any DNA fragments which are free from origin cell and able to circulate into bloodstream with or without microvescicles. CfDNAs are now being increasingly studied in different human diseases, like cancer or inflammatory diseases. However, to date it is unclear how IBD etiology is linked to cfDNAs in plasma. Extrachromosomal circular DNA (eccDNA) are non-plasmidic, nuclear, circular and closed DNA molecules found in all eukaryotes tested. CfDNAs appear to play an important role in autoimmune diseases, inflammatory processes, and cancer; recently, interest has also grown in IBD, and their role in the pathogenesis of IBD has been suggested. We now suggest that eccDNAs also plays a role in IBD. In this review we have comprehensively collected available knowledge in literature regarding cfDNA, eccDNA, and structures involving them such as neutrophil extracellular traps and exosomes, and their role in IBD. Finally, we focused on old and novel potential molecular therapies and drug delivery systems, such as nanoparticles, for IBD treatment.

Chicken infectious anemia (CIA) poses a significant threat to the chicken industry in China. De-spite its non-specific symptoms, the disease is often overlooked. This study aimed to conduct a comprehensive analysis of the etiology and pathology of CIA in Guangxi Province, China. Three strains of the chicken infectious anemia virus (CIAV) were isolated from liver samples of diseased 20-week-old chickens. The complete genomes of these strains were sequenced, and experiments on specific pathogen-free (SPF) chicks revealed that the GX21121 strain exhibited high virulence. Histopathological examination of the deceased chicks showed liver cell necrosis, fibrous-serous exudation, inflammatory cell infiltration, hemorrhage in liver tissues, as well as congestion in lung and renal tissues. Phylogenetic analysis of the genome revealed that the three strains had a close genetic relationship to the Heilongjiang wild-type (GenBank KY486144). The genetic evolu-tion of their VP1 genes indicated that all three CIAV isolates belonged to genotype IIIc. In sum-mary, this study demonstrated the genomic diversity of three CIAV strains in adult layer hens. The isolation and characterization of the GX21121 strain as a highly virulent isolate provide val-uable information for further investigations into the etiology, molecular epidemiology, and viral evolution of CIAV.

With the exponential growth of machine learning and development of Artificial Neural Network (ANNs) in recent years, there is great opportunity to leverage this approach and accelarate bio-logical discoveries through applications on the analysis of bioinformatics data. Various types of datasets including for example protein or gene interaction networks, molecular structures and cellular signalling pathways, have already been used for machine learning by training ANNs for inference and pattern classification. However, unlike regular data structures that are commonly used in the computer science and engineering fields, bioinformatics datasets present challenges that require unique algorithmic approaches. The recent development of the geometric and deep learning approach within the machine learning field, is very promising towards accelerating analysis complex bioinformatics datasets. The principles of ANNs and their importance for bio-informatics machine learning is demonstrated herein, through presentation of the undelying mathematical and statistical foundations from group theory, symmetry, linear algebra. Further-more, the structure and functions of ANN algorithms that form the core principles of artificial intelligence are explained, in relation to the bioinformatics data domain. Overall, the manuscript provides guidance for researchers to understand the principles required for practicing machine learning and artificial intelligence, with the special considerations towards bioinformatics applications.

Motivation: The visualization of biological data is a fundamental technique that enables researchers to understand and explain biology. Some of these visualizations have become iconic, for instance: tree views for taxonomy, cartoon rendering of 3D protein structures, or tracks to represent features in a gene or protein, for instance in a genome browser. Nightingale provides visualizations in the context of proteins and protein features. Results: Nightingale is a library of re-usable data visualization web components that are currently used by UniProt and InterPro, among other projects. The components can be used to display protein sequence features, variants, interaction data, 3D structure, etc. These components are flexible, allowing users to easily view multiple data sources within the same context, as well as compose these components to create a customized view. Availability: Nightingale examples and documentation are freely available at https://ebi-webcomponents.github.io/nightingale/. It is distributed under the MIT license, and its source code can be found at https://github.com/ebi-webcomponents/nightingale.

Preeclampsia (PE) is a serious complication of pregnancy with a pathogenesis that is not fully understood, though it involves the impaired invasion of extravillous trophoblasts (EVTs) into the decidual layer during implantation. Recently, others have found that diverse cell types, including EVTs, produce exosomes filled with molecular cargo (mainly proteins and RNAs) that can be transported to other cells and organs both locally and long-range. The cargo delivered by exosomes can signal to and modify the receiving cells and their environment, and EVT-derived exosomes may influence the pathogenesis of PE. Because the risk of PE is actually decreased by cigarette smoking, we considered the possibility that nicotine, a critical component of tobacco smoke, might protect against PE by modifying the content of exosomes from EVTs. In this study, we applied nicotine stimulation to cultured EVTs and subjected their secreted exosomes to proteomic analysis. We identified many proteins whose abundance in exosomes was modified by nicotine treatment of the donor EVTs, and we used bioinformatic annotation and network analysis to select five key hub proteins with potential roles in the pathogenesis or prevention of PE.

Here, one hundred patients (50 smokers and 50 non-smokers) clinically diagnosed with COVID-19 were studied. Yet, bioinformatics was used to predict epitopes on Tobacco mosaic virus coat protein (TMV-CP) to produce antibodies towards SARS-CoV-2. Death was three times higher in non-smokers than in smokers. However, biochemical parameters did not separate the groups. Bioinformatics analysis predicted the presence of B-cell epitopes in TMV-CP, suggesting the production of antibodies anti-TMV-CP in smoker patients. Smokers may develop severe forms of COVID-19, but survival was superior in the evaluated group than in non-smokers. Anti-TMV-CP antibodies, potentially present in smokers, might act as a pro-immune agent against SARS-CoV-2 at earlier stages of infection. These data are helpful for future studies assessing COVID-19 in smokers.

Chemoinformatics is the combination of physical chemistry theory with computer and information science techniques "in silico techniques" to a variety of descriptive and prescriptive chemistry issues, including applications to biology, drug discovery, and related molecular areas. On the other hand, machine learning has become a vital tool for drug designers to mine chemical information from enormous compound databases in order to build medications with important biological features. In Egypt, chemoinformatics were used recently in different applications. For example; the development of new antimicrobial agents like tetracycline analog B (Iodocycline). Another area of application in bioinformatics included heart disease classification based on hybrid ensemble stacking technique. In the era of COVID-19, hybrid approach for COVID-19 detection from chest radiography was applied. Furthermore, bio/chemoinformatics methods were used to compare antimicrobial medications in order to choose an effective nose-to-brain formulation targeting meningitis infection by using differences in the drugs' primary structural, topological, and electronic characteristics. An example for this included cefotaxime and ceftriaxone that were examined at three levels: at formulation level, by comparing the loading in gelatin and tripalmitin matrices as basis for the production of nanoparticulate systems, at biopharmaceutical level, through interaction with mucin and the P-gp efflux pumps, and at therapeutic level, through studying the interaction with S. pneumoniae bacterial receptors.

In the last two decades many detailed full transcriptomic studies on complex biological samples have been published and included in large gene expression repositories. These studies primarily provide bulk expression signal for each sample, including multiple cell-types mixed within the global signal. The cellular heterogeneity in these mixtures does not allow the activity of specific genes in specific cell types to be identified. Therefore, inferring relative cellular composition is a very powerful tool to achieve a more accurate molecular profiling of complex biological samples. In recent decades, computational techniques have been developed to solve this problem by applying cell deconvolution methods, designed to decompose cell mixtures into their cellular components and calculate the relative proportions of these elements. Some of them, only calculate the cell proportions (supervised methods), while other deconvolution algorithms can also identify the gene signatures specific for each cell type (unsupervised methods). In these work, five deconvolution methods (CIBERSORT, FARDEEP, DECONICA, LINSEED and ABIS) were implemented and used to analyze blood and immune cells, and also cancer cells, in complex mixture samples (using three bulk expression datasets). Our study provides three analytical tools (corrplots, cell-signature plots and bar-mixture plots) that allow a thorough comparative analysis of the cell mixture data. The work indicates that CIBERSORT is a robust method optimized for the identification of immune cell-types, but not as efficient in the identification of cancer cells. We also observed that LINSEED is a very powerful unsupervised method that provides precise and specific gene signatures for each of the five main immune cell-types analyzed: neutrophils and monocytes (in the myeloid lineage), B-cells, NK cells and T-cells (in the lymphoid lineage).

Selectins are responsible for the early stages of cell migration as they control cell adhesion. They make the microenvironment permissive for metastatic events by promoting the activation of other cell adhesion molecules (CAMs). We employed several robust bioinformatics tools to evaluate gene sequence; single nucleotide polymorphisms (SNPs) in intronic and UTR regions; and missense SNPs with amino acid change in L-selectin (SELL), E-selectin (SELE), P-selectin (SELP) and PSGL-1 (SELPLG). We demonstrated that gene polymorphisms rs2229569, rs1131498, rs4987360, rs4987301 and rs2205849; polymorphisms rs3917777, rs2205894 and rs2205893 of SELP gene; rs7138370, rs7300972 and rs2228315 variants of SELPLG gene; and rs1534904 and rs5368 polymorphisms of SELE gene may produce important alterations in the DNA structure and consequent alterations in the morphology and function of the corresponding proteins. These Selectin polymorphic variants deserve further investigation in cancer patients, as they may become useful clinical markers for risk determination, diagnostic and prognostic biomarkers or even targets for targeted therapies.

Untargeted metabolomics is a powerful tool for measuring and understanding complex biological chemistries. However, employment, bioinformatics and downstream analysis of mass spectrometry (MS) data can be daunting for inexperienced users. Numerous open-source and free-to-use data processing and analysis tools exist for various untargeted MS approaches, including liquid chro-matography (LC), but choosing the ‘correct’ pipeline isn’t straight-forward. This tutorial, in con-junction with a user-friendly online guide presents a workflow for connecting these tools to process, analyse and annotate various untargeted MS datasets. The workflow is intended to guide explor-atory analysis in order to inform decision-making regarding costly and time-consuming down-stream targeted MS approaches. We provide practical advice concerning experimental design, organisation of data and downstream analysis, and offer details on sharing and storing valuable MS data for posterity. The workflow is editable and modular, allowing flexibility for updated/ changing methodologies and increased clarity and detail as user participation becomes more common. Hence, the authors welcome contributions and improvements to the workflow via the online repository. We believe that this workflow will streamline and condense complex mass-spectrometry approaches into easier, more manageable, analyses thereby generating opportunities for researchers previously discouraged by inaccessible and overly complicated software.

Lupinus mutabilis Sweet (Fabaceae), “tarwi” or “chocho”, is an important grain legume in the Andean region. In Peru, studies on tarwi have been mainly focused on morphological features, however, the have not been molecularly characterized. Currently, it is possible to explore genetic parameters of plants with reliable and modern methods like genotyping-by-sequencing (GBS). We here for the first time used single nucleotide polymorphisms (SNPs) markers to infer the genetic diversity and population structure of 89 accessions of tarwi from nine Andean regions of Peru. A total of 5922 SNPs distributed along all chromosomes of tarwi were identified. STRUCTURE analysis revealed that this crop is grouped into two clusters. A dendrogram was generated using the UPGMA clustering algorithm and, similar to the principal coordinate analysis (PCoA), it showed two groups that correspond to the geographic origin of the tarwi samples. AMOVA showed a reduced variation between clusters (7.59 %) and indicated that variability within populations is 92.41 %. Population divergence (Fst) between clusters 1 and 2 revealed low genetic difference (0.019). We also detected a negative Fis for both populations, demonstrating that, similar to other Lupinus species, tarwi also depends on cross-pollination. SNPs markers were powerful and effective for the genotyping process in this germplasm. We hope that this information is the beginning of the path towards a modern genetic improvement and conservation strategies of this important Andean legume.

Senescence-accelerated mouse prone 8 (SAMP8) mice exhibit cognitive defects and neuron loss with aging, and are used to study anti-aging effects of Dendrobium nobile alkaloids (DNLA). SAMP8 mice were orally given DNLA from ages 6 to 10 months. At 10 months of age, behavioral tests and neuron damage were evaluated. Protein was extracted and subjected to phosphorylated proteomic analysis. The cognitive deficits and neuron loss in hippocampus and cortex of aged SAMP8 mice were improved by DNLA. Hippocampal proteomic analysis showed differentially expressed protein/genes in SAMP8 compared to age-matched senescence-accelerated resistant mice, including altered tubulin binding, microtubule binding, etc. via Gene Oncology. KEGG revealed endocytosis, mRNA surveillance, tight junction, protein processing in endoplasmic reticulum, aldosterone synthesis and secretion, and glucagon signaling pathway changes. Upregulated protein/genes in the hippocampus of SAMP8 mice, such as Lmtk3, Usp10, Dzip1, Csnk2b, and Rtn1, were attenuated by DNLA; whereas downregulated protein/genes, such as Kctd16, Psd3, Bsn, Atxn2l, and Kif1a, were rescued by DNLA. The aberrant protein/gene expressions of SAMP8 mice were correlated with transcriptome changes of Alzheimer’s disease in the GEO database, and were attenuated by DNLA. Thus, DNLA improved cognitive dysfunction and ameliorated neuronal injury in aged SAMP8 mice, and attenuated aberrant protein/gene expressions.

Hashimoto’s thyroiditis (HT) is a common autoimmune disease, with its prevalence rapidly increasing. Both genetic and environmental risk factors contribute to the development of HT. Recently, viral infection has been suggested to act as a trigger of HT by eliciting the host immune response and subsequent autoreactivity. We analyzed features of HT through bioinformatics analysis so as to identify markers of HT development. We accessed public microarray data of HT patients from the Gene Expression Omnibus (GEO) and obtained differentially expressed genes (DEGs) under HT. Gene Ontology (GO) and KEGG pathway enrichment analyses were performed for functional clustering of our protein-protein interaction (PPI) network. Utilizing ranked gene lists, we performed Gene Set Enrichment Analysis (GSEA) using the clusterprofiler R package. By comparing the expression signatures of the huge perturbation database with the queried rank-ordered gene list, connectivity map (CMap) analysis was performed to screen potential therapeutic targets and agents. The gene expression profile of the HT group was in line with general characteristics of HT. Biological processes related to the immune response and viral infection pathways were obtained for the upregulated DEGs. GSEA results revealed activation of autoimmune disease-related pathways and several viral infection pathways. Autoimmune disease and viral infection pathways were highly interconnected by common genes, while the HLA genes which are shared by both were significantly upregulated. CMap analysis suggested that perturbagens, including SRRM1, NLK, and CCDC92, have the potential to reverse the HT expression profile. Several lines of evidence suggested that viral infection and the host immune response are activated during HT. Viral infection is suspected to act as a key trigger of HT by causing autoimmunity. SRRM1, an alternative splicing factor, which responds to viral activity, might serve as potential marker of HT.

We attempted to understand the cervical cancer patient samples through Whole Exome Sequencing. We derived the variants from raw reads via our in-house benchmarked pipeline and validated the variants by IGV. This is the first cervical cancer exome data from the Indian cohort. Background: Cervical cancer (CC) is caused mainly by persistent infections of high-risk HPV, reduced parity, and factors like a decrease in average socioeconomic levels. We keep this because no cervical cancer exome data is available from the Indian cohort. Methods: The CC patient clinical samples were initially subjected to preparation using Qiagen DNA extraction, quantified, and the library preparation using the Agilent target enrichment system. Further, exome capture by Illumina platform (100X), quality check, alignment, and variant calling followed by the downstream analysis and finally visualized by IGV. Results: We observed a large number of SNVs or mutations from an Indian perspective, such as KMT2C, OR4M1, PDPR2P, EPHB1, FAS, OPCML, MGST1, C1QTNF9, HS6ST3, OR4K2, PRPSAP2, KCNJ12, FIGNL1, SFXN1, BAGE2, ARVCF, NAMPT variants of significance, unknown significance and possibly significant are reported. Conclusion: For the first time, KMT2C is observed as a novel potent mutation and pathogenic, showing the variant position at (7:152265091, T>A, SNV 62478356) from the Indian context in CC. Further, we visualized and validated the mutations using the Integrative Genomic Viewer (IGV) browser. Finally, we discuss the inherent challenges through KMT2C mutations.

The dynamic interplay of signaling networks in most major cellular processes is characterized by the orchestration of reversible protein phosphorylation. Consequently, analytic methods like quantitative phospho-peptidomics has been pushed forward from a highly specialized edge-technique to a powerful and versatile platform for comprehensively analyzing the phosphorylation profile of living organisms. Despite enormous progress in instrumentation and bioinformatics, a major problem remains a high number of missing values caused by the experimental procedure due to either a random phospho-peptide enrichment selectivity or borderline signal intensities, which both cause the exclusion for fragmentation using the commonly applied data dependent acquisition (DDA) mode. Consequently, an incomplete dataset reduces confidence in the subsequent statistical bioinformatic processing. Here, we successfully applied data independent acquisition (DIA) by using the filamentous fungus Magnaporthe oryzae as model organism and could prove that while maintaining data quality (such as phosphosite and peptide sequence confidence), the data completeness increases dramatically. Since the method presented here reduces the LC-MS/MS analysis from 3 h to 1 h and increases the number of phoshosites identified up to 10-fold in contrast to published studies in fungi, we pushed the phospho-proteomic technique beyond its current limits and could provide a sophisticated resource for investigation of signaling processes in filamentous fungi.

: Cancer metastasis is the lethal developmental step in cancer, responsible for the majority of cancer deaths. To metastasize, cancer cells must acquire the ability to disseminate systemically and to escape an activated immune response. Here, we endeavoured to investigate if metastatic dissemination reflects acquisition of genomic traits that are selected for. We acquired mutation and copy number data from 8,332 tumours representing 19 cancer types acquired from The Cancer Genome Atlas and the Hartwig Medical Foundation. A total of 827,344 non-synonymous mutations across 8,332 tumour samples representing 19 cancer types were timed as early or late relative to copy number alterations, and potential driver events were annotated. We found that metastatic cancers had significantly higher proportion of clonal mutations and a general enrichment of early mutations in p53 and RTK/KRAS pathways. However, while individual pathways demonstrated a clear time-separated preference for specific events, the relative timing did not vary between primary and metastatic cancers. These results indicate that the selective pressure that drives cancer development does not change dramatically between primary and metastatic cancer on a genomic level, and is mainly focused on alterations that increase proliferation.

Sulfate transporters (SULTRs) are responsible for the uptake of the sulfate (SO42−) ions in the rhizosphere by the roots and their distribution in plant organs. In this study, SULTR family members in the genome of the two oilseed crops, Camelina sativa, and Brassica napus, were identified and characterized based on their sequence structure, duplication events, phylogenetic relationships, phosphorylation sites, and expression levels. Herein, 36 and 45 putative SULTR genes were recognized from the genome of C. sativa, and B. napus, respectively. SULTR proteins were predicted as basophilic proteins with low hydrophilicity in both studied species. According to phylogenetic relationships, we divided SULTRs into five groups, in which SULTRs 3 showed highest variation. Besides, several duplication events were observed between SULTRs. The first duplication event was predicted approximately five million years ago between three SULTRs 3.1 in C. sativa. Two subunits were indicated in the 3D structure of SULTRs that the active binding sites differed between C. sativa and B. napus. According to available RNA-seq data, SULTRs showed diverse expression in tissues and response to stimuli. SULTRs 3 showed an expression in all tissues. SULTRs 3.1 were more upregulated in response to abiotic stresses in C. sativa, while SULTRs 3.3, and SULTRs 2.1 showed an upregulation in B. napus. Furthermore, SULTRs 3 and SULTRs 4.1 showed an upregulation in response to biotic stresses in B. napus. Based on the distribution of cis-regulatory elements in the promoter region, we speculated that SULTRs might be controlled by phytohormones such as ABA, and MeJA. Therefore, it seems that SULTR genes in C. sativa have been more influenced by evolutionary processes and have acquired more diversity.

Single-cell RNA sequencing data facilitates investigation of cell heterogeneity and subpopulations as well as differentially abundant states however modern single-cell RNA sequencing datasets are growing in size and complexity requiring advances in the bioinformatic methods that analyze them. Many methods exist for each step of analysis including read alignment, normalization, quality control, batch effect correction, imputation and dimensionality reduction. With so many options to choose from at each step of the analysis, benchmarking and a synthesis of the literature on the methods available is necessary to inform biological researchers on the most optimal workflow for their data. Here, recent key methods of analysis are highlighted with a focus on methods that facilitate identification of cell subpopulations and differentially abundant cell states. With a constantly expanding toolset for each step in single-cell RNA sequencing dataset analysis, biological researchers should stay informed to utilize the most applicable methods for their own analyses.

Our research team compared the performance of matrix-assisted laser desorption/ionization followed by a time of flight (MALDI-TOF) mass spectrometry and genomic DNA extraction followed by sequencing, assembly and alignment for phylogenetic assessment. We performed these comparisons to determine our methodology’s overall efficacy and accuracy for environmental bacteria. In addition, we collected samples from various contaminated rivers in the Dominican Republic. For both methods, we analyzed these results and reported the main differences between each method.

The ability to directly look into genome sequences has opened great opportunities in plant breeding. Yet, the assembly of full-length chromosomes remains one of the most difficult problems in modern genomics. Genetic maps are commonly used in de novo genome assembly and are constructed on the basis of a statistical analysis of the number of recombination rather than physical distance in base pairs. This may affect the accuracy of the ordering and orientation of scaffolds within the chromosome, especially in the region of recombination suppression. Here we report the use of Tyr-FISH for validation of genetic and pseudochromosome maps. For probe design we developed a pipeline that is based on selection of a unique sequence with minimal potential fluorescent background arising from non-specific in situ hybridization. In total of 24 unique genes were located on physical chromosomes 2 and 6. The order of markers has been corrected by integration genetic and cytogenetic maps. Tyr-FISH mapping showed that the order of 23.1% (chromosome 2) and 27.3% (chromosome 6) of the tested genes differed between physical chromosomes and pseudochromosomes. Also the position of the mlh1 gene, which was not on the genetic map, was defined on physical chromosome 2. The results suggest the possibility of Tyr-FISH mapping of any gene that may not be neither on the genetic map nor in the assembly. Hence, Tyr-FISH provides valuable information for the improvement of the genome assembly.

Abstract: Since ancient times aging has also been regarded as a disease, and humankind has always strived to extend the natural lifespan. Analyzing the genes involved in aging and disease allows for finding important indicators and biological markers for pathologies and possible therapeutic targets. An example of the use of omics technologies is the research regarding aging and the rare and fatal premature aging syndrome progeria (Hutchinson-Gilford progeria syndrome, HGPS). In our study, we focused on the in silico analysis of differentially expressed genes (DEGs) in progeria and aging, using a publicly available RNA Seq dataset (GEO dataset GSE113957) and a variety of bioinformatics tools. We identified several genes that appear to be involved both in natural aging and progeria (KRT8, KRT18, ACKR4, CCL2, UCP2, ADAMTS15, ACTN4P1, WNT16, IGFBP2). Further analyzing these genes and the pathways involved confirmed their possible roles in aging, suggesting the need for further in vitro and in vivo research. The graphical abstract illustrates the analysis workflow we used and will introduce in the following as an example to demonstrate the power of omics and bioinformatics.

The evolution and the emergence of new mutations of viruses affect their transmissibility and/or pathogenicity features, depending on different evolutionary scenarios of virus adaptation to the host. A typical trade-off scenario of SARS-CoV-2 evolution has been proposed, which leads to the appearance of an Omicron strain with lowered lethality, yet enhanced transmissibility. This direction of evolution might be partly explained by virus adaptation to therapeutic agents and enhanced escape from vaccine-induced and natural immunity formed by other SARS-CoV-2 strains. Omicron’s high mutation rate in the Spike protein, as well as its previously described high genome mutation rate (Kandeel et al., 2021), revealed a gap between it and other SARS-CoV-2 strains, indicating the absence of a transitional evolutionary form to the Omicron strain. Therefore, Omicron has emerged as a new serotype and divergent from the evolutionary lineage of other SARS-CoV-2 strains. Omicron is a rapidly evolving variant of high concern, whose new subvariants continue to manifest. Its further understanding and the further monitoring of key mutations that provide virus immune escape and/or high affinity towards the receptor could be useful for vaccine and therapeutic development in order to control the evolutionary direction of the COVID-19 pandemic.

The generation of different types of defective viral genomes (DVG) is an unavoidable consequence of the error-prone replication of RNA viruses. In recent years, a particular class of DVGs, those containing long deletions or genome rearrangements, has gain interest due to their potential therapeutic and biotechnological applications. Identifying such DVGs in high-throughput sequencing data has become an interesting computational problem. Up to nowadays, several algorithms have been proposed, though all incur in false positives, a problem of practical interest if such DVGs have to be synthetized and tested in the laboratory. Here we develop a novel software, DVGfinder, that wraps the two most commonly used algorithms into a pipeline that predicts DVGs. Using a gradient boosting classifier machine learning algorithm, we evaluate the performance of DVGfinder compared to previous algorithms and found that it outcompetes their precision and sensitivity in simulated datasets. DVGfinder generates user-friendly output files in HTML format that can assist users to identify DVGs based on their associated probability of being true positives.

HIV-associated neurocognitive disorder (HAND) is an array of neurocognitive changes associated with HIV infection, and the roles of microRNAs in HAND are not completely revealed yet. Based on published data and publicly available databases, we constructed an integrated miRNA-mRNA network involved in HAND. Bioinformatics analyses, including gene ontology, network analysis, and KEGG pathway analysis, were applied for further study of the network and the genes of the network. The axon guidance KEGG pathway, three genes NTNG1, EFNB2, CXCL12, and 17 miRNAs which regulates them, are spotlighted in our study. This study provides new perspectives to the knowledge of miRNAs’ roles in the process of HAND, and our findings provided potential therapeutic targets and clues of HAND.

Background: Recent researches have investigated the biological importance of RNA N1-methyladenosine (m1A) modifications in oncogenesis and progression of head and neck squamous cell carcinoma (HNSCC). However, whether m1A modifications also have latent effect in tumor microenvironment (TME) generation and immune regulation in HNSCC is unknown. Methods: We evaluated the m1A modification patterns and related to these modification patterns with TME cell infiltration features in 1041 HNSCC samples by bioinformatics approach. Results: The m1A score is an independent prognostic indicator. HNSCC patients with low m1A score group with poor overall survival (OS) was mainly characterized by stroma activation, lack of sufficient immune infiltration, and exhibited an immune- desert TME phenotype. Low m1A scores were also correlated with increased tumor mutation burden (TMB), and HNSCC patients with high TMB and low m1A scores had the worst OS. In addition, anti-CTLA-4 combined with anti-PD1 treatment was more effective in the high m1A score subgroup than in the low m1A score subgroup. Conclusions: This study revealed that m1A modifications play a non-negligible role in developing the TME versatility and complexity of HNSCC. Assessing m1A modification patterns in HNSCC helps improve our comprehension of its TME infiltration profile and guides more effective immunotherapeutic approaches.

Various kinds of primary metabolisms in plants are modulated through sulfate assimilation that the uptake of this inorganic compound can be regulated via the sulfate transporters, such as sulfotransfer-ases (SOTs), engaged in the sulfur metabolism. In the current study a genome-wide approach has been utilized for recognition and characterization of SOT family genes in the significant nutritional crop po-tato (Solanum tuberosum L.). As a result, 29 StSOT genes were identified in the potato genome, which were mapped onto the nine S. tuberosum chromosomes. The protein motifs structure demonstrated two highly conserved 5' PSB region and 3' PB motif that are essential for sulfotransferase and catalytic ac-tivities. The protein-protein interaction networks also significantly demonstrated an interesting collabo-ration between SOTs and the other genes, such as PRTase, APS-kinase, protein phosphatase and APRs, in sulfur compounds biosynthesis and regulation of the flavonoid and brassinosteroid metabolic pro-cesses, which clearly detected the importance of sulfotransferases for potato proper growth/development and stress dealing. Notably, the homology modeling of StSOT proteins and dock-ing analysis of their ligand-binding sites revealed the presence of some stress-responsive residues, such as proline, glycine, serine and lysine, in their active sites. The expression assay of StSOT genes via the potato RNA-seq data clearly suggested the engagements of these gene family members in plants growth and extension as well as responses to various hormones and biotic/abiotic stimulus circum-stances. Our predictions can be informative for the functional characterization of the SOT genes in po-tato and may the other nutritional crops.

The SARS-CoV-2 polybasic furin cleavage site is still a missing link. Remarkably, the two arginine residues of this protease recognition site are encoded by the CGG codon, which is rare in Betacoronavirus. However, the arginine pair is common at viral furin cleavage sites, but are not CGG-CGG encoded. The question is: Is this genetic footprint unique to the SARS-CoV-2? To address the issue, using Perl scripts, here I dissect in detail the NCBI Virus database in order to report the arginine dimers of the Betacoronavirus proteins. The main result reveals that a group of Middle East respiratory syndrome-related coronavirus (MERS-CoV) (isolates: camel/Nigeria/NVx/2016, host: Camelus dromedarius) also have the CGG-CGG arginine pair in the spike protein polybasic furin cleavage region. In addition, CGG-CGG encoded arginine pairs were found in the orf1ab polyprotein from HKU9 and HKU14 Betacoronavirus, as well as, in the nucleocapsid phosphoprotein from few SARS-CoV-2 isolates. To quantify the probability of finding the arginine CGG-CGG codon pair in Betacoronavirus, the likelihood ratio (LR) and a Markov model were defined. In conclusion, it is highly unlikely to find this genetic marker in betacoronaviruses wildlife, but they are there. Collectively, results shed light on recombination as origin of the virus CGG-CGG arginine pair in the S1/S2 cleavage site.

Interactome analyses have traditionally been applied to yeast, human and other model organisms due to the availability of protein-protein interactions data for these species. Recently these techniques have been applied to more diverse species using computational interaction prediction from genome sequence and other data types. This review describes the various types of computational interactome networks that can be created and how they have been used in diverse eukaryotic species, highlighting some of the key interactome studies in non-model organisms.

Identification of the differentially-expressed genes is important for clarification of the complex molecular mechanisms under drought conditions. In this experiment, transcriptome profiles of sensitive and tolerant tomato genotypes under drought stress were analyzed. Three up-regulated genes were selected, included CAB3 (Chlorophyll a-b binding protein3), SAMDC (S-adenosylmethionine decarboxylase proenzyme), and ACS9 (1-aminocyclopropane-1-carboxylate synthase 9). After bioinformatics analysis, tomato genotypes were subjected to drought stress and gene expression was determined using Real-Time-PCR. Physiological parameters of genotypes were also measured by spectrophotometer-based methods. According to the results, these three genes play a key role in stress tolerance. Expression of the CAB3 gene in both sensitive and tolerant genotypes was not significantly different compared to the control, but the SAMDC gene decreased in both genotypes and the ACS9 gene decreased in sensitive genotype and increased in tolerant genotype. The physiological analysis also showed that under stress conditions, the photosynthetic system of the plant was disrupted and the chlorophyll content was reduced, but, proline content and antioxidant enzymes activity increased, in which their quantity in the tolerant genotype was significantly higher than sensitive. Under drought stress, due to damage to the lipid membrane, Malondialdehyde content also increased, in which the sensitive genotype was more affected.

Gene expression data is key for the functional annotation of single nucleotide polymorphisms (SNPs) identified in genome-wide association studies (GWAS). Expression and splicing quantitative trait loci (e/sQTLs) in normal colon tissue, such as those from the University of Barcelona and University of Virginia RNA sequencing project (BarcUVa-Seq) and the Genotype-Tissue Expression project (GTEx), are required to gain biological insight of colon-related diseases risk loci. Moreover, transcriptome-wide association studies (TWAS) rely on reference gene expression imputation panels in the tissue of interest to nominate susceptibility genes. Also, it is of high interest to study the relationships between genes in a network framework. For facilitating these analyses, we have updated and expanded the scope of the Colon Transcriptome Explorer (CoTrEx) to the version 2.0. This web-based resource provides exhaustive visualization and analysis of transcriptome-wide gene expression profiles of normal colon tissue from BarcUVa-Seq and GTEx. In addition to the integration of new datasets, CoTrEx 2.0 provides additional e/sQTLs sets, as well as gene expression prediction models and regulatory and co-expression networks. It is freely available at https://barcuvaseq.org/cotrex/. Overall, it is of high interest for researchers aiming to investigate the genetic susceptibility to colon-related complex traits and diseases.

In 1984, Susumu Ohno hypothesized that the nylon-degrading enzyme NylB arose de novo via a frameshift mutation within a hypothetical precursor protein (PR.C). However, Ohno never tested his hypothesis or provided supporting biological evidence. For decades, Ohno’s famous frame-shift hypothesis has been uncritically accepted as the correct explanation for the origin of NylB and has been used to illustrate how simple it is for a totally new enzyme to arise spontaneously. In this paper we test Ohno’s hypothesis in light of data not available in 1984. We searched multiple protein databases and found that the NylB protein is widely occurring, has thousands of homologs, and is found in diverse organisms and diverse habitats. Conserved domain searches showed that the NylB sequence is homologous to beta lactamases - a family of highly conserved enzymes. However, our searches showed that there is no evidence for the existence of Ohno’s hypothetical PR.C protein, nor any credible homolog. Our results effectively falsify Ohno's frameshift hypothesis. We extended this analysis to other nylonases and found all the nylonases we examined had large numbers of homologs throughout the biosphere. This falsifies the long-held assumption that all nylonases evolved after the invention of nylon in 1935.

Distinguishing circular RNA (circRNA) reads from reads derived from the linear host transcript is a challenging task because of sequence overlap. We developed a computational approach, CiLiQuant, that determines the relative circular and linear abundance of transcripts and gene loci using backsplice and forward splice junction reads generated by existing mapping and circRNA discovery tools.

Next-generation sequencing (NGS) has actualized human papillomavirus (HPV) virome profiling for in-depth investigation of viral evolution and pathogenesis. However, viral computational analysis remains a bottleneck due to semantic discrepancies between computational tools and curated reference genomes. To address this, we developed and tested automated workflows for HPV taxonomic profiling and visualization using a customized Papillomavirus database in CLC Microbial Genomics Module. HPV genomes from Papilloma Virus Episteme were customized and incorporated into CLC “ready-to-use” workflows for stepwise data processing to include: 1) Taxonomic Analysis, 2) Estimate Alpha/Beta Diversities, and 3) Map Reads to Reference. Low-grade (n = 95) and high-grade (n = 60) Pap smears were tested with ensuing collective runtimes: Taxonomic Analysis (36 min); Alpha/Beta Diversities (5 sec); Map Reads (45 min). Tabular output conversion to visualizations entailed 1-2 keystrokes. Biodiversity analysis between low- (LSIL) and high-grade squamous intraepithelial lesions (HSIL) revealed loss of species richness and gain of dominance by HPV-16 in HSIL. Integrating clinically relevant, taxonomized HPV reference genomes within automated workflows proved to be an ultra-fast method of virome profiling. The entire process named “HPV DeepSeq” provides a simple, accurate and practical means of NGS data analysis for a broad range of applications in viral research.

Members of the DUF34 (domain of unknown function 34) family, also known as the NIF3 protein superfamily, are ubiquitous across superkingdoms. Proteins of this family have been widely annotated as “GTP cyclohydrolase I type 2” through electronic propagation based on one study. Here, the annotation status of this protein family was examined through comprehensive literature review and integrative bioinformatic analyses that revealed varied pleiotropic associations and phenotypes. This analysis combined with functional complementation studies strongly challenges the current annotation and suggests that DUF34 family members may serve as metal ion insertases, chaperones, or metallocofactor maturases. This general molecular function could explain how DUF34 subgroups participate in highly diversified pathways such as cell differentiation, metal ion homeostasis, pathogen virulence, redox and universal stress responses.

With a large number of annotated non-coding RNAs (ncRNAs), repetitive sequences are found to constitute functional components (termed as repetitive elements) in ncRNAs that perform specific biological functions. Bioinformatics analysis is a powerful tool for improving our understanding of the role of repetitive elements in ncRNAs. This chapter summarizes recent findings that reveal the role of repetitive elements in ncRNAs. Furthermore, relevant bioinformatics approaches are systematically reviewed, which promises to provide valuable resources for studying the functional impact of repetitive elements on ncRNAs.

The furin cleavage site, with an arginine doublet (RR), is one of the clues of the SARS-CoV-2 origin. This furin-RR is encoded by the CGG-CGG sequence. Because arginine can be encoded by six codons, in a previous work we found that in SARS-CoV-2, CGG was the minority arginine codon (3%). Also, analyzing the RR doublet from a large sample of furin cleavage sites of several kinds of viruses, we found that none of them were encoded by CGG-CGG. Here, we come back to the core of the matter, but from the perspective that in the human genome, in contrast, CGG is the majoroty arginine codon (21%). Here, we highlighted that the 6 arginine codons provide genetic markers to a traceability on the RR origin in the furin site, as well as, to weigh the probability of the theories about the origin of the virus.

Artificial intelligence AI or machine learning has proven to be a potential activity in the health and biomedical sciences. Previous research it has found that AI can learn new data and transform it into the useful knowledge. In the field of pharmacology, the aim is to design more efficient and novel vaccines using this method which are also cost effective. The underlying fact is to predict the molecular mechanism and structure for increased likelihood of developing new drugs. Clinical, electronic and high resolution imaging datasets can be used as inputs to aid the drug development niche. Moreover, the use of comprehensive target activity has been performed for repurposing a drug molecule by extending target profiles of drugs which also include off targets with therapeutic potential providing a new indication.

Systems biology has been established for more than a decade in the post-genomic era. With the help of the computational and mathematical tools, systems biology reconstitutes the entire scenario of the cell, tissue and even organism from the pieces data generated in the past decades. However, the modern biology is mainly focusing on the structure and function of the biomolecule, cell, tissue or organ, which are far from the essence of the life because of missing thermodynamic information. It is doubtable that the current systems biology-based omics is no-how to fully understand the dynamic courses of the structure, function and information in life. For this reason, we promote a novel concept of aquamoleculomics, in which the biological structure and function as well as thermodynamic characteristics and bioinformation of the aquamolecule complexes are included in this theoretical model of systems biology. Water is mother of life, matter and matrix of organism. Indeed, the fundamental roles of H₂O molecules in biological processes might be dramatically underestimated. Extremely speaking, H₂O networks in the living system might be engaged in all the biological processes including building all the biological structures, the residential places of the motherhood molecules as the honeycombs of honeybees.

Drug repurposing, the practice of utilizing existing drugs for novel clinical indications, has tremendous potential for improving human health outcomes and increasing therapeutic development efficiency. The goal of multidisease multitarget drug repurposing, also known as shotgun drug repurposing, is to develop platforms that assess the therapeutic potential of each existing drug for every clinical indication. Our Computational Analysis of Novel Drug Opportunities (CANDO) platform for shotgun multitarget repurposing implements several pipelines via large scale modelling and simulation of interactions between comprehensive libraries of drugs/compounds and protein structures. In these pipelines, each drug is described by an interaction signature that is then compared to all other signatures that are then sorted and ranked based on similarity. Pipelines within the platform are benchmarked based on their ability to recover known drugs for all indications in our library, and predictions are generated based on the hypothesis that (novel) drugs with similar signatures may be repurposed for the same indication(s). The drug-protein interactions in the platform used to create the drug-proteome signatures may be determined by any screening or docking method but the primary approach used thus far has been an in house similarity docking protocol. In this study, we calculated drug-proteome interaction signatures using the publicly available molecular docking method Autodock Vina and created hybrid decision tree pipelines that combined our original bio- and cheminformatic approach with the goal of assessing and benchmarking their drug repurposing capabilities and performance. The hybrid decision tree pipeline outperformed the corresponding two docking-based pipelines it was synthesized from, yielding an average indication accuracy of 13.3% at the top10 cutoff (the most stringent), relative to 10.9% and 7.1% for its constituent pipelines, and a random control accuracy of 2.2%. We demonstrate that docking based virtual screening pipelines have unique performance characteristics and that the CANDO shotgun repurposing paradigm is not dependent on a specific docking method. Our results also provide further evidence that multiple CANDO pipelines can be synthesized to enhance drug repurposing predictive capability relative to their constituent pipelines. Overall, this study indicates that pipelines consisting of varied docking based signature generation methods can capture unique and useful signal for accurate comparison of drug-proteome interaction signatures, leading to improvements in the benchmarking and predictive performance of the CANDO shotgun drug repurposing platform.

Rapid development of high-throughput omics technologies generates an increasing interests in algorithms for cutoff point identification. Existing cutoff methods and tools identify cutoff points based on association of continuous variables with another variable, such as phenotype, disease state or treatment group. These approaches are not applicable for descriptive studies in which continuous variables are reported without known association with any biologically meaningful variables. The most common shape of the ranked distribution of continuous variables in high-throughput descriptive studies corresponds to a biphasic exponential/super-exponential curve, where the first phase includes big number of variables with values slowly growing with rank and the second phase includes smaller number of variables rapidly growing with rank. This study describes an easy algorithm to identify the boundary between these phases to be used as a cutoff point. The major assumption of that approach is that small number of variables with high values dominate biological system and determine its major processes and functions. This approach was tested on three different datasets: genes in the human cerebral cortex, mammalian genes sensitive to chemical exposures, and proteins expressed in human heart. In every case, the described cutoff identification method produced shortlists of variables (genes, proteins) highly relevant for dominant functions/pathways of the analyzed biological systems. Thus, our described method for cutoff identification may be used to prioritize variables for a focused functional analysis, in situations where other methods of dichotomization of data are inaccessible.

The diagnosis of papillary thyroid carcinoma has always been a concerned and challenging issue and it is very important and meaningful to have a definite diagnosis before the operation. In this study, we tried to use an artificial intelligence algorithm instead of medical statistics to analyze the genetic fingerprint from gene chip results to identify papillary thyroid carcinoma. We trained 20 artificial neural network models with differential genes and other important genes related to the cell metabolic cycle as the list of input features, and apply them to the diagnosis of papillary thyroid cancer in the independent validation data set. The results showed that when we used the DEGs and all genes lists as input features the models got the best diagnostic performance with AUC=98.97% and 99.37% and the accuracy were both 96%. This study revealed that the proposed artificial neural network models constructed with genetic fingerprints could achieve a prediction of papillary thyroid carcinoma. Such models can support clinicians to make more accurate clinical diagnoses. At the same time, it provides a novel idea for the application of artificial intelligence in clinical medicine.

The exponential growth of biomedical data in recent years urged the application of numerous machine learning techniques to address emerging problems in biology and clinical research. By enabling automatic feature extraction, selection and generation of predictive models, these methods can be used to efficiently study complex biological systems. Machine learning techniques are frequently integrated with bioinformatic methods, as well as curated databases and biological networks, to enhance training and validation, identify the best interpretable features, and enable feature and model investigation. Here, we review recently developed methods that incorporate machine learning within the same framework with techniques from molecular evolution, protein structure analysis, systems biology and disease genomics. We outline the challenges posed for machine learning, and in particular, deep learning in biomedicine and suggest unique opportunities for machine learning techniques integrated with established bioinformatics approaches to overcome some of these challenges.

The SARS-CoV-2 high infectivity is due to the functional polybasic furin cleavage site in the S protein. How it was acquired is unknown. There are two challenges to face: (i) an evolutionary model, to fit the origin of the coronavirus; and (ii) a molecular mechanism for the site acquisition. Here we show genomic fingerprints which are specific of Pangolin-CoVs, Bat-SARS-like (CoVZC45, CoVZXC21), bat RatG13 and human SARS-CoV-2 coronaviruses. This, along with phylogenetic analysis, we found that these species have the same evolutionary origin in the bat, including a genetic recombination of S gene between Pangolin-CoV (2017) and RatG13 ancestors. However, this does not explain why SARS-CoV-2 is the only of them with the furin site, which consists in four amino acid (PRRA) motif. The Arginine doublet is encoded by CGGCGG codons. Surprisingly, none of the Arginine doublet of other furin site of viral proteins from several type of viruses, are encoded by the CGGCGG codons. This makes it difficult to consider a virus recombination as mechanism for the PRRA acquisition. The origin of SARS-CoV-2, is the origin of the recognition cleavage site. The bat coronavirus RaTG13 appears to be the closest relative of the SARS-CoV-2, but was isolated in 2013. So, new RatG13 samples would provide insights into the acquisition of the polybasic motif.

The clinical interaction between stroke and periodontitis has been consistently studied and confirmed. Hence, forecasting potentially new protein interactions in this association using bioinformatic strategies presents potential interest. In this exploratory study, we conducted a protein-protein network interaction (PPI) search with documented encoded proteins for both stroke and periodontitis. Genes of interest were collected via GWAS database. The STRING database was used to predict the PPI networks, first in a sensitivity purpose (confidence cut-off of 0.7), and then with a highest confidence cut-off (0.9). Genes over-representation was inspected in the final network. As a result, we foresee a prospective protein network of interaction between stroke and periodontitis. Inflammation, pro-coagulant/pro-thrombotic state and ultimately atheroma plaque rupture is the main biological mechanism derived from the network. These pilot results may pave the way to future molecular and therapeutic studies to further comprehend the mechanisms between these two conditions.

The usage of Quantum Similarity through the equation Z = {∀θ ∈ Z → ∃s ∈ S ∧ ∃t ∈ T : θ = (s, t)}, represents a way to analyze the way communication works in our DNA. Being able to create the object set reference for z being (s, t) in our DNA strands, we are able to set logical tags and representations of our DNA in a completely computational form. This will allow us to have a better understanding of the sequences that happen in our DNA. With this approach, we can also utilize mathematical formulas such as the Euler–Mascheroni constant, regression analysis, and computational proofs to answer important questions on Quantum biology, Quantum similarity, and Theoretical Physics.

High-quality chromosome-scale haplotype sequences— of diploid genomes, polyploid genomes and metagenomes — provide important insights into genetic variation associated with disease and biodiversity. However, whole-genome short read sequencing does not yield haplotype information that spans whole chromosomes directly. Computational assembly of shorter haplotype fragments is required for haplotype reconstruction, which can be challenging owing to limited fragment lengths and high haplotype and repeat variability across genomes. Recent advancements in long-read and chromosome-scale sequencing technologies, alongside computational innovations, are improving the reconstruction of haplotypes at the level of whole chromosomes. Here, we review recent methodological progress in these areas and discuss perspectives that could enable routine high-quality haplotype reconstruction in clinical and evolutionary studies.

Over the past decade, the One Strain Many Compounds (OSMAC) approach has been established for silent gene cluster activation and elicitation of secondary metabolite production, but so far the full secondary metabolome of a biosynthetically promising bacterium has not been elucidated yet. Here, we investigate the ability of seven categories of OSMAC conditions to enhance the diversity of new mass features from bacterial strains with little literature coverage but high biosynthetic potential. The strains Bacillus. amyloliquefaciens DSM7, Corallococcus. coralloides DSM2259, Pyxidicoccus. fallax HKI727, Rhodococcus. jostii DSM44719, and Streptomyces. griseochromogenes DSM40499 were selected after genome mining with antiSMASH. After cultivation under OSMAC conditions, the generated extracts were subjected to LC-MS and MZmine analysis to determine new mass features and evaluate the tested culture conditions. 4 predicted compounds, bacillibactin, desferrioxamine B, myxochelin A, and surfactin, were identified and up to 147 new mass features were detected in the generated extracts, which greatly surpasses the number of predicted gene clusters. Among the new mass features are bioactive compounds that were so far unreported for the strains such as cyclo-(Tyr-Pro) from DSM7 and nocardamin from DSM2259. Furthermore, the tested culture conditions were evaluated regarding their suitability for the generation of new mass features from the selected strains and promising new starting points for further screenings are postulated. Especially culture conditions with little prior literature coverage are responsible for the activation of predicted gene clusters

To fully appreciate genetics, one must understand the link between genotype (DNA sequence) and phenotype (observable characteristics). Advances in high-throughput genomic sequencing technologies and applications, so-called “-omics”, have made genetic sequencing readily available across fields in biology from applications in non-traditional study organisms to precision medicine. Thus, understanding these tools is critical for any biologist, especially those early in their career. This comprehensive review discusses the chronological development of different sequencing methods, the bioinformatics steps to analyzing this data, and social and ethical issues raised by these techniques that must be discussed and evaluated.

Hepatitis C Virus (HCV) is the main causative factor for liver cirrhosis and the development of liver cancer, with a confirmed ~ 180 million infections worldwide. E2 is an HCV structural protein responsible for virus entry to the host cell. Heat Shock Protein A5 (HSPA5), also termed BiP and GRP78, is the master regulator of the unfolded protein response mechanism, where it mainly localizes in the lumen of the Endoplasmic Reticulum (ER) in normal conditions. Under the stress of HCV infection or carcinogenesis, HSPA5 is upregulated. Consequently, HSPA5 escapes the ER retention localization and translocates to the cytoplasm and plasma membrane. Pep42, a cyclic peptide that was reported to target explicitly cell-surface HSPA5 in vivo. Owing to the high sequence and structural conservation between the C554-C566 region of HCV E2 and the Pep42, then we propose that the HCV E2 C554-C566 region could be the recognition site. The motivation of this work is to predict the possible binding mode between HCV E2 and HSPA5 by implementing molecular docking to test such proposed binding. Docking results reveal the high potent binding of the HCV E2 C554-C566 region to HSPA5 substrate-binding domain β (SBDβ). Moreover, the full-length HCV E2 also exhibits high binding potency to HSPA5 SBDβ. Defining the binding mode between HCV E2 and HSPA5 is of significance, so one can interfere with such binding and reducing the viral infection.

The subtilase family is the second largest family of serine proteases. Some fungi including Trichoderma species can capture and kill nematodes by secreting hydrolytic enzymes or toxins, among which serine proteases are important enzymes that allow fungi to infect nematodes. Subtilase can degrade nematode and insect body walls. In this study, subtilase family genes were identified from the Trichoderma harzianum genome database, and bioinformatics analysis of the characteristics and evolutionary status of these genes, along with structural and functional analyses of their proteins, was performed. Gene structure analysis revealed that all the 41 subtilase genes contained introns, while some did not have upstream or downstream regions. Chromosome localisation showed that subtilase family members were unevenly distributed in 22 Trichoderma chromosomes, and 3 clusters were present, indicating that they may be hot spots of subtilase genes. Conserved motif analyses showed these proteins contained a commonly conserved motif, and motifs belonging to the same subfamily remained highly similar. The upstream region of the subtilase genes were enriched with different type and numbers of cis-elements, indicating that subtilase genes are likely to play a role in the response to diverse stresses. Transcription of 31 genes was increased after 5 days of infection with nematodes, whereas that of 10 genes decreased. In these subtilase genes, ThSBT4, ThSBT5, ThSBT12, ThSBT27, ThSBT34, ThSBT35, ThSBT38, and ThSBT40 showed significantly upregulated expression with a log2 fold change value of more than 4, and ThSBT35 showed the highest peak. These results laid a theoretical foundation for further research on the function of the subtilase genes and the mechanism of the resistance response.

Latest evidence revealed a possible association between Parkinson’s disease (PD) and periodontitis. We explored the causal relationship of this association through two-sample Mendelian randomization (MR) in European ancestry populations. To this end, we used openly accessible data of genome-wide association studies (GWAS) on PD and periodontitis. As instrumental variables for periodontitis, seventeen single-nucleotide polymorphisms (SNPs) from a GWAS of periodontitis (1817 periodontitis cases vs. 2215 controls) and forty-five SNPs from a GWAS of PD (20,184 cases and 397,324 controls). Eight non-overlapping SNPs of periodontitis from an additional GWAS assisted in the validation of association being studied. Multiple approaches of MR were carried-out. There was no evidence of genetic liability of periodontitis being associated with a higher risk of PD (B= -0.0003, Standard Error [SE] 0.0003, P = 0.26). The eight independent SNPs (B= -0.0000, SE 0.0001, P = 0.99) validated this outcome. We found no association of genetically primed PD towards periodontitis (B= -0.0001, SE 0.0001, P = 0.19). This MR study found no conclusive evidence to support a bidirectional causal genetic liability between PD and periodontitis. Further GWAS studies are needed to confirm the consistency of these results.