Climate refugia: joint inference from fossil records, species distribution models and phylogeography

Climate refugia, locations where taxa survive periods of regionally adverse climate, are thought to be critical for maintaining biodiversity through the glacial–interglacial climate changes of the Quaternary. A critical research need is to better integrate and reconcile the three major lines of evidence used to infer the existence of past refugia – fossil records, species distribution models and phylogeographic surveys – in order to characterize the complex spatiotemporal trajectories of species and populations in and out of refugia. Here we review the complementary strengths, limitations and new advances for these three approaches. We provide case studies to illustrate their combined application, and point the way towards new opportunities for synthesizing these disparate lines of evidence. Case studies with European beech, Qinghai spruce and Douglas-fir illustrate how the combination of these three approaches successfully resolves complex species histories not attainable from any one approach. Promising new statistical techniques can capitalize on the strengths of each method and provide a robust quantitative reconstruction of species history. Studying past refugia can help identify contemporary refugia and clarify their conservation significance, in particular by elucidating the fine-scale processes and the particular geographic locations that buffer species against rapidly changing climate.

Honey Bee Apis mellifera Parasites in the Absence of Nosema ceranae Fungi and Varroa destructor Mites

Few areas of the world have western honey bee (Apis mellifera) colonies that are free of invasive parasites Nosema ceranae (fungi) and Varroa destructor (mites). Particularly detrimental is V. destructor; in addition to feeding on host haemolymph, these mites are important vectors of several viruses that are further implicated as contributors to honey bee mortality around the world. Thus, the biogeography and attendant consequences of viral communities in the absence of V. destructor are of significant interest. The island of Newfoundland, Province of Newfoundland and Labrador, Canada, is free of V. destructor; the absence of N. ceranae has not been confirmed. Of 55 Newfoundland colonies inspected visually for their strength and six signs of disease, only K-wing had prevalence above 5% (40/55 colonies = 72.7%). Similar to an earlier study, screenings again confirmed the absence of V. destructor, small hive beetles Aethina tumida (Murray), tracheal mites Acarapis woodi (Rennie), and Tropilaelaps spp. ectoparasitic mites. Of a subset of 23 colonies screened molecularly for viruses, none had Israeli acute paralysis virus, Kashmir bee virus, or sacbrood virus. Sixteen of 23 colonies (70.0%) were positive for black queen cell virus, and 21 (91.3%) had some evidence for deformed wing virus. No N. ceranae was detected in molecular screens of 55 colonies, although it is possible extremely low intensity infections exist; the more familiar N. apis was found in 53 colonies (96.4%). Under these conditions, K-wing was associated (positively) with colony strength; however, viruses and N. apis were not. Furthermore, black queen cell virus was positively and negatively associated with K-wing and deformed wing virus, respectively. Newfoundland honey bee colonies are thus free of several invasive parasites that plague operations in other parts of the world, and they provide a unique research arena to study independent pathology of the parasites that are present.

Life history, ecology and the biogeography of strong genetic breaks among 15 species of PaciWc rockWsh, Sebastes

Strong genetic change over short spatial scales is surprising among marine species with high dispersal potential. Concordant breaks among several species signals a role for geographic barriers to dispersal. Along the coast of California, such breaks have not been seen across the biogeographic barrier of Point Conception, but other potential geographic boundaries have been surveyed less often.;We tested for strong-population structure in 11 species of Sebastes sampled across two regions containing potential dispersal barriers, and conducted a meta-analysis including four additional species. We show two strong breaks north of Monterey Bay, spanning an oceanographic gradient and an upwelling jet. Moderate genetic structure is just as common in the north as it is in the south, across the biogeographic break at Point Conception. Gene Xow is generally higher among deep-water species, but these conclusions are confounded by phylogeny. Species in the subgenus Sebastosomus have higher structure than those in the subgenus;Pteropodus, despite having larvae with longer pelagic phases. DiVerences in settlement behavior in the face of ocean currents might help explain these diVerences. Across similar species across the same coastal environment, we document a wide variety of patterns in gene Xow, suggesting that interaction of individual species traits such as settlement behavior with environmental factors such as;oceanography can strongly impact population structure

Phylogeography of silver pheasant (Lophura nycthemera L.) across China: aggregate effects of refugia, introgression, and riverine barriers

The role of Pleistocene glacial cycles in forming the contemporary genetic structure of organisms has been well studied in China with a particular focus on the Tibetan Plateau. However, China has a complex topography and diversity of local climates, and how glacial cycles may have shaped the subtropical and tropical biota of the region remains mostly unaddressed. To investigate the factors that affected the phylogeography and population history of a widely distributed and nondeciduous forest species, we analysed morphological characters, mitochondrial DNA sequences and nuclear microsatellite loci in the Silver Pheasant (Lophura nycthemera). In a pattern generally consistent with phenotypic clusters, but not nominal subspecies, deeply divergent mitochondrial lineages restricted to different geographic regions were detected. Coalescent simulations indicated that the time of main divergence events corresponded to major glacial periods in the Pleistocene and gene flow was only partially lowered by drainage barriers between some populations. Intraspecific cytonuclear discordance was revealed in mitochondrial lineages from Hainan Island and the Sichuan Basin with evidence of nuclear gene flow from neighbouring populations into the latter. Unexpectedly, hybridization was revealed in Yingjiang between the Silver Pheasant and Kalij Pheasant (Lophura leucomelanos) with wide genetic introgression at both the mtDNA and nuclear levels. Our results highlight a novel phylogeographic pattern in a subtropical area generated from the combined effects of climate oscillation, partial drainage barriers and interspecific hybridization. Cytonuclear discordance combined with morphological differentiation implies that complex historical factors shaped the divergence process in this biodiversity hot spot area of southern China.

Evidence for selection on coloration in a Panamanian poison frog: a coalescent-based approach

Aim The strawberry poison frog, Oophaga pumilio, has undergone a remarkable radiation of colour morphs in the Bocas del Toro archipelago in Panama. This species shows extreme variation in colour and pattern between populations that have been geographically isolated for < 10,000 years. While previous research has suggested the involvement of divergent selection, to date no quantitative test has examined this hypothesis. Location Bocas del Toro archipelago, Panama. Methods We use a combination of population genetics, phylogeography and phenotypic analyses to test for divergent selection in coloration in O. pumilio. Tissue samples of 88 individuals from 15 distinct populations were collected. Using these data, we developed a gene tree using the mitochondrial DNA (mtDNA) d-loop region. Using parameters derived from our mtDNA phylogeny, we predicted the coalescence of a hypothetical nuclear gene underlying coloration. We collected spectral reflectance and body size measurements on 94 individuals from four of the populations and performed a quantitative analysis of phenotypic divergence. Results The mtDNA d-loop tree revealed considerable polyphyly across populations. Coalescent reconstructions of gene trees within population trees revealed incomplete genotypic sorting among populations. The quantitative analysis of phenotypic divergence revealed complete lineage sorting by colour, but not by body size: populations showed non-overlapping variation in spectral reflectance measures of body coloration, while variation in body size did not separate populations. Simulations of the coalescent using parameter values derived from our empirical analyses demonstrated that the level of sorting among populations seen in colour cannot reasonably be attributed to drift. Main conclusions These results imply that divergence in colour, but not body size, is occurring at a faster rate than expected under neutral processes. Our study provides the first quantitative support for the claim that strong diversifying selection underlies colour variation in the strawberry poison frog.

Macroevolutionary patterns in the diversification of parrots: effects of climate change, geological events and key innovations

Aim Parrots are thought to have originated on Gondwana during the Cretaceous. The initial split within crown group parrots separated the New Zealand taxa from the remaining extant species and was considered to coincide with the separation of New Zealand from Gondwana 82-85 Ma, assuming that the diversification of parrots was mainly shaped by vicariance. However, the distribution patterns of several extant parrot groups cannot be explained without invoking transoceanic dispersal, challenging this assumption. Here, we present a temporal and spatial framework for the diversification of parrots using external avian fossils as calibration points in order to evaluate the relative importance of the influences of past climate change, plate tectonics and ecological opportunity. Location Australasian, African, Indo-Malayan and Neotropical regions. Methods Phylogenetic relationships were investigated using partial sequences of the nuclear genes c-mos, RAG-1 and Zenk of 75 parrot and 21 other avian taxa. Divergence dates and confidence intervals were estimated using a Bayesian relaxed molecular clock approach. Biogeographic patterns were evaluated taking temporal connectivity between areas into account. We tested whether diversification remained constant over time and if some parrot groups were more species-rich than expected given their age. Results Crown group diversification of parrots started only about 58 Ma, in the Palaeogene, significantly later than previously thought. The Australasian lories and possibly also the Neotropical Arini were found to be unexpectedly species-rich. Diversification rates probably increased around the Eocene/Oligocene boundary and in the middle Miocene, during two periods of major global climatic aberrations characterized by global cooling. Main conclusions The diversification of parrots was shaped by climatic and geological events as well as by key innovations. Initial vicariance events caused by continental break-up were followed by transoceanic dispersal and local radiations. Habitat shifts caused by climate change and mountain orogenesis may have acted as a catalyst to the diversification by providing new ecological opportunities and challenges as well as by causing isolation as a result of habitat fragmentation. The lories constitute the only highly nectarivorous parrot clade, and their diet shift, associated with morphological innovation, may have acted as an evolutionary key innovation, allowing them to explore underutilized niches and promoting their diversification.

Diversity versus disparity and the role of ecological opportunity in a continental bird radiation

Aim The Neotropical parrots (Arini) are an unusually diverse group which colonized South America in the Oligocene. The newly invaded Neotropics may have functioned as an underused adaptive zone and provided novel ecological opportunities that facilitated diversification. Alternatively, diversification may have been driven by ecological changes caused by Andean uplift and/or climate change from the Miocene onwards. Our aim was to find out whether Arini diversified in a classical adaptive radiation after their colonization of South America, or whether their diversification occurred later and was influenced by more recent environmental change. Location Neotropics. Methods We generated a time-calibrated phylogeny of more than 80% of all Arini species in order to analyse lineage diversification. This chronogram was also used as the basis for the reconstruction of morphological evolution within Arini using a multivariate ratio analysis of three size measurements. Results We found a concentration of size evolution and partitioning of size niches in the early history of Arini consistent with the process of adaptive radia- tion, but there were no signs of an early burst of speciation or a decrease in speci- ation rates through time. Although we detected no overall temporal shifts in diversification rates, we discovered two young, unexpectedly species-rich clades. Main conclusions Arini show signs of an early adaptive radiation, but we found no evidence of the slowdown in speciation rate generally considered a feature of island or lake radiations. Historical processes and environmental change from the Miocene onwards may have kept diversification rates roughly constant ever since the colonization of the Neotropics. Thus, Arini may not yet have reached equilibrium diversity. The lack of diversity-dependent speciation might be a general feature of adaptive radiations on a continental scale, and diversification processes on continents might therefore not be as ecologically limited as in isolated lakes or on oceanic islands.

East Asian ethnolinguistic phylogeography

A polyphyletic understanding of Asian linguistic diversity was first propagated in 1823. Since 1901, various scholars have proposed larger linguistic phyla uniting two or more recognised Asian language families. The most recent proposal in this tradition, Starosta’s 2001 East Asian phylum, comprising the Trans-Himalayan, Hmong-Mien, Austroasiatic, Austronesian and Kradai language families, is reassessed in light of linguistic and non-linguistic evidence. Ethnolinguistically informed inferences based on Asian Y chromosomal phylogeography lead to a reconstruction of various episodes of ethnolinguistic prehistory which lie beyond the linguistic event horizon, i.e. at a time depth empirically inaccessible to historical linguistics. The Father Tongue correlation in population genetics, the evidence for refugia during the Last Glacial Maximum and the hypothesis of language families having arisen as the result of demographic bottlenecks in prehistory are shown to be crucial to an understanding of the ethnogenesis of East Asian linguistic phyla. The prehistory of several neighbouring Asian language families is discussed, and the Centripetal Migration model is opposed to the Farming Language Dispersal theory.

A prehistoric thoroughfare between the Ganges and the Himalayas

The greater Himalayan region, including the Tibetan plateau in the north and the Gangetic plain in the south, served as the principal prehistoric thoroughfare for the peopling of East and Southeast Asia. The descendants of ancient migrants through this region ultimately settled lands as far away as New Zealand, Madagascar and the Americas. Several of the keys to understanding the ethnogenesis of human diversity in Asia include the Father Tongue correlation, possible refugia during the Last Glacial Maximum and the hypothesis that language families may have arisen as the result of demographic bottlenecks in prehistory. Ethnolinguistically informed inferences based on Asian Y chromosomal phylogeography permit a reconstruction of episodes of ethnolinguistic prehistory which lie beyond the linguistic event horizon, i.e. beyond the time depth empirically accessible to historical linguistics. The origins of the language families which make up the hypothetical Uralo-Siberian and East Asian linguistic phyla are argued to have lain in the northeastern corner of the Indian subcontinent. Several other Asian language families are shown to be tied to the subcontinent. The Centripetal Migration model, which assumes that migrations in quest of a better life unfolded in both centrifugal and centripetal directions with respect to technologically more advanced centres of civilisation, is opposed to the Farming Language Dispersal theory, which assumes that all linguistic dispersals were driven by agricultural centrifugal migration.