Natural selection drives extremely rapid evolution in antiviral RNAi genes

RNA interference (RNAi) is perhaps best known as a laboratory tool. However, RNAi-related pathways represent an antiviral component of innate immunity in both plants and animals. Since viruses can protect themselves by suppressing RNAi, interaction between RNA viruses and host RNAi may represent an ancient coevolutionary "arms race." This could lead to strong directional selection on RNAi genes, but to date their evolution has not been studied. By comparing DNA sequences from different species of Drosophila, we show that the rate of amino acid evolution is substantially elevated in genes related to antiviral RNAi function (Dcr2, R2D2, and Ago2). They are among the fastest evolving 3% of all Drosophila genes; they evolve significantly faster than other components of innate immunity and faster than paralogous genes that mediate "housekeeping" functions. Based on DNA polymorphism data from three species of Drosophila, McDonald-Kreitman tests showed that this rapid evolution is due to strong positive selection. Furthermore, Dcr2 and Ago2 display reduced genetic diversity, indicative of a recent selective sweep in both genes. Together, these data show rapid adaptive evolution of the antiviral RNAi pathway in Drosophila. This is a signature of host-pathogen arms races and implies that the ancient battle between RNA viruses and host antiviral RNAi genes is active and significant in shaping RNAi function.     

Positive selection drives lactoferrin evolution in mammals

Lactoferrin (LF) is a member of the transferrin family that is abundantly expressed and secreted by glandular epithelial cells. The biological functions of LF involve in iron homeostasis regulation of the body and antibacterial activity. Previous studies demonstrated that it had a high cationic N-terminal domain that could interact with glycosaminoglycans, lipopolysaccharides and the bacterial virulence protein. Two anti-microbial peptides, lactoferricin (LFcin) and lactoferrampin (LFampin), were also isolated and identified in N-terminal of LF. Although the antibacterial mechanism was carefully studied, little was known about the molecular evolution of LF. In this study, we estimated the nonsynonymous-to-synonymous substitution ratios ( w = d_N \mathord

Deduced amino acid sequence of Escherichia coli adenosine deaminase reveals evolutionarily conserved amino acid residues: implications for catalytic function

The goal of the research reported here is to identify evolutionarily conserved amino acid residues associated with enzymatic deamination of adenosine. To do this, we isolated molecular clones of the Escherichia coli adenosine deaminase gene by functional complementation of adenosine deaminase deficient bacteria and deduced the amino acid sequence of the enzyme from the nucleotide sequence of the gene. Nucleotide sequence analysis revealed the presence of a 996-nucleotide open reading frame encoding a protein of 332 amino acids having a molecular weight of 36,345. The deduced amino acid sequence of the E. coli enzyme has approximately 33% identity with those of the mammalian adenosine deaminases. With conservative amino acid substitutions the overall sequence homology approaches 50%, suggesting that the structures and functions of the mammalian and bacterial enzymes are similar. Additional amino acid sequence analysis revealed specific residues that are conserved among all three adenosine deaminases and four AMP deaminases for which sequence information is currently available. In view of previously published enzymological data and the conserved amino acid residues identified in this study, we propose a model to account for the enzyme-catalyzed hydrolytic deamination of adenosine. Potential catalytic roles are assigned to the conserved His 214, Cys 262, Asp 295, and Asp 296 residues of mammalian adenosine deaminases and the corresponding conserved amino acid residues in bacterial adenosine deaminase and the eukaryotic AMP deaminases.     

Molecular characterization of HymA, an evolutionarily highly conserved and highly expressed protein of Aspergillus nidulans

Aspergillus nidulans reproduces asexually via uninucleate, haploid spores, which are produced on morphologically differentiated aerial structures, called conidiophores. These consist of four distinct cell types, a foot with a terminally swollen stalk, metulae, phialides and conidiospores. The molecular mechanisms underlying the morphological changes that occur during conidiophore development have been studied by mutant analysis. We have isolated the hymA mutant, in which conidiophore development is affected at the metula stage. In the mutant metulae do not differentiate properly but come to resemble hyphae (hym?=?hypha-like metulae). In this paper we have analyzed the corresponding gene. It encodes a highly expressed 44?kDa protein which resides in the cytoplasm and has homologues in yeast, plants, fly, worm, fish, mice and man. We constructed hym deletion strains of Saccharomyces cerevisiae and of A. nidulans and found that the gene is essential in S. cerevisiae but is dispensable in the filamentous fungus. A cellular function for the Hym protein has not yet been defined in any organism. To demonstrate functional conservation we constructed a chimeric protein comprised of the N-terminal half of the A.?nidulans and the C-terminal half of the mouse homologue MO25. This hybrid protein could fully substitute for HymA function in A. nidulans. In addition, the mouse protein itself partially rescued the hymA mutation in the fungus. HymA is thus highly conserved in evolution and probably serves similar functions. The fact that hymA is required for conidiophore development in A. nidulans suggests that homologous genes in other organisms might also be involved in morphogenesis.     

Single amino acid substitutions at conserved residues of human interferon-alpha can effect antiviral specific activity

Site-specific in vitro mutagenesis was used to direct various amino acid substitutions at conserved positions within the sequence of human interferon-alpha 1 (IFN-alpha 1). The antiviral specific activity of IFN-alpha 1, expressed in M13 as a fusion protein [IFN-alpha 1 (phi WT)], could be altered by single amino acid substitutions. The substitution of glycine for tyrosine at position 123 results in a loss of more than 99% of the antiviral specific activity on human cells, but causes no significant change in the antiviral specific activity on primary bovine cells. The tyrosine at position 123 is thus implicated in determining human cell specificity. Based on analysis of IFN-alpha 2, IFN-alpha 1 contains two dulsulphide bridges between cysteine residues 29 and 139 and cysteine residues 1 and 99. IFN-alpha 1 also contains a fifth cysteine residue at position 86. IFN-alpha 1 (phi WT) carrying three serine for cysteine substitutions at positions 1, 86 and 99 retains 23% of the antiviral specific activity of IFN-alpha 1 (phi WT) on human cells. However, the antiviral activity on bovine cells is not significantly affected by this modification. The presence of the disulphide bridge between residues 1 and 99 thus appears to be required for full antiviral activity on human but not bovine cells. A single serine for cysteine substitution at position 29 reduces the antiviral specific activity on both human and bovine cells by some 95%. This data shows that the disulphide bridge between residues 29 and 139 is critical for the antiviral activity of IFN-alpha's.     

Functional significance of conserved amino acid residues.

A systemic study of single amino acid substitutions in bacteriophage T4 lysozyme permitted a test of the concept that conserved amino acid residues are more functionally important than nonconserved residues. Substitutions of amino acid residues that are conserved among five bacteriophage-encoded lysozymes were found to lead more frequently to loss of function than substitutions of nonconserved residues. Of 163 residues tested, only 74 (45%) are sensitive to at least one substitution; however, all 14 residues that are fully conserved are sensitive to substitutions.     

Molecular characterization of HymA, an evolutionarily highly conserved and highly expressed protein of Aspergillus nidulans

Aspergillus nidulans reproduces asexually via uninucleate, haploid spores, which are produced on morphologically differentiated aerial structures, called conidiophores. These consist of four distinct cell types, a foot with a terminally swollen stalk, metulae, phialides and conidiospores. The molecular mechanisms underlying the morphological changes that occur during conidiophore development have been studied by mutant analysis. We have isolated the hymA mutant, in which conidiophore development is affected at the metula stage. In the mutant metulae do not differentiate properly but come to resemble hyphae (hym = hypha-like metulae). In this paper we have analyzed the corresponding gene. It encodes a highly expressed 44 kDa protein which resides in the cytoplasm and has homologues in yeast, plants, fly, worm, fish, mice and man. We constructed hym deletion strains of Saccharomyces cerevisiae and of A. nidulans and found that the gene is essential in S. cerevisiae but is dispensable in the filamentous fungus. A cellular function for the Hym protein has not yet been defined in any organism. To demonstrate functional conservation we constructed a chimeric protein comprised of the N-terminal half of the A.␣nidulans and the C-terminal half of the mouse homologue MO25. This hybrid protein could fully substitute for HymA function in A. nidulans. In addition, the mouse protein itself partially rescued the hymA mutation in the fungus. HymA is thus highly conserved in evolution and probably serves similar functions. The fact that hymA is required for conidiophore development in A. nidulans suggests that homologous genes in other organisms might also be involved in morphogenesis. Received: 11 February 1998 / Accepted: 14 September 1998     

Identification of structural and catalytic classes of highly conserved amino acid residues in lysine 2,3-aminomutase

Lysine 2,3-aminomutase (LAM) from Clostridium subterminale SB4 catalyzes the interconversion of (S)-lysine and (S)-beta-lysine by a radical mechanism involving coenzymatic actions of S-adenosylmethionine (SAM), a [4Fe-4S] cluster, and pyridoxal 5'-phosphate (PLP). The enzyme contains a number of conserved acidic residues and a cysteine- and arginine-rich motif, which binds iron and sulfide in the [4Fe-4S] cluster. The results of activity and iron, sulfide, and PLP analysis of variants resulting from site-specific mutations of the conserved acidic residues and the arginine residues in the iron-sulfide binding motif indicate two classes of conserved residues of each type. Mutation of the conserved residues Arg134, Asp293, and Asp330 abolishes all enzymatic activity. On the basis of the X-ray crystal structure, these residues bind the epsilon-aminium and alpha-carboxylate groups of (S)-lysine. However, among these residues, only Asp293 appears to be important for stabilizing the [4Fe-4S] cluster. Members of a second group of conserved residues appear to stabilize the structure of LAM. Mutations of arginine 130, 135, and 136 and acidic residues Glu86, Asp165, Glu236, and Asp172 dramatically decrease iron and sulfide contents in the purified variants. Mutation of Asp96 significantly decreases iron and sulfide content. Arg130 or Asp172 variants display no detectable activity, whereas variants mutated at the other positions display low to very low activities. Structural roles are assigned to this latter class of conserved amino acids. In particular, a network of hydrogen bonded interactions of Arg130, Glu86, Arg135, and the main chain carbonyl groups of Cys132 and Leu55 appears to stabilize the [4Fe-4S] cluster.     

Positive selection drives the evolution of the Acp29AB accessory gland protein in Drosophila.

Nucleotide sequence variation at the Acp29AB gene region has been surveyed in Drosophila melanogaster from Spain (12 lines), Ivory Coast (14 lines), and Malawi (13 lines) and in one line of D. simulans. The approximately 1.7-kb region studied encompasses the Acp29AB gene that codes for a male accessory gland protein and its flanking regions. Seventy-seven nucleotide and 8 length polymorphisms were detected. Nonsynonymous polymorphism was an order of magnitude lower than synonymous polymorphism, but still high relative to other non-sex-related genes. In D. melanogaster variation at this region revealed no major genetic differentiation between East and West African populations, while differentiation was highly significant between the European and the two African populations. Comparison of polymorphism and divergence at synonymous and nonsynonymous sites showed an excess of fixed nonsynonymous changes, which indicates that the evolution of the Acp29AB protein has been driven by directional selection at least after the split of the D. melanogaster and D. simulans lineages. The pattern of variation in extant populations of D. melanogaster favors a scenario where the fixation of advantageous replacement substitutions occurred in the early stages of speciation and balancing selection is maintaining variation in this species.     

Structure-function studies of nerve growth factor: functional importance of highly conserved amino acid residues.

Selected amino acid residues in chicken nerve growth factor (NGF) were replaced by site-directed mutagenesis. Mutated NGF sequences were transiently expressed in COS cells and the yield of NGF protein in conditioned medium was quantified by Western blotting. Binding of each mutant to NGF receptors on PC12 cells was evaluated in a competition assay. The biological activity was determined by measuring stimulation of neurite outgrowth from chick sympathetic ganglia. The residues homologous to the proposed receptor binding site of insulin (Ser18, Met19, Val21, Asp23) were substituted by Ala. Replacement of Ser18, Met19 and Asp23 did not affect NGF activity. Modification of Val21 notably reduced both receptor binding and biological activity, suggesting that this residue is important to retain a fully active NGF. The highly conserved Tyr51 and Arg99 were converted into Phe and Lys respectively, without changing the biological properties of the molecule. However, binding and biological activity were greatly impaired after the simultaneous replacement of both Arg99 and Arg102 by Gly. The three conserved Trp residues at positions 20, 75 and 98 were substituted by Phe. The Trp mutated proteins retained 15-60% of receptor binding and 40-80% of biological activity, indicating that the Trp residues are not essential for NGF activity. However, replacement of Trp20 significantly reduced the amount of NGF in the medium, suggesting that this residue may be important for protein stability.     

Positive selection drives the evolution of rhino, a member of the heterochromatin protein 1 family in Drosophila

Heterochromatin comprises a significant component of many eukaryotic genomes. In comparison to euchromatin, heterochromatin is gene poor, transposon rich, and late replicating. It serves many important biological roles, from gene silencing to accurate chromosome segregation, yet little is known about the evolutionary constraints that shape heterochromatin. A complementary approach to the traditional one of directly studying heterochromatic DNA sequence is to study the evolution of proteins that bind and define heterochromatin. One of the best markers for heterochromatin is the heterochromatin protein 1 (HP1), which is an essential, nonhistone chromosomal protein. Here we investigate the molecular evolution of five HP1 paralogs present in Drosophila melanogaster. Three of these paralogs have ubiquitous expression patterns in adult Drosophila tissues, whereas HP1D/rhino and HP1E are expressed predominantly in ovaries and testes respectively. The HP1 paralogs also have distinct localization preferences in Drosophila cells. Thus, Rhino localizes to the heterochromatic compartment in Drosophila tissue culture cells, but in a pattern distinct from HP1A and lysine-9 dimethylated H3. Using molecular evolution and population genetic analyses, we find that rhino has been subject to positive selection in all three domains of the protein: the N-terminal chromo domain, the C-terminal chromo-shadow domain, and the hinge region that connects these two modules. Maximum likelihood analysis of rhino sequences from 20 species of Drosophila reveals that a small number of residues of the chromo and shadow domains have been subject to repeated positive selection. The rapid and positive selection of rhino is highly unusual for a gene encoding a chromosomal protein and suggests that rhino is involved in a genetic conflict that affects the germline, belying the notion that heterochromatin is simply a passive recipient of "junk DNA" in eukaryotic genomes.     

Post-duplication charge evolution of phosphoglucose isomerases in teleost fishes through weak selection on many amino acid sites

Background  

The partitioning of ancestral functions among duplicated genes by neutral evolution, or subfunctionalization, has been considered the primary process for the evolution of novel proteins (neofunctionalization). Nonetheless, how a subfunctionalized protein can evolve into a more adaptive protein is poorly understood, mainly due to the limitations of current analytical methods, which can detect only strong selection for amino acid substitutions involved in adaptive molecular evolution. In this study, we employed a comparative evolutionary approach to this question, focusing on differences in the structural properties of a protein, specifically the electric charge, encoded by fish-specific duplicated phosphoglucose isomerase (Pgi) genes.     

Highly conserved configuration of catalytic amino acid residues among calicivirus-encoded proteases

A common feature of caliciviruses is the proteolytic processing of the viral polyprotein catalyzed by the viral 3C-like protease encoded in open reading frame 1 (ORF1). Here we report the identification and structural characterization of the protease domains and amino acid residues in sapovirus (SaV) and feline calicivirus (FCV). The in vitro expression and processing of a panel of truncated ORF1 polyproteins and corresponding mutant forms showed that the functional protease domain is 146 amino acids (aa) in SaV and 154 aa in FCV. Site-directed mutagenesis of the protease domains identified four amino acid residues essential to protease activities: H(31), E(52), C(116), and H(131) in SaV and H(39), E(60), C(122), and H(137) in FCV. A computer-assisted structural analysis showed that despite high levels of diversity in the primary structures of the protease domains in the family Caliciviridae, the configurations of the H, E, C, and H residues are highly conserved, with these residues positioned closely along the inner surface of the potential binding cleft for the substrate. These results strongly suggest that the H, E, C, and H residues are involved in the formation of a conserved catalytic surface of the SaV and FCV 3C-like proteases.     

Identification of highly conserved hydrophobic amino acid motif in deduced amino acid sequence of Elymus sibiricus L. mitochondrial S13 ribosomal protein

The nucleotide sequence of mitochondrial ribosomal protein rps13 gene from wild perennial grass Elymus sibiricus is presented. It was determined by the method of PCR amplification with specific oligonucleotide primers and the direct sequencing of the amplification product. The sequence of E. sibiricus mitochondrial gene for S13 predicts a hydrophobic ribosomal protein of 116 amino acids that shows strong similarity to those of wheat (99.7% identity) and maize (98%). The deduced amino acid sequence of S13 protein from E. sibiricus and homologous plant's (Zea mays, Daucus carota, Nicotiana tabacum, Marchantia polymorpha) and nonplant's (Escherichia coli) proteins shows the presence of hydrophobic amino acids' motif -L-X10-L-X10-M-X10-L-X10-L-. Slightly modified it can be found in many other ribosomal proteins. This conserved motif is presumed to be particularly important for association of the ribosomal S13 protein with other proteins in the small subunit of the mitochondrial ribosome.     

Simulating evolution by gene duplication of protein features that require multiple amino acid residues

Gene duplication is thought to be a major source of evolutionary innovation because it allows one copy of a gene to mutate and explore genetic space while the other copy continues to fulfill the original function. Models of the process often implicitly assume that a single mutation to the duplicated gene can confer a new selectable property. Yet some protein features, such as disulfide bonds or ligand binding sites, require the participation of two or more amino acid residues, which could require several mutations. Here we model the evolution of such protein features by what we consider to be the conceptually simplest route-point mutation in duplicated genes. We show that for very large population sizes N, where at steady state in the absence of selection the population would be expected to contain one or more duplicated alleles coding for the feature, the time to fixation in the population hovers near the inverse of the point mutation rate, and varies sluggishly with the lambda(th) root of 1/N, where lambda is the number of nucleotide positions that must be mutated to produce the feature. At smaller population sizes, the time to fixation varies linearly with 1/N and exceeds the inverse of the point mutation rate. We conclude that, in general, to be fixed in 10(8) generations, the production of novel protein features that require the participation of two or more amino acid residues simply by multiple point mutations in duplicated genes would entail population sizes of no less than 10(9).     

An evolutionarily conserved positively charged amino acid in the putative membrane-spanning domain of the foamy virus envelope protein controls fusion activity

Foamy viruses (FVs) are highly fusogenic, and their replication induces massive syncytium formation in infected cell cultures which is believed to be mediated by expression of the envelope (Env) protein. The FV Env is essential for virus particle egress. The unusually long putative membrane-spanning domain (MSD) of the transmembrane subunit carries dispersed charged amino acids and has an important function for particle envelopment. To better understand the capsid-envelope interaction and Env-mediated cell fusion, we generated a variety of FV MSD mutations. C-terminal deletions revealed the cytoplasmic domain to be dispensable but the full-length MSD to be required for fusogenic activity. The N-terminal 15 amino acids of the MSD were found to be sufficient for membrane anchorage and promotion of FV particle release. Expression of wild-type Env protein rarely induced syncytia due to intracellular retention. Coexpression with FV Gag-Pol resulted in particle export and a dramatic increase in fusion activity. A nonconservative mutation of K(959) in the middle of the putative MSD resulted in increased fusogenic activity of Env in the absence of Gag-Pol due to enhanced cell surface expression as well as structural changes in the mutant proteins. Coexpression with Gag-Pol resulted in a further increase in the fusion activity of mutant FV Env proteins. Our results suggest that an interaction between the viral capsid and Env is required for FV-induced giant-cell formation and that the positive charge in the MSD is an important determinant controlling intracellular transport and fusogenic activity of the FV Env protein.     

Tight junction protein Par6 interacts with an evolutionarily conserved region in the amino terminus of PALS1/stardust

Tight junctions are the structures in mammalian epithelial cells that separate the apical and basolateral membranes and may also be important in the establishment of cell polarity. Two evolutionarily conserved multiprotein complexes, Crumbs-PALS1 (Stardust)-PATJ and Cdc42-Par6-Par3-atypical protein kinase C, have been implicated in the assembly of tight junctions and in polarization of Drosophila melanogaster epithelia. These two complexes have been linked physically and functionally by an interaction between PALS1 and Par6. Here we identify an evolutionarily conserved region in the amino terminus of PALS1 as the Par6 binding site and identify valine and aspartic acid residues in this region as essential for interacting with the PDZ domain of Par6. We have also characterized, in more detail, the amino terminus of Drosophila Stardust and demonstrate that the interaction mechanism between Stardust and Drosophila Par6 is evolutionarily conserved. Par6 interferes with PATJ in binding PALS1, and these two interactions do not appear to function synergistically. Taken together, these results define the molecular mechanisms linking two conserved polarity complexes.     

Structural organization of essential iron-sulfur clusters in the evolutionarily highly conserved ATP-binding cassette protein ABCE1

The ABC protein ABCE1, formerly named RNase L inhibitor RLI1, is one of the most conserved proteins in evolution and is expressed in all organisms except eubacteria. Because of its fundamental role in translation initiation and/or ribosome biosynthesis, ABCE1 is essential for life. Its molecular mechanism has, however, not been elucidated. In addition to two ABC ATPase domains, ABCE1 contains a unique N-terminal region with eight conserved cysteines, predicted to coordinate iron-sulfur clusters. Here we present detailed information on the type and on the structural organization of the Fe-S clusters in ABCE1. Based on biophysical, biochemical, and yeast genetic analyses, ABCE1 harbors two essential diamagnetic [4Fe-4S](2+) clusters with different electronic environments, one ferredoxin-like (CPX(n)CX(2)CX(2)C; Cys at positions 4-7) and one unique ABCE1-type cluster (CXPX(2)CX(3)CX(n)CP; Cys at positions 1, 2, 3, and 8). Strikingly, only seven of the eight conserved cysteines coordinating the Fe-S clusters are essential for cell viability. Mutagenesis of the cysteine at position 6 yielded a functional ABCE1 with the ferredoxin-like Fe-S cluster in a paramagnetic [3Fe-4S](+) state. Notably, a lethal mutation of the cysteine at position 4 can be rescued by ligand swapping with an adjacent, extra cysteine conserved among all eukaryotes.