Notes on the statistical power of the binary state speciation and extinction (BiSSE) model

The Binary State Speciation and Extinction (BiSSE) method is one of the most popular tools for investigating the rates of diversification and character evolution. Yet, based on previous simulation studies, it is commonly held that the BiSSE method requires phylogenetic trees of fairly large sample sizes (>300 taxa) in order to distinguish between the different models of speciation, extinction, or transition rate asymmetry. Here, the power of the BiSSE method is reevaluated by simulating trees of both small and large sample sizes (30, 60, 90, and 300 taxa) under various asymmetry models and root state assumptions. Results show that the power of the BiSSE method can be much higher, also in trees of small sample size, for detecting differences in speciation rate asymmetry than anticipated earlier. This, however, is not a consequence of any conceptual or mathematical flaw in the method per se but rather of assumptions about the character state at the root of the simulated trees and thus the underlying macroevolutionary model, which led to biased results and conclusions in earlier power assessments. As such, these earlier simulation studies used to determine the power of BiSSE were not incorrect but biased, leading to an overestimation of type-II statistical error for detecting differences in speciation rate but not for extinction and transition rates.P17124-B0(VLID)241678

Multiple independent origins of auto-pollination in tropical orchids (Bulbophyllum) in light of the hypothesis of selfing as an evolutionary dead end

Background: The transition from outcrossing to selfing has long been portrayed as an ‘evolutionary dead end because, first, reversals are unlikely and, second, selfing lineages suffer from higher rates of extinction owing to a reduced potential for adaptation and the accumulation of deleterious mutations. We tested these two predictions in a clade of Madagascan Bulbophyllum orchids (30 spp.), including eight species where auto-pollinating morphs (i.e., selfers, without a ‘rostellum) co-exist with their pollinator-dependent conspecifics (i.e., outcrossers, possessing a rostellum). Specifically, we addressed this issue on the basis of a time-calibrated phylogeny by means of ancestral character reconstructions and within the state-dependent evolution framework of BiSSE (Binary State Speciation and Extinction), which allowed jointly estimating rates of transition, speciation, and extinction between outcrossing and selfing. Results: The eight species capable of selfing occurred in scattered positions across the phylogeny, with two likely originating in the Pliocene (ca. 4.43.1 Ma), one in the Early Pleistocene (ca. 2.4 Ma), and five since the mid-Pleistocene (ca. 1.3 Ma). We infer that this scattered phylogenetic distribution of selfing is best described by models including up to eight independent outcrossing-to-selfing transitions and very low rates of speciation (and either moderate or zero rates of extinction) associated with selfing. Conclusions: The frequent and irreversible outcrossing-to-selfing transitions in Madagascan Bulbophyllum are clearly congruent with the first prediction of the dead end hypothesis. The inability of our study to conclusively reject or support the likewise predicted higher extinction rate in selfing lineages might be explained by a combination of methodological limitations (low statistical power of our BiSSE approach to reliably estimate extinction in small-sized trees) and evolutionary processes (insufficient time elapsed for selfers to go extinct). We suggest that, in these tropical orchids, a simple genetic basis of selfing (via loss of the ‘rostellum) is needed to explain the strikingly recurrent transitions to selfing, perhaps reflecting rapid response to parallel and novel selective environments over Late Quaternary ( 1.3 Ma) time scales.P20726-B03P17124-B0(VLID)243455

Performance comparison of two reduced-representation based genome-wide marker-discovery strategies in a multi-taxon phylogeographic framework

Multi-locus genetic data are pivotal in phylogenetics. Today, high-throughput sequencing (HTS) allows scientists to generate an unprecedented amount of such data from any organism. However, HTS is resource intense and may not be accessible to wide parts of the scientific community. In phylogeography, the use of HTS has concentrated on a few taxonomic groups, and the amount of data used to resolve a phylogeographic pattern often seems arbitrary. We explore the performance of two genetic marker sampling strategies and the effect of marker quantity in a comparative phylogeographic framework focusing on six species (arthropods and plants). The same analyses were applied to data inferred from amplified fragment length polymorphism fingerprinting (AFLP), a cheap, non-HTS based technique that is able to straightforwardly produce several hundred markers, and from restriction site associated DNA sequencing (RADseq), a more expensive, HTS-based technique that produces thousands of single nucleotide polymorphisms. We show that in four of six study species, AFLP leads to results comparable with those of RADseq. While we do not aim to contest the advantages of HTS techniques, we also show that AFLP is a robust technique to delimit evolutionary entities in both plants and animals. The demonstrated similarity of results from the two techniques also strengthens biological conclusions that were based on AFLP data in the past, an important finding given the wide utilization of AFLP over the last decades. We emphasize that whenever the delimitation of evolutionary entities is the central goal, as it is in many fields of biodiversity research, AFLP is still an adequate technique.Te present study was co-funded by the Austrian Science Fund (FWF, project P25955 “Origin of steppe fora and fauna in inner-Alpine dry valleys” to P.S.), and the Tiroler Wissenschafsfonds (TWF, UNI-0404/2066,“Comparing information efciency of high- versus low-resolution genome scans for phylogeographic studies” to P.K.). Te computational results presented have been achieved using the HPC infrastructure LEO of the University of Innsbruck

Histological and micro-CT evidence of stigmatic rostellum receptivity promoting auto-pollination in the Madagascan orchid Bulbophyllum bicoloratum

Background: The rostellum, a projecting part of the gynostemium in orchid flowers, separates the anther(s) from the stigma and thus commonly prevents auto-pollination. Nonetheless, as a modified (usually distal) portion of the median stigma lobe, the rostellum has been frequently invoked of having re-gained a stigmatic function in rare cases of orchid auto-pollination. Here it is shown that a newly discovered selfing variant of Madagascan Bulbophyllum bicoloratum has evolved a modified rostellum allowing the penetration of pollen tubes from in situ pollinia. Methods: Gynostemium micro-morphology and anatomy of selfing and outcrossing variants of B. bicoloratum was studied by using light and scanning electron microscopy and histological sections. Pollen tube growth in the selfing variant was further observed via X-ray computed microtomography (micro-CT), providing 3D reconstructions of floral tissues at a micron scale. Findings: Selfing variants possess a suberect (‘displaced) rostellum rather than the conventional, erect type. Very early in anthesis, the pollinia of selfers are released from the anther and slide down onto the suberect rostellum, where pollen tube growth preferentially occurs through the non-vascularized, i.e. rear (adaxial) and (semi-) lateral parts. This penetrated tissue is comprised of a thin layer of elongate and loosely arranged cells, embedded in stigmatic exudates, as also observed in the stigmatic cavity of both selfing and outcrossing variants. Conclusions: Our results provide the first solid evidence of a stigmatic function for the rostellum in orchid flowers, thereby demonstrating for the first time the feasibility of the micro-CT technique for accurately visualizing pollen tube growth in flowering plants. Rostellum receptivity in B. bicoloratum probably uniquely evolved as an adaptation for reproductive assurance from an outcrossing ancestor possessing an erect (non-receptive) rostellum. These findings open up new avenues in the investigation of an organ that apparently re-gained its ‘primordial function of being penetrated by pollen tubes.P20726-B03P17124-B0(VLID)170467

Integrating phylogenomics, phylogenetics, morphometrics, relative genome size and ecological niche modelling disentangles the diversification of Eurasian Euphorbia seguieriana s. l. (Euphorbiaceae)

Next generation sequencing has revolutionised biology. Restriction-associated DNA sequencing (RADseq) has primarily been used to study infraspecific relationships but has also been applied in multi-species phylogenomic analyses. In this study, we used a combination of phylogenomic (with RADseq data) and phylogenetic (with sequences of the nuclear internal transcribed spacer, ITS) methods to explore relationships within the taxonomically intricate Euphorbia seguieriana s. l., one of the most widespread Euphorbia taxa inhabiting zonal and extrazonal steppes from Iberia to Central Asia. In the inferred phylogenies the southeastern Balkan and Anatolian populations were clearly separated, supporting the distinction of E. niciciana from E. seguieriana at the species level. Within E. seguieriana, the populations from the Caucasus, Iran, and easternmost Anatolia were sister to all other populations based on RADseq, making necessary the description of a new, morphologically divergent subspecies, E. seguieriana subsp. armeniaca. Conversely, additional studies are needed to understand the status of E. seguieriana subsp. hohenackeri, which is sympatric with E. seguieriana subsp. armeniaca. Niche analyses indicated that differences in the climatic niche between E. niciciana and E. seguieriana are relatively small compared with the climatic differences between the regions over which they are distributed. Contrary to previous believes, E. niciciana and E. seguieriana are allopatric and have likely diverged during the Pleistocene in two different glacial refugia as suggested by distribution modelling. Euphorbia niciciana nowadays has a submediterranean distribution, occupying habitats that are slightly warmer, moister, and less seasonal in temperature but more seasonal in precipitation than E. seguieriana, a characteristic species of continental steppes. Using flow cytometry, we demonstrate that the relative genome sizes of E. niciciana and E. seguieriana differ significantly. Additionally, multivariate morphometric analyses of 56 morphological characters indicated clear morphological divergence of the two species. Importantly, we also provide a revised taxonomic treatment including formal nomenclatural changes, an identification key and species descriptions. Our study demonstrates that an integrative approach, combining modern phylogenomic methods with traditional phylogenetic, cytogenetic, environmental and morphological analyses can result in satisfactorily resolved relationships in intricate groups of closely related species. Finally, phylogenetic inference using ITS sequences is still a useful tool for resolving relationships among the taxa at the species level, but the phylogenomic approach based on RADseq data certainly provides better resolution both among and within species.Austrian Science Fund Tiroler Wissenschaftsfond

Genomic Signature of Shifts in Selection in a Subalpine Ant and Its Physiological Adaptations

Understanding how organisms adapt to extreme environments is fundamental and can provide insightful case studies for both evolutionary biology and climate-change biology. Here, we take advantage of the vast diversity of lifestyles in ants to identify genomic signatures of adaptation to extreme habitats such as high altitude. We hypothesized two parallel patterns would occur in a genome adapting to an extreme habitat: 1) strong positive selection on genes related to adaptation and 2) a relaxation of previous purifying selection. We tested this hypothesis by sequencing the high-elevation specialist Tetramorium alpestre and four other phylogenetically related species. In support of our hypothesis, we recorded a strong shift of selective forces in T. alpestre, in particular a stronger magnitude of diversifying and relaxed selection when compared with all other ants. We further disentangled candidate molecular adaptations in both gene expression and protein-coding sequence that were identified by our genome-wide analyses. In particular, we demonstrate that T. alpestre has 1) a higher level of expression for stv and other heat-shock proteins in chill-shock tests and 2) enzymatic enhancement of Hex-T1, a rate-limiting regulatory enzyme that controls the entry of glucose into the glycolytic pathway. Together, our analyses highlight the adaptive molecular changes that support colonization of high-altitude environments.Research was supported by the Austrian Science Fund (FWF, P23409 and P30861)

Genomic Signature of Shifts in Selection in a Subalpine Ant and Its Physiological Adaptations

Understanding how organisms adapt to extreme environments is fundamental and can provide insightful case studies for both evolutionary biology and climate-change biology. Here, we take advantage of the vast diversity of lifestyles in ants to identify genomic signatures of adaptation to extreme habitats such as high altitude. We hypothesized two parallel patterns would occur in a genome adapting to an extreme habitat: 1) strong positive selection on genes related to adaptation and 2) a relaxation of previous purifying selection. We tested this hypothesis by sequencing the high-elevation specialist Tetramorium alpestre and four other phylogenetically related species. In support of our hypothesis, we recorded a strong shift of selective forces in T. alpestre, in particular a stronger magnitude of diversifying and relaxed selection when compared with all other ants. We further disentangled candidate molecular adaptations in both gene expression and protein-coding sequence that were identified by our genome-wide analyses. In particular, we demonstrate that T. alpestre has 1) a higher level of expression for stv and other heat-shock proteins in chill-shock tests and 2) enzymatic enhancement of Hex-T1, a rate-limiting regulatory enzyme that controls the entry of glucose into the glycolytic pathway. Together, our analyses highlight the adaptive molecular changes that support colonization of high-altitude environments

Bio18_Pleistocene

Oscillayers (Bio18) for time periods 20 kyr to 2.570 My

Bio17_Pleistocene

Oscillayers (Bio17) for time periods 20 kyr to 2.570 My

Bio5_Pliocene

Oscillayers (Bio5) for time periods 2.580 Myr to 5.40 My