Correlated evolution of dispersal traits and habitat preference in the melicgrasses

PREMISE

Seed dispersal is a critical process impacting individual plants and their communities. Plants have evolved numerous strategies and structures to disperse their seeds, but the evolutionary drivers of this diversity remain poorly understood in most lineages. We tested the hypothesis that the evolution of wind dispersal traits within the melicgrasses (Poaceae: Meliceae Link ex Endl.) was correlated with occupation of open and disturbed habitats.

METHODS

To evaluate wind dispersal potential, we collected seed dispersal structures (diaspores) from 24 melicgrass species and measured falling velocity and estimated dispersal distances. Species' affinity for open and disturbed habitats were recorded using georeferenced occurrence records and land cover maps. To test whether habitat preference and dispersal traits were correlated, we used phylogenetically informed multilevel models.

RESULTS

Melicgrasses display several distinct morphologies associated with wind dispersal, suggesting likely convergence. Open habitat taxa had slower-falling diaspores, consistent with increased wind dispersal potential. However, their shorter stature meant that dispersal distances, at a given wind speed, were not higher than those of their forest-occupying relatives. Species with affinities for disturbed sites had slower-falling diaspores and greater wind dispersal distances, largely explained by lighter diaspores.

CONCLUSIONS

Our results are consistent with the hypothesized evolutionary relationship between habitat preference and dispersal strategy. However, phylogenetic inertia and other plant functions (e.g., water conservation) likely shaped dispersal trait evolution in melicgrasses. It remains unclear if dispersal trait changes were precipitated by or predated changing habitat preferences. Nevertheless, our study provides promising results and a framework for disentangling dispersal strategy evolution.

On the origin of bird's nest fungi: Phylogenomic analyses of fungi in the Nidulariaceae (Agaricales, Basidiomycota)

Nidulariaceae, also known as bird's nest fungi, is an understudied group of mushroom-forming fungi. The common name is derived from their nest-like morphology. Bird's nest fungi are ubiquitous wood decomposers or saprobes on dung. Recent studies showed that species in the Nidulariaceae form a monophyletic group with five sub-clades. However, phylogenetic relationships among genera and placement of Nidulariaceae are still unclear. We present phylogenomic analyses of bird's nest fungi and related Agaricales fungi to gain insight into the evolution of Nidulariaceae. A species tree with 17 newly generated genomes of bird's nest fungi and representatives from all major clades of Agaricales was constructed using 1044 single-copy genes to explore the intergeneric relationships and pinpoint the placement of Nidulariaceae within Agaricales. We corroborated the hypothesis that bird's nest fungi are sister to Squamanitaceae, which includes mushroom-shaped fungi with a stipe and pileus that are saprobes and mycoparasites. Lastly, stochastic character mapping of discrete traits on phylogenies (SIMMAP) suggests that the ancestor of bird's nest fungi likely possessed an evanescent, globose peridium without strings attaching to the spore packets (funiculi). This analysis suggests that the funiculus was gained twice and that the persistent, cupulate peridium form was gained at least four times and lost once. However, alternative coding schemes and datasets with a wider array of Agaricales produced conflicting results during ancestral state reconstruction, indicating that there is some uncertainty in the number of peridium transitions and that taxon sampling may significantly alter ancestral state reconstructions. Overall, our results suggest that several key morphological characters of Nidulariaceae have been subject to homoplasy.

The evolution of ephemeral flora in Xinjiang, China: insights from plastid phylogenomic analyses of Brassicaceae

BACKGROUND

The ephemeral flora of northern Xinjiang, China, plays an important role in the desert ecosystems. However, the evolutionary history of this flora remains unclear. To gain new insights into its origin and evolutionary dynamics, we comprehensively sampled ephemeral plants of Brassicaceae, one of the essential plant groups of the ephemeral flora.

RESULTS

We reconstructed a phylogenetic tree using plastid genomes and estimated their divergence times. Our results indicate that ephemeral species began to colonize the arid areas in north Xinjiang during the Early Miocene and there was a greater dispersal of ephemeral species from the surrounding areas into the ephemeral community of north Xinjiang during the Middle and Late Miocene, in contrast to the Early Miocene or Pliocene periods.

CONCLUSIONS

Our findings, together with previous studies, suggest that the ephemeral flora originated in the Early Miocene, and species assembly became rapid from the Middle Miocene onwards, possibly attributable to global climate changes and regional geological events.

Identifying and addressing methodological incongruence in phylogenomics: A review

The availability of phylogenetic data has greatly expanded in recent years. As a result, a new era in phylogenetic analysis is dawning-one in which the methods we use to analyse and assess our data are the bottleneck to producing valuable phylogenetic hypotheses, rather than the need to acquire more data. This makes the ability to accurately appraise and evaluate new methods of phylogenetic analysis and phylogenetic artefact identification more important than ever. Incongruence in phylogenetic reconstructions based on different datasets may be due to two major sources: biological and methodological. Biological sources comprise processes like horizontal gene transfer, hybridization and incomplete lineage sorting, while methodological ones contain falsely assigned data or violations of the assumptions of the underlying model. While the former provides interesting insights into the evolutionary history of the investigated groups, the latter should be avoided or minimized as best as possible. However, errors introduced by methodology must first be excluded or minimized to be able to conclude that biological sources are the cause. Fortunately, a variety of useful tools exist to help detect such misassignments and model violations and to apply ameliorating measurements. Still, the number of methods and their theoretical underpinning can be overwhelming and opaque. Here, we present a practical and comprehensive review of recent developments in techniques to detect artefacts arising from model violations and poorly assigned data. The advantages and disadvantages of the different methods to detect such misleading signals in phylogenetic reconstructions are also discussed. As there is no one-size-fits-all solution, this review can serve as a guide in choosing the most appropriate detection methods depending on both the actual dataset and the computational power available to the researcher. Ultimately, this informed selection will have a positive impact on the broader field, allowing us to better understand the evolutionary history of the group of interest.

Pedomorphosis in the ancestry of marsupial mammals

Within mammals, different reproductive strategies (e.g., egg laying, live birth of extremely underdeveloped young, and live birth of well-developed young) have been linked to divergent evolutionary histories. How and when developmental variation across mammals arose is unclear. While egg laying is unquestionably considered the ancestral state for all mammals, many long-standing biases treat the extreme underdeveloped state of marsupial young as the ancestral state for therian mammals (clade including both marsupials and placentals), with the well-developed young of placentals often considered the derived mode of development. Here, we quantify mammalian cranial morphological development and estimate ancestral patterns of cranial shape development using geometric morphometric analysis of the largest comparative ontogenetic dataset of mammals to date (165 specimens, 22 species). We identify a conserved region of cranial morphospace for fetal specimens, after which cranial morphology diversified through ontogeny in a cone-shaped pattern. This cone-shaped pattern of development distinctively reflected the upper half of the developmental hourglass model. Moreover, cranial morphological variation was found to be significantly associated with the level of development (position on the altricial-precocial spectrum) exhibited at birth. Estimation of ancestral state allometry (size-related shape change) reconstructs marsupials as pedomorphic relative to the ancestral therian mammal. In contrast, the estimated allometries for the ancestral placental and ancestral therian were indistinguishable. Thus, from our results, we hypothesize that placental mammal cranial development most closely reflects that of the ancestral therian mammal, while marsupial cranial development represents a more derived mode of mammalian development, in stark contrast to many interpretations of mammalian evolution.

Whole-genome analyses disentangle reticulate evolution of primroses in a biodiversity hotspot

Biodiversity hotspots, such as the Caucasus mountains, provide unprecedented opportunities for understanding the evolutionary processes that shape species diversity and richness. Therefore, we investigated the evolution of Primula sect. Primula, a clade with a high degree of endemism in the Caucasus. We performed phylogenetic and network analyses of whole-genome resequencing data from the entire nuclear genome, the entire chloroplast genome, and the entire heterostyly supergene. The different characteristics of the genomic partitions and the resulting phylogenetic incongruences enabled us to disentangle evolutionary histories resulting from tokogenetic vs cladogenetic processes. We provide the first phylogeny inferred from the heterostyly supergene that includes all species of Primula sect. Primula. Our results identified recurrent admixture at deep nodes between lineages in the Caucasus as the cause of non-monophyly in Primula. Biogeographic analyses support the 'out-of-the-Caucasus' hypothesis, emphasizing the importance of this hotspot as a cradle for biodiversity. Our findings provide novel insights into causal processes of phylogenetic discordance, demonstrating that genome-wide analyses from partitions with contrasting genetic characteristics and broad geographic sampling are crucial for disentangling the diversification of species-rich clades in biodiversity hotspots.

A phylogenomically informed five-order system for the closest relatives of land plants

The evolution of streptophytes had a profound impact on life on Earth. They brought forth those photosynthetic eukaryotes that today dominate the macroscopic flora: the land plants (Embryophyta).¹ There is convincing evidence that the unicellular/filamentous Zygnematophyceae-and not the morphologically more elaborate Coleochaetophyceae or Charophyceae-are the closest algal relatives of land plants.^2-6 Despite the species richness (>4,000), wide distribution, and key evolutionary position of the zygnematophytes, their internal phylogeny remains largely unresolved.^7,8 There are also putative zygnematophytes with interesting body plan modifications (e.g., filamentous growth) whose phylogenetic affiliations remain unknown. Here, we studied a filamentous green alga (strain MZCH580) from an Austrian peat bog with central or parietal chloroplasts that lack discernible pyrenoids. It represents Mougeotiopsis calospora PALLA, an enigmatic alga that was described more than 120 years ago⁹ but never subjected to molecular analyses. We generated transcriptomic data of M. calospora strain MZCH580 and conducted comprehensive phylogenomic analyses (326 nuclear loci) for 46 taxonomically diverse zygnematophytes. Strain MZCH580 falls in a deep-branching zygnematophycean clade together with some unicellular species and thus represents a formerly unknown zygnematophycean lineage with filamentous growth. Our well-supported phylogenomic tree lets us propose a new five-order system for the Zygnematophyceae and provides evidence for at least five independent origins of true filamentous growth in the closest algal relatives of land plants. This phylogeny provides a robust and comprehensive framework for performing comparative analyses and inferring the evolution of cellular traits and body plans in the closest relatives of land plants.

The diverse terminology of reptile eggshell microstructure and its effect on phylogenetic comparative analyses

Reptile eggshell ensures water and gas exchange during incubation and plays a key role in reproductive success. The diversity of reptilian incubation and life history strategies has led to many clade-specific structural adaptations of their eggshell, which have been studied in extant taxa (i.e. birds, crocodilians, turtles, and lepidosaurs). Most studies on non-avian eggshells were performed over 30 years ago and categorized reptile eggshells into two main types: "hard" and "soft" - sometimes with a third intermediate category, "semi-rigid." In recent years, however, debate over the evolution of eggshell structure of major reptile clades has revealed how definitions of hard and soft eggshells influence inferred deep-time evolutionary patterns. Here, we review the diversity of extant and fossil eggshell with a focus on major reptile clades, and the criteria that have been used to define hard, soft, and semi-rigid eggshells. We show that all scoring approaches that retain these categories discretize continuous quantitative traits (e.g. eggshell thickness) and do not consider independent variation of other functionally important microstructural traits (e.g. degree of calcification, shell unit inner structure). We demonstrate the effect of three published approaches to discretizing eggshell type into hard, semi-rigid, and soft on ancestral state reconstructions using 200+ species representing all major extant and extinct reptile clades. These approaches result in different ancestral states for all major clades including Archosauria and Dinosauria, despite a difference in scoring for only 1-4% of the sample. Proposed scenarios of reptile eggshell evolution are highly conditioned by sampling, tree calibration, and lack of congruence between definitions of eggshell type. We conclude that the traditional "soft/hard/semi-rigid" classification of reptilian eggshells should be abandoned and provide guidelines for future descriptions focusing on specific functionally relevant characteristics (e.g. inner structures of shell units, pores, and membrane elements), analyses of these traits in a phylogenetic context, and sampling of previously undescribed taxa, including fossil eggs.

Phylogeny and evolution of morphological structures in a highly diverse lineage of fruiting-body-forming amoebae, order Trichiales (Myxomycetes, Amoebozoa)

Early phylogenetic studies refuted most previous assumptions concerning the evolution of the morphological traits in the fruiting bodies of the order Trichiales and did not detect discernible evolutionary patterns, yet they were based on a limited number of species. We infer a new Trichiales phylogeny based on three independently inherited genetic regions (nuclear and mitochondrial), with a fair taxonomic sampling encompassing its broad diversity. Besides, we study the evolutionary history of some key morphological characters. According to the new phylogeny, most fruiting body traits in Trichiales systematics do not represent exclusive synapomorphies or autapomorphies for most monophyletic groups. Instead, the evolution of the features derived from the peridium, stalk, capillitium, and spores showed intricate patterns, and character state transitions occurred rather within- than between clades. Thus, we should consider other evolutionary scenarios instead of assuming the homology of some characters. According to these results, we propose a new classification of Trichiales, including the creation of a new genus, Gulielmina, the resurrection of the family Dictydiaethaliaceae and the genus Ophiotheca, and the proporsal of 13 new combinations for species of the genera Arcyria (1), Hemitrichia (2), Ophiotheca (2), Oligonema (4), Gulielmina (3), and Perichaena (1).

Modeling Pulsed Evolution and Time-independent Variation Improves the Confidence Level of Ancestral and Hidden State Predictions

Ancestral state reconstruction is not only a fundamental tool for studying trait evolution, but also very useful for predicting the unknown trait values (hidden states) of extant species. A well-known problem in ancestral and hidden state predictions is that the uncertainty associated with predictions can be so large that predictions themselves are of little use. Therefore, for meaningful interpretation of predicted traits and hypothesis testing, it is prudent to accurately assess the uncertainty of the predictions. Commonly used constant-rate Brownian motion (BM) model fails to capture the complexity of tempo and mode of trait evolution in nature, making predictions under the BM model vulnerable to lack-of-fit errors from model misspecification. Using empirical data (mammalian body size and bacterial genome size), we show that the distribution of residual Z-scores under the BM model is neither homoscedastic nor normal as expected. Consequently, the 95% confidence intervals (CIs) of predicted traits are so unreliable that the actual coverage probability ranges from 33% (strongly permissive) to 100% (strongly conservative). Alternative methods such as BayesTraits and StableTraits that allow variable rates in evolution improve the predictions but are computationally expensive. Here we develop RasperGade, a method of ancestral and hidden state prediction that uses the Levy process to explicitly model gradual evolution, pulsed evolution and time-independent variation. Using the same empirical data, we show that RasperGade outperforms both BayesTraits and StableTraits in providing reliable confidence estimates and is orders-of-magnitude faster. Our results suggest that, when predicting the ancestral and hidden states of continuous traits, the rate variation should always be assessed and the quality of confidence estimates should always be examined.

Prey life-history influences the evolution of egg mass and indirectly reproductive investment in a group of free-living insect predators

The balance between risk and benefit of exploiting resources drives life-history evolution in organisms. Predators are naturally recognized as major drivers of the life-history evolution of their prey. Although prey may also influence the life-history evolution of their predators in the context of an evolutionary arms race, there is far more evidence of the role of predators than of prey.The goal of this study was to investigate the role of prey in life-history evolution of predators using ladybird beetle predators of aphids and coccids. These particular ladybirds and their prey were chosen because literature shows that the pace of life of aphids is faster than that of coccids and this difference is reflected in the life histories of the ladybirds that specialize on feeding on aphids or coccids.Thirty-four species of ladybird predators of aphids and eight of coccids belonging to five different tribes were collected and reared in the laboratory. The females were weighed as well as their eggs, and their reproductive investment estimated as the number of ovarioles. Phylogenetic relatedness was controlled for in the statistical analyses.Controlling for female mass revealed that ladybird predators of aphids lay bigger eggs than ladybird predators of coccids. This difference is not influenced by phylogenetic relatedness but only by the type of prey eaten. We suggest that ladybird predators of coccids lay smaller eggs because neonate larvae do not have to search, catch, and subdue prey. Both types of ladybirds have a similar reproductive investment relative to their body mass when phylogeny is controlled for.Recognizing the influence of prey on the life-history evolution of predators is important for understanding food web dynamics. From an applied perspective, this fine evolutionary tuning of prey-predator relationships should be used to guide and increase the efficiency of biological control programs.

Exploring the role of climatic niche changes in the evolution of the southern South American genus Baripus (Coleoptera: Carabidae): optimization of non-hereditary climatic variables and phylogenetic signal measurement

Baripus is a ground beetle genus endemic to southern South America, currently distributed across grassland and shrub habitats in mountain and lowland regions. The species of this genus are known to have been affected by the Andean orogeny and the climate changes that occurred during this process. In this study, we seek to understand how the orogeny of the Andes may have led to changes in the climatic niches of the species of Baripus over time. We integrated former ecological and historical biogeographic hypotheses, exploring the use of parsimony optimization of phylogenetically structured climate variables and ancestral character state reconstruction methods. We then performed regression analyses of the optimized climatic niche variables within the phylogenetic tree of Baripus. We were able to infer significant climatic niche constraints, and niche changes that provide new insights to the existing knowledge, supporting former ecological and biogeographic hypotheses for this genus. Such trends in climatic niche could be explained by the rain shadow effect caused by the Andean uplift as well as with other climate shifts associated with temperature and precipitation swings that occurred in this region from the Middle Miocene to the Pliocene.

Repeated dietary shifts in elapid snakes (Squamata: Elapidae) revealed by ancestral state reconstruction

Identifying the traits of ancestral organisms can reveal patterns and drivers of organismal diversification. Unfortunately, reconstructing complex multistate traits (such as diet) remains challenging. Adopting a ‘reconstruct, then aggregate’ approach in a maximum likelihood framework, we reconstructed ancestral diets for 298 species of elapid snakes. We tested whether different prey types were correlated with one another, tested for one-way contingency between prey type pairs, and examined the relationship between snake body size and dietary composition. We demonstrate that the evolution of diet was characterized by niche conservation punctuated by repeated dietary shifts. The ancestor of elapids most likely fed on reptiles and possibly amphibians, with deviations from this ancestral diet occurring repeatedly due to shifts into marine environments and changes in body size. Moreover, we demonstrate important patterns of prey use, including one-way dependency—most obviously the inclusion of eggs being dependent on a diet that already included the producers of those eggs. Despite imperfect dietary data, our approach produced a robust overview of dietary evolution. Given the paucity of natural history information for many organisms, our approach has the potential to increase the number of lineages to which ancestral state reconstructions of multistate traits can be robustly applied.

An Integrative Bioinformatic Analysis for Keratinase Detection in Marine-Derived Streptomyces

Keratinases present promising biotechnological applications, due to their ability to degrade keratin. Streptomyces appears as one of the main sources of these enzymes, but complete genome sequences of keratinolytic bacteria are still limited. This article reports the complete genomes of three marine-derived streptomycetes that show different levels of feather keratin degradation, with high (strain G11C), low (strain CHD11), and no (strain Vc74B-19) keratinolytic activity. A multi-step bioinformatics approach is described to explore genes encoding putative keratinases in these genomes. Despite their differential keratinolytic activity, multiplatform annotation reveals similar quantities of ORFs encoding putative proteases in strains G11C, CHD11, and Vc74B-19. Comparative genomics classified these putative proteases into 140 orthologous groups and 17 unassigned orthogroup peptidases belonging to strain G11C. Similarity network analysis revealed three network communities of putative peptidases related to known keratinases of the peptidase families S01, S08, and M04. When combined with the prediction of cellular localization and phylogenetic reconstruction, seven putative keratinases from the highly keratinolytic strain Streptomyces sp. G11C are identified. To our knowledge, this is the first multi-step bioinformatics analysis that complements comparative genomics with phylogeny and cellular localization prediction, for the prediction of genes encoding putative keratinases in streptomycetes.

A phylogenomic supermatrix of Galliformes (Landfowl) reveals biased branch lengths

Building taxon-rich phylogenies is foundational for macroevolutionary studies. One approach to improve taxon sampling beyond individual studies is to build supermatricies of publicly available data, incorporating taxa sampled across different studies and utilizing different loci. Most existing supermatrix studies have focused on loci commonly sequenced with Sanger technology ("legacy" markers, such as mitochondrial data and small numbers of nuclear loci). However, incorporating phylogenomic studies into supermatrices allows problem nodes to be targeted and resolved with considerable amounts of data, while improving taxon sampling with legacy data. Here we estimate phylogeny from a galliform supermatrix which includes well-known model and agricultural species such as the chicken and turkey. We assembled a supermatrix comprising 4500 ultra-conserved elements (UCEs) collected as part of recent phylogenomic studies in this group and legacy mitochondrial and nuclear (intron and exon) sequences. Our resulting phylogeny included 88% of extant species and recovered well-accepted relationships with strong support. However, branch lengths, which are particularly important in down-stream macroevolutionary studies, appeared vastly skewed. Taxa represented only by rapidly evolving mitochondrial data had high proportions of missing data and exhibited long terminal branches. Conversely, taxa sampled for slowly evolving UCEs with low proportions of missing data exhibited substantially shorter terminal branches. We explored several branch length re-estimation methods with particular attention to terminal branches and conclude that re-estimation using well-sampled mitochondrial sequences may be a pragmatic approach to obtain trees suitable for macroevolutionary analysis.

Genomic footprints of repeated evolution of CAM photosynthesis in a Neotropical species radiation

The adaptive radiation of Bromeliaceae (pineapple family) is one of the most diverse among Neotropical flowering plants. Diversification in this group was facilitated by shifts in several adaptive traits or "key innovations" including the transition from C₃ to CAM photosynthesis associated with xeric (heat/drought) adaptation. We used phylogenomic approaches, complemented by differential gene expression (RNA-seq) and targeted metabolite profiling, to address the mechanisms of C₃ /CAM evolution in the extremely species-rich bromeliad genus, Tillandsia, and related taxa. Evolutionary analyses of whole-genome sequencing and RNA-seq data suggest that evolution of CAM is associated with coincident changes to different pathways mediating xeric adaptation in this group. At the molecular level, C₃ /CAM shifts were accompanied by gene expansion of XAP5 CIRCADIAN TIMEKEEPER homologs, a regulator involved in sugar- and light-dependent regulation of growth and development. Our analyses also support the re-programming of abscisic acid-related gene expression via differential expression of ABF2/ABF3 transcription factor homologs, and adaptive sequence evolution of an ENO2/LOS2 enolase homolog, effectively tying carbohydrate flux to abscisic acid-mediated abiotic stress response. By pinpointing different regulators of overlapping molecular responses, our results suggest plausible mechanistic explanations for the repeated evolution of correlated adaptive traits seen in a textbook example of an adaptive radiation.

Evolutionary bedfellows: Reconstructing the ancestral state of autotomy and regeneration

Some form of regeneration occurs in all lifeforms and extends from single-cell organisms to humans. The degree to which regenerative ability is distributed across different taxa, however, is harder to ascertain given the potential for phylogenetic constraint or inertia, and adaptive processes to shape this pattern. Here, we examine the phylogenetic history of regeneration in two groups where the trait has been well-studied: arthropods and reptiles. Because autotomy is often present alongside regeneration in these groups, we performed ancestral state reconstructions for both traits to more precisely assess the timing of their origins and the degree to which these traits coevolve. Using an ancestral trait reconstruction, we find that autotomy and regeneration were present at the base of the arthropod and reptile trees. We also find that when autotomy is lost it does not re-evolve easily. Lastly, we find that the distribution of regeneration is intimately connected to autotomy with the association being stronger in reptiles than in arthropods. Although these patterns suggest that decoupling autotomy and regeneration at a broad phylogenetic scale may be difficult, the available data provides useful insight into their entanglement. Ultimately, our reconstructions provide the important groundwork to explore how selection may have played a role during the loss of regeneration in specific lineages.

Phylogeny and character evolution in the Dacrymycetes, and systematics of Unilacrymaceae and Dacryonaemataceae fam. nov

We present a multilocus phylogeny of the class Dacrymycetes, based on data from the 18S, ITS, 28S, RPB1, RPB2, TEF-1α, 12S, and ATP6 DNA regions, with c. 90 species including the types of most currently accepted genera. A variety of methodological approaches was used to infer phylogenetic relationships among the Dacrymycetes, from a supermatrix strategy using maximum likelihood and Bayesian inference on a concatenated dataset, to coalescence-based calculations, such as quartet-based summary methods of independent single-locus trees, and Bayesian integration of single-locus trees into a species tree under the multispecies coalescent. We evaluate for the first time the taxonomic usefulness of some cytological phenotypic characters, i.e., vacuolar contents (vacuolar bodies and lipid bodies), number of nuclei of recently discharged basidiospores, and pigments, with especial emphasis on carotenoids. These characters, along with several others traditionally used for the taxonomy of this group (basidium shape, presence and morphology of clamp connections, morphology of the terminal cells of cortical/marginal hyphae, presence and degree of ramification of the hyphidia), are mapped on the resulting phylogenies and their evolution through the class Dacrymycetes discussed. Our analyses reveal five lineages that putatively represent five different families, four of which are accepted and named. Three out of these four lineages correspond to previously circumscribed and published families (Cerinomycetaceae, Dacrymycetaceae, and Unilacrymaceae), and one is proposed as the new family Dacryonaemataceae. Provisionally, only a single order, Dacrymycetales, is accepted within the class. Furthermore, the systematics of the two smallest families, Dacryonaemataceae and Unilacrymaceae, are investigated to the species level, using coalescence-based species delimitation on multilocus DNA data, and a detailed morphological study including morphometric analyses of the basidiospores. Three species are accepted in Dacryonaema, the type, Da. rufum, the newly combined Da. macnabbii (basionym Dacrymyces macnabbii), and a new species named Da. macrosporum. Two species are accepted in Unilacryma, the new U. bispora, and the type, U. unispora, the latter treated in a broad sense pending improved sampling across the Holarctic.

Accuracy of ancestral state reconstruction for non-neutral traits

The assumptions underpinning ancestral state reconstruction are violated in many evolutionary systems, especially for traits under directional selection. However, the accuracy of ancestral state reconstruction for non-neutral traits is poorly understood. To investigate the accuracy of ancestral state reconstruction methods, trees and binary characters were simulated under the BiSSE (Binary State Speciation and Extinction) model using a wide range of character-state-dependent rates of speciation, extinction and character-state transition. We used maximum parsimony (MP), BiSSE and two-state Markov (Mk2) models to reconstruct ancestral states. Under each method, error rates increased with node depth, true number of state transitions, and rates of state transition and extinction; exceeding 30% for the deepest 10% of nodes and highest rates of extinction and character-state transition. Where rates of character-state transition were asymmetrical, error rates were greater when the rate away from the ancestral state was largest. Preferential extinction of species with the ancestral character state also led to higher error rates. BiSSE outperformed Mk2 in all scenarios where either speciation or extinction was state dependent and outperformed MP under most conditions. MP outperformed Mk2 in most scenarios except when the rates of character-state transition and/or extinction were highly asymmetrical and the ancestral state was unfavoured.

Sequence selection by FitSS4ASR alleviates ancestral sequence reconstruction as exemplified for geranylgeranylglyceryl phosphate synthase

For evolutionary studies, but also for protein engineering, ancestral sequence reconstruction (ASR) has become an indispensable tool. The first step of every ASR protocol is the preparation of a representative sequence set containing at most a few hundred recent homologs whose composition determines decisively the outcome of a reconstruction. A common approach for sequence selection consists of several rounds of manual recompilation that is driven by embedded phylogenetic analyses of the varied sequence sets. For ASR of a geranylgeranylglyceryl phosphate synthase, we additionally utilized FitSS4ASR, which replaces this time-consuming protocol with an efficient and more rational approach. FitSS4ASR applies orthogonal filters to a set of homologs to eliminate outlier sequences and those bearing only a weak phylogenetic signal. To demonstrate the usefulness of FitSS4ASR, we determined experimentally the oligomerization state of eight predecessors, which is a delicate and taxon-specific property. Corresponding ancestors deduced in a manual approach and by means of FitSS4ASR had the same dimeric or hexameric conformation; this concordance testifies to the efficiency of FitSS4ASR for sequence selection. FitSS4ASR-based results of two other ASR experiments were added to the Supporting Information. Program and documentation are available at https://gitlab.bioinf.ur.de/hek61586/FitSS4ASR.

Phylogenetic comparative analysis supports aposematic colouration-body size association in millipede assassins (Hemiptera: Reduviidae: Ectrichodiinae)

The diversity of colour patterns and its importance in interactions with the environment make colouration in animals an intriguing research focus. Aposematic colouration is positively correlated with body size in certain groups of animals, suggesting that warning colours are more effective or that crypsis is harder to achieve in larger animals. Surprisingly, this relationship has not been recovered in studies investigating insects, which may have been confounded by a focus on aposematic taxa that are also gregarious. Millipede assassin bugs (Hemiptera: Reduviidae: Ectrichodiinae) comprise species with cryptic and aposematic colour patterns across a range of body sizes, are typically solitary as adults and are thus an excellent model for investigating a possible association between colouration and body size. Here, we use a comprehensive phylogeny for Ectrichodiinae, ancestral state reconstruction of colouration, and phylogenetic comparative methods to test for a colouration-body size association. The ancestor of Ectrichodiinae is reconstructed as cryptically coloured, with multiple subsequent transitions between aposematic and cryptic colouration. Aposematic colouration is positively associated with male body length and supports the hypothesis that selection on Ectrichodiinae body size may influence evolutionary transitions between aposematic and cryptic colouration or alternatively that selection for aposematic colouration influences body size evolution.

Challenges of comprehensive taxon sampling in comparative biology: Wrestling with rosids

Using phylogenetic approaches to test hypotheses on a large scale, in terms of both species sampling and associated species traits and occurrence data-and doing this with rigor despite all the attendant challenges-is critical for addressing many broad questions in evolution and ecology. However, application of such approaches to empirical systems is hampered by a lingering series of theoretical and practical bottlenecks. The community is still wrestling with the challenges of how to develop species-level, comprehensively sampled phylogenies and associated geographic and phenotypic resources that enable global-scale analyses. We illustrate difficulties and opportunities using the rosids as a case study, arguing that assembly of biodiversity data that is scale-appropriate-and therefore comprehensive and global in scope-is required to test global-scale hypotheses. Synthesizing comprehensive biodiversity data sets in clades such as the rosids will be key to understanding the origin and present-day evolutionary and ecological dynamics of the angiosperms.

Genomic Signature of an Avian Lilliput Effect across the K-Pg Extinction

Survivorship following major mass extinctions may be associated with a decrease in body size-a phenomenon called the Lilliput Effect. Body size is a strong predictor of many life history traits (LHTs), and is known to influence demography and intrinsic biological processes. Pronounced changes in organismal size throughout Earth history are therefore likely to be associated with concomitant genome-wide changes in evolutionary rates. Here, we report pronounced heterogeneity in rates of molecular evolution (varying up to $\sim$20-fold) across a large-scale avian phylogenomic data set and show that nucleotide substitution rates are strongly correlated with body size and metabolic rate. We also identify potential body size reductions associated with the Cretaceous-Paleogene (K-Pg) transition, consistent with a Lilliput Effect in the wake of that mass extinction event. We posit that selection for reduced body size across the K-Pg extinction horizon may have resulted in transient increases in substitution rate along the deepest branches of the extant avian tree of life. This "hidden" rate acceleration may result in both strict and relaxed molecular clocks over-estimating the age of the avian crown group through the relationship between life history and demographic parameters that scale with molecular substitution rate. If reductions in body size (and/or selection for related demographic parameters like short generation times) are a common property of lineages surviving mass extinctions, this phenomenon may help resolve persistent divergence time debates across the tree of life. Furthermore, our results suggest that selection for certain LHTs may be associated with deterministic molecular evolutionary outcomes.

Know your farmer: Ancient origins and multiple independent domestications of ambrosia beetle fungal cultivars

Bark and ambrosia beetles are highly specialized weevils (Curculionidae) that have established diverse symbioses with fungi, most often from the order Ophiostomatales (Ascomycota, Sordariomycetes). The two types of beetles are distinguished by their feeding habits and intimacy of interactions with their symbiotic fungi. The tree tissue diet of bark beetles is facilitated by fungi, while ambrosia beetles feed solely on fungi that they farm. The farming life history strategy requires domestication of a fungus, which the beetles consume as their sole food source. Ambrosia beetles in the subfamily Platypodinae originated in the mid-Cretaceous (119-88 Ma) and are the oldest known group of farming insects. However, attempts to resolve phylogenetic relationships and the timing of domestication events for fungal cultivars have been largely inconclusive. We sequenced the genomes of 12 ambrosia beetle fungal cultivars and bark beetle associates, including the devastating laurel wilt pathogen, Raffaelea lauricola, to estimate a robust phylogeny of the Ophiostomatales. We find evidence for contemporaneous diversification of the beetles and their associated fungi, followed by three independent domestication events of the ambrosia fungi genus Raffaelea. We estimate the first domestication of an Ophiostomatales fungus occurred ~86 Ma, 25 million years earlier than prior estimates and in close agreement with the estimated age of farming in the Platypodinae (96 Ma). Comparisons of the timing of fungal domestication events with the timing of beetle radiations support the hypothesis that the first large beetle radiations may have spread domesticated "ambrosia" fungi to other fungi-associated beetle groups, perhaps facilitating the evolution of new farming lineages.

Piphillin: Improved Prediction of Metagenomic Content by Direct Inference from Human Microbiomes

Functional analysis of a clinical microbiome facilitates the elucidation of mechanisms by which microbiome perturbation can cause a phenotypic change in the patient. The direct approach for the analysis of the functional capacity of the microbiome is via shotgun metagenomics. An inexpensive method to estimate the functional capacity of a microbial community is through collecting 16S rRNA gene profiles then indirectly inferring the abundance of functional genes. This inference approach has been implemented in the PICRUSt and Tax4Fun software tools. However, those tools have important limitations since they rely on outdated functional databases and uncertain phylogenetic trees and require very specific data pre-processing protocols. Here we introduce Piphillin, a straightforward algorithm independent of any proposed phylogenetic tree, leveraging contemporary functional databases and not obliged to any singular data pre-processing protocol. When all three inference tools were evaluated against actual shotgun metagenomics, Piphillin was superior in predicting gene composition in human clinical samples compared to both PICRUSt and Tax4Fun (p<0.01 and p<0.001, respectively) and Piphillin’s ability to predict disease associations with specific gene orthologs exhibited a 15% increase in balanced accuracy compared to PICRUSt. From laboratory animal samples, no performance advantage was observed for any one of the tools over the others and for environmental samples all produced unsatisfactory predictions. Our results demonstrate that functional inference using the direct method implemented in Piphillin is preferable for clinical biospecimens. Piphillin is publicly available for academic use at http://secondgenome.com/Piphillin.

The Impact of Reconstruction Methods, Phylogenetic Uncertainty and Branch Lengths on Inference of Chromosome Number Evolution in American Daisies (Melampodium, Asteraceae)

Chromosome number change (polyploidy and dysploidy) plays an important role in plant diversification and speciation. Investigating chromosome number evolution commonly entails ancestral state reconstruction performed within a phylogenetic framework, which is, however, prone to uncertainty, whose effects on evolutionary inferences are insufficiently understood. Using the chromosomally diverse plant genus Melampodium (Asteraceae) as model group, we assess the impact of reconstruction method (maximum parsimony, maximum likelihood, Bayesian methods), branch length model (phylograms versus chronograms) and phylogenetic uncertainty (topological and branch length uncertainty) on the inference of chromosome number evolution. We also address the suitability of the maximum clade credibility (MCC) tree as single representative topology for chromosome number reconstruction. Each of the listed factors causes considerable incongruence among chromosome number reconstructions. Discrepancies between inferences on the MCC tree from those made by integrating over a set of trees are moderate for ancestral chromosome numbers, but severe for the difference of chromosome gains and losses, a measure of the directionality of dysploidy. Therefore, reliance on single trees, such as the MCC tree, is strongly discouraged and model averaging, taking both phylogenetic and model uncertainty into account, is recommended. For studying chromosome number evolution, dedicated models implemented in the program ChromEvol and ordered maximum parsimony may be most appropriate. Chromosome number evolution in Melampodium follows a pattern of bidirectional dysploidy (starting from x = 11 to x = 9 and x = 14, respectively) with no prevailing direction.

Increased diversification rates follow shifts to bisexuality in liverworts

Shifts in sexual systems are one of the key drivers of species diversification. In contrast to angiosperms, unisexuality prevails in bryophytes. Here, we test the hypotheses that bisexuality evolved from an ancestral unisexual condition and is a key innovation in liverworts. We investigate whether shifts in sexual systems influence diversification using hidden state speciation and extinction analysis (HiSSE). This new method compares the effects of the variable of interest to the best-fitting latent variable, yielding robust and conservative tests. We find that the transitions in sexual systems are significantly biased toward unisexuality, even though bisexuality is coupled with increased diversification. Sexual systems are strongly conserved deep within the liverwort tree but become much more labile toward the present. Bisexuality appears to be a key innovation in liverworts. Its effects on diversification are presumably mediated by the interplay of high fertilization rates, massive spore production and long-distance dispersal, which may separately or together have facilitated liverwort speciation, suppressed their extinction, or both. Importantly, shifts in liverwort sexual systems have the opposite effect when compared to angiosperms, leading to contrasting diversification patterns between the two groups. The high prevalence of unisexuality among liverworts suggests, however, a strong selection for sexual dimorphism.

Rates of morphological evolution are heterogeneous in Early Cretaceous birds

The Early Cretaceous is a critical interval in the early history of birds. Exceptional fossils indicate that important evolutionary novelties such as a pygostyle and a keeled sternum had already arisen in Early Cretaceous taxa, bridging much of the morphological gap between Archaeopteryx and crown birds. However, detailed features of basal bird evolution remain obscure because of both the small sample of fossil taxa previously considered and a lack of quantitative studies assessing rates of morphological evolution. Here we apply a recently available phylogenetic method and associated sensitivity tests to a large data matrix of morphological characters to quantify rates of morphological evolution in Early Cretaceous birds. Our results reveal that although rates were highly heterogeneous between different Early Cretaceous avian lineages, consistent patterns of significantly high or low rates were harder to pinpoint. Nevertheless, evidence for accelerated evolutionary rates is strongest at the point when Ornithuromorpha (the clade comprises all extant birds and descendants from their most recent common ancestors) split from Enantiornithes (a diverse clade that went extinct at the end-Cretaceous), consistent with the hypothesis that this key split opened up new niches and ultimately led to greater diversity for these two dominant clades of Mesozoic birds.

Integrating fossils, phylogenies, and niche models into biogeography to reveal ancient evolutionary history: the case of Hypericum (hypericaceae)

In disciplines such as macroevolution that are not amenable to experimentation, scientists usually rely on current observations to test hypotheses about historical events, assuming that “the present is the key to the past.” Biogeographers, for example, used this assumption to reconstruct ancestral ranges from the distribution of extant species. Yet, under scenarios of high extinction rates, the biodiversity we observe today might not be representative of the historical diversity and this could result in incorrect biogeographic reconstructions. Here, we introduce a new approach to incorporate into biogeographic inference the temporal, spatial, and environmental information provided by the fossil record, as a direct evidence of the extinct biodiversity fraction. First, inferences of ancestral ranges for those nodes in the phylogeny calibrated with the fossil record are constrained to include the geographic distribution of the fossil. Second, we use fossil distribution and past climate data to reconstruct the climatic preferences and potential distribution of ancestral lineages over time, and use this information to build a biogeographic model that takes into account “ecological connectivity” through time. To show the power of this approach, we reconstruct the biogeographic history of the large angiosperm genus Hypericum, which has a fossil record extending back to the Early Cenozoic. Unlike previous reconstructions based on extant species distributions, our results reveal that Hypericum stem lineages were already distributed in the Holarctic before diversification of its crown-group, and that the geographic distribution of the genus has been relatively stable throughout the climatic oscillations of the Cenozoic. Geographical movement was mediated by the existence of climatic corridors, like Beringia, whereas the equatorial tropical belt acted as a climatic barrier, preventing Hypericum lineages to reach the southern temperate regions. Our study shows that an integrative approach to historical biogeography—that combines sources of evidence as diverse as paleontology, ecology, and phylogenetics—could help us obtain more accurate reconstructions of ancient evolutionary history. It also reveals the confounding effect different rates of extinction across regions have in biogeography, sometimes leading to ancestral areas being erroneously inferred as recent colonization events.

Molecular phylogeny of Harpactorini (Insecta: Reduviidae): correlation of novel predation strategy with accelerated evolution of predatory leg morphology

Much research and discussion have focused on the effects of key innovations on lineage diversification, whereas little has been done to investigate their role in morphological evolution using phylogenetic approaches. Here we present the first comprehensive molecular phylogeny of the Harpactorini (Insecta: Reduviidae), the largest assassin bug tribe, sampling 229 terminal taxa and using five gene segments (28S D2, D3-D5, 16S, COI, and Deformed). Employing comparative phylogenetic methods, we demonstrate the correlation of a putative key innovation, the sticky trap predation strategy, with accelerated rates of morphological evolution of the predatory fore leg in assassin bugs. We show that bugs exhibiting sticky trap predation have evolved more slender and longer fore femora than non-sticky bugs. Using phylogenetically independent contrast analyses, we document correlated evolution between femoral thickness and length. We argue that the novel sticky trap predation strategy may allow sticky bugs to alleviate functional constraints on the fore femur and thus to attain a higher rate of evolution than other Harpactorini or Reduviidae. We discuss the possibility that sticky bugs represent a case of adaptive radiation. We also test historical supra-generic groups within the Harpactorini, and show that most of them are not monophyletic. We confirm the paraphyly of Harpactorini with respect to Rhaphidosomini.

Evolution of life history and behavior in Hominidae: towards phylogenetic reconstruction of the chimpanzee-human last common ancestor

The origin of the fundamental behavioral differences between humans and our closest living relatives is one of the central issues of evolutionary anthropology. The prominent, chimpanzee-based referential model of early hominin behavior has recently been challenged on the basis of broad multispecies comparisons and newly discovered fossil evidence. Here, we argue that while behavioral data on extant great apes are extremely relevant for reconstruction of ancestral behaviors, these behaviors should be reconstructed trait by trait using formal phylogenetic methods. Using the widely accepted hominoid phylogenetic tree, we perform a series of character optimization analyses using 65 selected life-history and behavioral characters for all extant hominid species. This analysis allows us to reconstruct the character states of the last common ancestors of Hominoidea, Hominidae, and the chimpanzee-human last common ancestor. Our analyses demonstrate that many fundamental behavioral and life-history attributes of hominids (including humans) are evidently ancient and likely inherited from the common ancestor of all hominids. However, numerous behaviors present in extant great apes represent their own terminal autapomorphies (both uniquely derived and homoplastic). Any evolutionary model that uses a single extant species to explain behavioral evolution of early hominins is therefore of limited use. In contrast, phylogenetic reconstruction of ancestral states is able to provide a detailed suite of behavioral, ecological and life-history characters for each hypothetical ancestor. The living great apes therefore play an important role for the confident identification of the traits found in the chimpanzee-human last common ancestor, some of which are likely to represent behaviors of the fossil hominins.

Niche conservatism constrains Australian honeyeater assemblages in stressful environments

The hypothesis of phylogenetic niche conservatism proposes that most extant members of a clade remain in ancestral environments because expansion into new ecological space imposes a selectional load on a population. A prediction that follows is that local assemblages contain increasingly phylogenetically clustered subsets of species with increasing difference from the ancestral environment of a clade. We test this in Australian Meliphagidae, a continental radiation of birds that originated in wet, subtropical environments, but subsequently spread to drier environments as Australia became more arid during the late Cenozoic. We find local assemblages are increasingly phylogenetically clustered along a gradient of decreasing precipitation. The pattern is less clear along a temperature gradient. We develop a novel phyloclimatespace to visualise the expansion of some lineages into drier habitats. Although few species extend into arid regions, those that do occupy larger ranges and thus local species richness does not decline predictably with precipitation.

Mutualism with sea anemones triggered the adaptive radiation of clownfishes

Background

Adaptive radiation is the process by which a single ancestral species diversifies into many descendants adapted to exploit a wide range of habitats. The appearance of ecological opportunities, or the colonisation or adaptation to novel ecological resources, has been documented to promote adaptive radiation in many classic examples. Mutualistic interactions allow species to access resources untapped by competitors, but evidence shows that the effect of mutualism on species diversification can greatly vary among mutualistic systems. Here, we test whether the development of obligate mutualism with sea anemones allowed the clownfishes to radiate adaptively across the Indian and western Pacific oceans reef habitats.

Results

We show that clownfishes morphological characters are linked with ecological niches associated with the sea anemones. This pattern is consistent with the ecological speciation hypothesis. Furthermore, the clownfishes show an increase in the rate of species diversification as well as rate of morphological evolution compared to their closest relatives without anemone mutualistic associations.

Conclusions

The effect of mutualism on species diversification has only been studied in a limited number of groups. We present a case of adaptive radiation where mutualistic interaction is the likely key innovation, providing new insights into the mechanisms involved in the buildup of biodiversity. Due to a lack of barriers to dispersal, ecological speciation is rare in marine environments. Particular life-history characteristics of clownfishes likely reinforced reproductive isolation between populations, allowing rapid species diversification.

Trophic specialization influences the rate of environmental niche evolution in damselfishes (Pomacentridae)

The rate of environmental niche evolution describes the capability of species to explore the available environmental space and is known to vary among species owing to lineage-specific factors. Trophic specialization is a main force driving species evolution and is responsible for classical examples of adaptive radiations in fishes. We investigate the effect of trophic specialization on the rate of environmental niche evolution in the damselfish, Pomacentridae, which is an important family of tropical reef fishes. First, phylogenetic niche conservatism is not detected in the family using a standard test of phylogenetic signal, and we demonstrate that the environmental niches of damselfishes that differ in trophic specialization are not equivalent while they still overlap at their mean values. Second, we estimate the relative rates of niche evolution on the phylogenetic tree and show the heterogeneity among rates of environmental niche evolution of the three trophic groups. We suggest that behavioural characteristics related to trophic specialization can constrain the evolution of the environmental niche and lead to conserved niches in specialist lineages. Our results show the extent of influence of several traits on the evolution of the environmental niche and shed new light on the evolution of damselfishes, which is a key lineage in current efforts to conserve biodiversity in coral reefs.