The relationships of passerines (such as the well-studied zebra finch) with non-passerine birds is one of the great enigmas of avian phylogenetic research, because decades of extensive morphological and molecular studies yielded highly inconsistent results between and within data sets. Here we show the first application of the virtually homoplasy-free retroposon insertions to this controversy. Our study examined ~200,000 retroposon-containing loci from various avian genomes and retrieved 51 markers resolving early bird phylogeny. Among these, we obtained statistically significant evidence that parrots are the closest and falcons the second-closest relatives of passerines, together constituting the Psittacopasserae and the Eufalconimorphae, respectively. Our new and robust phylogenetic framework has substantial implications for the interpretation of various conclusions drawn from passerines as model organisms. This includes insights of relevance to human neuroscience, as vocal learning (that is, birdsong) probably evolved in the psittacopasseran ancestor, >30 million years earlier than previously assumed.

Graphical AbstractHighlights d Interstitial annelid taxa are secondarily derived contra the archiannelid hypothesis d Several annelid clades adapted to the space between the sand grains by progenesis d Other interstitial annelids evolved by miniaturization of ancestral adult stages d Miniaturization is as important as progenesis in the adaptation to the interstitium

Low character variation among onychophoran species has been an obstacle for taxonomic and phylogenetic studies in the past, however we have identified a number of new and informative characters using morphological, molecular, and chromosomal techniques. Our analyses involved a detailed examination of Epiperipatus biolleyi from Costa Rica, Eoperipatus sp. from Thailand, and a new onychophoran species and genus from Costa Rica, Principapillatus hitoyensis gen. et sp. nov.. Scanning electron microscopy on embryos and specimens of varying age revealed novel morphological characters and character states, including the distribution of different receptor types along the antennae, the arrangement and form of papillae on the head, body and legs, the presence and shape of interpedal structures and fields of modified scales on the ventral body surface, the arrangement of lips around the mouth, the number, position and structure of crural tubercles and anal gland openings, and the presence and shape of embryonic foot projections. Karyotypic analyses revealed differences in the number and size of chromosomes among the species studied. The results of our phylogenetic analyses using mitochondrial COI and 12S rRNA gene sequences are in line with morphological and karyotype data. However, our data show a large number of unexplored, albeit informative, characters in the Peripatidae. We suggest that analysing these characters in additional species would help unravel species diversity and phylogeny in the Onychophora, and that inconsistencies among most diagnostic features used for the peripatid genera in the literature could be addressed by identifying a suite of characters common to all peripatids.

No abstract

Pax family genes encode a class of transcription factors that regulate various developmental processes. To shed light on the evolutionary history of these genes in Panarthropoda (Onychophora + Tardigrada + Arthropoda), we analyzed the Pax repertoire in the embryonic and adult transcriptomes of the onychophoran Euperipatoides rowelli. Our data revealed homologs of all five major bilaterian Pax subfamilies in this species, including Pax2/5/8, Pax4/6, Pox-neuro, Pax1/9/Pox-meso, and Pax3/7. In addition, we identified a new Pax member, pax-α, which does not fall into any other known Pax subfamily but instead clusters in the heterogenic Pax-α/β clade containing deuterostome, ecdysozoan, and lophotrochozoan gene sequences. These findings suggest that the last common bilaterian ancestor possessed six rather than five Pax genes, which have been retained in the panarthropod lineage. The expression data of Pax orthologs in the onychophoran embryo revealed distinctive patterns, some of which might be related to their ancestral roles in the last common panarthropod ancestor, whereas others might be specific to the onychophoran lineage. The derived roles include, for example, an involvement of pax2/5/8, pox-neuro, and pax3/7 in onychophoran nephridiogenesis, and an additional function of pax2/5/8 in the formation of the ventral and preventral organs. Furthermore, our transcriptomic analyses suggest that at least some Pax genes, including pax6 and pax-α, are expressed in the adult onychophoran head, although the corresponding functions remain to be clarified. The remarkable diversity of the Pax expression patterns highlights the functional and evolutionary plasticity of these genes in panarthropods.

Arthropods typically show two types of segmentation: the embryonic parasegments and the adult segments that lie out of register with each other. Such a dual nature of body segmentation has not been described from Onychophora, one of the closest arthropod relatives. Hence, it is unclear whether onychophorans have segments, parasegments, or both, and which of these features was present in the last common ancestor of Onychophora and Arthropoda. To address this issue, we analysed the expression patterns of the “segment polarity genes” engrailed, cubitus interruptus, wingless and hedgehog in embryos of the onychophoran Euperipatoides rowelli. Our data revealed that these genes are expressed in repeated sets with a specific anterior-to-posterior order along the body in embryos of E. rowelli. In contrast to arthropods, the expression occurs after the segmental boundaries have formed. Moreover, the initial segmental furrow retains its position within the engrailed domain throughout development, whereas no new furrow is formed posterior to this domain. This suggests that no re-segmentation of the embryo occurs in E. rowelli. Irrespective of whether or not there is a morphological or genetic manifestation of parasegments in Onychophora, our data clearly show that parasegments, even if present, cannot be regarded as the initial metameric units of the onychophoran embryo, because the expression of key genes that define the parasegmental boundaries in arthropods occurs after the segmental boundaries have formed. This is in contrast to arthropods, in which parasegments rather than segments are the initial metameric units of the embryo. Our data further revealed that the expression patterns of “segment polarity genes” correspond to organogenesis rather than segment formation. This is in line with the concept of segmentation as a result of concerted evolution of individual periodic structures rather than with the interpretation of ‘segments’ as holistic units.

Convergent phenotypic evolution is often caused by recurrent changes at particular nodes in the underlying gene regulatory networks (GRNs). The genes at such evolutionary ‘hotspots’ are thought to maximally affect the phenotype with minimal pleiotropic consequences. This has led to the suggestion that if a GRN is understood in sufficient detail, the path of evolution may be predictable. The repeated evolutionary loss of larval trichomes among Drosophila species is caused by the loss of shavenbaby (svb) expression. svb is also required for development of leg trichomes, but the evolutionary gain of trichomes in the ‘naked valley’ on T2 femurs in Drosophila melanogaster is caused by reduced microRNA-92a (miR-92a) expression rather than changes in svb. We compared the expression and function of components between the larval and leg trichome GRNs to investigate why the genetic basis of trichome pattern evolution differs in these developmental contexts. We found key differences between the two networks in both the genes employed, and in the regulation and function of common genes. These differences in the GRNs reveal why mutations in svb are unlikely to contribute to leg trichome evolution and how instead miR-92a represents the key evolutionary switch in this context. Our work shows that variability in GRNs across different developmental contexts, as well as whether a morphological feature is lost versus gained, influence the nodes at which a GRN evolves to cause morphological change. Therefore, our findings have important implications for understanding the pathways and predictability of evolution.

Convergent phenotypic evolution is often caused by recurrent changes at particular nodes in the underlying gene regulatory networks (GRNs). The genes at such evolutionary 'hotspots' are thought to maximally affect the phenotype with minimal pleiotropic consequences. This has led to the suggestion that if a GRN is understood in sufficient detail, the path of evolution may be predictable. The repeated loss of larval trichomes among Drosophila species is caused by the loss of shavenbaby (svb) expression. svb is also required for development of leg trichomes, but the evolutionary gain of trichomes in the 'naked valley' on T2 femurs in Drosophila melanogaster is caused by the loss of microRNA-92a (miR-92a) expression rather than changes in svb. We compared the architectures of the larval and leg trichome GRNs to investigate why the genetic basis of trichome pattern evolution differs in these developmental contexts. We found key differences between these two networks in both the genes employed, and in the regulation and function of common genes. These differences in the GRNs reveal why mutations in svb are unlikely to contribute to leg trichome evolution and how instead miR-92a represents the key evolutionary switch in this context. Our work shows that differences in the components and wiring of GRNs in different developmental contexts, as well as whether a morphological feature is lost versus gained, influence the nodes at which a GRN evolves to cause morphological change. Therefore our findings have important implications for understanding the pathways and predictability of evolution. Significance StatementA major goal of biology is to identify the genetic cause of organismal diversity. Convergent evolution of traits is often caused by changes in the same genes -evolutionary 'hotspots'. shavenbaby is a 'hotspot' for larval trichome loss in Drosophila, however microRNA-92a underlies the gain of leg trichomes. To understand this difference in the genetics of phenotypic evolution, we compared the underlying gene regulatory networks (GRNs). We found that differences in GRN architecture in different developmental contexts, and whether a trait is lost or gained, influence the pathway of evolution. Therefore hotspots in one context may not readily evolve in a different context. This has important implications for understanding the genetic basis of phenotypic change and the predictability of evolution.