Elevation as a selective force on mitochondrial respiratory chain complexes of the Phrynocephalus lizards in the Tibetan plateau

Animals living in extremely high elevations have to adapt to low temperatures and low oxygen availability (hypoxia), but the underlying genetic mechanisms associated with these adaptations are still unclear. The mitochondrial respiratory chain can provide >95% of the ATP in animal cells, and its efficiency is influenced by temperature and oxygen availability. Therefore, the respiratory chain complexes (RCCs) could be important molecular targets for positive selection associated with respiratory adaptation in high-altitude environments. Here, we investigated positive selection in 5 RCCs and their assembly factors by analyzing sequences of 106 genes obtained through RNA-seq of all 15 Chinese Phrynocephalus lizard species, which are distributed from lowlands to the Tibetan plateau (average elevation >4,500 m). Our results indicate that evidence of positive selection on RCC genes is not significantly different from assembly factors, and we found no difference in selective pressures among the 5 complexes. We specifically looked for positive selection in lineages where changes in habitat elevation happened. The group of lineages evolving from low to high altitude show stronger signals of positive selection than lineages evolving from high to low elevations. Lineages evolving from low to high elevation also have more shared codons under positive selection, though the changes are not equivalent at the amino acid level. This study advances our understanding of the genetic basis of animal respiratory metabolism evolution in extreme high environments and provides candidate genes for further confirmation with functional analyses.     

Both temperature fluctuations and East Asian monsoons have driven plant diversification in the karst ecosystems from southern China

Karst ecosystems in southern China are species‐rich and have high levels of endemism, yet little is known regarding the evolutionary processes responsible for the origin and diversification of karst biodiversity. The genus Primulina (Gesneriaceae) comprises ca. 170 species endemic to southern China with high levels of ecological (edaphic) specialization, providing an exceptional model to study the plant diversification in karsts. We used molecular data from nine chloroplast and 11 nuclear regions and macroevolutionary analyses to assess the origin and cause of species diversification due to palaeoenvironmental changes and edaphic specialization in Primulina. We found that speciation was positively associated with changes in past temperatures and East Asian monsoons through the evolutionary history of Primulina. Climatic change around the mid‐Miocene triggered an early burst followed by a slowdown of diversification rate towards the present with the climate cooling. We detected different speciation rates among edaphic types, and transitions among soil types were infrequently and did not impact the overall speciation rate. Our findings suggest that both global temperature changes and East Asian monsoons have played crucial roles in floristic diversification within the karst ecosystems in southern China, such that speciation was higher when climate was warmer and wetter. This is the first study to directly demonstrate that past monsoon activity is positively correlated with speciation rate in East Asia. This case study could motivate further investigations to assess the impacts of past environmental changes on the origin and diversification of biodiversity in global karst ecosystems, most of which are under threat.     

Adaptive Evolution in the GAF Domain of Phytochromes in Gymnosperms

The GAF domain of phytochrome is essential for photoconversion and signal transduction. In gymnosperms, it exists in all members of the phytochrome family that experience gene duplication. Maximum-likelihood models of codon substitution can provide a framework for constructing likelihood ratio tests of changes in selective pressure and make clear predictions about patterns of genetic change following gene duplication. In this study, 68 gymnosperm GAF sequences were analyzed to identify lineages and sites under positive selection. Our results indicate that (1) positive selection at a few sites (3.6%), rather than relaxation of selective constraints, has played a major role in the evolution of the gymnosperm GAF domain; (2) strong positive selective pressure tends to occur in the recent PHYP lineages of cogeneric species, but is absent in old lineages consisting of distantly related species; and (3) the selective pressure indicated by the ω ratio varies greatly among lineages and sites in the GAF domain.     

ASPM and the Evolution of Cerebral Cortical Size in a Community of New World Monkeys

The ASPM (abnormal spindle-like microcephaly associated) gene has been proposed as a major determinant of cerebral cortical size among primates, including humans. Yet the specific functions of ASPM and its connection to human intelligence remain controversial. This debate is limited in part by a taxonomic focus on Old World monkeys and apes. Here we expand the comparative context of ASPM sequence analyses with a study of New World monkeys, a radiation of primates in which enlarged brain size has evolved in parallel in spider monkeys (genus Ateles) and capuchins (genus Cebus). The primate community of Costa Rica is perhaps a model system because it allows for independent pairwise comparisons of smaller- and larger-brained species within two taxonomic families. Accordingly, we analyzed the complete sequence of exon 18 of ASPM in Ateles geoffroyi, Alouatta palliata, Cebus capucinus, and Saimiri oerstedii. As the analysis of multiple species in a genus improves phylogenetic reconstruction, we also analyzed eleven published sequences from other New World monkeys. Our exon-wide, lineage-specific analysis of eleven genera and the ratio of rates of nonsynonymous to synonymous substitutions (d_N/d_S) on ASPM revealed no detectable evidence for positive selection in the lineages leading to Ateles or Cebus, as indicated by d_N/d_S ratios of <1.0 (0.6502 and 0.4268, respectively). Our results suggest that a multitude of interacting genes have driven the evolution of larger brains among primates, with different genes involved in this process in different encephalized lineages, or at least with evidence for positive selection not readily apparent for the same genes in all lineages. The primate community of Costa Rica may serve as a model system for future studies that aim to elucidate the molecular mechanisms underlying cognitive capacity and cortical size.     

Gene flow,linked selection,and divergent sorting of ancient polymorphism shape genomic divergence landscape in a group of edaphic specialists

Interpreting the formation of genomic variation landscape, especially genomic regions with elevated differentiation (i.e. islands), is fundamental to a better understanding of the genomic consequences of adaptation and speciation. Edaphic islands provide excellent systems for understanding the interplay of gene flow and selection in driving population divergence and speciation. However, discerning the relative contribution of these factors that modify patterns of genomic variation remains difficult. We analysed 132 genomes from five recently divergent species in Primulina genus, with four species distributed in Karst limestone habitats and the fifth one growing in Danxia habitats. We demonstrated that both gene flow and linked selection have contributed to genome-wide variation landscape, where genomic regions with elevated differentiation (i.e., islands) were largely derived by divergent sorting of ancient polymorphism. Specifically, we identified several lineage-specific genomic islands that might have facilitated adaptation of P. suichuanensis to Danxia habitats. Our study is amongst the first cases disentangling evolutionary processes that shape genomic variation of plant specialists, and demonstrates the important role of ancient polymorphism in the formation of genomic islands that potentially mediate adaptation and speciation of endemic plants in special soil habitats.     

Whole-gene positive selection, elevated synonymous substitution rates, duplication, and indel evolution of the chloroplast clpP1 gene

Background

Synonymous DNA substitution rates in the plant chloroplast genome are generally relatively slow and lineage dependent. Non-synonymous rates are usually even slower due to purifying selection acting on the genes. Positive selection is expected to speed up non-synonymous substitution rates, whereas synonymous rates are expected to be unaffected. Until recently, positive selection has seldom been observed in chloroplast genes, and large-scale structural rearrangements leading to gene duplications are hitherto supposed to be rare.

Methodology/Principle Findings

We found high substitution rates in the exons of the plastid clpP1 gene in Oenothera (the Evening Primrose family) and three separate lineages in the tribe Sileneae (Caryophyllaceae, the Carnation family). Introns have been lost in some of the lineages, but where present, the intron sequences have substitution rates similar to those found in other introns of their genomes. The elevated substitution rates of clpP1 are associated with statistically significant whole-gene positive selection in three branches of the phylogeny. In two of the lineages we found multiple copies of the gene. Neighboring genes present in the duplicated fragments do not show signs of elevated substitution rates or positive selection. Although non-synonymous substitutions account for most of the increase in substitution rates, synonymous rates are also markedly elevated in some lineages. Whereas plant clpP1 genes experiencing negative (purifying) selection are characterized by having very conserved lengths, genes under positive selection often have large insertions of more or less repetitive amino acid sequence motifs.

Conclusions/Significance

We found positive selection of the clpP1 gene in various plant lineages to correlated with repeated duplication of the clpP1 gene and surrounding regions, repetitive amino acid sequences, and increase in synonymous substitution rates. The present study sheds light on the controversial issue of whether negative or positive selection is to be expected after gene duplications by providing evidence for the latter alternative. The observed increase in synonymous substitution rates in some of the lineages indicates that the detection of positive selection may be obscured under such circumstances. Future studies are required to explore the functional significance of the large inserted repeated amino acid motifs, as well as the possibility that synonymous substitution rates may be affected by positive selection.     

Molecular evolution of a family of resistance gene analogs of nucleotide-binding site sequences in <Emphasis Type="Italic">Solanum lycopersicum</Emphasis>

Nucleotide-binding site-leucine-rich repeats (NBS-LRR) gene families are one of the major plant resistance genes. Genomic NBS evolution was studied in many plant species for diverse arrays of NBS gene families. In this study, we focused on one family of NBS sequences in an attempt to understand how closely related NBS sequences evolved in the light of selection in domesticated plant species. A phylogenetic analysis revealed five major clades (A–E) and five subclades (A1–A5) within clade A of cloned NBS sequences. Positive selection was only detected in newly evolved NBS lineages in subclades of clade A. Positively selected codon sites were found among NBS sequences of clade A. A sliding-window analysis revealed that regions with Ka/Ks ratios of >1 were in the inter-motifs when paired clades were compared, but regions with Ka/Ks ratios of >1 were found across NBS sequences when subclades of clade A were compared. Our results based on a family of closely related NBS sequences showed that positive selection was first exerted on specific lineages across all NBS sequences after selective constraints. Subsequently, sequences with mutations in commonly conserved motifs were scrutinized by purifying selection. In the long term, conserved high frequency alleles in commonly conserved motifs and changes in inter-motifs were maintained in the investigated family of NBS sequences. Moreover, codons identified to be under positive selection in the inter-motifs were mainly located in regions involved in functions of ATP binding or hydrolysis.     

The evolutionary rate of leaf osmotic strength drives diversification of Primulina species in karst regions

Leaf water storage capacity and osmotic strength are important traits enabling species to adapt to environments that are often moisture limited. However, whether these drought tolerance traits are correlated with the species diversification rate (DR) of plant lineages is yet to be determined. In this study, we selected a species-rich genus (Primulina) of plants widely distributed in karst regions in which species frequently experience variable periods of drought. We measured water storage capacity-related traits (including leaf thickness and water content per mass) and saturated osmotic strength in the leaves of 58 Primulina species growing in a common garden. Subsequently, using phylogenetic methods, we examined the relationships between the rate of species diversification and the drought tolerance traits and between the species DR and evolutionary rates of these traits. We found that neither water storage capacity nor saturated osmotic strength showed significant relationships with the rate of species diversification. However, the evolutionary rate of saturated osmotic strength showed a significant correlation with the species DR, although no comparable significant relationship was detected regarding the evolutionary rate of water storage capacity. Our study indicates that the diversification among Primulina species has typically been accompanied by an extensive divergence of leaf osmotic strength but not a divergence in leaf water storage capacity. These findings will enhance our current understanding of how drought tolerance influences the diversification of plant species in karst regions.     

Phylogenetic Analysis of Burkholderia Species by Multilocus Sequence Analysis

Burkholderia comprises more than 60 species of environmental, clinical, and agro-biotechnological relevance. Previous phylogenetic analyses of 16S rRNA, recA, gyrB, rpoB, and acdS gene sequences as well as genome sequence comparisons of different Burkholderia species have revealed two major species clusters. In this study, we undertook a multilocus sequence analysis of 77 type and reference strains of Burkholderia using atpD, gltB, lepA, and recA genes in combination with the 16S rRNA gene sequence and employed maximum likelihood and neighbor-joining criteria to test this further. The phylogenetic analysis revealed, with high supporting values, distinct lineages within the genus Burkholderia. The two large groups were named A and B, whereas the B. rhizoxinica/B. endofungorum, and B. andropogonis groups consisted of two and one species, respectively. The group A encompasses several plant-associated and saprophytic bacterial species. The group B comprises the B. cepacia complex (opportunistic human pathogens), the B. pseudomallei subgroup, which includes both human and animal pathogens, and an assemblage of plant pathogenic species. The distinct lineages present in Burkholderia suggest that each group might represent a different genus. However, it will be necessary to analyze the full set of Burkholderia species and explore whether enough phenotypic features exist among the different clusters to propose that these groups should be considered separate genera.     

Molecular Adaptation of rbcL in the Heterophyllous Aquatic Plant Potamogeton

Background

Heterophyllous aquatic plants show marked phenotypic plasticity. They adapt to environmental changes by producing different leaf types: submerged, floating and terrestrial leaves. By contrast, homophyllous plants produce only submerged leaves and grow entirely underwater. Heterophylly and submerged homophylly evolved under selective pressure modifying the species-specific optima for photosynthesis, but little is known about the evolutionary outcome of habit. Recent evolutionary analyses suggested that rbcL, a chloroplast gene that encodes a catalytic subunit of RuBisCO, evolves under positive selection in most land plant lineages. To examine the adaptive evolutionary process linked to heterophylly or homophylly, we analyzed positive selection in the rbcL sequences of ecologically diverse aquatic plants, Japanese Potamogeton.

Principal Findings

Phylogenetic and maximum likelihood analyses of codon substitution models indicated that Potamogeton rbcL has evolved under positive Darwinian selection. The positive selection has operated specifically in heterophyllous lineages but not in homophyllous ones in the branch-site models. This suggests that the selective pressure on this chloroplast gene was higher for heterophyllous lineages than for homophyllous lineages. The replacement of 12 amino acids occurred at structurally important sites in the quaternary structure of RbcL, two of which (residue 225 and 281) were identified as potentially under positive selection.

Conclusions/Significance

Our analysis did not show an exact relationship between the amino acid replacements and heterophylly or homophylly but revealed that lineage-specific positive selection acted on the Potamogeton rbcL. The contrasting ecological conditions between heterophyllous and homophyllous plants have imposed different selective pressures on the photosynthetic system. The increased amino acid replacement in RbcL may reflect the continuous fine-tuning of RuBisCO under varying ecological conditions.     

Rapid Evolution, Genetic Variations, and Functional Association of the Human Spermatogenesis-Related Gene NYD-SP12

NYD-SP12 is a recently identified spermatogenesis-related gene with a pivotal role in human testis development. In this study, we analyzed between-species divergence and within-species variation of NYD-SP12 in seven representative primate species, four worldwide human populations, and 124 human clinical subjects. Our results indicate that NYD-SP12 evolves rapidly in both the human and the chimpanzee lineages, which is likely caused by Darwinian positive selection and/or sexual selection. We observed significant interpopulation divergence among human populations, which might be due to the varied demographic histories. In the association analysis, we demonstrated significant frequency discrepancy of a synonymous sequence polymorphism among the clinical groups with different sperm traits. [Reviewing Editor: Dr. Manyuan Long]     

广西报春苣苔属一新种及其系统位置(英文)

报春苣苔是一类生长在石灰岩基质上的多年生草本植物,目前已知的种类160余种,主要分布在我国南部和越南北部的石灰岩地区。该属植物许多种类分布极其狭域,往往仅零星分布在某个或少数几个地点。近年来,该属植物许多新种类不断地被发现和报道,预示着该属植物可能有更加丰富的多样性。该研究报道了在广西博白县发现的该属植物一新种。从形态上判断,该新种属于广义报春苣苔属,其具有的一些特别的形态特征,如花冠下唇比上唇长得多,上唇内部具有两个近圆形斑点(其它多数种类具有条形斑点),而易于与本属其它物种区分。而且,基于叶绿体mat K基因序列重建的系统发育关系表明,该新种的确属于报春苣苔属且广义报春苣苔属是单系类群。该新种与钟冠报春(P.swinglei)和疏花报春(P.laxiflora)最为近缘,其mat K基因与后两者分别有5个和11个碱基差异,暗示这些物种间经历了较长时间的独立演化或mat K基因在这些物种中具有较高的碱基替换速率。另外,该新种只生长在丹霞岩石上,而其近缘种以及同属大多数其它种类都生长在石灰岩基质上,暗示该新种获得了特殊的适应能力以及与其它近缘种间存在生态隔离。该新种以及近年来本属许多其它新种类的发现,进一步增加了对该属植物多样性的认识,以及该属植物的演化与多样性分化之谜。这些物种所具有的高度区域或局域特有性和稀有性,以及形态上的特殊性和遗传上的独立性,进一步强化了对该属植物开展保护的重要性。     

A cornucopia of cryptic species - a DNA barcode analysis of the gobiid fish genus Trimma (Percomorpha,Gobiiformes)

A genetic analysis of partial mitochondrial 5’ cytochrome c oxidase I gene (DNA barcode) sequences of 473 specimens assigned to 52 morphological species (including four known, but not formally named, species) of the gobiid genus Trimma revealed the presence of 94 genetic lineages. Each lineage was separated by > 2% sequence divergence. Thus there were an additional 42 haplogroups recognizable as provisional candidate species given that a value of > 2% difference is typical of different species of fishes. Such a high degree of apparently different cryptic species is, in our experience, virtually unprecedented among vertebrates. These results have precipitated further morphological research in a few cases, which has uncovered subtle differences independently corroborating the genetic results. However, such efforts are limited by the dearth of traditional systematists available to undertake the necessary time-consuming, and highly detailed, morphological research. In some cases, the genetic results we present are consistent with, and confirm, minor taxonomic distinctions based on morphology and/or colour pattern. In other instances, what had been recognized as a single species consists of several genetic lineages - up to eight in, for example, what we have identified based on morphology as Trimma okinawae. The increase from 52 to 94 potential species in our sampling raises the predicted total number of species in this genus from about 110 to nearly 200 (versus the 73 valid described species currently recognized).     

Gene-specific signatures of elevated non-synonymous substitution rates correlate poorly across the Plasmodium genus

Background

Comparative genome analyses of parasites allow large scale investigation of selective pressures shaping their evolution. An acute limitation to such analysis of Plasmodium falciparum is that there is only very partial low-coverage genome sequence of the most closely related species, the chimpanzee parasite P. reichenowi. However, if orthologous genes have been under similar selective pressures throughout the Plasmodium genus then positive selection on the P. falciparum lineage might be predicted to some extent by analysis of other lineages.

Principal Findings

Here, three independent pairs of closely related species in different sub-generic clades (P. falciparum and P. reichenowi; P. vivax and P. knowlesi; P. yoelii and P. berghei) were compared for a set of 43 candidate ligand genes considered likely to be under positive directional selection and a set of 102 control genes for which there was no selective hypothesis. The ratios of non-synonymous to synonymous substitutions (dN/dS) were significantly elevated in the candidate ligand genes compared to control genes in each of the three clades. However, the rank order correlation of dN/dS ratios for individual candidate genes was very low, less than the correlation for the control genes.

Significance

The inability to predict positive selection on a gene in one lineage by identifying elevated dN/dS ratios in the orthologue within another lineage needs to be noted, as it reflects that adaptive mutations are generally rare events that lead to fixation in individual lineages. Thus it is essential to complete the genome sequences of particular species of phylogenetic importance, such as P. reichenowi.     

Ecological genomics in Xanthomonas: the nature of genetic adaptation with homologous recombination and host shifts

Background

Comparative genomics provides insights into the diversification of bacterial species. Bacterial speciation usually takes place with lasting homologous recombination, which not only acts as a cohering force between diverging lineages but brings advantageous alleles favored by natural selection, and results in ecologically distinct species, e.g., frequent host shift in Xanthomonas pathogenic to various plants.

Results

Using whole-genome sequences, we examined the genetic divergence in Xanthomonas campestris that infected Brassicaceae, and X. citri, pathogenic to a wider host range. Genetic differentiation between two incipient races of X. citri pv. mangiferaeindicae was attributable to a DNA fragment introduced by phages. In contrast to most portions of the genome that had nearly equivalent levels of genetic divergence between subspecies as a result of the accumulation of point mutations, 10% of the core genome involving with homologous recombination contributed to the diversification in Xanthomonas, as revealed by the correlation between homologous recombination and genomic divergence. Interestingly, 179 genes were under positive selection; 98 (54.7%) of these genes were involved in homologous recombination, indicating that foreign genetic fragments may have caused the adaptive diversification, especially in lineages with nutritional transitions. Homologous recombination may have provided genetic materials for the natural selection, and host shifts likely triggered ecological adaptation in Xanthomonas. To a certain extent, we observed positive selection nevertheless contributed to ecological divergence beyond host shifting.

Conclusion

Altogether, mediated with lasting gene flow, species formation in Xanthomonas was likely governed by natural selection that played a key role in helping the deviating populations to explore novel niches (hosts) or respond to environmental cues, subsequently triggering species diversification.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1369-8) contains supplementary material, which is available to authorized users.     

Degree of Functional Divergence in Duplicates Is Associated with Distinct Roles in Plant Evolution

Gene duplication is a major mechanism to create new genes. After gene duplication, some duplicated genes undergo functionalization, whereas others largely maintain redundant functions. Duplicated genes comprise various degrees of functional diversification in plants. However, the evolutionary fate of high and low diversified duplicates is unclear at genomic scale. To infer high and low diversified duplicates in Arabidopsis thaliana genome, we generated a prediction method for predicting whether a pair of duplicate genes was subjected to high or low diversification based on the phenotypes of knock-out mutants. Among 4,017 pairs of recently duplicated A. thaliana genes, 1,052 and 600 are high and low diversified duplicate pairs, respectively. The predictions were validated based on the phenotypes of generated knock-down transgenic plants. We determined that the high diversified duplicates resulting from tandem duplications tend to have lineage-specific functions, whereas the low diversified duplicates produced by whole-genome duplications are related to essential signaling pathways. To assess the evolutionary impact of high and low diversified duplicates in closely related species, we compared the retention rates and selection pressures on the orthologs of A. thaliana duplicates in two closely related species. Interestingly, high diversified duplicates resulting from tandem duplications tend to be retained in multiple lineages under positive selection. Low diversified duplicates by whole-genome duplications tend to be retained in multiple lineages under purifying selection. Taken together, the functional diversities determined by different duplication mechanisms had distinct effects on plant evolution.     

Dinoflagellate phylogeny as inferred from heat shock protein 90 and ribosomal gene sequences

Background

Interrelationships among dinoflagellates in molecular phylogenies are largely unresolved, especially in the deepest branches. Ribosomal DNA (rDNA) sequences provide phylogenetic signals only at the tips of the dinoflagellate tree. Two reasons for the poor resolution of deep dinoflagellate relationships using rDNA sequences are (1) most sites are relatively conserved and (2) there are different evolutionary rates among sites in different lineages. Therefore, alternative molecular markers are required to address the deeper phylogenetic relationships among dinoflagellates. Preliminary evidence indicates that the heat shock protein 90 gene (Hsp90) will provide an informative marker, mainly because this gene is relatively long and appears to have relatively uniform rates of evolution in different lineages.

Methodology/Principal Findings

We more than doubled the previous dataset of Hsp90 sequences from dinoflagellates by generating additional sequences from 17 different species, representing seven different orders. In order to concatenate the Hsp90 data with rDNA sequences, we supplemented the Hsp90 sequences with three new SSU rDNA sequences and five new LSU rDNA sequences. The new Hsp90 sequences were generated, in part, from four additional heterotrophic dinoflagellates and the type species for six different genera. Molecular phylogenetic analyses resulted in a paraphyletic assemblage near the base of the dinoflagellate tree consisting of only athecate species. However, Noctiluca was never part of this assemblage and branched in a position that was nested within other lineages of dinokaryotes. The phylogenetic trees inferred from Hsp90 sequences were consistent with trees inferred from rDNA sequences in that the backbone of the dinoflagellate clade was largely unresolved.

Conclusions/Significance

The sequence conservation in both Hsp90 and rDNA sequences and the poor resolution of the deepest nodes suggests that dinoflagellates reflect an explosive radiation in morphological diversity in their recent evolutionary past. Nonetheless, the more comprehensive analysis of Hsp90 sequences enabled us to infer phylogenetic interrelationships of dinoflagellates more rigorously. For instance, the phylogenetic position of Noctiluca, which possesses several unusual features, was incongruent with previous phylogenetic studies. Therefore, the generation of additional dinoflagellate Hsp90 sequences is expected to refine the stem group of athecate species observed here and contribute to future multi-gene analyses of dinoflagellate interrelationships.     

Positive selection in AvrP4 avirulence gene homologues across the genus Melampsora

Pathogen genes involved in interactions with their plant hosts are expected to evolve under positive Darwinian selection or balancing selection. In this study a single copy avirulence gene, AvrP4, in the plant pathogen Melampsora lini, was used to investigate the evolution of such a gene across species. Partial translation elongation factor 1-alpha sequences were obtained to establish phylogenetic relationships among the Melampsora species. We amplified AvrP4 homologues from species pathogenic on hosts from different plant families and orders, across the inferred phylogeny. Translations of the AvrP4 sequences revealed a predicted signal peptide and towards the C-terminus of the protein, six identically spaced cysteines were identified in all sequences. Maximum likelihood analysis of synonymous versus non-synonymous substitution rates indicated that positive selection played a role in the evolution of the gene during the diversification of the genus. Fourteen codons under significant positive selection reside in the C-terminal 28 amino acid region, suggesting that this region interacts with host molecules in most sequenced accessions. Selection pressures on the gene may be either due to the pathogenicity or avirulence function of the gene or both.