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 127 records · Page 7Linked to original sources

Application of emerging technologies in the antiviral field.

Viral diseases pose a serious threat to global public health, agriculture, and biosecurity. Conventional antiviral strategies are often limited by an incomplete understanding of disease mechanisms, poor targeting precision, and slow response times. Emerging technologies are now reshaping the landscape of antiviral research. This review examines the roles of four key frontiers, including organoid models, gene editing, AI-driven molecular design, and synthetic biology. Organoids provide physiologically relevant platforms that model virus-host interactions and disease progression. Viral infections remain a major challenge to human and animal health, agriculture, and biosecurity. Progress in antiviral research is constrained by the complexity of viral pathogenesis, the diversity and rapid evolution of viruses, and the limited translational relevance of some traditional model systems. Recent advances in organoid technology, gene editing, artificial intelligence, and synthetic biology are expanding the toolkit available for antiviral research and development. In this review, we discuss how these four technological frontiers contribute to disease modeling, target discovery, molecular design, and translational innovation. Organoids, in particular, provide physiologically relevant systems for investigating viral infection, tissue tropism, host responses, and pathogenesis. Gene editing tools, such as CRISPR, enable precise manipulation of host and viral genomes, facilitating the development of resistant organisms and next-generation vaccine platforms. AI technologies, including AlphaFold for structure prediction and platforms for de novo protein design, address long-standing bottlenecks in structural biology and offer powerful means to engineer antiviral proteins, antibodies, and vaccine antigens. Synthetic biology, guided by the Design-Build-Test-Learn cycle, integrates computational design, genetic assembly, and functional validation into a cohesive pipeline. Together, these technologies form a synergistic workflow that spans disease modeling, target discovery, molecular design, construction, testing, and iterative optimization. This integrated approach is shifting antiviral development from traditional empirical methods toward more precise, intelligent strategies. The review also highlights ongoing challenges in integration and scalability, stressing that high-quality biological datasets and stronger interdisciplinary collaboration are essential for realizing translational potential. By presenting a cohesive view of these converging methodologies, this review offers a framework to guide the intelligent evolution of antiviral strategies in both human and animal health.

Antiviral↗

Efficiently detecting polymorphisms during the fragment assembly process.

MOTIVATION: Current genomic sequence assemblers assume that the input data is derived from a single, homogeneous source. However, recent whole-genome shotgun sequencing projects have violated this assumption, resulting in input fragments covering the same region of the genome whose sequences differ due to polymorphic variation in the population. While single-nucleotide polymorphisms (SNPs) do not pose a significant problem to state-of-the-art assembly methods, these methods do not handle insertion/deletion (indel) polymorphisms of more than a few bases. RESULTS: This paper describes an efficient method for detecting sequence discrepencies due to polymorphism that avoids resorting to global use of more costly, less stringent affine sequence alignments. Instead, the algorithm uses graph-based methods to determine the small set of fragments involved in each polymorphism and performs more sophisticated alignments only among fragments in that set. Results from the incorporation of this method into the Celera Assembler are reported for the D. melanogaster, H. sapiens, and M. musculus genomes.

Algorithms↗

Genomic exploration of the hemiascomycetous yeasts: 3. Methods and strategies used for sequence analysis and annotation.

The primary analysis of the sequences for our Hemiascomycete random sequence tag (RST) project was performed using a combination of classical methods for sequence comparison and contig assembly, and of specifically written scripts and computer visualization routines. Comparisons were performed first against DNA and protein sequences from Saccharomyces cerevisiae, then against protein sequences from other completely sequenced organisms and, finally, against protein sequences from all other organisms. Blast alignments were individually inspected to help recognize genes within our random genomic sequences despite the fact that only parts of them were available. For each yeast species, validated alignments were used to infer the proper genetic code, to determine codon usage preferences and to calculate their degree of sequence divergence with S. cerevisiae. The quality of each genomic library was monitored from contig analysis of the DNA sequences. Annotated sequences were submitted to the EMBL database, and the general annotation tables produced served as a basis for our comparative description of the evolution, redundancy and function of the Hemiascomycete genomes described in other articles of this issue.

Amino Acid Sequence↗

Comprehensive EST analysis of tomato and comparative genomics of fruit ripening.

A large tomato expressed sequence tag (EST) dataset (152 635 total) was analyzed to gain insights into differential gene expression among diverse plant tissues representing a range of developmental programs and biological responses. These ESTs were clustered and assembled to a total of 31 012 unique gene sequences. To better understand tomato gene expression at a plant system level and to identify differentially expressed and tissue-specific genes, we developed and implemented a digital expression analysis protocol. By clustering genes according to their relative abundance in the various EST libraries, expression patterns of genes across various tissues were generated and genes with similar patterns were grouped. In addition, tissues themselves were clustered for relatedness based on relative gene expression as a means of validating the integrity of the EST data as representative of relative gene expression. Arabidopsis and grape EST collections were also characterized to facilitate cross-species comparisons where possible. Tomato fruit digital expression data was specifically compared with publicly available grape EST data to gain insight into molecular manifestation of ripening processes across diverse taxa and resulted in identification of common transcription factors not previously associated with ripening.

Arabidopsis↗

A genome-wide survey of structural variation between human and chimpanzee.

Structural changes (deletions, insertions, and inversions) between human and chimpanzee genomes have likely had a significant impact on lineage-specific evolution because of their potential for dramatic and irreversible mutation. The low-quality nature of the current chimpanzee genome assembly precludes the reliable identification of many of these differences. To circumvent this, we applied a method to optimally map chimpanzee fosmid paired-end sequences against the human genome to systematically identify sites of structural variation > or = 12 kb between the two species. Our analysis yielded a total of 651 putative sites of chimpanzee deletion (n = 293), insertions (n = 184), and rearrangements consistent with local inversions between the two genomes (n = 174). We validated a subset (19/23) of insertion and deletions using PCR and Southern blot assays, confirming the accuracy of our method. The events are distributed throughout the genome on all chromosomes but are highly correlated with sites of segmental duplication in human and chimpanzee. These structural variants encompass at least 24 Mb of DNA and overlap with > 245 genes. Seventeen of these genes contain exons missing in the chimpanzee genomic sequence and also show a significant reduction in gene expression in chimpanzee. Compared with the pioneering work of Yunis, Prakash, Dutrillaux, and Lejeune, this analysis expands the number of potential rearrangements between chimpanzees and humans 50-fold. Furthermore, this work prioritizes regions for further finishing in the chimpanzee genome and provides a resource for interrogating functional differences between humans and chimpanzees.

Animals↗

Sequence assembly validation by multiple restriction digest fragment coverage analysis.

DNA sequence analysis depends on the accurate assembly of fragment reads for the determination of a consensus sequence. This report examines the possibility of analyzing multiple, independent restriction digests as a method for testing the fidelity of sequence assembly. A dynamic programming algorithm to determine the maximum likelihood alignment of error prone electrophoretic mobility data to the expected fragment mobilities given the consensus sequence and restriction enzymes is derived and used to assess the likelihood of detecting rearrangements in genomic sequencing projects. The method is shown to reliably detect errors in sequence fragment assembly without the necessity of making reference to an overlying physical map. An html form-based interface is available at http:/(/)www.ibc.wustl.edu/services/validate. html.

Algorithms↗

A computational approach for identifying pathogenicity islands in prokaryotic genomes.

BACKGROUND: Pathogenicity islands (PAIs), distinct genomic segments of pathogens encoding virulence factors, represent a subgroup of genomic islands (GIs) that have been acquired by horizontal gene transfer event. Up to now, computational approaches for identifying PAIs have been focused on the detection of genomic regions which only differ from the rest of the genome in their base composition and codon usage. These approaches often lead to the identification of genomic islands, rather than PAIs. RESULTS: We present a computational method for detecting potential PAIs in complete prokaryotic genomes by combining sequence similarities and abnormalities in genomic composition. We first collected 207 GenBank accessions containing either part or all of the reported PAI loci. In sequenced genomes, strips of PAI-homologs were defined based on the proximity of the homologs of genes in the same PAI accession. An algorithm reminiscent of sequence-assembly procedure was then devised to merge overlapping or adjacent genomic strips into a large genomic region. Among the defined genomic regions, PAI-like regions were identified by the presence of homolog(s) of virulence genes. Also, GIs were postulated by calculating G+C content anomalies and codon usage bias. Of 148 prokaryotic genomes examined, 23 pathogenic and 6 non-pathogenic bacteria contained 77 candidate PAIs that partly or entirely overlap GIs. CONCLUSION: Supporting the validity of our method, included in the list of candidate PAIs were thirty four PAIs previously identified from genome sequencing papers. Furthermore, in some instances, our method was able to detect entire PAIs for those only partial sequences are available. Our method was proven to be an efficient method for demarcating the potential PAIs in our study. Also, the function(s) and origin(s) of a candidate PAI can be inferred by investigating the PAI queries comprising it. Identification and analysis of potential PAIs in prokaryotic genomes will broaden our knowledge on the structure and properties of PAIs and the evolution of bacterial pathogenesis.

Bacteria↗

Discovery of human inversion polymorphisms by comparative analysis of human and chimpanzee DNA sequence assemblies.

With a draft genome-sequence assembly for the chimpanzee available, it is now possible to perform genome-wide analyses to identify, at a submicroscopic level, structural rearrangements that have occurred between chimpanzees and humans. The goal of this study was to investigate chromosomal regions that are inverted between the chimpanzee and human genomes. Using the net alignments for the builds of the human and chimpanzee genome assemblies, we identified a total of 1,576 putative regions of inverted orientation, covering more than 154 mega-bases of DNA. The DNA segments are distributed throughout the genome and range from 23 base pairs to 62 mega-bases in length. For the 66 inversions more than 25 kilobases (kb) in length, 75% were flanked on one or both sides by (often unrelated) segmental duplications. Using PCR and fluorescence in situ hybridization we experimentally validated 23 of 27 (85%) semi-randomly chosen regions; the largest novel inversion confirmed was 4.3 mega-bases at human Chromosome 7p14. Gorilla was used as an out-group to assign ancestral status to the variants. All experimentally validated inversion regions were then assayed against a panel of human samples and three of the 23 (13%) regions were found to be polymorphic in the human genome. These polymorphic inversions include 730 kb (at 7p22), 13 kb (at 7q11), and 1 kb (at 16q24) fragments with a 5%, 30%, and 48% minor allele frequency, respectively. Our results suggest that inversions are an important source of variation in primate genome evolution. The finding of at least three novel inversion polymorphisms in humans indicates this type of structural variation may be a more common feature of our genome than previously realized.

Animals↗

Infectious salmon anemia virus (ISAV) genomic segment 3 encodes the viral nucleoprotein (NP), an RNA-binding protein with two monopartite nuclear localization signals (NLS).

Infectious salmon anemia virus (ISAV) is the type species of the genus Isavirus belonging to the Orthomyxoviridae, and causes serious disease in Atlantic salmon (Salmo salar). This study presents the expression and functional analysis of the ISAV genome segment 3, and provides further evidence that it encodes the viral nucleoprotein (NP). The encoded protein was expressed in a baculovirus system, and Western blot analysis showed that it corresponds to the 66-71 kDa structural protein previously found in purified ISAV preparations. RNA-binding activity was established by the interaction of viral and recombinant NP with single-stranded RNA transcribed in vitro. Immunofluorescence studies of infected cells showed the ISAV NP to be an early protein. It locates to the nucleus of infected cells before it is transported to the cytoplasm prior to virus assembly. A similar localization pattern was observed in cells transfected with the NP gene, confirming that the encoded protein has an intrinsic ability to be imported into the nucleus. Two monopartite nuclear localization signals (NLS) at amino acids (230)RPKR(233) and (473)KPKK(476) were identified by computer analysis, and validated by site-directed mutagenesis. In contrast to other orthomyxovirus-NPs, that have several NLSs that function independent of each other, both NLSs had to be present for the ISAV NP protein to be transported into the nucleus, indicating that these motifs cooperate to target the protein to the nucleus.

Amino Acid Sequence↗

Hybrid genome assembly and phenotypic assays reveal carbohydrate metabolism diversity in Lacticaseibacillus strains.

Investigation of carbohydrate metabolism in lactic acid bacteria is essential for the rational selection of strains for fermentation processes, particularly in emerging applications involving non-conventional substrates or building of synthetic microbial consortia. However, establishing robust genotype-phenotype relationships remains challenging, as gene presence alone often fails to explain observed metabolic traits without considering the genomic context and regulatory architecture. In the present study, we combined hybrid genome assembly (Illumina and Oxford Nanopore) with high-throughput phenotype profiling (Biolog GENIII and PM2A) to investigate carbohydrate utilization in five Lacticaseibacillus strains. Phenotypic assays revealed clear intra- and inter-specific variability in substrate utilization. We therefore investigated whether such differences could be attributed to the organization and regulatory context of carbohydrate-associated loci, rather than to gene presence alone. Functional annotation based on COG and CAZyme databases revealed candidate genomic regions potentially involved in carbohydrate metabolism. Comparative analysis between predicted and experimentally observed substrate usage highlighted specific loci associated with carbohydrate utilization profile. The trehalose (tre) operon was conserved across all strains, while at least two distinct cellobiose-associated loci were detected in each genome. Despite the presence of these loci, L. paracasei strains were unable to metabolize cellobiose, a phenotype likely linked to the presence of a downstream TetR-type transcriptional repressor within the cellobiose (cel) operon. Additionally, a genomic region uniquely found in L. rhamnosus strains was associated with gentiobiose utilization, consistent with phenotypic observations. Overall, these findings highlight the importance of integrating phenotypic validation with complete genome context to support the identification of candidate structural and regulatory determinants of carbohydrate utilization in lactic acid bacteria. KEY POINTS: • Phenotype microarrays reveal metabolic traits of interest in isolated strains. • Regulatory context is key to understanding carbohydrate metabolism differences. • Basis of subspecies-dependent cellobiose metabolism in L. paracasei is provided.

Carbohydrate Metabolism↗

The genome sequence of silkworm, Bombyx mori.

We performed threefold shotgun sequencing of the silkworm (Bombyx mori) genome to obtain a draft sequence and establish a basic resource for comprehensive genome analysis. By using the newly developed RAMEN assembler, the sequence data derived from whole-genome shotgun (WGS) sequencing were assembled into 49,345 scaffolds that span a total length of 514 Mb including gaps and 387 Mb without gaps. Because the genome size of the silkworm is estimated to be 530 Mb, almost 97% of the genome has been organized in scaffolds, of which 75% has been sequenced. By carrying out a BLAST search for 50 characteristic Bombyx genes and 11,202 non-redundant expressed sequence tags (ESTs) in a Bombyx EST database against the WGS sequence data, we evaluated the validity of the sequence for elucidating the majority of silkworm genes. Analysis of the WGS data revealed that the silkworm genome contains many repetitive sequences with an average length of <500 bp. These repetitive sequences appear to have been derived from truncated transposons, which are interspersed at 2.5- to 3-kb intervals throughout the genome. This pattern suggests that silkworm may have an active mechanism that promotes removal of transposons from the genome. We also found evidence for insertions of mitochondrial DNA fragments at 9 sites. A search for Bombyx orthologs to Drosophila genes controlling sex determination in the WGS data revealed 11 Bombyx genes and suggested that the sex-determining systems differ profoundly between the two species.

Animals↗

Annotating large genomes with exact word matches.

We have developed a tool for rapidly determining the number of exact matches of any word within large, internally repetitive genomes or sets of genomes. Thus we can readily annotate any sequence, including the entire human genome, with the counts of its constituent words. We create a Burrows-Wheeler transform of the genome, which together with auxiliary data structures facilitating counting, can reside in about one gigabyte of RAM. Our original interest was motivated by oligonucleotide probe design, and we describe a general protocol for defining unique hybridization probes. But our method also has applications for the analysis of genome structure and assembly. We demonstrate the identification of chromosome-specific repeats, and outline a general procedure for finding undiscovered repeats. We also illustrate the changing contents of the human genome assemblies by comparing the annotations built from different genome freezes.

Algorithms↗

The planarian Schmidtea mediterranea as a model for epigenetic germ cell specification: analysis of ESTs from the hermaphroditic strain.

Freshwater planarians have prodigious regenerative abilities that enable them to form complete organisms from tiny body fragments. This plasticity is also exhibited by the planarian germ cell lineage. Unlike many model organisms in which germ cells are specified by localized determinants, planarian germ cells appear to be specified epigenetically, arising postembryonically from stem cells. The planarian Schmidtea mediterranea is well suited for investigating the mechanisms underlying epigenetic germ cell specification. Two strains of S. mediterranea exist: a hermaphroditic strain that reproduces sexually and an asexual strain that reproduces by means of transverse fission. To date, expressed sequence tags (ESTs) have been generated only from the asexual strain. To develop molecular reagents for studying epigenetic germ cell specification, we have sequenced 27,161 ESTs from two developmental stages of the hermaphroditic strain of S. mediterranea; this collection of ESTs represents approximately 10,000 unique transcripts. blast analysis of the assembled ESTs showed that 66% share similarity to sequences in public databases. We annotated the assembled ESTs using Gene Ontology terms as well as conserved protein domains and organized them in a relational database. To validate experimentally the Gene Ontology annotations, we used whole-mount in situ hybridization to examine the expression patterns of transcripts assigned to the biological process "reproduction." Of the 53 genes in this category, 87% were expressed in the reproductive organs. In addition to its utility for studying germ cell development, this EST collection will be an important resource for annotating the planarian genome and studying this animal's amazing regenerative abilities.

Animals↗

The extended tubulin superfamily.

Although most eukaryotic cells can express multiple isotypes of alphabeta-tubulin, the significance of this diversity has not always been apparent. Recent data indicate that particular alphabeta-tubulin isotypes, both genome encoded and those derived by post-translational modification, can directly influence microtubule structure and function--thus validating ideas originally proposed in the multitubulin hypothesis over 25 years ago. It has also become increasingly evident over the past year that some (but intriguingly not all) eukaryotes encode several other tubulin proteins, and to date five further members of the tubulin superfamily, gamma, delta, epsilon, zeta and eta, have been identified. Although the role of gamma-tubulin in the nucleation of microtubule assembly is now well established, far less is known about the functions of delta-, epsilon-, zeta- and eta-tubulin. Recent work has expanded our knowledge of the functions and localisation of these newer members of the tubulin superfamily, and the emerging data suggesting a restricted evolutionary distribution of these 'new' tubulin proteins, conforms to established knowledge of microtubule cell biology. On the basis of current evidence, we predict that delta-, epsilon-, zeta- and eta-tubulin all have functions associated with the centriole or basal body of eukaryotic cells and organisms.

Amino Acid Sequence↗

Base-calling of automated sequencer traces using phred. II. Error probabilities.

Elimination of the data processing bottleneck in high-throughput sequencing will require both improved accuracy of data processing software and reliable measures of that accuracy. We have developed and implemented in our base-calling program phred the ability to estimate a probability of error for each base-call, as a function of certain parameters computed from the trace data. These error probabilities are shown here to be valid (correspond to actual error rates) and to have high power to discriminate correct base-calls from incorrect ones, for read data collected under several different chemistries and electrophoretic conditions. They play a critical role in our assembly program phrap and our finishing program consed.

Base Sequence↗

META-DIFF: a k-mer-based pipeline that detects differentially abundant sequences in metagenomics whole genome sequencing.

Traditional case-control metagenomic studies are constrained by their dependence on taxonomic and functional databases. Because annotation occurs before differential analysis, they are limited to known elements and keep function and taxonomy separate. Although binning strategies have emerged to reconstruct genomes and mitigate this issue, they still require an assembly step, preventing the use of all available sequencing data. Here, we introduce META-DIFF, a pipeline based on differentially abundant k-mers independently of any prior annotation. From those k-mers, it reconstructs longer sequences and provides biological context, as well as the best set of unitigs to discriminate between conditions. Across both taxonomy-centric and functionally-centric benchmarks, it showed robust performance and displayed great reproducibility. It also behaved more conservatively than did other univariate methodologies, i.e. it maintained a high precision at the expense of recall, particularly in conditions of low fold-change and limited sequencing depth. The efficacy of META-DIFF was further validated through its application to a real-world colorectal cancer dataset, which produced both confirmatory and novel results compared with those of previous publications. The pipeline is able to exploit all reads and identify differentially abundant elements, including unknown DNA, prior to annotation. With the guidelines provided, META-DIFF provides users with great exploratory power to unravel microbiome changes.

Metagenomics↗

MaGIC: a program to generate targeted marker sets for genome-wide association studies.

High-throughput genotyping technologies such as DNA pooling and DNA microarrays mean that whole-genome screens are now practical for complex disease gene discovery using association studies. Because it is currently impractical to use all available markers, a subset is typically selected on the basis of required saturation density. Restricting markers to those within annotated genomic features of interest (e.g., genes or exons) or within feature-rich regions, reduces workload and cost while retaining much information. We have designed a program (MaGIC) that exploits genome assembly data to create lists of markers correlated with other genomic features. Marker lists are generated at a user-defined spacing and can target features with a user-defined density. Maps are in base pairs or linkage disequilibrium units (LDUs) as derived from the International HapMap data, which is useful for association studies and fine-mapping. Markers may be selected on the basis of heterozygosity and source database, and single nucleotide polymorphism (SNP) markers may additionally be selected on the basis of validation status. The import function means the method can be used for any genomic features such as housekeeping genes, long interspersed elements (LINES), or Alu repeats in humans, and is also functional for other species with equivalent data. The program and source code is freely available at http://cogent.iop.kcl.ac.uk/MaGIC.cogx.

Algorithms↗

Bacterial secrets of secretion: EuroConference on the biology of type IV secretion processes.

Type IV secretion systems (TFSS) mediate secretion or direct cell-to-cell transfer of virulence factors (proteins or protein-DNA complexes) from many Gram-negative animal, human and plant pathogens, such as Agrobacterium tumefaciens, Bartonella tribocorum, Bordetella pertussis, Brucella suis, Helicobacter pylori, Legionella pneumophila and Rickettsia prowazekii, into eukaryotic cells. Bacterial conjugation is also classified as a TFSS-like process mediating the spread of broad-host-range plasmids between Gram-negative bacteria such as RP4 and R388, which carry antibiotic resistance genes. Genetic, biochemical, cell biological and structural biology experiments led to significant progress in the understanding of several aspects of TFSS processes. X-ray crystallography revealed that homologues of the A. tumefaciens inner membrane-associated proteins VirB11 and VirD4 from H. pylori and R388, respectively, may form channels for substrate translocation or assembly of the transmembrane TFSS machinery. Biochemical and cell biological experiments revealed interactions between components of the periplasmic core components VirB8, VirB9 and VirB10, which may form the translocation channel. Analysis of A. tumefaciens virulence proteins VirE2 and VirF suggested that the periplasmic translocation route of the pertussis toxin from B. pertussis may be more generally valid than previously anticipated. Secretion and modification of toxins from H. pylori and L. pneumophila profoundly affect host cell metabolism, thus entering the discipline of cellular microbiology. Finally, results from genome sequencing projects revealed the presence of up to three TFSS in a single organism, and the analysis of their interplay and adaptation to different functions will be a future challenge. TFSS-carrying plasmids were discovered in different ecosystems, suggesting that genetic exchange may speed up their evolution and adaptation to different cell-cell interactions.

Animals↗