History, ocean channels, and distance determine phylogeographic patterns in three widespread Philippine fruit bats (Pteropodidae)

Phylogenomics and Spatiotemporal Dynamics of Bovine Leukemia Virus Focusing on Asian Native Cattle: Insights Into the Early Origin and Global Dissemination

Bovine leukemia virus (BLV), the causative agent of enzootic bovine leukosis, is currently one of the most important pathogens affecting the cattle industry worldwide. Determining where and in which host it originated, and how it dispersed across continents will provide valuable insights into its historical emergence as the cattle pathogen. Various species in the Bos genus were domesticated in Asia, where they also diversified. As native cattle (taurine cattle, zebu cattle, yak, and water buffalo) are indigenous and adapted to local environments, we hypothesized that Asian native cattle could have harbored BLV and, therefore, that they were important for virus emergence, maintenance, and spread. In this study, phylogeographic and ancestral trait analyses—including sequences obtained from Asian native cattle—were used to reconstruct the evolutionary history of BLV. It was shown that, since its probable emergence in Asia, the virus spread to South America and Europe via international trade of live cattle. It was inferred that zebu cattle were the hosts for the early origin of BLV, while taurine cattle played the significant role in the transmission worldwide. In addition, the results of positive selection analysis indicate that yak had a substantially minor role in the transmission of this virus. In this study, endogenous deltaretrovirus sequences in bats, collected in Asian countries, were also analyzed on whether these sequences were present in the bat genome. Endogenous deltaretrovirus sequences were detected from bat species endemic to specific regions and geographically isolated for a long time. Endogenous deltaretrovirus sequences from these geographically isolated species represent ancient exogenous deltaretroviruses distributions. The phylogenetic analysis revealed that these newly obtained endogenous deltaretrovirus sequences were closely related to those of BLV from Asian native cattle, indicating that BLV-related ancient deltaretroviruses circulated in Asia long before the emergence of BLV. Together, our analyses provide evidence for origin and spatiotemporal dynamics of BLV.

Colonization and diversification of the white-browed shortwing (Aves: Muscicapidae: Brachypteryx montana) in the Philippines

Molecular phylogenetic approaches have greatly improved our knowledge of the pattern and process of biological diversification across the globe; however, many regions remain poorly documented, even for well-studied vertebrate taxa. The Philippine archipelago, one of the least-studied 'biodiversity hotspots', is an ideal natural laboratory for investigating the factors driving diversification in an insular and geologically dynamic setting. We investigated the history and geography of diversification of the Philippine populations of a widespread montane bird, the White-browed Shortwing (Brachypteryx montana). Leveraging dense archipelago-wide sampling, we generated a multi-locus genetic dataset (one nuclear and two mtDNA markers), which we analyzed using phylogenetic, population genetic, and coalescent-based methods. Our results demonstrate that Philippine shortwings (1) likely colonized the Philippines from the Sunda Shelf to Mindanao in the late Miocene or Pliocene, (2) diversified across inter-island barriers into three divergent lineages during the Pliocene and early Pleistocene, (3) have not diversified within the largest island, Luzon, contrary to patterns observed in other montane taxa, and (4) colonized Palawan from the oceanic Philippines rather than from Borneo, challenging the assumption of Palawan functioning exclusively as a biogeographic extension of the Sunda Shelf. Additionally, our finding that divergent (c. 4.0 mya) lineages are coexisting in secondary sympatry on Mindanao without apparent gene flow suggests that the speciation process is likely complete for these shortwing lineages. Overall, these investigations provide insight into how topography and island boundaries influence diversification within remote oceanic archipelagos and echo the results of many other studies in demonstrating that taxonomic diversity continues to be underestimated in the Philippines.

Comparative phylogeography in the Atlantic forest and Brazilian savannas: pleistocene fluctuations and dispersal shape spatial patterns in two bumblebees

Background

Bombus morio and B. pauloensis are sympatric widespread bumblebee species that occupy two major Brazilian biomes, the Atlantic forest and the savannas of the Cerrado. Differences in dispersion capacity, which is greater in B. morio, likely influence their phylogeographic patterns. This study asks which processes best explain the patterns of genetic variation observed in B. morio and B. pauloensis, shedding light on the phenomena that shaped the range of local populations and the spatial distribution of intra-specific lineages.

Results

Results suggest that Pleistocene climatic oscillations directly influenced the population structure of both species. Correlative species distribution models predict that the warmer conditions of the Last Interglacial contributed to population contraction, while demographic expansion happened during the Last Glacial Maximum. These results are consistent with physiological data suggesting that bumblebees are well adapted to colder conditions. Intra-specific mitochondrial genealogies are not congruent between the two species, which may be explained by their documented differences in dispersal ability.

Conclusions

While populations of the high-dispersal B. morio are morphologically and genetically homogeneous across the species range, B. pauloensis encompasses multiple (three) mitochondrial lineages, and show clear genetic, geographic, and morphological differences. Because the lineages of B. pauloensis are currently exposed to distinct climatic conditions (and elevations), parapatric diversification may occur within this taxon. The eastern portion of the state of São Paulo, the most urbanized area in Brazil, represents the center of genetic diversity for B. pauloensis.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0803-0) contains supplementary material, which is available to authorized users.

Genetic Differentiation in Insular Lowland Rainforests: Insights from Historical Demographic Patterns in Philippine Birds

Phylogeographic studies of Philippine birds support that deep genetic structure occurs across continuous lowland forests within islands, despite the lack of obvious contemporary isolation mechanisms. To examine the pattern and tempo of diversification within Philippine island forests, and test if common mechanisms are responsible for observed differentiation, we focused on three co-distributed lowland bird taxa endemic to Greater Luzon and Greater Negros-Panay: Blue-headed Fantail (Rhipidura cyaniceps), White-browed Shama (Copsychus luzoniensis), and Lemon-throated Leaf-Warbler (Phylloscopus cebuensis). Each species has two described subspecies within Greater Luzon, and a single described subspecies on Greater Negros/Panay. Each of the three focal species showed a common geographic pattern of two monophyletic groups in Greater Luzon sister to a third monophyletic group found in Greater Negros-Panay, suggesting that common or similar biogeographic processes may have produced similar distributions. However, studied species displayed variable levels of mitochondrial DNA differentiation between clades, and genetic differentiation within Luzon was not necessarily concordant with described subspecies boundaries. Population genetic parameters for the three species suggested both rapid population growth from small numbers and geographic expansion across Luzon Island. Estimates of the timing of population expansion further supported that these events occurred asynchronously throughout the Pleistocene in the focal species, demanding particular explanations for differentiation, and support that co-distribution may be secondarily congruent.

Archipelago colonization by ecologically dissimilar amphibians: evaluating the expectation of common evolutionary history of geographical diffusion in co-distributed rainforest tree frogs in islands of Southeast Asia

Widespread, co-distributed species with limited relative dispersal abilities represent compelling focal taxa for comparative phylogeography. Forest vertebrates in island archipelagos often exhibit pronounced population structure resulting from limited dispersal abilities or capacity to overcome marine barriers to dispersal. The exceptionally diverse Old World tree frogs of the family Rhacophoridae have colonized the forested island archipelagos of Southeast Asia on multiple occasions, entering the islands of Indonesia and the Philippines via a "stepping stone" mode of dispersal along elongate island chains, separated by a series of marine channels. Here we evaluate the prediction that two tightly co-distributed Philippine rhacophorids colonized the archipelago during concomitant timescales and in the same, linear, "island-hopping" progression. We use a new multilocus dataset, utilize dense genetic sampling from the eastern arc of the Philippines, and we take a model-based phylogeographic approach to examining the two species for similar topological patterns of diversification, genetic structure, and timescales of diversification. Our results support some common mechanistic predictions (a general south-to-north polarity of colonization) but not others (timescale for colonization and manner and degree of lineage diversification), suggesting differing biogeographic scenarios of geographical diffusion through the archipelago and unique and idiosyncratic ecological capacities and evolutionary histories of each species.

Evidence for climate-driven diversification? A caution for interpreting ABC inferences of simultaneous historical events

Approximate Bayesian computation (ABC) is rapidly gaining popularity in population genetics. One example, msBayes, infers the distribution of divergence times among pairs of taxa, allowing phylogeographers to test hypotheses about historical causes of diversification in co-distributed groups of organisms. Using msBayes, we infer the distribution of divergence times among 22 pairs of populations of vertebrates distributed across the Philippine Archipelago. Our objective was to test whether sea-level oscillations during the Pleistocene caused diversification across the islands. To guide interpretation of our results, we perform a suite of simulation-based power analyses. Our empirical results strongly support a recent simultaneous divergence event for all 22 taxon pairs, consistent with the prediction of the Pleistocene-driven diversification hypothesis. However, our empirical estimates are sensitive to changes in prior distributions, and our simulations reveal low power of the method to detect random variation in divergence times and bias toward supporting clustered divergences. Our results demonstrate that analyses exploring power and prior sensitivity should accompany ABC model selection inferences. The problems we identify are potentially mitigable with uniform priors over divergence models (rather than classes of models) and more flexible prior distributions on demographic and divergence-time parameters.

Phylogeography of the prehensile-tailed skink Corucia zebrata on the Solomon Archipelago

The biogeography of islands is often strongly influenced by prior geological events. Corucia zebrata (Squamata: Scincidae) is endemic to the geologically complex Solomon Archipelago in Northern Melanesia. We examined the level of divergence for different island populations of C. zebrata and discussed these patterns in light of Pleistocene land bridges, island isolation, and island age. Corucia zebrata was sampled from 14 locations across the Solomon Archipelago and sequenced at two mitochondrial genes (ND2 and ND4; 1697 bp in total) and four nuclear loci (rhodopsin, an unknown intron, AKAP9, and PTPN12). Measures of genetic divergence, analyses of genetic variation, and Bayesian phylogenetic inference were used and the data assessed in light of geological information. Populations of C. zebrata on separate islands were found to be genetically different from each other, with reciprocal monophyly on mitochondrial DNA. Populations on islands previously connected by Pleistocene land bridges were marginally less divergent from each other than from populations on other nearby but isolated islands. There are indications that C. zebrata has radiated across the eastern islands of the archipelago within the last 1–4 million years. Nuclear loci were not sufficiently informative to yield further information about the phylogeography of C. zebrata on the Solomon Archipelago. Analyses of the mitochondrial data suggest that dispersal between islands has been very limited and that there are barriers to gene flow within the major islands. Islands that have been isolated during the Pleistocene glacial cycles are somewhat divergent in their mitochondrial genotypes, however, isolation by distance (IBD) and recent colonization of isolated but geologically younger islands appear to have had stronger effects on the phylogeography of C. zebrata than the Pleistocene glacial cycles. This contrasts with patterns reported for avian taxa, and highlights the fact that biogeographic regions for island species cannot be directly extrapolated among taxa of differing dispersal ability.

Molecular systematics of the Philippine forest skinks (Squamata: Scincidae: Sphenomorphus): testing morphological hypotheses of interspecific relationships

Skinks of the genus Sphenomorphus are the most diverse clade of squamates in the Philippine Archipelago. Morphological examination of these species has defined six phenotypic groups that are commonly used in characterizations of taxonomic hypotheses. We used a molecular phylogeny based on four mitochondrial and two nuclear genes to assess the group's biogeographical history in the archipelago and examine the phylogenetic validity of the currently recognized Philippine species groups. We re‐examined traditional characters used to define species groups and used multivariate statistics to quantitatively evaluate group structure in morphometric space. Clustering analyses of phenotypic similarity indicate that some (but not all) members of previously defined species groups are phenotypically most similar to other members of the same group. However, when species group membership was mapped on our partitioned Bayesian phylogenetic hypothesis, only one species group corresponds to a clade; all other species group arrangements are strongly rejected by our phylogeny. Our results demonstrate that (1) previously recognized species group relationships were misled by phenotypic convergence; (2) Sphenomorphus is widely paraphyletic; and (3) multiple lineages have independently invaded the Philippines. Based on this new perspective on the phylogenetic relationships of Philippine Sphenomorphus, we revise the archipelago's diverse assemblage of species at the generic level, and resurrect and/or expand four previously recognized genera, and describe two new genera to accommodate the diversity of Philippine skinks of the Sphenomorphus group.

© 2011 The Linnean Society of London, Zoological Journal of the Linnean Society, 2011, 163, 1217–1243.

Species interactions during diversification and community assembly in an island radiation of shrews

BACKGROUND

Closely related, ecologically similar species often have adjacent distributions, suggesting competitive exclusion may contribute to the structure of some natural communities. In systems such as island archipelagos, where speciation is often tightly associated with dispersal over oceanic barriers, competitive exclusion may prevent population establishment following inter-island dispersal and subsequent cladogenesis.

METHODOLOGY/PRINCIPAL FINDINGS

Using a combination of tools, we test the hypothesis that the distributions of shrew (Crocidura) species in the Philippines are the result of competitive exclusion preventing secondary invasion of occupied islands. We first compare ecological niche models between two widespread, allopatric species and find statistical support for their ecological similarity, implying that competition for habitat between these species is possible. We then examine dispersion patterns among sympatric species and find some signal for overdispersion of body size, but not for phylogenetic branch length. Finally, we simulate the process of inter-island colonization under a stochastic model of dispersal lacking ecological forces. Results are dependent on the geographic scope and colonization probability employed. However, some combinations suggest that the number of inter-island dispersal events necessary to populate the archipelago may be much higher than the minimum number of colonization events necessary to explain current estimates of species richness and phylogenetic relationships. If our model is appropriate, these results imply that alternative factors, such as competitive exclusion, may have influenced the process of inter-island colonization and subsequent cladogenesis.

CONCLUSIONS/SIGNIFICANCE

We interpret the combined results as providing tenuous evidence that similarity in body size may prevent co-occurrence in Philippine shrews and that competitive exclusion among ecologically similar species, rather than an inability to disperse among islands, may have limited diversification in this group, and, possibly other clades endemic to island archipelagos.

Exploring demographic, physical, and historical explanations for the genetic structure of two lineages of Greater Antillean bats

Observed patterns of genetic structure result from the interactions of demographic, physical, and historical influences on gene flow. The particular strength of various factors in governing gene flow, however, may differ between species in biologically relevant ways. We investigated the role of demographic factors (population size and sex-biased dispersal) and physical features (geographic distance, island size and climatological winds) on patterns of genetic structure and gene flow for two lineages of Greater Antillean bats. We used microsatellite genetic data to estimate demographic characteristics, infer population genetic structure, and estimate gene flow among island populations of Erophylla sezekorni/E. bombifrons and Macrotus waterhousii (Chiroptera: Phyllostomidae). Using a landscape genetics approach, we asked if geographic distance, island size, or climatological winds mediate historical gene flow in this system. Samples from 13 islands spanning Erophylla's range clustered into five genetically distinct populations. Samples of M. waterhousii from eight islands represented eight genetically distinct populations. While we found evidence that a majority of historical gene flow between genetic populations was asymmetric for both lineages, we were not able to entirely rule out incomplete lineage sorting in generating this pattern. We found no evidence of contemporary gene flow except between two genetic populations of Erophylla. Both lineages exhibited significant isolation by geographic distance. Patterns of genetic structure and gene flow, however, were not explained by differences in relative effective population sizes, island area, sex-biased dispersal (tested only for Erophylla), or surface-level climatological winds. Gene flow among islands appears to be highly restricted, particularly for M. waterhousii, and we suggest that this species deserves increased taxonomic attention and conservation concern.

Increased population sampling confirms low genetic divergence among Pteropus (Chiroptera: Pteropodidae) fruit bats of Madagascar and other western Indian Ocean islands

Fruit bats of the genus Pteropus occur throughout the Austral-Asian region west to islands off the eastern coast of Africa. Recent phylogenetic analyses of Pteropus from the western Indian Ocean found low sequence divergence and poor phylogenetic resolution among several morphologically defined species. We reexamine the phylogenetic relationships of these taxa by using multiple individuals per species. In addition, we estimate population genetic structure in two well-sampled taxa occurring on Madagascar and the Comoro Islands (P. rufus and P. seychellensis comorensis). Despite finding a similar pattern of low sequence divergence among species, increased sampling provides insight into the phylogeographic history of western Indian Ocean Pteropus, uncovering high levels of gene flow within species.

The phylogenetic relationships of cynopterine fruit bats (Chiroptera: Pteropodidae: Cynopterinae)

The subfamily Cynopterinae comprises ca. 24 species of pteropodid bats (Family Pteropodidae) distributed exclusively in South and Southeast Asia. Although some studies have supported monophyly of the subfamily, molecular analyses have produced contradictory results and there has been little agreement on relationships of cynopterines to other megabat groups. However, no previous studies have included a complete sampling of cynopterine genera. Here we describe a phylogenetic analysis of Cynopterinae based on more than 6000 bp from six different genes sampled in representatives of all 14 recognized genera. Our results support the monophyly of Cynopterinae but refute a close relationship of cynopterines with Nyctimeninae. Within Cynopterinae, our analyses consistently recovered two monophyletic clades, which we recommend be recognized formally as tribes: Cynopterini and Balionycterini. Biogeographic analyses indicate a Sundaland origin of the Cynopterinae and divergence date estimates suggest different timing of diversification of the two major cynopterine clades.

The role of tropical dry forest as a long-term barrier to dispersal: a comparative phylogeographical analysis of dry forest tolerant and intolerant frogs

We used a comparative phylogeographical approach to investigate the origins of the disjunct wet forest biota of the Golfo Dulce region along the Pacific slope of Costa Rica. This region is isolated by Pacific dry forests north and south and isolated from Caribbean wet forests by mountains. We studied three sympatric lowland frog species in the Craugastor fitzingeri species group that prefer wet forest but differ in their response to dry habitats. In dry forest, C. fitzingeri can survive along streams while C. crassidigitus and C. talamancae are entirely absent. We collected samples from across the ranges of all three species, and obtained mitochondrial DNA sequence data from the COI and cytochrome b genes. We observed significant phylogeographical structure in C. crassidigitus and C. talamancae, but much less in C. fitzingeri, demonstrating that mountain barriers and dry forest habitat have reduced mitochondrial gene flow in the strictly wet-forest species. Additionally, we discovered that the Golfo Dulce and Central Panama populations of C. crassidigitus appear to have diverged in the Pliocene or earlier, suggesting that the dry forest separating these populations is old. Our phylogenetic analysis of 12 of approximately 16 species of the C. fitzingeri species group suggests that the three lowland species are each other's closest relatives. Because of this shared phylogenetic history, we attribute the striking differences in phylogeographical structure to the different ecologies of the frogs. In summary, we find that what appear to be minor differences in the natural history of these three closely related species may profoundly impact the potential for dispersal, range size, and cladogenesis.