Molecular species delimitation of a symbiotic fig-pollinating wasp species complex reveals extreme deviation from reciprocal partner specificity

Background: Symbiotic relationships have contributed to major evolutionary innovations, the maintenance of fundamental ecosystem functions, and the generation and maintenance of biodiversity. However, the exact nature of host/symbiont associations, which has important consequences for their dynamics, is often poorly known due to limited understanding of symbiont taxonomy and species diversity. Among classical symbioses, figs and their pollinating wasps constitute a highly diverse keystone resource in tropical forest and savannah environments. Historically, they were considered to exemplify extreme reciprocal partner specificity (one-to-one host-symbiont species relationships), but recent work has revealed several more complex cases. However, there is a striking lack of studies with the specific aims of assessing symbiont diversity and how this varies across the geographic range of the host. Results: Here, we use molecular methods to investigate cryptic diversity in the pollinating wasps of a widespread Australian fig species. Standard barcoding genes and methods were not conclusive, but incorporation of phylogenetic analyses and a recently developed nuclear barcoding gene (ITS2), gave strong support for five pollinator species. Each pollinator species was most common in a different geographic region, emphasising the importance of wide geographic sampling to uncover diversity, and the scope for divergence in coevolutionary trajectories across the host plant range. In addition, most regions had multiple coexisting pollinators, raising the question of how they coexist in apparently similar or identical resource niches. Conclusion: Our study offers a striking example of extreme deviation from reciprocal partner specificity over the full geographical range of a fig-wasp system. It also suggests that superficially identical species may be able to co-exist in a mutualistic setting albeit at different frequencies in relation to their fig host’s range. We show that comprehensive sampling and molecular taxonomic techniques may be required to uncover the true structure of cryptic biodiversity underpinning intimate ecological interactions.

Characterisation of microsatellite markers for fig-pollinating wasps in the Pleistodontes imperialis species complex

We characterised a set of nine polymorphic microsatellite loci for Pleistodontes imperialis sp. 1, the pollinator wasp of Port Jackson fig (Ficus rubiginosa) in south-eastern Australia. Characterisation was performed on 30 female individuals collected from a population in Sydney, Australia. The average number of alleles per locus was 7.33, and eight loci were not in Hardy–Weinberg equilibrium. This was expected as fig wasps are known to be highly inbred. A test of genetic differentiation between two natural populations of P. imperialis sp. 1 (Sydney and Newcastle, Australia – some 120 km apart) yielded a very low FST value of 0.012, suggesting considerable gene flow. Bayesian clustering analysis using TESS 2.3.1, which does not assume Hardy–Weinberg equilibrium, however, indicated potential spatial substructuring between the Sydney and Newcastle populations, as well as within the Sydney population. The described loci were also characterised for two other species in the P. imperialis complex: P. imperialis sp. 2 (Townsville, Australia) and P. imperialis sp. 4 (Brisbane, Australia). Seven and six of the nine loci were polymorphic for P. imperialis sp. 2 and P.imperialis sp. 4, respectively.

Amplified fragment length polymorphism genetic fingerprinting challenges the taxonomic status of the near-endemic species Chara curta Nolte ex Kutz. (Characeae)

Amplified fragment length polymorphism (AFLP) genetic fingerprinting of 14 accessions of Chara curta and Chara aspera Willd., sampled across a range of habitats and morphologies in Britain, suggests that these taxa are part of the variation within a single species complex. Two primer combinations generating 397 fragments (97% of which were polymorphic), analysed by Jaccard's similarity coefficient and principal co-ordinate analysis, did not recover groups which reflect the current taxonomy. By contrast with the genetic study, a Gower general similarity coefficient and principal co-ordinate analysis of 52 morphological characters recovered the currently recognized species groups. A Mantel test showed no significant correlation between the genetic data and the morphological data, supporting the hypothesis that phenotypic variability in Chara L. is either to some extent environmentally induced or represents developmental stages. Implications for the conservation status of C. curta in Britain are discussed. (c) 2007 The Linnean Society of London, Botanical Journal of the Linnean Society, 2007, 155, 467-476.

Molecular genetic differentiation in earthworms inhabiting a heterogeneous Pb-polluted landscape

A Pb-mine site situated on acidic soil, but comprising of Ca-enriched islands around derelict buildings was used to study the spatial pattern of genetic diversity in Lumbricus rubellus. Two distinct genetic lineages ('A' and 'B'), differentiated at both the mitochondrial (mtDNA COII) and nuclear level (AFLPs) were revealed with a mean inter-lineage mtDNA sequence divergence of approximately 13%, indicative of a cryptic species complex. AFLP analysis indicates that lineage A individuals within one central 'ecological island' site are uniquely clustered, with little genetic overlap with lineage A individuals at the two peripheral sites. FTIR microspectroscopy of Pb-sequestering chloragocytes revealed different phosphate profiles in residents of adjacent acidic and calcareous islands. Bioinformatics found over-representation of Ca pathway genes in ESTPb libraries. Subsequent sequencing of a Ca-transport gene, SERCA, revealed mutations in the protein's cytosolic domain. We recommend the mandatory genotyping of all individuals prior to field-based ecotoxicological assays, particularly those using discriminating genomic technologies.

Genomic and genetic analyses of diversity and plant interactions of Pseudomonas fluorescens

Background: Pseudomonas fluorescens are common soil bacteria that can improve plant health through nutrient cycling, pathogen antagonism and induction of plant defenses. The genome sequences of strains SBW25 and Pf0-1 were determined and compared to each other and with P. fluorescens Pf-5. A functional genomic in vivo expression technology (IVET) screen provided insight into genes used by P. fluorescens in its natural environment and an improved understanding of the ecological significance of diversity within this species. Results: Comparisons of three P. fluorescens genomes (SBW25, Pf0-1, Pf-5) revealed considerable divergence: 61% of genes are shared, the majority located near the replication origin. Phylogenetic and average amino acid identity analyses showed a low overall relationship. A functional screen of SBW25 defined 125 plant-induced genes including a range of functions specific to the plant environment. Orthologues of 83 of these exist in Pf0-1 and Pf-5, with 73 shared by both strains. The P. fluorescens genomes carry numerous complex repetitive DNA sequences, some resembling Miniature Inverted-repeat Transposable Elements (MITEs). In SBW25, repeat density and distribution revealed 'repeat deserts' lacking repeats, covering approximately 40% of the genome. Conclusions: P. fluorescens genomes are highly diverse. Strain-specific regions around the replication terminus suggest genome compartmentalization. The genomic heterogeneity among the three strains is reminiscent of a species complex rather than a single species. That 42% of plant-inducible genes were not shared by all strains reinforces this conclusion and shows that ecological success requires specialized and core functions. The diversity also indicates the significant size of genetic information within the Pseudomonas pan genome.

The impact of genome defense on mobile elements in Microbotryum

Repeat induced point mutation (RIP), a mechanism causing hypermutation of repetitive DNA sequences in fungi, has been described as a ‘genome defense’ which functions to inactivate mobile elements and inhibit their deleterious effects on genome stability. Here we address the interactions between RIP and transposable elements in the Microbotryum violaceum species complex. Ten strains of M. violaceum, most of which belong to different species of the fungus, were all found to contain intragenomic populations of copia-like retrotransposons. Intragenomic DNA sequence variation among the copia-like elements was analyzed for evidence of RIP. Among species with RIP, there was no significant correlation between the frequency of RIP-induced mutations and inferred transposition rate based on diversity. Two strains of M. violaceum, from two different plant species but belonging to the same fungal lineage, contained copia-like elements with very low diversity, as would result from a high transposition rate, and these were also unique in showing no evidence of the hypermutation patterns indicative of the RIP genome defense. In this species, evidence of RIP was also absent from a Class II helitron-like transposable element. However, unexpectedly the absolute repetitive element load was lower than in other strains.

Diversity in a honey bee pathogen: first report of a third master variant of the Deformed Wing Virus quasispecies

Treatment of emerging RNA viruses is hampered by the high mutation and replication rates that enable these viruses to operate as a quasispecies. Declining honey bee populations have been attributed to the ectoparasitic mite Varroa destructor and its affiliation with Deformed Wing Virus (DWV). In the current study we use next-generation sequencing to investigate the DWV quasispecies in an apiary known to suffer from overwintering colony losses. We show that the DWV species complex is made up of three master variants. Our results indicate that a new DWV Type C variant is distinct from the previously described types A and B, but together they form a distinct clade compared with other members of the Iflaviridae. The molecular clock estimation predicts that Type C diverged from the other variants ~319 years ago. The discovery of a new master variant of DWV has important implications for the positive identification of the true pathogen within global honey bee populations.

Plant-pollinator biodiversity and pollination services in a complex Mediterranean landscape

Mediterranean landscapes comprise a complex mosaic of different habitats that vary in the diversity of their floral communities, pollinator communities and pollination services. Using the Greek Island of Lesvos as a model system, we assess the biodiversity value of six common habitats and measure ecosystemic 'health' using pollen grain deposition in three core flowering plants as a measure of pollination services. Three fire-driven habitats were assessed: freshly burnt areas, fully regenerated pine forests and intermediate age scrub; in addition we examined oak woodlands, actively managed olive groves and groves that had been abandoned from agriculture. Oak woodlands, pine forests and managed olive groves had the highest diversity of bees. The habitat characteristics responsible for structuring bee communities were: floral diversity, floral abundance, nectar energy availability and the variety of nectar resources present. Pollination services in two of our plant species, which were pollinated by a limited sub-set of the pollinator community, indicated that pollination levels were highest in the burnt and mature pine habitats. The third species, which was open to all flower visitors, indicated that oak woodlands had the highest levels of pollination from generalist species. Pollination was always more effective in managed olive groves than in abandoned groves. However, the two most common species of bee, the honeybee and a bumblebee, were not the primary pollinators within these habitats. We conclude that the three habitats of greatest overall value for plant-pollinator communities and provision of the healthiest pollination services are pine forests, oak woodland and managed olive groves. We indicate how the highest value habitats may be maintained in a complex landscape to safeguard and enhance pollination function within these habitats and potentially in adjoining agricultural areas. (c) 2005 Elsevier Ltd. All rights reserved.

Genetic analysis of complex demographic scenarios: Spatially expanding populations of the cane toad, Bufo marinus

Inferring the spatial expansion dynamics of invading species from molecular data is notoriously difficult due to the complexity of the processes involved. For these demographic scenarios, genetic data obtained from highly variable markers may be profitably combined with specific sampling schemes and information from other sources using a Bayesian approach. The geographic range of the introduced toad Bufo marinus is still expanding in eastern and northern Australia, in each case from isolates established around 1960. A large amount of demographic and historical information is available on both expansion areas. In each area, samples were collected along a transect representing populations of different ages and genotyped at 10 microsatellite loci. Five demographic models of expansion, differing in the dispersal pattern for migrants and founders and in the number of founders, were considered. Because the demographic history is complex, we used an approximate Bayesian method, based on a rejection-regression algorithm. to formally test the relative likelihoods of the five models of expansion and to infer demographic parameters. A stepwise migration-foundation model with founder events was statistically better supported than other four models in both expansion areas. Posterior distributions supported different dynamics of expansion in the studied areas. Populations in the eastern expansion area have a lower stable effective population size and have been founded by a smaller number of individuals than those in the northern expansion area. Once demographically stabilized, populations exchange a substantial number of effective migrants per generation in both expansion areas, and such exchanges are larger in northern than in eastern Australia. The effective number of migrants appears to be considerably lower than that of founders in both expansion areas. We found our inferences to be relatively robust to various assumptions on marker. demographic, and historical features. The method presented here is the only robust, model-based method available so far, which allows inferring complex population dynamics over a short time scale. It also provides the basis for investigating the interplay between population dynamics, drift, and selection in invasive species.

Temporal and spatial variability of mobile fauna on a submarine cliff and boulder scree complex: a community in flux

The mobile component of a community inhabiting a submarine boulder scree/cliff was investigated at Lough Hyne, Ireland at dawn, midday, dusk and night over a 1-week period. Line transects (50 m) were placed in the infralittoral (6 m) and circumlittoral (18 m) zones and also the interface between these two zones (12 m). The dominant mobile fauna of this cliff consisted of echinoderms (6 species), crustaceans (10 species) and fish (23 species). A different component community was identified at each time/depth interval using Multi-Dimensional Scaling (MDS) even though both species diversity (Shannon-Wiener indices) and richness (number of species) remained constant. These changes in community composition provided indirect evidence for migration by these mobile organisms. However, little evidence was found for migration between different zones with the exception of the several wrasse species. These species were observed to spend the daytime foraging in the deeper zone, but returned to the upper zone at night presumably for protection from predators. For the majority of species, migration was considered to occur to cryptic habitats such as holes and crevices. The number of organisms declined during the night, although crustacean numbers peaked, while fish and echinoderms were most abundant during day, possibly due to predator-prey interactions. This submarine community is in a state of flux, whereby, community characteristics, including trophic and energetic relationships, varied over small temporal (daily) and spatial (m) scales.

Time-resolved gas-phase kinetic and quantum chemical studies of reactions of silylene with chlorine-containing species. 1. HCl

Time-resolved kinetic studies of the reaction of silylene, SiH2, generated by laser flash photolysis of phenylsilane, have been carried out to obtain rate constants for its bimolecular reaction with HCL The reaction was studied in the gas phase at 10 Torr total pressure in SF6 bath gas, at five temperatures in the range of 296-611 K. The second-order rate constants fitted the Arrhenius equation: log(k/cm(3) molecule(-1) s(-1)) = (-11.51 +/- 0.06) + (1.92 +/- 0.47 kJ mol(-1))/RTIn10 Experiments at other pressures showed that these rate constants were unaffected by pressure in the range of 10-100 Torr, but showed small decreases in value of no more than 20% ( +/- 10%) at I Toff, at both the highest and lowest temperatures. The data are consistent with formation of an initial weakly bound donor-acceptor complex, which reacts by two parallel pathways. The first is by chlorine-to-silicon H-shift to make vibrationally excited chlorosilane, SiH3Cl*, which yields HSiCl by H-2 elimination from silicon. In the second pathway, the complex proceeds via H-2 elimination (4-center process) to make chlorosilylene, HSiCl, directly. This interpretation is supported by ab initio quantum calculations carried out at the G3 level which reveal the direct H-2 elimination route for the first time. RRKM modeling predicts the approximate magnitude of the pressure effect but is unable to determine the proportions of each pathway. The experimental data agree with the only previous measurements at room temperature. Comparisons with other reactions of SiH2 are also drawn.

Time-resolved gas-phase kinetic and quantum chemical studies of reactions of silylene with chlorine-containing species. 2. CH3Cl

Time-resolved kinetic studies of the reaction of silylene, SiH2, generated by laser flash photolysis of both silacyclopent-3-ene and phenylsilane, have been carried out to obtain second-order rate constants for its reaction with CH3Cl. The reaction was studied in the gas phase at six temperatures in the range 294-606 K. The second-order rate constants gave a curved Arrhenius plot with a minimum value at T approximate to 370 K. The reaction showed no pressure dependence in the presence of up to 100 Torr SF6. The rate constants, however, showed a weak dependence on laser pulse energy. This suggests an interpretation requiring more than one contributing reaction pathway to SiH2 removal. Apart from a direct reaction of SiH2 with CH3Cl, reaction of SiH2 with CH3 (formed by photodissociation of CH3Cl) seems probable, with contributions of up to 30% to the rates. Ab initio calculations (G3 level) show that the initial step of reaction of SiH2 with CH3Cl is formation of a zwitterionic complex (ylid), but a high-energy barrier rules out the subsequent insertion step. On the other hand, the Cl-abstraction reaction leading to CH3 + ClSiH2 has a low barrier, and therefore, this seems the most likely candidate for the main reaction pathway of SiH2 with CH3Cl. RRKM calculations on the abstraction pathway show that this process alone cannot account for the observed temperature dependence of the rate constants. The data are discussed in light of studies of other silylene reactions with haloalkanes.

Dicopper(II) complexes of a new di-para-xylyldioxatetraazamacrocycle and cascade species with dicarboxylate anions: thermodynamics and structural properties

The new dioxatetraazamacrocycle (L-1) was synthesized by a 2 + 2 condensation and characterized. Stability constants of its copper(II) complexes were determined by spectrophotometry in DMSO at 298.2 K in 0. 10 mol dm(-3) KClO4. Mainly dinuclear complexes are formed and the presence of mononuclear species is dependent on the counterion (Cl- or ClO4-). The association constants of the dinuclear copper(II) complexes with dicarboxylate anions [oxalate (oxa(2-)), malonate (mal(2-)), succinate (suc(2-)), and glutarate (glu(2-))] were also determined by spectrophotometry at 298.2 K in DMSO, and it was found that values decrease with an increase of the alkyl chain between the carboxylate groups. X-Band EPR spectra of the dicopper(II) complexes and of their cascade species in frozen DMSO exhibit dipole-dipole coupling, and their simulation, together with their UV-vis spectra, showed that the copper centres of the complexes in solution had square pyramidal geometries though with different distortions. From the experimental data, it was also possible to predict the Cu...Cu distances, the minimum being found at 6.4 angstrom for the (Cu2LCl4)-Cl-1 complex and then successively this distance slightly increases when the chloride anions are replaced by dicarboxylate anions, from 6.6 angstrom for oxa(2-) to 7.8 for glu(2-). The crystal structures of the dinuclear copper cascade species with oxa(2-) and suc(2-) were determined and showed one anion bridging both copper centres and Cu...Cu distances of 5.485(7) angstrom and 6.442(8) angstrom, respectively.

A square-planar nickel(II) monoradical complex with a bis(salicylidene)diamine ligand

The new square-planar Ni-II-N2O2 complex [Ni(L-Me)] (1(Me)), where L-Me, stands for the dianionic phenolato form of N,N'bis(3,5-di-tert-butyl-salicylidene)-4,5-dimethyl-1,2-phenyl- enediamine ((LH2)-L-Me), has been synthesised and fully characterised. X-ray crystallography was also used for the characterisation. The electrochemical one-electron oxidation of 1(Me) produces the thermally stable (within the temperature range 10-295 K) cationic species (1(Me))(+). The UV/Vis and X-band EPR experimental data, supported by DFT calculations, indicate that (1(Me))(+), is best described as a Ni-II monoradical complex and, thus, does NOT exist in a Ni-III ground state, in contrast to its demethylated counterpart [Ni(L-H)](+) (1(H))(+) below 170 K.