Selective Pressures on <Emphasis Type="Italic">Drosophila</Emphasis> Chemosensory Receptor Genes

The evolution and patterns of selection of genes encoding 10 Drosophila odorant receptors (Or) and the sex pheromone receptor Gr68a were investigated by comparing orthologous sequences across five to eight ecologically diverse species of Drosophila. Using maximum likelihood estimates of dN/dS ratios we show that all 11 genes sampled are under purifying selection, indicating functional constraint. Four of these genes (Or33c, Or42a, Or85e, and Gr68a) may be under positive selection, and if so, there is good evidence that 12 specific amino acid sites may be under positive selection. All of these sites are predicted to be located either in loop regions or just inside membrane spanning regions, and interestingly one of the two sites in Gr68a is in a similar position to a previously described polymorphism in Gr5a that causes a shift in sensitivity to its ligand trehalose. For three Ors, possible evidence for positive selection was detected along a lineage. These include Or22a in the lineage leading to D. mauritiana and Or22b in the lineage leading to D. simulans. This is of interest in light of previous data showing a change in ligand response profile for these species in the sensory neuron (ab3A) which expresses both Or22a and Or22b in D. melanogaster. In summary, while the main chemosensory function and/or structural integrity of these 10 Or genes and Gr68a are evolutionarily preserved, positive selection appears to be acting on some of these genes, at specific sites and along certain lineages, and provides testable hypotheses for further functional experimentation. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users. [Reviewing Editor: Dr. David Pollock]     

Polymorphism and selection in the major histocompatibility complex <Emphasis Type="Italic">DRA</Emphasis> and <Emphasis Type="Italic">DQA</Emphasis> genes in the family Equidae

The major histocompatibility complex genes coding for antigen binding and presenting molecules are the most polymorphic genes in the vertebrate genome. We studied the DRA and DQA gene polymorphism of the family Equidae. In addition to 11 previously reported DRA and 24 DQA alleles, six new DRA sequences and 13 new DQA alleles were identified in the genus Equus. Phylogenetic analysis of both DRA and DQA sequences provided evidence for trans-species polymorphism in the family Equidae. The phylogenetic trees differed from species relationships defined by standard taxonomy of Equidae and from trees based on mitochondrial or neutral gene sequence data. Analysis of selection showed differences between the less variable DRA and more variable DQA genes. DRA alleles were more often shared by more species. The DQA sequences analysed showed strong amongst-species positive selection; the selected amino acid positions mostly corresponded to selected positions in rodent and human DQA genes.     

Rapid Evolution by Positive Selection and Gene Gain and Loss: PLA<Subscript>2</Subscript> Venom Genes in Closely Related <Emphasis Type="Italic">Sistrurus</Emphasis> Rattlesnakes with Divergent Diets

Rapid evolution of snake venom genes by positive selection has been reported previously but key features of this process such as the targets of selection, rates of gene turnover, and functional diversity of toxins generated remain unclear. This is especially true for closely related species with divergent diets. We describe the evolution of PLA₂ gene sequences isolated from genomic DNA from four taxa of Sistrurus rattlesnakes which feed on different prey. We identified four to seven distinct PLA₂ sequences in each taxon and phylogenetic analyses suggest that these sequences represent a rapidly evolving gene family consisting of both paralogous and homologous loci with high rates of gene gain and loss. Strong positive selection was implicated as a driving force in the evolution of these protein coding sequences. Exons coding for amino acids that make up mature proteins have levels of variation two to three times greater than those of the surrounding noncoding intronic sequences. Maximum likelihood models of coding sequence evolution reveal that a high proportion (∼30%) of all codons in the mature protein fall into a class of codons with an estimated d _N /d _S (ω) ratio of at least 2.8. An analysis of selection on individual codons identified nine residues as being under strong (p < 0.01) positive selection, with a disproportionately high proportion of these residues found in two functional regions of the PLA₂ protein (surface residues and putative anticoagulant region). This is direct evidence that diversifying selection has led to high levels of functional diversity due to structural differences in proteins among these snakes. Overall, our results demonstrate that both gene gain and loss and protein sequence evolution via positive selection are important evolutionary forces driving adaptive divergence in venom proteins in closely related species of venomous snakes.     

Polymorphism and signatures of selection in the multimammate rat DQB gene

In order to test if DQB is a good candidate marker to investigate the relationship between major histocompatibility complex genes and pathogens in natural populations of Mastomys natalensis, we assessed the polymorphism and evolutionary history of this gene. Twenty-four individuals were genotyped for exon 2 of DQB using capillary electrophoresis single-strand conformation polymorphism, cloning, and sequencing. We found 21 different alleles. Four individuals show three alleles implying a duplication event in the history of this gene. Each distinct sequence translates to give a distinct amino acid sequence and there are strong signals of positive selection on peptide binding sites. Signals of recombination were found in the sequences suggesting that recombination has played a role in generating allelic diversity. Although trans-taxon polymorphism is present at the interspecific level in DQB exon 2 sequences of Mus species, we did not find any evidence of allele sharing among Muridae genera. This indicates a temporal limit of DQB allele sharing in Muridae of less than 8 Mya. In conclusion, although DQB seems to be a good marker to investigate pathogen-driven selection, the polymorphism of gene copy number may restrict its utility in natural populations.     

The comparative amino acid sequences, substrate specificities and gene or cDNA nucleotide sequences of some prokaryote and eukaryote amidinotransferases: implications for evolution

The amino acid sequences of the amidinotransferases and the nucleotide sequences of their genes or cDNA from four Streptomyces species (seven genes) and from the kidneys of rat, pig, human and human pancreas were compared. The overall amino acid and nucleotide sequences of the prokaryotes and eukaryotes were very similar and further, three regions were identified that were highly identical. Evidence is presented that there is virtually zero chance that the overall and high identity regions of the amino acid sequence similarities and the overall nucleotide sequence similarities between Streptomyces and mammals represent random match. Both rat and lamprey amidinotransferases were able to use inosamine phosphate, the amidine group acceptor of Streptomyces. We have concluded that the structure and function of the amidinotransferases and their genes has been highly conserved through evolution from prokaryotes to eukaryotes. The evolution has occurred with: (1) a high degree of retention of nucleotide and amino acid sequences; (2) a high degree of retention of the primitive Streptomyces guanine+cytosine (G+C) third codon position composition in certain high identity regions of the eukaryote cDNA; (3) a decrease in the specificities for the amidine group acceptors; and (4) most of the mutations silent in the regions suggested to code for active sites in the enzymes.     

Isolation and diversity analysis of resistance gene analogues (RGAs) from cultivated and wild strawberries

Degenerate oligonucleotide primers, designed based on conserved regions of Nucleotide Binding Site (NBS) domains from previously cloned plant resistance genes, were used to isolate Resistance Gene Analogues (RGAs) from wild and cultivated strawberries. Seven distinct families of RGAs of the NBS-LRR type were identified from two related wild species, Fragaria vesca and F. chiloensis, and six different Fragaria × ananassa cultivars. With one exception (GAV-3), the deduced amino acid sequences of strawberry RGAs showed strong similarity to TIR (Toll Interleukin I Receptor)-type R genes from Arabidopsis, tobacco and flax, suggesting the existence of common ancestors. GAV-3 seemed to be more closely related to the non-TIR type. Further studies showed that the recombination level and the ratio of non-synonymous to synonymous substitutions within families were low. These data suggest that NBS-encoding sequences of RGAs in strawberry are subject to a gradual accumulation of mutations leading to purifying selection, rather than to a diversifying process. The present paper is the first report on RGAs in strawberry.Communicated by M.-A. Grandbastien     

Distinct Evolutionary Patterns Between Two Duplicated Color Vision Genes Within Cyprinid Fishes

We investigated the molecular evolution of duplicated color vision genes (LWS-1 and SWS2) within cyprinid fish, focusing on the most cavefish-rich genus—Sinocyclocheilus. Maximum likelihood-based codon substitution approaches were used to analyze the evolution of vision genes. We found that the duplicated color vision genes had unequal evolutionary rates, which may lead to a related function divergence. Divergence of LWS-1 was strongly influenced by positive selection causing an accelerated rate of substitution in the proportion of pocket-forming residues. The SWS2 pigment experienced divergent selection between lineages, and no positively selected site was found. A duplicate copy of LWS-1 of some cyprinine species had become a pseudogene, but all SWS2 sequences remained intact in the regions examined in the cyprinid fishes examined in this study. The pseudogenization events did not occur randomly in the two copies of LWS-1 within Sinocyclocheilus species. Some cave species of Sinocyclocheilus with numerous morphological specializations that seem to be highly adapted for caves, retain both intact copies of color vision genes in their genome. We found some novel amino acid substitutions at key sites, which might represent interesting target sites for future mutagenesis experiments. Our data add to the increasing evidence that duplicate genes experience lower selective constraints and in some cases positive selection following gene duplication. Some of these observations are unexpected and may provide insights into the effect of caves on the evolution of color vision genes in fishes.     

Evolutionary history of an MHC gene in two leporid species: characterisation of <Emphasis Type="Italic">Mhc-DQA</Emphasis> in the European brown hare and comparison with the European rabbit

We surveyed the genetic diversity of the expressed major histocompatibility complex class II DQA locus in natural populations of European brown hares, Lepus europaeus, from Austria and Belgium (267 individuals in total). Based on cDNA sequences, we designed hare-specific primers to amplify the highly variable second exon of the DQA gene. Using cloning–sequencing methodology and capillary electrophoresis single-strand conformation polymorphism, we found ten alleles of the DQA exon 2 locus across these two European regions, of which eight are described for the first time. To search for signals of selection and recombination in the evolution of the DQA gene within the leporids, we augmented our sample with orthologous DQA alleles from the European rabbit, Oryctolagus cuniculus, in order to carry out a species level, species pairwise comparison. We found evidence of recombination in the history of the DQA sequences in leporids with some recombinant alleles bridging the species divide. In both species, selection on peptide binding site codons can be detected, though stronger for the rabbit. This result suggests that there may be a differential selection pressure in the deeper evolutionary history of these two species due to differences in several demographic and ecological traits likely subjecting them to differential selection by parasites. Finally, evolutionary relationships show a widespread and statistically significant intermingling of alleles from the two species. The many macroparasites shared between hares and rabbits may explain this pattern of trans-species polymorphism. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users. The nucleotide sequence data reported in this paper have been submitted to Genbank and have been assigned the accession numbers FJ225335–FJ225346.     

Local Patterns of Nucleotide Polymorphism Are Highly Variable in the Selfing Species Arabidopsis thaliana

Neighboring genes predictably share similar evolutionary histories to an extent delineated by recombination. This correlation should extend across multiple linked genes in a selfing species such as Arabidopsis thaliana due to its low effective recombination rate. To test this prediction, we performed a molecular population genetics analysis of nucleotide polymorphism and divergence in chromosomal regions surrounding four low-diversity loci. Three of these loci, At1g67140, At3g03700, and TERMINAL FLOWER1 (TFL1), have been previously implicated as targets of selection and we would predict stronger correlations in polymorphism between neighboring loci due to genetic hitchhiking around these loci. The remaining locus, At1g04300, was identified in a study of linkage disequilibrium surrounding the CRYPTOCHROME2 (CRY2) locus. Although we found broad valleys of reduced nucleotide variation around two of our focal genes, At1g67140 and At3g03700, all chromosomal regions exhibited extreme variation in the patterns of polymorphism and evolution between neighboring loci. Although three of our four regions contained potential targets of selection, application of the composite-likelihood-ratio test of selection in conjunction with a goodness-of-fit test supports the selection hypothesis only for the region containing At3g03700. The degree of discordance in evolutionary histories between linked loci within each region generally correlated with estimates of recombination and linkage disequilibrium for that region, with the exception of the region containing At1g04300. We discuss the implications of these data for future population genetics analyses and genomics studies in A. thaliana. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

棉花DUR3基因的鉴定及进化分析

植物DUR3同源蛋白属于钠离子/溶质共运蛋白家族的尿素高亲和力运输蛋白，在植物体对外源尿素的主动吸收及内源尿素的再分配过程中具有重要作用。为明确棉花DUR3基因的结构和进化情况，基于生物信息学的方法，从全基因组水平鉴定陆地棉和雷蒙德氏棉的DUR3基因，并对基因结构、跨膜结构域、基序分布、进化关系等进行分析。结果表明：(1)从陆地棉A亚组和D亚组染色体各鉴定出1个DUR3基因，从雷蒙德氏棉基因组鉴定出1个DUR3基因。这3个棉花DUR3同源蛋白同其他植物DUR3同源蛋白一样，具有15个跨膜结构域，具有3个位置一致、高度保守的基序。(2)基因结构分析表明，双子叶植物DUR3基因的外显子个数明显多于单子叶植物，这3个棉花DUR3基因的外显子个数亦是如此。(3)根据物种间种属亲缘关系，对不同物种DUR3氨基酸序列构建的进化树显示，棉花的同双子叶植物的聚在一起。(4)DUR3直系同源基因和旁系同源基因的K_a/K_s比值普遍均大于1,说明这些基因在进化过程中主要受到正向选择的作用。该研究结果为深入研究棉花DUR3同源蛋白提供了理论基础。     

The sequence and organization of the core histone H3 and H4 genes in the early branching amitochondriate protistTrichomonas vaginalis

Among the unicellular protists, several of which are parasitic, some of the most divergent eukaryotic species are found. The evolutionary distances between protists are so large that even slowly evolving proteins like histones are strongly divergent. In this study we isolated cDNA and genomic histone H3 and H4 clones fromTrichomonas vaginalis. Two histone H3 and three histone H4 genes were detected on three genomic clones with one complete H3 and two complete H4 sequences. H3 and H4 genes were divergently transcribed with very short intergenic regions of only 194 bp, which containedT. vaginalis-specific as well as histone-specific putative promoter elements. Southern blot analysis showed that there may be several more histone gene pairs. The two complete histone H4 genes were different on the nucleotide level but encoded the same amino acid sequence. Comparison of the amino acid sequences of theT. vaginalis H3 and H4 histones with sequences from animals, fungi, and plants as well as other protists revealed a significant divergence not only from the sequences in multicellular organisms but especially from the sequences in other protists likeEntamoeba histolytica, Trypanosoma cruzi, andLeishmania infantum.     

Molecular evolution of paclitaxel biosynthetic genes <Emphasis Type="Italic">TS</Emphasis> and <Emphasis Type="Italic">DBAT</Emphasis> of <Emphasis Type="Italic">Taxus</Emphasis> species

Evolutionary patterns of sequence divergence were analyzed in genes from the conifer genus Taxus (yew), encoding paclitaxel biosynthetic enzymes taxadiene synthase (TS) and 10-deacetylbaccatin III-10β-O-acetyltransferase (DBAT). N-terminal fragments of TS, full-length DBAT and internal transcribed spacer (ITS) were amplified from 15 closely related Taxus species and sequenced. Premature stop codons were not found in TS and DBAT sequences. Codon usage bias was not found, suggesting that synonymous mutations are selectively neutral. TS and DBAT gene trees are not consistent with the ITS tree, where species formed monophyletic clades. In fact, for both genes, alleles were sometimes shared across species and parallel amino acid substitutions were identified. While both TS and DBAT are, overall, under purifying selection, we identified a number of amino acids of TS under positive selection based on inference using maximum likelihood models. Positively selected amino acids in the N-terminal region of TS suggest that this region might be more important for enzyme function than previously thought. Moreover, we identify lineages with significantly elevated rates of amino acid substitution using a genetic algorithm. These findings demonstrate that the pattern of adaptive paclitaxel biosynthetic enzyme evolution can be documented between closely related Taxus species, where species-specific taxane metabolism has evolved recently.     

Positive Selection in Prion Protein

The prion diseases, such as Creutzfeldt-Jakob disease of humans and bovine spongiform encephalopathy, involve the aberrant metabolism and accumulation of prion protein PrP. There are three contradictory hypotheses about evolution of prion protein gene PRNP. Population genetic studies have proposed that PRNP could be under balancing selection, strong purifying selection, or mainly positive selection. We made use of the maximum likelihood tests for detection of positive selection at the amino acid level and present availability of PRNP coding sequences to contribute to these disagreements. Positive selection could occur at amino acids residing in active sites, and at amino acids involved in protein-protein interactions. Thus we tested a hypothesis that positive selection at the amino acid level in PrP might have taken place in human and related species from the superordinal group Euarchonta, as well as in bovine and related species from the superordinal clade Laurasiatheria. Our study and the present experimental evidences indicate that positive selection at the amino acid level might have taken place in the PrP signal sequences and conformationally plastic PrP regions, as well as at the protein X binding sites. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users. Prof. Vera Gamulin passed away.     

The Chlamydomonas chloroplast clpP gene contains translated large insertion sequences and is essential for cell growth

Sequence determination of the chloroplast clpP gene from two distantly related Chlamydomonas species (C. reinhardtii and C. eugametos) revealed the presence of translated large insertion sequences (IS1 and IS2) that divide the clpP gene into two or three sequence domains (SDs) and are not found in homologous genes in other organisms. These insertion sequences do not resemble RNA introns, and are not spliced out at the mRNA level. Instead, each insertion sequence forms a continuous open reading frame with its upstream and downstream sequence domains. IS1 specifies a potential polypeptide sequence of 286 and 318 amino acid residues in C. reinhardtii and C. eugametos, respectively. IS2 encodes a 456 amino acid polypeptide and is present only in C. eugametos. The two Chlamydomonas IS1 sequences show substantial similarity; however, there is no significant sequence similarity either between IS1 and IS2 or between these insertion sequences and any other known protein coding sequences. The C. reinhardtii clpP gene was further shown to be essential for cell growth, as demonstrated through targeted gene disruption by particle gun-mediated chloroplast transformation. Only heteroplasmic transformants could be obtained, even under mixotrophic growth conditions. The heteroplasmic transformants were stable only under selection pressure for the disrupted clpP, rapidly segregated into wild-type cells when the selection pressure was removed, and grew significantly more slowly than wildtype cells under phototrophic conditions.     

Evidence for Recombination in the Microcystin Synthetase (<Emphasis Type="Italic">mcy</Emphasis>) Genes ofToxic Cyanobacteria <Emphasis Type="Italic">Microcystis</Emphasis><Emphasis Type="Bold">spp</Emphasis>

Recombination has been suggested to be an important factor for the genetic variation of bacterial genes, but few studies have dealt with intragenic recombination between the same or closely related species of cyanobacteria. Here we provide strong evidence for recombination in the microcystin synthetase (mcy) gene cluster of the toxic cyanobacteria Microcystis spp. This gene cluster contains 10 genes (mcyA to J) that encode a mixed polyketide synthase (PKS)/nonribosomal peptide synthetase (NRPS) complex. mcy gene sequences were determined for four selected regions (within mcyA, D, G, and J) within the mcy gene cluster from 1 Canadian and 10 Asian toxic Microcystis and compared with previously published mcy sequences. Split decomposition analysis indicated a reticulate phylogeny of mcyA, and several potential recombination tracts of mcyA were identified by the RDP analysis and a runs test implemented in GENECONV. In contrast, no recombination was detected in the mcyD, G, and J sequences. However, discrepancies among the four mcy gene genealogies were evident from the results of independent split decomposition analyses, which were further supported by incongruence length difference (ILD) tests. Taken together, these findings suggest that both intragenic and intergenic recombination within the mcy gene cluster contributes to the genetic diversity of the mcy genes of Microcystis spp.This article contains online supplementary material.     

Natural selection in the <Emphasis Type="Italic">TLR-</Emphasis>related genes in the course of primate evolution

The innate immune system constitutes the front line of host defense against pathogens. Toll-like receptors (TLRs) recognize molecules derived from pathogens and play crucial roles in the innate immune system. Here, we provide evidence that the TLR-related genes have come under natural selection pressure in the course of primate evolution. We compared the nucleotide sequences of 16 TLR-related genes, including TLRs (TLR1–10), MYD88, TILAP, TICAM1, TICAM2, MD2, and CD14, among seven primate species. Analysis of the non-synonymous/synonymous substitution ratio revealed the presence of both strictly conserved and rapidly evolving regions in the TLR-related genes. The genomic segments encoding the intracellular Toll/interleukin 1 receptor domains, which exhibited lower rates of non-synonymous substitution, have undergone purifying selection. In contrast, TLR4, which carried a high proportion of non-synonymous substitutions in the part of extracellular domain spanning 200 amino acids, was found to have been the suggestive target of positive Darwinian selection in primate evolution. However, sequence analyses from 25 primate species, including eight hominoids, six Old World monkeys, eight New World monkeys, and three prosimians, showed no evidence that the pressure of positive Darwinian selection has shaped the pattern of sequence variations in TLR4 among New World monkeys and prosimians. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Variation in recombination rates across Vitis species

Recombination rate data are presented for three populations of grape based on framework genetic linkage maps developed with simple-sequence repeat markers. These linkage maps were constructed from different Vitis species and represent three genetic backgrounds. The first population is pure Vitis vinifera, derived from a cross of the European cultivars Riesling and Cabernet Sauvignon. The second is an interspecific cross between two commercially used rootstock cultivars of different North American Vitis species parentage, Ramsey (Vitis champinii) and Riparia Gloire (Vitis riparia). The third population, D8909-15 (Vitis rupestris × (Vitis arizonica/Vitis girdiana form)) × F8909-17 (V. rupestris × (V. arizonica/Vitis candicans form)), is an F1 from two half-sibs. Genome-wide and chromosome-wide recombination rates varied across the three populations and among the six Vitis parents. Global recombination rates in the parents of the third F1 population, with a complex Vitis background, were significantly reduced. In the first and third populations, the recombination rate was significantly greater in the male parent. Specific genome locations with frequent heterogeneity in recombination were identified, suggesting that recombination rates are not equal across the Vitis genome. The identification of regions with suppressed or high recombination will aid grape breeders and geneticists who rely on recombination events to introgress disease resistance genes from the genomes of wild Vitis species, develop fine-scale genetic maps, and clone disease resistance genes. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Genome wide evolutionary analyses reveal serotype specific patterns of positive selection in selected Salmonella serotypes

Background  

The bacterium Salmonella enterica includes a diversity of serotypes that cause disease in humans and different animal species. Some Salmonella serotypes show a broad host range, some are host restricted and exclusively associated with one particular host, and some are associated with one particular host species, but able to cause disease in other host species and are thus considered "host adapted". Five Salmonella genome sequences, representing a broad host range serotype (Typhimurium), two host restricted serotypes (Typhi [two genomes] and Paratyphi) and one host adapted serotype (Choleraesuis) were used to identify core genome genes that show evidence for recombination and positive selection.     

Inter- and Intralocus Recombination Drive MHC Class IIB Gene Diversification in a Teleost, the Three-Spined Stickleback Gasterosteus aculeatus

The mutational mechanism underlying the striking diversity in MHC (major histocompatibility complex) genes in vertebrates is still controversial. In order to evaluate the role of inter- and intragenic recombination in MHC gene diversification, we examined patterns of nucleotide polymorphism across an exon/intron boundary in a sample of 31 MHC class IIB sequences of three-spined stickleback (Gasterosteus aculeatus). MHC class IIB genes of G. aculeatus were previously shown to be under diversifying (positive) selection in mate choice and pathogen selection experiments. Based on recoding of alignment gaps, complete intron 2 sequences were grouped into three clusters using maximum-parsimony analysis. Two of these groups had >90% bootstrap support and were tentatively assigned single locus status. Intron nucleotide diversity within and among loci was low (p-distance within and among groups = 0.016 and 0.019, respectively) and fourfold lower than the rate of silent mutations in exon 2, suggesting that noncoding regions are homogenized by frequent interlocus recombination. A substitution analysis using GENECONV revealed as many intergenic conversion events as intragenic ones. Recombination between loci may explain the occurrence of sequence variants that are particularly divergent, as is the case in three-spined stickleback, with nucleotide diversity attaining d_N = 0.39 (peptide-binding residues only). For both MHC class II loci we also estimated the amount of intragenic recombination as population rate (4N_er) under the coalescent and found it to be approximately three times higher compared to point mutations (Watterson estimate per gene, 4N_eμ). Nonindependence of molecular evolution across loci and frequent recombination suggest that MHC class II genes of bony fish may follow different evolutionary dynamics than those of mammals. Our finding of widespread recombination suggests that phylogenies of MHC genes should not be based on coding segments but rather on noncoding introns. [Reviewing Editor: Dr. Richard Kliman]     

Molecular population genetic analysis of the streptokinase gene of Streptococcus pyogenes: mosaic alleles generated by recombination

To understand the mechanisms governing molecular evolution of the streptokinase gene (skn), a 384 bp DNA fragment encoding two variable regions of the molecule was characterized in 47 isolates of Streptococcus pyogenes. The results reveal that alleles of the streptokinase gene have a mosaic structure, and provide strong evidence for intragenic recombination. Moreover, organisms that are well differentiated in overall chromosomal character have identical skn alleles, which suggests that horizontal gene transfer and recombination have participated in the evolution of this locus. No simple relationship between skn allele and serum opacity factor production or specific disease was identified. The predicted amino acid sequences of highly divergent skn alleles are strikingly similar in hydrophobicity and hydrophobicity profiles, distribution of amphipathic and flexible regions, surface probability plots, and antigenic indices, indicating that despite extensive nucleotide polymorphism in the two skn variable regions, selective pressure has constrained overall structural divergence. These results add to an important emerging theme that intragenic recombination plays a critical role in diversifying genes coding for streptococcal virulence factors.