Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “protein function annotation”

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,639 records · Page 91Linked to original sources

Alternative splicing and differential gene expression in colon cancer detected by a whole genome exon array.

BACKGROUND: Alternative splicing is a mechanism for increasing protein diversity by excluding or including exons during post-transcriptional processing. Alternatively spliced proteins are particularly relevant in oncology since they may contribute to the etiology of cancer, provide selective drug targets, or serve as a marker set for cancer diagnosis. While conventional identification of splice variants generally targets individual genes, we present here a new exon-centric array (GeneChip Human Exon 1.0 ST) that allows genome-wide identification of differential splice variation, and concurrently provides a flexible and inclusive analysis of gene expression. RESULTS: We analyzed 20 paired tumor-normal colon cancer samples using a microarray designed to detect over one million putative exons that can be virtually assembled into potential gene-level transcripts according to various levels of prior supporting evidence. Analysis of high confidence (empirically supported) transcripts identified 160 differentially expressed genes, with 42 genes occupying a network impacting cell proliferation and another twenty nine genes with unknown functions. A more speculative analysis, including transcripts based solely on computational prediction, produced another 160 differentially expressed genes, three-fourths of which have no previous annotation. We also present a comparison of gene signal estimations from the Exon 1.0 ST and the U133 Plus 2.0 arrays. Novel splicing events were predicted by experimental algorithms that compare the relative contribution of each exon to the cognate transcript intensity in each tissue. The resulting candidate splice variants were validated with RT-PCR. We found nine genes that were differentially spliced between colon tumors and normal colon tissues, several of which have not been previously implicated in cancer. Top scoring candidates from our analysis were also found to substantially overlap with EST-based bioinformatic predictions of alternative splicing in cancer. CONCLUSION: Differential expression of high confidence transcripts correlated extremely well with known cancer genes and pathways, suggesting that the more speculative transcripts, largely based solely on computational prediction and mostly with no previous annotation, might be novel targets in colon cancer. Five of the identified splicing events affect mediators of cytoskeletal organization (ACTN1, VCL, CALD1, CTTN, TPM1), two affect extracellular matrix proteins (FN1, COL6A3) and another participates in integrin signaling (SLC3A2). Altogether they form a pattern of colon-cancer specific alterations that may particularly impact cell motility.

Algorithms↗

The putative malate/lactate dehydrogenase from Pseudomonas putida is an NADPH-dependent delta1-piperideine-2-carboxylate/delta1-pyrroline-2-carboxylate reductase involved in the catabolism of D-lysine and D-proline.

A Pseudomonas putida ATCC12633 gene, dpkA, encoding a putative protein annotated as malate/L-lactate dehydrogenase in various sequence data bases was disrupted by homologous recombination. The resultant dpkA(-) mutant was deprived of the ability to use D-lysine and also D-proline as a sole carbon source. The dpkA gene was cloned and overexpressed in Escherichia coli, and the gene product was characterized. The enzyme showed neither malate dehydrogenase nor lactate dehydrogenase activity but catalyzed the NADPH-dependent reduction of such cyclic imines as Delta(1)-piperideine-2-carboxylate and Delta(1)-pyrroline-2-carboxylate to form L-pipecolate and L-proline, respectively. NADH also served as a hydrogen donor for both substrates, although the reaction rates were less than 1% of those with NADPH. The reverse reactions were also catalyzed by the enzyme but at much lower rates. Thus, the enzyme has dual metabolic functions, and we named the enzyme Delta(1)-piperideine-2-carboxylate/Delta(1)-pyrroline-2-carboxylate reductase, the first member of a novel subclass in a large family of NAD(P)-dependent oxidoreductases.

Bacterial Proteins↗

Evolutionary origins of the endocannabinoid system.

Endocannabinoid system evolution was estimated by searching for functional orthologs in the genomes of twelve phylogenetically diverse organisms: Homo sapiens, Mus musculus, Takifugu rubripes, Ciona intestinalis, Caenorhabditis elegans, Drosophila melanogaster, Saccharomyces cerevisiae, Arabidopsis thaliana, Plasmodium falciparum, Tetrahymena thermophila, Archaeoglobus fulgidus, and Mycobacterium tuberculosis. Sequences similar to human endocannabinoid exon sequences were derived from filtered BLAST searches, and subjected to phylogenetic testing with ClustalX and tree building programs. Monophyletic clades that agreed with broader phylogenetic evidence (i.e., gene trees displaying topographical congruence with species trees) were considered orthologs. The capacity of orthologs to function as endocannabinoid proteins was predicted with pattern profilers (Pfam, Prosite, TMHMM, and pSORT), and by examining queried sequences for amino acid motifs known to serve critical roles in endocannabinoid protein function (obtained from a database of site-directed mutagenesis studies). This novel transfer of functional information onto gene trees enabled us to better predict the functional origins of the endocannabinoid system. Within this limited number of twelve organisms, the endocannabinoid genes exhibited heterogeneous evolutionary trajectories, with functional orthologs limited to mammals (TRPV1 and GPR55), or vertebrates (CB2 and DAGLbeta), or chordates (MAGL and COX2), or animals (DAGLalpha and CB1-like receptors), or opisthokonta (animals and fungi, NAPE-PLD), or eukaryotes (FAAH). Our methods identified fewer orthologs than did automated annotation systems, such as HomoloGene. Phylogenetic profiles, nonorthologous gene displacement, functional convergence, and coevolution are discussed.

Animals↗

A chromosome-scale assembly for the genome of southern corn rootworm, Diabrotica undecimpunctata.

Diabrotica undecimpunctata ssp. howardi, the southern corn rootworm or eastern 12-spotted cucumber beetle, is a generalist insect herbivore that causes damage and yield loss to several crops in North America including maize. Unresolved phylogenetic relationships within and among D. undecimpunctata subspecies are impacting current quarantine policies. We report the chromosome-level haploid genome assembly, icDiaUnde3, constructed using HiFi and Hi-C read data from a single male D. undecimpunctata collected and identified as subspecies howardi based on geographic location and morphology. The primary 1.74 Gbp assembly is scaffolded into 11 chromosome-length scaffolds representing 9 autosomes, a single X chromosome and a supernumerary (B) chromosome (scaffold N50 = 162.8 Mb and L50 = 5). Ab initio and evidence-based structural reference sequence (RefSeq) annotations predicted 18,959 protein-coding genes, in which 99.2% of the 1,367 Benchmark Universal Single-Copy Orthologs from Insecta were complete. Repeat elements occupy 1.26 Gbp (72.33%) of the icDiaUnde3 assembly, with nearly 36% predicted to be retroelements. Alignment of whole chromosomes from icDiaUnde3 with those previously assembled from Diabrotica spp. predicted 2 and 6 autosomal inversions with D. balteata and D. virgifera virgifera, respectively. The mitochondrial genome had an annotated gene order and orientation conserved among beetles. The icDiaUnde3 reference genome assembly is a vital resource for taxonomic, comparative, and functional studies to enhance sustainable crop production.

agriculture↗

Complete reannotation of the Arabidopsis genome: methods, tools, protocols and the final release.

BACKGROUND: Since the initial publication of its complete genome sequence, Arabidopsis thaliana has become more important than ever as a model for plant research. However, the initial genome annotation was submitted by multiple centers using inconsistent methods, making the data difficult to use for many applications. RESULTS: Over the course of three years, TIGR has completed its effort to standardize the structural and functional annotation of the Arabidopsis genome. Using both manual and automated methods, Arabidopsis gene structures were refined and gene products were renamed and assigned to Gene Ontology categories. We present an overview of the methods employed, tools developed, and protocols followed, summarizing the contents of each data release with special emphasis on our final annotation release (version 5). CONCLUSION: Over the entire period, several thousand new genes and pseudogenes were added to the annotation. Approximately one third of the originally annotated gene models were significantly refined yielding improved gene structure annotations, and every protein-coding gene was manually inspected and classified using Gene Ontology terms.

Alternative Splicing↗

Plant metabolomics: towards biological function and mechanism.

Metabolite profiling is a fast growing technology and is useful for phenotyping and diagnostic analyses of plants. It is also rapidly becoming a key tool in functional annotation of genes and in the comprehensive understanding of the cellular response to biological conditions. Metabolomics approaches have recently been used to assess the natural variance in metabolite content between individual plants, an approach with great potential for the improvement of the compositional quality of crops. Here, we assess the contribution of metabolite profiling to these areas.

Botany↗

The BioTools Suite. A comprehensive suite of platform-independent bioinformatics tools.

The BioTools Suite is a set of three comprehensive, platform-independent software packages (PepTool, GeneTool, and ChromaTool) developed for sequence assembly and analysis. In addition to supporting a large number of standard bioinformatics functions, these programs also incorporate a number of useful innovations including uniform graphical-user interface (GUI) design, direct internet connectivity, a novel approach to feature annotation, and a variety of enhanced algorithms for large scale proteome and genome analysis. This article describes the key features, recent changes, and general operation of all three programs.

Computational Biology↗

Rapid automatic detection and alignment of repeats in protein sequences.

Many large proteins have evolved by internal duplication and many internal sequence repeats correspond to functional and structural units. We have developed an automatic algorithm, RADAR, for segmenting a query sequence into repeats. The segmentation procedure has three steps: (i) repeat length is determined by the spacing between suboptimal self-alignment traces; (ii) repeat borders are optimized to yield a maximal integer number of repeats, and (iii) distant repeats are validated by iterative profile alignment. The method identifies short composition biased as well as gapped approximate repeats and complex repeat architectures involving many different types of repeats in the query sequence. No manual intervention and no prior assumptions on the number and length of repeats are required. Comparison to the Pfam-A database indicates good coverage, accurate alignments, and reasonable repeat borders. Screening the Swissprot database revealed 3,000 repeats not annotated in existing domain databases. A number of these repeats had been described in the literature but most were novel. This illustrates how in times when curated databases grapple with ever increasing backlogs, automatic (re)analysis of sequences provides an efficient way to capture this important information.

Algorithms↗

In search of differentially expressed genes and proteins.

A great challenge for modern cell biology is the successful examination of the co-expression of thousands of genes under physiological or pathological conditions and how the expression patterns define the different states of a single cell, tissue or a microorganism. Gene expression can be analyzed today on a large scale by advanced technical approaches for differential screening of proteins and mRNAs. The identification of differentially expressed mRNAs has been successfully applied to understand gene function and the underlying molecular mechanism(-s) of differentiation, development and disease state. Analysis of gene expression by the systematic mapping of thousands of proteins present in a cell or tissue can be achieved by the use of two-dimensional (2D) gel electrophoresis, quantitative computer image analysis, and protein identification techniques. In this article, we comment on some of these techniques and try to stress their advantages and drawbacks. We show how data from RNA/DNA mapping, sequence information from genome projects and protein pattern profiling can be linked with each other and annotated. These comprehensive approaches permit the study of differential gene and protein expressions in cells or tissues.

Animals↗

Gammaherpesvirus lytic gene expression as characterized by DNA array.

Gammaherpesviruses are associated with a number of diseases including lymphomas and other malignancies. Murine gammaherpesvirus 68 (MHV-68) constitutes the most amenable animal model for this family of pathogens. However experimental characterization of gammaherpesvirus gene expression, at either the protein or RNA level, lags behind that of other, better-studied alpha- and beta-herpesviruses. We have developed a cDNA array to globally characterize MHV-68 gene expression profiles, thus providing an experimental supplement to a genome that is chiefly annotated by homology. Viral genes started to be transcribed as early as 3 h postinfection (p.i.), and this was followed by a rapid escalation of gene expression that could be seen at 5 h p.i. Individual genes showed their own transcription profiles, and most genes were still being expressed at 18 h p.i. Open reading frames (ORFs) M3 (chemokine-binding protein), 52, and M9 (capsid protein) were particularly noticeable due to their very high levels of expression. Hierarchical cluster analysis of transcription profiles revealed four main groups of genes and allowed functional predictions to be made by comparing expression profiles of uncharacterized genes to those of genes of known function. Each gene was also categorized according to kinetic class by blocking de novo protein synthesis and viral DNA replication in vitro. One gene, ORF 73, was found to be expressed with alpha-kinetics, 30 genes were found to be expressed with beta-kinetics, and 42 genes were found to be expressed with gamma-kinetics. This fundamental characterization furthers the development of this model and provides an experimental basis for continued investigation of gammaherpesvirus pathology.

Animals↗

Multiomic study of cutaneous T-cell lymphoma reveals single-cell clonal evolution in progression and therapy resistance.

Cutaneous T-cell lymphoma (CTCL) remains a challenging disease due to its significant heterogeneity, therapy resistance, and relentless progression. Multiomics technologies offer the potential to provide uniquely precise views of disease progression and response to therapy. Here, we present a comprehensive multiomics view of CTCL clonal evolution, incorporating exome, whole-genome, epigenome, bulk, single-cell T-cell receptor, and single-cell RNA sequencing of 99 clinically annotated serial skin, peripheral blood, and lymph node samples from 34 patients with CTCL. We leveraged this extensive data set to define the molecular underpinnings of CTCL progression in individual patients at single-cell resolution with the goal of identifying clinically useful biomarkers and therapeutic targets. Our studies identified recurrent progression-associated clonal genomic alterations; we highlight mutation of CCR4, phosphoinositide 3-kinase inhibitor signaling, and programmed cell death protein 1 (PD-1) checkpoint pathways as evasion tactics deployed by malignant T cells. We identified a gain-of-function mutation in STAT3 (D661Y) and demonstrated, using cleavage under targets and release using nuclease (CUT&RUN) and RNA sequencing, that it enhances binding to and transcription of genes in Rho GTPase pathways. With our previous work implicating this pathway in histone deacetylase inhibitor-resistant CTCL, these data provide further support for a previously unrecognized role for Rho GTPase pathway dysregulation in CTCL progression. Recurrent progression-associated mutations were common in the epigenetic modifier EZH2, suggesting that EZH2 inhibition may benefit patients with CTCL. Our findings support an approach in which genomic analysis is widely used for improved disease monitoring, biomarker-informed clinical trial design, and genome-guided therapeutic decision-making. Moreover, these molecular changes present new opportunities for therapeutic targeting in this challenging and incurable cancer.

Multiomics↗

Array2BIO: from microarray expression data to functional annotation of co-regulated genes.

BACKGROUND: There are several isolated tools for partial analysis of microarray expression data. To provide an integrative, easy-to-use and automated toolkit for the analysis of Affymetrix microarray expression data we have developed Array2BIO, an application that couples several analytical methods into a single web based utility. RESULTS: Array2BIO converts raw intensities into probe expression values, automatically maps those to genes, and subsequently identifies groups of co-expressed genes using two complementary approaches: (1) comparative analysis of signal versus control and (2) clustering analysis of gene expression across different conditions. The identified genes are assigned to functional categories based on Gene Ontology classification and KEGG protein interaction pathways. Array2BIO reliably handles low-expressor genes and provides a set of statistical methods for quantifying expression levels, including Benjamini-Hochberg and Bonferroni multiple testing corrections. An automated interface with the ECR Browser provides evolutionary conservation analysis for the identified gene loci while the interconnection with Crème allows prediction of gene regulatory elements that underlie observed expression patterns. CONCLUSION: We have developed Array2BIO - a web based tool for rapid comprehensive analysis of Affymetrix microarray expression data, which also allows users to link expression data to Dcode.org comparative genomics tools and integrates a system for translating co-expression data into mechanisms of gene co-regulation. Array2BIO is publicly available at http://array2bio.dcode.org.

Algorithms↗

I-conotoxin superfamily revisited.

The I-conotoxin superfamily (I-Ctx) is known to have four disulfide bonds with the cysteine arrangement C-C-CC-CC-C-C, and the members inhibit or modify ion channels of nerve cells. Recently, Olivera and co-workers (FEBS J. 2005; 272: 4178-4188) have suggested that the previously described I-Ctx should now be divided into two different gene superfamilies, namely, I1 and I2, in view of their having two different types of signal peptides and exhibiting distinct functions. We have revisited the 28 entries presently grouped as I-Ctx in UniProt Swiss-Prot knowledgebase, and on the basis of in silico analysis have divided them into I1 and I2 superfamilies. The sequence analysis has provided a framework for in silico annotation enabling us to carry out computer-based functional characterization of the UniProtKB/TrEMBL entry Q59AA4 from Conus miles and to predict it as a member of the I2 superfamily. Furthermore, we have predicted the mature toxin of this entry and have proposed that it may be an inhibitor of voltage-gated potassium channels.

Amino Acid Sequence↗

EFICAz: a comprehensive approach for accurate genome-scale enzyme function inference.

EFICAz (Enzyme Function Inference by Combined Approach) is an automatic engine for large-scale enzyme function inference that combines predictions from four different methods developed and optimized to achieve high prediction accuracy: (i) recognition of functionally discriminating residues (FDRs) in enzyme families obtained by a Conservation-controlled HMM Iterative procedure for Enzyme Family classification (CHIEFc), (ii) pairwise sequence comparison using a family specific Sequence Identity Threshold, (iii) recognition of FDRs in Multiple Pfam enzyme families, and (iv) recognition of multiple Prosite patterns of high specificity. For FDR (i.e. conserved positions in an enzyme family that discriminate between true and false members of the family) identification, we have developed an Evolutionary Footprinting method that uses evolutionary information from homofunctional and heterofunctional multiple sequence alignments associated with an enzyme family. The FDRs show a significant correlation with annotated active site residues. In a jackknife test, EFICAz shows high accuracy (92%) and sensitivity (82%) for predicting four EC digits in testing sequences that are <40% identical to any member of the corresponding training set. Applied to Escherichia coli genome, EFICAz assigns more detailed enzymatic function than KEGG, and generates numerous novel predictions.

Amino Acid Sequence↗

Chromosomal, in silico and in vitro expression analysis of cardiovascular-based genes encoding zinc finger proteins.

Three hundred and sixty expressed sequence tags (ESTs) from human heart cDNA libraries corresponding to one hundred and twenty six unique zinc finger proteins (ZFPs) were annotated and classified into seven types of ZFPs as reported previously. Among these 126 cvbZFPs (cardiovascular-based ZFPs), the C(2)H(2)-type and the C(2)C(2)-type are the two major ZFP types which account for more than 80% of ZFP genes present in the cardiovascular system. The expression patterns of 11 randomly selected ZFP genes (at least one for each type) in normal fetal, adult and hypertrophic adult hearts, respectively, were determined using reverse transcriptase-polymerase chain reaction (RT-PCR) analysis. The results suggest that ZFPs may be involved in the processes of either developmental control (downregulated or upregulated expression) or basic cellular functional regulation (constant expression). Interestingly, PAF-1 (peroxisome assembly factor-1), a C(3)HC(4)-type ZFP (RING domain-containing ZFP) showing a downregulated expression pattern in normal tissues was found to be upregulated in hypertrophic adult heart, suggesting a possible role for this fetal gene in the pathogenesis of cardiac hypertrophy. In silico Northern analysis of 15 tissues showed that over 90% of cvbZFPs demonstrate widespread tissue distribution, suggesting the vast majority of ZFPs are functionally shared among tissues. The potential importance of transcriptional repressors in cardiovascular development and disease, such as HFHZ, was supported by the observation that one-third (39 of 126) of cvbZFPs possess this function. Of these, 26 are C(2)H(2)-type and the remaining 13 included 8 C(2)C(2)-type, 1 C(3)HC(4)-type, 1 C(2)HC(4)C(HD)-type, 2 C(3)H-type and 1 combination type. Of particular interest was the observation that ZFPs which contain a KRAB domain are the major subtype present (51. 3% of the total repressors in cvbZFPs). Chromosomal distribution analysis showed that mapping loci of cvbZFP genes are concentrated on chromosomes 1, 3, 6, 8, 10, 11, 12, 19 and X. In particular, chromosome 19 appears to be enriched in ZFP genes with C(2)H(2)-type as the predominant type present. Overall, this report provides a fundamental initial step toward understanding the potential role of ZFPs in regulating cadiac development and disease.

Adult↗

Mining Bacillus subtilis chromosome heterogeneities using hidden Markov models.

We present here the use of a new statistical segmentation method on the Bacillus subtilis chromosome sequence. Maximum likelihood parameter estimation of a hidden Markov model, based on the expectation-maximization algorithm, enables one to segment the DNA sequence according to its local composition. This approach is not based on sliding windows; it enables different compositional classes to be separated without prior knowledge of their content, size and localization. We compared these compositional classes, obtained from the sequence, with the annotated DNA physical map, sequence homologies and repeat regions. The first heterogeneity revealed discriminates between the two coding strands and the non-coding regions. Other main heterogeneities arise; some are related to horizontal gene transfer, some to t-enriched composition of hydrophobic protein coding strands, and others to the codon usage fitness of highly expressed genes. Concerning potential and established gene transfers, we found 9 of the 10 known prophages, plus 14 new regions of atypical composition. Some of them are surrounded by repeats, most of their genes have unknown function or possess homology to genes involved in secondary catabolism, metal and antibiotic resistance. Surprisingly, we notice that all of these detected regions are a + t-richer than the host genome, raising the question of their remote sources.

Bacillus subtilis↗

Evaluation of monocot and eudicot divergence using the sugarcane transcriptome.

Over 40,000 sugarcane (Saccharum officinarum) consensus sequences assembled from 237,954 expressed sequence tags were compared with the protein and DNA sequences from other angiosperms, including the genomes of Arabidopsis and rice (Oryza sativa). Approximately two-thirds of the sugarcane transcriptome have similar sequences in Arabidopsis. These sequences may represent a core set of proteins or protein domains that are conserved among monocots and eudicots and probably encode for essential angiosperm functions. The remaining sequences represent putative monocot-specific genetic material, one-half of which were found only in sugarcane. These monocot-specific cDNAs represent either novelties or, in many cases, fast-evolving sequences that diverged substantially from their eudicot homologs. The wide comparative genome analysis presented here provides information on the evolutionary changes that underlie the divergence of monocots and eudicots. Our comparative analysis also led to the identification of several not yet annotated putative genes and possible gene loss events in Arabidopsis.

Arabidopsis↗

Protein encoding genes in an ancient plant: analysis of codon usage, retained genes and splice sites in a moss, Physcomitrella patens.

BACKGROUND: The moss Physcomitrella patens is an emerging plant model system due to its high rate of homologous recombination, haploidy, simple body plan, physiological properties as well as phylogenetic position. Available EST data was clustered and assembled, and provided the basis for a genome-wide analysis of protein encoding genes. RESULTS: We have clustered and assembled Physcomitrella patens EST and CDS data in order to represent the transcriptome of this non-seed plant. Clustering of the publicly available data and subsequent prediction resulted in a total of 19,081 non-redundant ORF. Of these putative transcripts, approximately 30% have a homolog in both rice and Arabidopsis transcriptome. More than 130 transcripts are not present in seed plants but can be found in other kingdoms. These potential "retained genes" might have been lost during seed plant evolution. Functional annotation of these genes reveals unequal distribution among taxonomic groups and intriguing putative functions such as cytotoxicity and nucleic acid repair. Whereas introns in the moss are larger on average than in the seed plant Arabidopsis thaliana, position and amount of introns are approximately the same. Contrary to Arabidopsis, where CDS contain on average 44% G/C, in Physcomitrella the average G/C content is 50%. Interestingly, moss orthologs of Arabidopsis genes show a significant drift of codon fraction usage, towards the seed plant. While averaged codon bias is the same in Physcomitrella and Arabidopsis, the distribution pattern is different, with 15% of moss genes being unbiased. Species-specific, sensitive and selective splice site prediction for Physcomitrella has been developed using a dataset of 368 donor and acceptor sites, utilizing a support vector machine. The prediction accuracy is better than those achieved with tools trained on Arabidopsis data. CONCLUSION: Analysis of the moss transcriptome displays differences in gene structure, codon and splice site usage in comparison with the seed plant Arabidopsis. Putative retained genes exhibit possible functions that might explain the peculiar physiological properties of mosses. Both the transcriptome representation (including a BLAST and retrieval service) and splice site prediction have been made available on http://www.cosmoss.org, setting the basis for assembly and annotation of the Physcomitrella genome, of which draft shotgun sequences will become available in 2005.

Alternative Splicing↗