Alignment of Biological Sequences with Jalview.

In this chapter, we introduce core functionality of the Jalview interactive platform for the creation, analysis, and publication of multiple sequence alignments. A workflow is described based on Jalview's core functions: from data import to figure generation, including import of alignment reliability scores from T-Coffee and use of Jalview from the command line. The accompanying notes provide background information on the underlying methods and discuss additional options for working with Jalview to perform multiple sequence alignment, functional site analysis, and publication of alignments on the web.

Learning anatomy through Thiel- vs. formalin-embalmed cadavers: Student perceptions of embalming methods and effect on functional anatomy knowledge.

Thiel-embalmed cadavers, which have been adopted for use in anatomy teaching in relatively few universities, show greater flexibility and color retention compared to formalin-embalmed cadavers, properties which might be considered advantageous for anatomy teaching. This study aimed to investigate student attitudes toward the dissection experience with Thiel- compared to formalin/ethanol-embalmed cadavers. It also aimed to determine if one embalming method is more advantageous in terms of learning functional anatomy through the comparison of student anterior forearm functional anatomy knowledge. Student opinions and functional anatomy knowledge were obtained through use of a questionnaire from students at two medical schools, one using Thiel-, and one using more traditional formalin/ethanol-embalmed cadavers. Both the Thiel group and the formalin group of students were surveyed shortly after completing an anterior forearm dissection session. Significant differences (P-values <0.01) in some attitudes were found toward the dissection experience between cohorts using Thiel- vs. formalin-embalmed cadavers. The Thiel group of students felt more confident about recognizing anatomy in the living individual, found it easier to identify and dissect anatomical structures, and indicated more active exploration of functional anatomy due to the retained flexibility of the cadaver. However, on testing, no significant difference in functional anatomy knowledge was found between the two cohorts. Overall, although Thiel embalming may provide an advantageous learning experience in some investigated areas, more research needs to be carried out, especially to establish whether student perception is based on reality, at least in terms of structure identification. Anat Sci Educ 11: 166-174. © 2017 American Association of Anatomists.

Global quantitative SILAC phosphoproteomics reveals differential phosphorylation is widespread between the procyclic and bloodstream form lifecycle stages of Trypanosoma brucei.

We report a global quantitative phosphoproteomic study of bloodstream and procyclic form Trypanosoma brucei using SILAC labeling of each lifecycle stage. Phosphopeptide enrichment by SCX and TiO2 led to the identification of a total of 10096 phosphorylation sites on 2551 protein groups and quantified the ratios of 8275 phosphorylation sites between the two lifecycle stages. More than 9300 of these sites (92%) have not previously been reported. Model-based gene enrichment analysis identified over representation of Gene Ontology terms relating to the flagella, protein kinase activity, and the regulation of gene expression. The quantitative data reveal that differential protein phosphorylation is widespread between bloodstream and procyclic form trypanosomes, with significant intraprotein differential phosphorylation. Despite a lack of dedicated tyrosine kinases, 234 phosphotyrosine residues were identified, and these were 3-4 fold over-represented among site changing >10-fold between the two lifecycle stages. A significant proportion of the T. brucei kinome was phosphorylated, with evidence that MAPK pathways are functional in both lifecycle stages. Regulation of gene expression in T. brucei is exclusively post-transcriptional, and the extensive phosphorylation of RNA binding proteins observed may be relevant to the control of mRNA stability in this organism.

Tumorigenic fragments of APC cause dominant defects in directional cell migration in multiple model systems.

Nonsense mutations that result in the expression of truncated, N-terminal, fragments of the adenomatous polyposis coli (APC) tumour suppressor protein are found in most sporadic and some hereditary colorectal cancers. These mutations can cause tumorigenesis by eliminating ß-catenin-binding sites from APC, which leads to upregulation of ß-catenin and thereby results in the induction of oncogenes such as MYC. Here we show that, in three distinct experimental model systems, expression of an N-terminal fragment of APC (N-APC) results in loss of directionality, but not speed, of cell motility independently of changes in ß-catenin regulation. We developed a system to culture and fluorescently label live pieces of gut tissue to record high-resolution three-dimensional time-lapse movies of cells in situ. This revealed an unexpected complexity of normal gut cell migration, a key process in gut epithelial maintenance, with cells moving with spatial and temporal discontinuity. Quantitative comparison of gut tissue from wild-type mice and APC heterozygotes (APC(Min/+); multiple intestinal neoplasia model) demonstrated that cells in precancerous epithelia lack directional preference when moving along the crypt-villus axis. This effect was reproduced in diverse experimental systems: in developing chicken embryos, mesoderm cells expressing N-APC failed to migrate normally; in amoeboid Dictyostelium, which lack endogenous APC, expressing an N-APC fragment maintained cell motility, but the cells failed to perform directional chemotaxis; and multicellular Dictyostelium slug aggregates similarly failed to perform phototaxis. We propose that N-terminal fragments of APC represent a gain-of-function mutation that causes cells within tissue to fail to migrate directionally in response to relevant guidance cues. Consistent with this idea, crypts in histologically normal tissues of APC(Min/+) intestines are overpopulated with cells, suggesting that a lack of migration might cause cell accumulation in a precancerous state.

Genome sequence and analysis of the tuber crop potato.

Potato (Solanum tuberosum L.) is the world's most important non-grain food crop and is central to global food security. It is clonally propagated, highly heterozygous, autotetraploid, and suffers acute inbreeding depression. Here we use a homozygous doubled-monoploid potato clone to sequence and assemble 86% of the 844-megabase genome. We predict 39,031 protein-coding genes and present evidence for at least two genome duplication events indicative of a palaeopolyploid origin. As the first genome sequence of an asterid, the potato genome reveals 2,642 genes specific to this large angiosperm clade. We also sequenced a heterozygous diploid clone and show that gene presence/absence variants and other potentially deleterious mutations occur frequently and are a likely cause of inbreeding depression. Gene family expansion, tissue-specific expression and recruitment of genes to new pathways contributed to the evolution of tuber development. The potato genome sequence provides a platform for genetic improvement of this vital crop.

Prophossi: automating expert validation of phosphopeptide-spectrum matches from tandem mass spectrometry.

MOTIVATION: Complex patterns of protein phosphorylation mediate many cellular processes. Tandem mass spectrometry (MS/MS) is a powerful tool for identifying these post-translational modifications. In high-throughput experiments, mass spectrometry database search engines, such as MASCOT provide a ranked list of peptide identifications based on hundreds of thousands of MS/MS spectra obtained in a mass spectrometry experiment. These search results are not in themselves sufficient for confident assignment of phosphorylation sites as identification of characteristic mass differences requires time-consuming manual assessment of the spectra by an experienced analyst. The time required for manual assessment has previously rendered high-throughput confident assignment of phosphorylation sites challenging. RESULTS: We have developed a knowledge base of criteria, which replicate expert assessment, allowing more than half of cases to be automatically validated and site assignments verified with a high degree of confidence. This was assessed by comparing automated spectral interpretation with careful manual examination of the assignments for 501 peptides above the 1% false discovery rate (FDR) threshold corresponding to 259 putative phosphorylation sites in 74 proteins of the Trypanosoma brucei proteome. Despite this stringent approach, we are able to validate 80 of the 91 phosphorylation sites (88%) positively identified by manual examination of the spectra used for the MASCOT searches with a FDR < 15%. CONCLUSIONS: High-throughput computational analysis can provide a viable second stage validation of primary mass spectrometry database search results. Such validation gives rapid access to a systems level overview of protein phosphorylation in the experiment under investigation. AVAILABILITY: A GPL licensed software implementation in Perl for analysis and spectrum annotation is available in the supplementary material and a web server can be assessed online at http://www.compbio.dundee.ac.uk/prophossi.

The phosphoproteome of bloodstream form Trypanosoma brucei, causative agent of African sleeping sickness.

The protozoan parasite Trypanosoma brucei is the causative agent of human African sleeping sickness and related animal diseases, and it has over 170 predicted protein kinases. Protein phosphorylation is a key regulatory mechanism for cellular function that, thus far, has been studied in T.brucei principally through putative kinase mRNA knockdown and observation of the resulting phenotype. However, despite the relatively large kinome of this organism and the demonstrated essentiality of several T. brucei kinases, very few specific phosphorylation sites have been determined in this organism. Using a gel-free, phosphopeptide enrichment-based proteomics approach we performed the first large scale phosphorylation site analyses for T.brucei. Serine, threonine, and tyrosine phosphorylation sites were determined for a cytosolic protein fraction of the bloodstream form of the parasite, resulting in the identification of 491 phosphoproteins based on the identification of 852 unique phosphopeptides and 1204 phosphorylation sites. The phosphoproteins detected in this study are predicted from their genome annotations to participate in a wide variety of biological processes, including signal transduction, processing of DNA and RNA, protein synthesis, and degradation and to a minor extent in metabolic pathways. The analysis of phosphopeptides and phosphorylation sites was facilitated by in-house developed software, and this automated approach was validated by manual annotation of spectra of the kinase subset of proteins. Analysis of the cytosolic bloodstream form T. brucei kinome revealed the presence of 44 phosphorylated protein kinases in our data set that could be classified into the major eukaryotic protein kinase groups by applying a multilevel hidden Markov model library of the kinase catalytic domain. Identification of the kinase phosphorylation sites showed conserved phosphorylation sequence motifs in several kinase activation segments, supporting the view that phosphorylation-based signaling is a general and fundamental regulatory process that extends to this highly divergent lower eukaryote.

Jalview Version 2--a multiple sequence alignment editor and analysis workbench.

UNLABELLED: Jalview Version 2 is a system for interactive WYSIWYG editing, analysis and annotation of multiple sequence alignments. Core features include keyboard and mouse-based editing, multiple views and alignment overviews, and linked structure display with Jmol. Jalview 2 is available in two forms: a lightweight Java applet for use in web applications, and a powerful desktop application that employs web services for sequence alignment, secondary structure prediction and the retrieval of alignments, sequences, annotation and structures from public databases and any DAS 1.53 compliant sequence or annotation server. AVAILABILITY: The Jalview 2 Desktop application and JalviewLite applet are made freely available under the GPL, and can be downloaded from www.jalview.org.

Kinomer v. 1.0: a database of systematically classified eukaryotic protein kinases.

The regulation of protein function through reversible phosphorylation by protein kinases and phosphatases is a general mechanism controlling virtually every cellular activity. Eukaryotic protein kinases can be classified into distinct, well-characterized groups based on amino acid sequence similarity and function. We recently reported a highly sensitive and accurate hidden Markov model-based method for the automatic detection and classification of protein kinases into these specific groups. The Kinomer v. 1.0 database presented here contains annotated classifications for the protein kinase complements of 43 eukaryotic genomes. These span the taxonomic range and include fungi (16 species), plants (6), diatoms (1), amoebas (2), protists (1) and animals (17). The kinomes are stored in a relational database and are accessible through a web interface on the basis of species, kinase group or a combination of both. In addition, the Kinomer v. 1.0 HMM library is made available for users to perform classification on arbitrary sequences. The Kinomer v. 1.0 database is a continually updated resource where direct comparison of kinase sequences across kinase groups and across species can give insights into kinase function and evolution. Kinomer v. 1.0 is available at http://www.compbio.dundee.ac.uk/kinomer/.

A novel short splice variant of the tumour suppressor LKB1 is required for spermiogenesis.

LKB1 was discovered as a tumour suppressor mutated in Peutz-Jeghers syndrome, and is a gene involved in cell polarity as well as an upstream protein kinase for members of the AMP-activated protein kinase family. We report that mammals express two splice variants caused by alternate usage of 3'-exons. LKB1(L) is the previously described form, while LKB1(S) is a novel form in which the last 63 residues are replaced by a unique 39-residue sequence lacking known phosphorylation (Ser(431)) and farnesylation (Cys(433)) sites. Both isoforms are widely expressed in rodent and human tissues, although LKB1(S) is particularly abundant in haploid spermatids in the testis. Male mice in which expression of Lkb1(S) is knocked out are sterile, with the number of mature spermatozoa in the epididymis being dramatically reduced, and those spermatozoa that are produced have heads with an abnormal morphology and are non-motile. These results identify a previously undetected variant of LKB1, and suggest that it has a crucial role in spermiogenesis and male fertility.

TarO: a target optimisation system for structural biology.

TarO (http://www.compbio.dundee.ac.uk/taro) offers a single point of reference for key bioinformatics analyses relevant to selecting proteins or domains for study by structural biology techniques. The protein sequence is analysed by 17 algorithms and compared to 8 databases. TarO gathers putative homologues, including orthologues, and then obtains predictions of properties for these sequences including crystallisation propensity, protein disorder and post-translational modifications. Analyses are run on a high-performance computing cluster, the results integrated, stored in a database and accessed through a web-based user interface. Output is in tabulated format and in the form of an annotated multiple sequence alignment (MSA) that may be edited interactively in the program Jalview. TarO also simplifies the gathering of additional annotations via the Distributed Annotation System, both from the MSA in Jalview and through links to Dasty2. Routes to other information gateways are included, for example to relevant pages from UniProt, COG and the Conserved Domains Database. Open access to TarO is available from a guest account with private accounts for academic use available on request. Future development of TarO will include further analysis steps and integration with the Protein Information Management System (PIMS), a sister project in the BBSRC 'Structural Proteomics of Rational Targets' initiative.

Draft genome of the filarial nematode parasite Brugia malayi.

Parasitic nematodes that cause elephantiasis and river blindness threaten hundreds of millions of people in the developing world. We have sequenced the approximately 90 megabase (Mb) genome of the human filarial parasite Brugia malayi and predict approximately 11,500 protein coding genes in 71 Mb of robustly assembled sequence. Comparative analysis with the free-living, model nematode Caenorhabditis elegans revealed that, despite these genes having maintained little conservation of local synteny during approximately 350 million years of evolution, they largely remain in linkage on chromosomal units. More than 100 conserved operons were identified. Analysis of the predicted proteome provides evidence for adaptations of B. malayi to niches in its human and vector hosts and insights into the molecular basis of a mutualistic relationship with its Wolbachia endosymbiont. These findings offer a foundation for rational drug design.

Identification of multiple distinct Snf2 subfamilies with conserved structural motifs.

The Snf2 family of helicase-related proteins includes the catalytic subunits of ATP-dependent chromatin remodelling complexes found in all eukaryotes. These act to regulate the structure and dynamic properties of chromatin and so influence a broad range of nuclear processes. We have exploited progress in genome sequencing to assemble a comprehensive catalogue of over 1300 Snf2 family members. Multiple sequence alignment of the helicase-related regions enables 24 distinct subfamilies to be identified, a considerable expansion over earlier surveys. Where information is known, there is a good correlation between biological or biochemical function and these assignments, suggesting Snf2 family motor domains are tuned for specific tasks. Scanning of complete genomes reveals all eukaryotes contain members of multiple subfamilies, whereas they are less common and not ubiquitous in eubacteria or archaea. The large sample of Snf2 proteins enables additional distinguishing conserved sequence blocks within the helicase-like motor to be identified. The establishment of a phylogeny for Snf2 proteins provides an opportunity to make informed assignments of function, and the identification of conserved motifs provides a framework for understanding the mechanisms by which these proteins function.

Retrotransposon populations of Vicia species with varying genome size.

The (non-LTR) LINE and Ty3-gypsy-type LTR retrotransposon populations of three Vicia species that differ in genome size (Vicia faba, Vicia melanops and Vicia sativa) have been characterised. In each species the LINE retrotransposons comprise a complex, very heterogeneous set of sequences, while the Ty3-gypsy elements are much more homogeneous. Copy numbers of all three retrotransposon groups (Ty1-copia, Ty3-gypsy and LINE) in these species have been estimated by random genomic sequencing and Southern hybridisation analysis. The Ty3-gypsy elements are extremely numerous in all species, accounting for 18-35% of their genomes. The Ty1-copia group elements are somewhat less abundant and LINE elements are present in still lower amounts. Collectively, 20-45% of the genomes of these three Vicia species are comprised of retrotransposons. These data show that the three retrotransposon groups have proliferated to different extents in members of the Vicia genus and high proliferation has been associated with homogenisation of the retrotransposon population.

A preliminary crystallographic analysis of the putative mevalonate diphosphate decarboxylase from Trypanosoma brucei.

Mevalonate diphosphate decarboxylase catalyses the last and least well characterized step in the mevalonate pathway for the biosynthesis of isopentenyl pyrophosphate, an isoprenoid precursor. A gene predicted to encode the enzyme from Trypanosoma brucei has been cloned, a highly efficient expression system established and a purification protocol determined. The enzyme gives monoclinic crystals in space group P2(1), with unit-cell parameters a = 51.5, b = 168.7, c = 54.9 A, beta = 118.8 degrees. A Matthews coefficient VM of 2.5 A3 Da(-1) corresponds to two monomers, each approximately 42 kDa (385 residues), in the asymmetric unit with 50% solvent content. These crystals are well ordered and data to high resolution have been recorded using synchrotron radiation.

GOtcha: a new method for prediction of protein function assessed by the annotation of seven genomes.

BACKGROUND: The function of a novel gene product is typically predicted by transitive assignment of annotation from similar sequences. We describe a novel method, GOtcha, for predicting gene product function by annotation with Gene Ontology (GO) terms. GOtcha predicts GO term associations with term-specific probability (P-score) measures of confidence. Term-specific probabilities are a novel feature of GOtcha and allow the identification of conflicts or uncertainty in annotation. RESULTS: The GOtcha method was applied to the recently sequenced genome for Plasmodium falciparum and six other genomes. GOtcha was compared quantitatively for retrieval of assigned GO terms against direct transitive assignment from the highest scoring annotated BLAST search hit (TOPBLAST). GOtcha exploits information deep into the 'twilight zone' of similarity search matches, making use of much information that is otherwise discarded by more simplistic approaches. At a P-score cutoff of 50%, GOtcha provided 60% better recovery of annotation terms and 20% higher selectivity than annotation with TOPBLAST at an E-value cutoff of 10(-4). CONCLUSIONS: The GOtcha method is a useful tool for genome annotators. It has identified both errors and omissions in the original Plasmodium falciparum annotation and is being adopted by many other genome sequencing projects.

ELM server: A new resource for investigating short functional sites in modular eukaryotic proteins.

Multidomain proteins predominate in eukaryotic proteomes. Individual functions assigned to different sequence segments combine to create a complex function for the whole protein. While on-line resources are available for revealing globular domains in sequences, there has hitherto been no comprehensive collection of small functional sites/motifs comparable to the globular domain resources, yet these are as important for the function of multidomain proteins. Short linear peptide motifs are used for cell compartment targeting, protein-protein interaction, regulation by phosphorylation, acetylation, glycosylation and a host of other post-translational modifications. ELM, the Eukaryotic Linear Motif server at http://elm.eu.org/, is a new bioinformatics resource for investigating candidate short non-globular functional motifs in eukaryotic proteins, aiming to fill the void in bioinformatics tools. Sequence comparisons with short motifs are difficult to evaluate because the usual significance assessments are inappropriate. Therefore the server is implemented with several logical filters to eliminate false positives. Current filters are for cell compartment, globular domain clash and taxonomic range. In favourable cases, the filters can reduce the number of retained matches by an order of magnitude or more.

Visual representation of database search results: the RHIMS Plot.

SUMMARY: An algorithm and software are described that provide a fast method to produce a novel, function-oriented visualization of the results of a sequence database search. Text mining of sequence annotations allows position specific plots of potential functional similarity to be compared in a simple compact representation. AVAILABILITY: The application can be accessed via a web server at http://www.compbio.dundee.ac.uk. The RHIMS software may be obtained by request to the authors.

Genome sequence of the human malaria parasite Plasmodium falciparum.

The parasite Plasmodium falciparum is responsible for hundreds of millions of cases of malaria, and kills more than one million African children annually. Here we report an analysis of the genome sequence of P. falciparum clone 3D7. The 23-megabase nuclear genome consists of 14 chromosomes, encodes about 5,300 genes, and is the most (A + T)-rich genome sequenced to date. Genes involved in antigenic variation are concentrated in the subtelomeric regions of the chromosomes. Compared to the genomes of free-living eukaryotic microbes, the genome of this intracellular parasite encodes fewer enzymes and transporters, but a large proportion of genes are devoted to immune evasion and host-parasite interactions. Many nuclear-encoded proteins are targeted to the apicoplast, an organelle involved in fatty-acid and isoprenoid metabolism. The genome sequence provides the foundation for future studies of this organism, and is being exploited in the search for new drugs and vaccines to fight malaria.