Biomedical subjects

Genomics

Find indexed PubMed genomics citations. Search gene expression, sequencing and genetic variation in titles, abstracts and supplied subjects, then open the PubMed record.

Follow this topic with RSS

At least 1,747 records · Page 97Linked to original sources

PubMed2002

On the sequencing of the human genome.

Two recent papers using different approaches reported draft sequences of the human genome. The international Human Genome Project (HGP) used the hierarchical shotgun approach, whereas Celera Genomics adopted the whole-genome shotgun (WGS) approach. Here, we analyze whether the latter paper provides a meaningful test of the WGS approach on a mammalian genome. In the Celera paper, the authors did not analyze their own WGS data. Instead, they decomposed the HGP's assembled sequence into a "perfect tiling path", combined it with their WGS data, and assembled the merged data set. To study the implications of this approach, we perform computational analysis and find that a perfect tiling path with 2-fold coverage is sufficient to recover virtually the entirety of a genome assembly. We also examine the manner in which the assembly was anchored to the human genome and conclude that the process primarily depended on the HGP's sequence-tagged site maps, BAC maps, and clone-based sequences. Our analysis indicates that the Celera paper provides neither a meaningful test of the WGS approach nor an independent sequence of the human genome. Our analysis does not imply that a WGS approach could not be successfully applied to assemble a draft sequence of a large mammalian genome, but merely that the Celera paper does not provide such evidence.

Chromosomes, Artificial, Bacterial↗

PubMed2026

Seqwin: ultrafast identification of signature sequences in microbial genomes.

MOTIVATION: Polymerase chain reaction (PCR) enables rapid, cost-effective diagnostics but requires prior identification of genomic regions that allow sensitive and specific detection of target microbial groups, herein referred to as microbial signature sequences. We introduce Seqwin, an open-source framework designed to automate microbial genome signature discovery. Tens of thousands of microbial genomes are now available for a single species, limiting the application of existing manual and automated approaches for identifying signatures. Modern approaches that are capable of leveraging all available microbial genomes will ensure sensitive and accurate DNA signature identification and enable robust pathogen detection for clinical, environmental, and public health applications. RESULTS: Seqwin builds weighted pan-genome minimizer graphs and uses a traversal algorithm to identify signature sequences that occur frequently in target genomes but remain rare in non-targets. Unlike earlier tools that depend on strict presence or absence of sequences, Seqwin accommodates natural sequence variation and scales to very large genome collections. When applied to genomes from C. difficile, M. tuberculosis, and S. enterica, Seqwin recovered more high-quality signatures than alternative methods with lower computational burden. Seqwin's analysis of nearly 15 000 S. enterica genomes yielded over 200 candidate signatures in three minutes. Seqwin provides an open-source solution for the long-standing need for scalable microbial signature discovery and diagnostic assay design. AVAILABILITY AND IMPLEMENTATION: Seqwin is available on GitHub (https://github.com/treangenlab/Seqwin) and can be installed via Bioconda (https://bioconda.github.io/recipes/seqwin/README.html). Benchmarking datasets, outputs, and scripts are available on Zenodo (https://doi.org/10.5281/zenodo.19874011).

Software↗

PubMed2005

GenColors: accelerated comparative analysis and annotation of prokaryotic genomes at various stages of completeness.

SUMMARY: GenColors is a new web-based software/database system aimed at an improved and accelerated annotation of prokaryotic genomes, considering information on related genomes and making extensive use of genome comparison. It offers a seamless integration of data from ongoing sequencing projects and annotated genomic sequences obtained from GenBank. The genome comparison tools determine, for example, best-bidirectional hits, gene conservation, syntenies and gene core sets. Swiss-Prot/TrEMBL hits allow annotations in an effective manner. To further support the annotation base-specific quality data can also be displayed if available. With GenColors dedicated genome browsers containing a group of related genomes can be easily set up and maintained. It has been efficiently used for Borrelia garinii and is currently applied to various ongoing genome projects. AVAILABILITY: Detailed information on GenColors is available at http://gencolors.imb-jena.de. Online usage of GenColors-based genome browsers is the preferred application mode. The system is also available upon request for local installation.

Borrelia↗

PubMed2025

Annotated genome of the Atlantic dog whelk, Nucella lapillus.

Nucella lapillus is an important player in rocky shore food chains and has been a focal organism of ecological and evolutionary studies for decades. Despite poor dispersal, they have a broad geographic range, which makes them an ideal species to examine isolation by distance and selection across environmental gradients. Here we present the fully annotated genome of N. lapillus generated with Oxford Nanopore Techonology sequencing at ∼37× coverage. The genome assembly is 2.32 Gbp and consists of 2,525 contigs, with an N50 length of 2 Mbp. Repeat annotation identified 2,491 families that cover 67.56% of the genome, which is similar to other gastropods. Despite its large size and high proportion of repeats, the genome is of high quality. Benchmarking Universal Single-Copy Ortholog (BUSCO) analysis revealed a score of 96.8%. Functional annotation of the genome produced 45,848 protein-coding genes with a 96.6% BUSCO score. Genomic resources for mollusks lag behind that of other phyla, perhaps because many of their innate characteristics complicate DNA extraction, sequencing, and assembly. This new N. lapillus genome will increase our genomic understanding of the second largest phylum (and the most diverse class within said phylum) and serve as a key resource to advance studies on the organismal biology and population genetics of this iconic species as well as the connection between genomic variation and community-level processes.

Animals↗

PubMed2006

MBGD: a platform for microbial comparative genomics based on the automated construction of orthologous groups.

The microbial genome database for comparative analysis (MBGD) is a comprehensive platform for microbial comparative genomics. The central function of MBGD is to create orthologous groups among multiple genomes from precomputed all-against-all similarity relationships using the DomClust algorithm. The database now contains >300 published genomes and the number continues to grow. For researchers who are interested in ongoing genome projects, we have now started a new service called 'My MBGD,' which allows users to add their own genome sequences to MBGD for the purpose of identifying orthologs among both the new and the existing genomes. Furthermore, in order to make available the rapidly accumulating information on closely related genome sequences, we enhanced the interface for pairwise genome comparisons using the CGAT interface, which allows users to see nucleotide sequence alignments of non-coding as well as coding regions. MBGD is available at http://mbgd.genome.ad.jp/.

Algorithms↗

PubMed2026

Genomic surveillance of enterovirus D68 circulating in 2025 reveals the emergence of a novel A2/B3 recombinant lineage.

Enterovirus D68 (EV-D68) has re-emerged over the past decade as a significant respiratory pathogen associated with severe respiratory disease and acute flaccid myelitis. Its circulation has typically followed a biennial pattern, with predominance in late summer and early fall, a pattern that was temporarily disrupted during the COVID-19 pandemic. Surveillance in 2025 revealed off-season circulation of EV-D68. This study describes the genomic characteristics of the 2025 EV-D68 viruses and the clinical features of affected patients. Between May and December 2025, remnant respiratory specimens positive for rhinovirus/enterovirus were screened for EV-D68 and subjected to whole-genome sequencing. Phylogenetic analyses were performed using maximum-likelihood methods. Recombination was assessed using subgenomic phylogenies, SimPlot similarity and BootScan analyses, and read-level inspection. Among 1,321 patients tested, 147 (11.1%) were EV-D68-positive, and 119 (81.0%) yielded complete genomes. EV-D68 positivity increased in July 2025, peaked in August (~21%), and remained elevated through September and October, exceeding levels observed in 2024. Patients had a median age of 36 years, with infections disproportionately affecting older adults. Phylogenetic analysis demonstrated exclusive circulation of subclade A2. Five genomes formed a distinct recombinant lineage (A2-Re). Subgenomic phylogenies showed clustering with A2 viruses in the P1 region and with B3 viruses in the P2-P3 regions. SimPlot and BootScan analyses identified a recombination breakpoint near the 2A/2B junction (~nt 3,700). The recombinant lineage was associated with temporally clustered cases in September-October. These findings demonstrate recombination between distinct EV-D68 subclades and underscore the importance of whole-genome surveillance for accurate viral characterization. Continued genomic monitoring is essential for detecting emerging variants with potential implications for transmissibility, pathogenicity, and public health preparedness.IMPORTANCEThis study highlights an increased off-season circulation of Enterovirus D68 (EV-D68) and a higher burden of disease in adults in 2025. The identification of a novel A2-B3 recombinant lineage provides evidence of ongoing viral evolution through recombination, a mechanism that may alter transmissibility, virulence, or immune responses. Detection of this lineage in temporally clustered cases suggests local transmission and underscores the potential for rapid spread of newly emerged variants. These findings emphasize the limitations of partial genomic approaches and the critical role of whole-genome sequencing in accurately characterizing circulating strains and identifying recombination events. Enhanced genomic surveillance is essential to detect emerging variants in real time, inform diagnostic assay performance, and support public health responses. Continued monitoring of EV-D68 evolution will be important for anticipating changes in disease burden, guiding clinical awareness, and strengthening preparedness for future outbreaks.

Humans↗

PubMed2025

Comparative genomics of carbapenem-resistant Acinetobacter baumannii isolated from pediatric patients in a tertiary care hospital.

Acinetobacter baumannii is a short gram-negative bacillus, notable for its intrinsic multidrug resistance and genomic plasticity, which facilitates the acquisition of additional resistance genes via mobile genetic elements. Due to its increasing carbapenem resistance, the World Health Organization has classified it as a critical priority pathogen. This study performed a comparative genomic analysis of 20 carbapenem-resistant A. baumannii clinical strains isolated from the Hospital Infantil de México Federico Gómez (CRAB-HIMFG), alongside 11 genomes from other Mexican strains. The pangenome was determined to be open, and core genome single-nucleotide polymorphism-based analysis grouped the CRAB-HIMFG strains within CC758/IC5 and CC92/IC2. A novel sequence type (ST) in the MLST-Pasteur scheme was identified, related to STPas156, and in the MLST-Oxford scheme, associated with STOxf758 and STOxf1054. Virulence and resistance genes comprised 0.61% to 2.23% of the pangenome. Oxacillinase genes and efflux pumps primarily mediated carbapenem resistance, while virulence genes included those encoding biofilm and type IV pili. Capsule typing revealed a correlation with established international clones, IC2 and IC5. Plasmids exhibited high diversity, harboring maintenance modules and toxin-antitoxin systems, with the dissemination of resistance genes linked to insertion sequences. Biofilm formation and twitching motility were not always expressed, as they depend on additional environmental factors. Our study shows that comparative genomics is an essential tool to analyze clinically and epidemiologically significant genomes, providing critical insights into gene distribution, genomic architecture, and horizontal gene transfer mechanisms in microbial populations.IMPORTANCEIn recent years, a reported increase in the mortality rate associated with infections caused by A. baumannii, along with a rise in carbapenem resistance, poses a serious clinical challenge. The WHO considered this microorganism critical for research into alternative therapies and epidemiological surveillance. Despite advances in bioinformatics, genomic studies have yet to fully elucidate the structural rearrangements and secretion systems of A. baumannii. This knowledge gap hinders our understanding of its remarkable genomic plasticity and its ability to acquire and spread resistance and virulence genes through horizontal gene transfer.

Acinetobacter baumannii↗

PubMed2005

Visualization of comparative genomic analyses by BLAST score ratio.

BACKGROUND: The first microbial genome sequence, Haemophilus influenzae, was published in 1995. Since then, more than 400 microbial genome sequences have been completed or commenced. This massive influx of data provides the opportunity to obtain biological insights through comparative genomics. However few tools are available for this scale of comparative analysis. RESULTS: The BLAST Score Ratio (BSR) approach, implemented in a Perl script, classifies all putative peptides within three genomes using a measure of similarity based on the ratio of BLAST scores. The output of the BSR analysis enables global visualization of the degree of proteome similarity between all three genomes. Additional output enables the genomic synteny (conserved gene order) between each genome pair to be assessed. Furthermore, we extend this synteny analysis by overlaying BSR data as a color dimension, enabling visualization of the degree of similarity of the peptides being compared. CONCLUSIONS: Combining the degree of similarity, synteny and annotation will allow rapid identification of conserved genomic regions as well as a number of common genomic rearrangements such as insertions, deletions and inversions. The script and example visualizations are available at: http://www.microbialgenomics.org/BSR/.

Algorithms↗

PubMed2026

[Applications and Challenges of Deep Learning in Human Genome Research].

In recent years, the advent of high-throughput omics technologies has fueled an explosive growth in human genomic data. Uncovering the latent functions within this vast data has become a significant challenge in functional genomics research. While traditional statistical methods have proved successful for analyzing smaller-scale datasets in the past, they exhibit clear limitations in analytical efficiency and integrating multi-dimensional data, struggling to meet the escalating demands of contemporary genomic analysis. The introduction of deep learning (DL) technologies offers a novel paradigm for this field. This review systematically examines the advances in applying deep learning to human genomics research. Studies demonstrate that when ample labeled data is available, discriminative DL computational methods-such as Convolutional Neural Networks (CNNs) and Long Short-Term Memory networks (LSTMs)-achieve high accuracy and efficiency in genomic variant discovery tasks. Furthermore, generative DL methods, particularly Large Language Models (LLMs) leveraging self-supervised pre-training strategies, effectively integrate complex genomic information and exhibit superior performance in functional genomic sequence annotation and gene regulation studies. This review also explores the application of LLMs in multi-omics data integration and prediction. Looking ahead, the continued accumulation of long-read sequencing and high-dimensional data is expected to enable DL technologies to integrate increasingly complex and heterogeneous genomic information, playing an increasingly crucial role in human genomics research.

Deep Learning↗

PubMed2026

Whole-Genome Analysis Reveals Antimicrobial Resistance and Population Structure of Environmental and Veterinary Acinetobacter baumannii.

Acinetobacter (A.) baumannii is an important multidrug-resistant pathogen increasingly recognized across animal and environmental settings, and carbapenem-resistant A. baumannii (CRAB) is classified as a critical-priority pathogen by the World Health Organization. This study investigated the antimicrobial resistance (AMR) and genomic characteristics of 122 A. baumannii isolates comprising 72 veterinary and 50 environmental isolates collected in Andhra Pradesh, India. Antimicrobial susceptibility testing, whole-genome sequencing (WGS), resistance and virulence gene profiling, multilocus sequence typing (MLST), core-genome analysis, single nucleotide polymorphism (SNP) phylogeny, and pan-genome analysis were performed. Overall, 58.2% of isolates were multidrug-resistant (MDR), and 41.8% were extensively drug-resistant (XDR). Sequence type (ST) 52 predominated among veterinary isolates, whereas ST2 was more frequent among environmental isolates. The presence of carbapenem-resistant isolates along with the ST2 lineage enhances the similarity to clinical A. baumannii. Several intrinsic resistance genes, including blaOXA-23, armA, aph(3″)-Ib, aph(6)-Id, tet(B), mph(E), and msr(E), were more prevalent in the ST2-associated population. Virulence-associated determinants were widely conserved. Core-genome MLST (cgMLST) and core-genome SNP (cgSNP) analyses identified highly related isolates within both lineages, while pairwise SNP differences were 0-7. Pan-genome analysis identified 4204 gene clusters and distinct accessory gene patterns between ST2 and ST52. These findings indicate that resistance gene distribution was closely associated with lineage structure and support integrated genomic surveillance of A. baumannii across animal and environmental reservoirs.

Acinetobacter baumannii↗

PubMed2002

Genetic differentiation for nuclear, mitochondrial and chloroplast genomes in common wild rice ( Oryza rufipogon Griff.) and cultivated rice ( Oryza sativa L.).

The genetic differentiation of nuclear, mitochondrial (mt) and chloroplast (cp) genomes was investigated by Southern and PCR analysis using 75 varieties of cultivated rice ( Oryza sativa L.) and 118 strains of common wild rice (CWR, Oryza rufipogon Griff.) from ten countries of Asia. The distinguishing differences between the Indica and Japonica cultivars were detected both in the nuclear genome and the cytoplasmic genome, confirming that the Indica-Japonica differentiation is of major importance for the three different classes of genome in cultivated rice. This differentiation was also detected in common wild rice with some differences among the genome compartments and the various regions. For nuclear DNA variation, both Indica-like and Japonica-like types were observed in the Chinese CWR, with the latter more-frequent than the former. No Japonica-like type was found in South Asia, and only two strains of the Japonica-like type were detected in Southeast Asia, thus the Indica-like type is the major type among South and Southeast Asian CWR. For mtDNA, only a few strains of the Japonica-like type were detected in CWR. For cpDNA, the Japonica type was predominant among the CWR strains from China, Bangladesh and Burma, while the Indica type was predominant among the CWR strains from Thailand, Malaysia, Cambodia and Sri Lanka, and both types were found in similar frequencies among the Indian CWR. Altogether, however, the degree of Indica-Japonica differentiation in common wild rice was much-less important than that in cultivated rice. Cluster analyses for nuclear and mitochondrial DNA variation revealed that some CWR strains showed large genetic distances from cultivated rice and formed clusters distinct from cultivated rice. Coincidence in the genetic differentiation between the three different classes of genome was much higher in cultivated rice than in CWR. Among the 75 cultivars, about 3/4 entries were "homoeotype" showing congruent results for nuclear, mt and cpDNA regarding the Indica-Japonica differentiation. In CWR, the proportions of homoeotypes were 5.7%, 15% and 48.8% in China, South Asia and Southeast Asia, respectively. Based on the average genetic distance among all the strains of CWR and cultivated rice for nuclear and mitochondrial genomes, the variability of the nuclear genome was found to be higher than that of the mitochondrial genome. The global pattern based on all genomes shows much-more diversification in CWR than that in cultivated rice.

Journal Article↗

PubMed1983

Localization of the regions of homology between the genomes of herpes simplex virus, type 1, and pseudorabies virus.

Only 8% of the sequences of the genomes of pseudorabies (PRV) and herpes simplex (type 1) (HSV) viruses are homologous. These homologous sequences have been shown previously to be distributed throughout most of the genomes of the two viruses. By means of blot hybridization of restriction fragments of HSV-1 DNA to cloned, nick-translated restriction fragments of PRV DNA, it was possible to compare the location on the genomes of these viruses of the homologous regions. The results showed that the genome of PRV is, for the most part, colinear with the IL arrangement of the genome of HSV-1. An inversion or translocation of sequences mapping on the PRV genome between 0.07 and 0.39 map units was observed on the genome of one of these viruses. A comparison of the map positions of five genes with known functions confirmed these findings. The genes coding for the major immediate-early protein, the major capsid protein, and the thymidine kinase occupy similar positions on the genome of PRV and on the genome of HSV-1 in the IL arrangement. However, the genes for DNA polymerase and for the major DNA binding protein appear to be inverted relative to one another on the genomes of the two viruses.

Base Sequence↗

PubMed1990

Sequence of the genome RNA of rubella virus: evidence for genetic rearrangement during togavirus evolution.

The nucleotide sequence of the rubella virus (RUB) genomic RNA was determined. The RUB genomic RNA is 9757 nucleotides in length [excluding the poly(A) tail] and has a G/C content of 69.5%, the highest of any RNA virus sequenced to data. The RUB genomic RNA contains two long open reading frames (ORFs), a 5'-proximal ORF of 6656 nucleotides and a 3'-proximal ORF of 3189 nucleotides which encodes the structural proteins. Thus, the genomic organization of RUB is similar to that of alphaviruses, the other genus of the Togavirus family, and the 5'-proximal ORF of RUB therefore putatively codes for the nonstructural proteins. Sequences homologous to three regions of nucleotide sequence highly conserved among alphaviruses (a stem-and-loop structure at the 5' end of the genome, a 51-nucleotide conserved sequence near the 5' end of the genome, and a 20-nucleotide conserved sequence at the subgenomic RNA start site) were found in the RUB genomic RNA. Amino acid sequence comparisons between the nonstructural ORF of RUB and alphaviruses revealed only one short (122 amino acids) region of significant homology, indicating that these viruses are only distantly related. This region of homology is located at the NH2 terminus of nsP3 in the alphavirus genome. The RUB nonstructural protein ORF contains two global amino acid motifs conserved in a large number of positive-polarity RNA viruses, a motif indicative of helicase activity and a motif indicative of replicase activity. The order of the helicase motif and the nsP3 homology region in the RUB genome is reversed with respect to the alphavirus genome indicating that a genetic rearrangement has occurred during the evolution of these viruses.

Amino Acid Sequence↗

PubMed2004

Insights in metabolism and toxin production from the complete genome sequence of Clostridium tetani.

The decryption of prokaryotic genome sequences progresses rapidly and provides the scientific community with an enormous amount of information. Clostridial genome sequencing projects have been finished only recently, starting with the genome of the solvent-producing Clostridium acetobutylicum in 2001. A lot of attention has been devoted to the genomes of pathogenic clostridia. In 2002, the genome sequence of C. perfringens, the causative agent of gas gangrene, has been released. Currently in the finishing stage and prior to publication are the genomes of the foodborne botulism-causing C. botulinum and of C. difficile, the causative agent of a wide spectrum of clinical manifestations such as antibiotic-associated diarrhea. Our team sequenced the genome of neuropathogenic C. tetani, a Gram-positive spore-forming bacterium predominantly found in the soil. In deep wound infections it occasionally causes spastic paralysis in humans and vertebrate animals, known as tetanus disease, by the secretion of potent neurotoxin, designated tetanus toxin. The toxin blocks the release of neurotransmitters from presynaptic membranes of interneurons of the spinal cord and the brainstem, thus preventing muscle relaxation. Fortunately, this disease is successfully controlled through immunization with tetanus toxoid, a formaldehyde-treated tetanus toxin, but nevertheless, an estimated 400,000 cases still occur each year, mainly of neonatal tetanus. The World Health Organization has stated that neonatal tetanus is the second leading cause of death from vaccine preventable diseases among children worldwide. This minireview focuses on an analysis of the genome sequence of C. tetani E88, a vaccine production strain, which is a toxigenic non-sporulating variant of strain Massachusetts. The genome consists of a 2,799,250 bp chromosome encoding 2618 open reading frames. The tetanus toxin is encoded on a 74,082 kb plasmid, containing 61 genes. Additional virulence-related factors as well as an insight into the metabolic strategy of C. tetani with regard to its pathogenic phenotype will be presented. The information from other clostridial genomes by means of comparative analysis will also be explored.

Journal Article↗

PubMed2005

The first complete chloroplast genome sequence of a lycophyte, Huperzia lucidula (Lycopodiaceae).

We used a unique combination of techniques to sequence the first complete chloroplast genome of a lycophyte, Huperzia lucidula. This plant belongs to a significant clade hypothesized to represent the sister group to all other vascular plants. We used fluorescence-activated cell sorting (FACS) to isolate the organelles, rolling circle amplification (RCA) to amplify the genome, and shotgun sequencing to 8x depth coverage to obtain the complete chloroplast genome sequence. The genome is 154,373 bp, containing inverted repeats of 15,314 bp each, a large single-copy region of 104,088 bp, and a small single-copy region of 19,657 bp. Gene order is more similar to those of mosses, liverworts, and hornworts than to gene order for other vascular plants. For example, the Huperzia chloroplast genome possesses the bryophyte gene order for a previously characterized 30 kb inversion, thus supporting the hypothesis that lycophytes are sister to all other extant vascular plants. The lycophyte chloroplast genome data also enable a better reconstruction of the basal tracheophyte genome, which is useful for inferring relationships among bryophyte lineages. Several unique characters are observed in Huperzia, such as movement of the gene ndhF from the small single copy region into the inverted repeat. We present several analyses of evolutionary relationships among land plants by using nucleotide data, inferred amino acid sequences, and by comparing gene arrangements from chloroplast genomes. The results, while still tentative pending the large number of chloroplast genomes from other key lineages that are soon to be sequenced, are intriguing in themselves, and contribute to a growing comparative database of genomic and morphological data across the green plants.

Chromosome Mapping↗

PubMed2026

Predicting genome-wide functional constraints with GPN-Star.

Genomic language models have emerged as a powerful approach for learning genome-wide functional constraints directly from DNA sequences1. However, standard genomic language models adapted from natural language processing often require large model sizes and computational resources, yet still fall short of classical evolutionary models in predictive tasks2-4. Here we introduce a genomic pretrained network with species tree and alignment representations (GPN-Star), which is a biologically grounded genomic language model featuring a phylogeny-aware architecture that leverages whole-genome alignments and species trees to model evolutionary relationships explicitly. Trained on alignments spanning vertebrate, mammal and primate evolutionary timescales, GPN-Star achieves state-of-the-art performance across a wide range of variant effect prediction tasks in both coding and non-coding regions of the human genome. Analyses across timescales show task-dependent advantages of modelling more recent versus deeper evolution. To demonstrate its potential to advance human genetics, we show that GPN-Star substantially outperforms previous methods in prioritizing pathogenic and fine-mapped genome-wide association study variants, yields strong enrichments of complex trait heritability and improves power in rare variant association testing5. Extending beyond humans, we train GPN-Star for five model organisms-Mus musculus, Gallus gallus, Drosophila melanogaster, Caenorhabditis elegans and Arabidopsis thaliana-demonstrating the robustness and generalizability of the framework. Taken together, these results position GPN-Star as a scalable, powerful and flexible tool for genome interpretation, well suited to leverage the growing abundance of comparative genomics data.

Journal Article↗

PubMed2004

A bayesian analysis of metazoan mitochondrial genome arrangements.

Genome arrangements are a potentially powerful source of information to infer evolutionary relationships among distantly related taxa. Mitochondrial genome arrangements may be especially informative about metazoan evolutionary relationships because (1) nearly all animals have the same set of definitively homologous mitochondrial genes, (2) mitochondrial genome rearrangement events are rare relative to changes in sequences, and (3) the number of possible mitochondrial genome arrangements is huge, making convergent evolution of genome arrangements appear highly unlikely. In previous studies, phylogenetic evidence in genome arrangement data is nearly always used in a qualitative fashion-the support in favor of clades with similar or identical genome arrangements is considered to be quite strong, but is not quantified. The purpose of this article is to quantify the uncertainty among the relationships of metazoan phyla on the basis of mitochondrial genome arrangements while incorporating prior knowledge of the monophyly of various groups from other sources. The work we present here differs from our previous work in the statistics literature in that (1) we incorporate prior information on classifications of metazoans at the phylum level, (2) we describe several advances in our computational approach, and (3) we analyze a much larger data set (87 taxa) that consists of each unique, complete mitochondrial genome arrangement with a full complement of 37 genes that were present in the NCBI (National Center for Biotechnology Information) database at a recent date. In addition, we analyze a subset of 28 of these 87 taxa for which the non-tRNA mitochondrial genomes are unique where the assumption of our inversion-only model of rearrangement is more plausible. We present summaries of Bayesian posterior distributions of tree topology on the basis of these two data sets.

Animals↗

PubMed1978

Nucleotide sequence relationships between the genomes of an endogenous and an exogenous avian tumor virus.

We have used mapping of large T1 oligonucleotides to examine the genome of Rous-associated virus-O (RAV-O), an endogenous virus of chickens, and to compare it with that of Prague strain Rous sarcoma virus, subgroup B, (Pr-RSV-B), an exogenous sarcoma virus. To extend the sensitivity of such comparisons, we have developed a system of nucleic acid hybridization and hybridization-competition combined with fingerprinting. This method allows us to estimate the relative degree of relatedness of various portions of the viral genomes. From the results of this study, we have concluded that the genomes of Pr-RSV-B and RAV-O are related in the following way. The 5'-terminal half of the genomes (corresponding to the gag and pol regions) is virtually identical, with only scattered single nucleotide differences. This region is followed by a region comprising 25 to 30% of the genome (the env region) which contains substantial nucleotide sequence differences, most or all of which are due to single base changes. The env-coding region can be further subdivided into three regions: a more variable region probably containing sequences coding for subgroup specificity, flanked by relatively common sequences on each side. To the 3' side of the env region, the RAV-O genome contains a very short sequence not found in Pr-RSV-B, whereas the Pr-RSV-B genome contains a much longer unrelated sequence. The central portion of this sequence comprises the src gene as defined by transformation-defective mutants. Particularly striking is the absence, in the RAV-O genome, of any nucleotide sequence related to the "c region" found very near the 3' end of all exogenous tumor viruses. Both the Pr-RSV-B and RAV-O genomes contain the identical terminally redundant sequence of 21 nucleotides near each end of the genome.

Avian Leukosis Virus↗