PubMed Health⌕ Search

SEARCH · PubMed Health

Results for “draft genome sequence”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

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

At least 685 records · Page 38Linked to original sources

From chicken coops to genome maps: generating phenotype from the molecular blueprint.

The tools of molecular and cellular biology can be used to precisely describe traits in terms of a sequence of nucleic acids when their molecular and cellular bases are well understood. The entire genome of elite production birds, however, cannot be written as a series of A's, T's, C's, and G's because the interaction between alleles at the same and different loci is too large and there is likely to be many genotypes that encode the same production trait phenotype. A first draft of the genetic map of the chicken is anticipated within the next few years, but a complete molecular description of the genome of birds with elite production characteristics is not anticipated in the near future. Quantitative genetics will remain the cornerstone of breeding programs for production traits. Novel sequences encoding traits such as enhanced nutritional capability (e.g., expression of phytase) and resistance to specific diseases could be introduced into lines of chickens using the tools of molecular and cellular biology. Cloning could be used by the poultry industry to disperse highly desirable genotypes without the need for grandparent and parent flocks for multiplication.

Animals↗

Whole genome sequence of a superbug-Escherichia coli strain KAB-AI-497 isolated from the vagina of a 20 year old pregnant woman with premature rupture of membrane (PROM) in a resource limited setting, Kabale Regional Referral Hospital, in Uganda.

OBJECTIVES: The objective of the study is to sequence the whole genome of multidrug resistant E. coli strain KAB-AI-497 that causes bacterial vaginosis and implicated in premature rupture of membrane in pregnant woman. DATA DESCRIPTION: The DNA of the E. coli strain KAB-AI-497 was extracted using the MagAttract HMW DNA Kit, and the extracted DNA was sequenced using an MGI DNBSEQ G99ARS platform. FastQC was used to perform quality control analysis and the reads were trimmed by Trimmomatic. De novo genome assembly was performed by SPAdes and it resulted to a draft assembled genome that has 5.1 Mb genome size, 153 contigs, and 50.5% GC content. Quality analysis of the assembled genome revealed it has 98.46% completeness and 0.97% contamination. The closest E. coli strain to this strain KAB-AI-497 in terms of similarity was Escherichia coli SMS-3-5 with an average nucleotide identity of 98.43% and genome coverage of 86.19%, which confirmed the species level identity of the strain. The assembled genome was annotated using the NCBI Prokaryotic Genome Annotation Pipeline which identified 4,726 protein coding genes in the strain genome. Furthermore, the annotation revealed the genome has resistant genes responsible for resistance against many antibiotic classes such as tetracycline, fluoroquinolone, and penicillin.

Escherichia coli↗

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↗

Genome rearrangements in mammalian evolution: lessons from human and mouse genomes.

Although analysis of genome rearrangements was pioneered by Dobzhansky and Sturtevant 65 years ago, we still know very little about the rearrangement events that produced the existing varieties of genomic architectures. The genomic sequences of human and mouse provide evidence for a larger number of rearrangements than previously thought and shed some light on previously unknown features of mammalian evolution. In particular, they reveal that a large number of microrearrangements is required to explain the differences in draft human and mouse sequences. Here we describe a new algorithm for constructing synteny blocks, study arrangements of synteny blocks in human and mouse, derive a most parsimonious human-mouse rearrangement scenario, and provide evidence that intrachromosomal rearrangements are more frequent than interchromosomal rearrangements. Our analysis is based on the human-mouse breakpoint graph, which reveals related breakpoints and allows one to find a most parsimonious scenario. Because these graphs provide important insights into rearrangement scenarios, we introduce a new visualization tool that allows one to view breakpoint graphs superimposed with genomic dot-plots.

Algorithms↗

Map error reduction: using genetic and sequence-based physical maps to order closely linked markers.

The Marshfield comprehensive genetic maps are frequently used for linkage and association studies, however, for some regions of these maps the marker order has low level of likelihood ratio support. In order to investigate the level of statistical support and the accuracy of the genetic maps compared to sequence-based physical maps, two approximately 30 cM autosomal regions were selected. The first region was selected from chromosome 3 and consisted predominately of draft sequence. The second region was selected from chromosome 21 and consisted of finished sequence data. The physical order of these markers was based upon their position on Celera (CEL) and Human Genome Project-Santa Cruz (HGP-sc) sequence-based physical maps. The chromosome 3 and 21 regions contained 100 and 61 markers, respectively, on the Marshfield genetic map. The genetic and physical map order was consistent for 88.9 and 89.2% of the markers in the region on chromosome 3 and 21, respectively. Using a novel scoring criterion to assess inconsistent marker order between genetic and physical maps, it was determined that the physical order was likely the correct order for 3.3 and 7.1% of the markers in the chromosome 3 and 21 regions, respectively. To increase the accuracy of the order of markers selected for fine mapping a method is presented which combines information from genetic and sequence-based physical maps.

Genetic Linkage↗

The complete human olfactory subgenome.

Olfactory receptors likely constitute the largest gene superfamily in the vertebrate genome. Here we present the nearly complete human olfactory subgenome elucidated by mining the genome draft with gene discovery algorithms. Over 900 olfactory receptor genes and pseudogenes (ORs) were identified, two-thirds of which were not annotated previously. The number of extrapolated ORs is in good agreement with previous theoretical predictions. The sequence of at least 63% of the ORs is disrupted by what appears to be a random process of pseudogene formation. ORs constitute 17 gene families, 4 of which contain more than 100 members each. "Fish-like" Class I ORs, previously considered a relic in higher tetrapods, constitute as much as 10% of the human repertoire, all in one large cluster on chromosome 11. Their lower pseudogene fraction suggests a functional significance. ORs are disposed on all human chromosomes except 20 and Y, and nearly 80% are found in clusters of 6-138 genes. A novel comparative cluster analysis was used to trace the evolutionary path that may have led to OR proliferation and diversification throughout the genome. The results of this analysis suggest the following genome expansion history: first, the generation of a "tetrapod-specific" Class II OR cluster on chromosome 11 by local duplication, then a single-step duplication of this cluster to chromosome 1, and finally an avalanche of duplication events out of chromosome 1 to most other chromosomes. The results of the data mining and characterization of ORs can be accessed at the Human Olfactory Receptor Data Exploratorium Web site (http://bioinfo.weizmann.ac.il/HORDE).

Chromosome Mapping↗

Chromosome-Level Assembly and Annotation of the Grey Reef Shark (Carcharhinus amblyrhynchos) Genome.

To date less than 5% of shark species have nuclear reference genomes, despite next-generation sequencing advances. Particularly for threatened shark species, there is a lack of reliable genomes which are crucial in facilitating research and conservation applications. We assembled the first nuclear reference genome of the endangered grey reef shark (Carcharhinus amblyrhynchos) using long-read PacBio HiFi and Omni-C sequencing to reach chromosome-level contiguity (36 pseudochromosomes; 2.9 Gbp) and high completeness (94% complete BUSCOs). BRAKER3 annotated 16,505 protein-coding genes after masking repetitive elements which accounted for 59% of the genome. We identified potential X and Y sex chromosomes on pseudochromosomes 36 and 57, respectively. The quality and completeness of the draft genome of C. amblyrhynchos will enable researchers to investigate genetic variations and adaptations specific to this species as well as across other Carcharhinus spp., opening new venues for comparative genomics and advancing conservation genetic applications.

Animals↗

Identification and phenotypic characterization of the cell-division protein CdpA.

Analysis of the automated computer annotation of the early draft phase genome of Lactobacillus acidophilus NCFM revealed the previously discovered S-layer gene slpA and an additional partial ORF with weak similarities to S-layer proteins. The entire gene was sequenced to reveal a 1799-bp gene coding for 599 amino acids with a calculated molecular mass of 64.8 kDa. No transcription or translation signals could be determined in close proximity to the 5'-region. However, a strong putative terminator with a free energy of -16.84 kcal/mol was identified directly downstream of the gene. A PSI-Blast analysis showed similarities to members of S-layer proteins, cell-wall associated proteinases and hexosyl-transferases. Calculation of an unrooted phylogenetic tree with other examples of S-layer proteins and proteinases placed the deduced protein separately from both groups. A derivative of L. acidophilus NCFM was constructed by targeted integration into the gene. SDS-PAGE analysis of non-covalently linked proteins of the cell wall of the mutant, compared to the wild type, revealed the loss of a cell-surface protein. Phenotypic analyses of the mutant revealed significant changes in cell morphology, altered responses to various environmental stresses, and lowered cell adhesion. Based on the in silico and functional analyses, we ascertained that this protein plays a role in cell-wall processing during the growth and cell-cell separation and designated the gene as cell-division protein, cdpA.

Bacterial Adhesion↗

The implications of electronic health record for personalized medicine.

The emerging concept of an electronic health record (EHR) targeted at a patient centric, cross-institutional and longitudinal information entity (possibly spanning the individuals lifetime) has great promise for personalized medicine. In fact, it is probably the only vehicle through which we may truly realize the personalization of medicine beyond population-based genetic profiles that are expected to become part of medication and treatment indications in the near future. The new EHR standards include mechanisms that integrate clinical data with genomic testing results obtained through applying research-type procedures, such as full DNA sequencing, to an individual patient. Although the most optimal process for the utilization of integrated clinical-genomic data in the EHR framework is still unclear, the new Health Level Seven (HL7) Clinical Genomics Draft Standard for Trial Use suggests using the 'encapsulate & bubble-up' approach, which includes two main phases: the encapsulation of raw genomic data and bubbling-up the most clinically significant portions of that data, while associating it with clinical phenotypes residing in the individual's EHR.

Computational Biology↗

Accelerating inference in genomic and proteomic foundation models via speculative decoding.

MOTIVATION: Genomic and protein foundation models (GFMs and PFMs) have demonstrated strong performance in learning the language of DNA and proteins, but their use in large-scale sequence generation is limited by the latency of autoregressive decoding. Because every token triggers a forward pass of a large Transformer, whose inference is relatively slow, long-sequence generation quickly becomes costly. RESULTS: In this work we adapt speculative decoding to a representative GFM: the DNA model DNAGPT and two representative PFMs: ProGen2 and ProtGPT2. We implement a probabilistic variant of speculative decoding, in which a lightweight draft model proposes short token spans and a larger target model verifies or corrects them in parallel, while preserving the target model's sampling distribution. Across all three models we systematically study the effect of speculation window length, temperature, draft architecture and prompt length, and we benchmark tokens per second over multiple runs per configuration. Speculative decoding yields consistent speedups over standard key-value cached decoding, with maximum observed speedup reaching 100% increase, while average gains across models ranging between 20% and 40% (e.g. 1.2×-1.4×), without changing the underlying target model predictions. Our results show that speculative decoding is a practical and model-agnostic strategy for accelerating genomic and proteomic sequence generation without sacrificing prediction quality. AVAILABILITY AND IMPLEMENTATION: All code and results are freely available at https://github.com/Georgakopoulos-Soares-lab/BioSpecDec.

Genomics↗

Sawfish: improving long-read structural variant discovery and genotyping with local haplotype modeling.

MOTIVATION: Structural variants (SVs) play an important role in evolutionary and functional genomics but are challenging to characterize. High-accuracy, long-read sequencing can substantially improve SV characterization when coupled with effective calling methods. While state-of-the-art long-read SV callers are highly accurate, further improvements are achievable by systematically modeling local haplotypes during SV discovery and genotyping. RESULTS: We describe sawfish, an SV caller for mapped high-quality long reads incorporating systematic SV haplotype modeling to improve accuracy and resolution. Assessment against the draft Genome in a Bottle (GIAB) SV benchmark from the T2T-HG002-Q100 diploid assembly shows that sawfish has the highest accuracy among state-of-the-art long-read SV callers across every tested SV size group. Additionally, sawfish maintains the highest accuracy at every tested depth level from 10- to 32-fold coverage, such that other callers required at least 30-fold coverage to match sawfish accuracy at 15-fold coverage. Sawfish also shows the highest accuracy in the GIAB challenging medically relevant genes benchmark, demonstrating improvements in both comprehensive and medically relevant contexts.When joint-genotyping seven samples from CEPH-1463, sawfish has over 9000 more pedigree-concordant calls than other state-of-the-art SV callers, with the highest proportion of concordant SVs (81%). Sawfish's quality model enables selection for an even higher proportion of concordant SVs (88%), while still calling nearly 5000 more pedigree-concordant SVs than other callers. These results demonstrate that sawfish improves on the state-of-the-art for long-read SV calling accuracy across both individual and joint-sample analyses. AVAILABILITY AND IMPLEMENTATION: Sawfish source code, pre-compiled Linux binaries, and documentation are released on GitHub: https://github.com/PacificBiosciences/sawfish.

Haplotypes↗

An SNP resource for rice genetics and breeding based on subspecies indica and japonica genome alignments.

Dense coverage of the rice genome with polymorphic DNA markers is an invaluable tool for DNA marker-assisted breeding, positional cloning, and a wide range of evolutionary studies. We have aligned drafts of two rice subspecies, indica and japonica, and analyzed levels and patterns of genetic diversity. After filtering multiple copy and low quality sequence, 408,898 candidate DNA polymorphisms (SNPs/INDELs) were discerned between the two subspecies. These filters have the consequence that our data set includes only a subset of the available SNPs (in particular excluding large numbers of SNPs that may occur between repetitive DNA alleles) but increase the likelihood that this subset is useful: Direct sequencing suggests that 79.8% +/- 7.5% of the in silico SNPs are real. The SNP sample in our database is not randomly distributed across the genome. In fact, 566 rice genomic regions had unusually high (328 contigs/48.6 Mb/13.6% of genome) or low (237 contigs/64.7 Mb/18.1% of genome) polymorphism rates. Many SNP-poor regions were substantially longer than most SNP-rich regions, covering up to 4 Mb, and possibly reflecting introgression between the respective gene pools that may have occurred hundreds of years ago. Although 46.2% +/- 8.3% of the SNPs differentiate other pairs of japonica and indica genotypes, SNP rates in rice were not predictive of evolutionary rates for corresponding genes in another grass species, sorghum. The data set is freely available at http://www.plantgenome.uga.edu/snp.

Breeding↗

Haplotype-specific sequence encoding the protein kinase, interferon-inducible double-stranded RNA-dependent activator in the human leukocyte antigen class II region.

The protein kinase, interferon-inducible double-stranded (ds)RNA-dependent activator (PRKRA) is a dsRNA-binding protein which activates a protein kinase participating in the antiviral activity of interferon. Our previous studies indicated that the nucleotide sequence encoding PRKRA, which appeared to be an intronless gene, was present in PAC HS265J14 containing the human leukocyte antigen (HLA) DR subregion. In this study, we further investigated and characterized the PRKRA gene on the human genome by means of Southern blotting and polymerase chain reaction with homozygous typing cell lines for HLA genes. Results indicated that the presence of PRKRA in the DR subregion was dependent on the DR53 group. Consistently, fluorescence in situ hybridization profiles with PRKRA as a probe showed that the hybridization signal on Chromosome (Chr) 6p21.3 was seen only in the samples carrying the DR haplotypes that belonged to the DR53 group. Interestingly, another hybridization signal, which was mapped on Chr 2q31.2-q32.1, was always detected in the samples examined, i.e., even in the samples negative for the DR53 group. The outcome of a sequence-database homology search further indicated that the PRKRA gene with introns appeared to be present in a recently opened draft-sequence, RP11-65L3 (GenBank accession number AC009948), which is located between D2S335 and D2S2257. Together, the data presented here indicate that the PRKRA gene in the DR subregion is a processed pseudogene (PRKRApsi), which could have been generated only on the DR53 common ancestor's genome, and that the master copy of PRKRApsi is most probably present on Chr 2q31.2-q32.1.

Alleles↗

Homology-driven assembly of a sequence-ready mouse BAC contig map spanning regions related to the 46-Mb gene-rich euchromatic segments of human chromosome 19.

Draft sequence derived from the 46-Mb gene-rich euchromatic portion of human chromosome 19 (HSA19) was utilized to generate a sequence-ready physical map spanning homologous regions of mouse chromosomes. Sequence similarity searches with the human sequence identified more than 1000 individual orthologous mouse genes from which 382 overgo probes were developed for hybridization. Using human gene order and spacing as a model, these probes were used to isolate and assemble bacterial artificial chromosome (BAC) clone contigs spanning homologous mouse regions. Each contig was verified, extended, and joined to neighboring contigs by restriction enzyme fingerprinting analysis. Approximately 3000 mouse BACs were analyzed and assembled into 44 contigs with a combined length of 41.4 Mb. These BAC contigs, covering 90% of HSA19-related mouse DNA, are distributed throughout 15 homology segments derived from different regions of mouse chromosomes 7, 8, 9, 10, and 17. The alignment of the HSA19 map with the ordered mouse BAC contigs revealed a number of structural differences in several overtly conserved homologous regions and more precisely defined the borders of the known regions of HSA19-syntenic homology. Our results demonstrate that given a human draft sequence, BAC contig maps can be constructed quickly for comparative sequencing without the need for preestablished mouse-specific genetic or physical markers and indicate that similar strategies can be applied with equal success to genomes of other vertebrate species.

Animals↗

Transgenic gene knock-outs: functional genomics and therapeutic target selection.

The completion of the first draft of the human genome presents both a tremendous opportunity and enormous challenge to the pharmaceutical industry since the whole community, with few exceptions, will soon have access to the same pool of candidate gene sequences from which to select future therapeutic targets. The commercial imperative to select and pursue therapeutically relevant genes from within the overall content of the genome will be particularly intense for those gene families that currently represent the chemically tractable or 'drugable' gene targets. As a consequence the emphasis within exploratory research has shifted towards the evaluation and adoption of technology platforms that can add additional value to the gene selection process, either through functional studies or direct/indirect measures of disease alignment e.g., genetics, differential gene expression, proteomics, tissue distribution, comparative species data etc. The selection of biological targets for the development of potential new medicines relies, in part, on the quality of the in vivo biological data that correlates a particular molecular target with the underlying pathophysiology of a disease. Within the pharmaceutical industry, studies employing transgenic animals and, in particular, animals with specific gene deletions are playing an increasingly important role in the therapeutic target gene selection, drug candidate selection and product development phases of the overall drug discovery process. The potential of phenotypic information from gene knock-outs to contribute to a high-throughput target selection/validation strategy has hitherto been limited by the resources required to rapidly generate and characterise a large number of knock-out transgenics in a timely fashion. The offerings of several companies that provide an opportunity to overcome these hurdles, albeit at a cost, are assessed with respect to the strategic business needs of the pharmaceutical industry.

Animals↗

Comparative genomic analysis of Streptococcus parasuis and Streptococcus suis reveals mobile element-associated enrichment of antimicrobial resistance and lack of detectable same-MGE colocalization with virulence-associated genes within stable species boundaries.

Streptococcus suis is a major porcine pathogen and a zoonotic agent that causes meningitis and septicemia in humans. Streptococcus parasuis, a recently recognized close relative, remains poorly characterized with regard to its clinical significance and genomic features. In this study, we generated a single-contig closed genome assembly with genome-wide DNA methylation profiles for S. parasuis strain A1, isolated from a diseased pig in Xinjiang, China, and complemented in silico genomic predictions with isolate-level experimental validation of antimicrobial resistance (AMR) genotypes, virulence genotypes, and phenotypic susceptibility for this reference strain. Using this high-quality genome as a reference anchor, we performed comparative genomic analyses across 195 streptococcal genomes, comprising 15 S. parasuis and 180 S. suis strains, to distinguish genome-level co-occurrence of resistance and virulence determinants from their physical colocalization on the same mobile genetic element (MGE).Species boundaries remained clearly delineated at the genomic level, with a median interspecies average nucleotide identity (ANI) of approximately 86.0%, compared with intraspecies ANI medians of 97.5% for S. parasuis and 96.2% for S. suis. Pangenome analysis identified 12,693 gene clusters, of which 1086 were core clusters, and functional annotation revealed significant differences in accessory gene repertoires between the two species. Within this stable genomic framework, S. parasuis genomes carried a higher AMR gene burden; strain A1 harbored 10 AMR genes, multiple virulence-associated genes, three genomic islands, and eight prophage regions. For strain A1, PCR validation confirmed six AMR genes and six virulence genes, and disk diffusion testing demonstrated a multidrug-resistant phenotype consistent with the genotypic profile.Among 235 predicted mobile elements, 19 harbored AMR genes and seven carried Virulence Factor Database (VFDB) homologs, but none carried both categories simultaneously. This finding reflects a lack of detectable same-MGE colocalization under the applied annotation and assembly framework; it should not be interpreted as evidence of biological physical decoupling. Under a random-placement model, the expected number of co-carrying regions was only 0.57, and the probability of observing zero co-carrying regions was P = 0.55. This negative result should be interpreted with caution, given the limited number of cargo-bearing regions and the predominantly draft status of most genomes. Furthermore, the A1 genome contained multiple restriction-modification systems, showed depletion of several methylation motif families in mobile regions, and had limited CRISPR spacer matching evidence, suggesting prior exposure to the relevant sequence space. None of the genomes met our predefined criteria for whole-genome convergence.Collectively, our results support a model in which S. parasuis accumulates AMR-related genes in a modular fashion via mobile elements within stable species boundaries, with no detectable same-MGE colocalization of AMR and virulence determinants under our analytical pipeline. These findings imply that AMR surveillance strategies for this species should prioritize tracking mobile genetic elements rather than inferring wholesale genomic convergence toward S. suis.

Streptococcus suis↗

Unraveling genomic regulatory networks in the simple chordate, Ciona intestinalis.

The draft genome of the primitive chordate, Ciona intestinalis, was published three years ago. Since then, significant progress has been made in utilizing Ciona's genomic and morphological simplicity to better understand conserved chordate developmental processes. Extensive annotation and sequencing of staged EST libraries make the Ciona genome one of the best annotated among those that are publicly available. The formation of the Ciona tadpole depends on simple, well-defined cellular lineages, and it is possible to trace the lineages of key chordate tissues such as the notochord and neural tube to the fertilized egg. Electroporation methods permit the targeted expression of regulatory genes and signaling molecules in defined cell lineages, as well as the rapid identification of regulatory DNAs underlying cell-specific gene expression. The recent sequencing of a second Ciona genome (C. savignyi) permits the use of simple alignment algorithms for the identification of conserved noncoding sequences, including microRNA genes and enhancers. Detailed expression profiles are now available for almost every gene that encodes a regulatory protein or cell-signaling molecule. The combination of gene-expression profiles, comparative genome analysis, and gene-disruption assays should permit the determination of high-resolution genomic regulatory networks underlying the specification of basic chordate tissues such as the heart, blood, notochord, and neural tube.

Animals↗