Assessing Bayesian Phylogenetic Information Content of Morphological Data Using Knowledge from Anatomy Ontologies

Morphology remains a primary source of phylogenetic information for many groups of organisms, and the only one for most fossil taxa. Organismal anatomy is not a collection of randomly assembled and independent "parts", but instead a set of dependent and hierarchically nested entities resulting from ontogeny and phylogeny. How do we make sense of these dependent and at times redundant characters? One promising approach is using ontologies---structured controlled vocabularies that summarize knowledge about different properties of anatomical entities, including developmental and structural dependencies. Here we assess whether the proximity of ontology-annotated characters within an ontology predicts evolutionary patterns. To do so, we measure phylogenetic information across characters and evaluate if it is hierarchically structured by ontological knowledge---in much the same way as phylogeny structures across-species diversity. We implement an approach to evaluate the Bayesian phylogenetic information (BPI) content and phylogenetic dissonance among ontology-annotated anatomical data subsets. We applied this to datasets representing two disparate animal groups: bees (Hexapoda: Hymenoptera: Apoidea, 209 chars) and characiform fishes (Actinopterygii: Ostariophysi: Characiformes, 463 chars). For bees, we find that BPI is not substantially structured by anatomy since dissonance is often high among morphologically related anatomical entities. For fishes, we find substantial information for two clusters of anatomical entities instantiating concepts from the jaws and branchial arch bones, but among-subset information decreases and dissonance increases substantially moving to higher-level subsets in the ontology. We further applied our approach to addressing particular evolutionary hypotheses with an example of morphological evolution in miniature fishes. While we show that ontology does indeed structure phylogenetic information, additional relationships and processes, such as convergence, likely play a substantial role in explaining BPI and dissonance, and merit future investigation. Our work demonstrates how complex morphological datasets can be interrogated with ontologies by allowing one to access how information is spread hierarchically across anatomical concepts, how congruent this information is, and what sorts of processes may structure it: phylogeny, development, or convergence.

randtip, a generalized framework to expand incomplete phylogenies using non-molecular phylogenetic information

The increasing availability of molecular information has lifted our understanding of species evolutionary relationships to unprecedent levels. However, current estimates of the world's biodiversity suggest that about a fifth of all extant species are yet to be described, and we still lack molecular information for many of the known species. Hence, evolutionary biologists will have to tackle phylogenetic uncertainty for a long time to come. This prospect has urged the development of software to expand phylogenies based on non-molecular phylogenetic information, and while the available tools provide some valuable features, major drawbacks persist and some of the proposed solutions are hardly generalizable to any group of organisms. Here, we present a completely generalized and flexible framework to expand incomplete molecular phylogenies. The framework is implemented in the R package "randtip", a toolkit of functions that was designed to randomly bind phylogenetically uncertain taxa in backbone phylogenies through a fully customizable and automatic procedure that uses taxonomic ranks as a major source of phylogenetic information. Although randtip is capable of automatically generating fully operative phylogenies for any group of organisms using just a list of species and a backbone tree, we stress that the "blind" expansion of phylogenies (using randtip or any other available software) often leads to suboptimal solutions. Thus, we discuss a variety of circumstances that may require customizing simulation parameters beyond default settings to optimally expand the trees, including a detailed step-by-step workflow. Phylogenetic uncertainty should be tackled with caution, assessing potential pitfalls and opportunities to optimize parameter space prior to launch any simulation. Used judiciously, our framework will help evolutionary biologists to efficiently expand incomplete molecular phylogenies and thereby account for phylogenetic uncertainty in quantitative analyses.

Macroevolutionary dynamics of climatic niche space

How and why lineages evolve along niche space as they diversify and adapt to different environments is fundamental to evolution. Progress has been hampered by the difficulties of linking a comprehensive empirical characterization of species niches with flexible evolutionary models that describe their evolution. Consequently, the relative influence of external episodic and biotic factors remains poorly understood. Here we characterize species' two-dimensional temperature and precipitation niche space occupied (i.e., species niche envelope) as complex geometries and assess their evolution across a large vertebrate radiation (all Aves) using a model that captures heterogeneous evolutionary rates on time-calibrated phylogenies. We find that extant birds coevolved from warm, mesic climatic niches into colder and drier environments and responded to the K-Pg boundary with a dramatic increase in disparity. Contrary to expectations of subsiding rates of niche evolution as lineages diversify, our results show that overall rates have increased steadily, with some lineages experiencing exceptionally high evolutionary rates, associated with colonization of novel niche spaces, and others showing niche stasis. Both competition- and environmental change-driven niche evolution transpire and result in highly heterogeneous rates near the present. Our findings share the limitations of all work based purely on extant taxa but highlight the growing ecological and conservation insights arising from the model-based integration of increasingly comprehensive and robust environmental and phylogenetic information.

Dynamics of inter-farm transmission of highly pathogenic avian influenza H5N6 integrating vehicle movements and phylogenetic information

Highly pathogenic avian influenza (HPAI) in poultry holdings commonly spreads through animal trade, and poultry production and health-associated vehicle (PPHaV) movement. To effectively control the spread of disease, it is essential that the contact structure via those movements among farms is thoroughly explored. However, few attempts have been made to scrutinize PPHaV movement compared to poultry trade. Therefore, our study aimed to elucidate the role of PPHaV movement on HPAI transmission. We performed network analysis using PPHaV movement data based on a global positioning system, with phylogenetic information of the isolates during the 2016–2017 HPAI H5N6 epidemic in the Republic of Korea. Moreover, the contribution of PPHaV movement to the spread of HPAI was estimated by Bayesian modeling. The network analysis revealed that there was the relationship between phylogenetic clusters and the contact network via PPHaV movement. Furthermore, the similarity of farm poultry species and the shared integrators between inter-linked infected premises (IPs) were associated with ties within the same phylogenetic clusters. Additionally, PPHaV movement among phylogenetically clustered IPs was estimated to contribute to approximately 30% of HPAI H5N6 infections in IPs on average. This study provides insight into how HPAI spread via PPHaV movement and scientific basis for control strategies.

Conserving evolutionarily distinct species is critical to safeguard human well-being

Although there is growing interest in safeguarding the Tree of Life to preserve the human benefits that are directly provided by biodiversity, their evolutionary distribution remains unknown, which has hampered our understanding of the potential of phylodiversity indicators to evince them. Here, I drew on a global review of plant benefits and comprehensive phylogenetic information to breakdown their evolutionary distribution and thereby show why the commonly used Phylogenetic Diversity and Evolutionary Distinctiveness indicators can unequivocally help to preserve these natural services. Beneficial species clumped within phylogenetically overdispersed genera and closely related species often contributed very few and redundant benefits, suggesting that multiple plant lineages are required to maintain a wide variety of services. Yet, a reduced number of species stood out as multi-beneficial and evolutionarily distinct plants relative to both the entire phylogeny and the subset of beneficial species, and they collectively contributed a higher-than-expected number of records for most types of benefits. In addition to providing a clear mechanistic understanding for the recently proved success of Phylogenetic Diversity in capturing plant benefits, these findings stress the decisive role that conservation programmes aimed at protecting evolutionarily distinct taxa will play in safeguarding the beneficial potential of biodiversity for the future.

Phylogeny-Guided Microbiome OTU-Specific Association Test (POST)

Background: The relationship between host conditions and microbiome profiles, typically characterized by operational taxonomic units (OTUs), contains important information about the microbial role in human health. Traditional association testing frameworks are challenged by the high-dimensionality and sparsity of typical microbiome profiles. Incorporating phylogenetic information is often used to address these challenges with the assumption that evolutionarily similar taxa tend to behave similarly. However, this assumption may not always be valid due to the complex effect of microbes, and phylogenetic information should be incorporated in a data-supervised fashion. Results: In this work, we propose a local collapsing test called Phylogeny-guided microbiome OTU-Specific association Test (POST). In POST, whether or not to borrow information and how much information to borrow from the neighboring OTUs in the phylogenic tree are supervised by phylogenetic distance and the outcome-OTU association. POST is constructed under the kernel machine framework to accommodate complex OTU effects and extends kernel machine microbiome tests from community-level to OTU-level. Using simulation studies, we showed that when the phylogenetic tree is informative, POST has better performance than existing OTU-level association tests. When the phylogenetic tree is not informative, POST achieves similar performance as existing methods. Finally, we show that POST can identify more outcome-associated OTUs that are of biological relevance in real data applications on bacterial vaginosis and on preterm birth. Conclusions: Using POST, we show that the power of detecting associated microbiome features can be enhanced by adaptively leveraging the phylogenetic information when testing for a target OTU. We developed an user friendly R package POSTm which is now available at CRAN (https://CRAN.R-project.org/package=POSTm) for public access.

1438Dynamics of inter-farm transmission of highly pathogenic avian influenza integrating vehicle movements and phylogenetic information

Background Highly pathogenic avian influenza (HPAI), a zoonotic infectious disease, has been considered a severe threat to public health. The fundamental prevention and control strategy against HPAI includes minimizing the outbreaks of poultry holdings where the virus primarily spreads through animal trade and poultry production associated vehicle movement (PPVM). However, very few attempts have been made to elucidate the association between PPVM and HPAI transmission compared to studies on poultry trade. Therefore, our study aimed to elucidate the role of PPVM on HPAI transmission. Methods We performed network analysis using PPVM data based on a global positioning system (GPS), with phylogenetic information of the HPAI virus for reliable estimation. Moreover, the contribution of PPVM to HPAI infection was estimated by Bayesian inference. Results The network analysis revealed that the connection via PPVM between the same genetic group of infected premises (IPs) was more prevalent than that of different genotype IPs. Moreover, the similarity of farm poultry species and the overlapped integrators between inter-linked IPs was associated with potential transmission route formation. Additionally, the contribution of PPVM among phylogenetically clustered IPs was estimated to have 28.25% of HPAI infections in IPs on average. Conclusions HPAI control strategies including targeted movement restriction and standstill should be established against the HPAI transmission via PPVM. Key messages This is a solid and novel study depicting the need for combining epidemiological analysis with data regarding molecular epidemiology of pathogens.

Measuring Phylogenetic Information of Incomplete Sequence Data

Widely used approaches for extracting phylogenetic information from aligned sets of molecular sequences rely upon probabilistic models of nucleotide substitution or amino-acid replacement. The phylogenetic information that can be extracted depends on the number of columns in the sequence alignment and will be decreased when the alignment contains gaps due to insertion or deletion events. Motivated by the measurement of information loss, we suggest assessment of the Effective Sequence Length (ESL) of an aligned data set. The ESL can differ from the actual number of columns in a sequence alignment because of the presence of alignment gaps. Furthermore, the estimation of phylogenetic information is affected by model misspecification. Inevitably, the actual process of molecular evolution differs from the probabilistic models employed to describe this process. This disparity means the amount of phylogenetic information in an actual sequence alignment will differ from the amount in a simulated data set of equal size, which motivated us to develop a new test for model adequacy. Via theory and empirical data analysis, we show how to disentangle the effects of gaps and model misspecification. By comparing the Fisher information of actual and simulated sequences, we identify which alignment sites and tree branches are most affected by gaps and model misspecification.

Morphology of the pupae of Physonota humilis Boheman and Physonota stigmatilis Boheman (Coleoptera: Chrysomelidae: Cassidinae: Ischyrosonychini)

Pupal morphology has been described for 11 species in six genera of the Neotropical tortoise beetle tribe, Ischyrosonychini Chapuis, 1875. This life stage may offer valuable phylogenetic information but more pupae need to be documented. The pupae of Physonota humilis Boheman, 1856 and P. stigmatilis Boheman, 1854 are described and illustrated for the first time. The pupa of P. humilis does not exhibit lateral scoli on the abdominal segments. Additionally, the body surfaces of P. humilis and P. stigmatilis are somewhat tuberculate, different from other described Physonota Boheman, 1854 pupae. A key to the described pupae of Physonota is provided and 10 phylogenetic character hypotheses are proposed.

Combining taxonomic, phylogenetic and functional diversity reveals new global priority areas for tetrapod conservation

We are in the midst of a sixth mass extinction but little is known about the global patterns of biodiversity when accounting for taxonomic, phylogenetic and functional information. Here, we present the first integrated analysis of global variation in taxonomic, functional diversity and phylogenetic diversity of more than 17,000 tetrapod species (terrestrial mammals, amphibians, reptiles and birds). We used a new metric (z-Diversity) able to synthetize taxonomic, functional and phylogenetic information across different sets of species to provide a comprehensive estimation of biodiversity. Our analyses reveal that hotspots of tetrapod diversity are clustered in specific regions of the world such as central Africa and the Indian peninsula, and that climate stability and energy availability have an overarching importance in explaining tetrapod spatial patterns. Future research might take advantage of these methods to perform an informed prioritization of protected areas.