As the volume of data relating to proteins increases, researchers rely more and more on the analysis of published data, thus increasing the importance of good access to these data that vary from the supplemental material of individual articles, all the way to major reference databases with professional staff and long-term funding. Specialist protein resources fill an important middle ground, providing interactive web interfaces to their databases for a focused topic or family of proteins, using specialized approaches that are not feasible in the major reference databases. Many are labors of love, run by a single lab with little or no dedicated funding and there are many challenges to building and maintaining them. This perspective arose from a meeting of several specialist protein resources and major reference databases held at the Wellcome Trust Genome Campus (Cambridge, UK) on August 11 and 12, 2014. During this meeting some common key challenges involved in creating and maintaining such resources were discussed, along with various approaches to address them. In laying out these challenges, we aim to inform users about how these issues impact our resources and illustrate ways in which our working together could enhance their accuracy, currency, and overall value.

Cardiolipin is a glycerophospholipid found predominantly in the mitochondrial membranes of eukaryotes and in bacterial membranes. Cardiolipin interacts with protein complexes and plays pivotal roles in cellular energy metabolism, membrane dynamics, and stress responses. We recently identified the mitochondrial phosphatase, PTPMT1, as the enzyme that converts phosphatidylglycerolphosphate (PGP) to phosphatidylglycerol, a critical step in the de novo biosynthesis of cardiolipin. Upon examination of PTPMT1 evolutionary distribution, we found a PTPMT1-like phosphatase in the bacterium Rhodopirellula baltica. The purified recombinant enzyme dephosphorylated PGP in vitro. Moreover, its expression restored cardiolipin deficiency and reversed growth impairment in a Saccharomyces cerevisiae mutant lacking the yeast PGP phosphatase, suggesting that it is a bona fide PTPMT1 ortholog. When ectopically expressed, this bacterial PGP phosphatase was localized in the mitochondria of yeast and mammalian cells. Together, our results demonstrate the conservation of function between bacterial and mammalian PTPMT1 orthologs.

The major human intestinal pathogen Giardia lamblia is a very early branching eukaryote with a minimal genome of broad evolutionary and biological interest.

Green fluorescent protein (GFP) has been found in a wide range of Cnidaria, a basal group of metazoans in which it is associated with pigmentation, fluorescence, and light absorbance. A GFP has been recently discovered in the pigmentless chordate Branchiostoma floridae (amphioxus) that shows intense fluorescence mainly in the head region.

Myriapods (e.g., centipedes and millipedes) display a simple homonomous body plan relative to other arthropods. All members of the class are terrestrial, but they attained terrestriality independently of insects. Myriapoda is the only arthropod class not represented by a sequenced genome. We present an analysis of the genome of the centipede Strigamia maritima. It retains a compact genome that has undergone less gene loss and shuffling than previously sequenced arthropods, and many orthologues of genes conserved from the bilaterian ancestor that have been lost in insects. Our analysis locates many genes in conserved macro-synteny contexts, and many small-scale examples of gene clustering. We describe several examples where S. maritima shows different solutions from insects to similar problems. The insect olfactory receptor gene family is absent from S. maritima, and olfaction in air is likely effected by expansion of other receptor gene families. For some genes S. maritima has evolved paralogues to generate coding sequence diversity, where insects use alternate splicing. This is most striking for the Dscam gene, which in Drosophila generates more than 100,000 alternate splice forms, but in S. maritima is encoded by over 100 paralogues. We see an intriguing linkage between the absence of any known photosensory proteins in a blind organism and the additional absence of canonical circadian clock genes. The phylogenetic position of myriapods allows us to identify where in arthropod phylogeny several particular molecular mechanisms and traits emerged. For example, we conclude that juvenile hormone signalling evolved with the emergence of the exoskeleton in the arthropods and that RR-1 containing cuticle proteins evolved in the lineage leading to Mandibulata. We also identify when various gene expansions and losses occurred. The genome of S. maritima offers us a unique glimpse into the ancestral arthropod genome, while also displaying many adaptations to its specific life history.

Insulin homeostasis in pancreatic β cells is now recognized as a critical element in the progression of obesity and type II diabetes (T2D). Proteins that interact with insulin to direct its sequential synthesis, folding, trafficking, and packaging into reserve granules in order to manage release in response to elevated glucose remain largely unknown. Using a conformation-based approach combined with mass spectrometry, we have generated the insulin biosynthetic interaction network (insulin BIN), a proteomic roadmap in the β cell that describes the sequential interacting partners of insulin along the secretory axis. The insulin BIN revealed an abundant C2 domain-containing transmembrane protein 24 (TMEM24) that manages glucose-stimulated insulin secretion from a reserve pool of granules, a critical event impaired in patients with T2D. The identification of TMEM24 in the context of a comprehensive set of sequential insulin-binding partners provides a molecular description of the insulin secretory pathway in β cells.

Specialized metabolic enzymes biosynthesize chemicals of ecological importance, often sharing a pedigree with primary metabolic enzymes. However, the lineage of the enzyme chalcone isomerase (CHI) remained unknown. In vascular plants, CHI-catalysed conversion of chalcones to chiral (S)-flavanones is a committed step in the production of plant flavonoids, compounds that contribute to attraction, defence and development. CHI operates near the diffusion limit with stereospecific control. Although associated primarily with plants, the CHI fold occurs in several other eukaryotic lineages and in some bacteria. Here we report crystal structures, ligand-binding properties and in vivo functional characterization of a non-catalytic CHI-fold family from plants. Arabidopsis thaliana contains five actively transcribed genes encoding CHI-fold proteins, three of which additionally encode amino-terminal chloroplast-transit sequences. These three CHI-fold proteins localize to plastids, the site of de novo fatty-acid biosynthesis in plant cells. Furthermore, their expression profiles correlate with those of core fatty-acid biosynthetic enzymes, with maximal expression occurring in seeds and coinciding with increased fatty-acid storage in the developing embryo. In vitro, these proteins are fatty-acid-binding proteins (FAPs). FAP knockout A. thaliana plants show elevated α-linolenic acid levels and marked reproductive defects, including aberrant seed formation. Notably, the FAP discovery defines the adaptive evolution of a stereospecific and catalytically 'perfected' enzyme from a non-enzymatic ancestor over a defined period of plant evolution.

Many cancer cells show distorted epigenetic landscapes. The Cancer Genome Atlas (TCGA) project profiles thousands of tumors, allowing the discovery of somatic alterations in the epigenetic machinery and the identification of potential cancer drivers among members of epigenetic protein families.

ATR and Chk1 are two key protein kinases in the replication checkpoint. Activation of ATR-Chk1 has been extensively investigated, but checkpoint termination and replication fork restart are less well understood. Here, we report that DNA damage not only activates Chk1, but also exposes a degron-like region at the carboxyl terminus of Chk1 to an Fbx6-containing SCF (Skp1-Cul1-F box) E3 ligase, which mediates the ubiquitination and degradation of Chk1 and, in turn, terminates the checkpoint. The protein levels of Chk1 and Fbx6 showed an inverse correlation in both cultured cancer cells and in human breast tumor tissues. Further, we show that low levels of Fbx6 and consequent impairment of replication stress-induced Chk1 degradation are associated with cancer cell resistance to the chemotherapeutic agent, camptothecin. We propose that Fbx6-dependent Chk1 degradation contributes to S phase checkpoint termination and that a defect in this mechanism might increase tumor cell resistance to certain anticancer drugs.

About 10% of all protein kinases are predicted to be enzymatically inactive pseudokinases, but the structural details of kinase inactivation have remained unclear. We present the first structure of a pseudokinase, VRK3, and that of its closest active relative, VRK2. Profound changes to the active site region underlie the loss of catalytic activity, and VRK3 cannot bind ATP because of residue substitutions in the binding pocket. However, VRK3 still shares striking structural similarity with VRK2, and appears to be locked in a pseudoactive conformation. VRK3 also conserves residue interactions that are surprising in the absence of enzymatic function; these appear to play important architectural roles required for the residual functions of VRK3. Remarkably, VRK3 has an "inverted" pattern of sequence conservation: although the active site is poorly conserved, portions of the molecular surface show very high conservation, suggesting that they form key interactions that explain the evolutionary retention of VRK3.

Eukarya, Archaea, and some Bacteria encode all or part of the essential mevalonate (MVA) metabolic pathway clinically modulated using statins. Curiously, two components of the MVA pathway are often absent from archaeal genomes. The search for these missing elements led to the discovery of isopentenyl phosphate kinase (IPK), one of two activities necessary to furnish the universal five-carbon isoprenoid building block, isopentenyl diphosphate (IPP). Unexpectedly, we now report functional IPKs also exist in Bacteria and Eukarya. Furthermore, amongst a subset of species within the bacterial phylum Chloroflexi, we identified a new enzyme catalyzing the missing decarboxylative step of the putative alternative MVA pathway. These results demonstrate, for the first time, a functioning alternative MVA pathway. Key to this pathway is the catalytic actions of a newly uncovered enzyme, mevalonate phosphate decarboxylase (MPD) and IPK. Together, these two discoveries suggest that unforeseen variation in isoprenoid metabolism may be widespread in nature. DOI: http://dx.doi.org/10.7554/eLife.00672.001.

To reconstruct the evolutionary origin of multicellular animals from their unicellular ancestors, the genome sequences of diverse unicellular relatives are essential. However, only the genome of the choanoflagellate Monosiga brevicollis has been reported to date. Here we completely sequence the genome of the filasterean Capsaspora owczarzaki, the closest known unicellular relative of metazoans besides choanoflagellates. Analyses of this genome alter our understanding of the molecular complexity of metazoans' unicellular ancestors showing that they had a richer repertoire of proteins involved in cell adhesion and transcriptional regulation than previously inferred only with the choanoflagellate genome. Some of these proteins were secondarily lost in choanoflagellates. In contrast, most intercellular signalling systems controlling development evolved later concomitant with the emergence of the first metazoans. We propose that the acquisition of these metazoan-specific developmental systems and the co-option of pre-existing genes drove the evolutionary transition from unicellular protists to metazoans.

Gastric cancer is the second leading cause of worldwide cancer mortality, yet the underlying genomic alterations remain poorly understood. Here we perform exome and transcriptome sequencing and SNP array assays to characterize 51 primary gastric tumours and 32 cell lines. Meta-analysis of exome data and previously published data sets reveals 24 significantly mutated genes in microsatellite stable (MSS) tumours and 16 in microsatellite instable (MSI) tumours. Over half the patients in our collection could potentially benefit from targeted therapies. We identify 55 splice site mutations accompanied by aberrant splicing products, in addition to mutation-independent differential isoform usage in tumours. ZAK kinase isoform TV1 is preferentially upregulated in gastric tumours and cell lines relative to normal samples. This pattern is also observed in colorectal, bladder and breast cancers. Overexpression of this particular isoform activates multiple cancer-related transcription factor reporters, while depletion of ZAK in gastric cell lines inhibits proliferation. These results reveal the spectrum of genomic and transcriptomic alterations in gastric cancer, and identify isoform-specific oncogenic properties of ZAK.

Recent advances in sequencing technologies have greatly increased the identification of mutations in cancer genomes. However, it remains a significant challenge to identify cancer-driving mutations, since most observed missense changes are neutral passenger mutations. Various computational methods have been developed to predict the effects of amino acid substitutions on protein function and classify mutations as deleterious or benign. These include approaches that rely on evolutionary conservation, structural constraints, or physicochemical attributes of amino acid substitutions. Here we review existing methods and further examine eight tools: SIFT, PolyPhen2, Condel, CHASM, mCluster, logRE, SNAP, and MutationAssessor, with respect to their coverage, accuracy, availability and dependence on other tools.

Metagenomics projects based on shotgun sequencing of populations of micro-organisms yield insight into protein families. We used sequence similarity clustering to explore proteins with a comprehensive dataset consisting of sequences from available databases together with 6.12 million proteins predicted from an assembly of 7.7 million Global Ocean Sampling (GOS) sequences. The GOS dataset covers nearly all known prokaryotic protein families. A total of 3,995 medium- and large-sized clusters consisting of only GOS sequences are identified, out of which 1,700 have no detectable homology to known families. The GOS-only clusters contain a higher than expected proportion of sequences of viral origin, thus reflecting a poor sampling of viral diversity until now. Protein domain distributions in the GOS dataset and current protein databases show distinct biases. Several protein domains that were previously categorized as kingdom specific are shown to have GOS examples in other kingdoms. About 6,000 sequences (ORFans) from the literature that heretofore lacked similarity to known proteins have matches in the GOS data. The GOS dataset is also used to improve remote homology detection. Overall, besides nearly doubling the number of current proteins, the predicted GOS proteins also add a great deal of diversity to known protein families and shed light on their evolution. These observations are illustrated using several protein families, including phosphatases, proteases, ultraviolet-irradiation DNA damage repair enzymes, glutamine synthetase, and RuBisCO. The diversity added by GOS data has implications for choosing targets for experimental structure characterization as part of structural genomics efforts. Our analysis indicates that new families are being discovered at a rate that is linear or almost linear with the addition of new sequences, implying that we are still far from discovering all protein families in nature.

The eukaryotic protein kinase (ePK) domain mediates the majority of signaling and coordination of complex events in eukaryotes. By contrast, most bacterial signaling is thought to occur through structurally unrelated histidine kinases, though some ePK-like kinases (ELKs) and small molecule kinases are known in bacteria. Our analysis of the Global Ocean Sampling (GOS) dataset reveals that ELKs are as prevalent as histidine kinases and may play an equally important role in prokaryotic behavior. By combining GOS and public databases, we show that the ePK is just one subset of a diverse superfamily of enzymes built on a common protein kinase-like (PKL) fold. We explored this huge phylogenetic and functional space to cast light on the ancient evolution of this superfamily, its mechanistic core, and the structural basis for its observed diversity. We cataloged 27,677 ePKs and 18,699 ELKs, and classified them into 20 highly distinct families whose known members suggest regulatory functions. GOS data more than tripled the count of ELK sequences and enabled the discovery of novel families and classification and analysis of all ELKs. Comparison between and within families revealed ten key residues that are highly conserved across families. However, all but one of the ten residues has been eliminated in one family or another, indicating great functional plasticity. We show that loss of a catalytic lysine in two families is compensated by distinct mechanisms both involving other key motifs. This diverse superfamily serves as a model for further structural and functional analysis of enzyme evolution.

This paper reports a preliminary in silico analysis of the sea urchin kinome. The predicted protein kinases in the sea urchin genome were identified, annotated and classified, according to both function and kinase domain taxonomy. The results show that the sea urchin kinome, consisting of 353 protein kinases, is closer to the Drosophila kinome (239) than the human kinome (518) with respect to total kinase number. However, the diversity of sea urchin kinases is surprisingly similar to humans, since the urchin kinome is missing only 4 of 186 human subfamilies, while Drosophila lacks 24. Thus, the sea urchin kinome combines the simplicity of a non-duplicated genome with the diversity of function and signaling previously considered to be vertebrate-specific. More than half of the sea urchin kinases are involved with signal transduction, and approximately 88% of the signaling kinases are expressed in the developing embryo. These results support the strength of this nonchordate deuterostome as a pivotal developmental and evolutionary model organism.

Dictyostelium discoideum is a widely studied model organism with both unicellular and multicellular forms in its developmental cycle. The Dictyostelium genome encodes 285 predicted protein kinases, similar to the count of the much more advanced Drosophila. It contains members of most kinase classes shared by fungi and metazoans, as well as many previously thought to be metazoan specific, indicating that they have been secondarily lost from the fungal lineage. This includes the entire tyrosine kinase-like (TKL) group, which is expanded in Dictyostelium and includes several novel receptor kinases. Dictyostelium lacks tyrosine kinase group kinases, and most tyrosine phosphorylation appears to be mediated by TKL kinases. About half of Dictyostelium kinases occur in subfamilies not present in yeast or metazoa, suggesting that protein kinases have played key roles in the adaptation of Dictyostelium to its habitat. This study offers insights into kinase evolution and provides a focus for signaling analysis in this system.

The ciliate Tetrahymena thermophila is a model organism for molecular and cellular biology. Like other ciliates, this species has separate germline and soma functions that are embodied by distinct nuclei within a single cell. The germline-like micronucleus (MIC) has its genome held in reserve for sexual reproduction. The soma-like macronucleus (MAC), which possesses a genome processed from that of the MIC, is the center of gene expression and does not directly contribute DNA to sexual progeny. We report here the shotgun sequencing, assembly, and analysis of the MAC genome of T. thermophila, which is approximately 104 Mb in length and composed of approximately 225 chromosomes. Overall, the gene set is robust, with more than 27,000 predicted protein-coding genes, 15,000 of which have strong matches to genes in other organisms. The functional diversity encoded by these genes is substantial and reflects the complexity of processes required for a free-living, predatory, single-celled organism. This is highlighted by the abundance of lineage-specific duplications of genes with predicted roles in sensing and responding to environmental conditions (e.g., kinases), using diverse resources (e.g., proteases and transporters), and generating structural complexity (e.g., kinesins and dyneins). In contrast to the other lineages of alveolates (apicomplexans and dinoflagellates), no compelling evidence could be found for plastid-derived genes in the genome. UGA, the only T. thermophila stop codon, is used in some genes to encode selenocysteine, thus making this organism the first known with the potential to translate all 64 codons in nuclear genes into amino acids. We present genomic evidence supporting the hypothesis that the excision of DNA from the MIC to generate the MAC specifically targets foreign DNA as a form of genome self-defense. The combination of the genome sequence, the functional diversity encoded therein, and the presence of some pathways missing from other model organisms makes T. thermophila an ideal model for functional genomic studies to address biological, biomedical, and biotechnological questions of fundamental importance.

Approximately one in five breast cancers are driven by amplification and overexpression of the human epidermal growth factor receptor 2 (HER2) receptor kinase, and HER2-enriched (HER2E) is one of four major transcriptional subtypes of breast cancer. We set out to understand the genomics of HER2 amplification independent of subtype, and the underlying drivers and biology of HER2E tumors.