An amphibian species pushed out of Britain by a moving hybrid zone

Classical theory states that hybrid zones will be stable in troughs of low population density where dispersal is hampered. Yet, evidence for moving hybrid zones is mounting. One possible reason that moving zones have been underappreciated is that they may drive themselves into oblivion and with just the superseding species remaining, morphological and genetic signals of past species replacement may be difficult to appreciate. Using genetic data (32 diagnostic single nucleotide polymorphisms) from a clinal hybrid zone of the common toad (Bufo bufo) and the spined toad (Bufo spinosus) in France for comparison, alleles of the latter species were documented in common toads in the south of Great Britain, at frequencies in excess of 10%. Because long distance dispersal across the Channel is unlikely, the conclusion reached was that the continental toad hybrid zone which previously extended into Britain, moved southwards and extirpated B. spinosus. Species distribution models for the mid‐Holocene and the present support that climate has locally changed in favour of B. bufo. The system bears resemblance with the demise of Homo neanderthalensis and the rise of Homo sapiens and provides an example that some paleoanthropologists demanded in support of a hominin “leaky replacement” scenario. The toad example is informative just because surviving pure B. spinosus and an extant slowly moving interspecific hybrid zone are available for comparison.     

Hybrid zone formation and contrasting outcomes of secondary contact over transects in common toads

Much progress in speciation research stems from documenting patterns of morphological and genetic variation in hybrid zones. Contrasting patterns of marker introgression in different sections of the contact can provide valuable insights on the relative importance of various evolutionary mechanisms maintaining species differences in the face of hybridization and gene flow and on hybrid zone temporal and spatial dynamics. We studied species interactions in the common toads Bufo bufo and B. spinosus in France and northwestern Italy using morphological and molecular data from the mitochondrial and nuclear genomes in an extensive survey, including two independent transects west and east of the Alps. At both, we found sharp, coincident and concordant nuclear genetic transitions. However, morphological clines were wider or absent and mtDNA introgression was asymmetric. We discuss alternative, nonexclusive hypotheses about evolutionary processes generating these patterns, including drift, selection, long‐distance dispersal and spatial shifts in hybrid zone location and structure. The distribution of intraspecific mtDNA lineages supports a scenario in which B. bufo held a local refugium during the last glacial maximum. Present‐day genetic profiles are best explained by an advance of B. spinosus from a nearby Iberian refugium, largely superseding the local B. bufo population, followed by an advance of B. bufo from the Balkans, with prongs north and south of the Alps, driving B. spinosus southwards. A pendulum moving hybrid zone, first northwards and then southwards, explains the wide areas of introgression at either side of the current position of the contact zones.     

Testing an hypothesis of hybrid zone movement for toads in France

Hybrid zone movement may result in substantial unidirectional introgression of selectively neutral material from the local to the advancing species, leaving a genetic footprint. This genetic footprint is represented by a trail of asymmetric tails and displaced cline centres in the wake of the moving hybrid zone. A peak of admixture linkage disequilibrium is predicted to exist ahead of the centre of the moving hybrid zone. We test these predictions of the movement hypothesis in a hybrid zone between common (Bufo bufo) and spined toads (B. spinosus), using 31 nuclear and one mtDNA SNPs along a transect in the northwest of France. Average effective selection in Bufo hybrids is low and clines vary in shape and centre. A weak pattern of asymmetric introgression is inferred from cline discordance of seven nuclear markers. The dominant direction of gene flow is from B. spinosus to B. bufo and is in support of southward movement of the hybrid zone. Conversely, a peak of admixture linkage disequilibrium north of the hybrid zone suggests northward movement. These contrasting results can be explained by reproductive isolation of the B. spinosus and B. bufo gene pools at the southern (B. spinosus) side of the hybrid zone. The joint occurrence of asymmetric introgression and admixture linkage disequilibrium can also be explained by the combination of low dispersal and random genetic drift due to low effective population sizes.     

Polymerase chain reaction primers for microsatellite loci in the Common Toad Bufo bufo

The Common Toad Bufo bufo is a wide‐ranging species, with a distribution encompassing much of Europe and Asia. Few molecular studies have been undertaken on this species, and only one polymorphic microsatellite locus has been identified. Therefore, little is known about the genetic variability within and between B. bufo populations. The value of such information is essential for monitoring the species and its environment. This paper reports the characterization of 15 B. bufo microsatellite loci using individuals from a Spanish and UK sample.     

Molecular phylogeny and biogeography of Malagasy frogs of the genus Gephyromantis

The genus Gephyromantis is a clade within the Malagasy-Comoroan family Mantellidae composed of rainforest frogs that live and breed to varying degrees independently from water. Based on DNA sequences of five mitochondrial and five nuclear genes we inferred the phylogeny of these frogs with full taxon coverage at the species level. Our preferred consensus tree from a partitioned Bayesian analysis of 5843 base pairs of 51 nominal and candidate species supports various major clades within the genus although the basal relationships among these remain unresolved. The data provide strong evidence for the monophyly of the subgenera Gephyromantis (after exclusion of Gephyromantis klemmeri), Laurentomantis, Vatomantis, and Phylacomantis. Species assigned to the subgenus Duboimantis belong to two strongly supported clades of uncertain relationships. G. klemmeri, previously in the subgenus Gephyromantis, was placed with high support sister to the Laurentomantis clade, and the Laurentomantis + G. klemmeri clade was sister to Vatomantis. A reconstruction of ancestral distribution areas indicates a diversification of several subgenera in the northern biogeographic regions of Madagascar and the dispersal out of northern Madagascar for several clades.     

Phylogeography and genetic relationships of North African <Emphasis Type="Italic">Bufo mauritanicus</Emphasis> Schlegel, 1841 estimated from mitochondrial DNA sequences

Genetic variation within the North African toad Bufo mauritanicus was estimated by sequencing partial 12S rRNA and 16S rRNA mitochondrial regions from widespread populations in Morocco, Algeria and Tunisia. Unlike many other wide ranging species from this area, B. mauritanicus demonstrated very low levels of intraspecific variation. The minimal intraspecific genetic variation may be due to a relatively recent, possibly post-glacial, expansion into its current range. Further phylogeographic studies of other North African species are needed to assess if this is a common biogeographical phenomenon. Phylogenetic analyses support immunological data that B. mauritanicus is part of a clade of predominantly sub-Saharan Bufo, recently assigned to a new genus Amietophrynus. Two different lineages within this clade, B. mauritanicus and the B. pardalis group, appear to have reverted from 20 chromosomes to the more typical 22 chromosomes found in most other Bufonids. However, the alternative hypothesis that the Bufo species with 20 chromosomes form a monophyletic lineage cannot be rejected.     

A phylogenetic analysis of species in the Bufo crucifer group (Anura: Bufonidae), based on indolealkylamines and proteins from skin secretions

This research tested the utility of two classes of skin secretion compounds to the phylogeny of the Bufo crucifer group. Skin secretions from specimens of nine populations of B. crucifer group were obtained and submitted to qualitative analysis. We observed a clear difference in the composition of the skin secretion molecules obtained from the species of Bufo studied. Fifty-nine molecules, 16 indolealkylamines and 43 proteins, were used as characters, and 39 of these were parsimonious informative. The tree topology of the skin secretion combined data showed areas of congruence and conflict when compared to an mtDNA phylogeny of the B. crucifer group. We used the Templeton test to evaluate the heterogeneity between the skin secretion and mtDNA data. Although not recommended, we performed a combined analysis with the two partitions. The skin secretion characters from the species of Bufo studied have phylogenetic signal. These data are indicative, at least as a preliminary study, of the phylogenetic relationships among the B. crucifer group taxa.     

Phylogeny of east asian bufonids inferred from mitochondrial DNA sequences (Anura: amphibia)

We investigated the relationships of Asian bufonids using partial sequences of mitochondrial DNA genes. Twenty-six samples representing 14 species of Bufo from China and Vietnam and 2 species of Torrentophryne from China were examined. Three samples of Bufo viridis from Armenia and Georgia were also sequenced to make a comparison to its sibling tetraploid species B. danatensis. Bufo americanus, from Canada, was used as the outgroup. Sequences from the 12S ribosomal RNA, 16S ribosomal RNA, cytochrome b, and the control region were analyzed using parsimony. East Asian bufonids were grouped into two major clades. One clade included B. andrewsi, B. bankorensis, B. gargarizans, B. tibetanus, B. tuberculatus, its sister clade B. cryptotympanicus, and the 2 species of Torrentophryne. The second clade consisted of B. galeatus, B. himalayanus, B. melanostictus, and a new species from Vietnam. The placement of three taxa (B. raddei, B. viridis, and its sister species, B. danatensis) was problematic. The genus Torrentophryne should be synonymized with Bufo to remove paraphyly. Because B. raddei does not belong to the clade that includes B. viridis and B. danatensis, it was removed from the viridis species group. The species status of B. bankorensis from Taiwan is evaluated.     

First report of a green toad (Bufo viridis sensu lato) in the Early Pleistocene of Spain: Palaeobiogeographical and palaeoecological implications

For the first time unequivocal fossil remains of a green toad (Bufo viridis s.l.) are described in the Iberian Peninsula. The fossils come from the Cueva Victoria site, a late Early Pleistocene (ca. 1.1–1.2 Ma) karstic filling in semi-arid southeastern Spain (Murcia region). By extension, other remains from two other Early Pleistocene Spanish localities, Barranco León D (ca. 1.3 Ma) and Almenara-Casablanca 3 (ca. 1.1 Ma), are cautiously attributed to the group B. viridis. The B. viridis group was previously reported with some uncertainty to the west of its current distribution area in Western Europe (Spain and France) in the Pliocene (Bufo cf. viridis) and less probably in the Early Miocene (Bufo aff. viridis). Since no osteological differences have been established between the recently described extant species of B. viridis s.l. (e.g. Bufo balearicus, Bufo siculus, Bufo boulengeri, B. viridis sensu stricto and Bufo variabilis) no precise palaeobiogeographical relationships can be drawn for the Spanish fossils. However, the occurrence of a third species of bufonid toad during the Pleistocene in the South of the Iberian Peninsula raises some interesting ecological questions in relation to the local disappearance of the green toad, which can be hypothetically linked to the intensification of the Pleistocene glacial/interglacial climate dynamic or to probable competition with another toad, Bufo calamita.     

Bradyrhizobium canariense and Bradyrhizobium japonicum are the two dominant rhizobium species in root nodules of lupin and serradella plants growing in Europe

Forty three Bradyrhizobium strains isolated in Poland from root nodules of lupin species (Lupinus albus, L. angustifolius and L. luteus), and pink serradella (Ornithopus sativus) were examined based on phylogenetic analyses of three housekeeping (atpD, glnII and recA) and nodulation (nodA) gene sequences. Additionally, seven strains originating from root-nodules of yellow serradella (O. compressus) from Asinara Island (Italy) were included in this study. Phylogenetic trees revealed that 15 serradella strains, including all yellow serradella isolates, and six lupin strains grouped in Bradyrhizobium canariense (BC) clade, whereas eight strains from pink serradella and 15 lupin strains were assigned to Bradyrhizobium japonicum (BJ1). Apparently, these species are the two dominant groups in soils of central Europe, in the nodules of lupin and serradella plants. Only three strains belonged to other chromosomal lineages: one formed a cluster that was sister to B. canariense, one strain grouped outside the branch formed by B. japonicum super-group, and one strain occupied a distant position in the genus Bradyrhizobium, clustering with strains of the Rhodopseudomonas genus. All strains in nodulation nodA gene tree grouped in a cluster referred to as Clade II, which is in line with earlier data on this clade dominance among Bradyrhizobium strains in Europe. The nodA tree revealed four well-supported subgroups within Clade II (II.1-II.4). Interestingly, all B. canariense strains clustered in subgroup II.1 whereas B. japonicum strains dominated subgroups II.2-II.4.     

A new species of ferret-badger, Genus Melogale, from Vietnam

Ferret-badgers, genus Melogale, are distributed in the Indochinese region, Java, Bali and NE-Borneo. There are currently four species described each having very similar phenotypes. In March 2005, a living ferret-badger of a different phenotype was confiscated by rangers from Cuc Phuong National Park, Vietnam. This individual died and the carcass was not preserved. In January 2006, a newly deceased individual of the same phenotype was found at the Endangered Primate Rescue Center, Cuc Phuong National Park. Due to several different characteristics these individuals vary greatly from the current species. Thus, we describe an additional species, M. cucphuongensis sp. nov. from northern Vietnam, which occurs sympatrically with M. moschata and M. personata, but differs from both species clearly in skull morphology and other features.Based on a 423 bp-long fragment of the mitochondrial cytochrome b gene, M. cucphuongensis sp. nov. is a member of the genus Melogale and represents a sister lineage to a clade consisting of M. personata and M. moschata.     

Phylogenetic relationships in European Ceriporiopsis species inferred from nuclear and mitochondrial ribosomal DNA sequences

The aim of this work was to clarify taxonomy and examine evolutionary relationships within European Ceriporiopsis species using a combined analysis of the large subunit (nLSU) nuclear rRNA and small subunit (mtSSU) mitochondrial rRNA gene sequences. Data from the ITS region were applied to enhance the view of the phylogenetic relationships among different species. The studied samples grouped into four complex clades, suggesting that the genus Ceriporiopsis is polyphyletic. The generic type Ceriporiopsis gilvescens formed a separate group together with Ceriporiopsis guidella and Phlebia spp. in the phlebioid clade. In this clade, the closely related species Ceriporiopsis resinascens and Ceriporiopsis pseudogilvescens grouped together with Ceriporiopsis aneirina. C. resinascens and C. pseudogilvescens have identical LSU and SSU sequences but differ in ITS. Ceriporiopsis pannocincta also fell in the phlebioid clade, but showed closer proximity to Gloeoporus dichrous than to C. gilvescens or C. aneirina–C. pseudogilvescens–C. resinascens group. Another clade was composed of a Ceriporiopsis balaenae–Ceriporiopsis consobrina group and was found to be closely related to Antrodiella and Frantisekia, with the overall clade highly reminiscent of the residual polyporoid clade. The monotypic genus Pouzaroporia, erected in the past for Ceriporiopsis subrufa due to its remarkable morphological differences, also fell within the residual polyporoid clade. Ceriporiopsis subvermispora held an isolated position from the other species of the genus. Therefore, the previously proposed name Gelatoporia subvermispora has been adopted for this species. Physisporinus rivulosus appeared unrelated to two other European Physisporinus species. Moreover, Ceriporiopsis (=Skeletocutis) jelicii grouped in a separate clade, distinct from Ceriporiopsis species. Finally, the ITS data demonstrated the proximity of some Ceriporiopsis species (Ceriporiopsis portcrosensis and Ceriporiopsis subsphaerospora) to Skeletocutis amorpha.     

Endosymbiotic bacteria nodulating a new endemic lupine Lupinus mariae-josephi from alkaline soils in Eastern Spain represent a new lineage within the Bradyrhizobium genus

Lupinus mariae-josephi is a recently described endemic Lupinus species from a small area in Eastern Spain where it thrives in soils with active lime and high pH. The L. mariae-josephi root symbionts were shown to be very slow-growing bacteria with different phenotypic and symbiotic characteristics from those of Bradyrhizobium strains nodulating other Lupinus. Their phylogenetic status was examined by multilocus sequence analyses of four housekeeping genes (16S rRNA, glnII, recA, and atpD) and showed the existence of a distinct evolutionary lineage for L. mariae-josephi that also included Bradyrhizobium jicamae. Within this lineage, the tested isolates clustered in three different sub-groups that might correspond to novel sister Bradyrhizobium species. These core gene analyses consistently showed that all the endosymbiotic bacteria isolated from other Lupinus species of the Iberian Peninsula were related to strains of the B. canariense or B. japonicum lineages and were separate from the L. mariae-josephi isolates. Phylogenetic analysis based on nodC symbiotic gene sequences showed that L. mariae-josephi bacteria also constituted a new symbiotic lineage distant from those previously defined in the genus Bradyrhizobium. In contrast, the nodC genes of isolates from other Lupinus spp. from the Iberian Peninsula were again clearly related to the B. canariense and B. japonicum bv. genistearum lineages. Speciation of L. mariae-josephi bradyrhizobia may result from the colonization of a singular habitat by their unique legume host.     

Genetic diversity of root nodulating bacteria associated with Retama sphaerocarpa in sites with different soil and environmental conditions

The genetic diversity of root nodulating bacteria isolated from Retama sphaerocarpa was studied using BOX-A1R PCR and phylogenetic analysis of the 16S rRNA region, as well as the housekeeping genes atpD, glnII and recA. A total of 193 isolates were obtained from eight different sites with different soil and environmental conditions in the Iberian Peninsula. These isolates corresponded to 31 different strains that successfully nodulated R. sphaerocarpa seedlings in reinoculation trials. About one-third of the strains clustered with B. canariense or B. cytisi within Bradyrhizobium group I. The remaining strains clustered with B. elkanii/B. pachyrhizi within Bradyrhizobium group II or in separate clades that could represent new lineages. Based on the 16S rRNA and combined atpD + glnII + recA sequences, two to three lineages of root nodulating bacteria were found at each sampling site, except for Collado Garcia where five species were detected. B. canariense and B. elkanii/B. pachyrhizi were the most abundant species, whereas the least abundant were those related to B. retamae and a putative new lineage. B. canariense was found only in soils with neutral and acid pH, whereas B. retamae was the dominant species in alkaline soils.     

The phylogeny of Ephemeroptera in Pterygota revealed by the mitochondrial genome of Siphluriscus chinensis (Hexapoda: Insecta)

The mayfly species Siphluriscus chinensis (Siphluriscidae) has valuable structures useful for phylogeny reconstruction, given its putative basal position within the Ephemeroptera. Here its nearly complete mitochondrial genome is sequenced. We built phylogenetic trees through multiple analytical strategies with some other insect mitogenomes. Structurally, the obtained mitochondrial genome of S. chinensis is 16,616 bp in length, ¹ containing 37 genes and an extra trnK-like (trnK2 (AAA)) gene. The 12 PCGs start with typical ATN codons, except the nad1 gene which starts with an unnormalized TTG. Like other known mayfly mitogenomes, the strand bias has negative AT-skew and negative GC-skew. Phylogenetically, our topologies suggest that Odonata is the basally diverged clade in Pterygota; Ephemeroptera is the sister group of the Neoptera; and S. chinensis is indeed the most basal mayfly branch.     

Phylogenetic selection of target species in Amaryllidaceae tribe Haemantheae for acetylcholinesterase inhibition and affinity to the serotonin reuptake transport protein

We present phylogenetic analyses of 37 taxa of Amaryllidaceae, tribe Haemantheae and Amaryllis belladonna L. as an outgroup, in order to provide a phylogenetic framework for the selection of candidate plants for lead discoveries in relation to Alzheimer's disease and depression. DNA sequences from the nuclear ribosomal internal transcribed spacer (ITS) and the plastid trnL-F regions were used. Maximum parsimony analyses provide increased support for the sister relationship of Haemanthus and Scadoxus. Within Haemanthus, a well supported clade (89% BS) corresponds to a summer rainfall group (mainly Eastern Cape) with white-pale pink flowers. A second strongly supported clade (100% BS) corresponds to a winter rainfall group (mainly Western Cape) with red-pale pink flowers. Haemanthus montanus, which is from the summer rainfall region, is sister to the winter rainfall group. Alkaloid profiles and bioactivity profiles were investigated for 16 taxa using gas chromatography-mass spectrometry (GC-MS) and assays measuring acetylcholinesterase (AChE) inhibition and affinity to the serotonin reuptake transport protein (SERT). No alkaloids were detected by GC-MS in extracts of the two species of Gethyllis included in the present study suggesting that Gethyllis (and possibly Apodolirion) species may not produce the alkaloids characteristic for the family. AChE inhibitory activity was found in all investigated clades except the Apodolirion-Gethyllis clade, which can be explained by the apparent lack of alkaloids in this clade. In spite of infra-specific variability of alkaloid profiles observed, dose-dependent SERT activity appears to be pronounced and restricted to the genus Haemanthus within tribe Haemantheae. Three of eight Haemanthus species tested had IC₅₀ < 10 μg/ml. Two of the most active extracts in the present study contained primarily montanine type alkaloids which have not been tested for SERT affinity previously. Simultaneous evaluation of bioactivity and alkaloid profiles in a phylogenetic framework can potentially be used to select candidate species for phytotherapy and drug discovery.     

Phylogenetic Framework for Trema (Celtidaceae)

We used ITS and trnL sequence data, analyzed separately and combined by MP, to explore species relationships and concepts in Trema (Celtidaceae), a pantropical genus of pioneer trees. Whether Trema is monophyletic or includes Parasponia is still unresolved. Three clades within Trema received moderate to high support, one from the New World and two from the Old World, but their relationships were not resolved. In the New World, specimens of T. micrantha formed two groups consistent with endocarp morphology. Group I, with smaller brown endocarps, is a highly supported clade sister to T. lamarckiana. Group II, with larger black endocarps, is poorly resolved with several subclades, including the highly supported T. integerrima clade. Both Old World clades contain Asian and African species, with three or more species in each region. Trema orientalis is not monophyletic: specimens from Africa formed a highly supported clade sister to T. africana, while those from Asia were sister to T. aspera from Australia.     

Historical biogeography and character evolution of Cistaceae (Malvales) based on analysis of plastid rbcL and trnL-trnF sequences

Cistaceae consist of eight genera and about 180 species. Some taxonomic limits and relationships within the family remain unresolved when relying exclusively on morphological data. In the present study, a phylogeny was reconstructed and divergence times were estimated for 47 species representing various groups in Cistaceae and using coding (rbcL) and spacer (trnL-trnF) sequences of plastid DNA. The firm set of morphological synapomorphies that indicates the monophyly of the family is supported by both Bayesian and parsimony analyses. Five major lineages can be distinguished within the Cistaceae: (1) an early-diverging lineage containing Fumana species; (2) the New World Lechea clade; (3) the Helianthemum s.l. clade, containing two sister groups, one of species from the New World (Crocanthemum, Hudsonia) and the other with species from the Old World (Helianthemum s. str.); (4) the Tuberaria clade; and (5) a cohesive complex consisting of Halimium and Cistus species. Evolutionary shifts in 12 key characters of Cistaceae are inferred based on the most plausible phylogenetic hypothesis. Reconstructing the evolution of ovule position supports anatropous ovules as the ancestral condition within the Cistaceae, which is currently found only in Fumana. The Cistus-Halimium assemblage is consistently obtained as a natural clade and further supported by a cytological synapomorphy (chromosome number n=9). Optimisation of ancestral distribution areas and estimates of divergence times reveal an early divergence (10.17-18.51 Ma) of the Mediterranean-European genera, which may be related to subtropical vegetation, as complemented by paleobotanical data. In addition, the occurrence of multiple, independent migration events from the Old World to America between the Middle Miocene (8.44-14.7 Ma; Lechea) and the Upper Miocene (5.15-9.20 Ma; Crocanthemum/Hudsonia), and to the Canary Islands in the Pleistocene is inferred. We argue that the Mediterranean basin has been the main centre of differentiation of Cistaceae.     

Differential sensitivity of four Lemnaceae species to zinc sulphate

Lemnaceae are currently the only freshwater plants required for regulatory toxicity testing of pesticides and other chemicals. Toxicological protocols allow for the use of different Lemnaceae species in tests. However, few studies have compared the relative sensitivity of individual duckweed species. Zinc is an essential plant nutrient but is also a common pollutant in aquatic environments and elevated levels are phytotoxic. This study shows that four species of Lemnaceae differ in their relative sensitivities to zinc sulphate, a commonly used reference chemical. Comparative zinc sensitivity, in order, from most tolerant to most sensitive was: Landoltia punctata > Lemna minor > Wolffia brasiliensis > Lemna gibba. Zinc sensitivity was also endpoint dependant. EC₅₀ values typically increased in order of: specific biomass growth rate < specific frond number growth rate < chlorophyll absorbance. However, specific frond number growth rate was the most sensitive endpoint for L. punctata. Unlike the other species, L. punctata displayed no significant colony disintegration. Lemna species and L. punctata appear to be employing distinct response strategies when exposed to zinc. L. gibba and L. minor produce and release young, single fronds which are severely affected by zinc. In contrast, L. punctata produces fewer fronds, which are not released and form large colonies of high biomass that are relatively zinc tolerant.     

Retama species growing in different ecological–climatic areas of northeastern Algeria have a narrow range of rhizobia that form a novel phylogenetic clade within the Bradyrhizobium genus

Sixty-seven isolates were isolated from nodules collected on roots of Mediterranean shrubby legumes Retama raetam and Retama sphaerocarpa growing in seven ecological–climatic areas of northeastern Algeria. Genetic diversity of the Retama isolates was analyzed based on genotyping by restriction fragment length polymorphism of PCR-amplified fragments of the 16S rRNA gene, the intergenic spacer (IGS) region between the 16S and 23S rRNA genes (IGS), and the symbiotic genes nifH and nodC. Eleven haplotypes assigned to the Bradyrhizobium genus were identified. Significant biogeographical differentiation of the rhizobial populations was found, but one haplotype was predominant and conserved across the sites. All isolates were able to cross-nodulate the two Retama species. Accordingly, no significant genetic differentiation of the rhizobial populations was found in relation to the host species of origin. Sequence analysis of the 16S rRNA gene grouped the isolates with Bradyrhizobium elkanii, but sequence analyses of IGS, the housekeeping genes (dnaK, glnII, recA), nifH, and nodC yielded convergent results showing that the Retama nodule isolates from the northeast of Algeria formed a single evolutionary lineage, which was well differentiated from the currently named species or well-delineated unnamed genospecies of bradyrhizobia. Therefore, this study showed that the Retama species native to northeastern Algeria were associated with a specific clade of bradyrhizobia. The Retama isolates formed three sub-groups based on IGS and housekeeping gene phylogenies, which might form three sister species within a novel bradyrhizobial clade.