A multilocus phylogeny of the Sulidae (Aves: Pelecaniformes)

Karyotype evolution of suliformes and description of a ♂Z1Z1Z2Z2/♀Z1Z2W multiple sex chromosome system in boobies (Sula spp.)

Our comprehension of avian karyotypes still needs to be improved, especially for Suliform birds. To enhance understanding of chromosomal evolution in this order, we conducted conventional and molecular cytogenetic analysis in five species, named Sula dactylatra, Sula leucogaster, Sula sula (Sulidae), Fregata magnificens (Fregatidae), and Nannopterum brasilianum (Phalacrocoracidae). The diploid chromosome number for S. dactylatra and S. leucogaster was established as 2n = 76 in males, and 2n = 75 in females, but S. sula displayed a karyotype of 2n = 76 chromosomes in males. The disparity in diploid chromosome numbers between male and female Sula is due to a multiple sex chromosome system of the Z1Z1Z2Z2/Z1Z2W type. We propose that the emergence of this multiple-sex chromosome system resulted from a Robertsonian translocation involving the W chromosome and the smallest microchromosome. Fregata magnificens exhibited a diploid number 76 (2n = 76), while N. brasilianum displayed a diploid number of 74 (2n = 74) in both sexes. The ribosomal cluster was located in one microchromosome pair in S. dactylatra, S. leucogaster, S. sula, and F. magnificens and in four pairs in N. brasilianum. Our findings provide evidence of a conserved multiple-sex chromosome system within the Sula genus, shedding light on the high karyotype diversity in Suliformes.

A new fossil subspecies of booby (Aves, Sulidae: Papasula) from Mauritius and Rodrigues, Mascarene Islands, with notes on P. abbotti from Assumption Island

A new subspecies of Papasula booby is described from fossil remains collected in the Mascarene Islands of Mauritius and Rodrigues, southwestern Indian Ocean. The Mascarene Booby Papasula abbotti nelsoni ssp. nov., larger than nominate Abbott's Booby P. a. abbotti from Christmas Island, northeast Indian Ocean, was approximately the same size as the extinct Hiva Oa Booby P. a. costelloi from the Marquesas in the Eastern Pacific. Mentioned in early accounts of Mauritius in 1668 and in 1725-26 and 1761 on Rodrigues, the Mascarene booby became extinct by the end of the 18th century. Members of another isolated but now extirpated population of Papasula abbotti from Assumption Island in the Seychelles Archipelago, the island from which the type specimen was collected, are identical to extant Christmas Island birds in size and colouration and discussed in detail herein.

Complete Mitochondrial Genome of Great Frigatebird (Fregata minor): Phylogenetic Position and Gene Rearrangement

The complete mitogenome sequence of the Great Frigatebird, Fregata minor was sequenced for the first time in this study. The mitogenome (16,899 bp) comprises of 13 protein-coding genes (PCGs), two ribosomal RNA (rRNA) genes, and 22 transfer RNA (tRNA) genes, and a control region (CR). The mitogenome was AT-rich (55.60%) with 11 overlapping and 18 intergenic spacer regions. Most of the PCGs were started by a typical ATG initiation codon except for cox1 and nad3. A maximum-likelihood phylogeny of concatenated PCGs resulted in a well-resolved phylogeny of all the species of Suliformes and illuminates the sister relationship of F. minor with F. magnificens. The present mitogenome-based phylogeny clearly enlightens the evolutionary position of Suliformes and Pelecaniformes species. Unique tandem repeats were identified in both F. minor and F. magnificens, which can be employed as a species-specific marker. To illuminate the population structure of this migratory seabirds, the present study advocate more sampling and the generation of additional molecular data to clarify their genetic diversity. The present study also rejects an earlier hypothesis on the mitochondrial gene order of Suliformes and corroborated the typical avian gene order in frigatebirds.

The Galápagos Islands: biogeographic patterns and geology

In the traditional biogeographic model, the Galápagos Islands appeared a few million years ago in a sea where no other islands existed and were colonized from areas outside the region. However, recent work has shown that the Galápagos hotspot is 139 million years old (Early Cretaceous), and so groups are likely to have survived at the hotspot by dispersal of populations onto new islands from older ones. This process of metapopulation dynamics means that species can persist indefinitely in an oceanic region, as long as new islands are being produced. Metapopulations can also undergo vicariance into two metapopulations, for example at active island arcs that are rifted by transform faults. We reviewed the geographic relationships of Galápagos groups and found 10 biogeographic patterns that are shared by at least two groups. Each of the patterns coincides spatially with a major tectonic structure; these structures include: the East Pacific Rise; west Pacific and American subduction zones; large igneous plateaus in the Pacific; Alisitos terrane (Baja California), Guerrero terrane (western Mexico); rifting of North and South America; formation of the Caribbean Plateau by the Galápagos hotspot, and its eastward movement; accretion of Galápagos hotspot tracks; Andean uplift; and displacement on the Romeral fault system. All these geological features were active in the Cretaceous, suggesting that geological change at that time caused vicariance in widespread ancestors. The present distributions are explicable if ancestors survived as metapopulations occupying both the Galápagos hotspot and other regions before differentiating, more or less in situ.

Redrawing Wallace’s Line based on the fauna of Christmas Island, eastern Indian Ocean

Based on a comprehensive literature survey, we determined the sources of the terrestrial vertebrate species on Christmas Island, asking where they originated relative to Wallace’s Line (the southern end of the divide lies 1100 km to the east, where the Lombok Strait adjoins the eastern Indian Ocean). The two bats, Pipistrellus murrayi and Pteropus natalis, are from the west. Concerning the endemic and ‘resident’ bird species, one is from the west (Collocalia natalis), four are from the east (Accipiter fasciatus, Egretta novaehollandiae, Falco cenchroides and Ninox natalis) and the other 15 are ambiguous or indeterminate. Most of the land-locked species are also from the east: rodents Rattus macleari and Rattus nativitatis, and squamates Cryptoblepharus egeriae, Emoia nativitatis and Lepidodactylus listeria. Additionally, two have westerly origins (Crocidura trichura and Cyrtodactylus sadleiri), one is ambiguous (Emoia atrocostata) and another is unknown (Ramphotyphlops exocoeti). West-directed surface currents that flow across the eastern Indian Ocean towards Christmas Island would have facilitated most of the land-animal colonizations. We therefore suggest that Wallace’s Line be redrawn such that the landmass is placed on the Australasian side of this fundamental biogeographical boundary.

Additive Traits Lead to Feeding Advantage and Reproductive Isolation, Promoting Homoploid Hybrid Speciation

Speciation through homoploid hybridization (HHS) is considered extremely rare in animals. This is mainly because the establishment of reproductive isolation as a product of hybridization is uncommon. Additionally, many traits are underpinned by polygeny and/or incomplete dominance, where the hybrid phenotype is an additive blend of parental characteristics. Phenotypically intermediate hybrids are usually at a fitness disadvantage compared with parental species and tend to vanish through backcrossing with parental population(s). It is therefore unknown whether the additive nature of hybrid traits in itself could lead successfully to HHS. Using a multi-marker genetic data set and a meta-analysis of diet and morphology, we investigated a potential case of HHS in the prions (Pachyptila spp.), seabirds distinguished by their bills, prey choice, and timing of breeding. Using approximate Bayesian computation, we show that the medium-billed Salvin's prion (Pachyptila salvini) could be a hybrid between the narrow-billed Antarctic prion (Pachyptila desolata) and broad-billed prion (Pachyptila vittata). Remarkably, P. salvini's intermediate bill width has given it a feeding advantage with respect to the other Pachyptila species, allowing it to consume a broader range of prey, potentially increasing its fitness. Available metadata showed that P. salvini is also intermediate in breeding phenology and, with no overlap in breeding times, it is effectively reproductively isolated from either parental species through allochrony. These results provide evidence for a case of HHS in nature, and show for the first time that additivity of divergent parental traits alone can lead directly to increased hybrid fitness and reproductive isolation.

Fish diversification at the pace of geomorphological changes: evolutionary history of western Iberian Leuciscinae (Teleostei: Leuciscidae) inferred from multilocus sequence data

The evolutionary history of western Iberian Leuciscinae, obligatory freshwater fish, is directly linked to the evolution of the hydrographic network of the Iberian Peninsula after its isolation from the rest of Europe, which involved dramatic rearrangements such as the transition from endorheic lakes to open basins draining to the Atlantic. Previous phylogenetic research on western Iberian leuciscines, using mainly mitochondrial DNA and more recently one or two nuclear genes, has found contradictory results and there remain many unresolved issues regarding species relationships, taxonomy, and evolutionary history. Moreover, there is a lack of integration between phylogenetic and divergence time estimates and information on the timing of geomorphological changes and paleobasin rearrangements in the Iberian Peninsula. This study presents the first comprehensive fossil-calibrated multilocus coalescent species tree of western Iberian Leuciscinae (including 14 species of Achondrostoma, Iberochondrostoma, Pseudochondrostoma and Squalius endemic to the Iberian Peninsula, seven of which endemic to Portugal) based on seven nuclear genes, and from which we infer their biogeographic history by comparing divergence time estimates to known dated geological events. The phylogenetic pattern suggests slow-paced evolution of leuciscines during the Early-Middle Miocene endorheic phase of the main Iberian river basins, with the shift to exorheism in the late Neogene-Quaternary allowing westward dispersals that resulted in many cladogenetic events and a high rate of endemism in western Iberia. The results of this study also: (i) confirm the paraphyly of S. pyrenaicus with respect to S. carolitertii, and thus the possible presence of a new taxon in the Portuguese Tagus currently assigned to S. pyrenaicus; (ii) support the taxonomic separation of the Guadiana and Sado populations of S. pyrenaicus; (iii) show the need for further population sampling and taxonomic research to clarify the phylogenetic status of A. arcasii from the Minho basin and of the I. lusitanicum populations in the Sado and Tagus basins; and (iv) indicate that A. occidentale, I. olisiponensis and P. duriensis are the most ancient lineages within their respective genera.

A phylogenetic test of sympatric speciation in the Hydrobatinae (Aves: Procellariiformes)

Phylogenetic relationships among species can provide insight into how new species arise. For example, careful consideration of both the phylogenetic and geographic distributions of species in a group can reveal the geographic models of speciation within the group. One such model, sympatric speciation, may be more common than previously thought. The Hydrobatinae (Aves: Procellariformes) is a diverse subfamily of Northern Hemisphere storm-petrels for which the taxonomy is unclear. Previous studies showed that Hydrobates (formally Oceanodroma) castro breeding in the Azores during the cool season is sister species to H. monteiroi, a hot season breeder at the same locations, which suggests sympatric speciation by allochrony. To test whether other species within the subfamily arose via sympatric speciation by allochrony, we sequenced the cytochrome b gene and five nuclear introns to estimate a phylogenetic tree using multispecies coalescent methods, and to test whether species breeding in the same geographic area are monophyletic. We found that speciation within the Hydrobatinae appears to have followed several geographic modes of divergence. Sympatric seasonal species in Japan likely did not arise through sympatric speciation, but allochrony may have played a role in the divergence of H. matsudairae, a cool season breeder, and H. monorhis, a hot season breeder. No other potential cases of sympatric speciation were discovered within the subfamily. Despite breeding in the same geographic area, hydrobatine storm-petrels breeding in Baja California (H. microsoma and H. melania) are each sister to a species breeding off the coast of Peru (H. tethys and H. markhami, respectively). In fact, antitropical sister species appear to have diverged at multiple times, suggesting allochronic divergence might be common. In addition, allopatry has likely played a role in divergence of H. furcata, a north Pacific breeder, and H. pelagius, a north Atlantic breeder. This study demonstrates that a variety of mechanisms of divergence have played a role in generating the diversity of the Hydrobatinae and supports the current taxonomy of the subfamily.

Sexual Size Dimorphism and Body Condition in the Australasian Gannet

Sexual size dimorphism is widespread throughout seabird taxa and several drivers leading to its evolution have been hypothesised. While the Australasian Gannet (Morus serrator) has previously been considered nominally monomorphic, recent studies have documented sexual segregation in diet and foraging areas, traits often associated with size dimorphism. The present study investigated the sex differences in body mass and structural size of this species at two colonies (Pope’s Eye, PE; Point Danger, PD) in northern Bass Strait, south-eastern Australia. Females were found to be 3.1% and 7.3% heavier (2.74 ± 0.03, n = 92; 2.67 ± 0.03 kg, n = 43) than males (2.66 ± 0.03, n = 92; 2.48 ± 0.03 kg, n = 43) at PE and PD, respectively. Females were also larger in wing ulna length (0.8% both colonies) but smaller in bill depth (PE: 2.2%; PD: 1.7%) than males. Despite this dimorphism, a discriminant function provided only mild accuracy in determining sex. A similar degree of dimorphism was also found within breeding pairs, however assortative mating was not apparent at either colony (R² < 0.04). Using hydrogen isotope dilution, a body condition index was developed from morphometrics to estimate total body fat (TBF) stores, where TBF(%) = 24.43+1.94*(body mass/wing ulna length) – 0.58*tarsus length (r² = 0.84, n = 15). This index was used to estimate body composition in all sampled individuals. There was no significant difference in TBF(%) between the sexes for any stage of breeding or in any year of the study at either colony suggesting that, despite a greater body mass, females were not in a better condition than males. While the driving mechanism for sexual dimorphism in this species is currently unknown, studies of other Sulids indicate segregation in foraging behaviour, habitat and diet may be a contributing factor.

Evolutionary factors affecting the cross-species utility of newly developed microsatellite markers in seabirds

Microsatellite loci are ideal for testing hypotheses relating to genetic segregation at fine spatio-temporal scales. They are also conserved among closely related species, making them potentially useful for clarifying interspecific relationships between recently diverged taxa. However, mutations at primer binding sites may lead to increased nonamplification, or disruptions that may result in decreased polymorphism in nontarget species. Furthermore, high mutation rates and constraints on allele size may also with evolutionary time, promote an increase in convergently evolved allele size classes, biasing measures of interspecific genetic differentiation. Here, we used next-generation sequencing to develop microsatellite markers from a shotgun genome sequence of the sub-Antarctic seabird, the thin-billed prion (Pachyptila belcheri), that we tested for cross-species amplification in other Pachyptila and related sub-Antarctic species. We found that heterozygosity decreased and the proportion of nonamplifying loci increased with phylogenetic distance from the target species. Surprisingly, we found that species trees estimated from interspecific FST provided better approximations of mtDNA relationships among the studied species than those estimated using DC , even though FST was more affected by null alleles. We observed a significantly nonlinear second order polynomial relationship between microsatellite and mtDNA distances. We propose that the loss of linearity with increasing mtDNA distance stems from an increasing proportion of homoplastic allele size classes that are identical in state, but not identical by descent. Therefore, despite high cross-species amplification success and high polymorphism among the closely related Pachyptila species, we caution against the use of microsatellites in phylogenetic inference among distantly related taxa.

Phylogeny and forelimb disparity in waterbirds

Previous work has shown that the relative proportions of wing components (i.e., humerus, ulna, carpometacarpus) in birds are related to function and ecology, but these have rarely been investigated in a phylogenetic context. Waterbirds including "Pelecaniformes," Ciconiiformes, Procellariiformes, Sphenisciformes, and Gaviiformes form a highly supported clade and developed a great diversity of wing forms and foraging ecologies. In this study, forelimb disparity in the waterbird clade was assessed in a phylogenetic context. Phylogenetic signal was assessed via Pagel's lambda, Blomberg's K, and permutation tests. We find that different waterbird clades are clearly separated based on forelimb component proportions, which are significantly correlated with phylogeny but not with flight style. Most of the traditional contents of "Pelecaniformes" (e.g., pelicans, cormorants, and boobies) cluster with Ciconiiformes (herons and storks) and occupy a reduced morphospace. These taxa are closely related phylogenetically but exhibit a wide range of ecologies and flight styles. Procellariiformes (e.g., petrels, albatross, and shearwaters) occupy a wide range of morphospace, characterized primarily by variation in the relative length of carpometacarpus and ulna. Gaviiformes (loons) surprisingly occupy a wing morphospace closest to diving petrels and penguins. Whether this result may reflect wing proportions plesiomorphic for the waterbird clade or a functional signal is unclear. A Bayesian approach detecting significant rate shifts across phylogeny recovered two such shifts. At the base of the two sister clades Sphenisciformes + Procellariiformes, a shift to an increase evolutionary rate of change is inferred for the ulna and carpometacarpus. Thus, changes in wing shape begin prior to the loss of flight in the wing-propelled diving clade. Several shifts to slower rate of change are recovered within stem penguins.

Phylogeography and genetic structure of two Patagonian shag species (Aves: Phalacrocoracidae)

We compared the phylogeographic and genetic structure of two sympatric shag species, Phalacrocorax magellanicus (rock shag) and Phalacrocorax atriceps (imperial shag), from Patagonia (southern South America). We used multilocus genotypes of nuclear DNA (microsatellite loci) from 324 individuals and mitochondrial DNA sequences (ATPase) from 177 individuals, to evaluate hypotheses related to the effect of physical and non-physical barriers on seabird evolution. Despite sharing many ecological traits, the focal species strongly differ in two key aspects: P. magellanicus has a strong tendency to remain at/near their breeding colonies during foraging trips and the non-breeding season, while P. atriceps exhibits the converse pattern. Both species showed similar mtDNA genetic structure, where colonies from the Atlantic Coast, Pacific Coast and Fuegian region were genetically divergent. We also found similarities in the results of Bayesian clustering analysis of microsatellites, with both species having four clusters. However population differentiation (e.g. Fst, Φst) was higher in P. magellanicus compared to P. atriceps, and average membership probabilities of individuals to specific clusters (Q-values) were also higher in the former. Phalacrocorax magellanicus has strong phylogeographic structure, consistent with the impact of Pleistocene glaciations, with diagnostic haplotypes associated with each of the three mentioned regions. The same pattern was not as evident for P. atriceps. Migration rate estimators were higher for P. atriceps than for P. magellanicus; however both species followed an n-island-like model of gene flow, this implies that dispersal occurs across the continental land mass that separates Atlantic and Pacific Oceans. Our results supported the hypothesis that non-physical barriers are important drivers of the genetic and phylogeographic structure in seabirds, and also that physical barriers constitute effective but not absolute impediments to gene flow.

Evidence for asymmetrical divergence-gene flow of nuclear loci, but not mitochondrial loci, between seabird sister species: blue-footed (Sula nebouxii) and Peruvian (S. variegata) boobies

Understanding the process of speciation requires understanding how gene flow influences divergence. Recent analyses indicate that divergence can take place despite gene flow and that the sex chromosomes can exhibit different levels of gene flow than autosomes and mitochondrial DNA. Using an eight marker dataset including autosomal, z-linked, and mitochondrial loci we tested the hypothesis that blue-footed (Sula nebouxii) and Peruvian (S. variegata) boobies diverged from their common ancestor with gene flow, paying specific attention to the differences in gene flow estimates from nuclear and mitochondrial markers. We found no gene flow at mitochondrial markers, but found evidence from the combined autosomal and z-linked dataset that blue-footed and Peruvian boobies experienced asymmetrical gene flow during or after their initial divergence, predominantly from Peruvian boobies into blue-footed boobies. This gene exchange may have occurred either sporadically between periods of allopatry, or regularly throughout the divergence process. Our results add to growing evidence that diverging species can remain distinct but exchange genes.

Beyond phylogeny: pelecaniform and ciconiiform birds, and long-term niche stability

Phylogenetic trees are a starting point for the study of further evolutionary and ecological questions. We show that for avian evolutionary relationships, improved taxon sampling, longer sequences and additional data sets are giving stability to the prediction of the grouping of pelecaniforms and ciconiiforms, thus allowing inferences to be made about long-term niche occupancy. Here we report the phylogeny of the pelecaniform birds and their water-carnivore allies using complete mitochondrial genomes, and show that the basic groupings agree with nuclear sequence phylogenies, even though many short branches are not yet fully resolved. In detail, we show that the Pelecaniformes (minus the tropicbird) and the Ciconiiformes (storks, herons and ibises) form a natural group within a seabird water-carnivore clade. We find pelicans are the closest relatives of the shoebill (in a clade with the hammerkop), and we confirm that tropicbirds are not pelecaniforms. In general, the group appears to be an adaptive radiation into an 'aquatic carnivore' niche that it has occupied for 60-70 million years. From an ecological and life history perspective, the combined pelecaniform-ciconiform group is more informative than focusing on differences in morphology. These findings allow a start to integrating molecular evolution and macroecology.

An improved phylogeny of the Andean tit-tyrants (Aves, Tyrannidae): more characters trump sophisticated analyses

The phylogeny of the flycatcher genus Anairetes was previously inferred using short fragments of mitochondrial DNA and parsimony and distance-based methods. The resulting topology spurred taxonomic revision and influenced understanding of Andean biogeography. More than a decade later, we revisit the phylogeny of Anairetes tit-tyrants using more mtDNA characters, seven unlinked loci (three mitochondrial genes, six nuclear loci), more closely related outgroup taxa, partitioned Bayesian analyses, and two coalescent species-tree approaches (Bayesian estimation of species trees, BEST; Bayesian evolutionary analysis by sampling trees, (*)BEAST). Of these improvements in data and analyses, the fourfold increase in mtDNA characters was both necessary and sufficient to incur a major shift in the topology and near-complete resolution. The species-tree analyses, while theoretically preferable to concatenation or single gene approaches, yielded topologies that were compatible with mtDNA but with weaker statistical resolution at nodes. The previous results that had led to taxonomic and biogeographic reappraisal were refuted, and the current results support the resurrection of the genus Uromyias as the sister clade to Anairetes. The sister relationship between these two genera corresponds to an ecological dichotomy between a depauperate humid cloud forest clade and a diverse dry-tolerant clade that has diversified along the latitudinal axis of the Andes. The species-tree results and the concatenation results each reaffirm the primacy of mtDNA to provide phylogenetic signal for avian phylogenies at the species and subspecies level. This is due in part to the abundance of informative characters in mtDNA, and in part to its lower effective population size that causes it to more faithfully track the species tree.

Eight independent nuclear genes support monophyly of the plovers: the role of mutational variance in gene trees

Molecular phylogenies of Charadriiformes based on mtDNA genes and one to three nuclear loci do not support the traditional placement of Pluvialis in the plovers (Charadriidae), assigning it instead to oystercatchers, stilts, and avocets (Haematopodidae and Recurvirostridae). To investigate this hypothesis of plover paraphyly, the relationships among Pluvialis and closely related families were revisited by sequencing two individuals of all taxa except Peltohyas for eight independent single copy nuclear protein-coding loci selected for their informativeness at this phylogenetic depth. The species tree estimated jointly with the gene trees in the coalescent programme (*)BEAST strongly supported plover monophyly, as did Bayesian analysis of the concatenated matrix. The data sets that supported plover paraphyly in Baker et al. (2007) and Fain and Houde (2007) reflect two to four independent gene histories, and thus discordance with the plover monophyly species tree might have arisen by chance through stochastic mutational variance. For the plovers we conclude there is no conclusive evidence of coalescent variance from ancient incomplete lineage sorting across the interior branch leading to Pluvialis in the species tree. Rather, earlier studies seem have been misled by faster evolving mtDNA genes with high mutational variance, and a few nuclear genes that had low resolving power at the Pluvialis sister group level. These findings are of general relevance in avian phylogenetics, as they show that careful attention needs to be paid to the number and the phylogenetic informativeness of genes required to obtain accurate estimates of the species tree, especially where there is mutational heterogeneity in gene trees.