Search PubMedSearch

SEARCH · Search PubMed

Results for “Python”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5Linked to original sources

Riboflavin-binding protein from reptiles: a comparison with avian riboflavin-binding proteins.

1. Riboflavin-binding protein (RBP) has been isolated for the first time from reptilian sources. 2. RBP from eggs of Python molurus (Indian python) and Chrysemys picta (painted turtle) has been isolated and compared to RBP from Gallus gallus domesticus (chicken), a well-characterized protein, and a newly isolated RBP from Cairina moschata (Muscovy duck). 3. Each of the proteins is phosphorylated and glycosylated. 4. The ratio of riboflavin binding to protein is 1:1 and the KD for each protein is between 1-3 nM. 5. The mol. wts, different for each species, range from 30,000-40,000, with the reptilian proteins being approx. 10,000 larger than the avian proteins.

Amino Acids

GeneFEAST: the pivotal, gene-centric step in functional enrichment analysis interpretation.

SUMMARY: GeneFEAST, implemented in Python, is a gene-centric functional enrichment analysis summarization and visualization tool that can be applied to large functional enrichment analysis (FEA) results arising from upstream FEA pipelines. It produces a systematic, navigable HTML report, making it easy to identify sets of genes putatively driving multiple enrichments and to explore gene-level quantitative data first used to identify input genes. Further, GeneFEAST can juxtapose FEA results from multiple studies, making it possible to highlight patterns of gene expression amongst genes that are differentially expressed in at least one of multiple conditions, and which give rise to shared enrichments under those conditions. Thus, GeneFEAST offers a novel, effective way to address the complexities of linking up many overlapping FEA results to their underlying genes and data, advancing gene-centric hypotheses, and providing pivotal information for downstream validation experiments. AVAILABILITY AND IMPLEMENTATION: GeneFEAST GitHub repository: https://github.com/avigailtaylor/GeneFEAST; Zenodo record: 10.5281/zenodo.14753734; Python Package Index: https://pypi.org/project/genefeast; Docker container: ghcr.io/avigailtaylor/genefeast.

Software

Predicting coarse-grained representations of biogeochemical cycles from metabarcoding data.

MOTIVATION: Taxonomic analysis of environmental microbial communities is now routinely performed thanks to advances in DNA sequencing. Determining the role of these communities in global biogeochemical cycles requires the identification of their metabolic functions, such as hydrogen oxidation, sulfur reduction, and carbon fixation. These functions can be directly inferred from metagenomics data, but in many environmental applications metabarcoding is still the method of choice. The reconstruction of metabolic functions from metabarcoding data and their integration into coarse-grained representations of biogeochemical cycles remains a difficult bioinformatics problem today. RESULTS: We developed a pipeline, called Tabigecy, which exploits taxonomic affiliations to predict metabolic functions constituting biogeochemical cycles. In a first step, Tabigecy uses the tool EsMeCaTa to predict consensus proteomes from input affiliations. To optimize this process, we generated a precomputed database containing information about 2404 taxa from UniProt. The consensus proteomes are searched using bigecyhmm, a newly developed Python package relying on Hidden Markov Models to identify key enzymes involved in metabolic function of biogeochemical cycles. The metabolic functions are then projected on coarse-grained representation of the cycles. We applied Tabigecy to two salt cavern datasets and validated its predictions with microbial activity and hydrochemistry measurements performed on the samples. The results highlight the utility of the approach to investigate the impact of microbial communities on biogeochemical processes. AVAILABILITY AND IMPLEMENTATION: The Tabigecy pipeline is available at https://github.com/ArnaudBelcour/tabigecy. The Python package bigecyhmm and the precomputed EsMeCaTa database are also separately available at https://github.com/ArnaudBelcour/bigecyhmm and https://doi.org/10.5281/zenodo.13354073, respectively.

Metagenomics

RabbitSketch: a high-performance sketching library for genome analysis.

SUMMARY: We present RabbitSketch, a highly optimized library of sketching algorithms such as MinHash, OrderMinHash, and HyperLogLog that can exploit the power of modern multi-core CPUs. It provides significant speedups compared to existing implementations, ranging from 2.30× to 49.55×, as well as flexible and easy-to-use interfaces for both Python and C++. As a result, the similarity analysis of 455GB genomic data can be completed in only 5 minutes using RabbitSketch with merely 20 lines of Python code. As a case study, we enhanced RabbitTClust by integrating RabbitSketch's Kssd algorithm, resulting in a 1.54× speedup with no loss in accuracy. AVAILABILITY AND IMPLEMENTATION: RabbitSketch is available at https://github.com/RabbitBio/RabbitSketch with an archived version at Zenodo: https://doi.org/10.5281/zenodo.14903962. Detailed API documentation is available at https://rabbitsketch.readthedocs.io/en/latest.

Software

FuNTB: a functional network clustering tool for the analysis of genome-wide genetic variants in Mycobacterium tuberculosis.

MOTIVATION: Tuberculosis (TB), caused by Mycobacterium tuberculosis (Mtb), still claims around 1.25 million lives each year. The growing threat of drug resistance-often driven by single‑nucleotide polymorphisms (SNPs) in Mtb genomes underscores the need for high‑quality genomic data and powerful bioinformatics tools. We present FuNTB, a python‑based pipeline that detects non‑synonymous SNPs in Mtb and builds functional network clusters to reveal genotype-phenotype relationships. RESULTS: FuNTB profiles non‑synonymous SNPs at the gene level across user‑defined phenotypes, pinpointing both shared and unique mutations. It ingests annotated Variant Call Format (VCF) files or MTBseq outputs and merges them with clinical metadata to produce network‑XML files compatible with Cytoscape and Gephi. When applied to the CRyPTIC Mtb collection, FuNTB rapidly recovered established resistance genes and surfaced novel candidates, validating its utility for mapping genotype-phenotype associations. AVAILABILITY AND IMPLEMENTATION: FuNTB is implemented in Python 3.8+ and is freely available under the MIT license at https://doi.org/10.5281/zenodo.15399917.

Mycobacterium tuberculosis

Fast and flexible minimizer digestion with digest.

SUMMARY: Minimizer digestion is an increasingly common component of bioinformatics tools, including tools for de Bruijn graph assembly and sequence classification. We describe a new open source tool and library to facilitate efficient digestion of genomic sequences. It can produce digests based on the related ideas of minimizers, modimizers or syncmers. Digest uses efficient data structures, scales well to many threads, and produces digests with expected spacings between digested elements. AVAILABILITY AND IMPLEMENTATION: Digest is implemented in C++17 with a Python API, and is available open-source at https://github.com/VeryAmazed/digest. The python library is available on Bioconda. Rust bindings are available as a public crate at https://crates.io/crates/digest-rs.

Software

Tsbrowse: an interactive browser for ancestral recombination graphs.

SUMMARY: Ancestral recombination graphs (ARGs) represent the interwoven paths of genetic ancestry of a set of recombining sequences. The ability to capture the evolutionary history of samples makes ARGs valuable in a wide range of applications in population and statistical genetics. ARG-based approaches are increasingly becoming a part of genetic data analysis pipelines due to breakthroughs enabling ARG inference at biobank-scale. However, there is a lack of visualization tools, which are crucial for validating inferences and generating hypotheses. We present tsbrowse, an open-source, web-based Python application for the interactive visualization of the fundamental building blocks of ARGs, i.e. nodes, edges and mutations. We demonstrate the application of tsbrowse to various data sources and scenarios, and highlight its key features of browsability along the genome, user interactivity, and scalability to very large sample sizes. AVAILABILITY AND IMPLEMENTATION: Tsbrowse is installed as a Python package from PyPI (https://pypi.org/project/tsbrowse/), while a development version is maintained at https://github.com/tskit-dev/tsbrowse. Documentation is available at https://tskit.dev/tsbrowse/docs/. Source code is archived on Zenodo with DOI, https://doi.org/10.5281/zenodo.15683039.

Software

ShortCake: an integrated platform for efficient and reproducible single-cell analysis.

SUMMARY: Recent advances in single-cell analysis have introduced new computational challenges. Researchers often need to use multiple analysis tools written in different programming languages while managing version conflicts between related packages within a single workflow. For the research community, minimizing the time spent on environment setup and installation issues is essential. We present ShortCake, a containerized platform that integrates a suite of single-cell analysis tools written in R and Python. ShortCake isolates competing Python tools into separate virtual environments that can be easily accessed within a Jupyter notebook. This enables users to effortlessly transition between various environments, including R, even within a single notebook. Additionally, ShortCake offers multiple "flavors," enabling users to select container images tailored to their specific needs. ShortCake provides a unified environment with fixed versions of various tools, thus streamlining workflows, reducing setup time, and improving reproducibility. AVAILABILITY AND IMPLEMENTATION: The ShortCake image is available on DockerHub (https://hub.docker.com/r/rnakato/shortcake) and Zenodo (DOIs: 10.5281/zenodo.17116765 and 10.5281/zenodo.17118158). The source code is available on GitHub (https://github.com/rnakato/ShortCake).

Single-Cell Analysis

Tomtom-lite: accelerating Tomtom enables large-scale and real-time motif similarity scoring.

SUMMARY: Pairwise sequence similarity is a core operation in genomic analysis, yet most attention has been given to sequences made up of discrete characters. With the growing prevalence of machine learning, calculating similarities for sequences of continuous representations, e.g. frequency-based position-weight matrices (PWMs) and attribution-based contribution-weight matrices, is taking on newfound importance. Tomtom has previously been proposed as an algorithm for identifying pairs of PWMs whose similarity is statistically significant, but the implementation remains inefficient for both real-time and large-scale analysis. Accordingly, we have re-implemented Tomtom as a numba-accelerated Python function that is natively multi-threaded, avoids cache misses, more efficiently caches intermediate values, and uses approximations at compute bottlenecks. Here, we provide a detailed description of the original Tomtom method and present results demonstrating that our re-implementation can achieve over a 1000-fold speedup compared with the original tool on reasonable tasks. AVAILABILITY AND IMPLEMENTATION: Our implementation of Tomtom is freely available as a Python package at https://github.com/jmschrei/memesuite-lite, which can be downloaded via pip install memelite or at https://zenodo.org/records/17008952.

Software

DNAvi: integration, statistics, and visualization of cell-free DNA fragment traces.

SUMMARY: DNAvi is a Python-based tool for rapid grouped analysis and visualization of cell-free DNA fragment size profiles directly from electrophoresis data, overcoming the need for sequencing in basic fragmentomic screenings. It enables normalization, statistical comparison, and publication-ready plotting of multiple samples, supporting quality control and exploratory fragmentomics in clinical and research workflows. AVAILABILITY AND IMPLEMENTATION: DNAvi is implemented in Python and freely available on GitHub at https://github.com/anjahess/DNAvi under a GNU General Public License v3.0, along with source code, documentation, and examples. An archived version is available under https://doi.org/10.5281/zenodo.18401705.

Software

ORFannotate: reproducible coding sequence annotation of transcriptome assemblies.

SUMMARY: Accurate annotation of coding sequences and translational features within transcript models is essential for interpreting assembled transcriptomes and their functional potential. Existing open reading frame (ORF) prediction tools typically operate on transcript FASTA files and do not reintegrate coding sequence (CDS) information back into transcript models, limiting their utility in long-read sequencing workflows where GTF/GFF annotations are the primary output. We present ORFannotate, a lightweight, GTF-native Python command-line tool that predicts ORFs from transcript annotations and reinserts precise, exon-aware CDS and UTR features into the original GTF/GFF file. In addition, ORFannotate provides biologically informative translational context by annotating Kozak sequence strength, detecting non-overlapping upstream ORFs (uORFs) with coding probabilities, characterising 5' and 3' untranslated regions (UTRs), and predicting nonsense-mediated decay (NMD) susceptibility. All annotations are consolidated in a transcript-level summary to support downstream analysis. By generating GTF files with accurate CDS annotations, ORFannotate facilitates reproducible analysis of both long- and short-read transcriptomes and integrates seamlessly with visualization tools, genome browsers, and comparative transcript analysis workflows. ORFannotate is fast, scalable and provides a practical solution for transcriptome annotation beyond coding potential prediction alone. AVAILABILITY AND IMPLEMENTATION: ORFannotate is implemented in Python and freely available under the GNU General Public License v3 (GPL-3.0) at: https://github.com/egustavsson/ORFannotate (DOI: https://doi.org/10.5281/zenodo.16812866).

Open Reading Frames

ssHiCstuff: a package for the design and analysis of ssDNA-specific Hi-C experiments.

MOTIVATION: Single-strand DNA-specific Hi-C (ssHi-C) is a recently developed technique enabling the capture of chromatin interactions involving single-stranded DNA (ssDNA), an intermediate of various DNA metabolic processes. ssHi-C entails the restoration of restriction sites in ssDNA regions of interest upon introduction of designer, internally barcoded "annealing oligonucleotides" prior to the restriction digestion step of Hi-C. The design of these "annealing oligonucleotides," as well as the analysis of the resulting ssHi-C data presents specific challenges, such as (i) differentiating ssDNA from dsDNA-derived contacts, (ii) tracking probe-specific interactions, and (iii) calibrating the amount of ssDNA contacts across biological samples. Dedicated computational tools are therefore needed to facilitate the design of, and extract biological information from, ssHi-C experiments. RESULTS: We present ssHiCstuff, a Rust- and Python-based package for the design of key reagents for ssHi-C experiments and for the analysis of ssHi-C data. ssHiCstuff provides (i) an automated annealing oligonucleotides design module, (ii) an end-to-end analyses pipeline, and (iii) a graphical user interface. ssHiCstuff simplifies the high-resolution analysis of ssDNA interactions at genome-wide scale. A graphical user interface (GUI) implemented in Python is also available for biologists without coding skills. AVAILABILITY: ssHiCstuff is freely available at https://github.com/Piazzalab/ssHiCstuff and https://zenodo.org/records/19677479 (https://doi.org/10.5281/zenodo.19677479) under the GPL 3.0 license. The annealing oligonucleotides design and the visualization modules are additionally freely available on a web browser at https://bioshiny.ens-lyon.fr/public/app/sshicstuff. A test dataset is available at https://zenodo.org/records/20035366 (https://doi.org/10.5281/zenodo.20035366).

DNA, Single-Stranded

Plotgardener App: a graphical interface for publication-ready genomic visualization.

SUMMARY: Plotgardener is an R package used for generating high-quality genomic visualizations. Despite its broad range of functions and versatility, its reliance on code presents a barrier for many potential users. To address this, we developed a macOS desktop application version of Plotgardener that enables users to create publication-ready genomic plots with no programming experience. The application employs a modular architecture comprising an Electron.js backend, a React frontend, and a Python parser that dynamically analyzes the Plotgardener package to ensure interface fields remain synchronized with package updates. By lowering the technical barrier to advanced genomic visualization, the Plotgardener desktop application broadens access to powerful visualization workflows for researchers and clinicians. AVAILABILITY: The current release of the Plotgardener App is an open source macOS desktop application built with Electron (Node.js), featuring a React frontend and a Python-based parser. The download link is available at https://phanstiellab.github.io/plotgardener/articles/guides/plotgardenerApp.html and on Zenodo (doi: https://doi.org/10.5281/zenodo.21684228). The source code is hosted on GitHub at https://github.com/rishabhsvemuri/ThePlotgardenerApp.

Genomics

MyESL: A Software for Evolutionary Sparse Learning in Molecular Phylogenetics and Genomics.

Evolutionary sparse learning uses supervised machine learning to build evolutionary models where genomic sites loci are parameters. It uses the Least Absolute Shrinkage and Selection Operator with bi-level sparsity to connect a specific phylogenetic hypothesis with sequence variation across genomic loci. The MyESL software addresses the need for open-source tools to perform evolutionary sparse learning analyses, offering features to preprocess input phylogenomic alignments, post-process output models to generate molecular evolutionary metrics, and make Least Absolute Shrinkage and Selection Operator regression adaptable and efficient for phylogenetic trees and alignments. The core of MyESL, which constructs models with logistic regressions using bi-level sparsity, is written in C++. Its input data preprocessing and result post-processing tools are developed in Python. Compared to other tools, MyESL is more computationally efficient and provides evolution-friendly inputs and outputs. These features have already enabled the use of MyESL in two phylogenomic applications, one to identify outlier sequences and fragile clades in inferred phylogenies and another to build genetic models of convergent traits. In addition to the use in a Python environment, MyESL is available as a standalone executable compatible across multiple platforms, which can be directly integrated into scripts and third-party software. The source code, executable, and documentation for MyESL are openly accessible at https://github.com/kumarlabgit/MyESL.

Phylogeny

[A histological study of striated muscle and nerve fiber in the human caused by a high voltage electric current].

A case of an electrocution is reported, involving a 37-year-old man who while engaged in an electrical work, received an electric shock by touching an electric wire of 11,000 volts. Both of his arms had to be amputated. A piece of muscle tissue was removed from the amputated arm to be examined by light and electron microscopy. Under light microscopy, the muscle tissue, after myoglobin staining, showed a particular pattern that closely resembled the pattern on the skin of a python. Thus the author has named the pattern seen in this case as the "Python-pattern". The electron microscopy inspection revealed on myoglobin staining excessive stretching and fragmentation of the myofibrils where the myoglobin staining was negative in the slide.

Adult

Death caused by a constricting snake--an infant death.

A 21-month-old male was found dead in his crib. An escaped pet python snake (Python reticulatus) was found near the child on a cribside shelf. Autopsy findings included petechiae associated with asphyxial death plus puncture wounds consistent with the reptile's teeth. This report details this 1982 death of a child in Reno, Nevada, and includes illustration of the reptile's bite marks plus description of snake jaw mechanics in swallowing.

Animals

Disease patterns in the Detroit Zoo: a study of reptilian and amphibian populations from 1973 through 1983.

A retrospective study was conducted to determine disease patterns in reptilian and amphibian populations at the Detroit Zoo from 1973 through 1983. In the reptilian population (mean +/- SD = 285.2 +/- 28), overall annual mortality rates were 1% to 40%. Mortality rates were highest in the fall months (20%) and lowest in the winter months (6%). The most frequently affected reptiles were iguana (Iguana iguana), reticulated python (Python reticulatus), rattlesnakes (Crotalus spp), common boa (Constrictor constrictor), and lizards (various genera of suborder Lacertilia). Of the 1,300 reptilian deaths from 1973 through 1983, 36.6% were caused by microbial agents, 12% by parasites, 11.6% by trauma, and 9.3% by nutritional deficiencies. The main microbial organisms that caused death were Aeromonas spp, Salmonella spp, Pseudomonas spp, Proteus spp, and Edwardsiella spp. The main parasites that caused death were Entamoeba spp and lungworms. Among amphibians, frogs and toads were the most frequently affected, and starvation and trauma were the most frequent causes of death.

Amphibians

Plasma concentrations of chloramphenicol in snakes.

Plasma chloramphenicol concentrations after a subcutaneous injection were studied in 87 snakes of 16 different species. The biological half-life of chloramphenicol varied from 3.3 hours in the indigo snake (Drymarchon corais couperi) to 22.1 hours in the midland water snake (Nerodia sipedon). A single dosage of 50 mg of chloramphenicol/kg of body weight produced plasma concentrations greater than 5 micrograms/ml for nearly 72 hours in 2 species of water snakes (Nerodia erythrogaster, Nerodia sipedon), for 24 hours in the Burmese python (Python molurus bivittatus), and for less than 12 hours in the gray rat snake, Indigo snake, and eastern king snake (Elaphe obsoleta spiloides, Drymarchon coraise couperi, and Lampropeltis getulus getulus). A dosage of 50 mg/kg administered to water snakes every 72 hours for 18 days maintained a minimum plasma concentration of chloramphenicol between 2 and 5 micrograms/ml.

Animals