Actomyosin organelle functions of SPIRE actin nucleators precede animal evolution.

An important question in cell biology is how cytoskeletal proteins evolved and drove the development of novel structures and functions. Here we address the origin of SPIRE actin nucleators. Mammalian SPIREs work with RAB GTPases, formin (FMN)-subgroup actin assembly proteins and class-5 myosin (MYO5) motors to transport organelles along actin filaments towards the cell membrane. However, the origin and extent of functional conservation of SPIRE among species is unknown. Our sequence searches show that SPIRE exist throughout holozoans (animals and their closest single-celled relatives), but not other eukaryotes. SPIRE from unicellular holozoans (choanoflagellate), interacts with RAB, FMN and MYO5 proteins, nucleates actin filaments and complements mammalian SPIRE function in organelle transport. Meanwhile SPIRE and MYO5 proteins colocalise to organelles in Salpingoeca rosetta choanoflagellates. Based on these observations we propose that SPIRE originated in unicellular ancestors of animals providing an actin-myosin driven exocytic transport mechanism that may have contributed to the evolution of complex multicellular animals.

Barley Ror1 encodes a class XI myosin required for mlo-based broad-spectrum resistance to the fungal powdery mildew pathogen.

Loss-of-function alleles of plant MLO genes confer broad-spectrum resistance to powdery mildews in many eudicot and monocot species. Although barley (Hordeum vulgare) mlo mutants have been used in agriculture for more than 40 years, understanding of the molecular principles underlying this type of disease resistance remains fragmentary. Forward genetic screens in barley have revealed mutations in two Required for mlo resistance (Ror) genes that partially impair immunity conferred by mlo mutants. While Ror2 encodes a soluble N-ethylmaleimide-sensitive factor-attached protein receptor (SNARE), the identity of Ror1, located at the pericentromeric region of barley chromosome 1H, remained elusive. We report the identification of Ror1 based on combined barley genomic sequence information and transcriptomic data from ror1 mutant plants. Ror1 encodes the barley class XI myosin Myo11A (HORVU.MOREX.r3.1HG0046420). Single amino acid substitutions of this myosin, deduced from non-functional ror1 mutant alleles, map to the nucleotide-binding region and the interface between the relay-helix and the converter domain of the motor protein. Ror1 myosin accumulates transiently in the course of powdery mildew infection. Functional fluorophore-labeled Ror1 variants associate with mobile intracellular compartments that partially colocalize with peroxisomes. Single-cell expression of the Ror1 tail region causes a dominant-negative effect that phenocopies ror1 loss-of-function mutants. We define a myosin motor for the establishment of mlo-mediated resistance, suggesting that motor protein-driven intracellular transport processes are critical for extracellular immunity, possibly through the targeted transfer of antifungal and/or cell wall cargoes to pathogen contact sites.

Design of typical genes for heterologous gene expression.

Heterologous protein expression is an important method for analysing cellular functions of proteins, in genetic circuit engineering and in overexpressing proteins for biopharmaceutical applications and structural biology research. The degeneracy of the genetic code, which enables a single protein to be encoded by a multitude of synonymous gene sequences, plays an important role in regulating protein expression, but substantial uncertainty exists concerning the details of this phenomenon. Here we analyse the influence of a profiled codon usage adaptation approach on protein expression levels in the eukaryotic model organism Saccharomyces cerevisiae. We selected green fluorescent protein (GFP) and human α-synuclein (αSyn) as representatives for stable and intrinsically disordered proteins and representing a benchmark and a challenging test case. A new approach was implemented to design typical genes resembling the codon usage of any subset of endogenous genes. Using this approach, synthetic genes for GFP and αSyn were generated, heterologously expressed and evaluated in yeast. We demonstrate that GFP is expressed at high levels, and that the toxic αSyn can be adapted to endogenous, low-level expression. The new software is publicly available as a web-application for performing host-specific protein adaptations to a set of the most commonly used model organisms ( https://odysseus.motorprotein.de ).

Proteogenomics analysis of CUG codon translation in the human pathogen Candida albicans.

BACKGROUND: Yeasts of the CTG-clade lineage, which includes the human-infecting Candida albicans, Candida parapsilosis and Candida tropicalis species, are characterized by an altered genetic code. Instead of translating CUG codons as leucine, as happens in most eukaryotes, these yeasts, whose ancestors are thought to have lost the relevant leucine-tRNA gene, translate CUG codons as serine using a serine-tRNA with a mutated anticodon, [Formula: see text]. Previously reported experiments have suggested that 3-5% of the CTG-clade CUG codons are mistranslated as leucine due to mischarging of the [Formula: see text]. The mistranslation was suggested to result in variable surface proteins explaining fast host adaptation and pathogenicity. RESULTS: In this study, we reassess this potential mistranslation by high-resolution mass spectrometry-based proteogenomics of multiple CTG-clade yeasts, including various C. albicans strains, isolated from colonized and from infected human body sites, and C. albicans grown in yeast and hyphal forms. Our data do not support a bias towards CUG codon mistranslation as leucine. Instead, our data suggest that (i) CUG codons are mistranslated at a frequency corresponding to the normal extent of ribosomal mistranslation with no preference for specific amino acids, (ii) CUG codons are as unambiguous (or ambiguous) as the related CUU leucine and UCC serine codons, (iii) tRNA anticodon loop variation across the CTG-clade yeasts does not result in any difference of the mistranslation level, and (iv) CUG codon unambiguity is independent of C. albicans' strain pathogenicity or growth form. CONCLUSIONS: Our findings imply that C. albicans does not decode CUG ambiguously. This suggests that the proposed misleucylation of the [Formula: see text] might be as prevalent as every other misacylation or mistranslation event and, if at all, be just one of many reasons causing phenotypic diversity.

Critical assessment of coiled-coil predictions based on protein structure data.

Coiled-coil regions were among the first protein motifs described structurally and theoretically. The simplicity of the motif promises that coiled-coil regions can be detected with reasonable accuracy and precision in any protein sequence. Here, we re-evaluated the most commonly used coiled-coil prediction tools with respect to the most comprehensive reference data set available, the entire Protein Data Bank, down to each amino acid and its secondary structure. Apart from the 30-fold difference in minimum and maximum number of coiled coils predicted the tools strongly vary in where they predict coiled-coil regions. Accordingly, there is a high number of false predictions and missed, true coiled-coil regions. The evaluation of the binary classification metrics in comparison with naïve coin-flip models and the calculation of the Matthews correlation coefficient, the most reliable performance metric for imbalanced data sets, suggests that the tested tools' performance is close to random. This implicates that the tools' predictions have only limited informative value. Coiled-coil predictions are often used to interpret biochemical data and are part of in-silico functional genome annotation. Our results indicate that these predictions should be treated very cautiously and need to be supported and validated by experimental evidence.

The Protein-Coding Human Genome: Annotating High-Hanging Fruits.

The major transcript variants of human protein-coding genes are annotated to a certain degree of accuracy combining manual curation, transcript data, and proteomics evidence. However, there is considerable disagreement on the annotation of about 2000 genes-they can be protein-coding, noncoding, or pseudogenes-and on the annotation of most of the predicted alternative transcripts. Pure transcriptome mapping approaches seem to be limited in discriminating functional expression from noise. These limitations have partially been overcome by dedicated algorithms to detect alternative spliced micro-exons and wobble splice variants. Recently, knowledge about splice mechanism and protein structure are incorporated into an algorithm to predict neighboring homologous exons, often spliced in a mutually exclusive manner. Predicted exons are evaluated by transcript data, structural compatibility, and evolutionary conservation, revealing hundreds of novel coding exons and splice mechanism re-assignments. The emerging human pan-genome is necessitating distinctive annotations incorporating differences between individuals and between populations.

Predicting Genes in Closely Related Species with Scipio and WebScipio.

Scipio and WebScipio are homology-based gene prediction software designed for annotating multigenic families and for transferring annotations from one species to closely related species. The strengths include the power to cope with sequencing-related problems such as sequencing errors and assemblies with short contigs but also the ability to correctly predict genes with unusually long introns and/or rather short exons. WebScipio is connected to diArk, the largest collection of eukaryotic genome assemblies, and thereby offers a very convenient way to correct existing annotations and to extend protein family datasets. WebScipio is also a key resource for researchers interested in mutually exclusive splicing, allowing to search for alternative exons not only in introns but also in up- and downstream regions in case of incompleteness of the search sequence. In this chapter, I describe how to use Scipio and WebScipio keeping a first-time user in mind.

Endogenous Stochastic Decoding of the CUG Codon by Competing Ser- and Leu-tRNAs in Ascoidea asiatica.

Although the "universal" genetic code is now known not to be universal, and stop codons can have multiple meanings, one regularity remains, namely that for a given sense codon there is a unique translation. Examining CUG usage in yeasts that have transferred CUG away from leucine, we here report the first example of dual coding: Ascoidea asiatica stochastically encodes CUG as both serine and leucine in approximately equal proportions. This is deleterious, as evidenced by CUG codons being rare, never at conserved serine or leucine residues, and predominantly in lowly expressed genes. Related yeasts solve the problem by loss of function of one of the two tRNAs. This dual coding is consistent with the tRNA-loss-driven codon reassignment hypothesis, and provides a unique example of a proteome that cannot be deterministically predicted. VIDEO ABSTRACT.

Identifying Sequenced Eukaryotic Genomes and Transcriptomes with diArk.

The diArk Eukaryotic Genome Database is a manually curated and updated repository of available eukaryotic genome and transcriptome assemblies. diArk is a key resource for researchers interested in comparative eukaryotic genomics, and the entry point to browsing sequenced eukaryotes in general and to find the most closely related species to the own organism of interest in particular. The exponentially increasing number of sequenced species demands sophisticated search and data presentation tools. In this chapter we describe how to navigate the diArk database keeping a first-time user in mind.

Waggawagga-CLI: A command-line tool for predicting stable single α-helices (SAH-domains), and the SAH-domain distribution across eukaryotes.

Stable single-alpha helices (SAH-domains) function as rigid connectors and constant force springs between structural domains, and can provide contact surfaces for protein-protein and protein-RNA interactions. SAH-domains mainly consist of charged amino acids and are monomeric and stable in polar solutions, characteristics which distinguish them from coiled-coil domains and intrinsically disordered regions. Although the number of reported SAH-domains is steadily increasing, genome-wide analyses of SAH-domains in eukaryotic genomes are still missing. Here, we present Waggawagga-CLI, a command-line tool for predicting and analysing SAH-domains in protein sequence datasets. Using Waggawagga-CLI we predicted SAH-domains in 24 datasets from eukaryotes across the tree of life. SAH-domains were predicted in 0.5 to 3.5% of the protein-coding content per species. SAH-domains are particularly present in longer proteins supporting their function as structural building block in multi-domain proteins. In human, SAH-domains are mainly used as alternative building blocks not being present in all transcripts of a gene. Gene ontology analysis showed that yeast proteins with SAH-domains are particular enriched in macromolecular complex subunit organization, cellular component biogenesis and RNA metabolic processes, and that they have a strong nuclear and ribonucleoprotein complex localization and function in ribosome and nucleic acid binding. Human proteins with SAH-domains have roles in all types of RNA processing and cytoskeleton organization, and are predicted to function in RNA binding, protein binding involved in cell and cell-cell adhesion, and cytoskeletal protein binding. Waggawagga-CLI allows the user to adjust the stabilizing and destabilizing contribution of amino acid interactions in i,i+3 and i,i+4 spacings, and provides extensive flexibility for user-designed analyses.

The landscape of human mutually exclusive splicing.

Mutually exclusive splicing of exons is a mechanism of functional gene and protein diversification with pivotal roles in organismal development and diseases such as Timothy syndrome, cardiomyopathy and cancer in humans. In order to obtain a first genomewide estimate of the extent and biological role of mutually exclusive splicing in humans, we predicted and subsequently validated mutually exclusive exons (MXEs) using 515 publically available RNA-Seq datasets. Here, we provide evidence for the expression of over 855 MXEs, 42% of which represent novel exons, increasing the annotated human mutually exclusive exome more than fivefold. The data provide strong evidence for the existence of large and multi-cluster MXEs in higher vertebrates and offer new insights into MXE evolution. More than 82% of the MXE clusters are conserved in mammals, and five clusters have homologous clusters in Drosophila Finally, MXEs are significantly enriched in pathogenic mutations and their spatio-temporal expression might predict human disease pathology.

Myosin repertoire expansion coincides with eukaryotic diversification in the Mesoproterozoic era.

BACKGROUND: The last eukaryotic common ancestor already had an amazingly complex cell possessing genomic and cellular features such as spliceosomal introns, mitochondria, cilia-dependent motility, and a cytoskeleton together with several intracellular transport systems. In contrast to the microtubule-based dyneins and kinesins, the actin-filament associated myosins are considerably divergent in extant eukaryotes and a unifying picture of their evolution has not yet emerged. RESULTS: Here, we manually assembled and annotated 7852 myosins from 929 eukaryotes providing an unprecedented dense sequence and taxonomic sampling. For classification we complemented phylogenetic analyses with gene structure comparisons resulting in 79 distinct myosin classes. The intron pattern analysis and the taxonomic distribution of the classes suggest two myosins in the last eukaryotic common ancestor, a class-1 prototype and another myosin, which is most likely the ancestor of all other myosin classes. The sparse distribution of class-2 and class-4 myosins outside their major lineages contradicts their presence in the last eukaryotic common ancestor but instead strongly suggests early eukaryote-eukaryote horizontal gene transfer. CONCLUSIONS: By correlating the evolution of myosin diversity with the history of Earth we found that myosin innovation occurred in independent major "burst" events in the major eukaryotic lineages. Most myosin inventions happened in the Mesoproterozoic era. In the late Neoproterozoic era, a process of extensive independent myosin loss began simultaneously with further eukaryotic diversification. Since the Cambrian explosion, myosin repertoire expansion is driven by lineage- and species-specific gene and genome duplications leading to subfunctionalization and fine-tuning of myosin functions.

Correction: Distribution and evolution of stable single α-helices (SAH domains) in myosin motor proteins.

[This corrects the article DOI: 10.1371/journal.pone.0174639.].

Distribution and evolution of stable single α-helices (SAH domains) in myosin motor proteins.

Stable single-alpha helices (SAHs) are versatile structural elements in many prokaryotic and eukaryotic proteins acting as semi-flexible linkers and constant force springs. This way SAH-domains function as part of the lever of many different myosins. Canonical myosin levers consist of one or several IQ-motifs to which light chains such as calmodulin bind. SAH-domains provide flexibility in length and stiffness to the myosin levers, and may be particularly suited for myosins working in crowded cellular environments. Although the function of the SAH-domains in human class-6 and class-10 myosins has well been characterised, the distribution of the SAH-domain in all myosin subfamilies and across the eukaryotic tree of life remained elusive. Here, we analysed the largest available myosin sequence dataset consisting of 7919 manually annotated myosin sequences from 938 species representing all major eukaryotic branches using the SAH-prediction algorithm of Waggawagga, a recently developed tool for the identification of SAH-domains. With this approach we identified SAH-domains in more than one third of the supposed 79 myosin subfamilies. Depending on the myosin class, the presence of SAH-domains can range from a few to almost all class members indicating complex patterns of independent and taxon-specific SAH-domain gain and loss.

Nuclear codon reassignments in the genomics era and mechanisms behind their evolution.

The canonical genetic code ubiquitously translates nucleotide into peptide sequence with several alterations known in viruses, bacteria, mitochondria, plastids, and single-celled eukaryotes. A new hypothesis to explain genetic code changes, termed tRNA loss driven codon reassignment, has been proposed recently when the polyphyly of the yeast codon reassignment events has been uncovered. According to this hypothesis, the driving force for genetic code changes are tRNA or translation termination factor loss-of-function mutations or loss-of-gene events. The free codon can subsequently be captured by all tRNAs that have an appropriately mutated anticodon and are efficiently charged. Thus, codon capture most likely happens by near-cognate tRNAs and tRNAs whose anticodons are not part of the recognition sites of the respective aminoacyl-tRNA-synthetases. This hypothesis comprehensively explains the CTG codon translation as alanine in Pachysolen yeast together with the long known translation of the same codon as serine in Candida albicans and related species, and can also be applied to most other known reassignments.

How tRNAs dictate nuclear codon reassignments: Only a few can capture non-cognate codons.

mRNA decoding by tRNAs and tRNA charging by aminoacyl-tRNA synthetases are biochemically separated processes that nevertheless in general involve the same nucleotides. The combination of charging and decoding determines the genetic code. Codon reassignment happens when a differently charged tRNA replaces a former cognate tRNA. The recent discovery of the polyphyly of the yeast CUG sense codon reassignment challenged previous mechanistic considerations and led to the proposal of the so-called tRNA loss driven codon reassignment hypothesis. Accordingly, codon capture is caused by loss of a tRNA or by mutations in the translation termination factor, subsequent reduction of the codon frequency through reduced translation fidelity and final appearance of a new cognate tRNA. Critical for codon capture are sequence and structure of the new tRNA, which must be compatible with recognition regions of aminoacyl-tRNA synthetases. The proposed hypothesis applies to all reported nuclear and organellar codon reassignments.

Fine-Tuning Motile Cilia and Flagella: Evolution of the Dynein Motor Proteins from Plants to Humans at High Resolution.

The flagellum is a key innovation linked to eukaryogenesis. It provides motility by regulated cycles of bending and bend propagation, which are thought to be controlled by a complex arrangement of seven distinct dyneins in repeated patterns of outer- (OAD) and inner-arm dynein (IAD) complexes. Electron tomography showed high similarity of this axonemal repeat pattern across ciliates, algae, and animals, but the diversity of dynein sequences across the eukaryotes has not yet comprehensively been resolved and correlated with structural data. To shed light on the evolution of the axoneme I performed an exhaustive analysis of dyneins using the available sequenced genome data. Evidence from motor domain phylogeny allowed expanding the current set of nine dynein subtypes by eight additional isoforms with, however, restricted taxonomic distributions. I confirmed the presence of the nine dyneins in all eukaryotic super-groups indicating their origin predating the last eukaryotic common ancestor. The comparison of the N-terminal tail domains revealed a most likely axonemal dynein origin of the new classes, a group of chimeric dyneins in plants/algae and Stramenopiles, and the unique domain architecture and origin of the outermost OADs present in green algae and ciliates but not animals. The correlation of sequence and structural data suggests the single-headed class-8 and class-9 dyneins to localize to the distal end of the axonemal repeat and the class-7 dyneins filling the region up to the proximal heterodimeric IAD. Tracing dynein gene duplications across the eukaryotes indicated ongoing diversification and fine-tuning of flagellar functions in extant taxa and species.

Coordinated recruitment of Spir actin nucleators and myosin V motors to Rab11 vesicle membranes.

There is growing evidence for a coupling of actin assembly and myosin motor activity in cells. However, mechanisms for recruitment of actin nucleators and motors on specific membrane compartments remain unclear. Here we report how Spir actin nucleators and myosin V motors coordinate their specific membrane recruitment. The myosin V globular tail domain (MyoV-GTD) interacts directly with an evolutionarily conserved Spir sequence motif. We determined crystal structures of MyoVa-GTD bound either to the Spir-2 motif or to Rab11 and show that a Spir-2:MyoVa:Rab11 complex can form. The ternary complex architecture explains how Rab11 vesicles support coordinated F-actin nucleation and myosin force generation for vesicle transport and tethering. New insights are also provided into how myosin activation can be coupled with the generation of actin tracks. Since MyoV binds several Rab GTPases, synchronized nucleator and motor targeting could provide a common mechanism to control force generation and motility in different cellular processes.

Axonopathy in the Central Nervous System Is the Hallmark of Mice with a Novel Intragenic Null Mutation of Dystonin.

Dystonia musculorum is a neurodegenerative disorder caused by a mutation in the dystonin gene. It has been described in mice and humans where it is called hereditary sensory autonomic neuropathy. Mutated mice show severe movement disorders and die at the age of 3-4 weeks. This study describes the discovery and molecular, clinical, as well as pathological characterization of a new spontaneously occurring mutation in the dystonin gene in C57BL/6N mice. The mutation represents a 40-kb intragenic deletion allele of the dystonin gene on chromosome 1 with exactly defined deletion borders. It was demonstrated by Western blot, mass spectrometry, and immunohistology that mice with a homozygous mutation were entirely devoid of the dystonin protein. Pathomorphological lesions were restricted to the brain stem and spinal cord and consisted of swollen, argyrophilic axons and dilated myelin sheaths in the white matter and, less frequently, total chromatolysis of neurons in the gray matter. Axonal damage was detected by amyloid precursor protein and nonphosphorylated neurofilament immunohistology. Axonopathy in the central nervous system (CNS) represents the hallmark of this disease. Mice with the dystonin mutation also showed suppurative inflammation in the respiratory tract, presumably due to brain stem lesion-associated food aspiration, whereas skeletal muscles showed no pathomorphological changes. This study describes a novel mutation in the dystonin gene in mice leading to axonopathy in the CNS. In further studies, this model may provide new insights into the pathogenesis of neurodegenerative diseases and may elucidate the complex interactions of dystonin with various other cellular proteins especially in the CNS.

A novel nuclear genetic code alteration in yeasts and the evolution of codon reassignment in eukaryotes.

The genetic code is the cellular translation table for the conversion of nucleotide sequences into amino acid sequences. Changes to the meaning of sense codons would introduce errors into almost every translated message and are expected to be highly detrimental. However, reassignment of single or multiple codons in mitochondria and nuclear genomes, although extremely rare, demonstrates that the code can evolve. Several models for the mechanism of alteration of nuclear genetic codes have been proposed (including "codon capture," "genome streamlining," and "ambiguous intermediate" theories), but with little resolution. Here, we report a novel sense codon reassignment in Pachysolen tannophilus, a yeast related to the Pichiaceae. By generating proteomics data and using tRNA sequence comparisons, we show that Pachysolen translates CUG codons as alanine and not as the more usual leucine. The Pachysolen tRNACAG is an anticodon-mutated tRNA(Ala) containing all major alanine tRNA recognition sites. The polyphyly of the CUG-decoding tRNAs in yeasts is best explained by a tRNA loss driven codon reassignment mechanism. Loss of the CUG-tRNA in the ancient yeast is followed by gradual decrease of respective codons and subsequent codon capture by tRNAs whose anticodon is not part of the aminoacyl-tRNA synthetase recognition region. Our hypothesis applies to all nuclear genetic code alterations and provides several testable predictions. We anticipate more codon reassignments to be uncovered in existing and upcoming genome projects.