Untangling the colonization history of the Australo-Pacific reed warblers, one of the world's great island radiations.

Island radiations, such as those of the Australo-Pacific, offer unique insight into diversification, extinction, and early speciation processes. Yet, their speciation and colonization histories are often obscured by conflicting genomic signals from incomplete lineage sorting or hybridization. Here, we integrated mitogenomes and genome-wide SNPs to unravel the evolutionary history of one of the world's most geographically widespread island radiations. The Australo-Pacific reed warblers (Acrocephalus luscinius complex) are a speciose lineage including five species that have become extinct since the 19th century and ten additional species of conservation concern. The radiation spans over 10,000 km across Australo-Papua, Micronesia and Polynesia, including the Mariana, Hawaii and Pitcairn Island archipelagos. Earlier mtDNA studies suggested a stepping-stone colonization process, resulting in archipelago-level secondary sympatry of divergent mtDNA lineages in the Mariana Islands and Marquesas. These studies hypothesised that morphologically similar species on neighbouring islands arose from ecological convergence. Using hDNA from historical museum specimens and modern genetic samples, we show that incomplete lineage sorting and/or gene flow have shaped the radiation of Australo-Pacific reed warblers rather than secondary sympatry. The nuclear genome reconstructs a simpler biogeographic history than mtDNA, showing close relationships between species in the Mariana Islands and Marquesas despite their paraphyletic mtDNA lineages. Gene flow likely involved early and late colonizing waves of the radiation before the loss of ancestral dispersive ability. Our results highlight how collection genomics can elucidate evolutionary history and inform conservation efforts for threatened species.

Runs of homozygosity reveal past bottlenecks and contemporary inbreeding across diverging populations of an island-colonizing bird.

Genomes retain evidence of the demographic history and evolutionary forces that have shaped populations and drive speciation. Across island systems, contemporary patterns of genetic diversity reflect population demography, including colonization events, bottlenecks, gene flow and genetic drift. Here, we investigate genome-wide diversity and the distribution of runs of homozygosity (ROH) using whole-genome resequencing of individuals (>22× coverage) from six populations across three archipelagos of Berthelot's pipit (Anthus berthelotii)-a passerine that has recently undergone island speciation. We show the most dramatic reduction in diversity occurs between the mainland sister species (the tawny pipit) and Berthelot's pipit and is lowest in the populations that have experienced sequential bottlenecks (i.e., the Madeiran and Selvagens populations). Pairwise sequential Markovian coalescent (PSMC) analyses estimated that Berthelot's pipit diverged from its sister species ~2 million years ago, with the Madeiran archipelago founded 50,000 years ago, and the Selvagens colonized 8000 years ago. We identify many long ROH (>1 Mb) in these most recently colonized populations. Population expansion within the last 100 years may have eroded long ROH in the Madeiran archipelago, resulting in a prevalence of short ROH (<1 Mb). However, the extensive long and short ROH detected in the Selvagens suggest strong recent inbreeding and bottleneck effects, with as much as 38% of the autosomes consisting of ROH >250 kb. These findings highlight the importance of demographic history, as well as selection and genetic drift, in shaping contemporary patterns of genomic diversity across diverging populations.

Conservation genomics and systematics of a near-extinct island radiation.

Conservation benefits from incorporating genomics to explore the impacts of population declines, inbreeding, loss of genetic variation and hybridization. Here we use the near-extinct Mariana Islands reedwarbler radiation to showcase how ancient DNA approaches can allow insights into the population dynamics of extinct species and threatened populations for which historical museum specimens or material with low DNA yield (e.g., scats, feathers) are the only sources for DNA. Despite their having paraphyletic mitochondrial DNA (mtDNA), nuclear single nucleotide polymorphisms (SNPs) support the distinctiveness of critically endangered Acrocephalus hiwae and the other three species in the radiation that went extinct between the 1960s and 1990s. Two extinct species, A. yamashinae and A. luscinius, were deeply divergent from each other and from a third less differentiated lineage containing A. hiwae and extinct A. nijoi. Both mtDNA and SNPs suggest that the two isolated populations of A. hiwae from Saipan and Alamagan Islands are sufficiently distinct to warrant subspecies recognition and separate conservation management. We detected no significant differences in genetic diversity or inbreeding between Saipan and Alamagan, nor strong signatures of geographical structuring within either island. However, the implications of possible signatures of inbreeding in both Saipan and Alamagan, and long-term population declines in A. hiwae that pre-date modern anthropogenic threats require further study with denser population sampling. Our study highlights the value that conservation genomics studies of island radiations have as windows onto the possible future for the world's biota as climate change and habitat destruction increasingly fragment their ranges and contribute to rapid declines in population abundances.

Ancestral transoceanic colonization and recent population reduction in a nonannual killifish from the Seychelles archipelago.

Whether freshwater fish colonize remote islands following tectonic or transoceanic dispersal remains an evolutionary puzzle. Integrating dating of known tectonic events with phylogenomics and current species distribution, we find that killifish species distribution is not explained by species dispersal by tectonic drift only. Investigating the colonization of a nonannual killifish (golden panchax, Pachypanchax playfairii) on the Seychelle islands, we found genetic support for transoceanic dispersal and experimentally discovered an adaptation to complete tolerance to seawater. At the macroevolutionary scale, despite their long-lasting isolation, nonannual golden panchax show stronger genome-wide purifying selection than annual killifishes from continental Africa. However, progressive decline in effective population size over a more recent timescale has probably led to the segregation of slightly deleterious mutations across golden panchax populations, which represents a potential threat for species preservation in the long term.

Speciation and secondary contact in a fossorial island endemic, the São Tomé caecilian.

A period of isolation in allopatry typically precedes local adaptation and subsequent divergence among lineages. Alternatively, locally adapted phenotypes may arise and persist in the face of gene flow, resulting in strong correlations between ecologically-relevant phenotypic variation and corresponding environmental gradients. Quantifying genetic, ecological, and phenotypic divergence in such lineages can provide insights into the abiotic and biotic mechanisms that structure populations and drive the accumulation of phenotypic and taxonomic diversity. Low-vagility organisms whose distributions span ephemeral geographic barriers present the ideal evolutionary context within which to address these questions. Here, we combine genetic (mtDNA and genome-wide SNPs) and phenotypic data to investigate the divergence history of caecilians (Amphibia: Gymnophiona) endemic to the oceanic island of São Tomé in the Gulf of Guinea archipelago. Consistent with a previous mtDNA study, we find two phenotypically and genetically distinct lineages that occur along a north-to-south axis with extensive admixture in the centre of the island. Demographic modelling supports divergence in allopatry (~300 kya) followed by secondary contact (~95 kya). Consequently, in contrast to a morphological study that interpreted latitudinal phenotypic variation in these caecilians as a cline within a single widespread species, our analyses suggest a history of allopatric lineage divergence and subsequent hybridization that may have blurred species boundaries. We propose that late Pleistocene volcanic activity favoured allopatric divergence between these lineages with local adaptation to climate maintaining a stable hybrid zone in the centre of São Tomé Island. Our study joins a growing number of systems demonstrating lineage divergence on volcanic islands with stark environmental transitions across small geographic distances.

Um período de isolamento em alopatria geralmente precede adaptação local e divergência subsequente entre linhagens evolutivas. Alternativamente, fenótipos adaptados localmente podem surgir e persistir apesar de fluxo gênico, resultando em fortes correlações entre variação fenotípica ecologicamente relevante e os gradientes ambientais correspondentes. Quantificar divergência genética, ecológica e fenotípica em tais linhagens pode ajudar a clarificar os mecanismos abióticos e bióticos que estruturam as populações e levam ao acúmulo de diversidade fenotípica e taxonômica. Organismos de baixa vagilidade, cujas áreas de distribuição incluem barreiras geográficas efêmeras, representam um contexto evolutivo ideal para abordar essas questões. Neste estudo, combinamos dados genéticos (mtDNA e SNPs genômicos) e fenotípicos para investigar a história de divergência de cecílias endêmicas da ilha oceânica de São Tomé, no arquipélago do Golfo da Guiné. Consistentemente com um estudo anterior de mtDNA, encontramos duas linhagens fenotipicamente e geneticamente distintas que ocorrem ao longo de um eixo norte-sul, com extensa mistura genética no centro da ilha. Modelagem demográfica suportou um cenário de divergência em alopatria (~ 300 mil anos atrás) seguida de contato secundário (~ 95 mil anos atrás). Ao contrário de um estudo morfológico que interpretou a variação fenotípica latitudinal nessas cecílias como uma clina dentro de uma única espécie amplamente difundida, nossas análises sugerem uma história de divergência de linhagens em alopatria e subsequente hibridização que pode ter confundido os limites das espécies. Propomos que atividade vulcânica durante o Pleistoceno tardio favoreceu divergência alopátrica entre essas linhagens, com adaptação local ao clima mantendo uma zona híbrida estável no centro da Ilha de São Tomé. Nosso estudo se une a um número crescente de sistemas que demonstram divergência entre linhagens em ilhas vulcânicas com transições ambientais marcantes ao longo de distâncias geográficas curtas.

Drosophila yakuba mayottensis, a new model for the study of incipient ecological speciation.

A full understanding of how ecological factors drive the fixation of genetic changes during speciation is obscured by the lack of appropriate models with clear natural history and powerful genetic toolkits. In a recent study, we described an early stage of ecological speciation in a population of the generalist species Drosophila yakuba (melanogaster subgroup) on the island of Mayotte (Indian Ocean). On this island, flies are strongly associated with the toxic fruits of noni (Morinda citrifolia) and show a partial degree of pre-zygotic reproductive isolation. Here, I mine the nuclear and mitochondrial genomes and provide a full morphological description of this population. Only 29 nuclear sites (< 4 × 10-7 of the genome) are fixed in this population and absent from 3 mainland populations and the closest relative D. santomea, but no mitochondrial or morphological character distinguish Mayotte flies from the mainland. This result indicates that physiological and behavioral traits may evolve faster than morphology at the early stages of speciation. Based on these differences, the Mayotte population is designated as a new subspecies, Drosophila yakuba mayottensis subsp. nov., and its strong potential in understanding the genetics of speciation and plant-insect interactions is discussed.

Recurrent specialization on a toxic fruit in an island Drosophila population.

Recurrent specialization on similar host plants offers a unique opportunity to unravel the evolutionary and genetic mechanisms underlying dietary shifts. Recent studies have focused on ecological races belonging to the same species, but it is hard in many cases to untangle the role of adaptive introgression versus distinct mutations in facilitating recurrent evolution. We discovered on the island of Mayotte a population of the generalist fly Drosophila yakuba that is strictly associated with noni (Morinda citrifolia). This case strongly resembles Drosophila sechellia, a genetically isolated insular relative of D. yakuba whose intensely studied specialization on toxic noni fruits has always been considered a unique event in insect evolution. Experiments revealed that unlike mainland D. yakuba strains, Mayotte flies showed strong olfactory attraction and significant toxin tolerance to noni. Island females strongly discriminated against mainland males, suggesting that dietary adaptation has been accompanied by partial reproductive isolation. Population genomic analysis indicated a recent colonization (∼29 kya), at a time when year-round noni fruits may have presented a predictable resource on the small island, with ongoing migration after colonization. This relatively recent time scale allowed us to search for putatively adaptive loci based on genetic variation. Strong signals of genetic differentiation were found for several detoxification genes, including a major toxin tolerance locus in D. sechellia Our results suggest that recurrent evolution on a toxic resource can involve similar historical events and common genetic bases, and they establish an important genetic system for the study of early stages of ecological specialization and speciation.

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

BACKGROUND AND AIMS: Genetic drift due to geographical isolation, gene flow and mutation rates together make it difficult to determine the evolutionary relationships of present-day species. In this study, population genetic data were used to model and decipher interspecific relationships, speciation patterns and gene flow between three species of spruce with similar morphology, Picea wilsonii, P. neoveitchii and P. morrisonicola. Picea wilsonii and P. neoveitchii occur from central to north-west China, where they have overlapping distributions. Picea morrisonicola, however, is restricted solely to the island of Taiwan and is isolated from the other two species by a long distance. METHODS: Sequence variations were examined in 18 DNA fragments for 22 populations, including three fragments from the chloroplast (cp) genome, two from the mitochondrial (mt) genome and 13 from the nuclear genome. KEY RESULTS: In both the cpDNA and the mtDNA, P. morrisonicola accumulated more species-specific mutations than the other two species. However, most nuclear haplotypes of P. morrisonicola were shared by P. wilsonii, or derived from the dominant haplotypes found in that species. Modelling of population genetic data supported the hypothesis that P. morrisonicola derived from P. wilsonii within the more recent past, most probably indicating progenitor-derivative speciation with a distinct bottleneck, although further gene flow from the progenitor to the derivative continued. In addition, the occurrence was detected of an obvious mtDNA introgression from P. neoveitchii to P. wilsonii despite their early divergence. CONCLUSIONS: The extent of mutation, introgression and lineage sorting taking place during interspecific divergence and demographic changes in the three species had varied greatly between the three genomes. The findings highlight the complex evolutionary histories of these three Asian spruce species.