Nucleotide diversity and molecular evolution of the WAG-2 gene in common wheat (Triticum aestivum L) and its relatives

In this work, we examined the genetic diversity and evolution of the WAG-2 gene based on new WAG-2 alleles isolated from wheat and its relatives. Only single nucleotide polymorphisms (SNP) and no insertions and deletions (indels) were found in exon sequences of WAG-2 from different species. More SNPs and indels occurred in introns than in exons. For exons, exons+introns and introns, the nucleotide polymorphism π decreased from diploid and tetraploid genotypes to hexaploid genotypes. This finding indicated that the diversity of WAG-2 in diploids was greater than in hexaploids because of the strong selection pressure on the latter. All dn/ds ratios were < 1.0, indicating that WAG-2 belongs to a conserved gene affected by negative selection. Thirty-nine of the 57 particular SNPs and eight of the 10 indels were detected in diploid species. The degree of divergence in intron length among WAG-2 clones and phylogenetic tree topology suggested the existence of three homoeologs in the A, B or D genome of common wheat. Wheat AG-like genes were divided into WAG-1 and WAG-2 clades. The latter clade contained WAG-2, OsMADS3 and ZMM2 genes, indicating functional homoeology among them.     

The analysis of nucleotide substitutions, gaps, and recombination events between RHD and RHCE genes through complete sequencing

We determined the entire nucleotide sequences of all introns within the RHD and RHCE genes by amplifying genomic DNA using long PCR methods. The RHD and RHCE genes were 57,295 and 57,831 bp in length, respectively. Aligning both genes revealed 138 gaps (insertions and deletions) below 100 bp, 1116 substitutions in all introns and all exons (coding region), and 5 gaps of over 100 bp. Homologies (%) between the RH genes were 93.8% over all introns and coding exons and 91.7% over all exons and introns. Various short tandem repeats (STRs) and many interspersed nuclear elements were identified in both genes. The proportions of Alu sequences in the RHD and RHCE genes were 25.9 and 25.7%, respectively and these Alu sequences were concentrated in several regions. We confirmed multiple recombinations in introns 1 and 2. Such multiple recombination, which probably arose due to the concentrations of Alu sequences and the high level of the homology (%), is one of most important factors in the formation and evolution of RH gene. The variability of the Rh system may be generated because of these features of RH genes. Apparent mutational hotspots and regions with low of K values (the numbers of substitutions per nucleotide site) caused by recombinations as well as true mutational hotspots may be found in human genome. Accordingly, in searching for and identifying single nucleotide polymorphisms (SNPs) especially in noncoding regions, apparent mutational hotspots and areas of low K values by recombination should be noted since the unequal distribution of SNPs will reduce the power of SNPs as genetic maker. Combining the complete sequences' data of both RH genes with serological findings will provide beneficial information with which to elucidate the mechanism of recombination, mutation, polymorphism, and evolution of other genes containing the RH gene as well as to analyze Rh variants and develop new methods of Rh genotyping.     

DNA polymorphism, haplotype structure and balancing selection in the Leavenworthia PgiC locus.

A study of DNA polymorphism and divergence was conducted for the cytosolic phosphoglucose isomerase (PGI:E.C.5.3.1.9) gene of five species of the mustard genus Leavenworthia: Leavenworthia stylosa, L. alabamica, L. crassa, L. uniflora, and L. torulosa. Sequences of an internal 2.3-kb PgiC gene region spanning exons 6-16 were obtained from 14 L. stylosa plants from two natural populations and from one to several plants for each of the other species. The level of nucleotide polymorphism in L. stylosa PgiC gene was quite high (pi = 0.051, theta = 0.052). Although recombination is estimated to be high in this locus, extensive haplotype structure was observed for the entire 2.3-kb region. The L. stylosa sequences fall into at least two groups, distinguished by the presence of several indels and nucleotide substitutions, and one of the three charge change nucleotide replacements within the region sequenced correlates with the haplotypes. The differences between the haplotypes are older than between the species, and the haplotypes are still segregating in at least two of five species studied. There is no evidence of recent or ancient population subdivision that could maintain distinct haplotypes. The age of the haplotypes and the results of Kelly's Z(nS) and Wall's B and Q tests with recombination suggest that the haplotypes are maintained due to balancing selection at or near this locus.     

Determination of the gene structure of human oligophrenin-1 and identification of three novel polymorphisms by screening of DNA from 164 patients with non-specific X-linked mental retardation

We have recently shown that mutations in oligophrenin-1 (OPHN1) are responsible for non-specific X-linked mental retardation (MRX). The structure of the gene encoding the OPHN1 protein was determined by isolation of genomic DNA clones from the human cosmid library. Genomic fragments containing exons were sequenced, and the sequences of the exons and flanking introns were defined. Knowledge of the genomic structure of the OPHN1 gene, which spans at least 500 kb and consists of 25 exons, will facilitate the search for additional mutations in OPHN1. OPHN1 was screened for mutations in 164 subjects with non-specific mental retardation. Three nucleotide substitutions were identified, one of which was a silent mutation in the codon threonine 301 at position 903 (G-->C). The other substitutions were located in exon 2, a G-->A substitution at position 133 (A45T), and in exon 10, a C-->T substitution at position 902 (T301M), but these are common polymorphisms rather than disease-causing mutations.     

Identification of 167 polymorphisms in 88 genes from candidate neurodegeneration pathways

Catalogs of intra-gene polymorphisms are needed to facilitate wide-ranging candidate gene-based association studies in common complex diseases. With this in mind, we have scanned multiple alignments of expressed sequence tags and of genomic DNA sequences (PCR products from four to eight unrelated individuals) to find polymorphisms in 195 genes putatively involved in neurodegenerative illness (including components of oxidative stress, excitotoxicity, inflammation, apoptosis and aging). This led to the discovery of 167 polymorphisms in 88 genes. These comprised 163 single nucleotide polymorphisms, one insertion/deletion, and three other variations involving more than one base pair. The polymorphisms were distributed in the exons (87), introns (70), and gene flanking regions (10). Of the exonic polymorphisms, 17 would give rise to non-synonymous amino acid substitutions. These findings now provide a valuable resource for association studies in neurodegenerative disorders such as Alzheimer's disease and Parkinson's disease.     

Discovery of intron polymorphisms in cultivated tomato using both tomato and Arabidopsis genomic information

A low level of genetic variation has limited the application of molecular markers for characterizing important traits in cultivated tomato. To detect polymorphisms in tomato conserved ortholog sets (COS), expressed sequence tags (ESTs) were searched against tomato and Arabidopsis genomic sequences to define the positions of introns. Introns were amplified from 12 different accessions of tomato by polymerase chain reaction and nucleotide sequences were determined by sequencing. Results indicated that there was a possibility of 71% to amplify introns from tomato genomic DNA through this approach. A total of 201 introns were sequenced from 86 COS unigenes. The intron positions and numbers were conserved between tomato and Arabidopsis, but average intron length was three times longer in tomato than in Arabidopsis. A total of 307 single nucleotide polymorphisms (SNPs) and 75 indels were detected in introns of 57 COS unigenes among 12 tomato lines. Within cultivated tomato germplasm 172 SNPs and 47 indels were detected in introns of 33 COS unigenes. In addition, 41 SNPs were identified in the exons of 27 COS unigenes. The frequency of SNPs was 2.4 times higher in introns than in exons in the 22 COS unigenes having both intronic and exonic polymorphisms. These results indicate that intronic regions may contain sufficient variation to develop sufficient marker resources for genome-wide analysis in cultivated tomato.     

Characterisation of the genetic effects of the ADFP gene and its association with production traits in dairy goats

Adipose differentiation-related protein (ADFP) is important for regulation of lipid metabolism and insulin secretion in beta-cells. In this study, we investigated polymorphisms within the caprine ADFP gene and determined its relationship with production traits. As there was no sequence information available for the caprine ADFP gene, we generated DNA sequence data and examined the genomic organisation. The caprine ADFP gene is organised into 7 exons and 6 introns that span approximately 8.7 kbp and is transcribed into mRNA containing 1353 bp of sequence coding for a protein of 450 amino acids. The protein sequences showed substantial similarity (71–99%) to orthologues from cattle, human and mouse. We identified polymorphisms in the sequences using DNA sequencing, PCR-RFLP and forced PCR-RFLP methods. Seven single nucleotide polymorphisms (SNPs) were identified using samples from 4 different goat populations consisting of 1408 healthy and unrelated individuals. Six haplotypes involving the 7 SNPs from the caprine ADFP gene were identified and their effects on production traits were analysed. Haplotype 6 had the highest haplotype frequency and was highly significantly associated with chest circumference and milk yield in the analysed populations. The results of this study suggest that the ADFP gene is a strong candidate gene affecting production traits and may be used for marker-assisted selection and management in Chinese dairy goat breeding programmes.     

Polymorphism at the ribosomal DNA spacers and its relation to breeding structure of the widespread mushroom Schizophyllum commune

The common split-gilled mushroom Schizophyllum commune is found throughout the world on woody substrates. This study addresses the dispersal and population structure of this fungal species by studying the phylogeny and evolutionary dynamics of ribosomal DNA (rDNA) spacer regions. Extensive sampling (n = 195) of sequences of the intergenic spacer region (IGS1) revealed a large number of unique haplotypes (n = 143). The phylogeny of these IGS1 sequences revealed strong geographic patterns and supported three evolutionarily distinct lineages within the global population. The same three geographic lineages were found in phylogenetic analysis of both other rDNA spacer regions (IGS2 and ITS). However, nested clade analysis of the IGS1 phylogeny suggested the population structure of S. commune has undergone recent changes, such as a long distance colonization of western North America from Europe as well as a recent range expansion in the Caribbean. Among all spacer regions, variation in length and nucleotide sequence was observed between but not within the tandem rDNA repeats (arrays). This pattern is consistent with strong within-array and weak among-array homogenizing forces. We present evidence for the suppression of recombination between rDNA arrays on homologous chromosomes that may account for this pattern of concerted evolution.     

Population genetic analysis of the N-acylsphingosine amidohydrolase gene associated with mental activity in humans

To understand the evolution of human mental activity, we performed population genetic analyses of nucleotide sequences ( approximately 11 kb) from a worldwide sample of 60 chromosomes of the N-acylsphingosine amidohydrolase (ASAH1) gene. ASAH1 hydrolyzes ceramides and regulates neuronal development, and its deficiency often results in mental retardation. In the region ( approximately 4.4 kb) encompassing exons 3 and 4 of this gene, two distinct lineages (V and M) have been segregating in the human population for 2.4 +/- 0.4 million years (MY). The persistence of these two lineages is attributed to ancient population structure of humans in Africa. However, all haplotypes belonging to the V lineage exhibit strong linkage disequilibrium, a high frequency (62%), and small nucleotide diversity (pi = 0.05%). These features indicate a signature of positive Darwinian selection for the V lineage. Compared with the orthologs in mammals and birds, it is only Val at amino acid site 72 that is found exclusively in the V lineage in humans, suggesting that this Val is a likely target of positive selection. Computer simulation confirms that demographic models of modern humans except for the ancient population structure cannot explain the presence of two distinct lineages, and neutrality is incompatible with the observed small genetic variation of the V lineage at ASAH1. On the basis of the above observations, it is argued that positive selection is possibly operating on ASAH1 in the modern human population.     

The evolution of an alpha-esterase pseudogene inactivated in the Drosophila melanogaster lineage

Previous analyses of the alpha-esterase cluster of Drosophila melanogaster revealed 10 active genes and the DmalphaE4a-Psi pseudogene. Here, we reconstruct the evolution of the pseudogene from the sequences of 12 alleles from widely scattered D. melanogaster populations and single alleles from Drosophila simulans and Drosophila yakuba. All of the DmalphaE4a-Psi alleles contain numerous inactivating mutations, suggesting that pseudogene alleles are fixed in natural populations. Several lines of evidence also suggest that DmalphaE4a is now evolving without selective constraint in the D. melanogaster lineage. There are three polymorphic indels which result in frameshifts; a key nucleotide of the intron splice acceptor is polymorphic; the neutral mutation parameter is the same for replacement and silent sites; one of the nonsilent polymorphisms results in a stop codon; only 1 of the 13 replacement polymorphisms is biochemically conservative; residues that are conserved among active esterases have different states in DmalphaE4a-Psi; and there are about half as many transitional polymorphisms as transversional ones. In contrast, the D. simulans and D. yakuba orthologs DsalphaE4a and DyalphaE4a do not have the inactivating mutations of DmalphaE4a-Psi and appear to be evolving under the purifying selection typical of protein- encoding genes. For instance, there have been more substitutions in the introns than in the exons, and more in silent sites than in replacement sites. Furthermore, most of the amino acid substitutions that have occurred between DyalphaE4a and DsalphaE4a are located in sites that typically vary among active alpha-esterases rather than those that are usually conserved. We argue that the original alphaE4a gene had a function which it has lost since the divergence of the D. melanogaster and D. simulans lineages.     

Phylogenetic relationships among bisexual and unisexual lineages of the weevil Scepticus insularis (Coleoptera: Curculionidae) on Hokkaido, Northern Japan

We examined the geographic distributions and phylogenetic relationships of bisexual and unisexual (parthenogenetic) forms of the weevil Scepticus insularis on Hokkaido Island, northern Japan. Unisexual beetles were widespread throughout Hokkaido, whereas bisexuals were found only in three remote areas. Bisexuals (females and males) and unisexual females occurred sympatrically in two areas. We determined nucleotide sequences for part of the mitochondrial NADH dehydrogenase subunit 2 (ND2; 633 bp) gene for 104 individuals, and for part of the nuclear internal transcribed spacer region 2 (ITS2; 360 bp) for 91 individuals. In an ND2 gene tree, haplotypes of S. insularis fell into two distinct clades (A and B), which were genetically differentiated from one another by 9.1% nucleotide sequence divergence. Haplotypes of females identified as unisexual were all in clade A, whereas those of females identified as bisexual belonged to clade B. Haplotypes of males were in clade B, except for two males having a clade-A haplotype. Circumstantial evidence suggests that these two males were produced by unisexual females. The ND2 tree suggests that the current unisexual form of S. insularis on Hokkaido was of a single origin. In contrast, a gene tree for ITS2 haplotypes show no clear divergence between the two modes of reproduction, with two major haplotypes shared by unisexual females, bisexual females, and males. This incongruence between the nuclear and mitochondrial phylogenies may be attributable to occasional gene flow between the unisexual and bisexual lineages through males occasionally produced by unisexual females, but our results do not exclude the possibility that the two lineages share polymorphic ancestral ITS2 haplotypes.     

中国狼不同地理种群遗传多样性及系统进化分析

为了解中国狼不同地理种群遗传多样性及系统发育情况,从中国境内狼的主要分布区青海、新疆、内蒙古和吉林4个地区采集样品,用分子生物学技术手段成功地获得44个个体线粒体DNA控制区第一高变区(HVRⅠ)序列和40个线粒体Cyt b部分序列。线粒体控制区HVRⅠ共检测到51个变异位点,位点变异率为8.76%;线粒体Cyt b部分序列发现31个变异位点,位点变异率为5.33%,未见插入及缺失现象,变异类型全部为碱基置换。共定义了16个线粒体HVRⅠ单倍型,其中吉林与内蒙种群存在共享单倍型,估计这两地间种群亲缘关系较近。4个地理种群中新疆种群拥有较高的遗传多样性(0.94)。中国狼种群总体平均核苷酸多态性为2.27%,与世界其他国家地区相比,中国狼种群拥有相对较高的遗传多样性。通过线粒体HVRⅠ单倍型构建的系统进化树可以看出,中国狼在进化上分为2大支,其中位于青藏高原的青海种群独立为一支,推测其可能长期作为独立种群进化。基于青海种群与新疆,内蒙种群的线粒体Cyt b遗传距离,推测中国狼2个世系可能在更新世冰川时期青藏高原受地质作用急速隆起后出现分歧,分歧时间大约在1.1 MY前。     

A recently silenced, duplicate PgiC locus in Clarkia

Previous electrophoretic analysis showed that 17 diploid species of the wildflower Clarkia (Onagraceae) have two cytosolic isozymes of phosphoglucose isomerase (PGIC; EC 5.3.1.9), whereas 15 other diploid species have a single PGIC. Molecular studies revealed that the two isozymes in the former species are encoded by duplicate genes, PgiC1 and PgiC2, whereas the single isozyme in the latter is always encoded by PgiC1. Phylogenetic analysis of the nucleotide sequences implied that PgiC2 was silenced four times independently in the genus. Here we describe a psi PgiC2 from C. mildrediae, a species in which only PgiC1 is expressed. The discovery of the psi PgiC2 is significant because it confirms a formal prediction of the phylogenetic analysis. The psi PgiC2 includes 5,039 nucleotides corresponding to 18 of the 23 exons of PgiC, as well as the intervening introns and 3' nontranslated region. The absence of an increase of nucleotide substitutions in its "exons" suggests that the gene was silenced recently. The present study appears to be the first to establish that a specific duplicate gene locus regularly expressed in a group of related plant species has been silenced in one of them. The multiple independent silencings of PgiC2 suggest that it remained functional but inessential in ancestral lineages. We discuss the possibility that PgiC2 may have been preserved in these lineages by selection against mutants causing defective PGIC1- PGIC2 heterodimers.      

Structure of the human ornithine transcarbamylase gene

Complementary and genomic DNA clones corresponding to the human ornithine transcarbamylase (OTC) [EC 2.1.3.3]mRNA have been isolated and analyzed. The OTC gene is about 73 kilobase pairs (kb) long and contains 10 exons interrupted by 9 introns of highly variable sizes. The smallest intron is 80 base pairs and the largest, 21.7 kb. The 5'- and 3'-flanking regions, entire exons and all the exon/intron boundaries were sequenced. The nucleotide and deduced amino acid sequences of isolated OTC cDNAs as well as the corresponding regions of the genomic DNA were compared with those of human OTC cDNA (Horwich, A.L., Fenton, W.A., Williams, K.R., Kalousek, F., Kraus, J.P., Doolittle, R.F., Koningsberg, W., & Rosenberg, L.E. (1984) Science 224, 1068-1074). We found 20 nucleotide substitutions among these sequences, of which 6 related to amino acid changes. The nature of these nucleotide substitutions is discussed.     

High diversity of the chicken growth hormone gene and effects on growth and carcass traits

The chicken growth hormone (cGH) gene plays a crucial role in controlling growth and metabolism, leading to potential correlations between cGH polymorphisms and economic traits. In this study, DNA from four divergent chicken breeds were screened for single nucleotide polymorphisms (SNPs) in the cGH gene using denaturing high-performance liquid chromatography and sequencing. A total of 46 SNPs were identified, of which 4 were in the 5' untranslated region, 1 in the 3' untranslated region, 5 in exons (two of which are nonsynonymous), with the remaining 36 in introns. The nucleotide diversity in the cGH gene ( theta = 2.7 x 10(-3)) was higher than that reported for other chicken genes, even within the same breeds. The associations of five of these SNPs and their haplotypes with chicken growth and carcass traits were determined using polymerase chain reaction-restriction fragment length polymorphism analysis in a F2 resource population cross of two of the four chicken breeds (White Recessive Rock and Xinghua). This analysis shows that, among other correlations, G+1705A was significantly associated with body weight at all ages measured, shank length at three of four ages measured, and average daily gain within weeks 0 to 4. Thus, this cGH polymorphism, or another polymorphism that is in linkage disequilibrium with G+1705A, appears to correspond to a significant growth-related quantitative trait locus difference between the two breeds used to construct the resource population.     

Excessive variation in Y chromosomal DNA in Rumex acetosa (Polygonaceae)

Rumex acetosa is one of the few angiosperms that possesses sex chromosomes. The same types of abundant repetitive sequences cover both heterochromatic Y chromosomes present in males. The aim of this study was to investigate genetic variation in paternally inherited Y chromosomal DNA and in maternally inherited cpDNA, and to find out whether the examined genomic regions are suited to a phylogeographic study in R. acetosa. DNA sequence polymorphisms present in the 850-bp heterochromatic segment on the Y chromosomes were compared to variation in the 409-bp long chloroplast section (trnL- trnF spacer) in R. acetosa originating from several European locations and from the Altai mountains in Russia. A great amount of genetic variation was detected within the Y chromosomal region while only four chloroplast genotypes were detected. Although the chloroplast haplotypes possessed some geographic pattern, no clear phylogeographic pattern was detected based on the variable Y chromosomes. The mean Y chromosomal nucleotide diversity among all samples equaled 6.6 %, and the mean proportion of polymorphic sites per individual equaled 8.2 % among SNP sites and 1.7 % among all sites investigated. The high number of substitutions detected in the Y chromosomal DNA shows that this heterochromatic sequence has a high mutation rate. The diversity pattern indicates that gene flow via pollen is extensive and it blurs any geographical pattern in the Y chromosomal variation. The high number of repeats and uncertainty concerning the extent of recombination between the two Y chromosomes impair the usability of the Y chromosomal segment for phylogeographic or population genetic studies.     

Genealogy of the "Green Revolution" gene in rice

During the "Green Revolution" of rice, high-yielding varieties (HYVs) were developed using a semi-dwarf gene (sd1 or OsGA20ox2). The presence or absence of the two mutant alleles (DGWG type in Dee-geo-woo-gen and JKK type in Jikkoku) were surveyed by PCR using 256 accessions of eight wild and two cultivate rice species. The DGWG allele was detected in a landrace (Oryza sativa) and two accessions of wild rice (O. rufipogon), all of which are from China, showing their limited distribution. Genealogical studies of the OsGA20ox2 gene showed that the 62 sequences of O. sativa and O. rufipogon included 20 distinct haplotypes, indicating that the species complex contained OsGA20ox2 genes from two different lineages. The silent site nucleotide diversities (pi and theta(w)) were extremely low in Japonica rice, suggesting a genetic bottleneck. The haplotype network showed that the DGWG and JKK alleles were derived in different lineages. The DGWG carrier (W1944) had unique polymorphisms in the surrounding region of the locus, suggesting that the DGWG allele has been preserved in the wild progenitor, rather than that the DGWG allele has been introgressed from HYVs to W1944. Although a semi-dwarfing plant is a weak competitor under saturated fields, the crossing experiment revealed that the DGWG variant might have been preserved as a hidden variation in the genetic background of wild rice, without expressing a short-stature.     

Nucleotide variability at the acetyl coenzyme A carboxylase gene and the signature of herbicide selection in the grass weed Alopecurus myosuroides (Huds.)

Acetyl coenzyme A carboxylase (ACCase) is the target of highly effective herbicides. We investigated the nucleotide variability of the ACCase gene in a sample of 18 black-grass (Alopecurus myosuroides [Huds.]) populations to search for the signature of herbicide selection. Sequencing 3,396 bp encompassing ACCase herbicide-binding domain in 86 individuals revealed 92 polymorphisms, which formed 72 haplotypes. The ratio of nonsynonymous versus synonymous substitutions was very low, in agreement with ACCase being a vital metabolic enzyme. Within black grass, most nonsynonymous substitutions were related to resistance to ACCase-inhibiting herbicides. Differentiation between populations was strong, in contrast to expectations for an allogamous, annual plant. Significant H tests revealed recent hitchhiking events within populations. These results were consistent with recent and local positive selection. We propose that, although they have only been used since at most 15 black-grass generations, ACCase-inhibiting herbicides have exerted a positive selection targeting resistant haplotypes that has been strong enough to have a marked effect upon ACCase nucleotide diversity. A minimum-spanning network of nonrecombinant haplotypes revealed multiple, independent apparitions of resistance-associated mutations. This study provides the first evidence for the signature of ongoing, recent, pesticide selection upon variation at the gene encoding the targeted enzyme in natural plant populations.     

Contrasting selection modes at the Adh1 locus in outcrossing Miscanthus sinensis vs. inbreeding Miscanthus condensatus (Poaceae)

We estimated DNA sequence variation of the Adh1 locus in the outcrossing Miscanthus sinensis (Poaceae) and its close selfing relative, M. condensatus. Tajima's test of selection is significantly negative for both overall exons and replacement sites in M. sinensis. Among its entire sample, nucleotide diversity of nonsynonymous sites is higher than that of synonymous sites. A McDonald and Kreitman test of neutrality indicates an excess of intraspecific replacement polymorphisms, suggesting possible directional selection toward advantageous mutants. However, frequent intragenic recombination suggests both purifying and positive selection is unlikely. Recent demographic expansions coupled with relaxation of purifying selection may have resulted in elevated genetic diversity at the Adh1 locus as well as the trnL-trnF intergenic spacer of cpDNA in this outcrossing species. In contrast, low levels of genetic diversity were detected at both the Adh1 locus and the cpDNA spacer in M. condensatus, consistent with bottlenecks associated with selfing in all populations. While Tajima's D and Fu and Li's F statistics did not reveal deviation from neutrality at the Adh1 locus in M. condensatus, 12 replacements vs. 10 synonymous changes were detected. Based on pairwise comparisons of the d(N)/d(S) ratio, lineages of closely related populations of the species distributed along saline habitats appeared to be under directional selection.