From DNA to fitness differences: sequences and structures of adaptive variants of Colias phosphoglucose isomerase (PGI)

Colias eurytheme butterflies display extensive allozyme polymorphism in the enzyme phosphoglucose isomerase (PGI). Earlier studies on biochemical and fitness effects of these genotypes found evidence of strong natural selection maintaining this polymorphism in the wild. Here we analyze the molecular features of this polymorphism by sequencing multiple alleles and modeling their structures. PGI is a dimer with rotational symmetry. Each monomer provides a critical residue to the other monomer's catalytic center. Sequenced alleles differ at multiple amino acid positions, including cryptic charge-neutral variation, but most consistent differences among the electromorph alleles are at the charge-changing amino acid sites. Principal candidate sites of selection, identified by structural and functional analyses and by their variants' population frequencies, occur in interpenetrating loops across the interface between monomers, where they may alter subunit interactions and catalytic center geometry. Comparison to a second (and basal) species, Colias meadii, also polymorphic for PGI under natural selection, reveals one fixed amino acid difference between their PGIs, which is located in the interpenetrating loop and accompanies functional differences among their variants. We also study nucleotide variability among the PGI alleles, comparing these data to similar data from another glycolytic enzyme gene, glyceraldehyde-3-phosphate dehydrogenase. Despite extensive nonsynonymous and synonymous polymorphism at PGI in each species, the only base changes fixed between species are the two causing the amino acid replacement; this absence of synonymous fixation yields a significant McDonald-Kreitman test. Analyses of these data suggest historical population expansion. Positive peaks of Tajima's D statistic, representing regions of neutral "hitchhiking," are found around the principal candidate sites of selection. This study provides novel views of molecular-structural mechanisms, and beginnings of historical evidence, for a long-persistent balanced enzyme polymorphism at PGI in these and perhaps other species.     

Comparison of Populations of Pratylenchus brachyurus Based on Isozyme Phenotypes

Enzymes from females of five Pratylenchus brachyurus populations and one P. scribneri population were analyzed by isoelectric focusing electrophoresis. Of the 18 enzyme systems investigated, only malate dehydrogenase (MDH), phosphoglucomutase (PGM), and phosphoglucose isomerase (PGI) were detected from all five P. brachyurus populations and P. scribneri. Faint bands were detected for isocitrate dehydrogenase and phosphogluconate dehydrogenase from one P. brachyurus population. Three distinct phenotypic groups were found in the MDH and PGM systems for P. brachyurus populations, but only a single electromorph was detected for PGI. Multiple electromorphs for MDH, PGM, and PGI were detected for P. scribneri; there was no similarity among these patterns with those from P. brachyurus. No phenotypic differences in PGI were observed between females and mixed juveniles of one population of P. brachyurus.     

Functional Effects of Pgi Allozymes in ESCHERICHIA COLI

Five alleles representing three electromorphs of phosphoglucose isomerase (PGI) have been transferred from natural isolates of E. coli into the genetic background of E. coli K12 and examined for their effect on growth rate in chemostats limited for glucose or fructose. With glucose limitation, all alleles are selectively neutral or nearly neutral within the limit of resolution of the technique, whether the genetic background is nonmutant or whether it contains a deletion of the locus of glucose-6-phosphate dehydrogenase, the enzyme that provides an alternative metabolic pathway for the substrate of PGI. With fructose limitation, one of the naturally occurring alleles has a small but reproducible detrimental effect on growth rate. A kinetic difference in this detrimental allozyme, apparently relating to an inhibition constant, has been observed in some, but not all, lots of substrate, and a similar difference has also been noted in one of the rare electromorphs that could not be transferred into E. coli K12. These results support a model of genetic variation in which the alleles are neutral or nearly neutral in the prevailing environment but have a potential for selection that can be expressed under the appropriate conditions of environment or genetic background. This hypothesis is discussed in the context of allozyme polymorphisms observed in other organisms.     

Hybrid fusions show that inter-monomer electron transfer robustly supports cytochrome bc1 function in vivo

Electronic connection between Q_o and Q_i quinone catalytic sites of dimeric cytochrome bc₁ is a central feature of the energy-conserving Q cycle. While both the intra- and inter-monomer electron transfers were shown to connect the sites in the enzyme, mechanistic and physiological significance of the latter remains unclear. Here, using a series of mutated hybrid cytochrome bc₁-like complexes, we show that inter-monomer electron transfer robustly sustains the function of the enzyme in vivo, even when the two subunits in a dimer come from different species. This indicates that minimal requirement for bioenergetic efficiency is to provide a chain of cofactors for uncompromised electron flux between the catalytic sites, while the details of protein scaffold are secondary.     

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.     

Low Levels of Polymorphisms and No Evidence for Diversifying Selection on the Plasmodium knowlesi Apical Membrane Antigen 1 Gene

Infection with Plasmodium knowlesi, a zoonotic primate malaria, is a growing human health problem in Southeast Asia. P. knowlesi is being used in malaria vaccine studies, and a number of proteins are being considered as candidate malaria vaccine antigens, including the Apical Membrane Antigen 1 (AMA1). In order to determine genetic diversity of the ama1 gene and to identify epitopes of AMA1 under strongest immune selection, the ama1 gene of 52 P. knowlesi isolates derived from human infections was sequenced. Sequence analysis of isolates from two geographically isolated regions in Sarawak showed that polymorphism in the protein is low compared to that of AMA1 of the major human malaria parasites, P. falciparum and P. vivax. Although the number of haplotypes was 27, the frequency of mutations at the majority of the polymorphic positions was low, and only six positions had a variance frequency higher than 10%. Only two positions had more than one alternative amino acid. Interestingly, three of the high-frequency polymorphic sites correspond to invariant sites in PfAMA1 or PvAMA1. Statistically significant differences in the quantity of three of the six high frequency mutations were observed between the two regions. These analyses suggest that the pkama1 gene is not under balancing selection, as observed for pfama1 and pvama1, and that the PkAMA1 protein is not a primary target for protective humoral immune responses in their reservoir macaque hosts, unlike PfAMA1 and PvAMA1 in humans. The low level of polymorphism justifies the development of a single allele PkAMA1-based vaccine.     

腾冲嗜热厌氧菌丙氨酸消旋酶底物通道氨基酸位点的功能

【目的】通过定点突变探究腾冲嗜热厌氧菌MB4中生物合成型丙氨酸消旋酶Tt Alr底物通道内氨基酸位点A172和S173的功能。【方法】利用定点突变PCR技术构建突变体,通过亲和层析法纯化酶蛋白,采用D-氨基酸氧化酶偶联法检测各突变蛋白的活性及其稳定性。【结果】通过定点突变PCR成功得到8个突变体,酶学特性分析发现,A172位点突变为丝氨酸(S)后酶蛋白的相对活性有所提升,但含有该位点突变的酶蛋白稳定性均大幅下降;S173位点突变为天门冬氨酸(D)后导致突变体蛋白的最适反应温度提升了15°C,半衰期大幅延长,但相对活性明显下降。【结论】丙氨酸消旋酶Tt Alr底物通道内A172和S173位点均是影响酶蛋白催化活性和稳定性的关键位点。     

A horizontally transferred nuclear gene is associated with microhabitat variation in a natural plant population

Horizontal gene transfer involves the non-sexual interspecific transmission of genetic material. Even if they are initially functional, horizontally transferred genes are expected to deteriorate into non-expressed pseudogenes, unless they become adaptively relevant in the recipient organism. However, little is known about the distributions of natural transgenes within wild species or the adaptive significance of natural transgenes within wild populations. Here, we examine the distribution of a natural plant-to-plant nuclear transgene in relation to environmental variation within a wild population. Festuca ovina is polymorphic for an extra (second) expressed copy of the nuclear gene (PgiC) encoding cytosolic phosphoglucose isomerase, with the extra PgiC locus having been acquired horizontally from the distantly related grass genus Poa. We investigated variation at PgiC in samples of F. ovina from a fine-scale, repeating patchwork of grassland microhabitats, replicated within spatially separated sites. Even after accounting for spatial effects, the distributions of F. ovina individuals carrying the additional PgiC locus, and one of the enzyme products encoded by the locus, are significantly associated with fine-scale habitat variation. Our results suggest that the PgiC transgene contributes, together with the unlinked ‘native’ PgiC locus, to local adaptation to a fine-scale mosaic of edaphic and biotic grassland microhabitats.     

A precise method for the quantitation of proteins taking into account their amino acid composition.

A method for the quantitation of protein in biological material is described which gives the same response for all proteins irrespective of their amino acid composition. The method is based on the ninhydrin reaction of amino acids released after total acid hydrolysis of 5- to 20-μl solutions containing 1 to 100 μg of protein. The ammonia is released from the hydrolysate by diffusion and the amino acids are quantitated without fractionation using the continuous-flow system of an amino acid analyzer. Calibration is obtained with solutions of known amino acid content. The protein of a sample is calculated by multiplying the nanomoles of total amino acids found by a conversion factor F. F is the weight in micrograms of 1 nmol of the specific mixture of amino acid residues that the protein of the sample is composed of F has to be determined once for all further quantitations of the same material by quantitative amino acid analysis following standard procedures. By this method as little as 30 ng of protein per aliquot of hydrolysate analyzed can be determined.     

Multiple electrophoretic variants of Cu,Zn superoxide dismutase as expression of the enzyme aging. Effects of H2O2, ascorbate and metal ions

Multiple electrophoretic bands, with R_F identical to the natural molecular variants, are produced by treatment of purified Cu, Zn Superoxide dismutase with either H₂O₂ or ascorbate plus Fe(III) EDTA. The ascorbate reaction is also due to H₂O₂ since it is inhibited by catalase. However while H₂O₂ inactivates the enzyme, the electromorphs produced by ascorbate-Fe(III) EDTA have only slightly less activity than the native enzyme and this property parallels the natural situation. It is concluded that oxidative aging can be responsible for the multiple molecular variants of the natural enzyme, under conditions where the oxidant attack is preferentially directed to amino acid side chains outside the active site. Such conditions may occur when a metal ion coordinated to the protein surface undergoes a redox cycle with biological reductants, like ascorbate.     

A unique enzyme of acetic acid bacteria, PQQ-dependent alcohol dehydrogenase, is also present in Frateuria aurantia

A membrane-bound, pyrroloquinoline quinone (PQQ)-dependent alcohol dehydrogenase (ADH) was purified from Frateuria aurantia LMG 1558^T. Although F. aurantia belongs to a group of γ-Proteobacteria, the characteristics of its PQQ-ADH were similar to the enzyme characteristics of the typical high-acetic acid-resistant bacterium Gluconacetobacter europaeus from the group of α-Proteobacteria. The PQQ-dependent ADH was solubilized from the membranes and purified after anionic, cationic, and affinity chromatography with specific activity of 117 U/mg. The purified enzyme was estimated to be composed of two subunits of ca. 72 and 45 kDa, as judged by SDS-polyacrylamide gel electrophoresis. The purified enzyme had maximum activity at pH 4.5 and showed the highest substrate specificity to ethanol, isoamyl alcohol, 1-butanol, and 1-propanol. The deduced sequences of cloned genes adhA and adhB encoding subunits I and II of PQQ-ADH showed 80 % amino acid (AA) identity to AdhA and 68 % AA identity to AdhB of Ga. europaeus V3 (LMG 18494). Because of the high similarity between genes encoding subunits I and II of PQQ-ADH and its homologous genes found in a distantly related taxonomic group of acetic acid bacteria, the results suggest the possibility of horizontal gene transfer between these two groups of genera.     

Electric charge divergence in proteins: insights into the evolution of their three-dimensional properties

Previous studies of protein evolution have identified important mutations in various proteins that affect a small number of residues, but dramatically alter protein function. However, the evolutionary process underlying the three-dimensional protein properties, which are determined by a much larger number of residues, remains unclear. Based on a comparative evolutionary analysis of teleost phosphoglucose isomerases (PGIs; EC 5.3.1.9), we previously demonstrated that the relatively weak selection on many amino acid sites has played an important role in the evolution of protein electric charge as a model of three-dimensional protein properties. To ascertain the generality of this finding, we sought further evidence of this type of protein evolution. For this purpose, we analyzed the vertebrate isoforms of fructose-1,6-bisphosphate aldolase (ALD; EC 4.1.2.13), for which electric charges are known to have diverged after gene duplication. The results showed that the divergence in electric charge between the ALD isoforms was also driven by weak selection on many amino acid sites, as in PGI, confirming the generality of earlier findings. To obtain further insights, ALD and PGI were compared to the proteins pancreatic ribonuclease (EC 3.1.27.5) and triose-phosphate isomerase (EC 5.3.1.1), for which electric charges likely evolved through a well-defined mode of molecular evolution; namely, strong selection on specific amino acid sites. Comparison of the number and composition of amino acids on the protein surface suggested that the absolute number of evolutionarily changeable amino acids in a protein affects the strength of selection pressure acting on individual amino acid sites.     

Inter-monomer electron transfer is too slow to compete with monomeric turnover in bc1 complex

The homodimeric bc₁ complexes are membrane proteins essential in respiration and photosynthesis. The ~ 11 Å distance between the two b_L-hemes of the dimer opens the possibility of electron transfer between them, but contradictory reports make such inter-monomer electron transfer controversial. We have constructed in Rhodobacter sphaeroides a heterodimeric expression system similar to those used before, in which the bc₁ complex can be mutated differentially in the two copies of cyt b to test for inter-monomer electron transfer, but found that genetic recombination by cross-over then occurs to produce wild-type homodimer. Selection pressure under photosynthetic growth always favored the homodimer over heterodimeric variants enforcing inter-monomer electron transfer, showing that the latter are not competitive. These results, together with kinetic analysis of myxothiazol titrations, demonstrate that inter-monomer electron transfer does not occur at rates competitive with monomeric turnover. We examine the results from other groups interpreted as demonstrating rapid inter-monomer electron transfer, conclude that similar mechanisms are likely to be in play, and suggest that such claims might need to be re-examined.     

Evolutionary Relationships and Sequence Variation of α-Amylase Variants Encoded by Duplicated Genes in the Amy Locus of Drosophila Melanogaster

To infer the genealogical relationships of α-amylase electromorphs of Drosophila melanogaster, we determined the nucleotide sequences of a collection of electromorphs sampled throughout the world. On average there were 1.0 amino acid substitutions between identical electromorphs and 3.9 between different electromorphs, respectively. We found that the evolution of AMY(1) through AMY(6) electromorphs occurred by sequential accumulation of single amino acid substitutions each causing one charge difference. The nucleotide diversities at synonymous sites within Amy(1),Amy(2),Amy(3),Amy(4) and Amy(6) were 0.0321, 0.0000, 0.0355, 0.0059 and 0.0030, respectively. We also obtained evidence of genetic exchanges, such as intrachromosomal recombination, interchromosomal recombination or gene conversion, between the two duplicated Amy genes as well as among the alleles.     

Genomic Characterization of Phenylalanine Ammonia Lyase Gene in Buckwheat

Phenylalanine Ammonia Lyase (PAL) gene which plays a key role in bio-synthesis of medicinally important compounds, Rutin/quercetin was sequence characterized for its efficient genomics application. These compounds possessing anti-diabetic and anti-cancer properties and are predominantly produced by Fagopyrum spp. In the present study, PAL gene was sequenced from three Fagopyrum spp. (F. tataricum, F. esculentum and F. dibotrys) and showed the presence of three SNPs and four insertion/deletions at intra and inter specific level. Among them, the potential SNP (position 949^th bp G>C) with Parsimony Informative Site was selected and successfully utilised to individuate the zygosity/allelic variation of 16 F. tataricum varieties. Insertion mutations were identified in coding region, which resulted the change of a stretch of 39 amino acids on the putative protein. Our Study revealed that autogamous species (F. tataricum) has lower frequency of observed SNPs as compared to allogamous species (F. dibotrys and F. esculentum). The identified SNPs in F. tataricum didn’t result to amino acid change, while in other two species it caused both conservative and non-conservative variations. Consistent pattern of SNPs across the species revealed their phylogenetic importance. We found two groups of F. tataricum and one of them was closely related with F. dibotrys. Sequence characterization information of PAL gene reported in present investigation can be utilized in genetic improvement of buckwheat in reference to its medicinal value.     

Molecular Interaction of Novel Compound 2-Methylheptyl Isonicotinate Produced by Streptomyces sp. 201 with Dihydrodipicolinate Synthase (DHDPS) Enzyme of Mycobacterium tuberculosis for its Antibacterial Activity

Antibiotic resistance is a growing problem in multi-drug-resistant tuberculosis which is caused by Mycobacterium tuberculosis (MTB). Hence there is an urgent need for designing or developing a novel or potent anti-tubercular agent. The Lysine/DAP biosynthetic pathway is a promising target because of its role in cell wall and amino acid biosynthesis. In our study we performed a molecular docking analysis of a novel antibacterial isolated from Streptomyces sp. 201 at three different binding site of dihydrodipicolinate synthase (DHDPS) enzyme of MTB. The molecular docking studies suggest that the novel molecule shows favourable interaction at the three different binding sites as compared to five experimentally known inhibitors of DHDPS.     

Selection shapes the patterns of codon usage in three closely related species of genus Misgurnus

Neutrality plots revealed that selection probably dominates codon bias, whereas mutation plays only a minor role, in shaping the codon bias in three loaches, Misgurnus anguillicaudatus, M. mohoity, and M. bipartitus. These three species also clearly showed similar tendencies in the preferential usage of codons. Nineteen, nine, and 14 preferred codon pairs and 179, 182, and 173 avoided codon pairs were also detected in M. anguillicaudatus, M. bipartitus, and M. mohoity, respectively, and the most frequently avoided type of cP3-cA1 dinucleotide in these species was nnUAnn. The expression-linked patterns of codon usage revealed that higher expression was associated with higher GC₃, lower ENC, and a smaller proportion of amino acids with high size/complexity (S/C) scores in these three species. These results elucidate selectively driven codon bias in Misgurnus species, and reveal the potential importance of expression-mediated selection in shaping the genome evolution of fish.     

The Evolution of Vp1 Gene in Enterovirus C Species Sub-Group That Contains Types CVA-21, CVA-24, EV-C95, EV-C96 and EV-C99

Genus Enterovirus (Family Picornaviridae,) consists of twelve species divided into genetically diverse types by their capsid protein VP1 coding sequences. Each enterovirus type can further be divided into intra-typic sub-clusters (genotypes). The aim of this study was to elucidate what leads to the emergence of novel enterovirus clades (types and genotypes). An evolutionary analysis was conducted for a sub-group of Enterovirus C species that contains types Coxsackievirus A21 (CVA-21), CVA-24, Enterovirus C95 (EV-C95), EV-C96 and EV-C99. VP1 gene datasets were collected and analysed to infer the phylogeny, rate of evolution, nucleotide and amino acid substitution patterns and signs of selection. In VP1 coding gene, high intra-typic sequence diversities and robust grouping into distinct genotypes within each type were detected. Within each type the majority of nucleotide substitutions were synonymous and the non-synonymous substitutions tended to cluster in distinct highly polymorphic sites. Signs of positive selection were detected in some of these highly polymorphic sites, while strong negative selection was indicated in most of the codons. Despite robust clustering to intra-typic genotypes, only few genotype-specific ‘signature’ amino acids were detected. In contrast, when different enterovirus types were compared, there was a clear tendency towards fixation of type-specific ‘signature’ amino acids. The results suggest that permanent fixation of type-specific amino acids is a hallmark associated with evolution of different enterovirus types, whereas neutral evolution and/or (frequency-dependent) positive selection in few highly polymorphic amino acid sites are the dominant forms of evolution when strains within an enterovirus type are compared.     

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.