ANNIE: integrated de novo protein sequence annotation

Identifying the causes and consequences of assembly gaps using a multiplatform genome assembly of a bird-of-paradise

Genome assemblies are currently being produced at an impressive rate by consortia and individual laboratories. The low costs and increasing efficiency of sequencing technologies now enable assembling genomes at unprecedented quality and contiguity. However, the difficulty in assembling repeat-rich and GC-rich regions (genomic "dark matter") limits insights into the evolution of genome structure and regulatory networks. Here, we compare the efficiency of currently available sequencing technologies (short/linked/long reads and proximity ligation maps) and combinations thereof in assembling genomic dark matter. By adopting different de novo assembly strategies, we compare individual draft assemblies to a curated multiplatform reference assembly and identify the genomic features that cause gaps within each assembly. We show that a multiplatform assembly implementing long-read, linked-read and proximity sequencing technologies performs best at recovering transposable elements, multicopy MHC genes, GC-rich microchromosomes and the repeat-rich W chromosome. Telomere-to-telomere assemblies are not a reality yet for most organisms, but by leveraging technology choice it is now possible to minimize genome assembly gaps for downstream analysis. We provide a roadmap to tailor sequencing projects for optimized completeness of both the coding and noncoding parts of nonmodel genomes.

Natural Defect of a Plant Rhabdovirus Glycoprotein Gene: A Case Study of Virus-Plant Coevolution

Seven isolates of a putative cytorhabdovirus (family Rhabdoviridae, order Mononegavirales) designated as citrus-associated rhabdovirus (CiaRV) were identified in citrus, passion fruit, and paper bush from the same geographical area in China. CiaRV, bean-associated cytorhabdovirus (Brazil), and papaya virus E (Ecuador) should be taxonomically classified in the species Papaya cytorhabdovirus. Due to natural mutations, the glycoprotein (G) and P4 genes were impaired in citrus-infecting isolates of CiaRV, resulting in an atypical rhabdovirus genome organization of 3' leader-N-P-P3-M-L-5' trailer. The P3 protein of CiaRV shared a common origin with begomoviral movement proteins (family Geminiviridae). Secondary structure analysis and trans-complementation of movement-deficient tomato mosaic virus and potato virus X mutants by CiaRV P3 supported its function in viral cell-to-cell trafficking. The wide geographical dispersal of CiaRV and related viruses suggests an efficient transmission mechanism, as well as an underlying risk to global agriculture. Both the natural phenomenon and experimental analyses demonstrated presence of the "degraded" type of CiaRV in citrus, in parallel to "undegraded" types in other host plant species. This case study shows a plant virus losing the function of an important but nonessential gene, likely due to host shift and adaption, which deepened our understanding of course of natural viral diversification.

Genome reconstruction of the non-culturable spinach downy mildew Peronospora effusa by metagenome filtering

Peronospora effusa (previously known as P. farinosa f. sp. spinaciae, and here referred to as Pfs) is an obligate biotrophic oomycete that causes downy mildew on spinach (Spinacia oleracea). To combat this destructive many disease resistant cultivars have been bred and used. However, new Pfs races rapidly break the employed resistance genes. To get insight into the gene repertoire of Pfs and identify infection-related genes, the genome of the first reference race, Pfs1, was sequenced, assembled, and annotated. Due to the obligate biotrophic nature of this pathogen, material for DNA isolation can only be collected from infected spinach leaves that, however, also contain many other microorganisms. The obtained sequences can, therefore, be considered a metagenome. To filter and obtain Pfs sequences we utilized the CAT tool to taxonomically annotate ORFs residing on long sequences of a genome pre-assembly. This study is the first to show that CAT filtering performs well on eukaryotic contigs. Based on the taxonomy, determined on multiple ORFs, contaminating long sequences and corresponding reads were removed from the metagenome. Filtered reads were re-assembled to provide a clean and improved Pfs genome sequence of 32.4 Mbp consisting of 8,635 scaffolds. Transcript sequencing of a range of infection time points aided the prediction of a total of 13,277 gene models, including 99 RxLR(-like) effector, and 14 putative Crinkler genes. Comparative analysis identified common features in the predicted secretomes of different obligate biotrophic oomycetes, regardless of their phylogenetic distance. Their secretomes are generally smaller, compared to hemi-biotrophic and necrotrophic oomycete species. We observe a reduction in proteins involved in cell wall degradation, in Nep1-like proteins (NLPs), proteins with PAN/apple domains, and host translocated effectors. The genome of Pfs1 will be instrumental in studying downy mildew virulence and for understanding the molecular adaptations by which new isolates break spinach resistance.

IMA Genome-F 11: Draft genome sequences of Fusarium xylarioides, Teratosphaeria gauchensis and T. zuluensis and genome annotation for Ceratocystis fimbriata

Draft genomes of the fungal species Fusarium xylarioides, Teratosphaeria gauchensis and T. zuluensis are presented. In addition an annotation of the genome of Ceratocystis fimbriata is presented. Overall these genomes provide a valuable resource for understanding the molecular processes underlying pathogenicity and potential management strategies of these economically important fungi.

Identification and Characterization of a Novel Emaravirus Associated With Jujube (Ziziphus jujuba Mill.) Yellow Mottle Disease

A previously unreported disease affecting jujube (Ziziphus jujuba Mill.) trees was observed in China (Liaoning province) in 2015 and named jujube yellow mottle disease (JYMD), due to prevalent symptoms on the leaves. Diseased plants produced also malformed and discolored fruits. In an attempt to identify the possible causal agent of JYMD, high-throughput sequencing of small RNA libraries was performed and a novel virus, tentatively named jujube yellow mottle-associated virus (JYMaV), was identified and characterized. Six genomic RNA segments of JYMaV were completely sequenced. Each one contains a single open reading frame in the viral complementary strand and two untranslated regions with complementary 5' and 3' terminal ends, thus showing typical features of other negative-stranded RNA viruses. RNA1 (7.1 kb), RNA2 (2.2 kb) and RNA3 (1.2 kb) encode putative proteins that, based on their conserved motifs, have been identified as the RNA dependent RNA polymerase, the glycoprotein and the nucleocapsid protein, respectively. These proteins share significant sequence identity (52.1-70.4%) with proteins encoded by raspberry leaf blotch virus (RLBV). RNA4 (1.5 kb) and RNA5 (1.2 kb) code for two putative 30 K movement proteins also related to the homologous RLBV protein. The functional role of the protein encoded by JYMaV RNA6 remains unknown. These data together with the phylogenetic relationships of JYMaV with other recognized emaraviruses support the proposal that JYMaV is the representative member of a novel species in the genus Emaravirus. In agreement with this proposal, virus-like particles and double-membrane-bound bodies, similar to those previously reported for other emaraviruses, were observed by transmission electron microscopy in extracts and tissues from symptomatic leaves, respectively. A specific RT-PCR-based detection method has been developed and used in a preliminary field survey that provided results strongly supporting the close association of JYMaV with the novel disease.

A DNA-Modified Live Vaccine Prime-Boost Strategy Broadens the T-Cell Response and Enhances the Antibody Response against the Porcine Reproductive and Respiratory Syndrome Virus

The Porcine Reproductive and Respiratory Syndrome Virus (PRRSV) induces reproductive disorders in sows and respiratory illnesses in growing pigs and is considered as one of the main pathogenic agents responsible for economic losses in the porcine industry worldwide. Modified live PRRSV vaccines (MLVs) are very effective vaccine types against homologous strains but they present only partial protection against heterologous viral variants. With the goal to induce broad and cross-protective immunity, we generated DNA vaccines encoding B and T antigens derived from a European subtype 1 strain that include T-cell epitope sequences known to be conserved across strains. These antigens were expressed either in a native form or in the form of vaccibodies targeted to the endocytic receptor XCR1 and CD11c expressed by different types of antigen-presenting cells (APCs). When delivered in skin with cationic nanoparticles and surface electroporation, multiple DNA vaccinations as a stand-alone regimen induced substantial antibody and T-cell responses, which were not promoted by targeting antigens to APCs. Interestingly, a DNA-MLV prime-boost strategy strongly enhanced the antibody response and broadened the T-cell responses over the one induced by MLV or DNA-only. The anti-nucleoprotein antibody response induced by the DNA-MLV prime-boost was clearly promoted by targeting the antigen to CD11c and XCR1, indicating a benefit of APC-targeting on the B-cell response. In conclusion, a DNA-MLV prime-boost strategy, by enhancing the potency and breadth of MLV vaccines, stands as a promising vaccine strategy to improve the control of PRRSV in infected herds.

Discovery of a novel splice variant of Fcar (CD89) unravels sequence segments necessary for efficient secretion: A story of bad signal peptides and good ones that nevertheless do not make it

The IgA receptor, Fcar (CD89) consists of 5 sequence segments: 2 segments (S1, S2) forming the potential signal peptide, 2 extracellular EC domains that include the IgA binding site, and the transmembrane and cytoplasmic tail (TM/C) region. Numerous Fcar splice variants have been reported with various combinations of the sequence segments mentioned above. Here, we report a novel splice variant termed variant APD isolated from a healthy volunteer that lacks only the IgA-binding EC1 domain. Despite possessing the complete signal peptide S1+S2, the variant APD is only found in the intracellular space whereas the wild-type variant 1 is efficiently secreted and variant 4 leaks to the extracellular space. Further mutational experiments involving signal peptide replacements, cleavage site modifications, and studies on alternative isoforms demonstrate that despite the completeness of the signal peptide motif, the presence of the EC1 domain is essential for efficient extracellular export.

A conserved region in the Closterovirus 1a polyprotein drives extensive remodeling of endoplasmic reticulum membranes and induces motile globules in Nicotiana benthamiana cells

In infected plant cells, closterovirus replicative polyproteins 1a and 1ab drive membrane remodeling and formation of multivesicular replication platforms. Polyprotein 1a contains a variable Central Region (CR) between the methyltransferase and helicase domains. In a previous study, we have found that transient expression of the Beet yellows virus CR-2 segment (aa 1305-1494) in Nicotiana benthamiana induces the formation of ~1µm mobile globules originating from the ER membranes. In the present study, sequence analysis has shown that a part of the CR named the "Zemlya region" (overlapping the CR-2), is conserved in all members of the Closterovirus genus and contains a predicted amphipathic helix (aa 1368-1385). By deletion analysis, the CR-2 region responsible for the induction of 1-μm globules has been mapped to aa 1368-1432. We suggest that the conserved membrane-modifying region of the BYV 1a may be involved in the biogenesis of closterovirus replication platforms.

Bioinformatic and mutational analysis of ophiovirus movement proteins, belonging to the 30K superfamily

Ophioviridae is a family of segmented, negative-sense, single-stranded RNA plant viruses. We showed that their cell-to-cell movement protein (MP) is an isolated member of the 30K MP superfamily with a unique structural organization. All 30K MPs share a core domain that contains a nearly-invariant signature aspartate. We examined its role in the MP of Citrus psorosis virus (CPsV) and Mirafiori lettuce big-vein virus (MiLBVV). Alanine substitution of this aspartate prevented plasmodesmata accumulation of MP(MiLBVV), while MP(CPsV) was not affected. The capacity of ophiovirus MPs to increase the plasmodesmata size exclusion limit and non-cell autonomous protein feature was abolished in both mutants. To investigate the role of the signature aspartate in cell-to-cell movement, we constructed a new movement-deficient Tobacco mosaic virus vector used for trans-complementation assays. We showed that both ophiovirus MP mutants lack the cell-to-cell movement capacity, confirming that this signature aspartate is essential for viral cell-to-cell movement.

Protein Sequence Annotation Tool (PSAT): a centralized web-based meta-server for high-throughput sequence annotations

Background

Here we introduce the Protein Sequence Annotation Tool (PSAT), a web-based, sequence annotation meta-server for performing integrated, high-throughput, genome-wide sequence analyses. Our goals in building PSAT were to (1) create an extensible platform for integration of multiple sequence-based bioinformatics tools, (2) enable functional annotations and enzyme predictions over large input protein fasta data sets, and (3) provide a web interface for convenient execution of the tools.

Results

In this paper, we demonstrate the utility of PSAT by annotating the predicted peptide gene products of Herbaspirillum sp. strain RV1423, importing the results of PSAT into EC2KEGG, and using the resulting functional comparisons to identify a putative catabolic pathway, thereby distinguishing RV1423 from a well annotated Herbaspirillum species. This analysis demonstrates that high-throughput enzyme predictions, provided by PSAT processing, can be used to identify metabolic potential in an otherwise poorly annotated genome.

Conclusions

PSAT is a meta server that combines the results from several sequence-based annotation and function prediction codes, and is available at http://psat.llnl.gov/psat/. PSAT stands apart from other sequence-based genome annotation systems in providing a high-throughput platform for rapid de novo enzyme predictions and sequence annotations over large input protein sequence data sets in FASTA. PSAT is most appropriately applied in annotation of large protein FASTA sets that may or may not be associated with a single genome.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-016-0887-y) contains supplementary material, which is available to authorized users.

HPMV: human protein mutation viewer - relating sequence mutations to protein sequence architecture and function changes

Next-generation sequencing advances are rapidly expanding the number of human mutations to be analyzed for causative roles in genetic disorders. Our Human Protein Mutation Viewer (HPMV) is intended to explore the biomolecular mechanistic significance of non-synonymous human mutations in protein-coding genomic regions. The tool helps to assess whether protein mutations affect the occurrence of sequence-architectural features (globular domains, targeting signals, post-translational modification sites, etc.). As input, HPMV accepts protein mutations - as UniProt accessions with mutations (e.g. HGVS nomenclature), genome coordinates, or FASTA sequences. As output, HPMV provides an interactive cartoon showing the mutations in relation to elements of the sequence architecture. A large variety of protein sequence architectural features were selected for their particular relevance to mutation interpretation. Clicking a sequence feature in the cartoon expands a tree view of additional information including multiple sequence alignments of conserved domains and a simple 3D viewer mapping the mutation to known PDB structures, if available. The cartoon is also correlated with a multiple sequence alignment of similar sequences from other organisms. In cases where a mutation is likely to have a straightforward interpretation (e.g. a point mutation disrupting a well-understood targeting signal), this interpretation is suggested. The interactive cartoon can be downloaded as standalone viewer in Java jar format to be saved and viewed later with only a standard Java runtime environment. The HPMV website is: http://hpmv.bii.a-star.edu.sg/ .

Genome Wide Re-Annotation of Caldicellulosiruptor saccharolyticus with New Insights into Genes Involved in Biomass Degradation and Hydrogen Production

Caldicellulosiruptor saccharolyticus has proven itself to be an excellent candidate for biological hydrogen (H₂) production, but still it has major drawbacks like sensitivity to high osmotic pressure and low volumetric H₂ productivity, which should be considered before it can be used industrially. A whole genome re-annotation work has been carried out as an attempt to update the incomplete genome information that causes gap in the knowledge especially in the area of metabolic engineering, to improve the H₂ producing capabilities of C. saccharolyticus. Whole genome re-annotation was performed through manual means for 2,682 Coding Sequences (CDSs). Bioinformatics tools based on sequence similarity, motif search, phylogenetic analysis and fold recognition were employed for re-annotation. Our methodology could successfully add functions for 409 hypothetical proteins (HPs), 46 proteins previously annotated as putative and assigned more accurate functions for the known protein sequences. Homology based gene annotation has been used as a standard method for assigning function to novel proteins, but over the past few years many non-homology based methods such as genomic context approaches for protein function prediction have been developed. Using non-homology based functional prediction methods, we were able to assign cellular processes or physical complexes for 249 hypothetical sequences. Our re-annotation pipeline highlights the addition of 231 new CDSs generated from MicroScope Platform, to the original genome with functional prediction for 49 of them. The re-annotation of HPs and new CDSs is stored in the relational database that is available on the MicroScope web-based platform. In parallel, a comparative genome analyses were performed among the members of genus Caldicellulosiruptor to understand the function and evolutionary processes. Further, with results from integrated re-annotation studies (homology and genomic context approach), we strongly suggest that Csac_0437 and Csac_0424 encode for glycoside hydrolases (GH) and are proposed to be involved in the decomposition of recalcitrant plant polysaccharides. Similarly, HPs: Csac_0732, Csac_1862, Csac_1294 and Csac_0668 are suggested to play a significant role in biohydrogen production. Function prediction of these HPs by using our integrated approach will considerably enhance the interpretation of large-scale experiments targeting this industrially important organism.

Sequence analysis reveals a conserved extension in the capping enzyme of the alphavirus supergroup, and a homologous domain in nodaviruses

Background

Members of the alphavirus supergroup include human pathogens such as chikungunya virus, hepatitis E virus and rubella virus. They encode a capping enzyme with methyltransferase-guanylyltransferase (MTase-GTase) activity, which is an attractive drug target owing to its unique mechanism. However, its experimental study has proven very difficult.

Results

We examined over 50 genera of viruses by sequence analyses. Earlier studies showed that the MTase-GTase contains a “Core” region conserved in sequence. We show that it is followed by a long extension, which we termed “Iceberg” region, whose secondary structure, but not sequence, is strikingly conserved throughout the alphavirus supergroup. Sequence analyses strongly suggest that the minimal capping domain corresponds to the Core and Iceberg regions combined, which is supported by earlier experimental data. The Iceberg region contains all known membrane association sites that contribute to the assembly of viral replication factories. We predict that it may also contain an overlooked, widely conserved membrane-binding amphipathic helix. Unexpectedly, we detected a sequence homolog of the alphavirus MTase-GTase in taxa related to nodaviruses and to chronic bee paralysis virus. The presence of a capping enzyme in nodaviruses is biologically consistent, since they have capped genomes but replicate in the cytoplasm, where no cellular capping enzyme is present. The putative MTase-GTase domain of nodaviruses also contains membrane-binding sites that may drive the assembly of viral replication factories, revealing an unsuspected parallel with the alphavirus supergroup.

Conclusions

Our work will guide the functional analysis of the alphaviral MTase-GTase and the production of domains for structure determination. The identification of a homologous domain in a simple model system, nodaviruses, which replicate in numerous eukaryotic cell systems (yeast, flies, worms, mammals, and plants), can further help crack the function and structure of the enzyme.

Reviewers

This article was reviewed by Valerian Dolja, Eugene Koonin and Sebastian Maurer-Stroh.

Electronic supplementary material

The online version of this article (doi:10.1186/s13062-015-0050-0) contains supplementary material, which is available to authorized users.

What's that gene (or protein)? Online resources for exploring functions of genes, transcripts, and proteins

The genomic era has enabled research projects that use approaches including genome-scale screens, microarray analysis, next-generation sequencing, and mass spectrometry-based proteomics to discover genes and proteins involved in biological processes. Such methods generate data sets of gene, transcript, or protein hits that researchers wish to explore to understand their properties and functions and thus their possible roles in biological systems of interest. Recent years have seen a profusion of Internet-based resources to aid this process. This review takes the viewpoint of the curious biologist wishing to explore the properties of protein-coding genes and their products, identified using genome-based technologies. Ten key questions are asked about each hit, addressing functions, phenotypes, expression, evolutionary conservation, disease association, protein structure, interactors, posttranslational modifications, and inhibitors. Answers are provided by presenting the latest publicly available resources, together with methods for hit-specific and data set-wide information retrieval, suited to any genome-based analytical technique and experimental species. The utility of these resources is demonstrated for 20 factors regulating cell proliferation. Results obtained using some of these are discussed in more depth using the p53 tumor suppressor as an example. This flexible and universally applicable approach for characterizing experimental hits helps researchers to maximize the potential of their projects for biological discovery.

On the necessity of dissecting sequence similarity scores into segment-specific contributions for inferring protein homology, function prediction and annotation

Background

Protein sequence similarities to any types of non-globular segments (coiled coils, low complexity regions, transmembrane regions, long loops, etc. where either positional sequence conservation is the result of a very simple, physically induced pattern or rather integral sequence properties are critical) are pertinent sources for mistaken homologies. Regretfully, these considerations regularly escape attention in large-scale annotation studies since, often, there is no substitute to manual handling of these cases. Quantitative criteria are required to suppress events of function annotation transfer as a result of false homology assignments.

Results

The sequence homology concept is based on the similarity comparison between the structural elements, the basic building blocks for conferring the overall fold of a protein. We propose to dissect the total similarity score into fold-critical and other, remaining contributions and suggest that, for a valid homology statement, the fold-relevant score contribution should at least be significant on its own. As part of the article, we provide the DissectHMMER software program for dissecting HMMER2/3 scores into segment-specific contributions. We show that DissectHMMER reproduces HMMER2/3 scores with sufficient accuracy and that it is useful in automated decisions about homology for instructive sequence examples. To generalize the dissection concept for cases without 3D structural information, we find that a dissection based on alignment quality is an appropriate surrogate. The approach was applied to a large-scale study of SMART and PFAM domains in the space of seed sequences and in the space of UniProt/SwissProt.

Conclusions

Sequence similarity core dissection with regard to fold-critical and other contributions systematically suppresses false hits and, additionally, recovers previously obscured homology relationships such as the one between aquaporins and formate/nitrite transporters that, so far, was only supported by structure comparison.

Transamidase subunit GAA1/GPAA1 is a M28 family metallo-peptide-synthetase that catalyzes the peptide bond formation between the substrate protein's omega-site and the GPI lipid anchor's phosphoethanolamine

The transamidase subunit GAA1/GPAA1 is predicted to be the enzyme that catalyzes the attachment of the glycosylphosphatidyl (GPI) lipid anchor to the carbonyl intermediate of the substrate protein at the ω-site. Its ~300-amino acid residue lumenal domain is a M28 family metallo-peptide-synthetase with an α/β hydrolase fold, including a central 8-strand β-sheet and a single metal (most likely zinc) ion coordinated by 3 conserved polar residues. Phosphoethanolamine is used as an adaptor to make the non-peptide GPI lipid anchor look chemically similar to the N terminus of a peptide.

Powerful sequence similarity search methods and in-depth manual analyses can identify remote homologs in many apparently "orphan" viral proteins

The genome sequences of new viruses often contain many "orphan" or "taxon-specific" proteins apparently lacking homologs. However, because viral proteins evolve very fast, commonly used sequence similarity detection methods such as BLAST may overlook homologs. We analyzed a data set of proteins from RNA viruses characterized as "genus specific" by BLAST. More powerful methods developed recently, such as HHblits or HHpred (available through web-based, user-friendly interfaces), could detect distant homologs of a quarter of these proteins, suggesting that these methods should be used to annotate viral genomes. In-depth manual analyses of a subset of the remaining sequences, guided by contextual information such as taxonomy, gene order, or domain cooccurrence, identified distant homologs of another third. Thus, a combination of powerful automated methods and manual analyses can uncover distant homologs of many proteins thought to be orphans. We expect these methodological results to be also applicable to cellular organisms, since they generally evolve much more slowly than RNA viruses. As an application, we reanalyzed the genome of a bee pathogen, Chronic bee paralysis virus (CBPV). We could identify homologs of most of its proteins thought to be orphans; in each case, identifying homologs provided functional clues. We discovered that CBPV encodes a domain homologous to the Alphavirus methyltransferase-guanylyltransferase; a putative membrane protein, SP24, with homologs in unrelated insect viruses and insect-transmitted plant viruses having different morphologies (cileviruses, higreviruses, blunerviruses, negeviruses); and a putative virion glycoprotein, ORF2, also found in negeviruses. SP24 and ORF2 are probably major structural components of the virions.

Experimental and bioinformatic evidence that raspberry leaf blotch emaravirus P4 is a movement protein of the 30K superfamily

Emaravirus is a recently described genus of negative-strand RNA plant viruses. Emaravirus P4 protein localizes to plasmodesmata, suggesting that it could be a viral movement protein (MP). In the current study, we showed that the P4 protein of raspberry leaf blotch emaravirus (RLBV) rescued the cell-to-cell movement of a defective potato virus X (PVX) that had a deletion mutation in the triple gene block 1 movement-associated protein. This demonstrated that RLBV P4 is a functional MP. Sequence analyses revealed that P4 is a distant member of the 30K superfamily of MPs. All MPs of this family contain two highly conserved regions predicted to form β-strands, namely β1 and β2. We explored by alanine mutagenesis the role of two residues of P4 (Ile106 and Asp127) located in each of these strands. We also made the equivalent substitutions in the 29K MP of tobacco rattle virus, another member of the 30K superfamily. All substitutions abolished the ability to complement PVX movement, except for the I106A substitution in the β1 region of P4. This region has been shown to mediate membrane association of 30K MPs; our results show that it is possible to make non-conservative substitutions of a well-conserved aliphatic residue within β1 without preventing the membrane association or movement function of P4.

Viral proteins originated de novo by overprinting can be identified by codon usage: application to the "gene nursery" of Deltaretroviruses

A well-known mechanism through which new protein-coding genes originate is by modification of pre-existing genes, e.g. by duplication or horizontal transfer. In contrast, many viruses generate protein-coding genes de novo, via the overprinting of a new reading frame onto an existing (“ancestral”) frame. This mechanism is thought to play an important role in viral pathogenicity, but has been poorly explored, perhaps because identifying the de novo frames is very challenging. Therefore, a new approach to detect them was needed. We assembled a reference set of overlapping genes for which we could reliably determine the ancestral frames, and found that their codon usage was significantly closer to that of the rest of the viral genome than the codon usage of de novo frames. Based on this observation, we designed a method that allowed the identification of de novo frames based on their codon usage with a very good specificity, but intermediate sensitivity. Using our method, we predicted that the Rex gene of deltaretroviruses has originated de novo by overprinting the Tax gene. Intriguingly, several genes in the same genomic region have also originated de novo and encode proteins that regulate the functions of Tax. Such “gene nurseries” may be common in viral genomes. Finally, our results confirm that the genomic GC content is not the only determinant of codon usage in viruses and suggest that a constraint linked to translation must influence codon usage.

How does novelty originate in nature? It is commonly thought that new genes are generated mainly by modifications of existing genes (the “tinkering” model). In contrast, we have shown recently that in viruses, numerous genes are generated entirely de novo (“from scratch”). The role of these genes remains underexplored, however, because they are difficult to identify. We have therefore developed a new method to detect genes originated de novo in viral genomes, based on the observation that each viral genome has a unique “signature”, which genes originated de novo do not share. We applied this method to analyze the genes of Human T-Lymphotropic Virus 1 (HTLV1), a relative of the HIV virus and also a major human pathogen that infects about twenty million people worldwide. The life cycle of HTLV1 is finely regulated – it can stay dormant for long periods and can provoke blood cancers (leukemias) after a very long incubation. We discovered that several of the genes of HTLV1 have originated de novo. These novel genes play a key role in regulating the life cycle of HTLV1, and presumably its pathogenicity. Our investigations suggest that such “gene nurseries” may be common in viruses.

Low-resolution structure of the soluble domain GPAA1 (yGPAA170-247) of the glycosylphosphatidylinositol transamidase subunit GPAA1 from Saccharomyces cerevisiae

The GPI (glycosylphosphatidylinositol) transamidase complex catalyses the attachment of GPI anchors to eukaryotic proteins in the lumen of ER (endoplasmic reticulum). The Saccharomyces cerevisiae GPI transamidase complex consists of the subunits yPIG-K (Gpi8p), yPIG-S (Gpi17p), yPIG-T (Gpi16p), yPIG-U (CDC91/GAB1) and yGPAA1. We present the production of the two recombinant proteins yGPAA1^70–247 and yGPAA1^70–339 of the luminal domain of S. cerevisiae GPAA1, covering the amino acids 70–247 and 70–339 respectively. The secondary structural content of the stable and monodisperse yGPAA1^70–247 has been determined to be 28% α-helix and 27% β-sheet. SAXS (small-angle X-ray scattering) data showed that yGPAA1^70–247 has an R_g (radius of gyration) of 2.72±0.025 nm and D_max (maximum dimension) of 9.14 nm. These data enabled the determination of the two domain low-resolution solution structure of yGPAA1^70–247. The large elliptical shape of yGPAA1^70–247 is connected via a short stalk to the smaller hook-like domain of 0.8 nm in length and 3.5 nm in width. The topological arrangement of yGPAA1^70–247 will be discussed together with the recently determined low-resolution structures of yPIG-K^24–337 and yPIG-S^38–467 from S. cerevisiae in the GPI transamidase complex.

Whole exome sequencing identifies a mutation for a novel form of corneal intraepithelial dyskeratosis

BACKGROUND

Corneal intraepithelial dyskeratosis is an extremely rare condition. The classical form, affecting Native American Haliwa-Saponi tribe members, is called hereditary benign intraepithelial dyskeratosis (HBID). Herein, we present a new form of corneal intraepithelial dyskeratosis for which we identified the causative gene by using deep sequencing technology.

METHODS AND RESULTS

A seven member Caucasian French family with two corneal intraepithelial dyskeratosis affected individuals (6-year-old proband and his mother) was ascertained. The proband presented with bilateral complete corneal opacification and dyskeratosis. Palmoplantar hyperkeratosis and laryngeal dyskeratosis were associated with the phenotype. Histopathology studies of cornea and vocal cord biopsies showed dyskeratotic keratinisation. Quantitative PCR ruled out 4q35 duplication, classically described in HBID cases. Next generation sequencing with mean coverage of 50× using the Illumina Hi Seq and whole exome capture processing was performed. Sequence reads were aligned, and screened for single nucleotide variants and insertion/deletion calls. In-house pipeline filtering analyses and comparisons with available databases were performed. A novel missense mutation M77T was discovered for the gene NLRP1 which maps to chromosome 17p13.2. This was a de novo mutation in the proband's mother, following segregation in the family, and not found in 738 control DNA samples. NLRP1 expression was determined in adult corneal epithelium. The amino acid change was found to destabilise significantly the protein structure.

CONCLUSIONS

We describe a new corneal intraepithelial dyskeratosis and how we identified its causative gene. The NLRP1 gene product is implicated in inflammation, autoimmune disorders, and caspase mediated apoptosis. NLRP1 polymorphisms are associated with various diseases.

Beware of moving targets: reference proteome content fluctuates substantially over the years

Reference proteomes are generated by increasingly sophisticated annotation pipelines as part of regular genome build releases; yet, the corresponding changes in reference proteomes' content are dramatic. In the history of the NCBI-curated human proteome, the total number of entries has remained roughly constant but approximately half of the proteins from the 2003 build 33 are no longer represented by entries in current releases, while about the same number of new proteins have been added (for sequence identity thresholds 50-90%). Although mostly hypothetical proteins are affected, there are also spectacular cases of entry removal/addition of well studied proteins. The changes between the 2003 and recent human proteomes are in a similar order of magnitude as the differences between recent human and chimpanzee proteome releases. As an application example, we show that the proteome fluctuations affect the interpretation (about 74% of hits) of organelle-specific mass-spectrometry data. Although proteome quality tends to improve with more recent releases as, for example, the fraction of proteins with functional annotation has increased over time, existing evidence implies that, apparently, the proteome content still remains incomplete, not just pertaining to isoforms/sequence variants but also to proteins and their families that are clearly distinct.

A DECADE AFTER THE FIRST FULL HUMAN GENOME SEQUENCING: WHEN WILL WE UNDERSTAND OUR OWN GENOME?

The contrast between the pomp of celebrating the first full human genome sequencing in 2000 and the cautious tone of recollections a decade thereafter could hardly be greater. The promises with regard to medical cures and biotechnology applications have been realized not even nearly to the expectations. Understanding the human genomes means knowing the genes' and proteins' functions and their interconnectedness via biomolecular mechanisms. This articles estimates how long will it take to achieve this goal if we extrapolate from the previous decade (indeed, a century!) and the possible disruptive trends in science, technology and society that may accelerate the pace of progress dramatically.

Evolution of viral proteins originated de novo by overprinting

New protein-coding genes can originate either through modification of existing genes or de novo. Recently, the importance of de novo origination has been recognized in eukaryotes, although eukaryotic genes originated de novo are relatively rare and difficult to identify. In contrast, viruses contain many de novo genes, namely those in which an existing gene has been “overprinted” by a new open reading frame, a process that generates a new protein-coding gene overlapping the ancestral gene. We analyzed the evolution of 12 experimentally validated viral genes that originated de novo and estimated their relative ages. We found that young de novo genes have a different codon usage from the rest of the genome. They evolve rapidly and are under positive or weak purifying selection. Thus, young de novo genes might have strain-specific functions, or no function, and would be difficult to detect using current genome annotation methods that rely on the sequence signature of purifying selection. In contrast to young de novo genes, older de novo genes have a codon usage that is similar to the rest of the genome. They evolve slowly and are under stronger purifying selection. Some of the oldest de novo genes evolve under stronger selection pressure than the ancestral gene they overlap, suggesting an evolutionary tug of war between the ancestral and the de novo gene.

Tachyon search speeds up retrieval of similar sequences by several orders of magnitude

Summary: The usage of current sequence search tools becomes increasingly slower as databases of protein sequences continue to grow exponentially. Tachyon, a new algorithm that identifies closely related protein sequences ~200 times faster than standard BLAST, circumvents this limitation with a reduced database and oligopeptide matching heuristic.

Availability and implementation: The tool is publicly accessible as a webserver at http://tachyon.bii.a-star.edu.sg and can also be accessed programmatically through SOAP.

Contact:sebastianms@bii.a-star.edu.sg

Supplementary information:Supplementary data are available at the Bioinformatics online.

Genome-wide analysis of cell wall-related genes in Tuber melanosporum

A genome-wide inventory of proteins involved in cell wall synthesis and remodeling has been obtained by taking advantage of the recently released genome sequence of the ectomycorrhizal Tuber melanosporum black truffle. Genes that encode cell wall biosynthetic enzymes, enzymes involved in cell wall polysaccharide synthesis or modification, GPI-anchored proteins and other cell wall proteins were identified in the black truffle genome. As a second step, array data were validated and the symbiotic stage was chosen as the main focus. Quantitative RT-PCR experiments were performed on 29 selected genes to verify their expression during ectomycorrhizal formation. The results confirmed the array data, and this suggests that cell wall-related genes are required for morphogenetic transition from mycelium growth to the ectomycorrhizal branched hyphae. Labeling experiments were also performed on T. melanosporum mycelium and ectomycorrhizae to localize cell wall components.

Mutations in ROGDI Cause Kohlschütter-Tönz Syndrome

Kohlschütter-Tönz syndrome (KTS) is an autosomal-recessive disease characterized by the combination of epilepsy, psychomotor regression, and amelogenesis imperfecta. The molecular basis has not yet been elucidated. Here, we report that KTS is caused by mutations in ROGDI. Using a combination of autozygosity mapping and exome sequencing, we identified a homozygous frameshift deletion, c.229_230del (p.Leu77Alafs(∗)64), in ROGDI in two affected individuals from a consanguineous family. Molecular studies in two additional KTS-affected individuals from two unrelated Austrian and Swiss families revealed homozygosity for nonsense mutation c.286C>T (p.Gln96(∗)) and compound heterozygosity for the splice-site mutations c.531+5G>C and c.532-2A>T in ROGDI, respectively. The latter mutation was also found to be heterozygous in the mother of the Swiss affected individual in whom KTS was reported for the first time in 1974. ROGDI is highly expressed throughout the brain and other organs, but its function is largely unknown. Possible interactions with DISC1, a protein involved in diverse cytoskeletal functions, have been suggested. Our finding that ROGDI mutations cause KTS indicates that the protein product of this gene plays an important role in neuronal development as well as amelogenesis.

Nuclear import of a lipid-modified transcription factor: mobilization of NFAT5 isoform a by osmotic stress

Lipid-modified transcription factors (TFs) are biomolecular oddities since their reduced mobility and membrane attachment appear to contradict nuclear import required for their gene-regulatory function. NFAT5 isoform a (selected from an in silico screen for predicted lipid-modified TFs) is shown to contribute about half of all endogenous expression of human NFAT5 isoforms in the isotonic state. Wild-type NFAT5a protein is indeed myristoylated and palmitoylated on its transport to the plasmalemma via the endoplasmic reticulum and the Golgi. In contrast, its lipid anchor-deficient mutants as well as isoforms NFAT5b/c are diffusely localized in the cytoplasm without preference to vesicular structures. Quantitative/live microscopy shows the plasmamembrane-bound fraction of NFAT5a moving into the nucleus upon osmotic stress despite the lipid anchoring. The mobilization mechanism is not based on proteolytic processing of the lipid-anchored N-terminus but appears to involve reversible palmitoylation. Thus, NFAT5a is an example of TFs immobilized with lipid anchors at cyotoplasmic membranes in the resting state and that, nevertheless, can translocate into the nucleus upon signal induction.

Computational approaches to selecting and optimising targets for structural biology

Selection of protein targets for study is central to structural biology and may be influenced by numerous factors. A key aim is to maximise returns for effort invested by identifying proteins with the balance of biophysical properties that are conducive to success at all stages (e.g. solubility, crystallisation) in the route towards a high resolution structural model. Selected targets can be optimised through construct design (e.g. to minimise protein disorder), switching to a homologous protein, and selection of experimental methodology (e.g. choice of expression system) to prime for efficient progress through the structural proteomics pipeline.

Here we discuss computational techniques in target selection and optimisation, with more detailed focus on tools developed within the Scottish Structural Proteomics Facility (SSPF); namely XANNpred, ParCrys, OB-Score (target selection) and TarO (target optimisation). TarO runs a large number of algorithms, searching for homologues and annotating the pool of possible alternative targets. This pool of putative homologues is presented in a ranked, tabulated format and results are also visualised as an automatically generated and annotated multiple sequence alignment. The target selection algorithms each predict the propensity of a selected protein target to progress through the experimental stages leading to diffracting crystals. This single predictor approach has advantages for target selection, when compared with an approach using two or more predictors that each predict for success at a single experimental stage. The tools described here helped SSPF achieve a high (21%) success rate in progressing cloned targets to diffraction-quality crystals.

Purification and crystallization of yeast glycosylphosphatidylinositol transamidase subunit PIG-S (PIG-S(71-467))

The transfer of glycosylphosphatidylinositol (GPI) anchors onto eukaryotic proteins is catalyzed by the transamidase complex, which is composed of at least five subunits (PIG-K, PIG-S, PIG-T, PIG-U and GPAA1). Here, the recombinant protein PIG-S(71-467) from Saccharomyces cerevisiae, including residues 71-467 of the entire 534-residue protein, was cloned, expressed and purified to homogeneity. The monodisperse protein was crystallized by the vapour-diffusion method. A diffraction data set was collected to 3.2 Å resolution with 91.6% completeness. The crystals belonged to space group C2, with unit-cell parameters a = 106.72, b = 59.33, c = 124.3 Å, β = 114.19°, and contained two molecules in the asymmetric unit.

F-BAR domain proteins: Families and function

The F-BAR domain is emerging as an important player in membrane remodeling pathways. F-BAR domain proteins couple membrane remodeling with actin dynamics associated with endocytic pathways and filopodium formation. Here, we provide a comprehensive analysis of F-BAR domain proteins in terms of their evolutionary relationships and protein function. F-BAR domain containing proteins can be categorized into five subfamilies based on their phylogeny which is consistent with the additional protein domains they possess, for example, RhoGAP domains, Cdc42 binding sites, SH3 domains and tyrosine kinase domains. We derive a protein-protein interaction network suggesting that dynamin1/2, N-WASP, Huntingtin, intersectin and Cdc42 are central nodes influencing F-BAR domain protein function.

Reconstruction of gene association network reveals a transmembrane protein required for adipogenesis and targeted by PPARγ

We have developed a method for reconstructing gene association networks and have applied this method to gene profiles from 3T3-L1 cells. Priorization of the candidate genes pinpointed a transcript annotated as APMAP (adipocyte plasma membrane-associated protein). Functional studies showed that APMAP is upregulated in murine and human adipogenic cell models as well as in a genetic mouse model of obesity. Silencing APMAP in 3T3-L1 cells strongly impaired the differentiation into adipocytes. Moreover, APMAP expression was strongly induced by the PPARγ ligand rosiglitazone in adipocytes in vitro and in vivo in adipose tissue. Using ChIP-qPCR and luciferase reporter assays, we show a functional PPARγ binding site. In addition, we provide evidence that the extracellular C-terminal domain of APMAP is required for the function of APMAP in adipocyte differentiation. Finally, we demonstrate that APMAP translocates from the endoplasmatic reticulum to the plasma membrane during adipocyte differentiation.

More than 1,001 problems with protein domain databases: transmembrane regions, signal peptides and the issue of sequence homology

Large-scale genome sequencing gained general importance for life science because functional annotation of otherwise experimentally uncharacterized sequences is made possible by the theory of biomolecular sequence homology. Historically, the paradigm of similarity of protein sequences implying common structure, function and ancestry was generalized based on studies of globular domains. Having the same fold imposes strict conditions over the packing in the hydrophobic core requiring similarity of hydrophobic patterns. The implications of sequence similarity among non-globular protein segments have not been studied to the same extent; nevertheless, homology considerations are silently extended for them. This appears especially detrimental in the case of transmembrane helices (TMs) and signal peptides (SPs) where sequence similarity is necessarily a consequence of physical requirements rather than common ancestry. Thus, matching of SPs/TMs creates the illusion of matching hydrophobic cores. Therefore, inclusion of SPs/TMs into domain models can give rise to wrong annotations. More than 1001 domains among the 10,340 models of Pfam release 23 and 18 domains of SMART version 6 (out of 809) contain SP/TM regions. As expected, fragment-mode HMM searches generate promiscuous hits limited to solely the SP/TM part among clearly unrelated proteins. More worryingly, we show explicit examples that the scores of clearly false-positive hits, even in global-mode searches, can be elevated into the significance range just by matching the hydrophobic runs. In the PIR iProClass database v3.74 using conservative criteria, we find that at least between 2.1% and 13.6% of its annotated Pfam hits appear unjustified for a set of validated domain models. Thus, false-positive domain hits enforced by SP/TM regions can lead to dramatic annotation errors where the hit has nothing in common with the problematic domain model except the SP/TM region itself. We suggest a workflow of flagging problematic hits arising from SP/TM-containing models for critical reconsideration by annotation users.

Sequence homology is a fundamental principle of biology. It implies common phylogenetic ancestry of genes and, subsequently, similarity of their protein products with regard to amino acid sequence, three-dimensional structure and molecular and cellular function. Originally an esoteric concept, homology with the proxy of sequence similarity is used to justify the transfer of functional annotation from well-studied protein examples to new sequences. Yet, functional annotation via sequence similarity seems to have hit a plateau in recent years since relentless annotation transfer led to error propagation across sequence databases; thus, leading experimental follow-up work astray. It must be emphasized that the trinity of sequence, 3D structural and functional similarity has only been proven for globular segments of proteins. For non-globular regions, similarity of sequence is not necessarily a result of divergent evolution from a common ancestor but the consequence of amino acid sequence bias. In our investigation, we found that protein domain databases contain many domain models with transmembrane regions and signal peptides, non-globular segments of proteins having hydrophobic bias. Many proteins have inherited completely wrong function assignments from these domain models. We fear that future function predictions will turn out futile if this issue is not immediately addressed.

Is LGI2 the candidate gene for partial epilepsy with pericentral spikes?

Partial epilepsy with pericentral spikes (PEPS) is a familial epilepsy with disease locus mapped to human chromosome region 4p15; yet, the causative gene is unknown. In this work, arguments based on protein sequence analysis and patient-specific chromosomal deletions are provided for LGI2 as the prime candidate gene for PEPS among the 52 genes known at the genome locus 4p15. Furthermore, we suggest that two reports of patients that were not classified as PEPS but show very similar phenotypes and deletions in the PEPS disease locus, could in fact describe the same disease. To test this hypothesis, patients with diagnosed PEPS or the described similar phenotypes could be screened for mutations in LGI2 and other shortlisted candidate genes. The linkage between PEPS and its disease causing gene(s) would allow diagnosis of the disease based on genetic screening as well as hereditary studies. Furthermore, previous knowledge on molecular disease mechanisms of related LGI proteins, for example LGI1 and autosomal dominant lateral temporal epilepsy, could be applied to deepen the understanding of the PEPS disease mechanism at the molecular level, which may facilitate therapeutic intervention in the future. Supplementary Table is available at http://www.worldscinet.com/jbcb/.

Parameterization of disorder predictors for large-scale applications requiring high specificity by using an extended benchmark dataset

BACKGROUND

Algorithms designed to predict protein disorder play an important role in structural and functional genomics, as disordered regions have been reported to participate in important cellular processes. Consequently, several methods with different underlying principles for disorder prediction have been independently developed by various groups. For assessing their usability in automated workflows, we are interested in identifying parameter settings and threshold selections, under which the performance of these predictors becomes directly comparable.

RESULTS

First, we derived a new benchmark set that accounts for different flavours of disorder complemented with a similar amount of order annotation derived for the same protein set. We show that, using the recommended default parameters, the programs tested are producing a wide range of predictions at different levels of specificity and sensitivity. We identify settings, in which the different predictors have the same false positive rate. We assess conditions when sets of predictors can be run together to derive consensus or complementary predictions. This is useful in the framework of proteome-wide applications where high specificity is required such as in our in-house sequence analysis pipeline and the ANNIE webserver.

CONCLUSIONS

This work identifies parameter settings and thresholds for a selection of disorder predictors to produce comparable results at a desired level of specificity over a newly derived benchmark dataset that accounts equally for ordered and disordered regions of different lengths.

Mapping the sequence mutations of the 2009 H1N1 influenza A virus neuraminidase relative to drug and antibody binding sites

In this work, we study the consequences of sequence variations of the "2009 H1N1" (swine or Mexican flu) influenza A virus strain neuraminidase for drug treatment and vaccination. We find that it is phylogenetically more closely related to European H1N1 swine flu and H5N1 avian flu rather than to the H1N1 counterparts in the Americas. Homology-based 3D structure modeling reveals that the novel mutations are preferentially located at the protein surface and do not interfere with the active site. The latter is the binding cavity for 3 currently used neuraminidase inhibitors: oseltamivir (Tamiflu^®), zanamivir (Relenza^®) and peramivir; thus, the drugs should remain effective for treatment. However, the antigenic regions of the neuraminidase relevant for vaccine development, serological typing and passive antibody treatment can differ from those of previous strains and already vary among patients.

This article was reviewed by Sandor Pongor and L. Aravind.