Whole-genome analyses resolve early branches in the tree of life of modern birds

Genomes of critically endangered saola are shaped by population structure and purging

The saola is one of the most elusive large mammals, standing at the brink of extinction. We constructed a reference genome and resequenced 26 saola individuals, confirming the saola as a basal member of the Bovini. Despite its small geographic range, we found that the saola is partitioned into two populations with high genetic differentiation (FST = 0.49). We estimate that these populations diverged and started declining 5,000-20,000 years ago, possibly due to climate changes and exacerbated by increasing human activities. The saola has long tracts without genomic diversity; however, most of these tracts are not shared by the two populations. Saolas carry a high genetic load, yet their gradual decline resulted in the purging of the most deleterious genetic variation. Finally, we find that combining the two populations, e.g., in an eventual captive breeding program, would mitigate the genetic load and increase the odds of species survival.

Changes in Evolutionary Developmental Control Points in the Amniote Limb May Explain Hyperphalangy

Amniotes show a great diversity of limb phenotypes, including limbs specialized for running, flying, swimming, and digging. Here, we have examined how this diversity is generated during limb development in 13 species using transcriptomics and in situ hybridization. The selected species show evolutionary changes in the number of phalanges and/or loss of claws. We first looked at genes that show cyclical expression during digit development. Significantly, we find that Gdf5 cycles more rapidly in digits developing more phalanges. We identified two novel cyclically expressed genes: Ackr3 and Wnt9a. We also identified a transition point at which phalanx formation stops and claw development begins. We found that this transition point is marked by the downregulation of multiple developmental genes in the phalanx-forming region, and upregulation of claw-related genes. The timing of this transition is conserved, taking place at the same developmental stage in all digits of all species examined-except in the clawless digits of the Chinese soft-shelled turtle, the crocodilians, and birds. We suggest a model based on transcriptional heterochrony, in which the frequency of phalanx formation and the timing of the phalanx-claw transition are evolutionary control points open to natural selection on the phenotype. Furthermore, our model suggests that relaxation of developmental constraints on the timing of the phalanx-claw transition allows the digits to develop more phalanges (hyperphalangy). This is seen in some turtles, crocodilians, and dolphins. More broadly, our findings are consistent with the hypothesis that "hotspots" in otherwise conserved developmental pathways may be targets for evolutionary tinkering.

Range size and abundance dynamics of Japanese breeding birds over 40 years suggest a potential crisis in warm areas

Understanding the current status of biodiversity is crucial to preventing its loss in a changing world. We examined changes in the geographical range size and abundance of 165 bird species breeding in Japan during the past 40 years, as well as temperature niche changes in the past 20 years. Higher temperatures were recorded within the ranges of non-native species than in those of native species, and we detected range-size expansion and increased abundance among non-native species. Although open-land species exhibited range reductions from the 1970s to the 1990s, many recovered and the ranges of only a few species declined after this period. Nevertheless, the abundance of open-land species did decline, despite range-size recovery; similar inconsistencies were detected for waterbirds and raptors. Analysis of long-term temperatures suggested that species left warmest areas within their distributions while maximum temperatures experienced by species during the survey years did not change systematically. Birds in warm regions may be facing a crisis, with attrition of native bird communities and expansion of non-native species. It is necessary to establish efficient measures to prevent further expansions of non-native species and conservation measures of native species within managed areas in warm regions with few intact habitats.

Towards the next generation of species delimitation methods: an overview of machine learning applications

Species delimitation is the process of distinguishing between populations of the same species and distinct species of a particular group of organisms. Various methods exist for inferring species limits, whether based on morphological, molecular, or other types of data. In the case of methods based on DNA sequences, most of them are rooted in the coalescent theory. However, coalescence-based models have limitations, for instance regarding complex evolutionary scenarios and large datasets. In this context, machine learning (ML) can be considered as a promising analytical tool, and provides an effective way to explore dataset structures when species-level divergences are hypothesized. In this review, we examine the use of ML in species delimitation and provide an overview and critical appraisal of existing workflows. We also provide simple explanations on how the main types of ML approaches operate, which should help uninitiated researchers and students interested in the field. Our review suggests that while current ML methods designed to infer species limits are analytically powerful, they also present specific limitations and should not be considered as definitive alternatives to coalescent methods for species delimitation. Future ML enterprises to delimit species should consider the constraints related to the use of simulated data, as in other model-based methods relying on simulations. Conversely, the flexibility of ML algorithms offers a significant advantage by enabling the analysis of diverse data types (e.g., genetic and phenotypic) and handling large datasets effectively. We also propose best practices for the use of ML methods in species delimitation, offering insights into potential future applications. We expect that the proposed guidelines will be useful for enhancing the accessibility, effectiveness, and objectivity of ML in species delimitation.

The telomere-to-telomere genome assembly and annotation of the rock carp (Procypris rabaudi)

Procypris rabaudi, commonly known as the rock carp, is an endemic economic fish in the middle-upper reaches of the Yangtze River. To enhance the understanding the biology of rack carps, a high-quality reference genome is required in different areas of study. Here, we generated the first telomere-to-telomere genome assembly and annotation of the rock carp, which spans 1.64 Gb with a contig N50 of 32.36 Mb. Hi-C assembly suggested that 99.83% sequences were positioned to 50 pseudo-chromosomes. Notably, 43 chromosomes were assembled without any gap. Through the integration of homologous-based predictions and RNA-sequencing technology, we identified 44,402 protein-coding genes, with 43,663 of them (98.3%) predicted to be functional. Furthermore, our assembled genome achieved 98.1% BUSCO completeness. This work improves the quality of the rock carp genome and provides valuable foundation for the future studies of genomics, biology and the fishery resources breeding.

Accurate, scalable, and fully automated inference of species trees from raw genome assemblies using ROADIES

Current genome sequencing initiatives across a wide range of life forms offer significant potential to enhance our understanding of evolutionary relationships and support transformative biological and medical applications. Species trees play a central role in many of these applications; however, despite the widespread availability of genome assemblies, accurate inference of species trees remains challenging due to the limited automation, substantial domain expertise, and computational resources required by conventional methods. To address this limitation, we present ROADIES, a fully automated pipeline to infer species trees starting from raw genome assemblies. In contrast to the prominent approach, ROADIES incorporates a unique strategy of randomly sampling segments of the input genomes to generate gene trees. This eliminates the need for predefining a set of loci, limiting the analyses to a fixed number of genes, and performing the cumbersome gene annotation and/or whole genome alignment steps. ROADIES also eliminates the need to infer orthology by leveraging existing discordance-aware methods that allow multicopy genes. Using the genomic datasets from large-scale sequencing efforts across four diverse life forms (placental mammals, pomace flies, birds, and budding yeasts), we show that ROADIES infers species trees that are comparable in quality to the state-of-the-art studies but in a fraction of the time and effort, including on challenging datasets with rampant gene tree discordance and complex polyploidy. With its speed, accuracy, and automation, ROADIES has the potential to vastly simplify species tree inference, making it accessible to a broader range of scientists and applications.

Phylogenomics establishes an Early Miocene reconstruction of reef vertebrate diversity

Oceans blanket more than two-thirds of Earth's surface, yet marine biodiversity is disproportionately concentrated in coral reefs. Investigating the origins of this exceptional diversity is crucial for predicting how reefs will respond to anthropogenic disturbances. Here, we use a genome-scale dataset to reconstruct the evolutionary history of the wrasses and parrotfishes (Labridae), which rank among the most species-rich and ecologically diverse lineages of reef fishes. We show that major labrid clades experienced pulses of evolutionary innovation and accelerated diversification during the Miocene approximately 20 to 15 million years ago that the origin of no single phenotypic trait can explain. These results draw parallels to the evolutionary histories of many clades after mass extinctions and corroborate recent fossil evidence for an Early Miocene extinction event in oceanic vertebrates and changes in coral reef faunal composition. Our data provide genomic evidence for a major Early Miocene reassembly of reef faunas.

Drivers of avian genomic change revealed by evolutionary rate decomposition

Modern birds have diversified into a striking array of forms, behaviours and ecological roles. Analyses of molecular evolutionary rates can reveal the links between genomic and phenotypic change1-4, but disentangling the drivers of rate variation at the whole-genome scale has been difficult. Using comprehensive estimates of traits and evolutionary rates across a family-level phylogeny of birds5,6, we find that genome-wide mutation rates across lineages are predominantly explained by clutch size and generation length, whereas rate variation across genes is driven by the content of guanine and cytosine. Here, to find the subsets of genes and lineages that dominate evolutionary rate variation in birds, we estimated the influence of individual lineages on decomposed axes of gene-specific evolutionary rates. We find that most of the rate variation occurs along recent branches of the tree, associated with present-day families of birds. Additional tests on axes of rate variation show rapid changes in microchromosomes immediately after the Cretaceous-Palaeogene transition. These apparent pulses of evolution are consistent with major changes in the genetic machineries for meiosis, heart performance, and RNA splicing, surveillance and translation, and correlate with the ecological diversity reflected in increased tarsus length. Collectively, our analyses paint a nuanced picture of avian evolution, revealing that the ancestors of the most diverse lineages of birds underwent major genomic changes related to mutation, gene usage and niche expansion in the early Palaeogene period.

The many origins of extremophile fishes

Extremophiles survive in environments that are considered uninhabitable for most living things. The evolution of extremophiles is of great interest because of how they may have contributed to the assembly of ecosystems, yet the evolutionary dynamics that drive extremophile evolution remain obscure. Here, we investigate the evolution of extremophiles in Zoarcoidea, a lineage of over 300 species of fishes that have colonized both poles, the deep sea, and hydrothermal vents. We show that a pulse of habitat invasion occurred across over 20 different zoarcoid lineages within the last 8 million years, far after the origin of their prototypical innovation for surviving in cold water: type III antifreeze protein. Instead, a secondary burst of anatomical, physiological and life history traits and a handful of founder events in extreme ecosystems appear to have propelled zoarcoid diversification. These results decentralize the role of prototypical changes to organismal biology in shaping extremophile radiations and provide a clear example of how a combination of ancient adaptations and recent contingency shapes the origination of lineages in challenging habitats.

Unique myoglobin adaptation to endothermy and flight since the origin of birds

Myoglobin (Mb) mediates oxygen diffusion and storage in muscle tissue and thus is important for the energy utilization and activity of animals. Birds generally have a high body temperature, and most species also possess the capability of powered flight. Both of these require high levels of aerobic metabolism. Within endothermic mammals, bats also independently evolved flight. Although the functional evolution of myoglobins in deep-diving amniote vertebrates has been well-studied, the functional evolution of myoglobin since the origins of both birds and bats is unclear. Here, with Mb-coding sequences from >200 extant amniote species, we reconstructed ancestral sequences to estimate the functional properties of myoglobin through amniote evolution. A dramatic change in net surface charge on myoglobin occurred during the origin of Aves, which might have been driven by positively selected amino acid substitutions that occurred on the lineage leading to all birds. However, in bats, no change in net surface charge occurred and instead, the Mb genes show evidence of strong purifying selection. The increased net surface charge on bird myoglobins implies an adaptation to flight-related endothermic and higher body temperatures, possibly by reducing harmful protein aggregations. Different from the findings of net surface charge, myoglobins of extant birds show lower stability compared with other amniotes, which probably accelerates the rate of oxygen utilization in muscles. In bats and other mammals, higher stability of Mb may be an alternative pathway for adaptation to endothermy, indicating divergent evolution of myoglobin in birds and bats.

Comparative Phylogenetics Reveal Clade-specific Drivers of Recombination Rate Evolution Across Vertebrates

Meiotic recombination is an integral cellular process, required for the production of viable gametes. Recombination rate is a fundamental genomic parameter, modulating genomic responses to selection. Our increasingly detailed understanding of its molecular underpinnings raises the prospect that we can gain insight into trait divergence by examining the molecular evolution of recombination genes from a pathway perspective, as in mammals, where protein-coding changes in later stages of the recombination pathway are connected to divergence in intra-clade recombination rate. Here, we leverage increased availability of avian and teleost genomes to reconstruct the evolution of the recombination pathway across two additional vertebrate clades: birds, which have higher and more variable rates of recombination and similar divergence times to mammals, and teleost fish, which have much deeper divergence times. Rates of molecular evolution of recombination genes are highly correlated between vertebrate clades and significantly elevated compared to control panels, suggesting that they experience similar selective pressures. Avian recombination genes are significantly more likely to exhibit signatures of positive selection than other clades, unrestricted to later stages of the pathway. Signatures of positive selection in genes linked to recombination rate variation in mammalian populations and those with signatures of positive selection across the avian phylogeny are highly correlated. In contrast, teleost fish recombination genes have significantly less evidence of positive selection despite high intra-clade recombination rate variability. Gaining clade-specific understanding of patterns of variation in recombination genes can elucidate drivers of recombination rate and thus, factors influencing genetic diversity, selection efficacy, and species divergence.

Ex-situ avian sex skews: determinants and implications for conservation

With over half of all avian species in decline globally, zoo-based recovery programs are increasingly relied upon to save species from extinction. The success of such programs not only rests with political will, but also in our understanding of species' breeding biology and how individuals and populations respond to changes in their environment. Sex skews, that is, an imbalance in the optimal number of males to females, is an underlying mechanism of population decline in some threatened species. Ex-situ (i.e., zoo-based) management practices will need to become more efficient to support the growing number of conservation reliant species and manage sex skews to amend, repair and restore population stability both in- and ex-situ. In this article, we analysed data from over 182,000 birds in global ex-situ collections. We interpreted sex ratio variation by observing the proportion of males within and between orders, International Union for Conservation of Nature (IUCN) threat status and housing inside and outside of a species' natural range. Overall, our results showed that male-biased sex skews are more prevalent ex-situ than they are in the wild and although they vary greatly at the institutional level, were closer to parity at a global level. The variation amongst threat status and housing outside of range were less significant. These findings have implications for the conservation management of threatened birds and the potential of ex-situ populations to function with maximum effect in an integrated management system.

How do big brains evolve?

In both birds and mammals, variation in brain size predominantly reflects variation in mass or volume of the pallium (neocortex) and, to a lesser extent, of the cerebellum, suggesting convergent coevolution of brains and cognition. When brain measures are based on neuron counts, however, a surprisingly different picture emerges: The number of neurons in the cerebellum surpasses those in the pallium of all mammals (including humans and other primates) and in many but not all birds studied to date. In particular, parrots and corvids, clades known for cognitive abilities that match those of primates, have brains that contain more pallial than cerebellar neurons. Birds and mammals may thus have followed different evolutionary routes of pallial-cerebellar coordination behind enhanced cognitive complexity.

Bayesian Selection of Relaxed-Clock Models: Distinguishing between Independent and Autocorrelated Rates

In Bayesian molecular-clock dating of species divergences, rate models are used to construct the prior on the molecular evolutionary rates for branches in the phylogeny, with independent and autocorrelated rate models being commonly used. The two classes of models, however, can result in markedly different divergence time estimates for the same data set, and thus selecting the best rate model appears important for obtaining reliable inferences of divergence times. However, the properties of Bayesian rate model selection are not well understood, in particular when the number of sequence partitions analyzed increases and when age calibrations (such as fossil calibrations) are misspecified. Furthermore, Bayesian rate model selection is computationally expensive as it requires the calculation of marginal likelihoods by Markov Chain Monte Carlo sampling, and therefore, methods that can speed up the model selection procedure without compromising its accuracy are desirable. In this study, we use a combination of computer simulations and real data analysis to investigate the statistical behavior of Bayesian rate model selection and we also explore approximations of the likelihood to improve computational efficiency in large phylogenomic data sets. Our simulations demonstrate that the posterior probability for the correct rate model converges to one as more molecular sequence partitions are analyzed and when no calibrations are used, as expected due to asymptotic Bayesian model selection theory. Furthermore, we also show the model selection procedure is robust to slight misspecification of calibrations, and reliable inference of the correct rate model is possible in this case. However, we show that when calibrations are seriously misspecified, calculated model probabilities are completely wrong and may converge to one for the wrong rate model. Finally, we demonstrate that approximating the phylogenetic likelihood under an arcsine branch-length transform can dramatically reduce the computational cost of rate model selection without compromising accuracy. We test the approximate procedure on two large phylogenies of primates (372 species) and flowering plants (644 species), replicating results obtained on smaller data sets using exact likelihood. Our findings and methodology can assist users in selecting the optimal rate model for estimating times and rates along the Tree of Life.

Robustness of divergence time estimation despite gene tree estimation error: a case study of fireflies (Coleoptera: Lampyridae)

Genomic data have become ubiquitous in phylogenomic studies, including divergence time estimation, but provide new challenges. These challenges include, among others, biological gene tree discordance, methodological gene tree estimation error, and computational limitations on performing full Bayesian inference under complex models. In this study, we use a recently published firefly (Coleoptera: Lampyridae) anchored hybrid enrichment data set (AHE; 436 loci for 88 Lampyridae species and 10 outgroup species) as a case study to explore gene tree estimation error and the robustness of divergence time estimation. First, we explored the amount of model violation using posterior predictive simulations because model violations are likely to bias phylogenetic inferences and produce gene tree estimation error. We specifically focused on missing data (either uniformly distributed or systematically) and the distribution of highly variable and conserved sites (either uniformly distributed or clustered). Our assessment of model adequacy showed that standard phylogenetic substitution models are not adequate for any of the 436 AHE loci. We tested if the model violations and alignment errors resulted indeed in gene tree estimation error by comparing the observed gene tree discordance to simulated gene tree discordance under the multispecies coalescent model. Thus, we show that the inferred gene tree discordance is not only due to biological mechanism but primarily due to inference errors. Lastly, we explored if divergence time estimation is robust despite the observed gene tree estimation error. We selected four subsets of the full AHE data set, concatenated each subset and performed a Bayesian relaxed clock divergence estimation in RevBayes. The estimated divergence times overlapped for all nodes that are shared between the topologies. Thus, divergence time estimation is robust using any well selected data subset as long as the topology inference is robust.

Cellular Scaling Rules for Brains of the Galliform Birds (Aves, Galliformes) Compared to Those of Songbirds and Parrots: Distantly Related Avian Lineages Have Starkly Different Neuronal Cerebrotypes

INTRODUCTION

Songbirds, especially corvids, and parrots are remarkably intelligent. Their cognitive skills are on par with primates and their brains contain primate-like numbers of neurons concentrated in high densities in the telencephalon. Much less is known about cognition and neuron counts in more basal bird lineages. Here, we focus on brain cellular composition of galliform birds, which have small brains relative to body size and a proportionally small telencephalon and are often perceived as cognitively inferior to most other birds.

METHODS

We use the isotropic fractionator to assess quantitatively the numbers and distributions of neurons and nonneuronal cells in 15 species of galliform birds and compare their cellular scaling rules with those of songbirds, parrots, marsupials, insectivores, rodents, and primates.

RESULTS

On average, the brains of galliforms contain about half the number of neurons found in parrot and songbird brains of the same mass. Moreover, in contrast to these birds, galliforms resemble mammals in having small telencephalic and dominant cerebellar neuronal fractions. Consequently, galliforms have much smaller absolute numbers of neurons in their forebrains than equivalently sized songbirds and parrots, which may limit their cognitive abilities. However, galliforms have similar neuronal densities and neuron counts in the brain and forebrain as equally sized non-primate mammals. Therefore, it is not surprising that cognitive abilities of galliforms are on par with non-primate mammals in many domains.

CONCLUSION

Comparisons performed in this study demonstrate that birds representing distantly related clades markedly differ in neuronal densities, neuron numbers, and the allocation of brain neurons to major brain divisions. In analogy with the concept of volumetric composition of the brain, known as the cerebrotype, we conclude that distantly related birds have distinct neuronal cerebrotypes.

Descriptive histological analysis of the upper, lower, and third eyelids and the conjunctiva-associated lymphoid tissue in birds of prey

BACKGROUD

In this study, we present data obtained using light microscopy for the histological analysis of the upper eyelid (palpebra dorsalis), lower eyelid (palpebra ventralis), and third eyelid (palpebra tertia) also known as the nictitating membrane. We characterized the organized conjunctiva-associated lymphoid tissue (CALT) in selected raptor species. The aim of this study is to compare the histological structures of these eyelids in owls and diurnal raptors to identify potential evolutionary links or independent adaptations to their environments.

MATERIALS AND METHODS

We examined 34 individuals from 18 species representing Accipitriformes, Falconiformes, and Strigiformes, sourced from the Wrocław Zoological Garden (Poland), private bird collections (Poland), and wild birds found dead in the natural environment (Poland). The study involved morphometric analysis of the length and thickness of the tarsal plate of the lower eyelid. Microscopic evaluation included histological staining, using Masson-Goldner trichrome, Mayer's hematoxylin & eosin, Movat pentachrome (modified Russell-Movat), and picro-Mallory trichrome.

RESULTS

The structure of the eyelids in the analyzed bird orders proved to be highly diverse in terms of the presence of common features. The third eyelid, as well as CALT, exhibited the most variations morphological structures among the analyzed species. Strigiformes emerged as the most distinctive group of raptors, characterized by the greatest differences in eyelid morphology. This group of birds is not only distinct from other raptors but also internally diverse, with many significant differences observed between individual owl species.

CONCLUSION

Despite some common features, the upper, lower, and third eyelids of raptors from the orders Accipitriformes, Falconiformes, and Strigiformes exhibit significant morphological variation. The third eyelid and conjunctiva-associated lymphoid tissue (CALT) display the most diverse structures among the analyzed species. Owls stand out as a group of raptors with the most distinct eyelid morphologies, both compared to other raptors and within their own group. The small number of birds may lead to difficulties in distinguishing individual variation from species-specific traits, as we cannot be certain whether the observed differences result from genetic or environmental factors specific to the individual birds or if they are traits typical for the species. To address this issue, further studies involving a larger number of individuals from the same species are necessary to more accurately determine whether the morphological traits described in this study are consistent within the species or if significant variation exists among individuals.

Minus the Error: Testing for Positive Selection in the Presence of Residual Alignment Errors

Positive selection is an evolutionary process which increases the frequency of advantageous mutations because they confer a fitness benefit. Inferring the past action of positive selection on protein-coding sequences is fundamental for deciphering phenotypic diversity and the emergence of novel traits. With the advent of genome-wide comparative genomic datasets, researchers can analyze selection not only at the level of individual genes but also globally, delivering systems-level insights into evolutionary dynamics. However, genome-scale datasets are generated with automated pipelines and imperfect curation that does not eliminate all sequencing, annotation, and alignment errors. Positive selection inference methods are highly sensitive to such errors. We present BUSTED-E: a method designed to detect positive selection for amino acid diversification while concurrently identifying some alignment errors. This method builds on the flexible branch-site random effects model (BUSTED) for fitting distributions of dN/dS, with a critical modification: it incorporates an "error-sink" component to represent an abiological evolutionary regime. Using several genome-scale biological datasets that were extensively filtered using state-of-the art automated alignment tools, we show that BUSTED-E identifies pervasive residual alignment errors, produces more realistic estimates of positive selection, reduces bias, and improves biological interpretation. The BUSTED-E model promises to be a more stringent filter to identify positive selection in genome-wide contexts, thus enabling further characterization and validation of the most biologically relevant cases.

wQFM-TREE: highly accurate and scalable quartet-based species tree inference from gene trees

Motivation

methods are becoming increasingly popular for species tree estimation from multi-locus data in the presence of gene tree discordance. Accurate Species TRee Algorithm (ASTRAL), a leading method in this class, solves the Maximum Quartet Support Species Tree problem within a constrained solution space, while heuristics like Weighted Quartet Fiduccia-Mattheyses (wQFM) and Weighted Quartet MaxCut (wQMC) use weighted quartets and a divide-and-conquer strategy. Recent studies showed wQFM to be more accurate than ASTRAL and wQMC, though its scalability is hindered by the computational demands of explicitly generating and weighting Θ ( n 4 ) quartets. Here, we introduce wQFM-TREE, a novel summary method that enhances wQFM by avoiding explicit quartet generation and weighting, enabling its application to large datasets.

Results

Extensive simulations under diverse and challenging model conditions, with hundreds or thousands of taxa and genes, consistently demonstrate that wQFM-TREE matches or improves upon the accuracy of ASTRAL. It outperformed ASTRAL in 25 of 27 model conditions (statistically significant in 20) involving 200-1000 taxa. Moreover, applying wQFM-TREE to re-analyze the green plant dataset from the One Thousand Plant Transcriptomes Initiative produced a tree highly congruent with established evolutionary relationships of plants. wQFM-TREE's remarkable accuracy and scalability make it a strong competitor to leading methods. Its algorithmic and combinatorial innovations also enhance quartet-based computations, advancing phylogenetic estimation.

Availability and implementation

wQFM-TREE is freely available in open source form at https://github.com/abdur-rafi/wQFM-TREE.

New insights on angiosperm crown age based on Bayesian node dating and skyline fossilized birth-death approaches

Despite considerable work in recent years, pinpointing the time when angiosperms originated has been challenging. However, the rapid development of molecular clock methodology has provided new tools to resolve this conundrum. In particular, the fossilized birth-death model establishes a rich interplay between molecules and stratigraphy by incorporating fossils explicitly into dating analyses. In this study, we apply Bayesian node dating and the skyline fossilized birth-death model, which differ in how the calibration is applied, to estimate the crown age of angiosperms. Node dating analyses with different calibration strategies show that the posterior distribution is strongly constrained by the effective prior at the node of crown angiosperms, dominated by the maximum age constraint. Using the skyline fossilized birth-death model, we reveal that assigning different priors for origin time resulted in similar crown ages for angiosperms. Moreover, the oldest fossils play a significant role in time estimates, and the dating results are robust to sampling assumptions of extant taxa. Our dating analyses indicate a largely Triassic crown age (255-202 Ma) for angiosperms, the period when mammals, dinosaurs, and squamate reptiles first appeared, and highlight the potential of morphological data to redefine the timeline of angiosperms.

An Extensive Survey of Vertebrate-specific, Nonvisual Opsins Identifies a Novel Subfamily, Q113-Bistable Opsin

A group of nonvisual opsins specific to vertebrates is essential to understand evolution of lateral eyes, one of the most prominent innovations in this lineage. Nevertheless, our knowledge of their evolutionary history remains limited. To develop an integrated view of their evolution, we surveyed these non-visual opsins (VA opsin, pinopsin, parapinopsin, parietopsin, and parapinopsin-like) in 451 vertebrate genomes. Through extensive manual curation, we completed a high-quality catalog. We could not find them in 202 mammals, supporting previous reports of their loss. VA opsins are highly conserved among nonmammals. In contrast, other opsin subfamilies experienced more dynamic molecular evolution with many secondary losses. In addition, we found a previously unreported opsin subfamily that we named Q113-Bistable (QB) opsin. We found its orthologs only in several lizards and the tuatara. Nevertheless, QB opsin pseudogenes were discovered in diverse taxa, including ray-finned fishes, indicating its ancient origin. QB opsin, parapinopsin, and parietopsin are extremely prone to be lost in the course of evolution, and loss events involving these opsins seem to occur concomitantly. Furthermore, we demonstrated the spectral properties of QB opsin as a UV-sensitive, bistable photo-pigment. This study provides the first integrated view of the entire evolutionary history of this group of opsins.

CASTER: Direct species tree inference from whole-genome alignments

Genomes contain mosaics of discordant evolutionary histories, challenging the accurate inference of the tree of life. Although genome-wide data are routinely used for discordance-aware phylogenomic analyses, because of modeling and scalability limitations, the current practice leaves out large chunks of genomes. As more high-quality genomes become available, we urgently need discordance-aware methods to infer the tree directly from a multiple genome alignment. In this study, we introduce Coalescence-Aware Alignment-Based Species Tree Estimator (CASTER), a theoretically justified site-based method that eliminates the need to predefine recombination-free loci. CASTER is scalable to hundreds of mammalian whole genomes. We demonstrate the accuracy and scalability of CASTER in simulations that include recombination and apply CASTER to several biological datasets, showing that its per-site scores can reveal both biological and artifactual patterns of discordance across the genome.

Species tree branch length estimation despite incomplete lineage sorting, duplication, and loss

Phylogenetic branch lengths are essential for many analyses, such as estimating divergence times, analyzing rate changes, and studying adaptation. However, true gene tree heterogeneity due to incomplete lineage sorting (ILS), gene duplication and loss (GDL), and horizontal gene transfer (HGT) can complicate the estimation of species tree branch lengths. While several tools exist for estimating the topology of a species tree addressing various causes of gene tree discordance, much less attention has been paid to branch length estimation on multi-locus datasets. For single-copy gene trees, some methods are available that summarize gene tree branch lengths onto a species tree, including coalescent-based methods that account for heterogeneity due to ILS. However, no such branch length estimation method exists for multi-copy gene family trees that have evolved with gene duplication and loss. To address this gap, we introduce the CASTLES-Pro algorithm for estimating species tree branch lengths while accounting for both GDL and ILS. CASTLES-Pro improves on the existing coalescent-based branch length estimation method CASTLES by increasing its accuracy for single-copy gene trees and extends it to handle multi-copy ones. Our simulation studies show that CASTLES-Pro is generally more accurate than alternatives, eliminating the systematic bias toward overestimating terminal branch lengths often observed when using concatenation. Moreover, while not theoretically designed for HGT, we show that CASTLES-Pro maintains relatively high accuracy under high rates of random HGT.

Code availability

CASTLES-Pro is implemented inside the software package ASTER, available at https://github.com/chaoszhang/ASTER .

Data availability

The datasets and scripts used in this study are available at https://github.com/ytabatabaee/CASTLES-Pro-paper .

Chromosome-level reference genome assembly of the gyrfalcon (Falco rusticolus) and population genomics offer insights into the falcon population in Mongolia

The taxonomic classification of a falcon population found in the Mongolian Altai region in Asia has been heavily debated for two centuries and previous studies have been inconclusive, hindering a more informed conservation approach. Here, we generated a chromosome-level gyrfalcon reference genome using the Vertebrate Genomes Project (VGP) assembly pipeline. Using whole genome sequences of 49 falcons from different species and populations, including "Altai" falcons, we analyzed their population structure, admixture patterns, and demographic history. We find that the Altai falcons are genomic mosaics of saker and gyrfalcon ancestries, and carry distinct W and mitochondrial haplotypes that cluster with the lanner falcon. The Altai maternally-inherited haplotypes diverged 422,000 years before present (290,000-550,000 YBP) from the ancestor of sakers and gyrfalcons, both of which, in turn, split 109,000 YBP (70,000-150,000 YBP). The Altai W chromosome has 31 coding variants in 29 genes that may possibly influence important structural, behavioral, and reproductive traits. These findings provide insights into the question of Altai falcons as a candidate distinct species.

ZW Sex Chromosome Differentiation in Palaeognathous Birds Is Associated with Mitochondrial Effective Population Size but Not Mitochondrial Genome Size or Mutation Rate

Eukaryotic genome size varies considerably, even among closely related species. The causes of this variation are unclear, but weak selection against supposedly costly "extra" genomic sequences has been central to the debate for over 50 years. The mutational hazard hypothesis, which focuses on the increased mutation rate to null alleles in superfluous sequences, is particularly influential, though challenging to test. This study examines the sex chromosomes and mitochondrial genomes of 15 flightless or semiflighted palaeognathous bird species. In this clade, the nonrecombining portion of the W chromosome has independently expanded stepwise in multiple lineages. Given the shared maternal inheritance of the W chromosome and mitochondria, theory predicts that mitochondrial effective population size (Ne) should decrease due to increased Hill-Robertson interference in lineages with expanded nonrecombining W regions. Our findings support the extent of the nonrecombining W region with three indicators of reduced selective efficiency: (i) the ratio of nonsynonymous to synonymous nucleotide changes in the mitochondrion, (ii) the probability of radical amino acid changes, and (iii) the number of ancient, W-linked genes lost through evolution. Next, we tested whether reduced Ne affects mitochondrial genome size, as predicted by weak selection against genome expansion. We find no support for a relationship between mitochondrial genome size and expanded nonrecombining W regions, nor with increased mitochondrial mutation rates (predicted to modulate selective costs). These results highlight the utility of nonrecombining regions and mitochondrial genomes for studying genome evolution and challenge the general idea of a negative relation between the efficacy of selection and genome size.

The Genomic Landscape, Causes, and Consequences of Extensive Phylogenomic Discordance in Murine Rodents

A species tree is a central concept in evolutionary biology whereby a single branching phylogeny reflects relationships among species. However, the phylogenies of different genomic regions often differ from the species tree. Although tree discordance is widespread in phylogenomic studies, we still lack a clear understanding of how variation in phylogenetic patterns is shaped by genome biology or the extent to which discordance may compromise comparative studies. We characterized patterns of phylogenomic discordance across the murine rodents-a large and ecologically diverse group that gave rise to the laboratory mouse and rat model systems. Combining recently published linked-read genome assemblies for seven murine species with other available rodent genomes, we first used ultraconserved elements (UCEs) to infer a robust time-calibrated species tree. We then used whole genomes to examine finer-scale patterns of discordance across ∼12 million years of divergence. We found that proximate chromosomal regions tended to have more similar phylogenetic histories. There was no clear relationship between local tree similarity and recombination rates in house mice, but we did observe a correlation between recombination rates and average similarity to the species tree. We also detected a strong influence of linked selection whereby purifying selection at UCEs led to appreciably less discordance. Finally, we show that assuming a single species tree can result in substantial deviation from the results with gene trees when testing for positive selection under different models. Collectively, our results highlight the complex relationship between phylogenetic inference and genome biology and underscore how failure to account for this complexity can mislead comparative genomic studies.

Genomic evidence for hybridization and introgression between blue peafowl and endangered green peafowl and molecular foundation of leucistic plumage of blue peafowl

INTRODUCTION

The blue peafowl (Pavo cristatus) and the green peafowl (Pavo muticus) have garnered significant public affection due to their stunning appearance, although the green peafowl is currently endangered. The causative mutation that causes the leucistic plumage of the blue peafowl (also called white peafowl) remains unknown.

RESULTS

In this study, we generated a chromosome-level reference genome of the blue peafowl with a contig N50 of 30.6 Mb, including the autosomes, Z and W sex chromosomes, and a complete mitochondria DNA sequence. Data from 77 peafowl whole genomes, 76 peafowl mitochondrial genomes, and 33 peafowl W chromosomes genomes provided the first substantial genetic evidence for recent hybridization between green peafowls and blue peafowls. We found 3 hybrid green peafowls in zoo samples rather than in the wild samples, with a blue peafowl genomic content of 16-34%. Maternal genetic analysis showed that 2 of the hybrid female green peafowls contained complete blue peafowl mitochondrial genomes and W chromosomes. Some animal protection agencies release captive green peafowls in order to maintain the wild population of green peafowls. Therefore, to better protect the endangered green peafowl, we suggest that purebred identification must be carried out before releasing green peafowls from zoos into the wild in order to prevent the hybrid green peafowl from contaminating the wild green peafowl. In addition, we also found that there were historical introgression events of green peafowl to blue peafowl in 4 zoo blue peafowl individuals. The introgressed genomic regions contain IGFBP1 and IGFBP3 genes that could affect blue peafowl body size. Finally, we identified that the nonsense mutation (g.4:12583552G>A) in the EDNRB2 gene is the genetic causative mutation for leucistic plumage of blue peafowl, preventing melanocytes from being transported into plumage, thereby inhibiting melanin deposition.

CONCLUSION

Our research provides both theoretical and empirical support for the conservation of the endangered green peafowl. The high-quality genome and genomic data also provide a valuable resource for blue peafowl genomics-assisted breeding.

Pronounced early differentiation underlies zebra finch gonadal germ cell development

The diversity of germ cell developmental strategies has been well documented across many vertebrate clades. However, much of our understanding of avian primordial germ cell (PGC) specification and differentiation has derived from only one species, the chicken (Gallus gallus). Of the three major classes of birds, chickens belong to Galloanserae, representing less than 4% of species, while nearly 95% of extant bird species belong to Neoaves. This represents a significant gap in our knowledge of germ cell development across avian species, hampering efforts to adapt genome editing and reproductive technologies developed in chicken to other birds. We therefore applied single-cell RNA sequencing to investigate inter-species differences in germ cell development between chicken and zebra finch (Taeniopygia castanotis), a Neoaves songbird species and a common model of vocal learning. Analysis of early embryonic male and female gonads revealed the presence of two distinct early germ cell types in zebra finch and only one in chicken. Both germ cell types expressed zebra finch Germline Restricted Chromosome (GRC) genes, present only in songbirds among birds. One of the zebra finch germ cell types expressed the canonical PGC markers, as did chicken, but with expression differences in several signaling pathways and biological processes. The second zebra finch germ cell cluster was marked by proliferation and fate determination markers, indicating beginning of differentiation. Notably, these two zebra finch germ cell populations were present in both male and female zebra finch gonads as early as HH25. Using additional chicken developmental stages, similar germ cell heterogeneity was identified in the more developed gonads of females, but not males. Overall, our study demonstrates a substantial heterochrony in zebra finch germ cell development compared to chicken, indicating a richer diversity of avian germ cell developmental strategies than previously known.

Whence the birds: 200 years of dinosaurs, avian antecedents

Among the most revolutionary insights emerging from 200 years of research on dinosaurs is that the clade Dinosauria is represented by approximately 11 000 living species of birds. Although the origin of birds among dinosaurs has been reviewed extensively, recent years have witnessed tremendous progress in our understanding of the deep evolutionary origins of numerous distinctive avian anatomical systems. These advances have been enabled by exciting new fossil discoveries, leading to an ever-expanding phylogenetic framework with which to pinpoint the origins of characteristic avian features. The present review focuses on four notable avian systems whose Mesozoic evolutionary history has been greatly clarified by recent discoveries: brain, kinetic palate, pectoral girdle and postcranial skeletal pneumaticity.

Complete Mitochondrial Genomes of Pluvialis fulva and Charadrius dubius with Phylogenetic Analysis of Charadriiformes

BACKGROUND

Plovers (Charadriidae), within the order of Charadriiformes, a group of modern birds distributed worldwide, are a frequent subject of molecular phylogenetic studies. While research on mitochondrial genome (mitogenome) variation within the family Charadriidae, especially intraspecific variation, is limited. Additionally, the monophyly of Charadrius and the phylogenetic placement of Pluvialis remain contentious. Nevertheless, recent studies utilizing complete mitogenomes from available databases to construct phylogenetic trees for Charadriidae and Charadriiformes remain scarce.

METHODS

This study aims to explore mitogenome variation within Charadrius dubius and clarify the phylogenetic placement of Pluvialis fulva. We sequenced the complete mitogenome of six C. dubius and one P. fulva, and all additional available mitogenomes were integrated within Charadriiformes. The average complete mitogenome length of C. dubius is 16,889 bp, and P. fulva is 16,859 bp.

RESULTS

Our results support the suggestion that the monophyly of Charadrius and P. fulva is nested within Charadriidae. The phylogenetic analysis of Charadriiformes based on mitogenomes strongly supports the recognition of three major shorebird clades: Charadrii, Lari and Scolopaci, with Lari and Scolopaci identified as sister clades.

CONCLUSIONS

Our study reinforces the credibility of the inferred evolutionary relationships within Charadriidae and Charadriiformes.

Museomics Sheds Light on Evolutionary Diversity in a Critically Endangered Cockatoo Species From Wallacea

Accurate identification of evolutionarily significant units of rare and threatened organisms provides a foundation for effective management and conservation. Up to seven subspecies of the critically endangered Yellow-crested Cockatoo (Cacatua sulphurea) have been described, four of which were commonly recognised pre-2014. In the absence of genotypic data, C. sulphurea subspecies delimitation has been based on morphology, behaviour and biogeography. To clarify genetic relationships and shed light on the diversification of this parrot radiation, whole genomes were sequenced for 16 museum specimens, covering the geographic range of the proposed seven subspecies as well as one C. galerita galerita. Combined with four museum-derived wild Cacatua sequences from NCBI, the results indicate there are three distinct C. sulphurea subspecies clusters centred in different biogeographic subregions of Wallacea (Timor; Sumba; as well as the Sulawesi Region and the main Lesser Sunda chain), separated by shallow genetic distances (da < 0.148%). The results raise questions about the recent species-level elevation of the phenotypically most distinct subspecies, C. s. citrinocristata, and about the origins of C. s. abbotti, the only subspecies west of Wallace's Line. Our analyses suggest C. s. abbotti is unlikely to be embedded within C. sulphurea, suggesting its origin on the remote Masalembu islands may be due to human translocation via historical trade routes. These genomic results inform the prioritisation and streamlining of conservation measures for the critically endangered C. sulphurea by identifying and delimiting likely conservation units.

Highlighting chromosomal rearrangements of five species of Galliformes (Domestic fowl, Common and Japanese quail, Barbary and Chukar partridge) and the Houbara bustard, an endangered Otidiformes: banding cytogenetic is a powerful tool

Birds are one of the most diverse groups among terrestrial vertebrates. They evolved from theropod dinosaurs, are closely related to the sauropsid group and separated from crocodiles about 240 million years ago. According to the IUCN, 12% of bird populations are threatened with potential extinction. Classical cytogenetics remains a powerful tool for comparing bird genomes and plays a crucial role in the preservation populations of endangered species. It thus makes it possible to detect chromosomal abnormalities responsible for early embryonic mortalities. Thus, in this work, we have provided new information on part of the evolutionary history by analysing high-resolution GTG-banded chromosomes to detect inter- and intrachromosomal rearrangements in six species. Indeed, the first eight autosomal pairs and the sex chromosomes of the domestic fowl Gallusgallusdomesticus Linnaeus, 1758 were compared with five species, four of which represent the order Galliformes (Common and Japanese quail, Gambras and Chukar partridge) and one Otidiformes species (Houbara bustard). Our findings suggest a high degree of conservation of the analysed ancestral chromosomes of the four Galliformes species, with the exception of (double, terminal, para and pericentric) inversions, deletion and the formation of neocentromeres (1, 2, 4, 7, 8, Z and W chromosomes). In addition to the detected rearrangements, reorganisation of the Houbara bustard chromosomes mainly included fusions and fissions involving both macro- and microchromosomes (especially on 2, 4 and Z chromosomes). We also found interchromosomal rearrangements involving shared microchromosomes (10, 11, 13, 14 and 19) between the two analysed avian orders. These rearrangements confirm that the structure of avian karyotypes will be more conserved at the interchromosomal but not at intrachromosomal scale. The appearance ofa small number of inter- and intrachromosomal rearrangements that occurred during evolution suggests a high degree of conservatism of genome organisation in these six species studied. A summary diagram of the rearrangements detected in this study is proposed to explain the chronology of the appearance of various evolutionary events starting from the ancestral karyotype.

How does the shape of the wing and hindlimb bones of aquatic birds relate to their locomotor abilities?

Aquatic birds represent diverse ecologies and locomotion types. Some became flightless or lost the ability for effective terrestrial locomotion, yet, certain species excel in water, on land, and in air, despite differing physical characteristics associated with each medium. In this exploratory study, we intend to quantitatively analyze the morphological variety of multiple limb bones of aquatic birds using 3D geometric morphometrics. Morphological variation is mainly driven by phylogeny, which also affects size and locomotion. However, the shape of the ulna, including the proportion and orientation of the epiphyses is influenced by size and aquatic propulsive techniques even when phylogeny is taken into consideration. Certain trends, possibly linked to functions, can be observed too in other bones, notably in cases where phylogenetic and functional signals are probably mixed when some taxa only englobe species with similar functional requirements: penguins exhibit the most distinctive wing bone morphologies, highly adapted to wing-propulsion; advanced foot-propellers exhibit femur morphology that reduces proximal mobility but supports stability; knee structures, like cnemial crests of varied sizes and orientations, are crucial for muscle attachments and efficient movement in water and on land; taxa relying on their feet in water but retaining terrestrial abilities share features enabling swimming and walking postures. Size-linked changes distinguish the wing bones of non-wing-propelled taxa. For hindlimbs, larger size relates to robust bones probably linked to terrestrial abilities, but robustness in femora can be connected to foot-propulsion. These results help us better understand birds' skeletal adaptation and can be useful inferring extinct species' ecology.

Dinosaur palaeoneurology: an evolving science

Our fascination with dinosaur brains and their capabilities essentially began with the first dinosaur discovery. The history of this study is a useful reflection of palaeoneurology as a whole and its relationship to a more inclusive evolutionary neuroscience. I argue that this relationship is imbued with high heuristic potential, but one whose realization requires overcoming certain constraints. These constraints include the need for a stable phylogenetic framework, methods for efficient and precise endocast construction, and fossil researchers who are steeped in a neuroscience perspective. The progress that has already been made in these areas sets the stage for a more mature palaeoneurology-not only one capable of being informed by neuroscience discoveries but one that drives such discoveries. I draw from work on the size, shape, behavioural correlates and developmental role of the dinosaur brain to outline current advances in dinosaur palaeoneurology. My examples largely are taken from theropods and centre on questions related to the origin of birds and their unique locomotory capabilities. The hope, however, is that these exemplify the potential for study in other dinosaur groups, and for utilizing the dinosaur-bird lineage as a parallel model on a par with mammals for studying encephalization.

Phylogenomic analysis of brachyuran crabs using transcriptome data reveals possible sources of conflicting phylogenetic relationships within the group

Despite extensive morphological and molecular studies, the phylogenetic interrelationships within the infraorder Brachyura and the phylogenetic positions of many taxa remain uncertain. Studies that used a limited number of molecular markers have often failed to provide sufficient resolution, and may be susceptible to stochastic errors and incomplete lineage sorting (ILS). Here we reconstructed the phylogenetic relationships within the Brachyura using transcriptome data of 56 brachyuran species, including 14 newly sequenced taxa. Five supermatrices were constructed in order to exclude different sources of systematic error. The results of the phylogenetic analyses indicate that Heterotremata is non-monophyletic, and that the two Old World primary freshwater crabs (Potamidae and Gecarcinucidae) and the Hymenosomatoidea form a clade that is sister to the Thoracotremata, and outside the Heterotremata. We also found that ILS is the main cause of the gene-tree discordance of these freshwater crabs. Divergence time estimations indicate that the Brachyura has an ancient origin, probably either in the Triassic or Jurassic, and that the majority of extant families and superfamilies first appeared during the Cretaceous, with a constant increase of diversity in Post-Cretaceous-Palaeogene times. The results support the hypothesis that the two Old World freshwater crab families included in this study (Potamidae and Gecarcinucidae) diverged from their marine ancestors around 120 Ma, in the Cretaceous. In addition, this work provides new insights that may aid in the reclassification of some of the more problematic brachyuran groups.

Warm Temperature is Associated With Reduced Body Mass and Diversification Rates While Increasing Extinction Risks in Cold-Adapted Seabirds

Anthropogenic rapid warming has caused decreases in richness and body mass of birds following the metabolic theory of ecology; yet, the pervasiveness of these shifts remains controversial among different taxa. Here, by combining phylogenetic methods and fossil data, we synthesized spatial patterns of richness and body mass for 328 seabird species belonging to two groups: Procellariimorphae (PM) and non-Procellariimorphae (NPM). We found that the relationship between body mass and richness, as well as diversification rate, exhibits distinct patterns in these two groups. Ancestral state reconstruction analyses indicate that smaller PM, as opposed to NPM seabirds, evolved in warmer waters from larger ancestors and exhibited a slower diversification rate. Different ancestral climatic origins explain the reduced influence of environmental factors on richness patterns among PM compared to NPM seabirds. Furthermore, whereas NPM seabirds in high latitudes face a high extinction risk, warmer sea temperatures positively correlate with a high extinction risk among PM seabirds. Our results indicate that PM seabirds, evolving from cold waters, have reduced body mass and diversification rate, making them more vulnerable to warmer temperature.

Terraces in species tree inference from gene trees

A terrace in a phylogenetic tree space is a region where all trees contain the same set of subtrees, due to certain patterns of missing data among the taxa sampled, resulting in an identical optimality score for a given data set. This was first investigated in the context of phylogenetic tree estimation from sequence alignments using maximum likelihood (ML) and maximum parsimony (MP). It was later extended to the species tree inference problem from a collection of gene trees, where a set of equally optimal species trees was referred to as a "pseudo" species tree terrace which does not consider the topological proximity of the trees in terms of the induced subtrees resulting from certain patterns of missing data. In this study, we mathematically characterize species tree terraces and investigate the mathematical properties and conditions that lead multiple species trees to induce/display an identical set of locus-specific subtrees owing to missing data. We report that species tree terraces are agnostic to gene tree heterogeneity. Therefore, we introduce and characterize a special type of gene tree topology-aware terrace which we call "peak terrace". Moreover, we empirically investigated various challenges and opportunities related to species tree terraces through extensive empirical studies using simulated and real biological data. We demonstrate the prevalence of species tree terraces and the resulting ambiguity created for tree search algorithms. Remarkably, our findings indicate that the identification of terraces could potentially lead to advances that enhance the accuracy of summary methods and provide reasonably accurate branch support.

Cellular, Molecular, and Genetic Mechanisms of Avian Beak Development and Evolution

Diverse research programs employing complementary strategies have been uncovering cellular, molecular, and genetic mechanisms essential to avian beak development and evolution. In reviewing these discoveries, I offer an interdisciplinary perspective on bird beaks that spans their derivation from jaws of dinosaurian reptiles, their anatomical and ecological diversification across major taxonomic groups, their common embryonic origins, their intrinsic patterning processes, and their structural integration. I describe how descriptive and experimental approaches, including gene expression and cell lineage analyses, tissue recombinations, surgical transplants, gain- and loss-of-function methods, geometric morphometrics, comparative genomics, and genome-wide association studies, have identified key constituent parts and putative genes regulating beak morphogenesis and evolution. I focus throughout on neural crest mesenchyme, which generates the beak skeleton and other components, and describe how these embryonic progenitor cells mediate species-specific pattern and link form and function as revealed by 20 years of research using chimeras between quail and duck embryos.

The Meaning and Measure of Concordance Factors in Phylogenomics

As phylogenomic datasets have grown in size, researchers have developed new ways to measure biological variation and to assess statistical support for specific branches. Larger datasets have more sites and loci and therefore less sampling variance. While we can more accurately measure the mean signal in these datasets, lower sampling variance is often reflected in uniformly high measures of branch support-such as the bootstrap and posterior probability-limiting their utility. Larger datasets have also revealed substantial biological variation in the topologies found across individual loci, such that the single species tree inferred by most phylogenetic methods represents a limited summary of the data for many purposes. In contrast to measures of statistical support, the degree of underlying topological variation among loci should be approximately constant regardless of the size of the dataset. "Concordance factors" (CFs) and similar statistics have therefore become increasingly important tools in phylogenetics. In this review, we explain why CFs should be thought of as descriptors of topological variation rather than as measures of statistical support, and argue that they provide important information about the predictive power of the species tree not contained in measures of support. We review a growing suite of statistics for measuring concordance, compare them in a common framework that reveals their interrelationships, and demonstrate how to calculate them using an example from birds. We also discuss how measures of topological variation might change in the future as we move beyond estimating a single "tree of life" toward estimating the myriad evolutionary histories underlying genomic variation.

Beyond genes-for-behaviour: The potential for genomics to resolve long-standing questions in avian brood parasitism

Behavioural ecology by definition of its founding 'Tinbergian framework' is an integrative field, however, it lags behind in incorporating genomic methods. 'Finding the gene/s for a behaviour' is still rarely feasible or cost-effective in the wild but as we show here, genomic data can be used to address broader questions. Here we use avian brood parasitism, a model system in behavioural ecology as a case study to highlight how behavioural ecologists could use the full potential of state-of-the-art genomic tools. Brood parasite-host interactions are one of the most easily observable and amenable natural laboratories of antagonistic coevolution, and as such have intrigued evolutionary biologists for decades. Using worked examples, we demonstrate how genomic data can be used to study the causes and mechanisms of (co)evolutionary adaptation and answer three key questions for the field: (i) Where and when should brood parasitism evolve?, (ii) When and how should hosts defend?, and (iii) Will coevolution persist with ecological change? In doing so, we discuss how behavioural and molecular ecologists can collaborate to integrate Tinbergen's questions and achieve the coherent science that he promoted to solve the mysteries of nature.

Insights into avian molecular cytogenetics-with reptilian comparisons

In last 100 years or so, much information has been accumulated on avian karyology, genetics, physiology, biochemistry and evolution. The chicken genome project generated genomic resources used in comparative studies, elucidating fundamental evolutionary processes, much of it funded by the economic importance of domestic fowl (which are also excellent model species in many areas). Studying karyotypes and whole genome sequences revealed population processes, evolutionary biology, and genome function, uncovering the role of repetitive sequences, transposable elements and gene family expansion. Knowledge of the function of many genes and non-expressed or identified regulatory components is however still lacking. Birds (Aves) are diverse, have striking adaptations for flight, migration and survival and inhabit all continents most islands. They also have a unique karyotype with ~ 10 macrochromosomes and ~ 30 microchromosomes that are smaller than other reptiles. Classified into Palaeognathae and Neognathae they are evolutionarily close, and a subset of reptiles. Here we overview avian molecular cytogenetics with reptilian comparisons, shedding light on their karyotypes and genome structure features. We consider avian evolution, then avian (followed by reptilian) karyotypes and genomic features. We consider synteny disruptions, centromere repositioning, and repetitive elements before turning to comparative avian and reptilian genomics. In this context, we review comparative cytogenetics and genome mapping in birds as well as Z- and W-chromosomes and sex determination. Finally, we give examples of pivotal research areas in avian and reptilian cytogenomics, particularly physical mapping and map integration of sex chromosomal genes, comparative genomics of chicken, turkey and zebra finch, California condor cytogenomics as well as some peculiar cytogenetic and evolutionary examples. We conclude that comparative molecular studies and improving resources continually contribute to new approaches in population biology, developmental biology, physiology, disease ecology, systematics, evolution and phylogenetic systematics orientation. This also produces genetic mapping information for chromosomes active in rearrangements during the course of evolution. Further insights into mutation, selection and adaptation of vertebrate genomes will benefit from these studies including physical and online resources for the further elaboration of comparative genomics approaches for many fundamental biological questions.

Whole Genomes Reveal Evolutionary Relationships and Mechanisms Underlying Gene-Tree Discordance in Neodiprion Sawflies

Rapidly evolving taxa are excellent models for understanding the mechanisms that give rise to biodiversity. However, developing an accurate historical framework for comparative analysis of such lineages remains a challenge due to ubiquitous incomplete lineage sorting (ILS) and introgression. Here, we use a whole-genome alignment, multiple locus-sampling strategies, and summary-tree and single nucleotide polymorphism-based species-tree methods to infer a species tree for eastern North American Neodiprion species, a clade of pine-feeding sawflies (Order: Hymenopteran; Family: Diprionidae). We recovered a well-supported species tree that-except for three uncertain relationships-was robust to different strategies for analyzing whole-genome data. Nevertheless, underlying gene-tree discordance was high. To understand this genealogical variation, we used multiple linear regression to model site concordance factors estimated in 50-kb windows as a function of several genomic predictor variables. We found that site concordance factors tended to be higher in regions of the genome with more parsimony-informative sites, fewer singletons, less missing data, lower GC content, more genes, lower recombination rates, and lower D-statistics (less introgression). Together, these results suggest that ILS, introgression, and genotyping error all shape the genomic landscape of gene-tree discordance in Neodiprion. More generally, our findings demonstrate how combining phylogenomic analysis with knowledge of local genomic features can reveal mechanisms that produce topological heterogeneity across genomes.

Phylogenomic Discordance is Driven by Wide-Spread Introgression and Incomplete Lineage Sorting During Rapid Species Diversification Within Rattlesnakes (Viperidae: Crotalus and Sistrurus)

-Phylogenomics allows us to uncover the historical signal of evolutionary processes through time and estimate phylogenetic networks accounting for these signals. Insight from genome-wide data further allows us to pinpoint the contributions to phylogenetic signal from hybridization, introgression, and ancestral polymorphism across the genome. Here, we focus on how these processes have contributed to phylogenetic discordance among rattlesnakes (genera Crotalus and Sistrurus), a group for which there are numerous conflicting phylogenetic hypotheses based on a diverse array of molecular datasets and analytical methods. We address the instability of the rattlesnake phylogeny using genomic data generated from transcriptomes sampled from nearly all known species. These genomic data, analyzed with coalescent and network-based approaches, reveal numerous instances of rapid speciation where individual gene trees conflict with the species tree. Moreover, the evolutionary history of rattlesnakes is dominated by incomplete speciation and frequent hybridization, both of which have likely influenced past interpretations of phylogeny. We present a new framework in which the evolutionary relationships of this group can only be understood in light of genome-wide data and network-based analytical methods. Our data suggest that network radiations, like those seen within the rattlesnakes, can only be understood in a phylogenomic context, necessitating similar approaches in our attempts to understand evolutionary history in other rapidly radiating species.

Characterization of genome-wide phylogenetic conflict uncovers evolutionary modes of carnivorous fungi

Mass extinction has often paved the way for rapid evolutionary radiation, resulting in the emergence of diverse taxa within specific lineages. The emergence and diversification of carnivorous nematode-trapping fungi (NTF) in Ascomycota have been linked to the Permian-Triassic (PT) extinction, but the processes underlying NTF radiation remain unclear. We conducted phylogenomic analyses using 23 genomes that represent three NTF lineages, each employing distinct nematode traps-mechanical traps (Drechslerella spp.), three-dimensional (3D) adhesive traps (Arthrobotrys spp.), and two-dimensional (2D) adhesive traps (Dactylellina spp.), and the genome of one non-NTF species as the outgroup. These analyses revealed multiple mechanisms that likely contributed to the tempo of the NTF evolution and rapid radiation. The species tree of NTFs based on 2,944 single-copy orthologous genes suggested that Drechslerella emerged earlier than Arthrobotrys and Dactylellina. Extensive genome-wide phylogenetic discordance was observed, mainly due to incomplete lineage sorting (ILS) between lineages. Two modes of non-vertical evolution (introgression and horizontal gene transfer) also contributed to phylogenetic discordance. The ILS genes that are associated with hyphal growth and trap morphogenesis (e.g., those associated with the cell membrane system and polarized cell division) exhibited signs of positive selection.IMPORTANCEBy conducting a comprehensive phylogenomic analysis of 23 genomes across three NTF lineages, the research reveals how diverse evolutionary mechanisms, including ILS and non-vertical evolution (introgression and horizontal gene transfer), contribute to the swift diversification of NTFs. These findings highlight the complex evolutionary dynamics that drive the rapid radiation of NTFs, providing valuable insights into the processes underlying their diversity and adaptation.

Palaeognaths Reveal Evolutionary Ancestry of the Avian Major Histocompatibility Complex Class II

The multigene family of the major histocompatibility complex (MHC) codes for the key antigen-presenting molecules of the vertebrate immune system. In birds, duplicated MHC class II (MHC-II) genes are highly homogenized by concerted evolution, and thus, identification of their orthologous relationships across long evolutionary timescales remains challenging. Relatively low evolutionary rate of avian MHC class IIA genes has been expected to provide a promising avenue to allow such inferences, but availability of MHC-IIA sequences in nonmodel bird species has been limited until recently. Here, taking advantage from accumulating genomic resources, we identified and analyzed MHC-IIA sequences from the most basal lineage of extant birds (Palaeognathae). Conserved region of the MHC-IIA membrane-proximal domain was used to search for orthologous relationships between palaeognath birds and nonavian reptiles. First, analyses of palaeognath sequences revealed the presence of a separate MHC-IIA gene lineage (DAA3) in kiwis, which did not cluster with previously described avian MHC-IIA lineages (DAA1 and DAA2). Next, phylogenetic reconstruction showed that kiwi DAA3 sequences form a single well-supported cluster with turtle MHC-IIA. High similarity of these sequences most likely reflects their remarkable evolutionary conservation and retention of ancient orthologous relationships, which can be traced back to basal archosauromorphs ca. 250 million years ago. Our analyses offer novel insights into macroevolutionary history of the MHC and reinforce the view that rapid accumulation of high-quality genome assemblies across divergent nonmodel species can substantially advance our understanding of gene evolution.

SEGUL: Ultrafast, memory-efficient and mobile-friendly software for manipulating and summarizing phylogenomic datasets

Phylogenetic studies now routinely require manipulating and summarizing thousands of data files. For most of these tasks, currently available software requires considerable computing resources and substantial knowledge of command-line applications. We develop an ultrafast and memory-efficient software, SEGUL, that performs common phylogenomic dataset manipulations and calculates statistics summarizing essential data features. Our software is available as standalone command-line interface (CLI) and graphical user interface (GUI) applications, and as a library for Rust, R and Python, with possible support of other languages. The CLI and library versions run native on Windows, Linux and macOS, including Apple ARM Macs. The GUI version extends support to include mobile iOS, iPadOS and Android operating systems. SEGUL leverages the high performance of the Rust programming language to offer fast execution times and low memory footprints regardless of dataset size and platform choice. The inclusion of a GUI minimizes bioinformatics barriers to phylogenomics while SEGUL's efficiency reduces economic barriers by allowing analysis on inexpensive hardware. Our support for mobile operating systems further enables teaching phylogenomics where access to computing power is limited.

Corvisyringophilus, a New Genus in the Family Syringophilidae (Acariformes: Prostigmata) and Its Phylogenetic Position among Primitive Genera

Syringophilidae is one of the most species-rich families in the superfamily Cheyletoidea, comprising approximately 420 species across 62 genera and two subfamilies. In this paper, we propose a new genus, Corvisyringophilus, and a new species, C. krummi gen. n. et sp. n., found in the wing covert quills of the Common Raven, Corvus corax Linnaeus, in Iceland. Corvisyringophilus is placed among the primitive genera of syringophilid mites, which possess the full complement of idiosomal and leg setae. Phylogenetic analysis based on morphological characters suggests that this genus forms a sister clade to Blaszakia Skoracki & Sikora, 2008, and Charadriphilus Bochkov & Mironov, 1998, which inhabit birds of the orders Musophagiformes and Charadriiformes, respectively. The study proposes that the current distribution patterns of quill mites, based on their morphological characteristics, may result from multiple host switching followed by co-speciation events, highlighting the complex evolutionary dynamics within this family.

Dating in the Dark: Elevated Substitution Rates in Cave Cockroaches (Blattodea: Nocticolidae) Have Negative Impacts on Molecular Date Estimates

Rates of nucleotide substitution vary substantially across the Tree of Life, with potentially confounding effects on phylogenetic and evolutionary analyses. A large acceleration in mitochondrial substitution rate occurs in the cockroach family Nocticolidae, which predominantly inhabit subterranean environments. To evaluate the impacts of this among-lineage rate heterogeneity on estimates of phylogenetic relationships and evolutionary timescales, we analyzed nuclear ultraconserved elements (UCEs) and mitochondrial genomes from nocticolids and other cockroaches. Substitution rates were substantially elevated in nocticolid lineages compared with other cockroaches, especially in mitochondrial protein-coding genes. This disparity in evolutionary rates is likely to have led to different evolutionary relationships being supported by phylogenetic analyses of mitochondrial genomes and UCE loci. Furthermore, Bayesian dating analyses using relaxed-clock models inferred much deeper divergence times compared with a flexible local clock. Our phylogenetic analysis of UCEs, which is the first genome-scale study to include all 13 major cockroach families, unites Corydiidae and Nocticolidae and places Anaplectidae as the sister lineage to the rest of Blattoidea. We uncover an extraordinary level of genetic divergence in Nocticolidae, including two highly distinct clades that separated ~115 million years ago despite both containing representatives of the genus Nocticola. The results of our study highlight the potential impacts of high among-lineage rate variation on estimates of phylogenetic relationships and evolutionary timescales.

Insertions and Deletions: Computational Methods, Evolutionary Dynamics, and Biological Applications

Insertions and deletions constitute the second most important source of natural genomic variation. Insertions and deletions make up to 25% of genomic variants in humans and are involved in complex evolutionary processes including genomic rearrangements, adaptation, and speciation. Recent advances in long-read sequencing technologies allow detailed inference of insertions and deletion variation in species and populations. Yet, despite their importance, evolutionary studies have traditionally ignored or mishandled insertions and deletions due to a lack of comprehensive methodologies and statistical models of insertions and deletion dynamics. Here, we discuss methods for describing insertions and deletion variation and modeling insertions and deletions over evolutionary time. We provide practical advice for tackling insertions and deletions in genomic sequences and illustrate our discussion with examples of insertions and deletion-induced effects in human and other natural populations and their contribution to evolutionary processes. We outline promising directions for future developments in statistical methodologies that would allow researchers to analyze insertions and deletion variation and their effects in large genomic data sets and to incorporate insertions and deletions in evolutionary inference.