Genetic diversity among Brazilian isolates of Beauveria bassiana: comparisons with non-Brazilian isolates and other Beauveria species

Aims:  The genetic diversity of Beauveria bassiana was investigated by comparing isolates of this species to each other (49 from different geographical regions of Brazil and 4 from USA) and to other Beauveria spp.
Methods and Results:  The isolates were examined by multilocus enzyme electrophoresis (MLEE), amplified fragment length polymorphism (AFLP), and rDNA sequencing. MLEE and AFLP revealed considerable genetic variability among B. bassiana isolates. Several isolates from South and Southeast Brazil had high similarity coefficients, providing evidence of at least one population with clonal structure. There were clear genomic differences between most Brazilian and USA B. bassiana isolates. A Mantel test using data generated by AFLP provided evidence that greater geographical distances were associated with higher genetic distances. AFLP and rDNA sequencing demonstrated notable genotypic variation between B. bassiana and other Beauveria spp.
Conclusion:  Geographical distance between populations apparently is an important factor influencing genotypic variability among B. bassiana populations in Brazil.
Significance and Impact of the Study:  This study characterized many B. bassiana isolates. The results indicate that certain Brazilian isolates are considerably different from others and possibly should be regarded as separate species from B. bassiana sensu latu . The information on genetic variation among the Brazilian isolates, therefore, will be important to comprehending the population structure of B. bassiana in Brazil.     

Genetic diversity, reproductive biology, and speciation in the entomopathogenic fungus Beauveria bassiana (Balsamo) Vuillemin.

Beauveria bassiana, a mitosporic fungus used for the biological control of many insect species, is recognized as a "species complex" comprising genetically diverse lineages. Being predominantly asexual, mating tests cannot be applied to delimit species in this species complex. Genetic tests offer an indirect means of identifying species among isolates. To this end, molecular genetic analysis of a sample of B. bassiana isolates with 2 subsamples, 1 representing a worldwide collection and another from a localized epizootic population was carried out. DNA markers generated through AFLPs (amplified fragment length polymorphisms) and SSCPs (single-strand conformation poly morphisms) and nucleotide sequence data of different allelic forms of 3 genes (large and small subunits of rRNA and beta-tubulin) were evaluated. The B. bassiana isolates from the worldwide sample showed 11% overall similarity and no closely clustered groups. Phylogenetic trees generated from the AFLP and SSCP data of this sample resolved the different isolates into distinct phylogenetic lineages. In the epizootic B. bassiana population, prevalence of recombination was evident from random association of alleles in multilocus tests and lack of phylogenetic concordance among 3 gene genealogies. Thus, the worldwide sample of B. bassiana exhibits a predominantly clonal structure, hinting at species divergence leading to cryptic speciation with recombination being customary among isolates sharing a close ecological niche.     

Antigenic relationships among pathogenic Beauveria bassiana with Engyodontium album (=B. alba) and non-pathogenic species of the genus Beauveria

Exoantigenic extracts of 15 isolates belonging to hyalohyphomycosis-causing Beauveria bassiana (1), and Engyodontium album (1), as well as other species of the genus Beauveria (one isolate each of B. brogniartii, B. densa, B. stephanoderis, B. velata, B. vermiconia and six isolates of unknown Beauveria species) were studied. Aqueous-merthiolated extracts derived from 10-day-old Sabouraud's dextrose agar slant cultures (25 °C) were concentrated (25X), and reacted against rabbit anti-B. bassiana serum in the presence of partially purified homologous antigen (20X) prepared from 5-week-old shaken cultures (30 °C), using a microimmunodiffusion procedure. Beauveria bassiana reference antigen and antiserum reacted to produce four bands of identity. With the exception of E. album, which was negative, extracts of the isolates of B. brogniartii, B. densa, B. stephanoderis, B. velata, B. vermiconia and the unknown Beauveria species all produced 2-4 lines of identity against the homologous anti-B. bassiana serum. These results suggested that all the species of the genus Beauveria tested were antigenically related to B. bassiana. Engyodontium album demonstrated antigenic distinctness, however, from B. bassiana and thus supported the validity of this taxon. This revised version was published online in August 2006 with corrections to the Cover Date.     

What phylogeny and gene genealogy analyses reveal about homoplasy in citrus microsatellite alleles

Sixty-five microsatellite alleles amplified from ancestral citrus accessions classified in three separate genera were evaluated for sequence polymorphism to establish the basis of inter- and intra-allelic genetic variation, evaluate the extent of size homoplasy, and determine an appropriate model (stepwise or infinite allele) for analysis of citrus microsatellite alleles. Sequences for each locus were aligned and subsequently used to determine relationships between alleles of different taxa via parsimony. Interallelic size variation at each SSR locus examined was due to changes in repeat copy number with one exception. Sequencing these alleles uncovered new distinct point mutations in the microsatellite region and the region flanking the microsatellite. Several of the point mutations were found to be genus, species, or allele specific, and some mutations were informative about the inferred evolutionary relationships among alleles. Overall, homoplasy was observed in alleles from all three loci, where the core microsatellite repeat was changed causing alleles of the same size class to be identical in state but not identical by descent. Because nearly all changes in allele size (with one exception) were due to expansion or contraction of the repeat motif, this suggests that a stepwise mutation model, which assumes homoplasy may occur, would be the most appropriate for analyzing Citrus SSR data. The collected data indicate that microsatellites can be a useful tool for evaluating Citrus species and two related genera since repeat motifs were reasonably well retained. However, this work also demonstrated that the number of microsatellite alleles is clearly an underestimate of the number of sequence variants present.     

Community composition, host range and genetic structure of the fungal entomopathogen Beauveria in adjoining agricultural and seminatural habitats

Although intensively investigated for biological control of insect pests, little is known about the ecology of the fungal entomopathogenic genus Beauveria in natural or agricultural habitats. In this study, we used molecular phylogenetic and genotypic information to infer species diversity, reproductive potential and genetic structure of Beauveria occurring within a single arable field and bordering hedgerow in Denmark. Isolates were sampled from cultivated field and hedgerow soils, from insects harbouring latent fungal infections, and from the phylloplanes of three plant species common in the hedgerow flora. A nuclear phylogeny of this local Beauveria assemblage resolved seven phylogenetic species, including (i) five phylogenetic species within Beauveria bassiana sensu stricto ; (ii) Clade C, a taxonomically uncharacterized species that is morphologically indistinguishable but phylogenetically distant from B. bassiana s.s. ; and (iii) Beauveria brongniartii. All seven species were present throughout the hedgerow habitat, including as infections in insects. Significantly, only B . bassiana s.s. phylogenetic species Eu_1 was isolated from tilled soils. Mating type polymerase chain reaction assays demonstrated that all five B. bassiana s.s. phylogenetic species possess bipolar outcrossing mating systems. Of these, only the Eu_1 population contained two mating types; however, a 31:2 skew in MAT1:MAT2 mating types suggests a low frequency of sexual reproduction in this population. The four remaining B. bassiana s.s. phylogenetic species were fixed for single mating types and these populations are evidently clonal. Multilocus microsatellite genotyping revealed polymorphism in all five phylogenetic species of B. bassiana s.s. ; however, all show evidence of clonal genetic structure.     

Insertions–Deletions in a Microsatellite Flanking Region May Be Resolved by Variation in Stuttering Patterns

Microsatellites or simple sequence repeats (SSRs) may display polymerase-chain-reaction-amplified fragment lengths mismatching the patterns expected from repeat copy number variation. We sequenced alleles of a nuclear dinucleotide SSR locus in two oak species which showed 2- and 1-bp length differences between alleles and three types of stuttering patterns in fragment length analysis. In accordance with the variation in stuttering, we identified three allele classes characterized by insertions–deletions in the flanking regions and overlapping repeat copy number ranges. Different alleles could thus only be safely separated when considering these stuttering patterns. Our results raise the question of how to adequately delimit alleles when such size homoplasy is present. We advise to thoroughly characterize SSR sequence variation during marker development and to carefully place primer sites along flanking regions to facilitate automated allele scoring and to minimize labor-intensive visual inspection.     

利用28S rDNAI型内含子研究琅琊山球孢白僵菌菌株的遗传多样性

通过四对引物对89株采自琅琊山的球孢白僵菌菌株进行PCR扩增,得到17类不同的单倍型,其中以BBBA型最多,占39．33％,而BABB型等其他9种单倍型较少,各占1．12％。该地区单倍型种类占已发现类型的57．7％,这说明在固定时间内同一地区的菌株也有较为丰富菌株类型,种群遗传结构是多样化的。具有较广寄主谱的BBBA型菌株是该地区的优势菌株,可作为生防菌株的重点筛选对象。     

利用28S rDNA Ⅰ型内含子研究琅琊山球孢白僵菌菌株的遗传多样性

通过四对引物对89株采自琅琊山的球孢白僵菌菌株进行PCR扩增,得到17类不同的单倍型,其中以BBBA型最多,占39.33%,而BABB型等其他9种单倍型较少,各占1.12%。该地区单倍型种类占已发现类型的57.7%,这说明在固定时间内同一地区的菌株也有较为丰富菌株类型,种群遗传结构是多样化的。具有较广寄主谱的BBBA型菌株是该地区的优势菌株,可作为生防菌株的重点筛选对象。     

Variability in tolerance to UV-B radiation among Beauveria spp. isolates

Solar radiation, particularly the UV-B component, negatively affects survival of entomopathogenic fungi in the field. In an effort to identify Beauveria spp. isolates with promise for use in biological control settings with high insolation, we examined 53 Beauveria bassiana isolates, 7 isolates of 4 other Beauveria spp. and Engyodontium albus (=Beauveria alba). The origins of these fungi varied widely as to host/substrate and country, but approximately 30% of these isolates were B. bassiana from ticks in Brazil. A preliminary trial with three B. bassiana isolates (Bb 19, CG 310 and CG 481) at several UV-B dosages indicated that 2h of weighted UV-B irradiance at 978mWm(-2) (providing a total dose of 7.04kJm(-2)) allowed separation of isolates into low, medium or high UV-B tolerance. This dose, therefore, was selected as a single dose to compare UV-B tolerances of all 60 Beauveria spp. isolates. There was high variability in tolerance to UV-B radiation among the B. bassiana isolates, ranging from virtually zero tolerance (e.g., Bb 03) to almost 80% tolerance (e.g., CG 228). In addition, surviving B. bassiana conidia demonstrated delayed germination; and this is likely to reduce virulence. Conidia of the other species were markedly more sensitive to UV-B, with E. albus (UFPE 3138) being the least UV-B tolerant. Among B. bassiana isolates originating from 0 degrees to 22 degrees latitudes, those from lower latitudes demonstrated statistically significant greater UV-B tolerances than those isolates from higher latitudes. Isolates from above 22 degrees , however, were unaffected by latitude of origin. A similar analysis based on host type did not indicate a correlation between original host and UV-B tolerance. The identification in this study of several B. bassiana isolates with relatively high UV-B tolerance will guide the selection of isolates for future arthropod microbial control experiments.     

Sporulation of Metarhizium anisopliae and Beauveria bassiana on Coptotermes formosanus and in vitro

The sporulation of 22 total isolates of Metarhizium anisopliae and Beauveria bassiana was quantified on cadavers of the Formosan subterranean termite, Coptotermes formosanus. Conidial production increased significantly over 11 days post-death. Effects of isolates of M. anisopliae and B. bassiana on in vivo sporulation were significant. Although the overall effects of fungal species on in vivo sporulation were not significant, the interactions between fungal species and certain times post-death were significant, indicating different sporulation patterns between the two fungal species. B. bassiana isolates could be categorized into a group with high total sporulation (day 11) and low quick sporulation (on days 2 and 3), while M. anisopliae isolates fell into another group with high quick sporulation and low total sporulation. This could give M. anisopliae an advantage over B. bassiana in termite microbial control due to termite defensive social behaviors. Conidial production was significantly higher in vitro than in vivo. In vitro and in vivo sporulation differed by as much as 89x and 232x among the selected isolates of B. bassiana and M. anisopliae, respectively. Correlation between in vivo and in vitro conidial production was positive and significant. This may allow preliminary in vitro screening of a large number of isolates for high in vivo sporulation.     

A Beauveria phylogeny inferred from nuclear ITS and EF1-alpha sequences: evidence for cryptic diversification and links to Cordyceps teleomorphs

Beauveria is a globally distributed genus of soil-borne entomopathogenic hyphomycetes of interest as a model system for the study of entomopathogenesis and the biological control of pest insects. Species recognition in Beauveria is difficult due to a lack of taxonomically informative morphology. This has impeded assessment of species diversity in this genus and investigation of their natural history. A gene-genealogical approach was used to investigate molecular phylogenetic diversity of Beauveria and several presumptively related Cordyceps species. Analyses were based on nuclear ribosomal internal transcribed spacer (ITS) and elongation factor 1-alpha (EF1-alpha) sequences for 86 exemplar isolates from diverse geographic origins, habitats and insect hosts. Phylogenetic trees were inferred using maximum parsimony and Bayesian likelihood methods. Six well supported clades within Beauveria, provisionally designated A-F, were resolved in the EF1-alpha and combined gene phylogenies. Beauveria bassiana, a ubiquitous species that is characterized morphologically by globose to subglobose conidia, was determined to be non-monophyletic and consists of two unrelated lineages, clades A and C. Clade A is globally distributed and includes the Asian teleomorph Cordyceps staphylinidaecola and its probable synonym C. bassiana. All isolates contained in Clade C are anamorphic and originate from Europe and North America. Clade B includes isolates of B. brongniartii, a Eurasian species complex characterized by ellipsoidal conidia. Clade D includes B. caledonica and B. vermiconia, which produce cylindrical and comma-shaped conidia, respectively. Clade E, from Asia, includes Beauveria anamorphs and a Cordyceps teleomorph that both produce ellipsoidal conidia. Clade F, the basal branch in the Beauveria phylogeny includes the South American species B. amorpha, which produces cylindrical conidia. Lineage diversity detected within clades A, B and C suggests that prevailing morphological species concepts underestimate species diversity within these groups. Continental endemism of lineages in B. bassiana s.l. (clades A and C) indicates that isolation by distance has been an important factor in the evolutionary diversification of these clades. Permutation tests indicate that host association is essentially random in both B. bassiana s.l. clades A and C, supporting past assumptions that this species is not host specific. In contrast, isolates in clades B and D occurred primarily on coleopteran hosts, although sampling in these clades was insufficient to assess host affliation at lower taxonomic ranks. The phylogenetic placement of Cordyceps staphylinidaecola/bassiana, and C. scarabaeicola within Beauveria corroborates prior reports of these anamorph-teleomorph connections. These results establish a phylogenetic framework for further taxonomic, phylogenetic and comparative biological investigations of Beauveria and their corresponding Cordyceps teleomorphs.     

RAPD-PCR characterization of Beauveria bassiana isolates from the rice water weevil Lissorhoptrus oryzophilus

B.E. Urtz and W.C. Rice. 1997. Random amplified polymorphic DNA-polymerase chain reaction (RAPD-PCR) was used to examine the genetic variability among Beauveria bassiana isolates from infected rice water weevil (RWW), Lissorhoptrus oryzophilus , collected in south-western Louisiana. Most of the B. bassiana isolates tightly clustered into one or the other of two groups that diverged at the 65% similarity level. Three soil isolates also clustered within the two groups while one soil isolate did not tightly cluster (mean similarity = 65%) with any of the RWW isolates. The results suggest that certain genotypes of B. bassiana commonly infect RWW while others do not.     

Diversity of Beauveria spp. isolates from pollen beetles Meligethes aeneus in Switzerland

Pollen beetles Meligethes aeneus were collected in oilseed rape fields at different sites in Switzerland in spring 2004-2005 and 32 isolates of the fungal genus Beauveria occurring as latent infections in the beetles were obtained and molecularly characterized. Three major clades, Beauveria bassiana sensu stricto (Clade A: n=13), Beauveriabrongniartii (Clade B: n=1) and Beauveria Clade C (n=18) were identified among the isolates based on sequences of the ITS region and the 5' end of EF1-α. B. bassiana s.s. was further separated in the two clades, Eu_1 (n=10) and Eu_4 (n=3). The intergenic region Bloc provided best resolution of the individual clades B. bassiana s.s. Eu_1, Eu_4 and B. brongniartii. No specific clade of Beauveria appeared to be associated with adult M. aeneus populations. However, data suggested high relative abundance of Beauveria Clade C among the fungal entomopathogens infecting M. aeneus. Characterization of the isolates by simple sequence repeats (SSR) revealed further genotypic diversity within the clades except B. bassiana s.s. Eu_4 which appeared to be clonal. However, the individual SSR markers were differentially amplifiable from isolates of the different clades. It is therefore important to identify the underlying phylogenetic affinity of Beauveria isolates to interpret results based on SSR markers. The data suggest that not all available SSR markers are suitable for reliable characterization of diversity within Beauveria Clade C.     

Genetic effects on microsatellite diversity in wild emmer wheat (Triticum dicoccoides) at the Yehudiyya microsite,Israel

This study investigated allele size constraints and clustering, and genetic effects on microsatellite (simple sequence repeat, SSR) diversity at 28 loci comprising seven types of tandem repeated dinucleotide motifs in a natural population of wild emmer wheat, Triticum dicoccoides, from a shade vs sun microsite in Yehudiyya, northeast of the Sea of Galilee, Israel. It was found that allele distribution at SSR loci is clustered and constrained with lower or higher boundary. This may imply that SSR have functional significance and natural constraints. Genetic factors, involving genome, chromosome, motif, and locus significantly affected SSR diversity. Genome B appeared to have a larger average repeat number (ARN), but lower variance in repeat number (sigma(ARN)(2)), and smaller number of alleles per locus than genome A. SSRs with compound motifs showed larger ARN than those with perfect motifs. The effects of replication slippage and recombinational effects (eg, unequal crossing over) on SSR diversity varied with SSR motifs. Ecological stresses (sun vs shade) may affect mutational mechanisms, influencing the level of SSR diversity by both processes.     

Phylogenetic origins of African and Neotropical Beauveria bassiana s.l. pathogens of the coffee berry borer, Hypothenemus hampei

A phylogenetic epidemiological study of Beauveria bassiana s.l. was conducted for African and Neotropical pathogens of the coffee berry borer (CBB), Hypothenemus hampei, based on inferences from two nuclear intergenic regions, EFutr and Bloc. CBB pathogens were distributed among four terminal clades, however, the majority of African and Neotropical isolates cluster in a well-supported monophyletic group, informally designated AFNEO_1. Although the relationship between African and Neotropical AFNEO_1 is unresolved, the majority of alleles detected were exclusive to either the African or the Neotropical populations. These fixed genetic differences suggest that their disjunction predates the world trade in coffee. Neotropical AFNEO_1 have a broad host range and CBB pathogens are intermixed phylogenetically with isolates from diverse indigenous insects. Several Neotropical AFNEO_1 isolates were isolated from coffee plants as epiphytes or endophytes, thus plants themselves may potentially serve as reservoirs of pathogens against their insect pests. Topological incongruence between the EFutr and Bloc phylogenies of Neotropical AFNEO_1 may signify that individuals within this population are recombining.     

Assessing genetic diversity of wheat (<Emphasis Type="Italic">Triticum aestivum</Emphasis> L.) germplasm using microsatellite markers

A set of 24 wheat microsatellite markers, representing at least one marker from each chromosome, was used for the assessment of genetic diversity in 998 accessions of hexaploid bread wheat (Triticum aestivum L.) which originated from 68 countries of five continents. A total of 470 alleles were detected with an average allele number of 18.1 per locus. The highest number of alleles per locus was detected in the B genome with 19.9, compared to 17.4 and 16.5 for genomes A and D, respectively. The lowest allele number per locus among the seven homoeologous groups was observed in group 4. Greater genetic variation exists in the non-centromeric regions than in the centromeric regions of chromosomes. Allele numbers increased with the repeat number of the microsatellites used and their relative distance from the centromere, and was not dependent on the motif of microsatellites. Gene diversity was correlated with the number of alleles. Gene diversity according to Nei for the 26 microsatellite loci varied from 0.43 to 0.94 with an average of 0.77, and was 0.78, 0.81 and 0.73 for three genomes A, B and D, respectively. Alleles for each locus were present in regular two or three base-pair steps, indicating that the genetic variation during the wheat evolution occurred step by step in a continuous manner. In most cases, allele frequencies showed a normal distribution. Comparative analysis of microsatellite diversity among the eight geographical regions revealed that the accessions from the Near East and the Middle East exhibited more genetic diversity than those from the other regions. Greater diversity was found in Southeast Europe than in North and Southwest Europe. The present study also indicates that microsatellite markers permit the fast and high throughput fingerprinting of large numbers of accessions from a germplasm collection in order to assess genetic diversity.     

Comparative analysis of the production of insecticidal and melanizing macromolecules by strains of Beauveria spp.: in vivo studies

Eleven strains of Beauveria bassiana, and a further five species of Beauveria sp., were tested by injection of 8x10(2) conidia into the haemocoel of the larvae of the lepidopteran Galleria mellonella with the aim of analysing their toxin producing activity in vivo. Although the virulent strains killed 100% of the insects at slightly different rates (4-6 days) there were significant differences in the pattern and intensity of host melanization caused by isolates. The majority of the isolates of Beauveria spp. induced a fast and intense melanization of the cuticle of the integument and of tracheal wall, which followed one of three patterns. Another small group of two B. bassiana strains, isolated from Ostrinia nubilalis, induced very weak or no melanization. Strains 618 and 101 of B. bassiana, were selected as models of "melanizing" and "non-melanizing" strains, respectively. Ultrastructural alterations of cells of hypodermal and tracheal epithelium and of haemocytes, assumed to be at least partially caused by fungal toxins, were revealed in larvae infected by both isolates. However, their effects on the fine structure of the hypodermis were different. Injection of sera obtained from haemolymph of insects infected with B. bassiana 618 showed that they have insecticidal, melanizing, and cytotoxic effects similar to those occurring during mycosis. Chromatographic studies and bioassays with fractions prepared from crude serum have allowed a partial identification of the toxic molecules secreted by the fungus in vivo. They are proteinaceous, as shown by protease treatments, thermolabile, negatively charged, and not glycosylated with alpha-d-mannose or alpha-d-glucose. If strain B. bassiana 618 produces melanizing macromolecules which are vivotoxins secreted during the mycosis, the mode of action of isolate 101 is different. Its capacity to kill the host depends on active mycelial development, and on the production of low molecular weight toxins.     

An improved method for Beauveria bassiana transformation using phosphinothricin acetlytransferase and green fluorescent protein fusion gene as a selectable and visible marker

An improved transformation method for the biocontrol agent, Beauveria bassiana, was developed. For convenience of transformation selection and detection, the coding regions of the genes for phosphinothricin acetyltransferase and green fluorescent protein were fused and an expression vector, pBFT, carrying this fusion was constructed. Under optimum conditions, over 60 transformants microg(-1) plasmid DNA were obtained. B. bassiana conidia frozen 1 month at -80 degrees C were fully competent for transformation. The method was significantly less laborious and more rapid than current methods for B. bassiana. The bar::egfp provides a selectable and visible marker which may expedite future genetic engineering of this fungus.     

Comparison of five tandem repeat loci between humans and chimpanzees.

Five tandem repeat loci were studied in humans and chimpanzees using VNTR probes derived from human DNA. Shared alleles were found at three loci and were often the modal allele in one species but never in both. There was no difference in the mean number of alleles per locus. However, these species exhibited substantially different levels of gene diversity, with chimpanzees monomorphic at two loci. Evidence of reduced variability in chimpanzees corroborates earlier comparisons using isozymes and plasma proteins. Molecular mechanisms, population dynamics, or both may be responsible for these differences. Equal numbers of alleles per locus may reflect high mutation rates. By one test, chimpanzees were out of equilibrium at one locus, which may reflect a typing error or population substructure. The long divergence time, and the high probability of backward mutations, precludes accurate estimation of genetic distance between these species.     

Expressed sequence tag-simple sequence repeat-based molecular variance in two Salicornia (Amaranthaceae) populations

Salicornia spp is one of the most salt-tolerant vascular plants and is native to salt marshes and estuaries. We developed expressed sequence tag derived-simple sequence repeat (EST-SSR) markers for estimating genetic diversity and marker-assisted Salicornia breeding. Six polymorphic EST-SSRs of 40 detected 27 alleles, ranging from three to five alleles per locus. The average number of alleles per locus was 4.33 and 4.17, and the major allele frequency at locus DY529765 was high, being 0.859 and 0.857 in S. bigelovii and S. europea, respectively. Gene diversity, heterozygosity and polymorphism information content were highest at locus DY529950 and similar in these two species. Gene diversity increased with increase in the number of alleles that had a low major allele frequency at a locus. Six polymorphic loci effectively discriminated 46 taxa into three clusters via different analyses. Significant deviation of F(ST) from zero in three suggested populations for six loci indicated population differentiation and limited gene flow among them. A reduced median network established that taxon SB65 is primitive. SMART (simple modular architecture research tool) analysis of peptide sequences of six EST-SSRs showed that loci DY529765, DY529950 and EC906203 contained transmembrane, TLC, AgrB and NTR domains and might be involved in salinity stress tolerance. These EST-SSRs are a valuable resource for marker development and may be useful in marker-assisted Salicornia breeding.