Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Genome alignment”

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 1,027 records · Page 57Linked to original sources

A computational and experimental approach toward a priori identification of alternatively spliced exons.

Alternative splicing is a powerful means of regulating gene expression and enhancing protein diversity. In fact, the majority of metazoan genes encode pre-mRNAs that are alternatively spliced to produce anywhere from two to tens of thousands of mRNA isoforms. Thus, an important part of determining the complete proteome of an organism is developing a catalog of all mRNA isoforms. Alternatively spliced exons are typically identified by aligning EST clusters to reference mRNAs or genomic DNA. However, this approach is not useful for genomes that lack robust EST coverage, and tools that enable accurate prediction of alternatively spliced exons would be extraordinarily useful. Here, we use comparative genomics to identify, and experimentally verify, potential alternative exons based solely on their high degree of conservation between Drosophila melanogaster and D. pseudoobscura. At least 40% of the exons that fit our prediction criteria are in fact alternatively spliced. Thus, comparative genomics can be used to accurately predict certain classes of alternative exons without relying on EST data.

Alternative Splicing↗

The gene CSTF2T, encoding the human variant CstF-64 polyadenylation protein tauCstF-64, lacks introns and may be associated with male sterility.

Messenger RNA polyadenylation in male germ cells does not seem to require the AAUAAA polyadenylation signal required in all other cell types. To account for this difference, we found a variant form of the polyadenylation protein, the 64,000 Mr protein of the cleavage stimulation factor (CstF-64), in mouse meiotic and postmeiotic germ cells. This protein is a candidate to alter polyadenylation in those cells. More recently, we reported the cloning from mouse pachytene spermatocytes of mouse tauCstF-64 (gene symbol Cstf2t), which is a homolog of CstF-64 fitting the criteria we expected for the variant CstF-64 protein. Here we report the cloning and mapping of the human ortholog of mouse tauCstF-64. The human tauCstF-64 cDNA (gene symbol CSTF2T) is 2324 bp in length and encodes a protein of 616 amino acids (64,442.90 Da). Although most highly related to mouse tauCstF-64 (89.8% identity), human tauCstF-64 is also related to the human and mouse somatic CstF-64 (74.9% and 73.4% identity, respectively). Alignment of human tauCstF-64 with human genome sequence from chromosome 10 shows that CSTF2T lacks introns. Radiation hybrid mapping places the human tauCstF-64 gene at 10q22-q23, which is the site of a translocation that has been associated with human neurological problems and male infertility.

Amino Acid Sequence↗

Isolation of monoclonal antibodies-escape variant of dengue virus serotype 1.

INTRODUCTION: During an outbreak from December 2004 to March 2005, 138 isolates of dengue virus were prospectively obtained from acute-phase serum samples of 1,067 patients with the provisional clinical diagnosis of acute dengue illness admitted to the adult wards of Hospital Tengku Ampuan Rahimah, Klang, Malaysia. Of the 138 dengue virus isolates, 87, 11, 24 and 3 were typed as dengue serotypes 1, 2, 3 and 4, respectively, by a commercial dengue virus typing kit using monoclonal antibodies (Mab). 13 dengue virus isolates could not be assigned to any specific serotype by serotyping Mab and molecular typing using dengue-type specific molecular typing primer pairs. We report the associated clinical features and limited molecular genetics of this Mab-escape dengue virus variant. METHODS: Limited molecular characterisation of the Mab-escape dengue virus variants with respect to a few concurrently isolated dengue serotype 1 virus was performed by reverse transcriptase polymerase chain reaction (RT-PCR), followed by nucleic acid sequencing of the 500-bp dengue virus partial genomic capsid-PreM fragment. RESULTS: The aligned nucleic acid sequence of RT-PCR products showed that these Mab-escape variants were of identical nucleic acid sequence, and shared the highest sequence homology (99 percent) with dengue virus serotype 1 (GeneBank accession No. AB178040) isolated from a Japanese patient in 2004. Though these Mab-escape dengue virus variants were noted to replicate to a 2-log higher titre than the current circulating dengue virus serotype 1, there was no significant difference between these variants and the currently circulating dengue virus serotype 1 with respect to disease severity (dengue fever versus dengue haemorrhagic fever) and clinical features. CONCLUSION: There was no significant difference in the proportion of patients developing dengue haemorrhagic fever following acute infection by Mab-escape dengue virus 1 variant in comparison with infection by the conventional dengue virus 1. Similarly, there was no significant difference in the pattern of clinical presentations following acute infection by the two different strains of virus.

Amino Acid Sequence↗

Evolution and classification of cystine knot-containing hormones and related extracellular signaling molecules.

The cystine knot three-dimensional structure is found in many extracellular molecules and is conserved among divergent species. The identification of proteins with a cystine knot structure is difficult by commonly used pairwise alignments because the sequence homology among these proteins is low. Taking advantage of complete genome sequences in diverse organisms, we used a complementary approach of pattern searches and pairwise alignments to screen the predicted protein sequences of five model species (human, fly, worm, slime mold, and yeast) and retrieved proteins with low sequence homology but containing a typical cystine knot signature. Sequence comparison between proteins known to have a cystine knot three-dimensional structure (transforming growth factor-beta, glycoprotein hormone, and platelet-derived growth factor subfamily members) identified new crucial amino acid residues (two hydrophilic amino acid residues flanking cysteine 5 of the cystine knot). In addition to the well known members of the cystine knot superfamily, novel subfamilies of proteins (mucins, norrie disease protein, von Willebrand factor, bone morphogenetic protein antagonists, and slit-like proteins) were identified as putative cystine knot-containing proteins. Phylogenetic analysis revealed the ancient evolution of these proteins and the relationship between hormones [e.g. transforming growth factor-beta (TGFbeta)] and extracellular matrix proteins (e.g. mucins). They are absent in the unicellular yeast genome but present in nematode, fly, and higher species, indicating that the cystine knot structure evolved in extracellular signaling molecules of multicellular organisms. All data retrieved by this study can be viewed at http://hormone.stanford.edu/.

Amino Acid Motifs↗

Payer perspectives on genomic testing in the United States: A systematic literature review.

PURPOSE: Health care stakeholders' perspectives on the value of genomic testing vary widely and directly affect the access and practice of genomic medicine. To our knowledge, a review of US health care payers' perspectives on genomic testing has not been performed. METHODS: We conducted a systematic literature review of US payers' perspectives on genomic testing in the MEDLINE, PubMed, and Cumulative Index to Nursing and Allied Health Literature (CINAHL) databases. Of the 161 nonduplicate records screened, we summarized findings from 20 included records, and using the framework method, common domains were recorded. RESULTS: Domains included clinical utility, coverage decision frameworks, potential harms, costs, paying for research, demand/pressure, the flexibility of outcomes considered, and personal utility. There was consensus on the definition of clinical utility as improved health outcomes, and the nuances of genomic testing were reported as challenging to fit within existing coverage decision frameworks. Perspectives varied on accepting broader outcomes or uses of genomic testing and whether costs influence coverage decisions. Study methodologies were heterogeneous. CONCLUSION: A deeper understanding of how payers approach genomic testing may allow comparison with other stakeholders' perspectives and may identify challenges, opportunities, and solutions to align a conceptual and evidentiary framework better to demonstrate the value of genomic testing.

Humans↗

On the inference of parsimonious indel evolutionary scenarios.

Given a multiple alignment of orthologous DNA sequences and a phylogenetic tree for these sequences, we investigate the problem of reconstructing a most parsimonious scenario of insertions and deletions capable of explaining the gaps observed in the alignment. This problem, called the Indel Parsimony Problem, is a crucial component of the problem of ancestral genome reconstruction, and its solution provides valuable information to many genome functional annotation approaches. We first show that the problem is NP-complete. Second, we provide an algorithm, based on the fractional relaxation of an integer linear programming formulation. The algorithm is fast in practice, and the solutions it produces are, in most cases, provably optimal. We describe a divide-and-conquer approach that makes it possible to solve very large instances on a simple desktop machine, while retaining guaranteed optimality. Our algorithms are tested and shown efficient and accurate on a set of 1.8 Mb mammalian orthologous sequences in the CFTR region.

Algorithms↗

Fundamentals of massive automatic pairwise alignments of protein sequences: theoretical significance of Z-value statistics.

MOTIVATION: Different automatic methods of sequence alignments are routinely used as a starting point for homology searches and function inference. Confidence in an alignment probability is one of the major fundamentals of massive automatic genome-scale pairwise comparisons, for clustering of putative orthologs and paralogs, sequenced genome annotation or multiple-genomic tree constructions. Extreme value distribution based on the Karlin-Altschul model, usually advised for large-scale comparisons are not always valid, particularly in the case of comparisons of non-biased with nucleotide-biased genomes (such that of Plasmodium falciparum). Z-values estimates based on Monte Carlo technics, can be calculated experimentally for any alignment output, whatever the method used. Empirically, a Z-value higher than approximately 8 is supposed reasonable to assess that an alignment score is significant, but this arbitrary figure was never theoretically justified. RESULTS: In this paper, we used the Bienaymé-Chebyshev inequality to demonstrate a theorem of the upper limit of an alignment score probability (or P-value). This theorem implies that a computed Z-value is a statistical test, a single-linkage clustering criterion and that 1/Z-value(2) is an upper limit to the probability of an alignment score whatever the actual probability law is. Therefore, this study provides the missing theoretical link between a Z-value cut-off used for an automatic clustering of putative orthologs and/or paralogs, and the corresponding statistical risk in such genome-scale comparisons (using non-biased or biased genomes).

Algorithms↗

Replication fork density increases during DNA synthesis in X. laevis egg extracts.

Duplication of the eukaryotic genome depends on the temporal and spatial organization of DNA replication during the cell cycle. To investigate the genomic organization of DNA replication in a higher eukaryote, multiple origins of replication must be simultaneously analyzed over large regions of the genome as DNA synthesis progresses through S phase of the cell cycle. We have employed a novel technique that allows for the quantitative analysis of DNA replication on a genome wide basis. The technique involves stretching and aligning individual DNA molecules on a glass surface. As a model system, Xenopus laevis egg extract was used to differentially label sperm chromatin at successive time points after the start of DNA synthesis. The differentially labeled DNA allows earlier and later replicating sequences to be distinguished, and hence the sites of DNA synthesis at any given time can be directly visualized. Genomic DNA was extracted, and measurements made on the linearized molecules provided a comprehensive analysis of the spatial and temporal organization of DNA replication in the X. laevis in vitro replication system. It was found that: (i) DNA synthesis initiates asynchronously at irregular intervals but continuously as DNA replication advances; and (ii) that the frequency of initiation (the number of activated origins per kilobase) increases as DNA synthesis nears completion. The implications of these findings for the regulation of DNA replication in early embryos is discussed.

Animals↗

Evolutionary implications of three novel members of the human sarcomeric myosin heavy chain gene family.

Sarcomeric myosin heavy chain (MyHC) is the major contractile protein of striated muscle. Six tandemly linked skeletal MyHC genes on chromosome 17 and two cardiac MyHC genes on chromosome 14 have been previously described in the human genome. We report the identification of three novel human sarcomeric MyHC genes on chromosomes 3, 7, and 20, which are notable for their atypical size and intron-exon structure. Two of the encoded proteins are structurally most like the slow-beta MyHC, whereas the third one is closest to the adult fast IIb isoform. Data from pairwise comparisons of aligned coding sequences imply the existence of ancestral genomes with four sarcomeric genes before the emergence of a dedicated smooth muscle MyHC gene. To further address the evolutionary relationships of the distinct sarcomeric and nonsarcomeric rod sequences, we have identified and further annotated human genomic DNA sequences corresponding to 14 class-II MyHCs. An extensive analysis provides a timeline for intron gain and loss, gene contraction and expansion, and gene conversion among genes encoding class-II myosins. One of the novel human genes is found to have introns at positions shared only with the molluscan catchin/MyHC gene, providing evidence for the structure of a pre-Cambrian ancestral gene.

Amino Acid Sequence↗

Panaln: indexing pangenome for read alignment.

MOTIVATION: Pangenome indexing is a critical supporting technology in biological sequence analysis such as read alignment applications. The need to accurately identify billions of small sequencing fragments carrying sequencing errors and genomic variants drives the development of scalable and efficient pangenome indexing approach. RESULTS: We propose a new wavelet tree-based approach, called Panaln, for indexing pangenome and introduce a batch computation approach for fast count query over Panaln. We present a simple and effective seeding strategy and develop a pangenome program that uses the seed-and-extend paradigm for read alignment. Experimental results on simulated and real data demonstrate that Panaln uses significantly less space for the compared pangenome methods with generally higher accuracy. We provide a scalable index construction by representing pangenome with a linear model. Additionally, Panaln brings enhanced accuracy compared to the popular single reference methods. AVAILABILITY AND IMPLEMENTATION: Package: https://anaconda.org/bioconda/panaln and source code: https://github.com/Lilu-guo/Panaln.

Software↗

New local potential useful for genome annotation and 3D modeling.

A new potential energy function representing the conformational preferences of sequentially local regions of a protein backbone is presented. This potential is derived from secondary structure probabilities such as those produced by neural network-based prediction methods. The potential is applied to the problem of remote homolog identification, in combination with a distance-dependent inter-residue potential and position-based scoring matrices. This fold recognition jury is implemented in a Java application called JThread. These methods are benchmarked on several test sets, including one released entirely after development and parameterization of JThread. In benchmark tests to identify known folds structurally similar to (but not identical with) the native structure of a sequence, JThread performs significantly better than PSI-BLAST, with 10% more structures identified correctly as the most likely structural match in a fold library, and 20% more structures correctly narrowed down to a set of five possible candidates. JThread also improves the average sequence alignment accuracy significantly, from 53% to 62% of residues aligned correctly. Reliable fold assignments and alignments are identified, making the method useful for genome annotation. JThread is applied to predicted open reading frames (ORFs) from the genomes of Mycoplasma genitalium and Drosophila melanogaster, identifying 20 new structural annotations in the former and 801 in the latter.

Animals↗

Rfam: annotating non-coding RNAs in complete genomes.

Rfam is a comprehensive collection of non-coding RNA (ncRNA) families, represented by multiple sequence alignments and profile stochastic context-free grammars. Rfam aims to facilitate the identification and classification of new members of known sequence families, and distributes annotation of ncRNAs in over 200 complete genome sequences. The data provide the first glimpses of conservation of multiple ncRNA families across a wide taxonomic range. A small number of large families are essential in all three kingdoms of life, with large numbers of smaller families specific to certain taxa. Recent improvements in the database are discussed, together with challenges for the future. Rfam is available on the Web at http://www.sanger.ac.uk/Software/Rfam/ and http://rfam.wustl.edu/.

Animals↗

Mining genomic databases to identify novel hydrogen producers.

The realization that fossil fuel reserves are limited and their adverse effect on the environment has forced us to look into alternative sources of energy. Hydrogen is a strong contender as a future fuel. Biological hydrogen production ranges from 0.37 to 3.3 moles H(2) per mole of glucose and, considering the high theoretical values of production (4.0 moles H(2) per mole of glucose), it is worth exploring approaches to increase hydrogen yields. Screening the untapped microbial population is a promising possibility. Sequence analysis and pathway alignment of hydrogen metabolism in complete and incomplete genomes has led to the identification of potential hydrogen producers.

Bacteria↗

Molecular characterization of the human protein kinase C theta gene locus (PRKCQ).

Members of the protein kinase C (PKC) family of serine/threonine kinases, in particular PKCtheta, play critical roles in the regulation of differentiation and proliferation of T lymphocytes. In this study the genomic structure of the human PRKCQ gene that encodes PKCtheta was determined. Two genomic P1 clones were isolated from human P1 libraries using the PKCtheta cDNA as a probe and have been used to confirm the assignment of the single PRKCQ locus to chromosome 10p15 by FISH analysis. The PRKCQ locus, the first mammalian PKC gene locus characterized so far, spans approximately 62 kb and is composed of 15 coding exons and 14 introns, varying in size between 98 and 16000 bp. All exon-intron boundaries have been determined by long-range PCR and subsequent DNA sequence analysis. Comparison with other known genomic PKC genes reveals a high degree of homology to the genomic organization of the Drosophila melanogaster dPRKC gene. Alignment of the intron positions in the PRKCQ gene with the intron locations in the dPRKC gene indicates that the sites of seven of the 14 PRKCQ introns are exactly conserved. Exons 5 (32 bp), 11 (174 bp) and 12 (92 bp) share highest similarity in size, organization and primary structure with their counterparts in the Drosophila gene. On the basis of this knowledge of the genomic PRKCQ locus, a directed search for potential genetic polymorphisms and/or genetic abnormalities involved in human genetic disease(s) can now be initiated.

Animals↗

Image analysis for comparative genomic hybridization based on a karyotyping program for windows.

OBJECTIVE: To describe the image processing and analysis techniques we developed for the quantitative analysis of comparative genomic hybridization (CGH). STUDY DESIGN: A system for CGH cytometry is based on a semiautomated karyotyping program using the Windows graphic user interface. RESULTS: After alignment and normalization of the fluorescence images of the test genome (FITC) and reference genome (TRITC), the chromosomes are segmented and arranged as a CGH karyogram. The CGH karyograms from different metaphases of one tumor sample are represented as a CGH sum karyogram. Mean DAPI, FITC, TRITC and RATIO images, as well as ratio profiles with the 95% confidence interval, can be displayed. CONCLUSION: Representation of CGH results in the form of pseudocolored mean ratio chromosomes enhances the visibility of the method. The ratio profile plus confidence interval facilitates the identification of preparation artifacts and the classification of chromosomal imbalances. The sum karyograms of several tumors can be combined into a CGH superkaryogram of a tumor subgroup; that helps identify recurrent DNA changes in a tumor and might lead to genetic tumor classification based on CGH.

DNA, Neoplasm↗

Computational analysis of full-length mouse cDNAs compared with human genome sequences.

Although the sequencing of the human genome is complete, identification of encoded genes and determination of their structures remain a major challenge. In this report, we introduce a method that effectively uses full-length mouse cDNAs to complement efforts in carrying out these difficult tasks. A total of 61,227 RIKEN mouse cDNAs (21,076 full-length and 40,151 EST sequences containing certain redundancies) were aligned with the draft human sequences. We found 35,141 non-redundant genomic regions that showed a significant alignment with the mouse cDNAs. We analyzed the structures and compositional properties of the regions detected by the full-length cDNAs, including cross-species comparisons, and noted a systematic bias of GENSCAN against exons of small size and/or low GC-content. Of the cDNAs locating the 35,141 genomic regions, 3,217 did not match any sequences of the known human genes or ESTs. Among those 3,217 cDNAs, 1,141 did not show any significant similarity to any protein sequence in the GenBank non-redundant protein database and thus are candidates for novel genes.

Algorithms↗

A comparative genetic map of the turkey genome.

Genetic markers (microsatellites and SNPs) were used to create and compare maps of the turkey and chicken genomes. A physical map of the chicken genome was built by comparing sequences of turkey markers with the chicken whole-genome sequence by BLAST analysis. A genetic linkage map of the turkey genome (Meleagris gallopavo) was developed by segregation analysis of genetic markers within the University of Minnesota/Nicholas Turkey Breeding Farms (UMN/NTBF) resource population. This linkage map of the turkey genome includes 314 loci arranged into 29 linkage groups. An additional 40 markers are tentatively placed within linkage groups based on two-point LOD scores and 16 markers remain unlinked. Total map distance contained within linkage groups is 2,011 cM with the longest linkage group (47 loci) measuring 413.3 cM. Average marker interval over the 29 linkage groups was 6.4 cM. All but one turkey linkage group could be aligned with the physical map of the chicken genome. The present genetic map of the turkey provides a comparative framework for future genomic studies.

Animals↗

A probabilistic model for the evolution of RNA structure.

BACKGROUND: For the purposes of finding and aligning noncoding RNA gene- and cis-regulatory elements in multiple-genome datasets, it is useful to be able to derive multi-sequence stochastic grammars (and hence multiple alignment algorithms) systematically, starting from hypotheses about the various kinds of random mutation event and their rates. RESULTS: Here, we consider a highly simplified evolutionary model for RNA, called "The TKF91 Structure Tree" (following Thorne, Kishino and Felsenstein's 1991 model of sequence evolution with indels), which we have implemented for pairwise alignment as proof of principle for such an approach. The model, its strengths and its weaknesses are discussed with reference to four examples of functional ncRNA sequences: a riboswitch (guanine), a zipcode (nanos), a splicing factor (U4) and a ribozyme (RNase P). As shown by our visualisations of posterior probability matrices, the selected examples illustrate three different signatures of natural selection that are highly characteristic of ncRNA: (i) co-ordinated basepair substitutions, (ii) co-ordinated basepair indels and (iii) whole-stem indels. CONCLUSIONS: Although all three types of mutation "event" are built into our model, events of type (i) and (ii) are found to be better modeled than events of type (iii). Nevertheless, we hypothesise from the model's performance on pairwise alignments that it would form an adequate basis for a prototype multiple alignment and genefinding tool.

Evolution, Molecular↗