Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “genome assembly validation”

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 73 records · Page 4Linked to original sources
PubMed2025

Whole-genome automated assembly pipeline for Chlamydia trachomatis strains from reference, in vitro and clinical samples using the integrated CtGAP pipeline.

Whole genome sequencing (WGS) is pivotal for the molecular characterization of Chlamydia trachomatis (Ct)-the leading bacterial cause of sexually transmitted infections and infectious blindness worldwide. Ct WGS can inform epidemiologic, public health and outbreak investigations of these human-restricted pathogens. However, challenges persist in generating high-quality genomes for downstream analyses given its obligate intracellular nature and difficulty with in vitro propagation. No single tool exists for the entirety of Ct genome assembly, necessitating the adaptation of multiple programs with varying success. Compounding this issue is the absence of reliable Ct reference strain genomes. We, therefore, developed CtGAP-Chlamydia trachomatisGenome Assembly Pipeline-as an integrated 'one-stop-shop' pipeline for assembly and characterization of Ct genome sequencing data from various sources including isolates, in vitro samples, clinical swabs and urine. CtGAP, written in Snakemake, enables read quality statistics output, adapter and quality trimming, host read removal, de novo and reference-guided assembly, contig scaffolding, selective ompA, multi-locus-sequence and plasmid typing, phylogenetic tree construction, and recombinant genome identification. Twenty Ct reference genomes were also generated. Successfully validated on a diverse collection of 363 samples containing Ct, CtGAP represents a novel pipeline requiring minimal bioinformatics expertise with easy adaptation for use with other bacterial species.

Chlamydia trachomatis↗
PubMed2026

Nerpa 2: probabilistic linking of biosynthetic gene clusters to nonribosomal peptides.

MOTIVATION: Nonribosomal peptides (NRPs) are bioactive microbial metabolites with high pharmaceutical potential. Although genome mining enables large-scale detection of biosynthetic gene clusters (BGCs) predicted to encode NRPs, reliably linking these clusters to their chemical products remains challenging due to the flexible and heterogeneous organization of NRP assembly pathways. RESULTS: We present Nerpa 2, a probabilistic framework for accurate and scalable linking of NRP BGCs to candidate chemical structures. The method represents assembly lines as hidden Markov models (HMMs) that capture uncertainty and alternative biosynthetic routes. On curated datasets of experimentally validated BGC-product pairs, our tool outperforms existing methods in linking accuracy and pathway reconstruction. When applied to large genome mining datasets, Nerpa 2 efficiently identifies BGCs likely associated with known compounds and highlights potential producers of novel chemistry. AVAILABILITY AND IMPLEMENTATION: Nerpa 2 is freely available at https://github.com/gurevichlab/nerpa.

Multigene Family↗
PubMed2026

Ecological Restoration of the Soil-Like Function in the Bauxite Residue: Natural Microbiomes Mediated Molecular Transformation of Dissolved Organic Matter.

Soilization of bauxite residues offers a scalable route for long-term carbon management and ecological restoration. However, the microbial processes that transform exogenous organic inputs into stable soil-like carbon pools remain poorly resolved. Here, we combined cross-ecosystem meta-analysis, machine-learning prediction, native synthetic community (SynCom) construction, 13C-labeled straw microcosms, field validation, Fourier transform ion cyclotron resonance mass spectrometry, and genome-resolved metagenomics to unravel microbiome-mediated carbon transformation at the dissolved organic matter (DOM) molecular scale. Our meta-analysis revealed that alkaline industrial wastes retained soil-like DOM signatures but were enriched in microbial humic- and protein-like components, indicating active yet incomplete carbon processing. Guided by these patterns, native SynCom inoculation increased 13C incorporation into total organic carbon (TOC) and dissolved organic carbon (DOC), enlarged biodegradable and adsorbable DOC fractions, and shifted DOM from recalcitrant aromatic pools toward oxygenated carbohydrate-, tannin-, and phenolic-like molecular classes. Genome-resolved analyses linked this transformation to complementary polymer degradation and nutrient-cycling functions across fungal and bacterial guilds, including enriched carbohydrate-active enzymes in straw-carbon-utilizing metagenome-assembled genomes. Null model and thermodynamic analyses further showed that microbial communities were constrained by homogeneous selection, whereas DOM molecules were diversified through variable selection and redox-dependent transformation. Field-scale validation confirmed that SynCom promoted TOC and DOC accumulation and humic-like, high-density DOM fractions under alkaline conditions. Together, these findings establish a mechanistic framework in which functional microbiomes couple plant carbon depolymerization, DOM molecular diversification, and mineral-interactive carbon stabilization, providing a microbiome-guided strategy for carbon sequestration and soilization in the bauxite residue.

Soil↗
PubMed2001

A BAC-based STS-content map spanning a 35-Mb region of human chromosome 1p35-p36.

We have devised a mapping method for rapid assembly and ordering of bacterial artificial chromosome (BAC) clones on a radiation hybrid (RH) panel, using sequence-tagged sites (STSs) and PCR. The protocol consists of two rounds of two-dimensional screening from a limited number of BACs to correspond each to an STS. In the first round, STSs are assembled in the RH bins and ordered according to PCR signals derived from 384-well microtiter plates (MTPs) in which BAC clones have been arrayed. In the second round, individual BAC clones are isolated from the MTPs to build a contig. We applied this method to a 35-Mb region spanning human chromosome 1p35-p36 and assembled 1366 BACs in 11 contigs, the longest being about 20 Mb. The working draft sequences of the human genome have been integrated into the contigs to validate the accuracy.

Chromosome Mapping↗
PubMed2005

Sequencing and comparing whole mitochondrial genomes of animals.

Comparing complete animal mitochondrial genome sequences is becoming increasingly common for phylogenetic reconstruction and as a model for genome evolution. Not only are they much more informative than shorter sequences of individual genes for inferring evolutionary relatedness, but these data also provide sets of genome-level characters, such as the relative arrangements of genes, which can be especially powerful. We describe here the protocols commonly used for physically isolating mitochondrial DNA (mtDNA), for amplifying these by polymerase chain reaction (PCR) or rolling circle amplification (RCA), for cloning, sequencing, assembly, validation, and gene annotation, and for comparing both sequences and gene arrangements. On several topics, we offer general observations based on our experiences with determining and comparing complete mitochondrial DNA sequences.

Animals↗
PubMed2004

Shotgun optical mapping of the entire Leishmania major Friedlin genome.

Leishmania is a group of protozoan parasites which causes a broad spectrum of diseases resulting in widespread human suffering and death, as well as economic loss from the infection of some domestic animals and wildlife. To further understand the fundamental genomic architecture of this parasite, and to accelerate the on-going sequencing project, a whole-genome XbaI restriction map was constructed using the optical mapping system. This map supplemented traditional physical maps that were generated by fingerprinting and hybridization of cosmid and P1 clone libraries. Thirty-six optical map contigs were constructed for the corresponding known 36 chromosomes of the Leishmania major Friedlin genome. The chromosome sizes ranged from 326.9 to 2821.3 kb, with a total genome size of 34.7 Mb; the average XbaI restriction fragment was 25.3 kb, and ranged from 15.7 to 77.8 kb on a per chromosomes basis. Comparison between the optical maps and the in silico maps of sequence drawn from completed, nearly finished, or large sequence contigs showed that optical maps served several useful functions within the path to create finished sequence by: guiding aspects of the sequence assembly, identifying misassemblies, detection of cosmid or PAC clones misplacements to chromosomes, and validation of sequence stemming from varying degrees of finishing. Our results also showed the potential use of optical maps as a means to detect and characterize map segmental duplication within genomes.

Animals↗
PubMed2025

Genome size estimation from long read overlaps.

MOTIVATION: Accurate genome size estimation is an important component of genomic analyses such as assembly and coverage calculation, though existing tools are primarily optimized for short-read data. RESULTS: We present LRGE, a novel tool that uses read-to-read overlap information to estimate genome size in a reference-free manner. LRGE calculates per-read genome size estimates by analysing the expected number of overlaps for each read, considering read lengths and a minimum overlap threshold. The final size is taken as the median of these estimates, ensuring robustness to outliers such as reads with no overlaps. Additionally, LRGE provides an expected confidence range for the estimate. We validate LRGE on a large, diverse bacterial dataset and confirm it generalizes to eukaryotic datasets. On bacterial genomes, LRGE outperforms k-mer-based methods in both accuracy and computational efficiency and produces genome size estimates comparable to those from assembly-based approaches, like Raven, while using significantly less computational resources. AVAILABILITY AND IMPLEMENTATION: Our method, LRGE (Long Read-based Genome size Estimation from overlaps), is implemented in Rust and is available as a precompiled binary for most architectures, a Bioconda package, a prebuilt container image, and a crates.io package as a binary (lrge) or library (liblrge). The source code is available at https://github.com/mbhall88/lrge and an archive at https://doi.org/10.5281/zenodo.17183812 under an MIT license.

Genome Size↗
PubMed2004

Large-scale validation of single nucleotide polymorphisms in gene regions.

Genome-wide association studies using large numbers of bi-allelic single nucleotide polymorphisms (SNPs) have been proposed as a potentially powerful method for identifying genes involved in common diseases. To assemble a SNP collection appropriate for large-scale association, we designed assays for 226,099 publicly available SNPs located primarily within known and predicted gene regions. Allele frequencies were estimated in a sample of 92 CEPH Caucasians using chip-based MALDI-TOF mass spectrometry with pooled DNA. Of the 204,200 designed assays that were functional, 125,799 SNPs were determined to be polymorphic (minor allele frequency > 0.02), of which 101,729 map uniquely to the human genome. Many of the commonly available RefSNP annotations were predictive of polymorphic status and could be used to improve the selection of SNPs from the public domain for genetic research. The set of uniquely mapping, polymorphic SNPs is located within 10 kb of 66% of known and predicted genes annotated in LocusLink, which could prove useful for large-scale disease association studies.

Databases, Genetic↗
PubMed2001

Integration of cytogenetic landmarks into the draft sequence of the human genome.

We have placed 7,600 cytogenetically defined landmarks on the draft sequence of the human genome to help with the characterization of genes altered by gross chromosomal aberrations that cause human disease. The landmarks are large-insert clones mapped to chromosome bands by fluorescence in situ hybridization. Each clone contains a sequence tag that is positioned on the genomic sequence. This genome-wide set of sequence-anchored clones allows structural and functional analyses of the genome. This resource represents the first comprehensive integration of cytogenetic, radiation hybrid, linkage and sequence maps of the human genome; provides an independent validation of the sequence map and framework for contig order and orientation; surveys the genome for large-scale duplications, which are likely to require special attention during sequence assembly; and allows a stringent assessment of sequence differences between the dark and light bands of chromosomes. It also provides insight into large-scale chromatin structure and the evolution of chromosomes and gene families and will accelerate our understanding of the molecular bases of human disease and cancer.

Chromosome Aberrations↗
PubMed2025

Infectious Subgenomic Amplicon Strategies for Japanese Encephalitis and West Nile Viruses.

Classical methods for constructing infectious cDNA clones of flaviviruses are often hindered by instability and toxicity. The Infectious-Subgenomic-Amplicons (ISA) method is an advancement which utilizes overlapping DNA fragments representing viral genomic sequence and in-cell recombination to bypass bacterial plasmid assembly. However, the ISA method has limitations due to the toxicity of some ISA DNA fragments in bacteria during synthetic production. We validated modified ISA strategies for producing toxic ISA Japanese encephalitis virus (JEV) and West Nile virus (WNV) DNA fragments. Three approaches were explored, including subdividing toxic DNA fragments into two sub-fragments for synthetic clonal production, using a low-copy bacterial plasmid, and subdividing the toxic DNA fragments into four short overlapping sub-fragments, each up to 1.8 kb. The latter novel approach in ISA applications enabled the synthesis of entirely bacteria-free ISA DNA fragments. Our results demonstrate that subdividing toxic fragments into sub-fragments smaller than 1.8 kb for synthesis is the efficient strategy, circumventing the need for bacterial plasmids and ensuring rapid production of synthetic flaviviruses. This method also shortens the production timeline. We also compared the efficacy of JEV and WNV ISA in zinc finger antiviral protein 1 (ZAP) wild-type and knockout cells and found that knockout cells may be more effective for ISA rescue of flaviviruses, including CpG-enriched strains for live attenuated vaccines. The validated modified ISA strategies provide an efficient approach for producing synthetic JEV and WNV. This will enable rapid research during outbreaks of emerging flaviviruses by facilitating the quick generation of new virus variants.

West Nile virus↗
PubMed2025

DeepES: deep learning-based enzyme screening to identify orphan enzyme genes.

MOTIVATION: Progress in sequencing technology has led to determination of large numbers of protein sequences, and large enzyme databases are now available. Although many computational tools for enzyme annotation were developed, sequence information is unavailable for many enzymes, known as orphan enzymes. These orphan enzymes hinder sequence similarity-based functional annotation, leading gaps in understanding the association between sequences and enzymatic reactions. RESULTS: Therefore, we developed DeepES, a deep learning-based tool for enzyme screening to identify orphan enzyme genes, focusing on biosynthetic gene clusters and reaction class. DeepES uses protein sequences as inputs and evaluates whether the input genes contain biosynthetic gene clusters of interest by integrating the outputs of the binary classifier for each reaction class. The validation results suggested that DeepES can capture functional similarity between protein sequences, and it can be implemented to explore orphan enzyme genes. By applying DeepES to 4744 metagenome-assembled genomes, we identified candidate genes for 236 orphan enzymes, including those involved in short-chain fatty acid production as a characteristic pathway in human gut bacteria. AVAILABILITY AND IMPLEMENTATION: DeepES is available at https://github.com/yamada-lab/DeepES. Model weights and the candidate genes are available at Zenodo (https://doi.org/10.5281/zenodo.11123900).

Deep Learning↗
PubMed2004

ApiEST-DB: analyzing clustered EST data of the apicomplexan parasites.

ApiEST-DB (http://www.cbil.upenn.edu/paradbs-servlet/) provides integrated access to publicly available EST data from protozoan parasites in the phylum Apicomplexa. The database currently incorporates a total of nearly 100,000 ESTs from several parasite species of clinical and/or veterinary interest, including Eimeria tenella, Neospora caninum, Plasmodium falciparum, Sarcocystis neurona and Toxoplasma gondii. To facilitate analysis of these data, EST sequences were clustered and assembled to form consensus sequences for each organism, and these assemblies were then subjected to automated annotation via similarity searches against protein and domain databases. The underlying relational database infrastructure, Genomics Unified Schema (GUS), enables complex biologically based queries, facilitating validation of gene models, identification of alternative splicing, detection of single nucleotide polymorphisms, identification of stage-specific genes and recognition of phylogenetically conserved and phylogenetically restricted sequences.

Animals↗
PubMed2025

A chromosome-level reference genome assembly of the Small snakehead (Channa asiatica).

The Small snakehead (Channa asiatica) is an economically important species in both aquaculture and ornamental trade, mainly distributed in South China and Southeast Asia. Despite its significance, limited genomic resources have impeded in-depth genetic studies and breeding programs. In this study, we used PacBio HiFi long-read sequencing, Illumina short-read sequencing, and Hi-C technologies to generate a high-quality chromosome-level genome of the C. asiatica. The final genome spans 659.44 Mb, with an impressive 98.18% anchored to 23 chromosomes. Notably, the contig N50 and scaffold N50 are 23.92 Mb and 29.61 Mb, validated by a BUSCO completeness score of 98.93%. Genome annotation identified 26,603 protein-coding genes, 99.29% of which were confirmed by BUSCO analysis, and 93.68% were functionally annotated. Approximately 27.72% of the genome sequences were classified as repeat elements. This high-fidelity genome assembly provides a robust foundation for advancing molecular breeding, comparative genomics, and evolutionary studies of C. asiatica and related species.

Animals↗
PubMed2006

Genome assembly comparison identifies structural variants in the human genome.

Numerous types of DNA variation exist, ranging from SNPs to larger structural alterations such as copy number variants (CNVs) and inversions. Alignment of DNA sequence from different sources has been used to identify SNPs and intermediate-sized variants (ISVs). However, only a small proportion of total heterogeneity is characterized, and little is known of the characteristics of most smaller-sized (<50 kb) variants. Here we show that genome assembly comparison is a robust approach for identification of all classes of genetic variation. Through comparison of two human assemblies (Celera's R27c compilation and the Build 35 reference sequence), we identified megabases of sequence (in the form of 13,534 putative non-SNP events) that were absent, inverted or polymorphic in one assembly. Database comparison and laboratory experimentation further demonstrated overlap or validation for 240 variable regions and confirmed >1.5 million SNPs. Some differences were simple insertions and deletions, but in regions containing CNVs, segmental duplication and repetitive DNA, they were more complex. Our results uncover substantial undescribed variation in humans, highlighting the need for comprehensive annotation strategies to fully interpret genome scanning and personalized sequencing projects.

Base Sequence↗
PubMed2002

Finishing a whole-genome shotgun: release 3 of the Drosophila melanogaster euchromatic genome sequence.

BACKGROUND: The Drosophila melanogaster genome was the first metazoan genome to have been sequenced by the whole-genome shotgun (WGS) method. Two issues relating to this achievement were widely debated in the genomics community: how correct is the sequence with respect to base-pair (bp) accuracy and frequency of assembly errors? And, how difficult is it to bring a WGS sequence to the accepted standard for finished sequence? We are now in a position to answer these questions. RESULTS: Our finishing process was designed to close gaps, improve sequence quality and validate the assembly. Sequence traces derived from the WGS and draft sequencing of individual bacterial artificial chromosomes (BACs) were assembled into BAC-sized segments. These segments were brought to high quality, and then joined to constitute the sequence of each chromosome arm. Overall assembly was verified by comparison to a physical map of fingerprinted BAC clones. In the current version of the 116.9 Mb euchromatic genome, called Release 3, the six euchromatic chromosome arms are represented by 13 scaffolds with a total of 37 sequence gaps. We compared Release 3 to Release 2; in autosomal regions of unique sequence, the error rate of Release 2 was one in 20,000 bp. CONCLUSIONS: The WGS strategy can efficiently produce a high-quality sequence of a metazoan genome while generating the reagents required for sequence finishing. However, the initial method of repeat assembly was flawed. The sequence we report here, Release 3, is a reliable resource for molecular genetic experimentation and computational analysis.

Animals↗
PubMed2004

Bottlenecks and roadblocks in high-throughput XAS for structural genomics.

Structural and functional characterization of the entire protein complement (the proteome) of an organism can provide an infrastructure upon which questions about biological pathways and systems biology can be framed. The technology necessary to perform this proteome-level structural and functional characterization is under development in numerous structural genomics and functional genomics initiatives. Given the ubiquity of metal active sites in a proteome, it seems appropriate to ask whether comprehensive local structural characterization of metal sites within a proteome (metalloproteomics) is either a valid or obtainable goal. With a proteome-wide knowledge of the active-site structures of all metalloproteins, one could start to ask how metal insertion, cluster assembly and metalloprotein expression are affected by growth conditions or developmental status etc. High-throughput X-ray absorption spectroscopy (HTXAS) is being developed as a technology for investigating the metalloproteome. In creating a pipeline from genome to metalloproteome, several bottlenecks to high-throughput determination of metal-site structures must be overcome. For example, automation of arraying small samples for XAS examination must be invented, automation of rapid data collection of multiple low-volume low-concentration samples must be developed, automation of data reduction and analysis must be perfected. Discussed here are the promises and the pitfalls of HTXAS development, including the results of initial feasibility experiments.

Automation↗
PubMed1997

Construction of a 1-Mb restriction-mapped cosmid contig containing the candidate region for the familial Mediterranean fever locus (MEFV) on chromosome 16p 13.3.

In this paper we describe the assembly and restriction map of a 1.05-Mb cosmid contig spanning the candidate region for familial Mediterranean fever (FMF), a recessively inherited disorder of inflammation localized to 16p13.3. Using a combination of cosmid walking and screening for P1, PAC, BAC, and YAC clones, we have generated a contig of genomic clones spanning approximately 1050 kb that contains the FMF critical region. The map consists of 179 cosmid, 15 P1, 10 PAC, 3 BAC, and 17 YAC clones, anchored by 27 STS markers. Eight additional STSs have been developed from the approximately 700 kb immediately centromeric to this genomic region. Five of the 35 STSs are microsatellites that have not been previously reported. NotI and EcoRI mapping of the overlapping cosmids, hybridization of restriction fragments from cosmids to one another, and STS analyses have been used to validate the assembly of the contig. Our contig totally subsumes the 250-kb interval recently reported, by founder haplotype analysis, to contain the FMF gene. Thus, our high-resolution clone map provides an ideal resource for transcriptional mapping toward the eventual identification of this disease gene.

Base Sequence↗
PubMed2004

Genomic annotation of the meningioma tumor suppressor locus on chromosome 1p34.

Meningioma is a frequently occurring tumor of the meninges surrounding the central nervous system. Loss of the short arm of chromosome 1 (1p) is the second most frequent chromosomal abnormality observed in these tumors. Previously, we identified a 3.7 megabase (Mb) region of consistent deletion on 1p33-p34 in a panel of 157 tumors. Loss of this region was associated with advanced disease and predictive for tumor relapse. In this report, a high-resolution integrated map of the region was constructed (CompView) to identify all markers in the smallest region of overlapping deletion (SRO). A regional somatic cell hybrid panel was used to more precisely localize those markers identified in CompView as within or overlapping the region. Additional deletion mapping using microsatellites localized to the region narrowed the SRO to approximately 2.8 Mb. The 88 markers remaining in the SRO were used to screen genomic databases to identify large-insert clones. Clones were assembled into a physical map of the region by PCR-based, sequence-tagged site (STS) content mapping. A sequence from clones was used to validate STS content by electronic PCR and to identify transcripts. A minimal tiling path of 43 clones was constructed across the SRO. Sequence data from the most current sequence assembly were used for further validation. A total of 59 genes were ordered within the SRO. In all, 17 of these were selected as likely candidates based on annotation using Gene Ontology Consortium terms, including the MUTYH, PRDX1, FOXD2, FOXE3, PTCH2, and RAD54L genes. This annotation of a putative tumor suppressor locus provides a resource for further analysis of meningioma candidate genes.

Chromosomes, Human, Pair 1↗