PubMed Health⌕ Search

SEARCH · PubMed Health

Results for “genome coverage”

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

Insights from human/mouse genome comparisons.

Large-scale public genomic sequencing efforts have provided a wealth of vertebrate sequence data poised to provide insights into mammalian biology. These include deep genomic sequence coverage of human, mouse, rat, zebrafish, and two pufferfish ( Fugu rubripes and Tetraodon nigroviridis) (Aparicio et al. 2002; Lander et al. 2001; Venter et al. 2001; Waterston et al. 2002). In addition, a high-priority has been placed on determining the genomic sequence of chimpanzee, dog, cow, frog, and chicken (Boguski 2002). While only recently available, whole genome sequence data have provided the unique opportunity to globally compare complete genome contents. Furthermore, the shared evolutionary ancestry of vertebrate species has allowed the development of comparative genomic approaches to identify ancient conserved sequences with functionality. Accordingly, this review focuses on the initial comparison of available mammalian genomes and describes various insights derived from such analysis.

Animals↗

Comprehensive genome analysis of 203 genomes provides structural genomics with new insights into protein family space.

We present an analysis of 203 completed genomes in the Gene3D resource (including 17 eukaryotes), which demonstrates that the number of protein families is continually expanding over time and that singleton-sequences appear to be an intrinsic part of the genomes. A significant proportion of the proteomes can be assigned to fewer than 6000 well-characterized domain families with the remaining domain-like regions belonging to a much larger number of small uncharacterized families that are largely species specific. Our comprehensive domain annotation of 203 genomes enables us to provide more accurate estimates of the number of multi-domain proteins found in the three kingdoms of life than previous calculations. We find that 67% of eukaryotic sequences are multi-domain compared with 56% of sequences in prokaryotes. By measuring the domain coverage of genome sequences, we show that the structural genomics initiatives should aim to provide structures for less than a thousand structurally uncharacterized Pfam families to achieve reasonable structural annotation of the genomes. However, in large families, additional structures should be determined as these would reveal more about the evolution of the family and enable a greater understanding of how function evolves.

Algorithms↗

Assessing Hardy-Weinberg equilibrium in T2T-aligned 1000 genomes project.

Quality control of markers in genome-wide association studies often includes testing for Hardy-Weinberg equilibrium (HWE). However, this is usually implemented in a homogeneous population without stratifying by sex. Previous work indicates sex-based selection at numerous autosomal loci in cohorts with active recruitment. Sex chromosome sequences can also interfere with autosomal SNPs. These motivate a re-examination of HWE in sex-aware analyses. Using the telomere-to-telomere (T2Tv2)-aligned high-coverage whole genome sequencing data from 2,490 individuals in the 1000 Genomes Project, we examined genome-wide sex-specific deviations from HWE across five super-populations. Our analyses were restricted to bi-allelic SNPs with non-missing genotypes and minor allele frequency (MAF) &#x2265;5% in both sexes of the five super-populations. We applied an allele-based framework to quantify both the magnitude and direction of Hardy-Weinberg disequilibrium (HWD), followed by a second-order omnibus meta-analysis that combined HWD results across populations and sexes. At a genome-wide significance threshold of p&#x2009;<&#x2009;5e-8, 0.9% of autosomal SNPs exhibited significant deviations from HWE. The majority of these deviations were associated with genomic features indicative of poor sequence quality. Restricting the analysis to reliable genomic regions substantially reduced the number of signals, yielding 255 autosomal SNPs and one non-pseudoautosomal chromosome X SNP. Among these, 140 autosomal SNPs displayed significant heterogeneity across populations but not across sexes. Notably, eight SNPs within a 15-bp region on chromosome 14q31.3 showed excess heterozygosity in both sexes of the African super-population (AFR). Finally, we developed a multivariate predictor of HWD based on sequence features, providing a practical tool that can be integrated into existing quality control pipelines for whole genome sequencing studies.

Journal Article↗

Targeted loss of heterozygosity in Candida albicans using CRISPR-Cas9 reveals the functional impact of allelic variation.

The diploid genome of the fungal pathogen Candida albicans is highly heterozygous, with most allele pairs diverging at either the coding or regulatory level. When faced with selection pressure like antifungal exposure, this hidden genetic diversity can provide a reservoir of adaptive mutations through loss of heterozygosity (LOH) events. Validating the potential phenotypic impact of LOH events observed in clinical or experimentally evolved strains can be difficult due to the challenge of precisely targeting one allele over the other. Here, we show that a CRISPR-Cas9 system can be used to overcome this challenge. By designing allele-specific guide RNA sequences, we can induce targeted, directed LOH events, which we validate by whole-genome long-read sequencing. Using this approach, we efficiently recapitulate a recently described LOH event that increases resistance to the antifungal fluconazole. Additionally, we find that the recombination tracts of these induced LOH events have similar lengths to those observed naturally. To facilitate future use of this method, we provide a database of allele-specific sgRNA sequences for Cas9 that provide near genome-wide coverage of heterozygous sites through either direct or indirect targeting. This approach will be useful in probing the adaptive role of LOH events in this important human pathogen.

Candida albicans↗

ZIPcnv: accurate and efficient inference of copy number variations from shallow whole-genome sequencing.

MOTIVATION: Shallow whole-genome sequencing (sWGS), a rapid and cost-effective sequencing technology, has gradually been widely adopted for CNV analyses. However, with genome&#x2011;wide coverage of only 0.1-5&#xd7;, sWGS data display a pronounced zero&#x2011;inflation phenomenon-a large fraction of loci has zero sequencing reads. Zero inflation causes read counts to fluctuate by several&#x2011;fold between adjacent windows. As a result, random upward blips in coverage can be misinterpreted as copy&#x2011;number gains (false positives), and true deletions often become indistinguishable from pervasive zero&#x2011;coverage noise. In addition, existing CNV detection tools developed for sWGS data often struggle to adapt across different CNV sizes. These combined effects severely constrain the accuracy of CNV inference. RESULTS: To address above challenges, we propose ZIPcnv, a novel CNV detection tool specifically designed for sWGS data. First, we apply a segment sliding window to smooth the raw read depth signal, which transforms the original zero-inflated statistical characteristics into approximately normal distribution characteristics. We then design a statistical process model that robustly detects persistent shifts under high background noise using a cumulative sum strategy, classifying genomic regions into candidate and non-candidate CNV regions. Finally, dynamic sliding windows are used for one-pass detection of CNVs of varying lengths, with window size adapting to the CNV region size. We evaluated the performance of ZIPcnv on simulated data and 190 real whole-genome sequencing samples. Experimental results show that ZIPcnv consistently outperforms currently popular CNV detection tools. AVAILABILITY AND IMPLEMENTATION: The ZIPcnv source code is freely available at https://github.com/Nevermore233/ZIPcnv.

DNA Copy Number Variations↗

Human genome organization.

Recent advances have been made in addressing three intriguing aspects of human genome organization: the organization of protein-coding sequences within chromosomes, the structural basis of the metaphase chromosomal banding pattern, and the function of non protein coding DNA. At the cytogenetic level, R band heterogeneity has been examined by fluorescence in situ hybridization using complex fractions of genomic DNA as probes. DNA fractionated according to GC content and CpG is island density both generated patterns related G and R bands and directly demonstrated regional variations in gene densities. A structural basis for metaphase bands has been proposed that is based on the differential size and packing of DNA loops and matrix attachment sites in G versus R bands. The model presents interesting opportunities for structure/function and organization investigations. At the molecular level, the human genome initiative has resulted in extensive genomic clone coverage for many large chromosomal regions, permitting detailed documentation of CpG islands, base composition and repeat sequence context, as well as fueling comprehensive gene searches. Sequence and functional characteristics are being examined at the kilobase level and are prompting new suggestions of roles for 'junk' DNA. Because of these developments, opportunities are now emerging for direct assessment of the molecular characteristics of individual metaphase bands.

Chromosome Banding↗

Evaluation of sequencing reads at scale using rdeval.

MOTIVATION: Large sequencing datasets are being produced and deposited into public archives at unprecedented rates. The availability of tools that can reliably and efficiently generate and store sequencing read summary statistics has become critical. RESULTS: As part of the effort by the Vertebrate Genomes Project (VGP) to generate high-quality reference genomes at scale, we sought to address the community's need for efficient sequence data evaluation by developing rdeval, a standalone tool to quickly compute and interactively display sequencing read metrics. Rdeval can either run on the fly or store key sequence data metrics in tiny read 'snapshot' files. Statistics can then be efficiently recalled from snapshots for additional processing. Rdeval can convert fa*[.gz] files to and from other popular formats including BAM and CRAM for better compression. Overall, while CRAM achieves the best compression, the gain compared to BAM is marginal, and BAM achieves the best compromise between data compression and access speed. Rdeval also generates a detailed visual report with multiple data analytics that can be exported in various formats. We showcase rdeval's functionalities using long-read data from different sequencing platforms and species, including human. For PacBio long-read sequencing, our analysis shows dramatic improvements in both read length and quality over time, as well as the benefit of increased coverage for genome assembly, though the magnitude varies by taxa. AVAILABILITY AND IMPLEMENTATION: Rdeval is implemented in C++ for data processing and in R for data visualization. Precompiled releases (Linux, MacOS, Windows) and commented source code for rdeval are available under MIT license at https://github.com/vgl-hub/rdeval. Documentation is available on ReadTheDocs (https://rdeval-documentation.readthedocs.io). Rdeval is also available in Bioconda and in Galaxy (https://usegalaxy.org). An automated test workflow ensures the consistency of software updates.

Software↗

Identification of genetic determinants that promote biofilm growth under heterotrophic conditions in Cupriavidus necator using transposon enrichment.

Cupriavidus necator is a metabolically versatile &#x3b2;-proteobacterium of growing interest for auto- and heterotrophic bioprocesses, yet the genetic determinants governing its biofilm formation remain largely uncharacterized, particularly under process-relevant heterotrophic conditions. Here, we applied a forward-genetics transposon-enrichment approach to identify loci which promote surface-associated growth. A high-density mini-Tn5 mutant library (26,185 insertion clones, exceeding the >17,000 required for genome-wide coverage) was cultivated as a biofilm in a microfluidic flow-cell system on fructose for 168&#x202f;h, and the surface-associated community was characterized by deep sequencing. Twelve genes showed significantly elevated insertion frequencies, several with documented links to biofilm formation in other bacteria, including the ferrous-iron uptake system (feoA/feoB), galU, and a GSDEF/EAL dual-domain protein. The gene B2043 (E6A55_RS29530), encoding this c-di-GMP-metabolizing protein, was selected for validation by markerless deletion. Under static conditions, the &#x394;B2043 mutant showed a 1.69&#x202f;&#xb1;&#x202f;0.06-fold increase in biofilm-associated biomass (p&#x202f;=&#x202f;5.16&#x202f;&#xd7;&#x202f;10-15). Under flow-through conditions, the mutant attached faster, entered exponential growth &#x223c;10&#x202f;h earlier, reached its biovolume plateau &#x223c;16&#x202f;h earlier than the wild-type, and formed distinct tower-like structures. These results identify B2043 as a negative regulator of biofilm formation acting predominantly during attachment, provide the first experimental evidence for c-di-GMP-dependent biofilm regulation in C. necator H16, and establish a functional-genomics framework - together with eleven further candidate loci - for engineering productive biofilms in this organism.

Biofilm formation↗

Genomic exploration of the hemiascomycetous yeasts: 12. Kluyveromyces marxianus var. marxianus.

As part of the comparative genomics project 'GENOLEVURES', we studied the Kluyveromyces marxianus var. marxianus strain CBS712 using a partial random sequencing strategy. With a 0.2 x genome equivalent coverage, we identified ca. 1300 novel genes encoding proteins, some containing spliceosomal introns with consensus splice sites identical to those of Saccharomyces cerevisiae, 28 tRNA genes, the whole rDNA repeat, and retrotransposons of the Ty1/2 family of S. cerevisiae with diverged Long Terminal Repeats. Functional classification of the K. marxianus genes, as well as the analysis of the paralogous gene families revealed few differences with respect to S. cerevisiae. Only 42 K. marxianus identified genes are without detectable homolog in the baker's yeast. However, we identified several genetic rearrangements between these two yeast species.

Ascomycota↗

Use of a dense single nucleotide polymorphism map for in silico mapping in the mouse.

Rapid expansion of available data, both phenotypic and genotypic, for multiple strains of mice has enabled the development of new methods to interrogate the mouse genome for functional genetic perturbations. In silico mapping provides an expedient way to associate the natural diversity of phenotypic traits with ancestrally inherited polymorphisms for the purpose of dissecting genetic traits. In mouse, the current single nucleotide polymorphism (SNP) data have lacked the density across the genome and coverage of enough strains to properly achieve this goal. To remedy this, 470,407 allele calls were produced for 10,990 evenly spaced SNP loci across 48 inbred mouse strains. Use of the SNP set with statistical models that considered unique patterns within blocks of three SNPs as an inferred haplotype could successfully map known single gene traits and a cloned quantitative trait gene. Application of this method to high-density lipoprotein and gallstone phenotypes reproduced previously characterized quantitative trait loci (QTL). The inferred haplotype data also facilitates the refinement of QTL regions such that candidate genes can be more easily identified and characterized as shown for adenylate cyclase 7.

Adenylyl Cyclases↗

MetaGLIMPSE: Meta-imputation of low-coverage sequencing data for modern and ancient genomes.

The advent of efficient and accurate imputation for low-coverage sequencing offers an unbiased alternative to SNP array imputation, increasing the accuracy of rare variant imputation across all populations. Since imputation accuracy generally increases with larger reference panels and closer ancestry match between target and reference samples, leveraging imputation from multiple reference panels improves imputation accuracy; however, individual reference panel genotypes are often privacy protected. Meta-imputation bypasses individual-level data by combining single-panel imputed genotypes through estimating panel- and marker-specific weights. We present a meta-imputation method, MetaGLIMPSE, that combines estimates from multiple reference panels for low-coverage sequencing imputation. Across all our scenarios, for both modern and ancient DNA samples, MetaGLIMPSE consistently outperforms the best single-panel imputation for coverages of 0.1&#xd7;-8&#xd7; and across all minor-allele frequencies, equaling the combined panel imputation for some parameters. Finally, MetaGLIMPSE is computationally efficient, meta-imputing 500 whole genomes in 16% of the time of GLIMPSE2.

Humans↗

Construction and characterization of a human chromosome 2-specific BAC library.

We have constructed a human chromosome 2-specific bacterial artificial chromosome (BAC) library using DNA from the somatic cell hybrid GM10826. The average size of the clones is about 63 kb. The coverage and distribution of the library were estimated by screening with known polymorphic genetic markers and fluorescence in situ hybridization (FISH). Twenty-one markers tested positive when DNA pools prepared from approximately one-sixth of the library were screened with 33 known markers. This is consistent with the theoretical calculation of 63% coverage at one genomic equivalent. This suggested that the coverage of the library is approximately 5-6x. FISH analysis with 54 BACs revealed single site hybridization to chromosome 2, and the clones were distributed randomly on the chromosome. We have also performed direct sequencing of the BAC insert ends to generate sequence-tagged sites suitable for mapping and chromosome walking. This is the first reported human chromosome 2-specific BAC library and should provide a resource for physical mapping and disease searching for this chromosome.

Base Sequence↗

Construction of a horse BAC library and cytogenetical assignment of 20 type I and type II markers.

A horse BAC library was constructed with about 40,000 clones and mean insert size of 110 kb representing a 1.5 genome equivalent coverage and a probability of finding a single sequence of 0.75. It was characterized by PCR screening of about 130 sequences of horse microsatellites and exonic gene sequences retrieved from databases. BACs containing 8 microsatellites and 12 genes were subsequently localized by fluorescent in situ hybridization (FISH) on chromosomes. Two linkage groups were newly assigned to chromosomes: LG2 to ECA3 and LG5 to ECA24, and five linkage groups were also oriented--LG3, LG4, LG5, LG8, and LG12--leaving only three groups unassigned. This work showed how this library makes an integrated map a realistic objective for the near future and how it can make comparative mapping more efficient in a search for candidate genes of interest.

Animals↗

Genomic organization of the human PAX3 gene: DNA sequence analysis of the region disrupted in alveolar rhabdomyosarcoma.

Mutations in the human PAX3 gene have previously been associated with two distinct diseases, Waardenburg syndrome and alveolar rhabdomyosarcoma. In this report we establish that the normal human PAX3 gene is encoded by 8 exons. Intron-exon boundary sequences were obtained for PAX3 exons 5, 6, 7, and 8 and together with previous work provide the complete genomic sequence organization for PAX3. Difficulties in obtaining overlapping genomic clone coverage of PAX3 were circumvented in part by RARE cleavage mapping, which showed that the entire PAX3 gene spans 100 kb of chromosome 2. Sequence analysis of the last intron of PAX3, which contains the previously mapped t(2;13)(q35;q14) translocation breakpoints of alveolar rhabdomyosarcoma, revealed the presence of a pair of inverted Alu repeats and a pair of inverted (GT)n-rich microsatellite repeats within a 5-kb region. This work establishes the complete structure of PAX3 and will permit high-resolution analyses of this locus for mutations associated with Waardenburg syndrome, alveolar rhabdomyosarcoma, and other phenotypes for which PAX3 may be a candidate locus.

Alleles↗

Resolving tylenchid evolutionary relationships through multiple gene analysis derived from EST data.

Sequence-based phylogenetic analyses typically are based on a small number of character sets and report gene trees which may not reflect the true species tree. We employed an EST mining strategy to suppress such incongruencies, and recovered the most robust phylogeny for five species of plant-parasitic nematode (Meloidogyne arenaria, M. chitwoodi, M. hapla, M. incognita, and M. javanica), three closely related tylenchid taxa (Heterodera glycines, Globodera pallida, and G. rostochiensis) and a distant taxon, Caenorhabditis elegans. Our multiple-gene approach is based on sampling more than 80,000 publicly available tylenchid EST sequences to identify phylum-wide orthologues. Bayesian inference, minimum evolution, maximum likelihood and protein distance methods were employed for phylogenetic reconstruction and hypothesis tests were constructed to elucidate differential selective pressures across the phylogeny for each gene. Our results place M. incognita and M. javanica as sister taxa, with M. arenaria as the next closely related nematode. Significant differences in selective pressure were revealed for some genes under some hypotheses, though all but one gene are exclusively under purifying selection, indicating conservation across the orthologous groups. This EST-based multi-gene analysis is a first step towards accomplishing genome-wide coverage for tylenchid evolutionary analyses.

Animals↗

Complete coverage of the Schizosaccharomyces pombe genome in yeast artificial chromosomes.

The genome of the fission yeast, Schizosaccharomyces pombe, consists of some 14 million base pairs of DNA contained in three chromosomes. On account of its excellent genetics we used it as a test system for a strategy designed to map mammalian chromosomes and genomes. Data obtained from hybridization fingerprinting established an ordered library of 1,248 yeast artificial chromosome clones with an average size of 535 kilobases. The clones fall into three contigs completely representing the three chromosomes of the organism. This work provides a high resolution physical and clone map of the genome, which has been related to available genetic and physical map information.

Base Composition↗

Identification of anthrax toxin genes in a Bacillus cereus associated with an illness resembling inhalation anthrax.

Bacillus anthracis is the etiologic agent of anthrax, an acute fatal disease among mammals. It was thought to differ from Bacillus cereus, an opportunistic pathogen and cause of food poisoning, by the presence of plasmids pXO1 and pXO2, which encode the lethal toxin complex and the poly-gamma-d-glutamic acid capsule, respectively. This work describes a non-B. anthracis isolate that possesses the anthrax toxin genes and is capable of causing a severe inhalation anthrax-like illness. Although initial phenotypic and 16S rRNA analysis identified this isolate as B. cereus, the rapid generation and analysis of a high-coverage draft genome sequence revealed the presence of a circular plasmid, named pBCXO1, with 99.6% similarity with the B. anthracis toxin-encoding plasmid, pXO1. Although homologues of the pXO2 encoded capsule genes were not found, a polysaccharide capsule cluster is encoded on a second, previously unidentified plasmid, pBC218. A/J mice challenged with B. cereus G9241 confirmed the virulence of this strain. These findings represent an example of how genomics could rapidly assist public health experts responding not only to clearly identified select agents but also to novel agents with similar pathogenic potentials. In this study, we combined a public health approach with genome analysis to provide insight into the correlation of phenotypic characteristics and their genetic basis.

Animals↗

Network methods for diagonal integration of unpaired single-cell multiomics data: a review.

MOTIVATION: Advances in single-cell sequencing have enabled multiomics profiling at unprecedented resolution; however, mass spectrometry-based single-cell proteomics (scMS) remains inherently destructive, precluding simultaneous transcriptomic capture. Unlike antibody-based methods such as CITE-seq, which permit paired profiling but are restricted to targeted protein panels, scMS provides unbiased, genome-scale coverage of the intracellular proteome yet necessitates post hoc integration of unpaired datasets. This diagonal integration challenge, where transcriptomes and proteomes are measured in separate cells lacking shared anchors, remains underserved by existing reviews, which focus predominantly on vertical integration strategies enabled by non-destructive assays. RESULTS: We survey the complete computational pipeline for constructing mechanistic proteogenomic networks from unpaired single-cell data, covering: (i) unimodal network inference such as knowledge-based approaches, probabilistic graphical models, temporal directionality inference, and generative and foundation model strategies that establish the transcriptomic scaffold; (ii) cross-modal integration architectures such as network propagation, graph neural networks (scMRDR, scmFormer, scCotag), and consensus frameworks designed explicitly for the unpaired proteomics setting; and (iii) benchmarking paradigms spanning network reconstruction (BEELINE, GRETA, CausalBench) and multi-task integration evaluation (scMultiBench, SCMMIB), with guidance on metric selection under network sparsity and class imbalance. We identify three principal axes of future development: generative proteomic translation from transcriptomic precursors, inductive prior embedding in next-generation architectures, and perturbation-based causal benchmarking. AVAILABILITY AND IMPLEMENTATION: This is a review article; no novel software is distributed. A curated benchmark resource table, methods starter guide, and per-method bottleneck annotations are provided in the Supplementary Material.

Multiomics↗