Japanese Hollies: Intolerant Hosts of Meloidogyne arenaria in Microplots

Japanese hollies were itttolerant of Meloidogyne arenaria in field microplot experiments. Ilex crenata var. rotundifolia was relatively more tolerant than I. crenata var. convexa or I. crenata var. helleri. When M. arenaria was added at various itfitial population densities to soil containing plants of "Helleri," "Convexa," and "Rotundifolia," respectively, 91, 75, and 25% were killed by the end of the third growing season. No control plants died during the same period. Initial numbers of M. arenaria larvae and eggs were the only population densities that were correlated (negatively), regardless of cultivar, with plant growth over the three growing seasons. A linear relation was found for initial density of M. arenaria and growth of I. crenata rotundifolia. Increasing nematode density by 10-fold suppressed the growth of this cultivar by 23%.     

Axenizing and Culturing Endomigratory Plant-Parasitic Nematodes using Pluronic F127, Including its Effects,on Population Dynamics of Pratylenchus penetrans

A non-chemical technique for surface sterilizing plant-parasitic nematodes for aseptic cultures is described. The method is most applicable to nematodes with active migratory infective stages and requires only a few starting specimens. Rate of achieving a primary aseptic culture with the technique ranged from 60%-100% depending on the conditions of the specimens collected for culturing. Aseptic cultures of species of Meloidogyne, Rotylenchuluz, Pratylenchus, and Radopholus initiated with the method remained contamination-free after 12 months of maintenance in tomato root explant or alfalfa callus cultures. Further studies of Pluronic F127, a polyol gel medium employed in the technique to confine the spread of contaminating bacteria or fungi associated with the nematodes, showed that the polyol gel was a suitable support medium for culturing corn root explant, alfalfa callus tissues, and consequently Pratylenchus species including P. agilis, P. brachyurus, P. scribneri, and P. penetrans. During the course of 10 months, P. penetrans reared in polyol-base medium followed a standard biological growth curve, multiplied to a higher population density, maintained a similar female-to-male ratio, and possessed a similar tendency to reside inside or outside host tissues as did P. penetrans reared in agar-base medium. The percentages of P. penetrans juveniles in the sub-populations residing outside or inside the host tissues reared in polyol-base medium also were similar to and fluctuated temporally in like manner as those reared in agar-base medium. Members of these sub-populations from the polyol- or agar-base were equally infective and reproductive after 9 months of culturing.     

Downregulation of net phosphorus-uptake capacity is inversely related to leaf phosphorus-resorption proficiency in four species from a phosphorus-impoverished environment

Background and Aims

Previous research has suggested a trade-off between the capacity of plants to downregulate their phosphorus (P) uptake capacity and their efficiency of P resorption from senescent leaves in species from P-impoverished environments.

Methods

To investigate this further, four Australian native species (Banksia attenuata, B. menziesii, Acacia truncata and A. xanthina) were grown in a greenhouse in nutrient solutions at a range of P concentrations [P]. Acacia plants received between 0 and 500 µm P; Banksia plants received between 0 and 10 µm P, to avoid major P-toxicity symptoms in these highly P-sensitive species.

Key Results

For both Acacia species, the net P-uptake rates measured at 10 µm P decreased steadily with increasing P supply during growth. In contrast, in B. attenuata, the net rate of P uptake from a solution with 10 µm P increased linearly with increasing P supply during growth. The P-uptake rate of B. menziesii showed no significant response to P supply in the growing medium. Leaf [P] of the four species supported this finding, with A. truncata and A. xanthina showing an increase up to a saturation value of 19 and 21 mg P g⁻¹ leaf dry mass, respectively (at 500 µm P), whereas B. attenuata and B. menziesii both exhibited a linear increase in leaf [P], reaching 10 and 13 mg P g⁻¹ leaf dry mass, respectively, without approaching a saturation point. The Banksia plants grown at 10 µm P showed mild symptoms of P toxicity, i.e. yellow spots on some leaves and drying and curling of the tips of the leaves. Leaf P-resorption efficiency was 69 % (B. attenuata), 73 % (B. menziesii), 34 % (A. truncata) and 36 % (A. xanthina). The P-resorption proficiency values were 0·08 mg P g⁻¹ leaf dry mass (B. attenuata and B. menziesii), 0·32 mg P g⁻¹ leaf dry mass (A. truncata) and 0·36 mg P g⁻¹ leaf dry mass (A. xanthina). Combining the present results with additional information on P-remobilization efficiency and the capacity to downregulate P-uptake capacity for two other Australian woody species, we found a strong negative correlation between these traits.

Conclusions

It is concluded that species that are adapted to extremely P-impoverished soils, such as many south-western Australian Proteaceae species, have developed extremely high P-resorption efficiencies, but lost their capacity to downregulate their P-uptake mechanisms. The results support the hypothesis that the ability to resorb P from senescing leaves is inversely related to the capacity to downregulate net P uptake, possibly because constitutive synthesis of P transporters is a prerequisite for proficient P remobilization from senescing tissues.     

A Classification of Tylenchulus semipenetrans Biotypes

The presence of two biotypes of the citrus nematode (Tylenchulus semipenetrans) in Italian citrus and olive orchards has been confirmed by comparing host specificity. Host reaction to California biotypes C1 and C3 and to three populations from Arizona, Texas, and Florida indicates that of these five United States biotypes, all except C3 consistently fit biotype C1. These findings, and the results of host-range studies in other countries, show that four biotypes of T. semipenetrans are distributed worldwide: the "Poncirus biotype," the "Citrus biotype," the "Mediterranean biotype," and the "Grass biotype."     

Effects of phosphorus supply on growth, phosphate concentration and cluster-root formation in three Lupinus species

Background and Aims

In some lupin species, phosphate deficiency induces cluster-root formation, which enhances P uptake by increasing root surface area and, more importantly, the release of root exudates which enhances P availability.

Methods

Three species of Lupinus, L. albus, L. atlanticus and L. micranthus, with inherently different relative growth rates were cultivated under hydroponics in a greenhouse at four phosphate concentrations (1, 10, 50 and 150 µm) to compare the role of internal P in regulating cluster-root formation.

Key Results

The highest growth rate was observed in L. atlanticus, followed by L. albus and L. micranthus. At 1 µm P, cluster-root formation was markedly induced in all three species. The highest P uptake and accumulation was observed in L. micranthus, followed by L. atlanticus and then L. albus. Inhibition of cluster-root formation was severe at 10 µm P in L. atlanticus, but occurred stepwise with increasing P concentration in the root medium in L. albus.

Conclusions

In L. atlanticus and L. albus cluster-root formation was suppressed by P treatments above 10 µm, indicating a P-inducible regulating system for cluster-root formation, as expected. By contrast, production of cluster roots in L. micranthus, in spite of a high internal P concentration, indicated a lower sensitivity to P status, which allowed P-toxicity symptoms to develop.     

Unidirectional hybridization and reproductive barriers between two heterostylous primrose species in north-west Yunnan,China

Background and Aims

Heteromorphy in flowers has a profound effect on breeding patterns within a species, but little is known about how it affects reproductive barriers between species. The heterostylous genus Primula is very diverse in the Himalaya region, but hybrids there have been little researched. This study examines in detail a natural hybrid zone between P. beesiana and P. bulleyana.

Methods

Chloroplast sequencing, AFLP (amplified fragment length polymorphism) markers and morphological comparisons were employed to characterize putative hybrids in the field, using synthetic F₁s from hand pollination as controls. Pollinator visits to parent species and hybrids were observed in the field. Hand pollinations were conducted to compare pollen tube growth, seed production and seed viability for crosses involving different morphs, species and directions of crossing.

Key Results

Molecular data revealed all hybrid derivatives examined to be backcrosses of first or later generations towards P. bulleyana: all had the chloroplast DNA (cpDNA) of this species. Some individuals had morphological traits suggesting they were hybrids, but they were genetically similar to P. bulleyana; they might have been advanced generation backcrosses. Viable F₁s could not be produced with P. bulleyana pollen on P. beesiana females, irrespective of the flower morphs used. Within-morph crosses for each species had very low (<10 %) seed viability, whereas crosses between pin P. bulleyana (female) and pin P. beesiana had a higher seed viability of 30 %. Thus genetic incompatibility mechanisms back up mechanical barriers to within-morph crosses in each species, but are not the same between the two species. The two species share their main pollinators, and pollinators were observed to fly between P. bulleyana and hybrids, suggesting that pollinator behaviour may not be an important isolating factor.

Conclusions

Hybridization is strongly asymmetric, with P. bulleyana the only possible mother and all detected hybrids being backcrosses in this direction. Partial ecological isolation and inhibition of heterospecific pollen, and possibly complete barriers to F₁ formation on P. beesiana, may be enough to make F₁ formation very rare in these species. Therefore, with no F₁ detected, this hybrid zone may have a finite life span as successive generations become more similar to P. bulleyana.     

Plastid DNA sequencing and nuclear SNP genotyping help resolve the puzzle of central American Platanus

Background and Aims

Recent research on the history of Platanus reveals that hybridization phenomena occurred in the central American species. This study has two goals: to help resolve the evolutive puzzle of central American Platanus, and to test the potential of real-time polymerase chain reaction (PCR) for detecting ancient hybridization.

Methods

Sequencing of a uniparental plastid DNA marker [psbA-trnH^(GUG) intergenic spacer] and qualitative and quantitative single nucleotide polymorphism (SNP) genotyping of biparental nuclear ribosomal DNA (nrDNA) markers [LEAFY intron 2 (LFY-i2) and internal transcribed spacer 2 (ITS2)] were used.

Key Results

Based on the SNP genotyping results, several Platanus accessions show the presence of hybridization/introgression, including some accessions of P. rzedowskii and of P. mexicana var. interior and one of P. mexicana var. mexicana from Oaxaca (= P. oaxacana). Based on haplotype analyses of the psbA-trnH spacer, five haplotypes were detected. The most common of these is present in taxa belonging to P. orientalis, P. racemosa sensu lato, some accessions of P. occidentalis sensu stricto (s.s.) from Texas, P. occidentalis var. palmeri, P. mexicana s.s. and P. rzedowskii. This is highly relevant to genetic relationships with the haplotypes present in P. occidentalis s.s. and P. mexicana var. interior.

Conclusions

Hybridization and introgression events between lineages ancestral to modern central and eastern North American Platanus species occurred. Plastid haplotypes and qualitative and quantitative SNP genotyping provide information critical for understanding the complex history of Mexican Platanus. Compared with the usual molecular techniques of sub-cloning, sequencing and genotyping, real-time PCR assay is a quick and sensitive technique for analysing complex evolutionary patterns.     

Changes in activity of metabolic and regulatory pathways during germination of S. coelicolor

Background

Bacterial spore germination is a developmental process during which all required metabolic pathways are restored to transfer cells from their dormant state into vegetative growth. Streptomyces are soil dwelling filamentous bacteria with complex life cycle, studied mostly for they ability to synthesize secondary metabolites including antibiotics.

Results

Here, we present a systematic approach that analyzes gene expression data obtained from 13 time points taken over 5.5 h of Streptomyces germination. Genes whose expression was significantly enhanced/diminished during the time-course were identified, and classified to metabolic and regulatory pathways. The classification into metabolic pathways revealed timing of the activation of specific pathways during the course of germination. The analysis also identified remarkable changes in the expression of specific sigma factors over the course of germination. Based on our knowledge of the targets of these factors, we speculate on their possible roles during germination. Among the factors whose expression was enhanced during the initial part of germination, SigE is though to manage cell wall reconstruction, SigR controls protein re-aggregation, and others (SigH, SigB, SigI, SigJ) control osmotic and oxidative stress responses.

Conclusions

From the results, we conclude that most of the metabolic pathway mRNAs required for the initial phases of germination were synthesized during the sporulation process and stably conserved in the spore. After rehydration in growth medium, the stored mRNAs are being degraded and resynthesized during first hour. From the analysis of sigma factors we conclude that conditions favoring germination evoke stress-like cell responses.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1173) contains supplementary material, which is available to authorized users.     

Host responses and metabolic profiles of wood components in Dutch elm hybrids with a contrasting tolerance to Dutch elm disease

Background and Aims

Changes occurring in the macromolecular traits of cell wall components in elm wood following attack by Ophiostoma novo-ulmi, the causative agent of Dutch elm disease (DED), are poorly understood. The purpose of this study was to compare host responses and the metabolic profiles of wood components for two Dutch elm (Ulmus) hybrids, ‘Groeneveld’ (a susceptible clone) and ‘Dodoens’ (a tolerant clone), that have contrasting survival strategies upon infection with the current prevalent strain of DED.

Methods

Ten-year-old plants of the hybrid elms were inoculated with O. novo-ulmi ssp. americana × novo-ulmi. Measurements were made of the content of main cell wall components and extractives, lignin monomer composition, macromolecular traits of cellulose and neutral saccharide composition.

Key Results

Upon infection, medium molecular weight macromolecules of cellulose were degraded in both the susceptible and tolerant elm hybrids, resulting in the occurrence of secondary cell wall ruptures and cracks in the vessels, but rarely in the fibres. The ¹³C nuclear magnetic resonance spectra revealed that loss of crystalline and non-crystalline cellulose regions occurred in parallel. The rate of cellulose degradation was influenced by the syringyl:guaiacyl ratio in lignin. Both hybrids commonly responded to the medium molecular weight cellulose degradation with the biosynthesis of high molecular weight macromolecules of cellulose, resulting in a significant increase in values for the degree of polymerization and polydispersity. Other responses of the hybrids included an increase in lignin content, a decrease in relative proportions of d-glucose, and an increase in proportions of d-xylose. Differential responses between the hybrids were found in the syringyl:guaiacyl ratio in lignin.

Conclusions

In susceptible ‘Groeneveld’ plants, syringyl-rich lignin provided a far greater degree of protection from cellulose degradation than in ‘Dodoens’, but only guaiacyl-rich lignin in ‘Dodoens’ plants was involved in successful defence against the fungus. This finding was confirmed by the associations of vanillin and vanillic acid with the DED-tolerant ‘Dodoens’ plants in a multivariate analysis of wood traits.     

Direct and indirect influences of arbuscular mycorrhizal fungi on phosphorus uptake by two root hemiparasitic Pedicularis species: do the fungal partners matter at low colonization levels?

Background and Aims

Because most parasitic plants do not form mycorrhizal associations, the nutritional roles of arbuscular mycorrhizal (AM) fungi in them have hardly been tested. Some facultative root hemiparasitic Pedicularis species form AM associations and hence are ideal for testing both direct and indirect effects of AM fungi on their nutrient acquisition. The aim of this study was to test the influence of AM inoculation on phosphorus (P) uptake by Pedicularis rex and P. tricolor.

Methods

³²P labelling was used in compartmented pots to assess the contribution of the AM pathway and the influence of AM inoculation on P uptake from a host plant into the root hemiparasites. Laboratory isolates of fungal species (Glomus mosseae and G. intraradices) and the host species (Hordeum vulgare ‘Fleet’) to which the two Pedicularis species showed obvious responses in haustorium formation and growth in previous studies were used.

Key Results

The AM colonization of both Pedicularis spp. was low (<15 % root length) and only a very small proportion of total plant P (<1 %) was delivered from the soil via the AM fungus. In a separate experiment, inoculation with AM fungi strongly interfered with P acquisition by both Pedicularis species from their host barley, almost certainly because the numbers of haustoria formed by the parasite were significantly reduced in AM plants.

Conclusions

Roles of AM fungi in nutrient acquisition by root parasitic plants were quantitatively demonstrated for the first time. Evidence was obtained for a novel mechanism of preventing root parasitic plants from overexploiting host resources through AM fungal-induced suppression of the absorptive structures in the parasites.     

Description of a New and Two Known Species of the Free-living Nematode Paroigolaimella Paramonov, 1952 (Diplogastridae) from India

This paper describes a new and two known species of Paroigolaimella collected from India. Paroigolaimella helalii n. sp. is characterized by having conspicuous sexual dimorphism in the stoma and pharynx, ovaries with a sphincter separating the mature oocyte from developing ones, a vagina leading to a strong ovijector, a pore-like vulva with cuticular flap; males with slender strongly arcuate spicules with dilated capitula; the gubernaculum slender with expanded plate-like distal end and nine pairs of genital papillae, and four to five pairs of copulatory muscle bands. P. coprophila (Sudhaus and Rehfeld, 1990) Sudhaus and Fürst von Lieven, 2003 collected from leaf litter from a farmyard has been redescribed with reassessment of its distinguishing characters from P. bernensis. P. bodamica (Micoletzky, 1922) n. comb. has been described and its status has been discussed with context to P. bernensis.     

Pan-genome dynamics of Pseudomonas gene complements enriched across hexachlorocyclohexane dumpsite

Background

Phylogenetic heterogeneity across Pseudomonas genus is complemented by its diverse genome architecture enriched by accessory genetic elements (plasmids, transposons, and integrons) conferring resistance across this genus. Here, we sequenced a stress tolerant genotype i.e. Pseudomonas sp. strain RL isolated from a hexachlorocyclohexane (HCH) contaminated pond (45 mg of total HCH g⁻¹ sediment) and further compared its gene repertoire with 17 reference ecotypes belonging to P. stutzeri, P. mendocina, P. aeruginosa, P. psychrotolerans and P. denitrificans, representing metabolically diverse ecosystems (i.e. marine, clinical, and soil/sludge). Metagenomic data from HCH contaminated pond sediment and similar HCH contaminated sites were further used to analyze the pan-genome dynamics of Pseudomonas genotypes enriched across increasing HCH gradient.

Results

Although strain RL demonstrated clear species demarcation (ANI ≤ 80.03%) from the rest of its phylogenetic relatives, it was found to be closest to P. stutzeri clade which was further complemented functionally. Comparative functional analysis elucidated strain specific enrichment of metabolic pathways like α-linoleic acid degradation and carbazole degradation in Pseudomonas sp. strain RL and P. stutzeri XLDN-R, respectively. Composition based methods (%codon bias and %G + C difference) further highlighted the significance of horizontal gene transfer (HGT) in evolution of nitrogen metabolism, two-component system (TCS) and methionine metabolism across the Pseudomonas genomes used in this study. An intact mobile class-I integron (3,552 bp) with a captured gene cassette encoding for dihydrofolate reductase (dhfra1) was detected in strain RL, distinctly demarcated from other integron harboring species (i.e. P. aeruginosa, P. stutzeri, and P. putida). Mobility of this integron was confirmed by its association with Tnp21-like transposon (95% identity) suggesting stress specific mobilization across HCH contaminated sites. Metagenomics data from pond sediment and recently surveyed HCH adulterated soils revealed the in situ enrichment of integron associated transposase gene (TnpA6100) across increasing HCH contamination (0.7 to 450 mg HCH g⁻¹ of soil).

Conclusions

Unlocking the potential of comparative genomics supplemented with metagenomics, we have attempted to resolve the environment and strain specific demarcations across 18 Pseudomonas gene complements. Pan-genome analyses of these strains indicate at astoundingly diverse metabolic strategies and provide genetic basis for the cosmopolitan existence of this taxon.

Electronic supplementary material

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

Population genetic evidence for speciation pattern and gene flow between Picea wilsonii,P. morrisonicola and P. neoveitchii

Background and Aims

Genetic drift due to geographical isolation, gene flow and mutation rates together make it difficult to determine the evolutionary relationships of present-day species. In this study, population genetic data were used to model and decipher interspecific relationships, speciation patterns and gene flow between three species of spruce with similar morphology, Picea wilsonii, P. neoveitchii and P. morrisonicola. Picea wilsonii and P. neoveitchii occur from central to north-west China, where they have overlapping distributions. Picea morrisonicola, however, is restricted solely to the island of Taiwan and is isolated from the other two species by a long distance.

Methods

Sequence variations were examined in 18 DNA fragments for 22 populations, including three fragments from the chloroplast (cp) genome, two from the mitochondrial (mt) genome and 13 from the nuclear genome.

Key Results

In both the cpDNA and the mtDNA, P. morrisonicola accumulated more species-specific mutations than the other two species. However, most nuclear haplotypes of P. morrisonicola were shared by P. wilsonii, or derived from the dominant haplotypes found in that species. Modelling of population genetic data supported the hypothesis that P. morrisonicola derived from P. wilsonii within the more recent past, most probably indicating progenitor–derivative speciation with a distinct bottleneck, although further gene flow from the progenitor to the derivative continued. In addition, the occurrence was detected of an obvious mtDNA introgression from P. neoveitchii to P. wilsonii despite their early divergence.

Conclusions

The extent of mutation, introgression and lineage sorting taking place during interspecific divergence and demographic changes in the three species had varied greatly between the three genomes. The findings highlight the complex evolutionary histories of these three Asian spruce species.     

Role of Nematodes and Soil-borne Fungi in Cotton Stunt

The nematodes, Pratylenchus brachyurus, Trichodorus christiei, and T. porosus and the soil-borne fungi, Rhizoctonia solani, Pythium debaryanum, P. irregulare, P. ultimum, and Fusarium spp. were the pathogens most frequently found in the roots and rhizosphere of field-grown cotton (Gossypium hirsutum) showing "stunt" symptoms. Field-plot application of the nematicide D-D (l,2-dichloropropane, 1,3-dichloropropene) at 373.4 liter/ha (40 gal/A) significantly increased plant growth and yield. A fungicidal mixture of Dexon (p-dimethylaminobenzenediazo sodium sulfonate at 23.5 kg/ha (2l lb/A) and Terraclor (pentachloronitrobenzene at 25.2 kg/ha (22.5 lb/A) was phytotoxic, but combined nematicide/fungicide treatments were not. Greenhouse temperature-tank experiments in soils from two locations showed significantly improved root and shoot growth following methyl bromide fumigation at both 25 C and 18 C and more severe "stunt" at the lower temperature.     

Transition State Analogues of Plasmodium falciparum and Human Orotate Phosphoribosyltransferases

The survival and proliferation of Plasmodium falciparum parasites and human cancer cells require de novo pyrimidine synthesis to supply RNA and DNA precursors. Orotate phosphoribosyltransferase (OPRT) is an indispensible component in this metabolic pathway and is a target for antimalarials and antitumor drugs. P. falciparum (Pf) and Homo sapiens (Hs) OPRTs are characterized by highly dissociative transition states with ribocation character. On the basis of the geometrical and electrostatic features of the PfOPRT and HsOPRT transition states, analogues were designed, synthesized, and tested as inhibitors. Iminoribitol mimics of the ribocation transition state in linkage to pyrimidine mimics using methylene or ethylene linkers gave dissociation constants (K_d) as low as 80 nm. Inhibitors with pyrrolidine groups as ribocation mimics displayed slightly weaker binding affinities for OPRTs. Interestingly, p-nitrophenyl riboside 5′-phosphate bound to OPRTs with K_d values near 40 nm. Analogues designed with a C5-pyrimidine carbon–carbon bond to ribocation mimics gave K_d values in the range of 80–500 nm. Acyclic inhibitors with achiral serinol groups as the ribocation mimics also displayed nanomolar inhibition against OPRTs. In comparison with the nucleoside derivatives, inhibition constants of their corresponding 5′-phosphorylated transition state analogues are largely unchanged, an unusual property for a nucleotide-binding site. In silico docking of the best inhibitor into the HsOPRT active site supported an extensive hydrogen bond network associated with the tight binding affinity. These OPRT transition state analogues identify crucial components of potent inhibitors targeting OPRT enzymes. Despite their tight binding to the targets, the inhibitors did not kill cultured P. falciparum.     

Phylogeography of Pinus subsection Australes in the Caribbean Basin

Background and Aims

Four species of Pinus subsection Australes occur in the Caribbean Basin: P. caribaea, P. cubensis, P. maestrensis and P. occidentalis. This study analyses the phylogeography of these species to assess possible colonization events from Central America to the islands and subsequent population expansions during glacial periods driven by both drier climate and larger emerged land areas.

Methods

Allele size data were obtained for plastid microsatellites for 314 individuals from 24 populations, covering the distribution range of subsection Australes in the Caribbean Basin.

Key Results

In total, 113 plastid haplotypes were identified. The highest genetic diversity was found in populations of P. caribaea. Overall, Caribbean Basin populations fit the isolation by distance model. Significant phylogeographical structure was found (R_ST = 0·671 > permuted R_ST = 0·101; P < 0·0001). The haplotype network and a Bayesian analysis of population structure (BAPS) indicated different Central American origins for P. caribaea var. bahamensis and P. caribaea var. caribaea plastids, with Central America populations in northern and south-eastern groups. Sudden expansion times for BAPS clusters were close to three glacial maxima.

Conclusions

Central America contains ancestral plastid haplotypes. Population expansion has played a major role in the distribution of genetic diversity in P. caribaea var. hondurensis. Two colonization events gave rise to the P. caribaea var. bahamensis and P. caribaea var. caribaea lineages. Plastid variation in the eastern species (P. cubensis, P. maestrensis and P. occidentalis) evolved independently from that in P. caribaea var. caribaea. Incomplete lineage sorting between P. cubensis and P. maestrensis is apparent. Inferred expansion times for P. caribaea var. bahamensis and for the eastern lineages correspond to glacial maxima, whereas those for P. caribaea var. hondurensis correspond to the beginning of the temperature decrease that led to Marine Isotope Stage 8.