Origin and Early Diversification of the Papain Family of Cysteine Peptidases.

Peptidases of the papain family play a key role in protein degradation, regulated proteolysis, and the host-pathogen arms race. Although the papain family has been the subject of many studies, knowledge about its diversity, origin, and evolution in Eukaryota, Bacteria, and Archaea is limited; thus, we aimed to address these long-standing knowledge gaps. We traced the origin and expansion of the papain family with a phylogenomic analysis, using sequence data from numerous prokaryotic and eukaryotic proteomes, transcriptomes, and genomes. We identified the full complement of the papain family in all prokaryotic and eukaryotic lineages. Analysis of the papain family provided strong evidence for its early diversification in the ancestor of eukaryotes. We found that the papain family has undergone complex and dynamic evolution through numerous gene duplications, which produced eight eukaryotic ancestral paralogous C1A lineages during eukaryogenesis. Different evolutionary forces operated on C1A peptidases, including gene duplication, horizontal gene transfer, and gene loss. This study challenges the current understanding of the origin and evolution of the papain family and provides valuable insights into their early diversification. The findings of this comprehensive study provide guidelines for future structural and functional studies of the papain family.

Tum1 is involved in the metabolism of sterol esters in Saccharomyces cerevisiae.

BACKGROUND: The only hitherto known biological role of yeast Saccharomyces cerevisiae Tum1 protein is in the tRNA thiolation pathway. The mammalian homologue of the yeast TUM1 gene, the thiosulfate sulfurtransferase (a.k.a. rhodanese) Tst, has been proposed as an obesity-resistance and antidiabetic gene. To assess the role of Tum1 in cell metabolism and the putative functional connection between lipid metabolism and tRNA modification, we analysed evolutionary conservation of the rhodanese protein superfamily, investigated the role of Tum1 in lipid metabolism, and examined the phenotype of yeast strains expressing the mouse homologue of Tum1, TST. RESULTS: We analysed evolutionary relationships in the rhodanese superfamily and established that its members are widespread in bacteria, archaea and in all major eukaryotic groups. We found that the amount of sterol esters was significantly higher in the deletion strain tum1Δ than in the wild-type strain. Expression of the mouse TST protein in the deletion strain did not rescue this phenotype. Moreover, although Tum1 deficiency in the thiolation pathway was complemented by re-introducing TUM1, it was not complemented by the introduction of the mouse homologue Tst. We further showed that the tRNA thiolation pathway is not involved in the regulation of sterol ester content in S. cerevisiae, as overexpression of the tEUUC, tKUUU and tQUUG tRNAs did not rescue the lipid phenotype in the tum1Δ deletion strain, and, additionally, deletion of the key gene for the tRNA thiolation pathway, UBA4, did not affect sterol ester content. CONCLUSIONS: The rhodanese superfamily of proteins is widespread in all organisms, and yeast TUM1 is a bona fide orthologue of mammalian Tst thiosulfate sulfurtransferase gene. However, the mouse TST protein cannot functionally replace yeast Tum1 protein, neither in its lipid metabolism-related function, nor in the tRNA thiolation pathway. We show here that Tum1 protein is involved in lipid metabolism by decreasing the sterol ester content in yeast cells, and that this function of Tum1 is not exerted through the tRNA thiolation pathway, but through another, currently unknown pathway.

Genesis and regulatory wiring of retroelement-derived domesticated genes: a phylogenomic perspective.

Molecular domestications of transposable elements have occurred repeatedly during the evolution of eukaryotes. Vertebrates, especially mammals, possess numerous single copy domesticated genes (DGs) that have originated from the intronless multicopy transposable elements. However, the origin and evolution of the retroelement-derived DGs (RDDGs) that originated from Metaviridae has been only partially elucidated, due to absence of genome data or to limited analysis of a single family of DGs. We traced the genesis and regulatory wiring of the Metaviridae-derived DGs through phylogenomic analysis, using whole-genome information from more than 90 chordate genomes. Phylogenomic analysis of these DGs in chordate genomes provided direct evidence that major diversification has occurred in the ancestor of placental mammals. Mammalian RDDGs have been shown to originate in several steps by independent domestication events and to diversify later by gene duplications. Analysis of syntenic loci has shown that diverse RDDGs and their chromosomal positions were fully established in the ancestor of placental mammals. By analysis of active Metaviridae lineages in amniotes, we have demonstrated that RDDGs originated from retroelement remains. The chromosomal gene movements of RDDGs were highly dynamic only in the ancestor of placental mammals. During the domestication process, de novo acquisition of regulatory regions is shown to be a prerequisite for the survival of the DGs. The origin and evolution of de novo acquired promoters and untranslated regions in diverse mammalian RDDGs have been explained by comparative analysis of orthologous gene loci. The origin of placental mammal-specific innovations and adaptations, such as placenta and newly evolved brain functions, was most probably connected to the regulatory wiring of DGs and their rapid fixation in the ancestor of placental mammals.

Serine pseudoproteases in physiology and disease.

Serine proteases (SPs) constitute a very important family of enzymes, both physiologically and pathologically. The effects produced by these proteins have been explained by their proteolytic activity. However, the discovery of pharmacologically active SP molecules that show no enzymatic activity, as the so-called pseudo SPs or SP homologs (SPHs), has exposed a profoundly neglected possibility of nonenzymatic functions of these SP molecules. In this review, the most thoroughly described SPHs are presented. The main physiological domains in which SPHs operate appear to be in reproduction, embryonic development, immune response, host defense, and hemostasis. Hitherto unexplained actions of SPs should therefore be considered also as the result of the ligand-like attributes of SPs. The gain of a novel function by an SPH is a consequence of specific amino acid replacements that have resulted in a novel interaction interface or a 'catalytic trap'. Unraveling the SP/SPH interactome will provide a description of previously unknown physiological functions of SPs/SPHs, aiding the creation of innovative medical approaches.

Conus consors snail venom proteomics proposes functions, pathways, and novel families involved in its venomic system.

For some decades, cone snail venoms have been providing peptides, generally termed conopeptides, that exhibit a large diversity of pharmacological properties. However, little attention has been devoted to the high molecular mass (HMM) proteins in venoms of mollusks. In order to shed more light on cone snail venom HMM components, the proteins of dissected and injected venom of a fish-hunting cone snail, Conus consors, were extensively assessed. HMM venom proteins were separated by two-dimensional polyacrylamide gel electrophoresis and analyzed by mass spectrometry (MS). The MS data were interpreted using UniProt database, EST libraries from C. consors venom duct and salivary gland, and their genomic information. Numerous protein families were discovered in the lumen of the venom duct and assigned a biological function, thus pointing to their potential role in venom production and maturation. Interestingly, the study also revealed original proteins defining new families of unknown function. Only two groups of HMM proteins passing the venom selection process, echotoxins and hyaluronidases, were clearly present in the injected venom. They are suggested to contribute to the envenomation process. This newly devised integrated HMM proteomic analysis is a big step toward identification of the protein arsenal used in a cone snail venom apparatus for venom production, maturation, and function.

Toward the massive genome of Proteus anguinus-illuminating longevity, regeneration, convergent evolution, and metabolic disorders.

Deciphering the genetic code of organisms with unusual phenotypes can help answer fundamental biological questions and provide insight into mechanisms relevant to human biomedical research. The cave salamander Proteus anguinus (Urodela: Proteidae), also known as the olm, is an example of a species with unique morphological and physiological adaptations to its subterranean environment, including regenerative abilities, resistance to prolonged starvation, and a life span of more than 100 years. However, the structure and sequence of the olm genome is still largely unknown owing to its enormous size, estimated at nearly 50 gigabases. An international Proteus Genome Research Consortium has been formed to decipher the olm genome. This perspective provides the scientific and biomedical rationale for exploring the olm genome and outlines potential outcomes, challenges, and methodological approaches required to analyze and annotate the genome of this unique amphibian.

Evolution of phospholipase A2 toxins in venomous animals.

Franc Gubensek devoted much of his research career to the phospholipases A2 (PLA2), which are the major pharmacologically active components of snake venoms. Our long collaboration started with an analysis of Vipera ammodytes ammodytoxin and ammodytin cDNAs and genes. These PLA2 genes provided us with an entry into the exciting area of molecular evolution. We studied the structures of the PLA2 genes, the evolution of multigene families encoding PLA2 toxins, and the horizontal transfer of unusual retroelements that we found in these genes. In the last decade a number of novel features have emerged concerning the evolution of PLA2s in venomous animals. The large amount of recently accumulated data has provided a timely opportunity to review current understanding of the evolution of PLA2 toxins in venomous animals.

The Evolution of αD-Conopeptides Targeting Neuronal Nicotinic Acetylcholine Receptors.

Venoms of the marine cone snails (Conus spp.) consist of numerous proteins and peptides showing a wide variety of biological activities such as on ion-channels and receptors. Peptides acting on neuronal nicotinic acetylcholine receptors belong to several peptide superfamilies including the recently described αD-conopeptides which are homodimers of identical peptides with 47-49 amino acids. Among the venom glands of 27 Conus species analyzed by cDNA cloning, precursors of αD-conopeptides were identified in four species only: C. betulinus, C. capitaneus, C. mustelinus, and C. vexillum. Phylogenetic analysis of the relationships among the αD-conopeptides revealed that they belong to clades, which are characterized by an AVV- and EMM-motif in the signal peptide sequence.

Analysis of the RNA virome of basal hexapods.

The diversity and evolution of RNA viruses has been well studied in arthropods and especially in insects. However, the diversity of RNA viruses in the basal hexapods has not been analysed yet. To better understand their diversity, evolutionary histories and genome organizations, we searched for RNA viruses in transcriptome and genome databases of basal hexapods. We discovered  40 novel RNA viruses, some of which are also present as endogenous viral elements derived from RNA viruses. Here, we demonstrated that basal hexapods host 14 RNA viral clades that have been recently identified in invertebrates. The following RNA viral clades are associated with basal hexapods: Reo, Partiti-Picobirna, Toti-Chryso, Mono-Chu, Bunya-Arena, Orthomyxo, Qinvirus, Picorna-Calici, Hepe-Virga, Narna-Levi, Tombus-Noda, Luteo-Sobemo, Permutotetra and Flavi. We have found representatives of the nine RNA viral clades that are present as endogenous genomic copies in the genomes of Machilis (Monocondylia) and Catajapyx (Diplura). Our study provided a first insight into the diversity of RNA viruses in basal hexapods and demonstrated that the basal hexapods possess quite high diversity of RNA viral clades.

Evolutionary History of Alpha Satellite DNA Repeats Dispersed within Human Genome Euchromatin.

Major human alpha satellite DNA repeats are preferentially assembled within (peri)centromeric regions but are also dispersed within euchromatin in the form of clustered or short single repeat arrays. To study the evolutionary history of single euchromatic human alpha satellite repeats (ARs), we analyzed their orthologous loci across the primate genomes. The continuous insertion of euchromatic ARs throughout the evolutionary history of primates starting with the ancestors of Simiformes (45-60 Ma) and continuing up to the ancestors of Homo is revealed. Once inserted, the euchromatic ARs were stably transmitted to the descendant species, some exhibiting copy number variation, whereas their sequence divergence followed the species phylogeny. Many euchromatic ARs have sequence characteristics of (peri)centromeric alpha repeats suggesting heterochromatin as a source of dispersed euchromatic ARs. The majority of euchromatic ARs are inserted in the vicinity of other repetitive elements such as L1, Alu, and ERV or are embedded within them. Irrespective of the insertion context, each AR insertion seems to be unique and once inserted, ARs do not seem to be subsequently spread to new genomic locations. In spite of association with (retro)transposable elements, there is no indication that such elements play a role in ARs proliferation. The presence of short duplications at most of ARs insertion sites suggests site-directed recombination between homologous motifs in ARs and in the target genomic sequence, probably mediated by extrachromosomal circular DNA, as a mechanism of spreading within euchromatin.

Widespread Evolution of Molecular Resistance to Snake Venom α-Neurotoxins in Vertebrates.

Venomous snakes are important subjects of study in evolution, ecology, and biomedicine. Many venomous snakes have alpha-neurotoxins (α-neurotoxins) in their venom. These toxins bind the alpha-1 nicotinic acetylcholine receptor (nAChR) at the neuromuscular junction, causing paralysis and asphyxia. Several venomous snakes and their predators have evolved resistance to α-neurotoxins. The resistance is conferred by steric hindrance from N-glycosylated asparagines at amino acids 187 or 189, by an arginine at position 187 that has been hypothesized to either electrostatically repulse positively charged neurotoxins or sterically interfere with α-neurotoxin binding, or proline replacements at positions 194 or 197 of the nAChR ligand-binding domain to inhibit α-neurotoxin binding through structural changes in the receptor. Here, we analyzed this domain in 148 vertebrate species, and assessed its amino acid sequences for resistance-associated mutations. Of these sequences, 89 were sequenced de novo. We find widespread convergent evolution of the N-glycosylation form of resistance in several taxa including venomous snakes and their lizard prey, but not in the snake-eating birds studied. We also document new lineages with the arginine form of inhibition. Using an in vivo assay in four species, we provide further evidence that N-glycosylation mutations reduce the toxicity of cobra venom. The nAChR is of crucial importance for normal neuromuscular function and is highly conserved throughout the vertebrates as a result. Our research shows that the evolution of α-neurotoxins in snakes may well have prompted arms races and mutations to this ancient receptor across a wide range of sympatric vertebrates. These findings underscore the inter-connectedness of the biosphere and the ripple effects that one adaption can have across global ecosystems.

Phylogenomic analysis of the cystatin superfamily in eukaryotes and prokaryotes.

BACKGROUND: The cystatin superfamily comprises cysteine protease inhibitors that play key regulatory roles in protein degradation processes. Although they have been the subject of many studies, little is known about their genesis, evolution and functional diversification. Our aim has been to obtain a comprehensive insight into their origin, distribution, diversity, evolution and classification in Eukaryota, Bacteria and Archaea. RESULTS: We have identified in silico the full complement of the cystatin superfamily in more than 2100 prokaryotic and eukaryotic genomes. The analysis of numerous eukaryotic genomes has provided strong evidence for the emergence of this superfamily in the ancestor of eukaryotes. The progenitor of this superfamily was most probably intracellular and lacked a signal peptide and disulfide bridges, much like the extant Giardia cystatin. A primordial gene duplication produced two ancestral eukaryotic lineages, cystatins and stefins. While stefins remain encoded by a single or a small number of genes throughout the eukaryotes, the cystatins have undergone a more complex and dynamic evolution through numerous gene and domain duplications. In the cystatin superfamily we discovered twenty vertebrate-specific and three angiosperm-specific orthologous families, indicating that functional diversification has occurred only in multicellular eukaryotes. In vertebrate orthologous families, the prevailing trends were loss of the ancestral inhibitory activity and acquisition of novel functions in innate immunity. Bacterial cystatins and stefins may be emergency inhibitors that enable survival of bacteria in the host, defending them from the host's proteolytic activity. CONCLUSION: This study challenges the current view on the classification, origin and evolution of the cystatin superfamily and provides valuable insights into their functional diversification. The findings of this comprehensive study provide guides for future structural and evolutionary studies of the cystatin superfamily as well as of other protease inhibitors and proteases.

Characterization of a lipid droplet protein from Yarrowia lipolytica that is required for its oleaginous phenotype.

Oleaginous microorganisms are characterized by their ability to store high amounts of triacylglycerol (TAG) in intracellular lipid droplets (LDs). In this work, we characterized a protein of the oleaginous yeast Yarrowia lipolytica that is associated with LD and plays a role in the regulation of TAG storage. This protein is required for the oleaginous phenotype of Y. lipolytica because deletion of the coding gene results in a strongly reduced TAG content of the mutant. Therefore, we named it Oleaginicity Inducing LD protein, Oil1. Furthermore, a mutant overexpressing OIL1 accumulates more TAG than the wild type and is delayed in TAG lipolysis when this process is stimulated. We found that Oil1p plays a role in protecting the TAG content of the LD from degradation through lipases under conditions where the cell aims at building up its TAG reserves. Heterologous expression studies showed that Oil1p rescued the phenotype of a Saccharomyces cerevisiae mutant deleted for the perilipin-like protein Pln1p and that its expression in COS-7 cells resulted in increased TAG accumulation, similar to the phenotype of a perilipin 1 expressing control strain. Despite this phenotypical parallels to mammalian perilipins, Oil1p is not a member of this protein family and its activity does not depend on phosphorylation. Rather, our results suggest that ubiquitination might contribute to the function of Oil1p in Y. lipolytica and that a different mechanism evolved in this species to regulate TAG homeostasis.

Strong and widespread action of site-specific positive selection in the snake venom Kunitz/BPTI protein family.

S1 family of serine peptidases is the largest family of peptidases. They are specifically inhibited by the Kunitz/BPTI inhibitors. Kunitz domain is characterized by the compact 3D structure with the most important inhibitory loops for the inhibition of S1 peptidases. In the present study we analysed the action of site-specific positive selection and its impact on the structurally and functionally important parts of the snake venom Kunitz/BPTI family of proteins. By using numerous models we demonstrated the presence of large numbers of site-specific positively selected sites that can reach between 30-50% of the Kunitz domain. The mapping of the positively selected sites on the 3D model of Kunitz/BPTI inhibitors has shown that these sites are located in the inhibitory loops 1 and 2, but also in the Kunitz scaffold. Amino acid replacements have been found exclusively on the surface, and the vast majority of replacements are causing the change of the charge. The consequence of these replacements is the change in the electrostatic potential on the surface of the Kunitz/BPTI proteins that may play an important role in the precise targeting of these inhibitors into the active site of S1 family of serine peptidases.

A genomic perspective on the chromodomain-containing retrotransposons: Chromoviruses.

Chromoviruses, chromodomain-containing retrotransposons, are the only Metaviridae (Ty3/gypsy group of retrotransposons) clade with a Eukaryota-wide distribution. They have a common evolutionary origin and are the most prolific and diverse Metaviridae clade. The fusion of a retrotransposon and a chromodomain, was most probably responsible for their extreme evolutionary success in Eukaryota. Analysis of the massive amount of genome sequence data for different eukaryotic lineages has provided an in depth insight into the diversity, evolution, neofunctionalization, high rate of genomic turnover and origin of chromoviruses in Eukaryota. This review attempts to summarise the unique aspects of chromoviruses from a genomic perspective.

Adaptive evolution in the snake venom Kunitz/BPTI protein family.

Snake venoms are rich sources of serine proteinase inhibitors that are members of the Kunitz/BPTI (bovine pancreatic trypsin inhibitor) family. However, only a few of their gene sequences have been determined from snakes. We therefore cloned the cDNAs for the trypsin and chymotrypsin inhibitors from a Vipera ammodytes venom gland cDNA library. Phylogenetic analysis of these and other snake Kunitz/BPTI homologs shows the presence of three clusters, where sequences cluster by functional role. Analysis of the nucleotide sequences from the snake Kunitz/BPTI family shows that positive Darwinian selection was operating on the highly conserved BPTI fold, indicating that this family evolved by gene duplication and rapid diversification.

What Can Domesticated Genes Tell Us about the Intron Gain in Mammals?

Domesticated genes, originating from retroelements or from DNA-transposons, constitute an ideal system for testing the hypothesis on the absence of intron gain in mammals. Since single-copy domesticated genes originated from the intronless multicopy transposable elements, the ancestral intron state for domesticated genes is zero. A phylogenomic approach has been used to analyse all domesticated genes in mammals and chordates that originated from the coding parts of transposable elements. A significant amount of intron gain was found only in domesticated genes of placental mammals, where more than 70 cases were identified. De novo gained introns show clear positional bias, since they are distributed mainly in 5' UTR and coding regions, while 3' UTR introns are very rare. In the coding regions of some domesticated genes up to 8 de novo gained introns have been found. Surprisingly, the majority of intron gains have occurred in the ancestor of placental mammals. Domesticated genes could constitute an excellent system on which to analyse the mechanisms of intron gain. This paper summarizes the current understanding of intron gain in mammals.

Variability of tetrodotoxin and of its analogues in the red-spotted newt, Notophthalmus viridescens (Amphibia: Urodela: Salamandridae).

Efts and adult specimens (n = 142) of the red-spotted newt Notophthalmus viridescens from various locations in Canada and USA were analyzed for the presence of tetrodotoxin (TTX) and of its analogues 6-epitetrodotoxin and 11-oxotetrodotoxin. Considerable individual variations in toxin levels were found within and among populations from New Hampshire, New York, Pennsylvania, and Virginia ranging from non-detectable to 69 µg TTX per g newt. TTX and its analogues were absent in efts and adults from various locations in the Canadian province Nova Scotia, the northernmost distribution of the newt, and in adults from Florida. Newts kept in captivity for several years and reared on toxin-free diet lost their toxicity. Bayesian and maximum likelihood phylogenetic analysis of specimens from the various populations using three phylogenetic markers (COI, ND2 and 16S RNA) revealed that populations from the northern states of the USA and Canada are genetically homogenous, whereas the newts from Florida exhibited a much higher level of genetic divergence. An exogenous source of TTX in the newts either via the food chain or by synthesis of symbiotic bacteria is suggested to explain the high variability and lack of TTX in certain populations.

Extensive intron gain in the ancestor of placental mammals.

BACKGROUND: Genome-wide studies of intron dynamics in mammalian orthologous genes have found convincing evidence for loss of introns but very little for intron turnover. Similarly, large-scale analysis of intron dynamics in a few vertebrate genomes has identified only intron losses and no gains, indicating that intron gain is an extremely rare event in vertebrate evolution. These studies suggest that the intron-rich genomes of vertebrates do not allow intron gain. The aim of this study was to search for evidence of de novo intron gain in domesticated genes from an analysis of their exon/intron structures. RESULTS: A phylogenomic approach has been used to analyse all domesticated genes in mammals and chordates that originated from the coding parts of transposable elements. Gain of introns in domesticated genes has been reconstructed on well established mammalian, vertebrate and chordate phylogenies, and examined as to where and when the gain events occurred. The locations, sizes and amounts of de novo introns gained in the domesticated genes during the evolution of mammals and chordates has been analyzed. A significant amount of intron gain was found only in domesticated genes of placental mammals, where more than 70 cases were identified. De novo gained introns show clear positional bias, since they are distributed mainly in 5' UTR and coding regions, while 3' UTR introns are very rare. In the coding regions of some domesticated genes up to 8 de novo gained introns have been found. Intron densities in Eutheria-specific domesticated genes and in older domesticated genes that originated early in vertebrates are lower than those for normal mammalian and vertebrate genes. Surprisingly, the majority of intron gains have occurred in the ancestor of placentals. CONCLUSIONS: This study provides the first evidence for numerous intron gains in the ancestor of placental mammals and demonstrates that adequate taxon sampling is crucial for reconstructing intron evolution. The findings of this comprehensive study slightly challenge the current view on the evolutionary stasis in intron dynamics during the last 100 - 200 My. Domesticated genes could constitute an excellent system on which to analyse the mechanisms of intron gain in placental mammals.

Venomic study on cone snails (Conus spp.) from South Africa.

From six Conus species (Conus coronatus, Conus lividus, Conus mozambicus f. lautus, Conus pictus, Conus sazanka, Conus tinianus) collected off the eastern coast of South Africa the venoms were analyzed using MALDI-TOF mass spectrometry. Between 56 and 151 molecular masses most in a range of 1000 to 2500 Da, were identified. Among the six venoms, between 0 and 27% (C. coronatus versus C. sazanka) of the peptide masses were found to be similar. In a study on venoms from 6 Conus species collected in the Philippines, the percentage of identical masses was between none and 9% only. The venoms from the South African Conus species antagonized the rat neuronal nicotinic acetylcholine receptors (nAChRs) α3ß2, α4ß2, and α7, except for C. coronatus venom that blocked the α4ß2 and α7 nAChRs only. HPLC-fractionation of C. tinianus venom led to the isolation of a peptide that is active on all three receptor subtypes. It consists of 16 amino acid residues cross-linked by two disulfide bridges as revealed by de novo sequencing using tandem mass spectrometry: GGCCSHPACQNNPDYC. Posttranslational modifications include C-terminal amidation and tyrosine sulfation. The new peptide is a member of the α-conotoxin family that are competitive antagonists of nAChRs. Phylogenetic analysis of the 16S RNA from numerous Conus species has clarified the evolutionary position of endemic South African Conus species and provided the first evidence for their close genetic relationships.