The phylogenomic forest of bird trees contains a hard polytomy at the root of Neoaves

Genome and life-history evolution link bird diversification to the end-Cretaceous mass extinction

Complex patterns of genome evolution associated with the end-Cretaceous [Cretaceous-Paleogene (K-Pg)] mass extinction limit our understanding of the early evolutionary history of modern birds. Here, we analyzed patterns of avian molecular evolution and identified distinct macroevolutionary regimes across exons, introns, untranslated regions, and mitochondrial genomes. Bird clades originating near the K-Pg boundary exhibited numerous shifts in the mode of molecular evolution, suggesting a burst of genomic heterogeneity at this point in Earth's history. These inferred shifts in substitution patterns were closely related to evolutionary shifts in developmental mode, adult body mass, and patterns of metabolic scaling. Our results suggest that the end-Cretaceous mass extinction triggered integrated patterns of evolution across avian genomes, physiology, and life history near the dawn of the modern bird radiation.

Dragons in the tropics - Phylogeography and speciation in Diporiphora lizards and common geographic breaks in co-distributed taxa

Co-distributed taxa can respond both similarly or differently to the same climatic and geological events, resulting in a range of phylogeographic patterns across the region. Using a nested approach on a taxonomically diverse yet morphologically conservative group of agamid lizards, we first aimed to evaluate more precisely the extent of phylogeographic structuring within the genus. Then, focusing on four lineages within the more widespread species, we assessed the impact of biogeographic barriers on phylogeographic structuring and demographic history of species, comparing to patterns previously observed in co-distributed taxa. These species occur in the Australian Monsoonal Tropics, a vast tropical savanna system with high richness and endemism associated with environmental heterogeneity and past climate fluctuations. The employment of genomic data helped to determine the relationships between specific taxa that were previously difficult to place. We found a local influence of biogeographic and climatic breaks on population dynamics, analogous to other species. We detected high levels of population structure in the West Kimberley and Arnhem Plateau, which are already known for high endemism. However, we also highlighted unique lineages in areas that have been overlooked until recently, in the South Kimberley and West Top End. Climatic and geographical features in the Arnhem Plateau act as a soft barrier between populations in the east and west regions of the Top End. These observations reflect patterns observed for other vertebrates across this rich biome, indicating how climatic variation, species' ecology, and landscape features interact to shape regional diversity and endemism.

Hybridization in birds-of-paradise: Widespread ancestral gene flow despite strong sexual selection in a lek-mating system

Sexual selection can directly contribute to reproductive isolation and is an important mechanism that can lead to speciation. Lek-mating is one of the most extreme forms of sexual selection, but surprisingly does not seem to preclude occasional hybridization in nature. However, hybridization among lekking species may still be trivial if selection against offspring with intermediate phenotypes prohibits introgression. Here we investigate this further by sequencing the genomes of nearly all bird-of-paradise (Paradisaeidae) species and 10 museum specimens of putative hybrid origin. We find that intergeneric hybridization indeed still takes place despite extreme differentiation in form, plumage, and behavior. In parallel, the genomes of contemporary species contain widespread signatures of past introgression, demonstrating that hybridization has repeatedly resulted in shared genetic variation despite strong sexual isolation. Our study raises important questions about extrinsic factors that modulate hybridization probability and the evolutionary consequences of introgressive hybridization between lekking species.

Evolving dispersal ability causes rapid adaptive radiation

Despite a long history of research since Darwin, the mechanism underlying rapid adaptive radiation remains poorly understood. All theories constructed to date require special assumptions, so none can comprehensively explain actual cases found in wide-ranging taxonomic groups. Here, we propose a simple theoretical solution to this problem. Namely, we extend the classical archipelago model of adaptive radiation into a more realistic model by adding one assumption, namely, the evolvability of dispersal ability, which is well supported empirically. Our individual-based simulations with evolvable dispersal ability showed that environmental heterogeneity among islands (or island-like habitats) led to an evolutionary decrease in dispersal ability. However, when islands are rather evenly distributed, as is often the case in actual archipelagos where adaptive radiation has been reported, the decline in dispersal ability that began in some island populations was quickly halted by the continuous influx of immigrants from other islands. The process of reduction in dispersal ability in these island populations was resumed almost synchronously when the dispersal ability began to decrease on the final island, which had maintained high dispersal ability and continued to release migrants for the longest duration. Then, a rapid loss of dispersal ability followed in all island populations. In short, the frequent simultaneous evolution of multiple allopatric incipient species was an inevitable consequence of the properties of ordinary archipelagos in our simulations. This study strongly suggests that the seemingly complex process of rapid radiation is driven by a simple mechanism of evolutionary reduction in dispersal ability.

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

Inference of species trees plays a crucial role in advancing our understanding of evolutionary relationships and has immense significance for diverse biological and medical applications. Extensive genome sequencing efforts are currently in progress across a broad spectrum of life forms, holding the potential to unravel the intricate branching patterns within the tree of life. However, estimating species trees starting from raw genome sequences is quite challenging, and the current cutting-edge methodologies require a series of error-prone steps that are neither entirely automated nor standardized. In this paper, we present ROADIES, a novel pipeline for species tree inference from raw genome assemblies that is fully automated, easy to use, scalable, free from reference bias, and provides flexibility to adjust the tradeoff between accuracy and runtime. The ROADIES pipeline eliminates the need to align whole genomes, choose a single reference species, or pre-select loci such as functional genes found using cumbersome annotation steps. Moreover, it leverages recent advances in phylogenetic inference to allow multi-copy genes, eliminating the need to detect orthology. Using the genomic datasets released from large-scale sequencing consortia across three diverse life forms (placental mammals, pomace flies, and birds), we show that ROADIES infers species trees that are comparable in quality with the state-of-the-art approaches but in a fraction of the time. By incorporating optimal approaches and automating all steps from assembled genomes to species and gene trees, ROADIES is poised to improve the accuracy, scalability, and reproducibility of phylogenomic analyses.

Complexity of avian evolution revealed by family-level genomes

Despite tremendous efforts in the past decades, relationships among main avian lineages remain heavily debated without a clear resolution. Discrepancies have been attributed to diversity of species sampled, phylogenetic method and the choice of genomic regions1-3. Here we address these issues by analysing the genomes of 363 bird species4 (218 taxonomic families, 92% of total). Using intergenic regions and coalescent methods, we present a well-supported tree but also a marked degree of discordance. The tree confirms that Neoaves experienced rapid radiation at or near the Cretaceous-Palaeogene boundary. Sufficient loci rather than extensive taxon sampling were more effective in resolving difficult nodes. Remaining recalcitrant nodes involve species that are a challenge to model due to either extreme DNA composition, variable substitution rates, incomplete lineage sorting or complex evolutionary events such as ancient hybridization. Assessment of the effects of different genomic partitions showed high heterogeneity across the genome. We discovered sharp increases in effective population size, substitution rates and relative brain size following the Cretaceous-Palaeogene extinction event, supporting the hypothesis that emerging ecological opportunities catalysed the diversification of modern birds. The resulting phylogenetic estimate offers fresh insights into the rapid radiation of modern birds and provides a taxon-rich backbone tree for future comparative studies.

Gene flow and an anomaly zone complicate phylogenomic inference in a rapidly radiated avian family (Prunellidae)

BACKGROUND

Resolving the phylogeny of rapidly radiating lineages presents a challenge when building the Tree of Life. An Old World avian family Prunellidae (Accentors) comprises twelve species that rapidly diversified at the Pliocene-Pleistocene boundary.

RESULTS

Here we investigate the phylogenetic relationships of all species of Prunellidae using a chromosome-level de novo assembly of Prunella strophiata and 36 high-coverage resequenced genomes. We use homologous alignments of thousands of exonic and intronic loci to build the coalescent and concatenated phylogenies and recover four different species trees. Topology tests show a large degree of gene tree-species tree discordance but only 40-54% of intronic gene trees and 36-75% of exonic genic trees can be explained by incomplete lineage sorting and gene tree estimation errors. Estimated branch lengths for three successive internal branches in the inferred species trees suggest the existence of an empirical anomaly zone. The most common topology recovered for species in this anomaly zone was not similar to any coalescent or concatenated inference phylogenies, suggesting presence of anomalous gene trees. However, this interpretation is complicated by the presence of gene flow because extensive introgression was detected among these species. When exploring tree topology distributions, introgression, and regional variation in recombination rate, we find that many autosomal regions contain signatures of introgression and thus may mislead phylogenetic inference. Conversely, the phylogenetic signal is concentrated to regions with low-recombination rate, such as the Z chromosome, which are also more resistant to interspecific introgression.

CONCLUSIONS

Collectively, our results suggest that phylogenomic inference should consider the underlying genomic architecture to maximize the consistency of phylogenomic signal.

The Chalcidoidea bush of life: evolutionary history of a massive radiation of minute wasps

Chalcidoidea are mostly parasitoid wasps that include as many as 500 000 estimated species. Capturing phylogenetic signal from such a massive radiation can be daunting. Chalcidoidea is an excellent example of a hyperdiverse group that has remained recalcitrant to phylogenetic resolution. We combined 1007 exons obtained with Anchored Hybrid Enrichment with 1048 ultra-conserved elements (UCEs) for 433 taxa including all extant families, >95% of all subfamilies, and 356 genera chosen to represent the vast diversity of the superfamily. Going back and forth between the molecular results and our collective knowledge of morphology and biology, we detected bias in the analyses that was driven by the saturation of nucleotide data. Our final results are based on a concatenated analysis of the least saturated exons and UCE datasets (2054 loci, 284 106 sites). Our analyses support an expected sister relationship with Mymarommatoidea. Seven previously recognized families were not monophyletic, so support for a new classification is discussed. Natural history in some cases would appear to be more informative than morphology, as illustrated by the elucidation of a clade of plant gall associates and a clade of taxa with planidial first-instar larvae. The phylogeny suggests a transition from smaller soft-bodied wasps to larger and more heavily sclerotized wasps, with egg parasitism as potentially ancestral for the entire superfamily. Deep divergences in Chalcidoidea coincide with an increase in insect families in the fossil record, and an early shift to phytophagy corresponds with the beginning of the "Angiosperm Terrestrial Revolution". Our dating analyses suggest a middle Jurassic origin of 174 Ma (167.3-180.5 Ma) and a crown age of 162.2 Ma (153.9-169.8 Ma) for Chalcidoidea. During the Cretaceous, Chalcidoidea may have undergone a rapid radiation in southern Gondwana with subsequent dispersals to the Northern Hemisphere. This scenario is discussed with regard to knowledge about the host taxa of chalcid wasps, their fossil record and Earth's palaeogeographic history.

A Guide to Phylogenomic Inference

Phylogenomics aims at reconstructing the evolutionary histories of organisms taking into account whole genomes or large fractions of genomes. Phylogenomics has significant applications in fields such as evolutionary biology, systematics, comparative genomics, and conservation genetics, providing valuable insights into the origins and relationships of species and contributing to our understanding of biological diversity and evolution. This chapter surveys phylogenetic concepts and methods aimed at both gene tree and species tree reconstruction while also addressing common pitfalls, providing references to relevant computer programs. A practical phylogenomic analysis example including bacterial genomes is presented at the end of the chapter.

A comprehensive phylogenomic study unveils evolutionary patterns and challenges in the mitochondrial genomes of Carcharhiniformes: A focus on Triakidae

The complex evolutionary patterns in the mitochondrial genome (mitogenome) of the most species-rich shark order, the Carcharhiniformes (ground sharks) has led to challenges in the phylogenomic reconstruction of the families and genera belonging to the order, particularly the family Triakidae (houndsharks). The current state of Triakidae phylogeny remains controversial, with arguments for both monophyly and paraphyly within the family. We hypothesize that this variability is triggered by the selection of different a priori partitioning schemes to account for site and gene heterogeneity within the mitogenome. Here we used an extensive statistical framework to select the a priori partitioning scheme for inference of the mitochondrial phylogenomic relationships within Carcharhiniformes, tested site heterogeneous CAT + GTR + G4 models and incorporated the multi-species coalescent model (MSCM) into our analyses to account for the influence of gene tree discordance on species tree inference. We included five newly assembled houndshark mitogenomes to increase resolution of Triakidae. During the assembly procedure, we uncovered a 714 bp-duplication in the mitogenome of Galeorhinus galeus. Phylogenetic reconstruction confirmed monophyly within Triakidae and the existence of two distinct clades of the expanded Mustelus genus. The latter alludes to potential evolutionary reversal of reproductive mode from placental to aplacental, suggesting that reproductive mode has played a role in the trajectory of adaptive divergence. These new sequences have the potential to contribute to population genomic investigations, species phylogeography delineation, environmental DNA metabarcoding databases and, ultimately, improved conservation strategies for these ecologically and economically important species.

The genome of the pygmy right whale illuminates the evolution of rorquals

BACKGROUND

Baleen whales are a clade of gigantic and highly specialized marine mammals. Their genomes have been used to investigate their complex evolutionary history and to decipher the molecular mechanisms that allowed them to reach these dimensions. However, many unanswered questions remain, especially about the early radiation of rorquals and how cancer resistance interplays with their huge number of cells. The pygmy right whale is the smallest and most elusive among the baleen whales. It reaches only a fraction of the body length compared to its relatives and it is the only living member of an otherwise extinct family. This placement makes the pygmy right whale genome an interesting target to update the complex phylogenetic past of baleen whales, because it splits up an otherwise long branch that leads to the radiation of rorquals. Apart from that, genomic data of this species might help to investigate cancer resistance in large whales, since these mechanisms are not as important for the pygmy right whale as in other giant rorquals and right whales.

RESULTS

Here, we present a first de novo genome of the species and test its potential in phylogenomics and cancer research. To do so, we constructed a multi-species coalescent tree from fragments of a whole-genome alignment and quantified the amount of introgression in the early evolution of rorquals. Furthermore, a genome-wide comparison of selection rates between large and small-bodied baleen whales revealed a small set of conserved candidate genes with potential connections to cancer resistance.

CONCLUSIONS

Our results suggest that the evolution of rorquals is best described as a hard polytomy with a rapid radiation and high levels of introgression. The lack of shared positive selected genes between different large-bodied whale species supports a previously proposed convergent evolution of gigantism and hence cancer resistance in baleen whales.

A Phylogenomic Assessment of Processes Underpinning Convergent Evolution in Open-Habitat Chats

Insights into the processes underpinning convergent evolution advance our understanding of the contributions of ancestral, introgressed, and novel genetic variation to phenotypic evolution. Phylogenomic analyses characterizing genome-wide gene tree heterogeneity can provide first clues about the extent of ILS and of introgression and thereby into the potential of these processes or (in their absence) the need to invoke novel mutations to underpin convergent evolution. Here, we were interested in understanding the processes involved in convergent evolution in open-habitat chats (wheatears of the genus Oenanthe and their relatives). To this end, based on whole-genome resequencing data from 50 taxa of 44 species, we established the species tree, characterized gene tree heterogeneity, and investigated the footprints of ILS and introgression within the latter. The species tree corroborates the pattern of abundant convergent evolution, especially in wheatears. The high levels of gene tree heterogeneity in wheatears are explained by ILS alone only for 30% of internal branches. For multiple branches with high gene tree heterogeneity, D-statistics and phylogenetic networks identified footprints of introgression. Finally, long branches without extensive ILS between clades sporting similar phenotypes provide suggestive evidence for the role of novel mutations in the evolution of these phenotypes. Together, our results suggest that convergent evolution in open-habitat chats involved diverse processes and highlight that phenotypic diversification is often complex and best depicted as a network of interacting lineages.

Ancient proteins resolve controversy over the identity of Genyornis eggshell

The realization that ancient biomolecules are preserved in "fossil" samples has revolutionized archaeological science. Protein sequences survive longer than DNA, but their phylogenetic resolution is inferior; therefore, careful assessment of the research questions is required. Here, we show the potential of ancient proteins preserved in Pleistocene eggshell in addressing a longstanding controversy in human and animal evolution: the identity of the extinct bird that laid large eggs which were exploited by Australia's indigenous people. The eggs had been originally attributed to the iconic extinct flightless bird Genyornis newtoni (†Dromornithidae, Galloanseres) and were subsequently dated to before 50 ± 5 ka by Miller et al. [Nat. Commun. 7, 10496 (2016)]. This was taken to represent the likely extinction date for this endemic megafaunal species and thus implied a role of humans in its demise. A contrasting hypothesis, according to which the eggs were laid by a large mound-builder megapode (Megapodiidae, Galliformes), would therefore acquit humans of their responsibility in the extinction of Genyornis. Ancient protein sequences were reconstructed and used to assess the evolutionary proximity of the undetermined eggshell to extant birds, rejecting the megapode hypothesis. Authentic ancient DNA could not be confirmed from these highly degraded samples, but morphometric data also support the attribution of the eggshell to Genyornis. When used in triangulation to address well-defined hypotheses, paleoproteomics is a powerful tool for reconstructing the evolutionary history in ancient samples. In addition to the clarification of phylogenetic placement, these data provide a more nuanced understanding of the modes of interactions between humans and their environment.

Comprehensive taxon sampling and vetted fossils help clarify the time tree of shorebirds (Aves, Charadriiformes)

Shorebirds (Charadriiformes) are a globally distributed clade of modern birds and, due to their ecological and morphological disparity, a frequent subject of comparative studies. While molecular phylogenies have been key to establishing the suprafamilial backbone of the charadriiform tree, a number of relationships at both deep and shallow taxonomic levels remain poorly resolved. The timescale of shorebird evolution also remains uncertain as a result of extensive disagreements among the published divergence dating studies, stemming largely from different choices of fossil calibrations. Here, we present the most comprehensive non-supertree phylogeny of shorebirds to date, based on a total-evidence dataset comprising 353 ingroup taxa (90% of all extant or recently extinct species), 27 loci (15 mitochondrial and 12 nuclear), and 69 morphological characters. We further clarify the timeline of charadriiform evolution by time-scaling this phylogeny using a set of 14 up-to-date and thoroughly vetted fossil calibrations. In addition, we assemble a taxonomically restricted 100-locus dataset specifically designed to resolve outstanding problems in higher-level charadriiform phylogeny. In terms of tree topology, our results are largely congruent with previous studies but indicate that some of the conflicts among earlier analyses reflect a genuine signal of pervasive gene tree discordance. Monophyly of the plovers (Charadriidae), the position of the ibisbill (Ibidorhyncha), and the relationships among the five subfamilies of the gulls (Laridae) could not be resolved even with greatly increased locus and taxon sampling. Moreover, several localized regions of uncertainty persist in shallower parts of the tree, including the interrelationships of the true auks (Alcinae) and anarhynchine plovers. Our node-dating and macroevolutionary rate analyses find support for a Paleocene origin of crown-group shorebirds, as well as exceptionally rapid recent radiations of Old World oystercatchers (Haematopodidae) and select genera of gulls. Our study underscores the challenges involved in estimating a comprehensively sampled and carefully calibrated time tree for a diverse avian clade, and highlights areas in need of further research.

The Application of Quantitative Metabolomics for the Taxonomic Differentiation of Birds

In the current pilot study, we propose the use of quantitative metabolomics to reconstruct the phylogeny of vertebrates, namely birds. We determined the concentrations of the 67 most abundant metabolites in the eye lenses of the following 14 species from 6 orders of the class Aves (Birds): the Black kite (Milvus migrans), Eurasian magpie (Pica pica), Northern raven (Corvus corax), Eurasian coot (Fulica atra), Godlewski's bunting (Emberiza godlewskii), Great crested grebe (Podiceps cristatus), Great tit (Parus major), Hawfinch (Coccothraustes coccothraustes), Hooded crow (Corvus cornix), House sparrow (Passer domesticus), Rock dove (Columba livia), Rook (Corvus frugilegus), Short-eared owl (Asio flammeus) and Ural owl (Strix uralensis). Further analysis shows that the statistical approaches generally used in metabolomics can be applied for differentiation between species, and the most fruitful results were obtained with hierarchical clustering analysis (HCA). We observed the grouping of conspecific samples independently of the sampling place and date. The HCA tree structure supports the key role of genomics in the formation of the lens metabolome, but it also indicates the influence of the species lifestyle. A combination of genomics-based and metabolomics-based phylogeny could potentially resolve arising issues and yield a more reliable tree of life.

Phylogenomic Coalescent Analyses of Avian Retroelements Infer Zero-Length Branches at the Base of Neoaves, Emergent Support for Controversial Clades, and Ancient Introgressive Hybridization in Afroaves

Retroelement insertions (RIs) are low-homoplasy characters that are ideal data for addressing deep evolutionary radiations, where gene tree reconstruction errors can severely hinder phylogenetic inference with DNA and protein sequence data. Phylogenomic studies of Neoaves, a large clade of birds (>9000 species) that first diversified near the Cretaceous−Paleogene boundary, have yielded an array of robustly supported, contradictory relationships among deep lineages. Here, we reanalyzed a large RI matrix for birds using recently proposed quartet-based coalescent methods that enable inference of large species trees including branch lengths in coalescent units, clade-support, statistical tests for gene flow, and combined analysis with DNA-sequence-based gene trees. Genome-scale coalescent analyses revealed extremely short branches at the base of Neoaves, meager branch support, and limited congruence with previous work at the most challenging nodes. Despite widespread topological conflicts with DNA-sequence-based trees, combined analyses of RIs with thousands of gene trees show emergent support for multiple higher-level clades (Columbea, Passerea, Columbimorphae, Otidimorphae, Phaethoquornithes). RIs express asymmetrical support for deep relationships within the subclade Afroaves that hints at ancient gene flow involving the owl lineage (Strigiformes). Because DNA-sequence data are challenged by gene tree-reconstruction error, analysis of RIs represents one approach for improving gene tree-based methods when divergences are deep, internodes are short, terminal branches are long, and introgressive hybridization further confounds species−tree inference.

Categorical edge-based analyses of phylogenomic data reveal conflicting signals for difficult relationships in the avian tree

Phylogenetic analyses fail to yield a satisfactory resolution of some relationships in the tree of life even with genome-scale datasets, so the failure is unlikely to reflect limitations in the amount of data. Gene tree conflicts are particularly notable in studies focused on these contentious nodes, and taxon sampling, different analytical methods, and/or data type effects can further confound analyses. Although many efforts have been made to incorporate biological conflicts, few studies have curated individual genes for their efficiency in phylogenomic studies. Here, we conduct an edge-based analysis of Neoavian evolution, examining the phylogenetic efficacy of two recent phylogenomic bird datasets and three datatypes (ultraconserved elements [UCEs], introns, and coding regions). We assess the potential causes for biases in signal-resolution for three difficult nodes: the earliest divergence of Neoaves, the position of the enigmatic Hoatzin (Opisthocomus hoazin), and the position of owls (Strigiformes). We observed extensive conflict among genes for all data types and datasets even after meticulous curation. Edge-based analyses (EBA) increased congruence and provided information about the impact of data type, GC content variation (GC_CV), and outlier genes on each of nodes we examined. First, outlier gene signals appeared to drive different patterns of support for the relationships among the earliest diverging Neoaves. Second, the placement of Hoatzin was highly variable, although our EBA did reveal a previously unappreciated data type effect with an impact on its position. It also revealed that the resolution with the most support here was Hoatzin + shorebirds. Finally, GC_CV, rather than data type (i.e., coding vs non-coding) per se, was correlated with a signal that supports monophyly of owls + Accipitriformes (hawks, eagles, and vultures). Eliminating high GC_CV loci increased the signal for owls + mousebirds. Categorical EBA was able to reveal the nature of each edge and provide a way to highlight especially problematic branches that warrant a further examination. The current study increases our understanding about the contentious parts of the avian tree, which show even greater conflicts than appreciated previously.

Speciation and population divergence in a mutualistic seed dispersing bird

Bird-mediated seed dispersal is crucial for the regeneration and viability of ecosystems, often resulting in complex mutualistic species networks. Yet, how this mutualism drives the evolution of seed dispersing birds is still poorly understood. In the present study we combine whole genome re-sequencing analyses and morphometric data to assess the evolutionary processes that shaped the diversification of the Eurasian nutcracker (Nucifraga), a seed disperser known for its mutualism with pines (Pinus). Our results show that the divergence and phylogeographic patterns of nutcrackers resemble those of other non-mutualistic passerine birds and suggest that their early diversification was shaped by similar biogeographic and climatic processes. The limited variation in foraging traits indicates that local adaptation to pines likely played a minor role. Our study shows that close mutualistic relationships between bird and plant species might not necessarily act as a primary driver of evolution and diversification in resource-specialized birds.

Genomic and phylogeographic analyses indicate that resource-specialization did not play a major role in the diversification and speciation of seed dispersing nutcrackers

Mendelian nightmares: the germline-restricted chromosome of songbirds

Germline-restricted chromosomes (GRCs) are accessory chromosomes that occur only in germ cells. They are eliminated from somatic cells through programmed DNA elimination during embryo development. GRCs have been observed in several unrelated animal taxa and show peculiar modes of non-Mendelian inheritance and within-individual elimination. Recent cytogenetic and phylogenomic evidence suggests that a GRC is present across the species-rich songbirds, but absent in non-passerine birds, implying that over half of all 10,500 bird species have extensive germline/soma genome differences. Here, we review recent insights gained from genomic, transcriptomic, and cytogenetic approaches with regard to the genetic content, phylogenetic distribution, and inheritance of the songbird GRC. While many questions remain unsolved in terms of GRC inheritance, elimination, and function, we discuss plausible scenarios and future directions for understanding this widespread form of programmed DNA elimination.

Phylogenomic discordance suggests polytomies along the backbone of the large genus Solanum

PREMISE

Evolutionary studies require solid phylogenetic frameworks, but increased volumes of phylogenomic data have revealed incongruent topologies among gene trees in many organisms both between and within genomes. Some of these incongruences indicate polytomies that may remain impossible to resolve. Here we investigate the degree of gene-tree discordance in Solanum, one of the largest flowering plant genera that includes the cultivated potato, tomato, and eggplant, as well as 24 minor crop plants.

METHODS

A densely sampled species-level phylogeny of Solanum is built using unpublished and publicly available Sanger sequences comprising 60% of all accepted species (742 spp.) and nine regions (ITS, waxy, and seven plastid markers). The robustness of this topology is tested by examining a full plastome dataset with 140 species and a nuclear target-capture dataset with 39 species of Solanum (Angiosperms353 probe set).

RESULTS

While the taxonomic framework of Solanum remained stable, gene tree conflicts and discordance between phylogenetic trees generated from the target-capture and plastome datasets were observed. The latter correspond to regions with short internodal branches, and network analysis and polytomy tests suggest the backbone is composed of three polytomies found at different evolutionary depths. The strongest area of discordance, near the crown node of Solanum, could potentially represent a hard polytomy.

CONCLUSIONS

We argue that incomplete lineage sorting due to rapid diversification is the most likely cause for these polytomies, and that embracing the uncertainty that underlies them is crucial to understand the evolution of large and rapidly radiating lineages.

OUP accepted manuscript

Phylogenetic analyses are increasingly being performed with data sets that incorporate hundreds of loci. Due to incomplete lineage sorting, hybridization, and horizontal gene transfer, the gene trees for these loci may often have topologies that differ from each other and from the species tree. The effect of these topological incongruences on divergence time estimation has not been fully investigated. Using a series of simulation experiments and empirical analyses, we demonstrate that when topological incongruence between gene trees and the species tree is not accounted for, the temporal duration of branches in regions of the species tree that are affected by incongruence is underestimated, whilst the duration of other branches is considerably overestimated. This effect becomes more pronounced with higher levels of topological incongruence. We show that this pattern results from the erroneous estimation of the number of substitutions along branches in the species tree, although the effect is modulated by the assumptions inherent to divergence time estimation, such as those relating to the fossil record or among-branch-substitution-rate variation. By only analyzing loci with gene trees that are topologically congruent with the species tree, or only taking into account the branches from each gene tree that are topologically congruent with the species tree, we demonstrate that the effects of topological incongruence can be ameliorated. Nonetheless, even when topologically congruent gene trees or topologically congruent branches are selected, error in divergence time estimates remains. This stems from temporal incongruences between divergence times in species trees and divergence times in gene trees, and more importantly, the difficulty of incorporating necessary assumptions for divergence time estimation. [Divergence time estimation; gene trees; species tree; topological incongruence.]

OUP accepted manuscript

Whole genome sequences are beginning to revolutionize our understanding of phylogenetic relationships. Yet, even whole genome sequences can fail to resolve the evolutionary history of the most rapidly radiating lineages, where incomplete lineage sorting, standing genetic variation, introgression, and other factors obscure the phylogenetic history of the group. To overcome such challenges, one emerging strategy is to integrate results across different methods. Most such approaches have been implemented on reduced representation genomic data sets, but whole genomes should provide the maximum possible evidence approach. Here, we test the ability of single nucleotide polymorphisms extracted from whole genome resequencing data, implemented in an integrative genomic approach, to resolve key nodes in the phylogeny of the mbuna, rock-dwelling cichlid fishes of Lake Malaŵi, which epitomize the phylogenetic intractability that often accompanies explosive lineage diversification. This monophyletic radiation has diversified at an unparalleled rate into several hundred species in less than 2 million years. Using an array of phylogenomic methods, we consistently recovered four major clades of mbuna, but a large basal polytomy among them. Although introgression between clades apparently contributed to the challenge of phylogenetic reconstruction, reduction of the data set to nonintrogressed sites still did not help to resolve the basal polytomy. On the other hand, relationships among six congeneric species pairs were resolved without ambiguity, even in one case where existing data led us to predict that resolution would be difficult. We conclude that the bursts of diversification at the earliest stages of the mbuna radiation may be phylogenetically unresolvable, but other regions of the tree are phylogenetically clearly supported. Integration of multiple phylogenomic approaches will continue to increase confidence in relationships inferred from these and other whole-genome data sets. [Incomplete lineage sorting; introgression; linkage disequilibrium; multispecies coalescence; rapid radiation; soft polytomy.]

Genome Stability Is in the Eye of the Beholder: CR1 Retrotransposon Activity Varies Significantly across Avian Diversity

Since the sequencing of the zebra finch genome it has become clear that avian genomes, while largely stable in terms of chromosome number and gene synteny, are more dynamic at an intrachromosomal level. A multitude of intrachromosomal rearrangements and significant variation in transposable element (TE) content have been noted across the avian tree. TEs are a source of genome plasticity, because their high similarity enables chromosomal rearrangements through nonallelic homologous recombination, and they have potential for exaptation as regulatory and coding sequences. Previous studies have investigated the activity of the dominant TE in birds, chicken repeat 1 (CR1) retrotransposons, either focusing on their expansion within single orders, or comparing passerines with nonpasserines. Here, we comprehensively investigate and compare the activity of CR1 expansion across orders of birds, finding levels of CR1 activity vary significantly both between and within orders. We describe high levels of TE expansion in genera which have speciated in the last 10 Myr including kiwis, geese, and Amazon parrots; low levels of TE expansion in songbirds across their diversification, and near inactivity of TEs in the cassowary and emu for millions of years. CR1s have remained active over long periods of time across most orders of neognaths, with activity at any one time dominated by one or two families of CR1s. Our findings of higher TE activity in species-rich clades and dominant families of TEs within lineages mirror past findings in mammals and indicate that genome evolution in amniotes relies on universal TE-driven processes.

Extensive Interspecific Gene Flow Shaped Complex Evolutionary History and Underestimated Species Diversity in Rapidly Radiated Dolphins

Recently diverged taxa are often characterized by high rates of hybridization, which can complicate phylogenetic reconstruction. For this reason, the phylogenetic relationships and evolutionary history of dolphins are still not very well resolved; the question of whether the genera Tursiops and Stenella are monophyletic is especially controversial. Here, we performed re-sequencing of six dolphin genomes and combined them with eight previously published dolphin SRA datasets and six whole-genome datasets to investigate the phylogenetic relationships of dolphins and test the monophyly hypothesis of Tursiops and Stenella. Phylogenetic reconstruction with the maximum likelihood and Bayesian methods of concatenated loci, as well as with coalescence analyses of sliding window trees, produced a concordant and well-supported tree. Our studies support the non-monophyletic status of Tursiops and Stenella because the species referred these genera do not form exclusive monophyletic clades. This suggests that the current taxonomy of both genera might not reflect their evolutionary history and may underestimate their diversity. A four-taxon D-statistic (ABBA-BABA) test, five-taxon DFOIL test, and tree-based PhyloNet analyses all showed extensive gene flow across dolphin species, which could explain the instability in resolving phylogenetic relationship of oceanic dolphins with different and limited markers. This study could be a good case to demonstrate how genomic data can reveal complex speciation and phylogeny in rapidly radiating animal groups.

Genome-wide ultraconserved elements resolve phylogenetic relationships and biogeographic history among Neotropical leaf-nosed bats in the genus Anoura (Phyllostomidae)

AnouraGray, 1838 are Neotropical nectarivorous bats and the most speciose genus within the phyllostomid subfamily Glossophaginae. However, Anoura species limits remain debated, and phylogenetic relationships remain poorly known, because previous studies used limited Anoura taxon sampling or focused primarily on higher-level relationships. Here, we conduct the first phylogenomic study of Anoura by analyzing 2039 genome-wide ultraconserved elements (UCEs) sequenced for 42 individuals from 8 Anoura species/lineages plus two outgroups. Overall, our results based on UCEs resolved relationships in the genus and supported (1) the monophyly of small-bodied Anoura species (previously genus Lonchoglossa); (2) monotypic status of A. caudifer; and (3) nested positions of "A. carishina", A. caudifer aequatoris, and A. geoffroyi peruana specimens within A. latidens, A. caudifer and A. geoffroyi, respectively (suggesting that these taxa are not distinct species). Additionally, (4) phylogenetic networks allowing reticulate edges did not explain gene tree discordance better than the species tree (without introgression), indicating that a coalescent model accounting for discordance solely through incomplete lineage sorting fit our data well. Sensitivity analyses indicated that our species tree results were not adversely affected by varying taxon sampling across loci. Tree calibration and Bayesian coalescent analyses dated the onset of diversification within Anoura to around ∼ 6-9 million years ago in the Miocene, with extant species diverging mainly within the past ∼ 4 million years. We inferred a historical biogeographical scenario for Anoura of parapatric speciation fragmenting the range of a wide-ranging ancestral lineage centered in the Central to Northern Andes, along with Pliocene-Pleistocene dispersal or founder event speciation in Amazonia and the Brazilian Atlantic forest during the last ∼ 2.5 million years.

What Have We Learned from the First 500 Avian Genomes?

The increased capacity of DNA sequencing has significantly advanced our understanding of the phylogeny of birds and the proximate and ultimate mechanisms molding their genomic diversity. In less than a decade, the number of available avian reference genomes has increased to over 500—approximately 5% of bird diversity—placing birds in a privileged position to advance the fields of phylogenomics and comparative, functional, and population genomics. Whole-genome sequence data, as well as indels and rare genomic changes, are further resolving the avian tree of life. The accumulation of bird genomes, increasingly with long-read sequence data, greatly improves the resolution of genomic features such as germline-restricted chromosomes and the W chromosome, and is facilitating the comparative integration of genotypes and phenotypes. Community-based initiatives such as the Bird 10,000 Genomes Project and Vertebrate Genome Project are playing a fundamental role in amplifying and coalescing a vibrant international program in avian comparative genomics.

Protein Structure, Models of Sequence Evolution, and Data Type Effects in Phylogenetic Analyses of Mitochondrial Data: A Case Study in Birds

Phylogenomic analyses have revolutionized the study of biodiversity, but they have revealed that estimated tree topologies can depend, at least in part, on the subset of the genome that is analyzed. For example, estimates of trees for avian orders differ if protein-coding or non-coding data are analyzed. The bird tree is a good study system because the historical signal for relationships among orders is very weak, which should permit subtle non-historical signals to be identified, while monophyly of orders is strongly corroborated, allowing identification of strong non-historical signals. Hydrophobic amino acids in mitochondrially-encoded proteins, which are expected to be found in transmembrane helices, have been hypothesized to be associated with non-historical signals. We tested this hypothesis by comparing the evolution of transmembrane helices and extramembrane segments of mitochondrial proteins from 420 bird species, sampled from most avian orders. We estimated amino acid exchangeabilities for both structural environments and assessed the performance of phylogenetic analysis using each data type. We compared those relative exchangeabilities with values calculated using a substitution matrix for transmembrane helices estimated using a variety of nuclear- and mitochondrially-encoded proteins, allowing us to compare the bird-specific mitochondrial models with a general model of transmembrane protein evolution. To complement our amino acid analyses, we examined the impact of protein structure on patterns of nucleotide evolution. Models of transmembrane and extramembrane sequence evolution for amino acids and nucleotides exhibited striking differences, but there was no evidence for strong topological data type effects. However, incorporating protein structure into analyses of mitochondrially-encoded proteins improved model fit. Thus, we believe that considering protein structure will improve analyses of mitogenomic data, both in birds and in other taxa.

Incorporating Topological and Age Uncertainty into Event-Based Biogeography of Sand Spiders Supports Paleo-Islands in Galapagos and Ancient Connections among Neotropical Dry Forests

Event-based biogeographic methods, such as dispersal-extinction-cladogenesis, have become increasingly popular for attempting to reconstruct the biogeographic history of organisms. Such methods employ distributional data of sampled species and a dated phylogenetic tree to estimate ancestral distribution ranges. Because the input tree is often a single consensus tree, uncertainty in topology and age estimates are rarely accounted for, even when they may affect the outcome of biogeographic estimates. Even when such uncertainties are taken into account for estimates of ancestral ranges, they are usually ignored when researchers compare competing biogeographic hypotheses. We explore the effect of incorporating this uncertainty in a biogeographic analysis of the 21 species of sand spiders (Sicariidae: Sicarius) from Neotropical xeric biomes, based on a total-evidence phylogeny including a complete sampling of the genus. Using a custom R script, we account for uncertainty in ages and topology by estimating ancestral ranges over a sample of trees from the posterior distribution of a Bayesian analysis, and for uncertainty in biogeographic estimates by using stochastic maps. This approach allows for counting biogeographic events such as dispersal among areas, counting lineages through time per area, and testing biogeographic hypotheses, while not overestimating the confidence in a single topology. Including uncertainty in ages indicates that Sicarius dispersed to the Galapagos Islands when the archipelago was formed by paleo-islands that are now submerged; model comparison strongly favors a scenario where dispersal took place before the current islands emerged. We also investigated past connections among currently disjunct Neotropical dry forests; failing to account for topological uncertainty underestimates possible connections among the Caatinga and Andean dry forests in favor of connections among Caatinga and Caribbean + Mesoamerican dry forests. Additionally, we find that biogeographic models including a founder-event speciation parameter (“+J”) are more prone to suffer from the overconfidence effects of estimating ancestral ranges using a single topology. This effect is alleviated by incorporating topological and age uncertainty while estimating stochastic maps, increasing the similarity in the inference of biogeographic events between models with or without a founder-event speciation parameter. We argue that incorporating phylogenetic uncertainty in biogeographic hypothesis-testing is valuable and should be a commonplace approach in the presence of rogue taxa or wide confidence intervals in age estimates, and especially when using models including founder-event speciation.

Comparing Ultraconserved Elements and Exons for Phylogenomic Analyses of Middle American Cichlids: When Data Agree to Disagree

Choosing among types of genomic markers to be used in a phylogenomic study can have a major influence on the cost, design, and results of a study. Yet few attempts have been made to compare categories of next-generation sequence markers limiting our ability to compare the suitability of these different genomic fragment types. Here, we explore properties of different genomic markers to find if they vary in the accuracy of component phylogenetic trees and to clarify the causes of conflict obtained from different data sets or inference methods. As a test case, we explore the causes of discordance between phylogenetic hypotheses obtained using a novel data set of ultraconserved elements (UCEs) and a recently published exon data set of the cichlid tribe Heroini. Resolving relationships among heroine cichlids has historically been difficult, and the processes of colonization and diversification in Middle America and the Greater Antilles are not yet well understood. Despite differences in informativeness and levels of gene tree discordance between UCEs and exons, the resulting phylogenomic hypotheses generally agree on most relationships. The independent data sets disagreed in areas with low phylogenetic signal that were overwhelmed by incomplete lineage sorting and nonphylogenetic signals. For UCEs, high levels of incomplete lineage sorting were found to be the major cause of gene tree discordance, whereas, for exons, nonphylogenetic signal is most likely caused by a reduced number of highly informative loci. This paucity of informative loci in exons might be due to heterogeneous substitution rates that are problematic to model (i.e., computationally restrictive) resulting in systematic errors that UCEs (being less informative individually but more uniform) are less prone to. These results generally demonstrate the robustness of phylogenomic methods to accommodate genomic markers with different biological and phylogenetic properties. However, we identify common and unique pitfalls of different categories of genomic fragments when inferring enigmatic phylogenetic relationships.

A re-assessment of the oldest therapsid Raranimus confirms its status as a basal member of the clade and fills Olson's gap

The non-mammalian therapsids comprise a paraphyletic assemblage of Permian-Jurassic synapsids closely related to mammals that includes six major clades of largely unresolved phylogenetic affinity. Understanding the early evolutionary radiation of therapsids is complicated by a gap in the fossil record during the Roadian (middle Permian) known as Olson's gap. Because of its early stratigraphic occurrence and its primitive features, Raranimus dashankouensis, from the Dashankou fauna (Rodian), Qingtoushan Formation (China), is currently considered the best candidate to fill this gap. However, it is known from only a single specimen, an isolated snout, which limits the amount of usable phylogenetic characters to reconstruct its affinities. In addition, understanding of the stratigraphy of the Qingtoushan Formation is poor. Here, we used CT scanning techniques to digitally reconstruct the bones and trigeminal canals of the snout of Raranimus in 3D. We confirm that Raranimus shares a high number of synapomorphies with more derived therapsids and is the only therapsid known so far to display a "pelycosaur"-like maxillary canal bearing a long caudal alveolar canal that gives off branches at regular intervals. This plesiomorphic feature supports the idea that Raranimus is basal to other therapsids.

Genomic phylogeography of the White-crowned Manakin Pseudopipra pipra (Aves: Pipridae) illuminates a continental-scale radiation out of the Andes

The complex landscape history of the Neotropics has generated opportunities for population isolation and diversification that place this region among the most species-rich in the world. Detailed phylogeographic studies are required to uncover the biogeographic histories of Neotropical taxa, to identify evolutionary correlates of diversity, and to reveal patterns of genetic connectivity, disjunction, and potential differentiation among lineages from different areas of endemism. The White-crowned Manakin (Pseudopipra pipra) is a small suboscine passerine bird that is broadly distributed through the subtropical rainforests of Central America, the lower montane cloud forests of the Andes from Colombia to central Peru, the lowlands of Amazonia and the Guianas, and the Atlantic forest of southeast Brazil. Pseudopipra is currently recognized as a single, polytypic biological species. We studied the effect of the Neotropical landscape on genetic and phenotypic differentiation within this species using genomic data derived from double digest restriction site associated DNA sequencing (ddRAD), and mitochondrial DNA. Most of the genetic breakpoints we identify among populations coincide with physical barriers to gene flow previously associated with avian areas of endemism. The phylogenetic relationships among these populations imply a novel pattern of Andean origination for this group, with subsequent diversification into the Amazonian lowlands. Our analysis of genomic admixture and gene flow reveals a complex history of introgression between some western Amazonian populations. These reticulate processes confound our application of standard concatenated and coalescent phylogenetic methods and raise the question of whether a lineage in the western Napo area of endemism should be considered a hybrid species. Lastly, analysis of variation in vocal and plumage phenotypes in the context of our phylogeny supports the hypothesis that Pseudopipra is a species-complex composed of at least 8, and perhaps up to 17 distinct species which have arisen in the last ∼2.5 Ma.

Data Types and the Phylogeny of Neoaves

The phylogeny of Neoaves, the largest clade of extant birds, has remained unclear despite intense study. The difficulty associated with resolving the early branches in Neoaves is likely driven by the rapid radiation of this group. However, conflicts among studies may be exacerbated by the data type analyzed. For example, analyses of coding exons typically yield trees that place Strisores (nightjars and allies) sister to the remaining Neoaves, while analyses of non-coding data typically yield trees where Mirandornites (flamingos and grebes) is the sister of the remaining Neoaves. Our understanding of data type effects is hampered by the fact that previous analyses have used different taxa, loci, and types of non-coding data. Herein, we provide strong corroboration of the data type effects hypothesis for Neoaves by comparing trees based on coding and non-coding data derived from the same taxa and gene regions. A simple analytical method known to minimize biases due to base composition (coding nucleotides as purines and pyrimidines) resulted in coding exon data with increased congruence to the non-coding topology using concatenated analyses. These results improve our understanding of the resolution of neoavian phylogeny and point to a challenge—data type effects—that is likely to be an important factor in phylogenetic analyses of birds (and many other taxonomic groups). Using our results, we provide a summary phylogeny that identifies well-corroborated relationships and highlights specific nodes where future efforts should focus.

Hybrid capture of 964 nuclear genes resolves evolutionary relationships in the mimosoid legumes and reveals the polytomous origins of a large pantropical radiation

PREMISE

Targeted enrichment methods facilitate sequencing of hundreds of nuclear loci to enhance phylogenetic resolution and elucidate why some parts of the "tree of life" are difficult (if not impossible) to resolve. The mimosoid legumes are a prominent pantropical clade of ~3300 species of woody angiosperms for which previous phylogenies have shown extensive lack of resolution, especially among the species-rich and taxonomically challenging ingoids.

METHODS

We generated transcriptomes to select low-copy nuclear genes, enrich these via hybrid capture for representative species of most mimosoid genera, and analyze the resulting data using de novo assembly and various phylogenomic tools for species tree inference. We also evaluate gene tree support and conflict for key internodes and use phylogenetic network analysis to investigate phylogenetic signal across the ingoids.

RESULTS

Our selection of 964 nuclear genes greatly improves phylogenetic resolution across the mimosoid phylogeny and shows that the ingoid clade can be resolved into several well-supported clades. However, nearly all loci show lack of phylogenetic signal for some of the deeper internodes within the ingoids.

CONCLUSIONS

Lack of resolution in the ingoid clade is most likely the result of hyperfast diversification, potentially causing a hard polytomy of six or seven lineages. The gene set for targeted sequencing presented here offers great potential to further enhance the phylogeny of mimosoids and the wider Caesalpinioideae with denser taxon sampling, to provide a framework for taxonomic reclassification, and to study the ingoid radiation.

The Origin of the Legumes is a Complex Paleopolyploid Phylogenomic Tangle Closely Associated with the Cretaceous-Paleogene (K-Pg) Mass Extinction Event

The consequences of the Cretaceous–Paleogene (K–Pg) boundary (KPB) mass extinction for the evolution of plant diversity remain poorly understood, even though evolutionary turnover of plant lineages at the KPB is central to understanding assembly of the Cenozoic biota. The apparent concentration of whole genome duplication (WGD) events around the KPB may have played a role in survival and subsequent diversification of plant lineages. To gain new insights into the origins of Cenozoic biodiversity, we examine the origin and early evolution of the globally diverse legume family (Leguminosae or Fabaceae). Legumes are ecologically (co-)dominant across many vegetation types, and the fossil record suggests that they rose to such prominence after the KPB in parallel with several well-studied animal clades including Placentalia and Neoaves. Furthermore, multiple WGD events are hypothesized to have occurred early in legume evolution. Using a recently inferred phylogenomic framework, we investigate the placement of WGDs during early legume evolution using gene tree reconciliation methods, gene count data and phylogenetic supernetwork reconstruction. Using 20 fossil calibrations we estimate a revised timeline of legume evolution based on 36 nuclear genes selected as informative and evolving in an approximately clock-like fashion. To establish the timing of WGDs we also date duplication nodes in gene trees. Results suggest either a pan-legume WGD event on the stem lineage of the family, or an allopolyploid event involving (some of) the earliest lineages within the crown group, with additional nested WGDs subtending subfamilies Papilionoideae and Detarioideae. Gene tree reconciliation methods that do not account for allopolyploidy may be misleading in inferring an earlier WGD event at the time of divergence of the two parental lineages of the polyploid, suggesting that the allopolyploid scenario is more likely. We show that the crown age of the legumes dates to the Maastrichtian or early Paleocene and that, apart from the Detarioideae WGD, paleopolyploidy occurred close to the KPB. We conclude that the early evolution of the legumes followed a complex history, in which multiple auto- and/or allopolyploidy events coincided with rapid diversification and in association with the mass extinction event at the KPB, ultimately underpinning the evolutionary success of the Leguminosae in the Cenozoic. [Allopolyploidy; Cretaceous–Paleogene (K–Pg) boundary; Fabaceae, Leguminosae; paleopolyploidy; phylogenomics; whole genome duplication events]

The Mitochondrial Genome of the Phytopathogenic Fungus Bipolaris sorokiniana and the Utility of Mitochondrial Genome to Infer Phylogeny of Dothideomycetes

A number of species in Bipolaris are important plant pathogens. Due to a limited number of synapomorphic characters, it is difficult to perform species identification and to estimate phylogeny of Bipolaris based solely on morphology. In this study, we sequenced the complete mitochondrial genome of Bipolaris sorokiniana, and presented the detailed annotation of the genome. The B. sorokiniana mitochondrial genome is 137,775 bp long, and contains two ribosomal RNA genes, 12 core protein-coding genes, 38 tRNA genes. In addition, two ribosomal protein genes (rps3 gene and rps5 gene) and the fungal mitochondrial RNase P gene (rnpB) are identified. The large genome size is mostly determined by the presence of numerous intronic and intergenic regions. A total of 28 introns are inserted in eight core protein-coding genes. Together with the published mitochondrial genome sequences, we conducted a preliminary phylogenetic inference of Dothideomycetes under various datasets and substitution models. The monophyly of Capnodiales, Botryosphaeriales and Pleosporales are consistently supported in all analyses. The Venturiaceae forms an independent lineage, with a distant phylogenetic relationship to Pleosporales. At the family level, the Mycosphaerellaceae, Botryosphaeriaceae. Phaeosphaeriaceae, and Pleosporaceae are recognized in the majority of trees.

Deep-Time Demographic Inference Suggests Ecological Release as Driver of Neoavian Adaptive Radiation

Assessing the applicability of theory to major adaptive radiations in deep time represents an extremely difficult problem in evolutionary biology. Neoaves, which includes 95% of living birds, is believed to have undergone a period of rapid diversification roughly coincident with the Cretaceous–Paleogene (K-Pg) boundary. We investigate whether basal neoavian lineages experienced an ecological release in response to ecological opportunity, as evidenced by density compensation. We estimated effective population sizes (Ne) of basal neoavian lineages by combining coalescent branch lengths (CBLs) and the numbers of generations between successive divergences. We used a modified version of Accurate Species TRee Algorithm (ASTRAL) to estimate CBLs directly from insertion–deletion (indel) data, as well as from gene trees using DNA sequence and/or indel data. We found that some divergences near the K-Pg boundary involved unexpectedly high gene tree discordance relative to the estimated number of generations between speciation events. The simplest explanation for this result is an increase in Ne, despite the caveats discussed herein. It appears that at least some early neoavian lineages, similar to the ancestor of the clade comprising doves, mesites, and sandgrouse, experienced ecological release near the time of the K-Pg mass extinction.

Avian Binocularity and Adaptation to Nocturnal Environments: Genomic Insights from a Highly Derived Visual Phenotype

Typical avian eyes are phenotypically engineered for photopic vision (daylight). In contrast, the highly derived eyes of the barn owl (Tyto alba) are adapted for scotopic vision (dim light). The dramatic modifications distinguishing barn owl eyes from other birds include: 1) shifts in frontal orientation to improve binocularity, 2) rod-dominated retina, and 3) enlarged corneas and lenses. Some of these features parallel mammalian eye patterns, which are hypothesized to have initially evolved in nocturnal environments. Here, we used an integrative approach combining phylogenomics and functional phenotypes of 211 eye-development genes across 48 avian genomes representing most avian orders, including the stem lineage of the scotopic-adapted barn owl. Overall, we identified 25 eye-development genes that coevolved under intensified or relaxed selection in the retina, lens, cornea, and optic nerves of the barn owl. The agtpbp1 gene, which is associated with the survival of photoreceptor populations, was pseudogenized in the barn owl genome. Our results further revealed that barn owl retinal genes responsible for the maintenance, proliferation, and differentiation of photoreceptors experienced an evolutionary relaxation. Signatures of relaxed selection were also observed in the lens and cornea morphology-associated genes, suggesting that adaptive evolution in these structures was essentially structural. Four eye-development genes (ephb1, phactr4, prph2, and rs1) evolved in positive association with the orbit convergence in birds and under relaxed selection in the barn owl lineage, likely contributing to an increased reliance on binocular vision in the barn owl. Moreover, we found evidence of coevolutionary interactions among genes that are expressed in the retina, lens, and optic nerve, suggesting synergetic adaptive events. Our study disentangles the genomic changes governing the binocularity and low-light perception adaptations of barn owls to nocturnal environments while revealing the molecular mechanisms contributing to the shift from the typical avian photopic vision to the more-novel scotopic-adapted eye.

Large-scale genomic sequence data resolve the deepest divergences in the legume phylogeny and support a near-simultaneous evolutionary origin of all six subfamilies

Phylogenomics is increasingly used to infer deep-branching relationships while revealing the complexity of evolutionary processes such as incomplete lineage sorting, hybridization/introgression and polyploidization. We investigate the deep-branching relationships among subfamilies of the Leguminosae (or Fabaceae), the third largest angiosperm family. Despite their ecological and economic importance, a robust phylogenetic framework for legumes based on genome-scale sequence data is lacking. We generated alignments of 72 chloroplast genes and 7621 homologous nuclear-encoded proteins, for 157 and 76 taxa, respectively. We analysed these with maximum likelihood, Bayesian inference, and a multispecies coalescent summary method, and evaluated support for alternative topologies across gene trees. We resolve the deepest divergences in the legume phylogeny despite lack of phylogenetic signal across all chloroplast genes and the majority of nuclear genes. Strongly supported conflict in the remainder of nuclear genes is suggestive of incomplete lineage sorting. All six subfamilies originated nearly simultaneously, suggesting that the prevailing view of some subfamilies as 'basal' or 'early-diverging' with respect to others should be abandoned, which has important implications for understanding the evolution of legume diversity and traits. Our study highlights the limits of phylogenetic resolution in relation to rapid successive speciation.

Tanopicobia gen. nov., a new genus of quill mites, its phylogenetic placement in the subfamily Picobiinae (Acariformes: Syringophilidae) and picobiine relationships with avian hosts

A new monotypic genus Tanopicobia gen. nov. is established for a new species Tanopicobia trachyphoni sp. nov., parasitizing Trachyphonus erythrocephalus Cabanis, 1878 (Piciformes: Lybiidae) from Tanzania. In phylogenetic analyses based on morphological data and constructed using the maximum parsimony approach, this taxon falls within the subfamily Picobiinae Johnston and Kethley, 1973 in the Neopicobia-species-group as closely related to the genus Pipicobia Glowska and Schmidt, 2014. Tanopicobia differs from Pipicobia by the following features in females: genital setae absent; setae ve are situated far and posteromedial to the level of setal bases vi; setae 3a are thick and knobbed. Additionally, a new generic key for subfamily Picobiinae is constructed and general host-parasite ecological and phylogenetic relationships are discussed. Picobiines are present in several lineages of neoavian birds, from basal Galloanseres to terminal Telluraves, which are infested by 70 (89.7% of all) species of these ectoparasites.

The uropygial gland of the monk parakeet Myiopsitta monachus: Histology, morphogenesis, and evolution within Psittaciformes (Aves)

We describe the morphology, histology, and histochemical characteristics of the uropygial gland (UG) of the monk parakeet Myiopsitta monachus. The UG has a heart-shape external appearance and adenomers extensively branched with a convoluted path, covered by a stratified epithelium formed by different cellular strata and divided into three zones (based on the epithelial height and lumen width), a cylindrical papilla with an internal structure of delicate type and two excretory pores surrounded by a feather tuft. Histochemical and lectin-histochemical techniques performed showed positivity against PAS, AB pH 2.5, AB-PAS, and some lectines, likely related to the granivorous feeding habits. Also, we describe the morphogenesis of the UG of the monk parakeet, which appears at embryological stage 34 as a pair of ectodermal invaginations. Heterochronic events in the onset development of the UG when compared with other birds could be recognized. Finally, to examine the phylogenetic occurrence of the UG within the Psittaciformes and infer its evolutionary history, we mapped its presence/absence over a molecular phylogeny. The reconstruction of the characters states at ancestral nodes revealed that the presence of the UG was the plesiomorphic feature for Psittaciformes and its loss evolved independently more than once.

Whole-Genome Analyses Resolve the Phylogeny of Flightless Birds (Palaeognathae) in the Presence of an Empirical Anomaly Zone

Palaeognathae represent one of the two basal lineages in modern birds, and comprise the volant (flighted) tinamous and the flightless ratites. Resolving palaeognath phylogenetic relationships has historically proved difficult, and short internal branches separating major palaeognath lineages in previous molecular phylogenies suggest that extensive incomplete lineage sorting (ILS) might have accompanied a rapid ancient divergence. Here, we investigate palaeognath relationships using genome-wide data sets of three types of noncoding nuclear markers, together totaling 20,850 loci and over 41 million base pairs of aligned sequence data. We recover a fully resolved topology placing rheas as the sister to kiwi and emu + cassowary that is congruent across marker types for two species tree methods (MP-EST and ASTRAL-II). This topology is corroborated by patterns of insertions for 4274 CR1 retroelements identified from multispecies whole-genome screening, and is robustly supported by phylogenomic subsampling analyses, with MP-EST demonstrating particularly consistent performance across subsampling replicates as compared to ASTRAL. In contrast, analyses of concatenated data supermatrices recover rheas as the sister to all other nonostrich palaeognaths, an alternative that lacks retroelement support and shows inconsistent behavior under subsampling approaches. While statistically supporting the species tree topology, conflicting patterns of retroelement insertions also occur and imply high amounts of ILS across short successive internal branches, consistent with observed patterns of gene tree heterogeneity. Coalescent simulations and topology tests indicate that the majority of observed topological incongruence among gene trees is consistent with coalescent variation rather than arising from gene tree estimation error alone, and estimated branch lengths for short successive internodes in the inferred species tree fall within the theoretical range encompassing the anomaly zone. Distributions of empirical gene trees confirm that the most common gene tree topology for each marker type differs from the species tree, signifying the existence of an empirical anomaly zone in palaeognaths.

Phylogenomic species delimitation in microendemic frogs of the Brazilian Atlantic Forest

The advent of next-generation sequencing allows researchers to use large-scale datasets for species delimitation analyses, yet one can envision an inflection point where the added accuracy of including more loci does not offset the increased computational burden. One alternative to including all loci could be to prioritize the analysis of loci for which there is an expectation of high informativeness. Here, we explore the issue of species delimitation and locus selection with montane species from two anuran genera that have been isolated in sky islands across the southern Brazilian Atlantic Forest: Melanophryniscus (Bufonidae) and Brachycephalus (Brachycephalidae). To delimit species, we obtained genetic data using target enrichment of ultraconserved elements from 32 populations (13 for Melanophryniscus and 19 for Brachycephalus), and we were able to create datasets that included over 800 loci with no missing data. We ranked loci according to their number of parsimony-informative sites, and we performed species delimitation analyses using BPP with the most informative 10, 20, 40, 80, 160, 320, and 640 loci. We identified three types of phylogenetic node: nodes with either consistently high or low support regardless of the number of loci or their informativeness and nodes that were initially poorly supported where support became stronger as we included more data. When viewed across all sensitivity analyses, our results suggest that the current species richness in both genera is likely underestimated. In addition, our results show the effects of different sampling strategies on species delimitation using phylogenomic datasets.

Phylogenomics Uncovers Confidence and Conflict in the Rapid Radiation of Australo-Papuan Rodents

The estimation of robust and accurate measures of branch support has proven challenging in the era of phylogenomics. In data sets of potentially millions of sites, bootstrap support for bifurcating relationships around very short internal branches can be inappropriately inflated. Such overestimation of branch support may be particularly problematic in rapid radiations, where phylogenetic signal is low and incomplete lineage sorting severe. Here, we explore this issue by comparing various branch support estimates under both concatenated and coalescent frameworks, in the recent radiation Australo-Papuan murine rodents (Muridae: Hydromyini). Using nucleotide sequence data from 1245 independent loci and several phylogenomic inference methods, we unequivocally resolve the majority of genus-level relationships within Hydromyini. However, at four nodes we recover inconsistency in branch support estimates both within and among concatenated and coalescent approaches. In most cases, concatenated likelihood approaches using standard fast bootstrap algorithms did not detect any uncertainty at these four nodes, regardless of partitioning strategy. However, we found this could be overcome with two-stage resampling, that is, across genes and sites within genes (using -bsam GENESITE in IQ-TREE). In addition, low confidence at recalcitrant nodes was recovered using UFBoot2, a recent revision to the bootstrap protocol in IQ-TREE, but this depended on partitioning strategy. Summary coalescent approaches also failed to detect uncertainty under some circumstances. For each of four recalcitrant nodes, an equivalent (or close to equivalent) number of genes were in strong support ($>$ 75% bootstrap) of both the primary and at least one alternative topological hypothesis, suggesting notable phylogenetic conflict among loci not detected using some standard branch support metrics. Recent debate has focused on the appropriateness of concatenated versus multigenealogical approaches to resolving species relationships, but less so on accurately estimating uncertainty in large data sets. Our results demonstrate the importance of employing multiple approaches when assessing confidence and highlight the need for greater attention to the development of robust measures of uncertainty in the era of phylogenomics.

Comparative Phylogenomics, a Stepping Stone for Bird Biodiversity Studies

Birds are a group with immense availability of genomic resources, and hundreds of forthcoming genomes at the doorstep. We review recent developments in whole genome sequencing, phylogenomics, and comparative genomics of birds. Short read based genome assemblies are common, largely due to efforts of the Bird 10K genome project (B10K). Chromosome-level assemblies are expected to increase due to improved long-read sequencing. The available genomic data has enabled the reconstruction of the bird tree of life with increasing confidence and resolution, but challenges remain in the early splits of Neoaves due to their explosive diversification after the Cretaceous-Paleogene (K-Pg) event. Continued genomic sampling of the bird tree of life will not just better reflect their evolutionary history but also shine new light onto the organization of phylogenetic signal and conflict across the genome. The comparatively simple architecture of avian genomes makes them a powerful system to study the molecular foundation of bird specific traits. Birds are on the verge of becoming an extremely resourceful system to study biodiversity from the nucleotide up.

A Robust Phylogenomic Time Tree for Biotechnologically and Medically Important Fungi in the Genera Aspergillus and Penicillium

Understanding the evolution of traits across technologically and medically significant fungi requires a robust phylogeny. Even though species in the Aspergillus and Penicillium genera (family Aspergillaceae, class Eurotiomycetes) are some of the most significant technologically and medically relevant fungi, we still lack a genome-scale phylogeny of the lineage or knowledge of the parts of the phylogeny that exhibit conflict among analyses. Here, we used a phylogenomic approach to infer evolutionary relationships among 81 genomes that span the diversity of Aspergillus and Penicillium species, to identify conflicts in the phylogeny, and to determine the likely underlying factors of the observed conflicts. Using a data matrix comprised of 1,668 genes, we found that while most branches of the phylogeny of the Aspergillaceae are robustly supported and recovered irrespective of method of analysis, a few exhibit various degrees of conflict among our analyses. Further examination of the observed conflict revealed that it largely stems from incomplete lineage sorting and hybridization or introgression. Our analyses provide a robust and comprehensive evolutionary genomic roadmap for this important lineage, which will facilitate the examination of the diverse technologically and medically relevant traits of these fungi in an evolutionary context.

The filamentous fungal family Aspergillaceae contains >1,000 known species, mostly in the genera Aspergillus and Penicillium. Several species are used in the food, biotechnology, and drug industries (e.g., Aspergillus oryzae and Penicillium camemberti), while others are dangerous human and plant pathogens (e.g., Aspergillus fumigatus and Penicillium digitatum). To infer a robust phylogeny and pinpoint poorly resolved branches and their likely underlying contributors, we used 81 genomes spanning the diversity of Aspergillus and Penicillium to construct a 1,668-gene data matrix. Phylogenies of the nucleotide and amino acid versions of this full data matrix as well as of several additional data matrices were generated using three different maximum likelihood schemes (i.e., gene-partitioned, unpartitioned, and coalescence) and using both site-homogenous and site-heterogeneous models (total of 64 species-level phylogenies). Examination of the topological agreement among these phylogenies and measures of internode certainty identified 11/78 (14.1%) bipartitions that were incongruent and pinpointed the likely underlying contributing factors, which included incomplete lineage sorting, hidden paralogy, hybridization or introgression, and reconstruction artifacts associated with poor taxon sampling. Relaxed molecular clock analyses suggest that Aspergillaceae likely originated in the lower Cretaceous and that the Aspergillus and Penicillium genera originated in the upper Cretaceous. Our results shed light on the ongoing debate on Aspergillus systematics and taxonomy and provide a robust evolutionary and temporal framework for comparative genomic analyses in Aspergillaceae. More broadly, our approach provides a general template for phylogenomic identification of resolved and contentious branches in densely genome-sequenced lineages across the tree of life.

Phylogenetic Signal of Indels and the Neoavian Radiation

The early radiation of Neoaves has been hypothesized to be an intractable “hard polytomy”. We explore the fundamental properties of insertion/deletion alleles (indels), an under-utilized form of genomic data with the potential to help solve this. We scored >5 million indels from >7000 pan-genomic intronic and ultraconserved element (UCE) loci in 48 representatives of all neoavian orders. We found that intronic and UCE indels exhibited less homoplasy than nucleotide (nt) data. Gene trees estimated using indel data were less resolved than those estimated using nt data. Nevertheless, Accurate Species TRee Algorithm (ASTRAL) species trees estimated using indels were generally similar to nt-based ASTRAL trees, albeit with lower support. However, the power of indel gene trees became clear when we combined them with nt gene trees, including a striking result for UCEs. The individual UCE indel and nt ASTRAL trees were incongruent with each other and with the intron ASTRAL trees; however, the combined indel+nt ASTRAL tree was much more congruent with the intronic trees. Finally, combining indel and nt data for both introns and UCEs provided sufficient power to reduce the scope of the polytomy that was previously proposed for several supraordinal lineages of Neoaves.

New Material of Paleocene-Eocene Pellornis (Aves: Gruiformes) Clarifies the Pattern and Timing of the Extant Gruiform Radiation

Pellornis mikkelseni is an early gruiform from the latest Paleocene-earliest Eocene Fur Formation of Denmark. At approximately 54 million years old, it is among the earliest clear records of the Gruiformes. The holotype specimen, and only material thus far recognised, was originally considered to comprise a partial postcranial skeleton. However, additional mechanical preparation of the nodule containing the holotype revealed that the skeleton is nearly complete and includes a well-preserved skull. In addition to extracting new information from the holotype, we identify and describe two additional specimens of P. mikkelseni which reveal further morphological details of the skeleton. Together, these specimens show that P. mikkelseni possessed a schizorhinal skull and shared many features with the well-known Paleogene Messelornithidae (“Messel rails”). To reassess the phylogenetic position of P. mikkelseni, we modified an existing morphological dataset by adding 20 characters, four extant gruiform taxa, six extinct gruiform taxa, and novel scorings based on the holotype and referred specimens. Phylogenetic analyses recover a clade containing P. mikkelseni, Messelornis, Songzia and crown Ralloidea, supporting P. mikkelseni as a crown gruiform. The phylogenetic position of P. mikkelseni illustrates that some recent divergence time analyses have underestimated the age of crown Gruiformes. Our results suggest a Paleocene origin for this important clade, bolstering evidence for a rapid early radiation of Neoaves following the end-Cretaceous mass extinction.