Using SINEs to probe ancient explosive speciation: "hidden" radiation of African cichlids?

Cichlid fishes of the east African Great Lakes represent a paradigm of adaptive radiation. We conducted a phylogenetic analysis of cichlids including pan-African and west African species by using insertion patterns of short interspersed elements (SINEs) at orthologous loci. The monophyly of the east African cichlids was consistently supported by seven independent insertions of SINE sequences that are uniquely shared by these species. In addition, data from four other loci indicated that the genera Tilapia (pan-African) and Steatocranus (west African) are the closest relatives to east African cichlids. However, relationships among Tilapia, Steatocranus, and the east African clade were ambiguous because of incongruencies among topologies suggested by insertion patterns of SINEs at six other loci. One plausible explanation for this phenomenon is incomplete lineage sorting of alleles containing or missing a SINE insertion at these loci during ancestral speciation. Such incomplete sorting may have taken place earlier than 14 MYA, followed by random and stochastic fixation of the alleles in subsequent lineages. These observations prompted us to consider the possibility that cichlid speciation occurred at an accelerated rate during this period when the African Great Lakes did not exist. The SINE method could be useful for detecting ancient exclusive speciation events that tend to remain hidden during conventional sequence analyses because of accumulated point mutations.     

Retroposition of the AFC Family of SINEs (Short Interspersed Repetitive Elements) Before and During the Adaptive Radiation of Cichlid Fishes in Lake Malawi and Related Inferences About Phylogeny

Lake Malawi is home to more than 450 species of endemic cichlids, which provide a spectacular example of adaptive radiation. To clarify the phylogenetic relationships among these fish, we examined the presence and absence of SINEs (short interspersed repetitive elements) at orthologous loci. We identified six loci at which a SINE sequence had apparently been specifically inserted by retroposition in the common ancestor of all the investigated species of endemic cichlids in Lake Malawi. At another locus, unique sharing of a SINE sequence was evident among all the investigated species of endemic non-Mbuna cichlids with the exception of Rhamphochromis sp. The relationships were in good agreement with those deduced in previous studies with various different markers, demonstrating that the SINE method is useful for the elucidation of phylogenetic relationships among cichlids in Lake Malawi. We also characterized a locus that exhibited transspecies polymorphism with respect to the presence or absence of the SINE sequence among non-Mbuna species. This result suggests that incomplete lineage sorting and/or interspecific hybridization might have occurred or be occurring among the species in this group, which might potentially cause misinterpretation of phylogenetic data, in particular when a single-locus marker, such as a sequence in the mitochondrial DNA, is used for analysis. Received: 15 December 2000 / Accepted: 30 January 2001     

SINEs of speciation: tracking lineages with retroposons

The value of short interspersed elements (SINEs) for diagnosing common ancestry is being expanded to examine the differential sorting of lineages through the course of speciation events. Because most SINEs are neutral markers of identical descent, are not precisely excised from the genome and have a known ancestral condition, they are advantageous for reconciling gene trees and species trees with minimal phylogenetic error. A population perspective on SINE evolution combined with coalescence theory provides a context for investigating the phenomenon of ancestral polymorphism and its role in producing incongruent SINE insertion patterns among multiple loci. Studies of human Alu repeats demonstrate the value of young polymorphic SINEs for assessing human genomic diversity and tracking ancient demographics of human populations, whereas incongruent insertion patterns revealed by older fixed SINE loci, such as those in African cichlid fishes, contain information that might help identify ancient radiations that are otherwise obscured by accumulated mutations in sequence data. Here, we review the utility of retroposons for inferring common ancestry, discuss limits to the method, and clarify confusion by providing examples from the literature that illustrate how discordant multi-locus insertion patterns of retroelements can indicate lineage-sorting events that should not be misinterpreted as phylogenetic noise.     

Extensive morphological convergence and rapid radiation in the evolutionary history of the family Geoemydidae (old world pond turtles) revealed by SINE insertion analysis

The family Geoemydidae is one of three in the superfamily Testudinoidea and is the most diversified family of extant turtle species. The phylogenetic relationships in this family and among related families have been vigorously investigated from both morphological and molecular viewpoints. The evolutionary history of Geoemydidae, however, remains controversial. Therefore, to elucidate the phylogenetic relationships of Geoemydidae and related species, we applied the SINE insertion method to investigate 49 informative SINE loci in 28 species. We detected four major evolutionary lineages (Testudinidae, Batagur group, Siebenrockiella group, and Geoemyda group) in the clade Testuguria (a clade of Geoemydidae + Testudinidae). All five specimens of Testudinidae form a monophyletic clade. The Batagur group comprises five batagurines. The Siebenrockiella group has one species, Siebenrockiella crassicollis. The Geoemyda group comprises 15 geoemydines (including three former batagurines, Mauremys reevesii, Mauremys sinensis, and Heosemys annandalii). Among these four groups, the SINE insertion patterns were inconsistent at four loci, suggesting that an ancestral species of Testuguria radiated and rapidly diverged into the four lineages during the initial stage of its evolution. Furthermore, within the Geoemyda group we identified three evolutionary lineages, namely Mauremys, Cuora, and Heosemys. The Heosemys lineage comprises Heosemys, Sacalia, Notochelys, and Melanochelys species, and its monophyly is a novel assemblage in Geoemydidae. Our SINE phylogenetic tree demonstrates extensive convergent morphological evolution between the Batagur group and the three species of the Geoemyda group, M. reevesii, M. sinensis, and H. annandalii.     

Inferring phylogeny despite incomplete lineage sorting

It is now well known that incomplete lineage sorting can cause serious difficulties for phylogenetic inference, but little attention has been paid to methods that attempt to overcome these difficulties by explicitly considering the processes that produce them. Here we explore approaches to phylogenetic inference designed to consider retention and sorting of ancestral polymorphism. We examine how the reconstructability of a species (or population) phylogeny is affected by (a) the number of loci used to estimate the phylogeny and (b) the number of individuals sampled per species. Even in difficult cases with considerable incomplete lineage sorting (times between divergences less than 1 N(e) generations), we found the reconstructed species trees matched the "true" species trees in at least three out of five partitions, as long as a reasonable number of individuals per species were sampled. We also studied the tradeoff between sampling more loci versus more individuals. Although increasing the number of loci gives more accurate trees for a given sampling effort with deeper species trees (e.g., total depth of 10 N(e) generations), sampling more individuals often gives better results than sampling more loci with shallower species trees (e.g., depth = 1 N(e)). Taken together, these results demonstrate that gene sequences retain enough signal to achieve an accurate estimate of phylogeny despite widespread incomplete lineage sorting. Continued improvement in our methods to reconstruct phylogeny near the species level will require a shift to a compound model that considers not only nucleotide or character state substitutions, but also the population genetics processes of lineage sorting. [Coalescence; divergence; population; speciation.].     

Proceedings of the SMBE Tri-National Young Investigators' Workshop 2005. Baleen whale phylogeny and a past extensive radiation event revealed by SINE insertion analysis

Baleen whales (suborder Mysticeti) comprise 11 extant species that are classified into four families. Although several phylogenetic hypotheses about these taxa have been proposed, their phylogenetic relationships remain confused. We addressed this problem using short interspersed repetitive element (SINE) insertion data, which now are regarded as almost ideal shared, derived characters at the molecular level. We reconstructed the phylogenetic relationships of baleen whales by characterizing 36 informative SINE loci. One of the intriguing conclusions is that balaenopterids and eschrichtiids radiated very rapidly during a very short evolutionary period. During this period, speciation occurred in balaenopterids and eschrichtiids while newly inserted SINE loci remains polymorphic. Later on, these SINEs were sorted incompletely into each lineage. Thus, there are now inconsistencies among species regarding the presence or absence of a given SINE. This is in sharp contrast to the phylogeny of toothed whales, for which no SINE inconsistencies have been found. Furthermore, we found monophyletic groupings between humpback and fin whales as well as between (sei+Bryde's) whales and blue whales, both of which have not previously been recognized. The comprehensive SINE insertion data, together with the mitochondrial DNA phylogeny that was recently completed (Sasaki, T., M. Nikaido, H. Healy et al. 2005. Mitochondrial phylogenetics and evolution of mysticete whales. Syst. Biol. 56:77-90; Rychel, A. L., T. W. Reeder, and A. Berta. 2004. Phylogeny of mysticete whales based on mitochondrial and nuclear data. Mol. Phylogenet. Evol. 32:892-901), provide a nearly complete picture of the evolutionary history of baleen whales.     

Phylogenetic, morphological, and chemotaxonomic incongruence in the North American endemic genus Echinacea

The study of recently formed species is important because it can help us to better understand organismal divergence and the speciation process. However, these species often present difficult challenges in the field of molecular phylogenetics because the processes that drive molecular divergence can lag behind phenotypic divergence. In the current study we show that species of the recently diverged North American endemic genus of purple coneflower, Echinacea, have low levels of molecular divergence. Data from three nuclear loci and two plastid loci provide neither resolved topologies nor congruent hypotheses about species-level relationships. This lack of phylogenetic resolution is likely due to the combined effects of incomplete lineage sorting, hybridization, and backcrossing following secondary contact. The poor resolution provided by molecular markers contrasts previous studies that found well-resolved and taxonomically supported relationships from metabolic and morphological data. These results suggest that phenotypic canalization, resulting in identifiable morphological species, has occurred rapidly within Echinacea. Conversely, molecular signals have been distorted by gene flow and incomplete lineage sorting. Here we explore the impact of natural history on the genetic organization and phylogenetic relationships of Echinacea.     

Incomplete Lineage Sorting and Hybridization Statistics for Large-Scale Retroposon Insertion Data

Ancient retroposon insertions can be used as virtually homoplasy-free markers to reconstruct the phylogenetic history of species. Inherited, orthologous insertions in related species offer reliable signals of a common origin of the given species. One prerequisite for such a phylogenetically informative insertion is that the inserted element was fixed in the ancestral population before speciation; if not, polymorphically inserted elements may lead to random distributions of presence/absence states during speciation and possibly to apparently conflicting reconstructions of their ancestry. Fortunately, such misleading fixed cases are relatively rare but nevertheless, need to be considered. Here, we present novel, comprehensive statistical models applicable for (1) analyzing any pattern of rare genomic changes, (2) testing and differentiating conflicting phylogenetic reconstructions based on rare genomic changes caused by incomplete lineage sorting or/and ancestral hybridization, and (3) differentiating between search strategies involving genome information from one or several lineages. When the new statistics are applied, in non-conflicting cases a minimum of three elements present in both of two species and absent in a third group are considered significant support (p<0.05) for the branching of the third from the other two, if all three of the given species are screened equally for genome or experimental data. Five elements are necessary for significant support (p<0.05) if a diagnostic locus derived from only one of three species is screened, and no conflicting markers are detected. Most potentially conflicting patterns can be evaluated for their significance and ancestral hybridization can be distinguished from incomplete lineage sorting by considering symmetric or asymmetric distribution of rare genomic changes among possible tree configurations. Additionally, we provide an R-application to make the new KKSC insertion significance test available for the scientific community at http://retrogenomics.uni-muenster.de:3838/KKSC_significance_test/.     

Retroposon insertion patterns of neoavian birds: strong evidence for an extensive incomplete lineage sorting era

More than 150 Ma, the avian lineage separated from that of other dinosaurs and later diversified into the more than 10,000 species extant today. The early neoavian bird radiations most likely occurred in the late Cretaceous (more than 65 Ma) but left behind few if any molecular signals of their archaic evolutionary past. Retroposed elements, once established in an ancestral population, are highly valuable, virtually homoplasy-free markers of species evolution; after applying stringent orthology criteria, their phylogenetically informative presence/absence patterns are free of random noise and independent of evolutionary rate or nucleotide composition effects. We screened for early neoavian orthologous retroposon insertions and identified six markers with conflicting presence/absence patterns, whereas six additional retroposons established before or after the presumed major neoavian radiation show consistent phylogenetic patterns. The exceptionally frequent conflicting retroposon presence/absence patterns of neoavian orders are strong indicators of an extensive incomplete lineage sorting era, potentially induced by an early rapid successive speciation of ancestral Neoaves.     

Species tree inference in a recent radiation of orioles (Genus Icterus): multiple markers and methods reveal cytonuclear discordance in the northern oriole group

Recent computational advances provide novel opportunities to infer species trees based on multiple independent loci. Thus, single gene trees no longer need suffice as proxies for species phylogenies. Several methods have been developed to deal with the challenges posed by incomplete and stochastic lineage sorting. In this study, we employed four Bayesian methods to infer the phylogeny of a clade of 11 recently diverged oriole species within the genus Icterus. We obtained well-resolved and mostly congruent phylogenies using a set of seven unlinked nuclear intron loci and sampling multiple individuals per species. Most notably, Bayesian concordance analysis generally agreed well with concatenation; the two methods agreed fully on eight of nine nodes. The coalescent-based method ^∗BEAST further supported six of these eight nodes. The fourth method used, BEST, failed to converge despite exhaustive efforts to optimize the tree search. Overall, the results obtained by new species tree methods and concatenation generally corroborate our findings from previous analyses and data sets. However, we found striking disagreement between mitochondrial and nuclear DNA involving relationships within the northern oriole group. Our results highlight the danger of reliance on mtDNA alone for phylogenetic inference. We demonstrate that in spite of low variability and incomplete lineage sorting, multiple nuclear loci can produce largely congruent phylogenies based on multiple species tree methods, even for very closely-related species.     

Shifting distributions and speciation: species divergence during rapid climate change

Questions about how shifting distributions contribute to species diversification remain virtually without answer, even though rapid climate change during the Pleistocene clearly impacted genetic variation within many species. One factor that has prevented this question from being adequately addressed is the lack of precision associated with estimates of species divergence made from a single genetic locus and without incorporating processes that are biologically important as populations diverge. Analysis of DNA sequences from multiple variable loci in a coalescent framework that (i) corrects for gene divergence pre-dating speciation, and (ii) derives divergence-time estimates without making a priori assumptions about the processes underlying patterns of incomplete lineage sorting between species (i.e. allows for the possibility of gene flow during speciation), is critical to overcoming the inherent logistical and analytical difficulties of inferring the timing and mode of speciation during the dynamic Pleistocene. Estimates of species divergence that ignore these processes, use single locus data, or do both can dramatically overestimate species divergence. For example, using a coalescent approach with data from six loci, the divergence between two species of montane Melanoplus grasshoppers is estimated at between 200,000 and 300,000 years before present, far more recently than divergence estimates made using single-locus data or without the incorporation of population-level processes. Melanoplus grasshoppers radiated in the sky islands of the Rocky Mountains, and the analysis of divergence between these species suggests that the isolation of populations in multiple glacial refugia was an important factor in promoting speciation. Furthermore, the low estimates of gene flow between the species indicate that reproductive isolation must have evolved rapidly for the incipient species boundaries to be maintained through the subsequent glacial periods and shifts in species distributions.     

Diversification of Hawaiian Cyrtandra (Gesneriaceae) under the influence of incomplete lineage sorting and hybridization

Cyrtandra (Gesneriaceae) is a genus of flowering plants with over 800 species distributed throughout Southeast Asia and the Pacific Islands. On the Hawaiian Islands, 60 named species and over 89 putative hybrids exist, most of which are identified on the basis of morphology. Despite many previous studies on the Hawaiian Cyrtandra lineage, questions regarding the reconciliation of morphology and genetics remain, many of which can be attributed to the relatively young age and evidence of hybridization between species. We utilized targeted enrichment, high‐throughput sequencing, and modern phylogenomics tools to test 31 Hawaiian Cyrtandra samples (22 species, two putative hybrids, four species with two samples each, one species with four samples) and two outgroups for species relationships and hybridization in the presence of incomplete lineage sorting (ILS). Both concatenated and species‐tree methods were used to reconstruct species relationships, and network analyses were conducted to test for hybridization. We expected to see high levels of ILS and putative hybrids intermediate to their parent species. Phylogenies reconstructed from the concatenated and species‐tree methods were highly incongruent, most likely due to high levels of incomplete lineage sorting. Network analyses inferred gene flow within this lineage, but not always between taxa that we expected. Multiple hybridizations were inferred, but many were on deeper branches of the island lineages suggesting a long history of hybridization. We demonstrated the utility of high‐throughput sequencing and a phylogenomic approach using 569 loci to understanding species relationships and gene flow in the presence of ILS.     

Genetic evidence on the demography of speciation in allopatric dolphin species

Under a neutral model, the stochastic lineage sorting that leads to gene monophyly proceeds slowly in large populations. Therefore, in many recent species with large population size, the genome will have mixed support for monophyly unless historical bottlenecks have accelerated coalescence. We use genealogical patterns in mitochondrial DNA and in introns of four nuclear loci to test for historical bottlenecks during the speciation and divergence of two temperate Lagenorhynchus dolphin species isolated by tropical Pacific waters (an antitropical distribution). Despite distinct morphologies, foraging behaviors, and mitochondrial DNAs, these dolphin species are polyphyletic at all four nuclear loci. The abundance of shared polymorphisms between these sister taxa is most consistent with the maintenance of large effective population sizes (5.09 x 10(4) to 10.9 x 10(4)) during 0.74-1.05 million years of divergence. A variety of population size histories are possible, however. We used gene tree coalescent probabilities to explore the rejection region for historical bottlenecks of different intensity given best estimates of effective population size under a strict isolation model of divergence. In L. obliquidens the data are incompatible with a colonization propagule of an effective size of 10 or fewer individuals. Although the ability to reject less extreme historical bottlenecks will require data from additional loci, the intermixed genealogical patterns observed between these dolphin sister species are highly probable only under an extended history of large population size. If similar demographic histories are inferred for other marine antitropical taxa, a parsimonious model for the Pleistocene origin of these distributions would not involve rare breaches of a constant dispersal barrier by small colonization propagules. Instead, a history of large population size in L. obliquidens and L. obscurus contributes to growing biological and environmental evidence that the equatorial barrier became permeable during glacial/interglacial cycles, leading to vicariant isolation of antitropical populations.     

Birds perching on bushes: Networks to visualize conflicting phylogenetic signals during early avian radiation

Hybridization is increasingly seen as an important source of adaptive genetic variation and biotic diversity. Recent phylogenetic studies on the early evolution of birds suggest that the early diversification of neoavian orders perhaps involved a period of extensive hybridization or incomplete lineage sorting. Phylogenetic error, saturation, long-branch attraction, and convergence make it difficult to detect ancient hybridization events and differentiate them from incomplete lineage sorting using sequence data. We used recently published retroposon marker data to visualize the early radiation of Neoaves within a phylogenetic network approach, and found that the most basal neoavian taxa indeed show a complex pattern of reticulated relationships. Moreover, the reticulation levels of different parts of the network are consistent with the insertion pattern of the retroposon elements. The use of network-based analyses on homoplasy-free data shows true conflicting signals and the taxa involved that are not represented in trees.     

Incomplete lineage sorting impacts the inference of macroevolutionary regimes from molecular phylogenies when concatenation is employed: An analysis based on Cetacea

Interest in methods that estimate speciation and extinction rates from molecular phylogenies has increased over the last decade. The application of such methods requires reliable estimates of tree topology and node ages, which are frequently obtained using standard phylogenetic inference combining concatenated loci and molecular dating. However, this practice disregards population‐level processes that generate gene tree/species tree discordance. We evaluated the impact of employing concatenation and coalescent‐based phylogeny inference in recovering the correct macroevolutionary regime using simulated data based on the well‐established diversification rate shift of delphinids in Cetacea. We found that under scenarios of strong incomplete lineage sorting, macroevolutionary analysis of phylogenies inferred by concatenating loci failed to recover the delphinid diversification shift, while the coalescent‐based tree consistently retrieved the correct rate regime. We suggest that ignoring microevolutionary processes reduces the power of methods that estimate macroevolutionary regimes from molecular data.     

Recurrent introgression of mitochondrial DNA among hares (Lepus spp.) revealed by species-tree inference and coalescent simulations

Understanding recent speciation history requires merging phylogenetic and population genetics approaches, taking into account the persistence of ancestral polymorphism and possible introgression. The emergence of a clear phylogeny of hares (genus Lepus) has been hampered by poor genomic sampling and possible occurrence of mitochondrial DNA (mtDNA) introgression from the arctic/boreal Lepus timidus into several European temperate and possibly American boreal species. However, no formal test of introgression, taking also incomplete lineage sorting into account, has been done. Here, to clarify the yet poorly resolved species phylogeny of hares and test hypotheses of mtDNA introgression, we sequenced 14 nuclear DNA and 2 mtDNA fragments (8205 and 1113 bp, respectively) in 50 specimens from 11 hare species from Eurasia, North America, and Africa. By applying an isolation-with-migration model to the nuclear data on subsets of species, we find evidence for very limited gene flow from L. timidus into most temperate European species, and not into the American boreal ones. Using a multilocus coalescent-based method, we infer the species phylogeny, which we find highly incongruent with mtDNA phylogeny using parametric bootstrap. Simulations of mtDNA evolution under the speciation history inferred from nuclear genes did not support the hypothesis of mtDNA introgression from L. timidus into the American L. townsendii but did suggest introgression from L. timidus into 4 temperate European species. One such event likely resulted in the complete replacement of the aboriginal mtDNA of L. castroviejoi and of its sister species L. corsicanus. It is remarkable that mtDNA introgression in hares is frequent, extensive, and always from the same donor arctic species. We discuss possible explanations for the phenomenon in relation to the dynamics of range expansions and species replacements during the climatic oscillations of the Pleistocene.     

Coalescent histories on phylogenetic networks and detection of hybridization despite incomplete lineage sorting

Analyses of the increasingly available genomic data continue to reveal the extent of hybridization and its role in the evolutionary diversification of various groups of species. We show, through extensive coalescent-based simulations of multilocus data sets on phylogenetic networks, how divergence times before and after hybridization events can result in incomplete lineage sorting with gene tree incongruence signatures identical to those exhibited by hybridization. Evolutionary analysis of such data under the assumption of a species tree model can miss all hybridization events, whereas analysis under the assumption of a species network model would grossly overestimate hybridization events. These issues necessitate a paradigm shift in evolutionary analysis under these scenarios, from a model that assumes a priori a single source of gene tree incongruence to one that integrates multiple sources in a unifying framework. We propose a framework of coalescence within the branches of a phylogenetic network and show how this framework can be used to detect hybridization despite incomplete lineage sorting. We apply the model to simulated data and show that the signature of hybridization can be revealed as long as the interval between the divergence times of the species involved in hybridization is not too small. We reanalyze a data set of 106 loci from 7 in-group Saccharomyces species for which a species tree with no hybridization has been reported in the literature. Our analysis supports the hypothesis that hybridization occurred during the evolution of this group, explaining a large amount of the incongruence in the data. Our findings show that an integrative approach to gene tree incongruence and its reconciliation is needed. Our framework will help in systematically analyzing genomic data for the occurrence of hybridization and elucidating its evolutionary role.     

First application of the SINE (short interspersed repetitive element) method to infer phylogenetic relationships in reptiles: an example from the turtle superfamily Testudinoidea

Although turtles (order Testudines) constitute one of the major reptile groups, their phylogenetic relationships remain largely unresolved. Hence, we attempted to elucidate their phylogeny using the SINE (short interspersed repetitive element) method, in which the sharing of a SINE at orthologous loci is indicative of synapomorphy. First, a detailed characterization of the tortoise polIII/SINE was conducted using 10 species from eight families of hidden-necked turtles (suborder Cryptodira). Our analysis of 382 SINE sequences newly isolated in the present study revealed two subgroups, namely Cry I and Cry II, which were distinguishable according to diagnostic nucleotides in the 3' region. Furthermore, four (IA-ID) and five (IIA-IIE) different SINE types were identified within Cry I and Cry II subgroups, respectively, based on features of insertions/deletions located in the middle of the SINE sequences. The relative frequency of occurrence of the subgroups and the types of SINEs in this family were highly variable among different lineages of turtles, suggesting active differential retroposition in each lineage. Further application of the SINE method using the most retrotranspositionally active types, namely IB and IC, challenged the established phylogenetic relationships of Bataguridae and its related families. The data for 11 orthologous loci demonstrated a close relationship between Bataguridae and Testudinidae, as well as the presence of the three clades within Bataguridae. Although the SINE method has been used to infer the phylogenies of a number of vertebrate groups, it has never been applied to reptiles. The present study represents the first application of this method to a phylogenetic analysis of this class of vertebrates, and it provides detailed information on the SINE subgroups and types. This information may be applied to the phylogenetic resolution of relevant turtle lineages.     

Contrasting patterns of genetic structure in two species of the coral trout Plectropomus (Serranidae) from east and west Australia: introgressive hybridisation or ancestral polymorphisms

Inter-specific genetic relationships among regional populations of two species of grouper (Plectropomus maculatus and Plectropomus leopardus) were examined using mitochondrial and nuclear markers. mtDNA revealed contrasting regional inter-specific patterns whilst nuclear markers revealed contrasting patterns among markers, irrespective of region. In eastern Australia (EA) the species form a single mtDNA lineage, but the two species are reciprocally monophyletic in Western Australia (WA). This supports previous evidence for hybridisation between these species on the east coast. WA P. leopardus forms a sister relationship with the EA P. leopardus-maculatus clade while WA P. maculatus is more basal and sister to the P. leopardus lineages, indicating mtDNA does not suffer from incomplete lineage sorting for these species. In contrast, one of three nuclear markers (locus 7-90TG) differentiated the species into two reciprocally monophyletic clades, with no evidence of hybridisation or ancestral polymorphism. The remaining two nuclear markers (2-22 and ETS-2) did not separate these two species, while distinguishing other plectropomid species, suggesting incomplete lineage sorting at these nuclear loci. These results together with coalescence analyses suggest that P. leopardus females have hybridised historically with P. maculatus males and that P. maculatus mitochondria were displaced through introgressive hybridisation and fixation in the P. maculatus founder population on the Great Barrier Reef. The contrasting regional patterns of mtDNA structure may be attributed to Quaternary sea-level changes and shelf width differences driving different reef configurations on each coast. These reef configurations have provided opportunities for local scale interaction and reproduction among species on the narrower EA continental shelves, but not on the broader WA continental shelves.     

A multilocus phylogeny of the Sulidae (Aves: Pelecaniformes)

Gene trees will often differ from the true species history, the species tree, as a result of processes such as incomplete lineage sorting. New methods such as Bayesian Estimation of the Species Tree (BEST) use the multispecies coalescent to model lineage sorting, and directly infer the species tree from multilocus DNA sequence data. The Sulidae (Aves: Pelecaniformes) is a family of ten booby and gannet species with a global distribution. We sequenced five nuclear intron loci and one mitochondrial locus to estimate a species tree for the Sulidae using both BEST and by concatenating nuclear loci. We also used fossil calibrated strict and relaxed molecular clocks in BEAST to estimate divergence times for major nodes in the sulid phylogeny. Individual gene trees showed little phylogenetic conflict but varied in resolution. With the exception of the mitochondrial gene tree, no gene tree was completely resolved. On the other hand, both the BEST and concatenated species trees were highly resolved, strongly supported, and topologically consistent with each other. The three sulid genera (Morus, Sula, Papasula) were monophyletic and the relationships within genera were mostly consistent with both a previously estimated mtDNA gene tree and the mtDNA gene tree estimated here. However, our species trees conflicted with the mtDNA gene trees in the relationships among the three genera. Most notably, we find that the endemic and endangered Abbott's booby (Papasula abbotti) is likely basal to all other members of the Sulidae and diverged from them approximately 22 million years ago.