PubMed Health⌕ Search

SEARCH · PubMed Health

Results for “functional inference”

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

UBE2A, which encodes a ubiquitin-conjugating enzyme, is mutated in a novel X-linked mental retardation syndrome.

We report a mutation of UBE2A/HR6A, which encodes a ubiquitin-conjugating enzyme (E2), a member of the ubiquitin proteasome pathway, as the cause of a novel X-linked mental retardation (XLMR) syndrome that affects three males in a two-generation family. A single-nucleotide substitution, c.382C-->T in UBE2A, led to a premature UAG stop codon (Q128X). As a consequence, the predicted polypeptide lacks the 25 C-terminal amino acid residues. The importance of this terminal sequence for UBE2 function is inferred by its conservation in vertebrates and in Drosophila. UBE2A mutations do not appear to significantly contribute to XLMR, since no UBE2A mutations were identified in 15 families with nonsyndromic and 4 families with syndromic idiopathic XLMR previously mapped to intervals encompassing this gene. This is the first description of a mutation in a ubiquitin-conjugating enzyme gene as the cause of a human disease.

Adult↗

Methods in comparative genomics: genome correspondence, gene identification and regulatory motif discovery.

In Kellis et al. (2003), we reported the genome sequences of S. paradoxus, S. mikatae, and S. bayanus and compared these three yeast species to their close relative, S. cerevisiae. Genomewide comparative analysis allowed the identification of functionally important sequences, both coding and noncoding. In this companion paper we describe the mathematical and algorithmic results underpinning the analysis of these genomes. (1) We present methods for the automatic determination of genome correspondence. The algorithms enabled the automatic identification of orthologs for more than 90% of genes and intergenic regions across the four species despite the large number of duplicated genes in the yeast genome. The remaining ambiguities in the gene correspondence revealed recent gene family expansions in regions of rapid genomic change. (2) We present methods for the identification of protein-coding genes based on their patterns of nucleotide conservation across related species. We observed the pressure to conserve the reading frame of functional proteins and developed a test for gene identification with high sensitivity and specificity. We used this test to revisit the genome of S. cerevisiae, reducing the overall gene count by 500 genes (10% of previously annotated genes) and refining the gene structure of hundreds of genes. (3) We present novel methods for the systematic de novo identification of regulatory motifs. The methods do not rely on previous knowledge of gene function and in that way differ from the current literature on computational motif discovery. Based on genomewide conservation patterns of known motifs, we developed three conservation criteria that we used to discover novel motifs. We used an enumeration approach to select strongly conserved motif cores, which we extended and collapsed into a small number of candidate regulatory motifs. These include most previously known regulatory motifs as well as several noteworthy novel motifs. The majority of discovered motifs are enriched in functionally related genes, allowing us to infer a candidate function for novel motifs. Our results demonstrate the power of comparative genomics to further our understanding of any species. Our methods are validated by the extensive experimental knowledge in yeast and will be invaluable in the study of complex genomes like that of the human.

Algorithms↗

Annotation of bacterial genomes using improved phylogenomic profiles.

MOTIVATION: Phylogenomic profiling is a large-scale comparative genomic method used to infer protein function from evolutionary information first described in a binary form by Pellegrini et al. (1999). Here, we propose improvements of this approach including the use of normalized Blastp bit scores, a normalization of the matrix of profiles to take into account the evolutionary distances between bacteria, the definition of a phylogenomic neighborhood based on continuous pairwise distances between genes and an original annotation procedure including the computation of a p-value for each functional assignment. RESULTS: The method presented here increases the number of Ecocyc enzymes identified as being evolutionarily related by about 25% with respect to the original binary form (absent/present) method. The fraction of 'false' positives is shown to be smaller than 20%. Based on their phylogenomic relationships, genes of unknown function can then be automatically related to annotated genes. Each gene annotation predicted is associated with a p-value, i.e. its probability to be obtained by chance. The validity of this method was extensively tested on a large set of genes of known function using the MultiFun database. We find that 50% of 3122 function attributions that can be made at a p-value level of 10(-11) correspond to the actual gene annotation. The method can be readily applied to any newly sequenced microbial genome. In contrast to earlier work on the same topic, our approach avoids the use of arbitrary cut-off values, and provides a reliability estimate of the functional predictions in form of p-values.

Algorithms↗

Support vector machine learning from heterogeneous data: an empirical analysis using protein sequence and structure.

MOTIVATION: Drawing inferences from large, heterogeneous sets of biological data requires a theoretical framework that is capable of representing, e.g. DNA and protein sequences, protein structures, microarray expression data, various types of interaction networks, etc. Recently, a class of algorithms known as kernel methods has emerged as a powerful framework for combining diverse types of data. The support vector machine (SVM) algorithm is the most popular kernel method, due to its theoretical underpinnings and strong empirical performance on a wide variety of classification tasks. Furthermore, several recently described extensions allow the SVM to assign relative weights to various datasets, depending upon their utilities in performing a given classification task. RESULTS: In this work, we empirically investigate the performance of the SVM on the task of inferring gene functional annotations from a combination of protein sequence and structure data. Our results suggest that the SVM is quite robust to noise in the input datasets. Consequently, in the presence of only two types of data, an SVM trained from an unweighted combination of datasets performs as well or better than a more sophisticated algorithm that assigns weights to individual data types. Indeed, for this simple case, we can demonstrate empirically that no solution is significantly better than the naive, unweighted average of the two datasets. On the other hand, when multiple noisy datasets are included in the experiment, then the naive approach fares worse than the weighted approach. Our results suggest that for many applications, a naive unweighted sum of kernels may be sufficient. AVAILABILITY: http://noble.gs.washington.edu/proj/seqstruct

Algorithms↗

Recombination and selection at Brassica self-incompatibility loci.

In Brassica species, self-incompatibility is controlled genetically by haplotypes involving two known genes, SLG and SRK, and possibly an as yet unknown gene controlling pollen incompatibility types. Alleles at the incompatibility loci are maintained by frequency-dependent selection, and diversity at SLG and SRK appears to be very ancient, with high diversity at silent and replacement sites, particularly in certain "hypervariable" portions of the genes. It is important to test whether recombination occurs in these genes before inferences about function of different parts of the genes can be made from patterns of diversity within their sequences. In addition, it has been suggested that, to maintain the relationship between alleles within a given S-haplotype, recombination is suppressed in the S-locus region. The high diversity makes many population genetic measures of recombination inapplicable. We have analyzed linkage disequilibrium within the SLG gene of two Brassica species, using published coding sequences. The results suggest that intragenic recombination has occurred in the evolutionary history of these alleles. This is supported by patterns of synonymous nucleotide diversity within both the SLG and SRK genes, and between domains of the SRK gene. Finally, clusters of linkage disequilibrium within the SLG gene suggest that hypervariable regions are under balancing selection, and are not merely regions of relaxed selective constraint.

Brassica↗

Identification and characterization of a novel calcineurin-binding protein in scallop testis.

Calcineurin has been inferred to function in meiosis and spermiogenesis in testis. Here, we identified a calcineurin-binding protein in scallop testis by Far-Western blot analysis using purified calcineurin as a probe. The molecular mass of the binding protein estimated on the blot was 75 kDa. The isolated cDNA clone encoded a novel 474-residue protein, named CaNBP75. The region between T6 and A210 of CaNBP75 was responsible for the interaction with calcineurin. CaNBP75 was predominantly expressed in testis and ovary of scallop. Thus, CaNBP75 may modulate the physiological function of calcineurin in the testis and ovary of scallop, such as in spermiogenesis or meiosis.

Animals↗

PRINTS and its automatic supplement, prePRINTS.

The PRINTS database houses a collection of protein fingerprints. These may be used to assign uncharacterised sequences to known families and hence to infer tentative functions. The September 2002 release (version 36.0) includes 1800 fingerprints, encoding approximately 11 000 motifs, covering a range of globular and membrane proteins, modular polypeptides and so on. In addition to its continued steady growth, we report here the development of an automatic supplement, prePRINTS, designed to increase the coverage of the resource and reduce some of the manual burdens inherent in its maintenance. The databases are accessible for interrogation and searching at http://www.bioinf.man.ac.uk/dbbrowser/PRINTS/.

Amino Acid Motifs↗

TreeDomViewer: a tool for the visualization of phylogeny and protein domain structure.

Phylogenetic analysis and examination of protein domains allow accurate genome annotation and are invaluable to study proteins and protein complex evolution. However, two sequences can be homologous without sharing statistically significant amino acid or nucleotide identity, presenting a challenging bioinformatics problem. We present TreeDomViewer, a visualization tool available as a web-based interface that combines phylogenetic tree description, multiple sequence alignment and InterProScan data of sequences and generates a phylogenetic tree projecting the corresponding protein domain information onto the multiple sequence alignment. Thereby it makes use of existing domain prediction tools such as InterProScan. TreeDomViewer adopts an evolutionary perspective on how domain structure of two or more sequences can be aligned and compared, to subsequently infer the function of an unknown homolog. This provides insight into the function assignment of, in terms of amino acid substitution, very divergent but yet closely related family members. Our tool produces an interactive scalar vector graphics image that provides orthological relationship and domain content of proteins of interest at one glance. In addition, PDF, JPEG or PNG formatted output is also provided. These features make TreeDomViewer a valuable addition to the annotation pipeline of unknown genes or gene products. TreeDomViewer is available at http://www.bioinformatics.nl/tools/treedom/.

Computer Graphics↗

Positional artifacts in microarrays: experimental verification and construction of COP, an automated detection tool.

Microarray technology is currently one of the most widely-used technologies in biology. Many studies focus on inferring the function of an unknown gene from its co-expressed genes. Here, we are able to show that there are two types of positional artifacts in microarray data introducing spurious correlations between genes. First, we find that genes that are close on the microarray chips tend to have higher correlations between their expression profiles. We call this the 'chip artifact'. Our calculations suggest that the carry-over during the printing process is one of the major sources of this type of artifact, which is later confirmed by our experiments. Based on our experiments, the measured intensity of a microarray spot contains 0.1% (for fully-hybridized spots) to 93% (for un-hybridized ones) of noise resulting from this artifact. Secondly, we, for the first time, show that genes that are close on the microtiter plates in microarray experiments also tend to have higher correlations. We call this the 'plate artifact'. Both types of artifacts exist with different severity in all cDNA microarray experiments that we analyzed. Therefore, we develop an automated web tool-COP (COrrelations by Positional artifacts) to detect these artifacts in microarray experiments. COP has been integrated with the microarray data normalization tool, ExpressYourself, which is available at http://bioinfo.mbb.yale.edu/ExpressYourself/. Together, the two can eliminate most of the common noises in microarray data.

Animals↗

MMDB: annotating protein sequences with Entrez's 3D-structure database.

Three-dimensional (3D) structure is now known for a large fraction of all protein families. Thus, it has become rather likely that one will find a homolog with known 3D structure when searching a sequence database with an arbitrary query sequence. Depending on the extent of similarity, such neighbor relationships may allow one to infer biological function and to identify functional sites such as binding motifs or catalytic centers. Entrez's 3D-structure database, the Molecular Modeling Database (MMDB), provides easy access to the richness of 3D structure data and its large potential for functional annotation. Entrez's search engine offers several tools to assist biologist users: (i) links between databases, such as between protein sequences and structures, (ii) pre-computed sequence and structure neighbors, (iii) visualization of structure and sequence/structure alignment. Here, we describe an annotation service that combines some of these tools automatically, Entrez's 'Related Structure' links. For all proteins in Entrez, similar sequences with known 3D structure are detected by BLAST and alignments are recorded. The 'Related Structure' service summarizes this information and presents 3D views mapping sequence residues onto all 3D structures available in MMDB (http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=structure).

Databases, Protein↗

Nucleotide bias causes a genomewide bias in the amino acid composition of proteins.

We analyzed the nucleotide contents of several completely sequenced genomes, and we show that nucleotide bias can have a dramatic effect on the amino acid composition of the encoded proteins. By surveying the genes in 21 completely sequenced eubacterial and archaeal genomes, along with the entire Saccharomyces cerevisiae genome and two Plasmodium falciparum chromosomes, we show that biased DNA encodes biased proteins on a genomewide scale. The predicted bias affects virtually all genes within the genome, and it could be clearly seen even when we limited the analysis to sets of homologous gene sequences. Parallel patterns of compositional bias were found within the archaea and the eubacteria. We also found a positive correlation between the degree of amino acid bias and the magnitude of protein sequence divergence. We conclude that mutational bias can have a major effect on the molecular evolution of proteins. These results could have important implications for the interpretation of protein-based molecular phylogenies and for the inference of functional protein adaptation from comparative sequence data.

Amino Acids↗

Misleading local sequence alignments: implications for comparative protein modelling.

Although it is well known that significant sequence similarity between proteins is reflected at the structural level, it is commonly assumed that any misaligned regions, as judged by the correct structure based alignment, are those where the local sequence identity is lower than the global. Recent studies have shown that this is not always the case and there can exist short stretches of high local identity which is not reflected in the structure based alignment. An analysis is presented of 290 pairs of homologous proteins with a view to quantifying the occurrence of these misleading local sequence alignments (MLSAs). It is found that such MLSAs are likely if the global sequence identity is less than 40% and can occur even when it is greater than 60%. The results have implications for automated homology modelling and also for the inference of function made by comparison.

Algorithms↗

Exaggerated response to cold stress in a congenic strain for the quantitative trait locus for blood pressure.

OBJECTIVE: Stroke-prone spontaneously hypertensive rats (SHRSP) are known to have sympathetic hyperactivity to various stimuli. In the search for 'intermediate phenotypes' inferring the function of hypertension genes, the present study assessed responsiveness to cold stress in a congenic strain derived from SHRSP/Izm and Wistar-Kyoto/Izm (WKY/Izm). DESIGN: A congenic strain, WKYpch1.0, was established by 10 generations of backcrossing to transfer the chromosomal fragment between D1Wox29 and D1Arb21 of SHRSP to WKY. This fragment covered the 100:1 confidence interval of the quantitative trait locus (QTL) for blood pressure identified in a previous study. Response to cold stress was studied by exposing rats to 4 degrees C for 4 h. Blood pressure was monitored with telemetry. Urine was collected during the exposure, and urinary concentrations of catecholamines were measured by high-performance liquid chromatography. RESULTS: Under the cold stress, urinary excretion of norepinephrine (NE) and vanillylmandelic acid (VMA), as well as the plasma level of NE, was significantly greater in WKYpch1.0 than in WKY. The increase in blood pressure during the cold stress was also greater in WKYpch1.0 than in WKY. Further, neonatal chemical sympathectomy using guanethidine abolished the exaggerated response in blood pressure and in urinary excretion of NE and VMA in WKYpch1.0. CONCLUSION: These results suggested that the QTL region on rat chromosome 1 harbored genes responsible for the exaggerated response of the sympathetic nervous system to the cold stress. The relationship of this with the pathogenesis of hypertension should be elucidated in future studies.

Animals↗

The effect of postural changes on the inferred pressures within the nucleus pulposus during lumbar discography.

By observing the variation of intradiscal pressure occurring at different body postures, it is possible to infer a functional hydrostatic behavior of a lumbar disc. Results from such observations on normal discs are already available. However, observations on degenerate discs are largely restricted to in vitro studies. The authors are now able to report a series of recordings taken from discographically degenerate lumbar discs in patients presenting with low-back pain. In this study of twenty patients, pressure observations were made on 21 normal and 19 abnormal discs. From the results that the authors have obtained, they can reaffirm that normal discs behave predictably and as previously described. The abnormal discs, however, did not behave as a single group. They showed patterns of pressure changes in different postures often dissimilar from that shown by the normal discs both in the absolute values recorded and the sequential changes that occur during the postural change. Unfortunately, the authors were unable to correlate either the extent or character of disc degeneration with the observed variation in pressure changes.

Back Pain↗

End-product control of expression of branched-chain amino acid biosynthesis genes in Streptomyces coelicolor A3(2): paradoxical relationships between DNA sequence and regulatory phenotype.

The branched-chain protein amino acids isoleucine, valine and leucine can provide precursors for synthesis of complex polyketide secondary metabolites in streptomycetes; therefore the regulation of their own synthesis is of interest. DNA sequences upstream of ilvBNC, ilvD, leuA, leuB, ilvE and leuCD in Streptomyces coelicolor A3(2) have been obtained in this laboratory or as part of the S. coelicolor genome sequencing project. Upstream of ilvB and leuA, typical features of classical attenuator systems can be discerned, in particular hypothetical short ORFs with runs of Ile/Val/Leu and Leu codons, respectively. No such features are apparent upstream of other genes or gene clusters present. All five upstream regions were fused to xylE (encoding catechol dioxygenase, CO) as a reporter gene in the SCP2*-based low-copy-number vector pIJ2839. All wild-type regions showed strong depression of CO activity in the presence of all three branched-chain amino acids whether or not the attenuation features were present. By site-directed mutagenesis, the Ile/Val/Leu and Leu triplets in the putative attenuator peptides for ilvB and leuA were replaced by ones for other amino acids. In the case of ilvB, this had no effect at all; for leuA, the wild-type regulatory phenotype persisted in at least some experiments. It was concluded that (i) an unknown regulatory mechanism must be operating in the ilv/leu system of S. coelicolor A3(2) in place of classical attenuation; and (ii) it is unsafe to infer the functioning of a regulatory mechanism from sequence homologies alone.

5' Untranslated Regions↗

Parallel construction of orthologous sequence-ready clone contig maps in multiple species.

Comparison is a fundamental tool for analyzing DNA sequence. Interspecies sequence comparison is particularly powerful for inferring genome function and is based on the simple premise that conserved sequences are likely to be important. Thus, the comparison of a genomic sequence with its orthologous counterpart from another species is increasingly becoming an integral component of genome analysis. In ideal situations, such comparisons are performed with orthologous sequences from multiple species. To facilitate multispecies comparative sequence analysis, a robust and scalable strategy for simultaneously constructing sequence-ready bacterial artificial chromosome (BAC) contig maps from targeted genomic regions has been developed. Central to this approach is the generation and utilization of "universal" oligonucleotide-based hybridization probes ("overgo" probes), which are designed from sequences that are highly conserved between distantly related species. Large collections of these probes are used en masse to screen BAC libraries from multiple species in parallel, with the isolated clones assembled into physical contig maps. To validate the effectiveness of this strategy, efforts were focused on the construction of BAC-based physical maps from multiple mammalian species (chimpanzee, baboon, cat, dog, cow, and pig). Using available human and mouse genomic sequence and a newly developed computer program to design the requisite probes, sequence-ready maps were constructed in all species for a series of targeted regions totaling approximately 16 Mb in the human genome. The described approach can be used to facilitate the multispecies comparative sequencing of targeted genomic regions and can be adapted for constructing BAC contig maps in other vertebrates.

Animals↗

Structure of sperm-specific protein SSP-19 from Caenorhabditis elegans.

Structural data are reported for SSP-19, a sperm-specific protein (SSP) family member from Caenorhabditis elegans. The SSP family [also known as the major sperm protein-like (MSP-like) family] contains proteins with only 107-109 amino acids, compared with 127 amino acids in the major sperm protein (MSP) family. MSP, the most abundant protein in nematode sperm, forms a dynamic actin-like cytoskeleton that provides the framework for the nematode sperm motility. In vivo, MSP dimers polymerize to form filaments that are constructed from two helical strands, which assemble into larger macromolecular structures. Little is known about the SSP family and a similar function is inferred from sequence and structural homology [Pfam (Protein Families Database of Alignments and HMMs) and SCOP (Structural Classification of Proteins) classification]. Despite the overall structural homology, the monomer-monomer interactions in SSP-19 are strikingly different from the interactions in the two MSP canonic domains described previously.

Amino Acid Sequence↗

EMatch: an efficient method for aligning atomic resolution subunits into intermediate-resolution cryo-EM maps of large macromolecular assemblies.

Structural analysis of biological machines is essential for inferring their function and mechanism. Nevertheless, owing to their large size and instability, deciphering the atomic structure of macromolecular assemblies is still considered as a challenging task that cannot keep up with the rapid advances in the protein-identification process. In contrast, structural data at lower resolution is becoming more and more available owing to recent advances in cryo-electron microscopy (cryo-EM) techniques. Once a cryo-EM map is acquired, one of the basic questions asked is what are the folds of the components in the assembly and what is their configuration. Here, a novel knowledge-based computational method, named EMatch, towards tackling this task for cryo-EM maps at 6-10 A resolution is presented. The method recognizes and locates possible atomic resolution structural homologues of protein domains in the assembly. The strengths of EMatch are demonstrated on a cryo-EM map of native GroEL at 6 A resolution.

Algorithms↗