A serpin in the cellulosome of the anaerobic fungus Piromyces sp. strain E2

A gene encoding a novel component of the cellulolytic complex (cellulosome) of the anaerobic fungus Piromyces sp. strain E2 was identified. The encoded 538 amino acid protein, named celpin, consists of a signal peptide, a positively charged domain of unknown function followed by two fungal dockerins, typical for components of the extracellular fungal cellulosome. The C-terminal end consists of a 380 amino acid serine proteinase inhibitor (or serpin) domain homologue, sharing 30 % identity and 50 % similarity to vertebrate and bacterial serpins. Detailed protein sequence analysis of the serpin domain revealed that it contained all features of a functional serpin. It possesses the conserved amino acids present in more than 70 % of known serpins, and it contained the consensus of inhibiting serpins. Because of the confined space of the fungal cellulosome inside plant tissue and the auto-proteolysis of plant material in the rumen, the fungal serpin is presumably involved in protection of the cellulosome against plant proteinases. The celpin protein of Piromyces sp. strain E2 is the first non-structural, non-hydrolytic fungal cellulosome component. Furthermore, the celpin protein of Piromyces sp. strain E2 is the first representative of a serine proteinase inhibitor of the fungal kingdom.     

Networks of coevolving sites in structural and functional domains of serpin proteins

Amino acids do not occur randomly in proteins; rather, their occurrence at any given site is strongly influenced by the amino acid composition at other sites, the structural and functional aspects of the region of the protein in which they occur, and the evolutionary history of the protein. The goal of our research study is to identify networks of coevolving sites within the serpin proteins (serine protease inhibitors) and classify them as being caused by structural-functional constraints or by evolutionary history. To address this, a matrix of pairwise normalized mutual information (NMI) values was computed among amino acid sites for the serpin proteins. The NMI matrix was partitioned into orthogonal patterns of amino acid variability by factor analysis. Each common factor pattern was interpreted as having phylogenetic and/or structural-functional explanations. In addition, we used a bootstrap factor analysis technique to limit the effects of phylogenetic history on our factor patterns. Our results show an extensive network of correlations among amino acid sites in key functional regions (reactive center loop, shutter, and breach). Additionally, we have discovered long-range coevolution for packed amino acids within the serpin protein core. Lastly, we have discovered a group of serpin sites which coevolve in the hydrophobic core region (s5B and s4B) and appear to represent sites important for formation of the "native" instead of the "latent" serpin structure. This research provides a better understanding on how protein structure evolves; in particular, it elucidates the selective forces creating coevolution among protein sites.     

Pattern of divergence of amino acid sequences encoded by paralogous genes in human and pufferfish

We used phylogenetic analyses of protein families containing two or more pairs of orthologues in the genomes of human and pufferfish (Takifugu rubripes) to test the hypothesis that these sequences show a strong signal of polyploidization events hypothesized to have occurred early in vertebrate history. In order to test for evidence of two distinct rounds of polyploidization (the 2R hypothesis), we compared the pattern of amino acid sequence divergence of proteins encoded by genes duplicated just prior to the most recent common ancestor of human and pufferfish with that of proteins encoded genes duplicated earlier. These sequence divergences were statistically indistinguishable, contrary to the prediction of the 2R hypothesis. The variance of amino acid sequence divergences between paralogues was significantly greater than expected from that of orthologues in the same families. Estimation of gene duplication times assuming a molecular clock provided earlier estimates than expected, suggesting that it may not be appropriate to time the duplication of paralogues using rate estimates derived from orthologous comparisons. Overall, the results indicate that amino acid sequences do not provide a strong signal supporting the hypothesis that gene duplications early in vertebrate history occurred by polyploidization. On the other hand, the data are easily explained under an alternative model that gene duplications occurred at different times in different vertebrate gene families.     

Sequence, genomic structure and tissue expression of carp (Cyprinus carpio L.) vertebrate ancient (VA) opsin

We report the isolation and characterisation of a novel opsin cDNA from the retina and pineal of the common carp (Cyprinus carpio L.). When a comparison of the amino acid sequences of salmon vertebrate ancient opsin (sVA) and the novel carp opsin are made, and the carboxyl terminus is omitted, the level of identity between these two opsins is 81% and represents the second example of the VA opsin family. We have therefore termed this C. carpio opsin as carp VA opsin (cVA opsin). We show that members of the VA opsin family may exist in two variants or isoforms based upon the length of the carboxyl terminus and propose that the mechanism of production of the short VA opsin isoform is alternative splicing of intron 4 of the VA opsin gene. The VA opsin gene consists of five exons, with intron 2 significantly shifted in a 3' direction relative to the corresponding intron in rod and cone opsins. The position (or lack) of intron 2 appears to be a diagnostic feature which separates the image forming rod and cone opsin families from the more recently discovered non-visual opsin families (pin-opsins (P), vertebrate ancient (VA), parapinopsin (PP)). Finally, we suggest that lamprey P opsin should be reassigned to the VA opsin family based upon its level of amino acid identity, genomic structure with respect to the position of intron 2 and nucleotide phylogeny.     

Toxin-resistant sodium channels: parallel adaptive evolution across a complete gene family

Approximately 75% of vertebrate proteins belong to protein families encoded by multiple evolutionarily related genes, a pattern that emerged as a result of gene and genome duplications over the course of vertebrate evolution. In families of genes with similar or related functions, adaptation to a strong selective agent should involve multiple adaptive changes across the entire gene family. However, we know of no evolutionary studies that have explicitly addressed this point. Here, we show how 4 taxonomically diverse species of pufferfishes (Tetraodontidae) each evolved resistance to the guanidinium toxins tetrodotoxin (TTX) and saxitoxin (STX) via parallel amino acid replacements across all 8 sodium channels present in teleost fish genomes. This resulted in diverse suites of coexisting sodium channel types that all confer varying degrees of toxin resistance, yet show remarkable convergence among genes and phylogenetically diverse species. Using site-directed mutagenesis and expression of a vertebrate sodium channel, we also demonstrate that resistance to TTX/STX is enhanced up to 15-fold by single, frequently observed replacements at 2 sites that have not previously been implicated in toxin binding but show similar or identical replacements in pufferfishes and in distantly related vertebrate and nonvertebrate animals. This study presents an example of natural selection acting upon a complete gene family, repeatedly arriving at a diverse but limited number of adaptive changes within the same genome. To be maximally informative, we suggest that future studies of molecular adaptation should consider all functionally similar paralogs of the affected gene family.     

Phosphoglucose isomerases of hagfish, zebrafish, gray mullet, toad, and snake, with reference to the evolution of the genes in vertebrates

Phosphoglucose isomerase (PGI) is a protein with multiple functions. To infer its structure changes and evolution in vertebrates, we cloned cDNAs encoding PGI genes from hagfish (Paramyxine yangi), gray mullet (Mugil cephalus), zebrafish (Danio rerio), toad (Bufo melanosticus), and snake (Boiga kraepelini). Only one PGI gene was cloned in each of hagfish, toad, and snake, but two PGI genes were found in zebrafish and gray mullet, respectively. The PGI of hagfish encodes 554 amino acids, in contrast to the PGIs of bonyfishes, toad, and snake which encode 553 amino acids and the PGIs of mammals which encode 558 amino acids. Among 558 aligned amino acid sites, there are 314 sites (56.27%) totally conserved. To see if diversifying selection acts on PGI amino acids of vertebrates, we calculated the pairwise ratio of nonsynonymous versus synonymous substitution per site (Ka/Ks) and the ratio of radical amino acid changes versus conservative amino acid changes per sites (dR/dC) between PGI sequences. The average pairwise ratio between nonsynonymous substitutions per nucleotide (Ka) and synonymous substitutions per nucleotide (Ks) among vertebrate PGI sequences equals 0.047 +/- 0.019. The average pairwise ratio between radical amino acid changes and conservative amino acid changes (dR/dC) among the vertebrate PGIs equal 0.938 +/- 0.158 for charge changes, 0.558 +/- 0.085 for polarity changes, and 0.465 +/- 0.0714 when both polarity and volume are considered. There is no amino acid within the vertebrate PGIs under diversifying selection as analyzed by the method of Yang et al. (2000b). The results suggest that the present vertebrate PGIs are at evolutionary stasis and are being subjected to intense purifying selection. The purifying selection is to maintain polarity and volume of the protein but not the charge groups of amino acids. Phylogenetic analysis reveals that vertebrate PGIs can be classified into three major groups: the mammalian, amphibian-reptilian, and teleostean PGIs. The gene tree suggests that the gene duplication event of PGI in bonyfishes occurred before diversification of Acanthopterygii but after the split of bonyfishes and tetrapods. The evolution of multiple functions of PGI is discussed.     

Comparative studies of vertebrate lipoprotein lipase: a key enzyme of very low density lipoprotein metabolism

Lipoprotein lipase (LIPL or LPL; E.C.3.1.1.34) serves a dual function as a triglyceride lipase of circulating chylomicrons and very-low-density lipoproteins (VLDL) and facilitates receptor-mediated lipoprotein uptake into heart, muscle and adipose tissue. Comparative LPL amino acid sequences and protein structures and LPL gene locations were examined using data from several vertebrate genome projects. Mammalian LPL genes usually contained 9 coding exons on the positive strand. Vertebrate LPL sequences shared 58-99% identity as compared with 33-49% sequence identities with other vascular triglyceride lipases, hepatic lipase (HL) and endothelial lipase (EL). Two human LPL N-glycosylation sites were conserved among seven predicted sites for the vertebrate LPL sequences examined. Sequence alignments, key amino acid residues and conserved predicted secondary and tertiary structures were also studied. A CpG island was identified within the 5'-untranslated region of the human LPL gene which may contribute to the higher than average (×4.5 times) level of expression reported. Phylogenetic analyses examined the relationships and potential evolutionary origins of vertebrate lipase genes, LPL, LIPG (encoding EL) and LIPC (encoding HL) which suggested that these have been derived from gene duplication events of an ancestral neutral lipase gene, prior to the appearance of fish during vertebrate evolution. Comparative divergence rates for these vertebrate sequences indicated that LPL is evolving more slowly (2-3 times) than for LIPC and LIPG genes and proteins.     

Complex history of a chromosomal paralogy region: insights from amphioxus aromatic amino acid hydroxylase genes and insulin-related genes

Aromatic amino acid hydroxylase (AAAH) genes and insulin-like genes form part of an extensive paralogy region shared by human chromosomes 11 and 12, thought to have arisen by tetraploidy in early vertebrate evolution. Cloning of a complementary DNA (cDNA) for an amphioxus (Branchiostoma floridae) hydroxylase gene (AmphiPAH) allowed us to investigate the ancestry of the human chromosome 11/12 paralogy region. Molecular phylogenetic evidence reveals that AmphiPAH is orthologous to vertebrate phenylalanine (PAH) genes; the implication is that all three vertebrate AAAH genes arose early in metazoan evolution, predating vertebrates. In contrast, our phylogenetic analysis of amphioxus and vertebrate insulin-related gene sequences is consistent with duplication of these genes during early chordate ancestry. The conclusion is that two tightly linked gene families on human chromosomes 11 and 12 were not duplicated coincidentally. We rationalize this paradox by invoking gene loss in the AAAH gene family and conclude that paralogous genes shared by paralogous chromosomes need not have identical evolutionary histories.     

Revising angiotensinogen from phylogenetic and genetic variants perspectives

Angiotensinogen (AGT) belongs to the serpin superfamily. It acts as the unique substrate of all angiotensin peptides, which generates a spectrum of angiotensin peptides in the renin-angiotensin system and regulates hypertension. This serpin belongs to the multiple member group V2 of the intron encoded vertebrate serpin classification. Despite huge advancements in the understanding of angiotensinogen based on biochemical properties and its roles in the RAS, phylogenetic history of AGT remains forgotten. To date, there is no comprehensive study illustrating the phylogenetic history of AGT. Herein, we investigated phylogenetic traits of AGT gene across vertebrates. Gene structures of AGT gene from selected ray-finned fishes varied in exon I and II with insertions of two novel introns in the core domain for ray-finned fishes at the position 77c and 233c. We that found AGT loci is conserved from lampreys to human and estimated to be older than 500 MY. By comparing AGT protein in 57 vertebrate genomes, we illustrated that the reactive center loop (RCL) of AGT protein became from inhibitory (in lampreys, GTEAKAETVVGIMPI†SMPPT) to non-inhibitory (in human, EREPTESTQQLNKPE†VLEVT) during period of 500 MY. We identified 690 AGT variants by analysis of 1092 human genomes with top three variation classes belongs to SNPs (89.7%), somatic SNVs (5.2%) and deletion (2.9%). There are 32 key residues out of 121 missense variants, which are deleterious for AGT protein, computed by combination of SIFT and PolyPhen V2 methods. These results may have clinical implications for understanding hypertension.     

Identification and cloning of caprine uterine serpin

The uterine serpins have been described in sheep, cattle, and pigs as a highly diverged group of the large superfamily of serpin proteins that typically function as serine proteinase inhibitors. Here, the range of species that possess and express a uterine serpin gene is extended to the goat. Sequencing of cDNA amplified from total RNA from a pregnant goat at day 25 of pregnancy resulted in a 1,292 bp full-length consensus cDNA sequence for caprine uterine serpin (CaUS). The predicted amino acid sequence of the caprine precursor showed 96%, 82%, 55%, and 56% identity to OvUS, BoUS, PoUS1, and PoUS2, respectively. The signal peptide extends from amino acids 1 to 25, resulting in a secreted protein of 404 amino acids and 46,227 Mr (excluding carbohydrate). Both the goat and sheep uterine serpins have a nine amino acid insert in the Helix I region that is not found in bovine or porcine uterine serpins. A total of 13 amino acids in CaUS are different than those for the nearest homologue, ovine uterine serpin. One of these is in the site of cleavage of the signal sequence, where a single nucleotide substitution (G --> C) changed the cysteine for the sheep, bovine, and porcine genes to a serine. In addition, the amino acid at the putative P1-P1' site (the scissile bond for antiproteinase activity) is a valine for CaUS, BoUS, PoUS1, and PoUS2 versus an alanine for OvUS. The hinge region of all five of the uterine serpins (P17-P9) is distinct from the consensus pattern for inhibitory sequences and it is unlikely, therefore, that the uterine serpins possess prototypical proteinase inhibitory activity. The goat uterine serpin was immunolocalized to the glandular epithelium of the endometrium from a pregnant nanny at day 25 of pregnancy. There was also immunoreactive product in scattered luminal epithelial cells. No immunoreaction product was detected in endometrium from a nanny at day 5 of the estrous cycle. Western blotting of uterine fluid collected from the pregnant uterine horn of a unilaterally-pregnant goat revealed the presence of a protein band at Mr approximately 56,000 that reacted with monoclonal antibody to OvUS. In conclusion, the range of species in which uterine serpins are present and expressed in the uterus includes the goat in addition to the previously described sheep, cow, and pig. In all of these species, the uterine serpin is derived primarily from glandular epithelium, is secreted into the uterine lumen, and contains sequence characteristics suggesting it is not an inhibitory serpin.     

Evidence of positive selection at codon sites localized in the C-terminal peptide of ORC6

Origin recognition complex 6 (Orc6) plays a central role in the initiation of DNA replication in all eukaryotic systems. The exact contribution of Orc6 to replication initiation has yet to be elucidated. Here, we analyzed the evolutionary dynamics of Orc6 in 15 vertebrates. Positive selection was detected in the region of exon 6 of the Orc6 gene. Site tests revealed a proportion of codon sites that displayed evidence of positive selection (ω > 1) within the coding sequences of the vertebrate Orc6 gene. Seven positively selected amino acid sites were identified and three were located in exon6. These results suggest that amino acid residues present in the middle region of the protein are more selectively constrained, whereas amino acid residues in the C-terminal peptide of the protein evolve at a faster rate, possibly because of heightened selective pressure during the course of evolution.     

Phylogeny of the Ascaridoidea (Nematoda: Ascaridida) based on three genes and morphology: hypotheses of structural and sequence evolution

Ascaridoid nematodes parasitize the gastrointestinal tract of vertebrate definitive hosts and are represented by more than 50 described genera. We used 582 nucleotides (83% of the coding sequence) of the mitochondrial gene cytochrome oxidase subunit 2, in combination with published small- and large-subunit nuclear rDNA sequences (2,557 characters) and morphological data (20 characters), to produce a phylogenetic hypothesis for representatives of this superfamily. This combined evidence phylogeny strongly supported clades that, with 1 exception, were consistent with Fagerholm's 1991 classification. Parsimony mapping of character states on the combined evidence tree was used to develop hypotheses for the evolution of morphological, life history, and amino acid characters. This analysis of character evolution revealed that certain key features that have been used by previous workers for developing taxonomic and evolutionary hypotheses represent plesiomorphic states. Cytochrome oxidase subunit 2 nucleotides show a strong compositional bias to A+T and a substitution bias to thymine. These biases are most apparent at third positions of codons and 4-fold degenerate sites, which is consistent with the nonrandom substitution pattern of A+T pressure. Despite nucleotide bias, cytochrome oxidase amino acid sequences show conservation and retention of critical functional residues, as inferred from comparisons to other organisms.     

Characterization, primary structure, and evolution of lamprey plasma albumin.

The most abundant protein found in blood plasma from the sea lamprey (Petromyzon marinus) has the hallmarks of a plasma albumin: namely, high abundance, solubility in distilled water, a small number of tryptophans, and a high content of cysteines and charged residues. As in other vertebrate albumins, not all the cysteines are disulfide bonded. An unusual feature of this protein is its molecular weight of 175,000, roughly 2.5 times the size of other vertebrate albumins. Its amino acid sequence, deduced from a series of overlapping cDNA clones, can be aligned with other members of the gene family including plasma albumin, alpha-fetoprotein, and vitamin-D binding protein, confirming that it is indeed an oversized albumin. An unusual feature of the sequence is a 28-amino acid stretch consisting of a serine-threonine repeat with the general motif (STTT). Lamprey albumin contains a 23-amino acid putative signal peptide and a 6-residue putative propeptide, which, when cleaved, yield a mature protein of 1,394 amino acids with a calculated molecular weight of 157,000. The sequence also includes nine potential N-linked glycosylation sites (Asn-X-Ser/Thr), consistent with observation that lamprey albumin is a glycoprotein. If all the potential glycosylation sites were occupied by clusters of 2,000 molecular weight each, the total molecular weight would be 175,000. Like other members of the gene family, lamprey albumin is composed of a series of 190-amino acid repeats, there being seven such domains all together. Quantitative amino acid sequence comparisons of lamprey albumin with the other members of the gene family indicate that it diverged from an ancestral albumin prior to the gene duplications leading to this diverse group. This notion is confirmed by the pattern of amino acid insertions and deletions observed in a consideration of all domains that compose this family. Furthermore, it suggests that the invention of albumin antedates the vertebrate radiation.     

Molecular insights into Cumacean family relationships (Crustacea, Cumacea)

Cumaceans are a diverse order of small, benthic marine crustaceans. Phylogenetic hypotheses for the eight currently recognized cumacean families have not been formally proposed. However, based on external morphological traits and Linnean classification, a few conflicting hypotheses of relatedness have been proposed. Family definitions rely on morphological characters that often overlap and diagnoses are based on a combination of non-unique characters. Morphological analysis does not provide a well-resolved phylogeny. In the present study, we use amino acid sequences from the mitochondrial cytochrome oxidase I gene to produce a molecular phylogenetic hypothesis for the families of Cumacea. Phylogenetic analyses at the amino acid level were performed under Bayesian, likelihood, and parsimony methods. Results strongly suggest that families lacking an articulated telson form a monophyletic group. This pleotelson clade, composed of the families Bodotriidae, Leuconidae, and Nannastacidae, is the most derived within the Cumacea. Within this group, the Bodotriidae resolve paraphyletically, with Leuconidae and Nannastacidae embedded within it. Comparison of the molecular phylogeny with that based on morphology suggests that many "diagnostic" characters are homoplasious.