Phylogenetic signal in AFLP data sets

One, two or three? Integrative species delimitation of short-range endemic Hemicycla species (Gastropoda: Helicidae) from the Canary Islands based on morphology, barcoding, AFLP and ddRADseq data

Hemicycla mascaensis and H. diegoi are short-range endemics that occur allopatrically in small areas in the Teno Mountains in the western part of Tenerife (Canary Islands). Both taxa have been recognised as distinct species based on differences in shell morphology and genital anatomy. Preliminary molecular analyses using mitochondrial markers suggested a potential paraphyly of H. diegoi with regard to H. mascaensis. We here use multilocus AFLP data and ddRADseq data as well as distribution data, data on shell morphology and genital anatomy to assess the status of these taxa using phylogenetic analyses, species tree reconstruction and molecular species delimitation based on the multispecies coalescent as implemented in BFD* and BPP in an integrative approach. Our analyses show that, based on the analysis of multilocus data, the two taxa are reciprocally monophyletic. Species delimitation methods, however, tend to recognise all investigated populations as distinct species, albeit neither lending unambiguous support to any of the species hypotheses. The comparison of the anatomy of distal genital organs further suggests differentiation within H. mascaensis. This highlights the need for a balanced weighting of arguments from different lines of evidence to determine species status and calls for cautious interpretations of the results of molecular species delimitation analyses, especially in organisms with low active dispersal capacities and expected distinct population structuring such as land snails. Taking all available evidence into account, we favour to recognise H. mascaensis and H. diegoi as distinct species, acknowledging, though, that the recognition of both taxa as subspecies (with possibly a third yet undescribed) would also be an option as morphological differentiation is within the limits of other land snail species that are traditionally subdivided into subspecies.

Phylogeography of western Mediterranean Cymbalaria (Plantaginaceae) reveals two independent long-distance dispersals and entails new taxonomic circumscriptions

The Balearic Islands, Corsica and Sardinia (BCS) constitute biodiversity hotspots in the western Mediterranean Basin. Oligocene connections and long distance dispersal events have been suggested to cause presence of BCS shared endemic species. One of them is Cymbalaria aequitriloba, which, together with three additional species, constitute a polyploid clade endemic to BCS. Combining amplified fragment length polymorphism (AFLP) fingerprinting, plastid DNA sequences and morphometrics, we inferred the phylogeography of the group and evaluated the species’ current taxonomic circumscriptions. Based on morphometric and AFLP data we propose a new circumscription for C. fragilis to additionally comprise a group of populations with intermediate morphological characters previously included in C. aequitriloba. Consequently, we suggest to change the IUCN category of C. fragilis from critically endangered (CR) to near threatened (NT). Both morphology and AFLP data support the current taxonomy of the single island endemics C. hepaticifolia and C. muelleri. The four species had a common origin in Corsica-Sardinia, and two long-distance dispersal events to the Balearic Islands were inferred. Finally, plastid DNA data suggest that interspecific gene flow took place where two species co-occur.

Use of AFLP for the study of eukaryotic pathogens affecting humans

Amplified fragment length polymorphism (AFLP) is a genotyping technique based on PCR amplification of specific restriction fragments from a particular genome. The methodology has been extensively used in plant biology to solve a variety of scientific questions, including taxonomy, molecular epidemiology, systematics, population genetics, among many others. The AFLP share advantages and disadvantages with other types of molecular markers, being particularly useful in organisms with no previous DNA sequence knowledge. In eukaryotic pathogens, the technique has not been extensively used, although it has the potential to solve many important issues as it allows the simultaneous examination of hundreds or even thousands of polymorphic sites in the genome of the organism. Here we describe the main applications published on the use of AFLP in eukaryotic pathogens, with emphasis in species of the groups fungi, protozoa and helminths, and discuss the role of this methodology in the context of new techniques derived from the advances of the next generation sequencing.

Genotypic diversity in root-endophytic fungi reflects efficient dispersal and environmental adaptation

Studying community structure and dynamics of plant-associated fungi is the basis for unravelling their interactions with hosts and ecosystem functions. A recent sampling revealed that only a few fungal groups, as defined by internal transcribed spacer region (ITS) sequence similarity, dominate culturable root endophytic communities of nonmycorrhizal Microthlaspi spp. plants across Europe. Strains of these fungi display a broad phenotypic and functional diversity, which suggests a genetic variability masked by ITS clustering into operational taxonomic units (OTUs). The aims of this study were to identify how genetic similarity patterns of these fungi change across environments and to evaluate their ability to disperse and adapt to ecological conditions. A first ITS-based haplotype analysis of ten widespread OTUs mostly showed a low to moderate genotypic differentiation, with the exception of a group identified as Cadophora sp. that was highly diverse. A multilocus phylogeny based on additional genetic loci (partial translation elongation factor 1α, beta-tubulin and actin) and amplified fragment length polymorphism profiling of 185 strains representative of the five dominant OTUs revealed a weak association of genetic differences with geography and environmental conditions, including bioclimatic and soil factors. Our findings suggest that dominant culturable root endophytic fungi have efficient dispersal capabilities, and that their distribution is little affected by environmental filtering. Other processes, such as inter- and intraspecific biotic interactions, may be more important for the local assembly of their communities.

Snails in the desert: Species diversification of Theba (Gastropoda: Helicidae) along the Atlantic coast of NW Africa

The spatial subdivision of species often plays a pivotal role in speciation. Across their entire range, species are rarely panmictic and crucial consequences of spatial subdivision are (1) random genetic drift including historical factors, (2) uniform selection, and (3) divergent selection. Each of these consequences may result in geographic variation and eventually reproductive isolation, but their relative importance in speciation is still unclear. In this study, we used a combination of genetic, morphological, and climatic data to obtain a comprehensive picture of differentiation among three closely related, parapatrically distributed taxa of the land snail genus Theba occurring along the Atlantic coasts of South Morocco and Western Sahara. We conducted Mantel and partial Mantel tests to relate phenotypic and genotypic variation of these species to geography and/or climate. As null hypothesis for an evolutionary scenario, we assumed nonadaptive speciation and expected a pattern of isolation by distance among taxa. Rejection of the null hypothesis would indicate isolation by environment due to adaptation. Generally, genetic drift plays an important role but is rarely considered as sole driver of speciation. It is the combination of drift and selection that predominantly drives speciation. This study, however, provides a potential example, in which nonadaptive speciation, that is, genetic drift, is apparently the main driver of shaping the diversity of Theba in NW Africa. Restriction of gene flow between populations caused by geographic isolation probably has played an important role. Climate oscillations during the Plio- and Pleistocene may have led to repeated ecological changes in NW Africa and disruptions of habitats promoting differentiation by geographic isolation. The inferred evolutionary scenario, however, did not fully explain the incongruence between the AFLP- and mtDNA-tree topologies. This incongruence might indicate past hybridization among the studied Theba forms.

Geography shapes the phylogeny of frailejones (Espeletiinae Cuatrec., Asteraceae): a remarkable example of recent rapid radiation in sky islands

Background

The páramo ecosystem, located above the timberline in the tropical Andes, has been the setting for some of the most dramatic plant radiations, and it is one of the world’s fastest evolving and most diverse high-altitude ecosystems. Today 144+ species of frailejones (subtribe Espeletiinae Cuatrec., Asteraceae) dominate the páramo. Frailejones have intrigued naturalists and botanists, not just for their appealing beauty and impressive morphological diversity, but also for their remarkable adaptations to the extremely harsh environmental conditions of the páramo. Previous attempts to reconstruct the evolutionary history of this group failed to resolve relationships among genera and species, and there is no agreement regarding the classification of the group. Thus, our goal was to reconstruct the phylogeny of the frailejones and to test the influence of the geography on it as a first step to understanding the patterns of radiation of these plants.

Methods

Field expeditions in 70 páramos of Colombia and Venezuela resulted in 555 collected samples from 110 species. Additional material was obtained from herbarium specimens. Sequence data included nrDNA (ITS and ETS) and cpDNA (rpl16), for an aligned total of 2,954 bp. Fragment analysis was performed with AFLP data using 28 primer combinations and yielding 1,665 fragments. Phylogenies based on sequence data were reconstructed under maximum parsimony, maximum likelihood and Bayesian inference. The AFLP dataset employed minimum evolution analyses. A Monte Carlo permutation test was used to infer the influence of the geography on the phylogeny.

Results

Phylogenies reconstructed suggest that most genera are paraphyletic, but the phylogenetic signal may be misled by hybridization and incomplete lineage sorting. A tree with all the available molecular data shows two large clades: one of primarily Venezuelan species that includes a few neighboring Colombian species; and a second clade of only Colombian species. Results from the Monte Carlo permutation test suggests a very strong influence of the geography on the phylogenetic relationships. Venezuelan páramos tend to hold taxa that are more distantly-related to each other than Colombian páramos, where taxa are more closely-related to each other.

Conclusions

Our data suggest the presence of two independent radiations: one in Venezuela and the other in Colombia. In addition, the current generic classification will need to be deeply revised. Analyses show a strong geographic structure in the phylogeny, with large clades grouped in hotspots of diversity at a regional scale, and in páramo localities at a local scale. Differences in the degrees of relatedness between sympatric species of Venezuelan and Colombian páramos may be explained because of the younger age of the latter páramos, and the lesser time for speciation of Espeletiinae in them.

A novel method to infer the origin of polyploids from Amplified Fragment Length Polymorphism data reveals that the alpine polyploid complex of Senecio carniolicus (Asteraceae) evolved mainly via autopolyploidy

Despite its evolutionary and ecological relevance, the mode of polyploid origin has been notoriously difficult to be reconstructed from molecular data. Here, we present a method to identify the putative parents of polyploids and thus to infer the mode of their origin (auto- vs. allopolyploidy) from Amplified Fragment Length Polymorphism (AFLP) data. To this end, we use Cohen's d of distances between in silico polyploids, generated within a priori defined scenarios of origin from a priori delimited putative parental entities (e.g. taxa, genetic lineages), and natural polyploids. Simulations show that the discriminatory power of the proposed method increases mainly with increasing divergence between the lower-ploid putative ancestors and less so with increasing delay of polyploidization relative to the time of divergence. We apply the new method to the Senecio carniolicus aggregate, distributed in the European Alps and comprising two diploid, one tetraploid and one hexaploid species. In the eastern part of its distribution, the S. carniolicus aggregate was inferred to comprise an autopolyploid series, whereas for western populations of the tetraploid species, an allopolyploid origin involving the two diploid species was the most likely scenario. Although this suggests that the tetraploid species has two independent origins, other evidence (ribotype distribution, morphology) is consistent with the hypothesis of an autopolyploid origin with subsequent introgression by the second diploid species. Altogether, identifying the best among alternative scenarios using Cohen's d can be straightforward, but particular scenarios, such as allopolyploid origin vs. autopolyploid origin with subsequent introgression, remain difficult to be distinguished.

Diversification of Bromelioideae (Bromeliaceae) in the Brazilian Atlantic rainforest: A case study in Aechmea subgenus Ortgiesia

Aechmea subgenus Ortgiesia comprises ca. 20 species distributed in Brazil, Argentina, Paraguay, and Uruguay, with a center of diversity in the Brazilian Atlantic rainforest. We examined interspecific relationships of Ortgiesia based on Amplified Fragment Length Polymorphisms (AFLP). Ninety-six accessions belonging to 14 species of Ortgiesia were sampled, and genotyped with 11 AFLP primer combinations. The neighbor joining (NJ) tree depicted two main genetic groups within Aechmea subgenus Ortgiesia, and four subgroups. The NJ tree showed short internal branches, indicating an overall shallow genetic divergence among Ortgiesia species as expected for the recently radiated subfamily Bromelioideae. Our results suggest that hybridization and/or incomplete lineage sorting may have hampered the reconstruction of interspecific relationships in Aechmea subgenus Ortgiesia. The mapping of petal color (yellow, blue, pink, or white), inflorescence type (simple or compound), and inflorescence shape (ellipsoid, subcylindric, cylindric, or pyramidal) against the NJ tree indicated that these characters are of limited taxonomic use in Aechmea subgenus Ortgiesia due to homoplasy. An analysis of the current distribution of Ortgiesia identified the southern region of the Brazilian Atlantic rainforest, between latitudes of 26° and 27°S, as the center of diversity for the subgenus.

Chemical and genetic characterization of Phlomis species and wild hybrids in Crete

The genus Phlomis is represented in the island of Crete (Greece, Eastern Mediterranean) by three species Phlomis cretica C. Presl., Phlomis fruticosa L., the island endemic Phlomis lanata Willd. and three hybrids Phlomis x cytherea Rech.f. (P. cretica x P. fruticosa), Phlomis x commixta Rech.f. (P. cretica x P. lanata) and Phlomis x sieberi Vierh. (P. fruticosa x P. lanata). This work describes (a) the profile of hybrids and parental species concerning their volatile compounds, (b) the suitability of ribosomal nuclear (ITS region), chloroplast (trnH-psbA), and AFLP markers to identify hybrids and (c) their competence to characterize the different chemotypes of both hybrids and their parental species. The cluster analysis and PCA constructed from chemical data (volatile oils) suggest that there are three groups of taxa. Group IA includes P. cretica and P. fruticosa, group IB includes P. x cytherea, whereas group II consists of P. x commixta, P. x sieberi and P. lanata. Volatile compounds detected only in the hybrids P. x sieberi and P. x commixta correspond to the 3% of the total compounds, value that is much higher in P. x cytherea (21%). Neighbor-joining, statistical parsimony analysis and the observations drawn from ribotypes spectrum of ITS markers divided Phlomis species in two groups, P. lanata and the complex P. cretica/P. fruticosa. In contrast to the ITS region, the plastid DNA marker follows a geographically related pattern. Neighbor-Net, PCA and Bayesian assignment analysis performed for AFLP markers separated the genotypes into three groups corresponding to populations of P. cretica, P. fruticosa, and P. lanata, respectively, while populations of P. x commixta, P. x cytherea, and P. x sieberi presented admixed ancestry. Most of the P. x cytherea samples were identified as F1 hybrids by Bayesian assignment test, while those of P. x commixta and P. x sieberi were identified as F2 hybrids. Overall, high chemical differentiation is revealed in one of the three hybrids, which is likely related with niche variation. Moreover, molecular markers show potential to identify Phlomis taxa.

Key processes for Cheirolophus (Asteraceae) diversification on oceanic islands inferred from AFLP data

The radiation of the genus Cheirolophus (Asteraceae) in Macaronesia constitutes a spectacular case of rapid diversification on oceanic islands. Twenty species - nine of them included in the IUCN Red List of Threatened Species - have been described to date inhabiting the Madeiran and Canarian archipelagos. A previous phylogenetic study revealed that the diversification of Cheirolophus in Macaronesia started less than 2 Ma. As a result of such an explosive speciation process, limited phylogenetic resolution was reported, mainly due to the low variability of the employed molecular markers. In the present study, we used highly polymorphic AFLP markers to i) evaluate species' boundaries, ii) infer their evolutionary relationships and iii) investigate the patterns of genetic diversity in relation to the potential processes likely involved in the radiation of Cheirolophus. One hundred and seventy-two individuals representing all Macaronesian Cheirolophus species were analysed using 249 AFLP loci. Our results suggest that geographic isolation played an important role in this radiation process. This was likely driven by the combination of poor gene flow capacity and a good ability for sporadic long-distance colonisations. In addition, we also found some traces of introgression and incipient ecological adaptation, which could have further enhanced the extraordinary diversification of Cheirolophus in Macaronesia. Last, we hypothesize that current threat categories assigned to Macaronesian Cheirolophus species do not reflect their respective evolutionary relevance, so future evaluations of their conservation status should take into account the results presented here.

Dealing with discordant genetic signal caused by hybridisation, incomplete lineage sorting and paucity of primary nucleotide homologies: a case study of closely related members of the genus Picris subsection Hieracioides (Compositae)

We investigated genetic variation and evolutionary history of closely related taxa of Picris subsect. Hieracioides with major focus on the widely distributed P. hieracioides and its closely related congeners, P. hispidissima, P. japonica, P. olympica, and P. nuristanica. Accessions from 140 sample sites of the investigated Picris taxa were analyzed on the infra- and the inter-specific level using nuclear (ITS1-5.8S-ITS2 region) and chloroplast (rpl32-trnL^(UAG) region) DNA sequences. Genetic patterns of P. hieracioides, P. hispidissima, and P. olympica were shown to be incongruent and, in several cases, both plastid and nuclear alleles transcended borders of the taxa and genetic lineages. The widespread P. hieracioides was genetically highly variable and non-monophyletic across both markers, with allele groups having particular geographic distributions. Generally, all gene trees and networks displayed only a limited and statistically rather unsupported resolution among ingroup taxa causing their phylogenetic relationships to remain rather unresolved. More light on these intricate evolutionary relationships was cast by the Bayesian coalescent-based analysis, although some relationships were still left unresolved. A combination of suite of phylogenetic analyses revealed the ingroup taxa to represent a complex of genetically closely related and morphologically similar entities that have undergone a highly dynamic and recent evolution. This has been especially affected by the extensive and recurrent gene flow among and within the studied taxa and/or by the maintenance of ancestral variation. Paucity of phylogenetically informative signal further hampers the reconstruction of relationships on the infra- as well as on the inter-specific level. In the present study, we have demonstrated that a combination of various phylogenetic analyses of datasets with extremely complex and incongruent phylogenetic signal may shed more light on the interrelationships and evolutionary history of analysed species groups.

Sequence-related amplified polymorphism (SRAP) markers: A potential resource for studies in plant molecular biology(1.)

In the past few decades, many investigations in the field of plant biology have employed selectively neutral, multilocus, dominant markers such as inter-simple sequence repeat (ISSR), random-amplified polymorphic DNA (RAPD), and amplified fragment length polymorphism (AFLP) to address hypotheses at lower taxonomic levels. More recently, sequence-related amplified polymorphism (SRAP) markers have been developed, which are used to amplify coding regions of DNA with primers targeting open reading frames. These markers have proven to be robust and highly variable, on par with AFLP, and are attained through a significantly less technically demanding process. SRAP markers have been used primarily for agronomic and horticultural purposes, developing quantitative trait loci in advanced hybrids and assessing genetic diversity of large germplasm collections. Here, we suggest that SRAP markers should be employed for research addressing hypotheses in plant systematics, biogeography, conservation, ecology, and beyond. We provide an overview of the SRAP literature to date, review descriptive statistics of SRAP markers in a subset of 171 publications, and present relevant case studies to demonstrate the applicability of SRAP markers to the diverse field of plant biology. Results of these selected works indicate that SRAP markers have the potential to enhance the current suite of molecular tools in a diversity of fields by providing an easy-to-use, highly variable marker with inherent biological significance.

A framework phylogeny of the American oak clade based on sequenced RAD data

Previous phylogenetic studies in oaks (Quercus, Fagaceae) have failed to resolve the backbone topology of the genus with strong support. Here, we utilize next-generation sequencing of restriction-site associated DNA (RAD-Seq) to resolve a framework phylogeny of a predominantly American clade of oaks whose crown age is estimated at 23-33 million years old. Using a recently developed analytical pipeline for RAD-Seq phylogenetics, we created a concatenated matrix of 1.40 E06 aligned nucleotides, constituting 27,727 sequence clusters. RAD-Seq data were readily combined across runs, with no difference in phylogenetic placement between technical replicates, which overlapped by only 43-64% in locus coverage. 17% (4,715) of the loci we analyzed could be mapped with high confidence to one or more expressed sequence tags in NCBI Genbank. A concatenated matrix of the loci that BLAST to at least one EST sequence provides approximately half as many variable or parsimony-informative characters as equal-sized datasets from the non-EST loci. The EST-associated matrix is more complete (fewer missing loci) and has slightly lower homoplasy than non-EST subsampled matrices of the same size, but there is no difference in phylogenetic support or relative attribution of base substitutions to internal versus terminal branches of the phylogeny. We introduce a partitioned RAD visualization method (implemented in the R package RADami; http://cran.r-project.org/web/packages/RADami) to investigate the possibility that suboptimal topologies supported by large numbers of loci--due, for example, to reticulate evolution or lineage sorting--are masked by the globally optimal tree. We find no evidence for strongly-supported alternative topologies in our study, suggesting that the phylogeny we recover is a robust estimate of large-scale phylogenetic patterns in the American oak clade. Our study is one of the first to demonstrate the utility of RAD-Seq data for inferring phylogeny in a 23-33 million year-old clade.

Impact of deep coalescence and recombination on the estimation of phylogenetic relationships among species using AFLP markers

Deep coalescence and the nongenealogical pattern of descent caused by recombination have emerged as a common problem for phylogenetic inference at the species level. Here we use computer simulations to assess whether AFLP-based phylogenies are robust to the uncertainties introduced by these factors. Our results indicate that phylogenetic signal can prevail even in the face of extensive deep coalescence allowing recovering the correct species tree topology. The impact of recombination on tree accuracy was related to total tree depth and species effective population size. The correct tree topology could be recovered upon many simulation settings due to a trade-off between the conflicting signals resulting from intra-locus recombination and the benefits of the joint consideration of unlinked loci that better matched overall the true species tree. Errors in tree topology were not only determined by deep coalescence, but also by the timing of divergence and the tree-building errors arising from an insufficient number of characters. DNA sequences generally outperformed AFLPs upon any simulated scenario, but this difference in performance was nearly negligible when a sufficient number of AFLP characters were sampled. Our simulations suggest that the impact of deep coalescence and intra-locus recombination on the reliability of AFLP trees could be minimal for effective population sizes equal to or lower than 10,000 (typical of many vertebrates and tree plants) given tree depths above 0.02 substitutions per site.

Outgroup effects on root position and tree topology in the AFLP phylogeny of a rapidly radiating lineage of cichlid fish

Phylogenetic analyses of rapid radiations are particularly challenging as short basal branches and incomplete lineage sorting complicate phylogenetic inference. Multilocus data of presence-absence polymorphisms such as obtained by AFLP genotyping overcome some of the difficulties, but also present their own intricacies. Here we analyze >1000 AFLP markers to address the evolutionary history of the Limnochromini, a cichlid fish lineage endemic to Lake Tanganyika, and to test for potential effects of outgroup composition on tree topology. The data support previous mitochondrial evidence on the tribe's taxonomy by confirming the polyphyly of the genus Limnochromis and - in contradiction to a recent taxonomic revision - nesting the genus Greenwoodochromis within the Limnochromini. Species relationships suggest that ecological segregation occurred during the rapid basal radiation of the Limnochromini. The large phylogenetic distance between candidate outgroup taxa and the Limnochromini radiation caused random outgroup effects. Bootstrap support for ingroup nodes was lower in outgroup-rooted than in midpoint-rooted trees, and root positions and ingroup tree topologies varied in response to the composition of the outgroup. These observations suggest that the predisposition for homoplastic evolution makes AFLP-based phylogenetic analyses particularly susceptible to random biases introduced by too-distant outgroup taxa.

Analyses of amplified fragment length polymorphisms (AFLP) indicate rapid radiation of Diospyros species (Ebenaceae) endemic to New Caledonia

Background

Radiation in some plant groups has occurred on islands and due to the characteristic rapid pace of phenotypic evolution, standard molecular markers often provide insufficient variation for phylogenetic reconstruction. To resolve relationships within a clade of 21 closely related New Caledonian Diospyros species and evaluate species boundaries we analysed genome-wide DNA variation via amplified fragment length polymorphisms (AFLP).

Results

A neighbour-joining (NJ) dendrogram based on Dice distances shows all species except D. minimifolia, D. parviflora and D. vieillardii to form unique clusters of genetically similar accessions. However, there was little variation between these species clusters, resulting in unresolved species relationships and a star-like general NJ topology. Correspondingly, analyses of molecular variance showed more variation within species than between them. A Bayesian analysis with BEAST produced a similar result. Another Bayesian method, this time a clustering method, Structure, demonstrated the presence of two groups, highly congruent with those observed in a principal coordinate analysis (PCO). Molecular divergence between the two groups is low and does not correspond to any hypothesised taxonomic, ecological or geographical patterns.

Conclusions

We hypothesise that such a pattern could have been produced by rapid and complex evolution involving a widespread progenitor for which an initial split into two groups was followed by subsequent fragmentation into many diverging populations, which was followed by range expansion of then divergent entities. Overall, this process resulted in an opportunistic pattern of phenotypic diversification. The time since divergence was probably insufficient for some species to become genetically well-differentiated, resulting in progenitor/derivative relationships being exhibited in a few cases. In other cases, our analyses may have revealed evidence for the existence of cryptic species, for which more study of morphology and ecology are now required.

Hybridization and adaptation to introduced balloon vines in an Australian soapberry bug

Contemporary adaptation of plant feeding insects to introduced hosts provides clear cases of ecologically based population divergence. In most cases the mechanisms permitting rapid differentiation are not well known. Here we study morphological and genetic variation associated with recent shifts by the Australian soapberry bug Leptocoris tagalicus onto two naturalized Neotropical balloon vines, Cardiospermum halicacabum and C. grandiflorum that differ in time since introduction. Our results show that these vines have much larger fruits than the native hosts (Whitewood tree -Atalaya hemiglauca- and Woolly Rambutan -Alectryon tomentosus-) and that bugs living on them have evolved significantly longer beaks and new allometries. Genetic analyses of mitochondrial haplotypes and amplified fragment length polymorphic (AFLP) markers indicate that the lineage of bugs on the annual vine C. halicacabum, the older introduction, is intermediate between the two subspecies of L. tagalicus found on native hosts. Moreover, where the annual vine and Whitewood tree co-occur, the morphology and genomic composition of the bugs are similar to those occurring in allopatry. These results show that hybridization provided the genetic elements underlying the strongly differentiated 'Halicacabum bugs'. In contrast, the bugs feeding on the recently introduced perennial balloon vine (C. grandiflorum) showed no evidence of admixture, and are genetically indistinguishable from the nearby populations on a native host.

Genomic distribution of AFLP markers relative to gene locations for different eukaryotic species

BACKGROUND

Amplified fragment length polymorphism (AFLP) markers are frequently used for a wide range of studies, such as genome-wide mapping, population genetic diversity estimation, hybridization and introgression studies, phylogenetic analyses, and detection of signatures of selection. An important issue to be addressed for some of these fields is the distribution of the markers across the genome, particularly in relation to gene sequences.

RESULTS

Using in-silico restriction fragment analysis of the genomes of nine eukaryotic species we characterise the distribution of AFLP fragments across the genome and, particularly, in relation to gene locations. First, we identify the physical position of markers across the chromosomes of all species. An observed accumulation of fragments around (peri) centromeric regions in some species is produced by repeated sequences, and this accumulation disappears when AFLP bands rather than fragments are considered. Second, we calculate the percentage of AFLP markers positioned within gene sequences. For the typical EcoRI/MseI enzyme pair, this ranges between 28 and 87% and is usually larger than that expected by chance because of the higher GC content of gene sequences relative to intergenic ones. In agreement with this, the use of enzyme pairs with GC-rich restriction sites substantially increases the above percentages. For example, using the enzyme system SacI/HpaII, 86% of AFLP markers are located within gene sequences in A. thaliana, and 100% of markers in Plasmodium falciparun. We further find that for a typical trait controlled by 50 genes of average size, if 1000 AFLPs are used in a study, the number of those within 1 kb distance from any of the genes would be only about 1-2, and only about 50% of the genes would have markers within that distance.

CONCLUSIONS

The high coverage of AFLP markers across the genomes and the high proportion of markers within or close to gene sequences make them suitable for genome scans and detecting large islands of differentiation in the genome. However, for specific traits, the percentage of AFLP markers close to genes can be rather small. Therefore, genome scans directed towards the search of markers closely linked to selected loci can be a difficult task in many instances.

Phylogeny, floral evolution, and inter-island dispersal in Hawaiian Clermontia (Campanulaceae) based on ISSR variation and plastid spacer sequences

Previous studies based on DNA restriction-site and sequence variation have shown that the Hawaiian lobeliads are monophyletic and that the two largest genera, Cyanea and Clermontia, diverged from each other ca. 9.7 Mya. Sequence divergence among species of Clermontia is quite limited, however, and extensive hybridization is suspected, which has interfered with production of a well-resolved molecular phylogeny for the genus. Clermontia is of considerable interest because several species posses petal-like sepals, raising the question of whether such a homeotic mutation has arisen once or several times. In addition, morphological and molecular studies have implied different patterns of inter-island dispersal within the genus. Here we use nuclear ISSRs (inter-simple sequence repeat polymorphisms) and five plastid non-coding sequences to derive biparental and maternal phylogenies for Clermontia. Our findings imply that (1) Clermontia is not monophyletic, with Cl. pyrularia nested within Cyanea and apparently an intergeneric hybrid; (2) the earliest divergent clades within Clermontia are native to Kauài, then Òahu, then Maui, supporting the progression rule of dispersal down the chain toward progressively younger islands, although that rule is violated in later-evolving taxa in the ISSR tree; (3) almost no sequence divergence among several Clermontia species in 4.5 kb of rapidly evolving plastid DNA; (4) several apparent cases of hybridization/introgression or incomplete lineage sorting (i.e., Cl. oblongifolia, peleana, persicifolia, pyrularia, samuelii, tuberculata), based on extensive conflict between the ISSR and plastid phylogenies; and (5) two origins and two losses of petaloid sepals, or--perhaps more plausibly--a single origin and two losses of this homeotic mutation, with its introgression into Cl. persicifolia. Our phylogenies are better resolved and geographically more informative than others based on ITS and 5S-NTS sequences and nuclear SNPs, but agree with them in supporting Clermontia's origin on Kauài or some older island and dispersal down the chain subsequently.

Automated masking of AFLP markers improves reliability of phylogenetic analyses

The amplified fragment length polymorphisms (AFLP) method has become an attractive tool in phylogenetics due to the ease with which large numbers of characters can be generated. In contrast to sequence-based phylogenetic approaches, AFLP data consist of anonymous multilocus markers. However, potential artificial amplifications or amplification failures of fragments contained in the AFLP data set will reduce AFLP reliability especially in phylogenetic inferences. In the present study, we introduce a new automated scoring approach, called “AMARE” (AFLP MAtrix REduction). The approach is based on replicates and makes marker selection dependent on marker reproducibility to control for scoring errors. To demonstrate the effectiveness of our approach we record error rate estimations, resolution scores, PCoA and stemminess calculations. As in general the true tree (i.e. the species phylogeny) is not known, we tested AMARE with empirical, already published AFLP data sets, and compared tree topologies of different AMARE generated character matrices to existing phylogenetic trees and/or other independent sources such as morphological and geographical data. It turns out that the selection of masked character matrices with highest resolution scores gave similar or even better phylogenetic results than the original AFLP data sets.

Radiating on oceanic islands: patterns and processes of speciation in the land snail genus Theba (Risso 1826)

Island radiations have played a major role in shaping our current understanding of allopatric, sympatric and parapatric speciation. However, the fact that species divergence correlates with island size emphasizes the importance of geographic isolation (allopatry) in speciation. Based on molecular and morphological data, we investigated the diversification of the land snail genus Theba on the two Canary Islands of Lanzarote and Fuerteventura. Due to the geological history of both islands, this study system provides ideal conditions to investigate the interplay of biogeography, dispersal ability and differentiation in generating species diversity. Our analyses demonstrated extensive cryptic diversification of Theba on these islands, probably driven mainly by non-adaptive allopatric differentiation and secondary gene flow. In a few cases, we observed a complete absence of gene flow among sympatrically distributed forms suggesting an advanced stage of speciation. On the Jandía peninsula genome scans suggested genotype-environment associations and potentially adaptive diversification of two closely related Theba species to different ecological environments. We found support for the idea that genetic differentiation was enhanced by divergent selection in different environments. The diversification of Theba on both islands is therefore best explained by a mixture of non-adaptive and adaptive speciation, promoted by ecological and geomorphological factors.

AFLP markers provide insights into the evolutionary relationships and diversification of New Caledonian Araucaria species (Araucariaceae)

PREMISE OF THE STUDY

Despite its small size, New Caledonia is characterized by a very diverse flora and striking environmental gradients, which make it an ideal setting to study species diversification. Thirteen of the 19 Araucaria species are endemic to the territory and form a monophyletic group, but patterns and processes that lead to such a high species richness are largely unexplored.

METHODS

We used 142 polymorphic AFLP markers and performed analyses based on Bayesian clustering algorithms, genetic distances, and cladistics on 71 samples representing all New Caledonian Araucaria species. We examined correlations between the inferred evolutionary relationships and shared morphological, ecological, or geographic parameters among species, to investigate evolutionary processes that may have driven speciation.

KEY RESULTS

We showed that genetic divergence among the present New Caledonian Araucaria species is low, suggesting recent diversification rather than pre-existence on Gondwana. We identified three genetic groups that included small-leaved, large-leaved, and coastal species, but detected no association with soil preference, ecological habitat, or rainfall. The observed patterns suggested that speciation events resulted from both differential adaptation and vicariance. Last, we hypothesize that speciation is ongoing and/or there are cryptic species in some genetically (sometimes also morphologically) divergent populations.

CONCLUSIONS

Further data are required to provide better resolution and understanding of the diversification of New Caledonian Araucaria species. Nevertheless, our study allowed insights into their evolutionary relationships and provides a framework for future investigations on the evolution of this emblematic group of plants in one of the world's biodiversity hotspots.

The relative contribution of band number to phylogenetic accuracy in AFLP data sets

We examined the effect of increasing the number of sampled amplified fragment length polymorphism (AFLP) bands to reconstruct an accurate and well-supported AFLP-based phylogeny. In silico AFLP was performed using simulated DNA sequences evolving along balanced and unbalanced model trees with recent, uniform and ancient radiations and average branch lengths (from the most internal node to the tip) ranging from 0.02 to 0.05 substitutions per site. Trees were estimated by minimum evolution (ME) and maximum parsimony (MP) methods from both DNA sequences and virtual AFLP fingerprints. The comparison of the true tree with the estimated AFLP trees suggests that moderate numbers of AFLP bands are necessary to recover the correct topology with high bootstrap support values (i.e. >70%). Fewer numbers of bands are necessary for shorter tree lengths and for balanced than for unbalanced tree topologies. However, branch length estimation was rather unreliable and did not improve substantially after a certain number of bands were sampled. These results hold for different levels of genome coverage and number of taxa analysed. In silico AFLP using bacterial genomic DNA sequences recovered a well-supported tree topology that mirrored an empirical phylogeny based on a set of 31 orthologous gene sequences when as few as 263 AFLP bands were scored. These results suggest that AFLPs may be an efficient alternative to traditional DNA sequencing for accurate topology reconstruction of shallow trees when not very short ancestral branches exist.

Polyphasic analysis of intraspecific diversity in Epicoccum nigrum warrants reclassification into separate species

BACKGROUND

Epicoccum nigrum Link (syn. E. purpurascens Ehrenb. ex Schlecht) is a saprophytic ascomycete distributed worldwide which colonizes a myriad of substrates. This fungus has been known as a biological control agent for plant pathogens and produces a variety of secondary metabolites with important biological activities as well as biotechnological application. E. nigrum produces darkly pigmented muriform conidia on short conidiophores on sporodochia and is a genotypically and phenotypically highly variable species. Since different isolates identified as E. nigrum have been evaluated as biological control agents and used for biocompound production, it is highly desirable that this species name refers to only one lineage. However, according to morphological and genetic variation, E. nigrum present two genotypes that may comprise more than one species.

METHODOLOGY/PRINCIPAL FINDINGS

We report the application of combined molecular (ITS and β-tubulin gene sequence analysis, PCR-RFLP and AFLP techniques), morphometric, physiological, genetic compatibility and recombination analysis to study the taxonomic relationships within an endophytic population that has been identified as E. nigrum. This combined analysis established two genotypes showing morphological, physiological and genetic divergence as well as genetic incompatibility characterized by colony inhibition, strongly indicating that these genotypes correspond to different species. Genotype 1 corresponds to E. nigrum while genotype 2 represents a new species, referred to in this study as Epicoccum sp.

CONCLUSIONS/SIGNIFICANCE

This research contributes to the knowledge of the Epicoccum genus and asserts that the classification of E. nigrum as a single variable species should be reassessed. In fact, based on the polyphasic approach we suggest the occurrence of cryptic species within E. nigrum and also that many of the sequences deposited as E. nigrum in GenBank and culture collection of microbial strains should be reclassified, including the reference strain CBS 161.73 sequenced in this work. In addition, this study provides valuable tools for differentiation of Epicoccum species.

Phylogenetic relationships in the genus Leontopodium (Asteraceae: Gnaphalieae) based on AFLP data

The genus Leontopodium comprises 30-41 species. The centre of diversity is the Sino-Himalayan region in south-western China, where about 15 species occur. The two species native to Europe, L. alpinum (known as the common 'Edelweiss') and L. nivale, are part of the cultural heritage of the people living there. Despite its importance, very little is known about the systematics of the genus. Because recent molecular studies have shown that species within this genus are closely related and difficult to distinguish with rDNA and cpDNA data, we used AFLPs to obtain a more detailed understanding of the phylogeny of the genus. Our main aims were as follows: (1) to clarify species relationships within the genus; and (2) to reveal information about the biogeography of the genus. We used AFLPs with six primer combinations to investigate 216 individuals in 38 populations of 16 different species. With AFLPs, we were able to recognize 10 different groups, all of which had strong bootstrap support. These results were also congruent with the morphology-based taxonomy of the genus. Most private and rare fragments were found in the Yunnan region (south-western China) relative to Europe and Mongolia/central China, suggesting a long-lasting in situ history of populations in the centre of diversity of the genus. Our results illustrate the utility of AFLPs to resolve phylogenetic relationships between these closely related species.

On the origin of Solanum nigrum: can networks help?

Black nightshades are a group of species best known for their 'poisonous' or noxious weedy reputation. It is not so well known that species of this group serve as emerging food source in many countries worldwide especially in the African continent. Despite the fact that the section has recently been studied extensively, taxonomy is still unsettled and debated because of inter- and intraspecific hybridization, phenotypic plasticity and polyploidization. In this study we analyze the genetic relationships among diploid, tetraploid and hexaploid species of sect. Solanum, which have possibly taken part in the formation of Solanum nigrum, utilizing multi-locus (SCoT, ISSR) markers combined with chloroplast trnL-F sequence data and morphological characters. We scored 51 morphological characters united with SCoT (171), ISSR (224) and trnL-F (1042), for simultaneous analysis of 49 terminals and 1488 characters. The topology of the tree is concordant with the results of the network analysis. In the phylogenetic networks, all the accessions of the diploid species shared a split with all of the polyploid species. This reflected a high portion of shared ISSR and SCoT bands between diploids and polyploids. In addition, a strong split divided the diploid species. The history of S. nigrum might be reticulate with hybrid speciation playing an important rule. Genetically differentiated diploids in few combinations have created a series of genetically distinct polyploid populations. The insufficient isolation that permitted further recombination between ancient polyploids and diploids have resulted in high level of genotypic and phenotypic polymorphism. This high level of novel genomic variability obviously enabled species to succeed in their new environment.

Species delimitation under the general lineage concept: an empirical example using wild North American hops (Cannabaceae: Humulus lupulus)

There is an emerging consensus that the intent of most species concepts is to identify evolutionarily distinct lineages. However, the criteria used to identify lineages differ among concepts depending on the perceived importance of various attributes of evolving populations. We have examined five different species criteria to ask whether the three taxonomic varieties of Humulus lupulus (hops) native to North America are distinct lineages. Three criteria (monophyly, absence of genetic intermediates, and diagnosability) focus on evolutionary patterns and two (intrinsic reproductive isolation and niche specialization) consider evolutionary processes. Phylogenetic analysis of amplified fragment length polymorphism (AFLP) data under a relaxed molecular clock, a stochastic Dollo substitution model, and parsimony identified all varieties as monophyletic, thus they satisfy the monophyly criterion for species delimitation. Principal coordinate analysis and a Bayesian assignment procedure revealed deep genetic subdivisions and little admixture between varieties, indicating an absence of genetic intermediates and compliance with the genotypic cluster species criterion. Diagnostic morphological and AFLP characters were found for all varieties, thus they meet the diagnosability criterion. Natural history information suggests that reproductive isolating barriers may have evolved in var. pubescens, potentially qualifying it as a species under a criterion of intrinsic reproductive isolation. Environmental niche modeling showed that the preferred habitat of var. neomexicanus is climatically unique, suggesting niche specialization and thus compliance with an ecological species criterion. Isolation by distance coupled with imperfect sampling can lead to erroneous lineage identification using some species criteria. Compliance with complementary pattern- and process-oriented criteria provides powerful corroboration for a species hypothesis and mitigates the necessity for comprehensive sampling of the entire species range, a practical impossibility in many systems. We hypothesize that var. pubescens maintains its genetic identity, despite substantial niche overlap with var. lupuloides, via the evolution of partial reproductive isolating mechanisms. Variety neomexicanus, conversely, will likely persist as a distinct lineage, regardless of limited gene flow with vars. lupuloides and pubescens because of ecological isolation--adaptation to the unique conditions of the Rocky Mountain cordillera. Thus, we support recognition of vars. neomexicanus and pubescens as species, but delay making a recommendation for var. lupuloides until sampling of genetic variation is complete or a stable biological process can be identified to explain its observed genetic divergence.

Genomic scans support repetitive continental colonization events during the rapid radiation of voles (Rodentia: Microtus): the utility of AFLPs versus mitochondrial and nuclear sequence markers

Single locus studies might not resolve phylogenetic relationships and the evolutionary history of taxa. The analysis of multiple markers promises higher resolution, and congruence among loci may indicate that the phylogenies represent the underlying species history. Here, we examine the utility of a genome-wide approach based on amplified fragment length polymorphisms (AFLP) and several DNA sequence markers in resolving phylogenetic signals in the rapidly radiating rodent genus Microtus which produced about 70 vole species within the last 1.2-2 myr. The current Holarctic distribution of Microtus is assumed to have resulted from three independent colonization events out of Asia to North America, Europe, and northern Asia without subsequent colonization, which would have led to deep splits between species from different continents. We investigated this hypothesis of three single colonization events by reconstructing the phylogenetic relationships among species from all three continents based on data from the first exon of the nuclear arginine vasopressin receptor 1a gene (EXON1), an adjacent noncoding region and the mitochondrial cytochrome b gene. The phylogenetic patterns obtained from these sequence markers are contrasted to genome-wide data on more than 1800 amplified fragment length polymorphisms (AFLP) analyzed for the same samples. Our results show that the single sequence markers partially resolve the phylogenetic relationships within Microtus, but with some incongruence mostly between EXON1 and the other loci. However, deeper nodes of the radiation are only weakly supported and neither the combination of the markers nor additional nuclear sequences improved the resolution significantly. AFLPs provided much stronger support for major continent-specific clades, and show also that reciprocal monophyly of American and European voles is incomplete. Our results demonstrate that Microtus voles colonized the American and European continents each repeatedly in several independent events on similar colonization routes during their radiation. More generally, this study supports the suitability of AFLPs as an alternative to sequence markers to resolve the evolutionary history of rapidly radiating taxa.

Phylogeography of 3 North Atlantic Wolffish species (Anarhichas spp.) with phylogenetic relationships within the family Anarhichadidae

Phylogenetic analyses of all 4 wolffish species (Atlantic, Spotted, Northern, and Bering wolffishes) and the Wolfeel were assessed with both mitochondrial (D-loop and ND1) and nuclear (amplified fragment length polymorphism) DNA to resolve relationships within the family Anarhichadidae. Species-specific mitochondrial DNA (mtDNA) mutation rates were estimated based on 2 possible dates of divergence between the Pacific and Atlantic lineages. Phylogeographic patterns within each of the 3 North Atlantic wolffishes were investigated with Markov chain Monte Carlo simulations based on mtDNA to determine whether population size changes occurred following the last glaciation and where wolffishes likely survived glaciation. All 3 species of North Atlantic wolffishes showed evidence of postglacial expansion but did not show evidence of persistence in multiple refugia in both the eastern and western Atlantic Ocean. Rather, the data supported persistence in a single refuge, with postglacial expansion into the rest of the range. Nucleotide diversity, in particular, was low in wolffishes compared with other marine fishes, possibly related to reductions in population sizes during the last glaciation.

Analyzing the relationship between sequence divergence and nodal support using Bayesian phylogenetic analyses

Determining the appropriate gene for phylogeny reconstruction can be a difficult process. Rapidly evolving genes tend to resolve recent relationships, but suffer from alignment issues and increased homoplasy among distantly related species. Conversely, slowly evolving genes generally perform best for deeper relationships, but lack sufficient variation to resolve recent relationships. We determine the relationship between sequence divergence and Bayesian phylogenetic reconstruction ability using both natural and simulated datasets. The natural data are based on 28 well-supported relationships within the subphylum Vertebrata. Sequences of 12 genes were acquired and Bayesian analyses were used to determine phylogenetic support for correct relationships. Simulated datasets were designed to determine whether an optimal range of sequence divergence exists across extreme phylogenetic conditions. Across all genes we found that an optimal range of divergence for resolving the correct relationships does exist, although this level of divergence expectedly depends on the distance metric. Simulated datasets show that an optimal range of sequence divergence exists across diverse topologies and models of evolution. We determine that a simple to measure property of genetic sequences (genetic distance) is related to phylogenic reconstruction ability in Bayesian analyses. This information should be useful for selecting the most informative gene to resolve any relationships, especially those that are difficult to resolve, as well as minimizing both cost and confounding information during project design.

Dissecting comimetic radiations in Heliconius reveals divergent histories of convergent butterflies

Mimicry among Heliconius butterflies provides a classic example of coevolution but unresolved relationships among mimetic subspecies have prevented examination of codiversification between species. We present amplified fragment length polymorphism and mtDNA datasets for the major comimetic races of Heliconius erato and H. melpomene. The AFLP data reveal unprecedented resolution, clustering samples by geography and race in both species. Our results show that, although H. erato and H. melpomene co-occur, mimic each other, and exhibit parallel shifts in color pattern, they experienced very different modes of diversification and geographic histories. Our results suggest that H. erato originated on the western side of South America whereas H. melpomene originated in the east. H. erato underwent rapid diversification and expansion with continued gene-flow following diversification, resulting in widely dispersed sister taxa. In contrast, H. melpomene underwent a slower pace of diversification with lower levels of gene flow, producing a stepwise directional expansion from west to east. Our results also suggest that each of the three main wing pattern phenotypes originated and/or was lost multiple times in each species. The rayed pattern is likely to be the ancestral phenotype in H. erato whereas postman or red patch is likely to be ancestral in H. melpomene. Finally, H. cydno and H. himera are monophyletic entities clearly nested within H. melpomene and H. erato, rather than being their respective sister species. Estimates of mtDNA divergence suggest a minimum age of 2.8 and 2.1 My for H. erato and H. melpomene, respectively, placing their origins in the late Pliocene.

Phylogenetic relationships and natural hybridization among the North American woody bamboos (Poaceae: Bambusoideae: Arundinaria)

In spite of the ecological and economic importance of temperate bamboos, relatively little is known about their population biology or evolutionary history. Recently, hybridization has emerged as a potential source of diversity in this group, as well as an underlying cause of taxonomic problems. As part of a broader phylogenetic study of the temperate bamboos, we report the results of an analysis of the North American Arundinaria gigantea species complex, including estimates of genetic variation and molecular evidence of natural hybridization among A. gigantea, A. tecta, and A. appalachiana. The study involved a comparative analysis of amplified fragment length polymorphisms (AFLPs) and chloroplast DNA sequences representing diversity within and among all three species plus individuals with intermediate or unusual morphological characteristics (putative hybrids). Molecular results support the recognition of three species previously defined on the basis of morphology, anatomy, and ecology, with most of the molecular variance accounted for by among-species variation. Molecular evidence also demonstrates that A. tecta and A. appalachiana are sister species, forming a clade that is significantly divergent from A. gigantea. The role of hybridization in the phylogenetic history of Arundinaria is discussed along with implications for the evolution and taxonomy of the temperate woody bamboos.

Evaluating the relationship between evolutionary divergence and phylogenetic accuracy in AFLP data sets

Using in silico amplified fragment length polymorphism (AFLP) fingerprints, we explore the relationship between sequence similarity and phylogeny accuracy to test when, in terms of genetic divergence, the quality of AFLP data becomes too low to be informative for a reliable phylogenetic reconstruction. We generated DNA sequences with known phylogenies using balanced and unbalanced trees with recent, uniform and ancient radiations, and average branch lengths (from the most internal node to the tip) ranging from 0.02 to 0.4 substitutions per site. The resulting sequences were used to emulate the AFLP procedure. Trees were estimated by maximum parsimony (MP), neighbor-joining (NJ), and minimum evolution (ME) methods from both DNA sequences and virtual AFLP fingerprints. The estimated trees were compared with the reference trees using a score that measures overall differences in both topology and relative branch length. As expected, the accuracy of AFLP-based phylogenies decreased dramatically in the more divergent data sets. Above a divergence of approximately 0.05, AFLP-based phylogenies were largely inaccurate irrespective of the distinct topology, radiation model, or phylogenetic method used. This value represents an upper bound of expected tree accuracy for data sets with a simple divergence history; AFLP data sets with a similar divergence but with unbalanced topologies and short ancestral branches produced much less accurate trees. The lack of homology of AFLP bands quickly increases with divergence and reaches its maximum value (100%) at a divergence of only 0.4. Low guanine-cytosine (GC) contents increase the number of nonhomologous bands in AFLP data sets and lead to less reliable trees. However, the effect of the lack of band homology on tree accuracy is surprisingly small relative to the negative impact due to the low information content of AFLP characters. Tree-building methods based on genetic distance displayed similar trends and outperformed parsimony at low but not at high divergences. However, the impact of using alternative phylogenetic methods on tree accuracy was generally small relative to the uncertainty arising from factors such as divergence, nonhomology of bands, or the low information content of AFLP characters. Nevertheless, our data suggest that under certain circumstances, AFLPs may be suitable to reconstruct deeper phylogenies than usually accepted.

Testing mitochondrial sequences and anonymous nuclear markers for phylogeny reconstruction in a rapidly radiating group: molecular systematics of the Delphininae (Cetacea: Odontoceti: Delphinidae)

BACKGROUND

Many molecular phylogenetic analyses rely on DNA sequence data obtained from single or multiple loci, particularly mitochondrial DNA loci. However, phylogenies for taxa that have undergone recent, rapid radiation events often remain unresolved. Alternative methodologies for discerning evolutionary relationships under these conditions are desirable. The dolphin subfamily Delphininae is a group that has likely resulted from a recent and rapid radiation. Despite several efforts, the evolutionary relationships among the species in the subfamily remain unclear.

RESULTS

Here, we compare a phylogeny estimated using mitochondrial DNA (mtDNA) control region sequences to a multi-locus phylogeny inferred from 418 polymorphic genomic markers obtained from amplified fragment length polymorphism (AFLP) analysis. The two sets of phylogenies are largely incongruent, primarily because the mtDNA tree provides very poor resolving power; very few species' nodes in the tree are supported by bootstrap resampling. The AFLP phylogeny is considerably better resolved and more congruent with relationships inferred from morphological data. Both phylogenies support paraphyly for the genera Stenella and Tursiops. The AFLP data indicate a close relationship between the two spotted dolphin species and recent ancestry between Stenella clymene and S. longirostris. The placement of the Lagenodelphis hosei lineage is ambiguous: phenetic analysis of the AFLP data is consistent with morphological expectations but the phylogenetic analysis is not.

CONCLUSION

For closely related, recently diverged taxa, a multi-locus genome-wide survey is likely the most comprehensive approach currently available for phylogenetic inference.

Linkage disequilibrium and population structure in wild and domesticated populations of Phaseolus vulgaris L

Together with the knowledge of the population structure, a critical aspect for the planning of association and/or population genomics studies is the level of linkage disequilibrium (LD) that characterizes the species and the population used for such an analysis. We have analyzed the population structure and LD in wild and domesticated populations of Phaseolus vulgaris L. using amplified fragment length polymorphism markers, most of which were genetically mapped in two recombinant inbred populations. Our results reflect the previous knowledge of the occurrence of two major wild gene pools of P. vulgaris, from which two independent domestication events originated, one in the Andes and one in Mesoamerica. The high level of LD in the whole sample was mostly due to the gene pool structure, with a much higher LD in domesticated compared to wild populations. In relation to association studies, our results also suggest that whole-genome-scan approaches are feasible in the common bean. Interestingly, an excess of inter-chromosomal LD was found in the domesticated populations, which suggests an important role for epistatic selection during domestication. Moreover, our results indicate the occurrence of a strong bottleneck in the Andean wild population before domestication, suggesting a Mesoamerican origin of P. vulgaris. Finally, our data support the occurrence of a single domestication event in Mesoamerica, and the same scenario in the Andes.

Reconstructing the species phylogeny of Pseudopanax (Araliaceae), a genus of hybridising trees

Pseudopanax (Araliaceae) comprises 12 tree species of diverse morphology and ecology, and is endemic to New Zealand. It is notable for the hybridisation that occurs between P. crassifolius and P. lessonii, which have very different leaves and habits. To provide context for the study of this hybridisation and other investigations, we examined the phylogeny of Pseudopanax using chloroplast DNA sequences (c.5900 base-pairs) and AFLP DNA-fingerprinting. Both approaches resolve two principal groups within Pseudopanax--the Arboreus group and the Crassifolius+Lessonii union--but how they are related to other genera remains unclear. There is, nevertheless, little compelling evidence against the monophyly of Pseudopanax, making unnecessary the recent re-segregation of the Arboreus group as Neopanax. The chloroplast data provided minimal additional resolution, although the position of P. linearis was discordant compared to other data. Analyses of the AFLP data strongly recovered each species, aside from the morphologically heterogeneous P. colensoi, and the two mainland species (P. arboreus and P. crassifolius) that each contained a nested island-endemic (P. kermadecensis and P. chathamicus, respectively). However, relationships amongst species within the two principal groups were poorly resolved. An example was the uncertainty of whether P. crassifolius grouped with P. lessonii and its allies, or the morphologically similar species it had been previously placed with. This in turn raises the issue of how hybridisation might affect phylogenies and the ability to reconstruct them, even when using multiple, independent markers.