Search

Recombination breaks down genetic linkage by reshuffling existing variants onto new genetic backgrounds. These dynamics are traditionally quantified by examining the correlations between alleles, and how they decay as a function of the recombination rate. However, the magnitudes of these correlations are strongly influenced by other evolutionary forces like natural selection and genetic drift, making it difficult to tease out the effects of recombination. Here, we introduce a theoretical framework for analyzing an alternative family of statistics that measure the homoplasy produced by recombination. We derive analytical expressions that predict how these statistics depend on the rates of recombination and recurrent mutation, the strength of negative selection and genetic drift, and the present-day frequencies of the mutant alleles. We find that the degree of homoplasy can strongly depend on this frequency scale, which reflects the underlying timescales over which these mutations occurred. We show how these scaling properties can be used to isolate the effects of recombination and discuss their implications for the rates of horizontal gene transfer in bacteria., Competing Interests: Conflicts of interest The author(s) declare no conflicts of interest., (© The Author(s) 2024. Published by Oxford University Press on behalf of The Genetics Society of America. All rights reserved. For commercial re-use, please contact reprints@oup.com for reprints and translation rights for reprints. All other permissions can be obtained through our RightsLink service via the Permissions link on the article page on our site—for further information please contact journals.permissions@oup.com.)

New Zealand scincid lizards, genus Oligosoma, represent a monophyletic radiation of a clade, Eugongylini, of species distributed geographically throughout the South Pacific with major radiations in Australia and New Caledonia. Viviparity has evolved independently on multiple occasions within these lineages. Studies of Australian species have revealed that placental specializations resulting in substantial placentotrophy have evolved in two lineages. The pattern of extraembryonic membrane development of oviparous species differs from viviparous species and identical placental architecture has evolved in both placentotrophic lineages. We analyzed extraembryonic membrane development in two New Zealand species, the sole oviparous species, Oligosoma suteri, and placental development of a representative viviparous species, Oligosoma polychroma, using histological techniques. We conclude that these two species share a basic pattern of extraembryonic membrane development with other squamates. Comparisons with Australian species indicate that morphogenesis of the yolk sac of O. suteri results in an elaborate structure previously known only in Oligosoma lichenigerum with a geographic distribution on Lord Howe Island and Norfolk Island. This finding supports a close relationship between these two taxa. We conclude also that the pattern of placental development of O. polychroma is identical to that of viviparous species of Australia. The terminal placental stage for each of these lineages includes a chorioallantoic placenta and an elaborate omphaloplacenta. This level of homoplasy in placental evolution is consistent with a hypothesis that selection favors regional differentiation of the maternal-embryonic interface and that the omphaloplacenta is an adaptation for histotrophic transport., (© 2024 Wiley Periodicals LLC.)

In addition to the type locality (the summit of Aprada-tepui, Bolívar State of Venezuela), the distribution of the egg-brooding frog Stefania satelles was long thought to include several isolated tabletop mountain (tepui) summits surrounding the large Chimantá Massif in Bolívar State (hence the Latin name "satelles"). However, multilocus molecular phylogenetic analyses have revealed that this taxon includes several undescribed morphologically cryptic species, and that S. satelles should be restricted to its type locality. Two tepui-summit species confused under that name in the literature remain to be named, and the present paper aims at describing these populations previously referred to as Stefania sp. 3 and S. sp. 5. Stefania sp. 3 is only known from the small summit of Angasima-tepui, while S. sp. 5 is only reported from the small summit of Upuigma-tepui, both mountains being located south of the Chimantá Massif. These new, phylogenetically distinct species are described based on external morphology and osteology and in comparison to close relatives in the S. ginesi clade, which consists exclusively of tepui summit species. Both new species have highly restricted geographic ranges (less than 3 km 2 ) and should be listed as Critically Endangered according to IUCN criteria., (© 2024. The Author(s).)

Leaf venation is a pivotal trait in the success of vascular plants. Whereas gymnosperms have single or sparsely branched parallel veins, angiosperms developed a hierarchical structure of veins that form a complex reticulum. Its physiological consequences are considered to have enabled angiosperms to dominate terrestrial ecosystems in the Late Cretaceous and Cenozoic. Although a hierarchical-reticulate venation also occurs in some groups of extinct seed plants, it is unclear whether these are stem relatives of angiosperms or have evolved these traits in parallel. Here, we re-examine the morphology of the enigmatic foliage taxon Furcula, a potential early Mesozoic angiosperm relative, and argue that its hierarchical vein network represents convergent evolution (in the Late Triassic) with flowering plants (which developed in the Early Cretaceous) based on details of vein architecture and the absence of angiosperm-like stomata and guard cells. We suggest that its nearest relatives are Peltaspermales similar to Scytophyllum and Vittaephyllum, the latter being a genus that originated during the Late Triassic (Carnian) and shares a hierarchical vein system with Furcula. We further suggest that the evolution of hierarchical venation systems in the early Permian, the Late Triassic, and the Early Cretaceous represent 'natural experiments' that might help resolve the selective pressures enabling this trait to evolve., (© 2024 The Authors New Phytologist © 2024 New Phytologist Foundation.)

Although molecular phylogenetics remains the most widely used method of inferring the evolutionary history of living groups, the last decade has seen a renewed interest in morphological phylogenetics, mostly driven by the promises that integrating the fossil record in phylogenetic trees offers to our understanding of macroevolutionary processes and dynamics and the possibility that the inclusion of fossil taxa could lead to more accurate phylogenetic hypotheses. The plant fossil record presents some challenges to its integration in a phylogenetic framework. Phylogenies including plant fossils often retrieve uncertain relationships with low support, or lack of resolution. This low support is due to the pervasiveness of morphological convergence among plant organs and the fragmentary nature of many plant fossils, and it is often perceived as a fundamental weakness reducing the utility of plant fossils in phylogenetics. Here I discuss the importance of uncertainty in morphological phylogenetics and how we can identify important information from different patterns and types of uncertainty. I also review a set of methodologies that can allow us to understand the causes underpinning uncertainty and how these practices can help us to further our knowledge of plant fossils. I also propose that a new visual language, including the use of networks instead of trees, represents an improvement on the old visualization based on consensus trees and more adequately serves phylogeneticists working with plant fossils. This set of methods and visualization tools represents an important way forward in a fundamental field for our understanding of the evolutionary history of plants., (© 2024 The Authors. American Journal of Botany published by Wiley Periodicals LLC on behalf of Botanical Society of America.)

Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.

Background: A general view in the study of pollination syndromes is that floral traits usually represent convergent floral adaptations to specific functional pollinator groups. However, the definition of convergence is elusive and contradictory in the literature. Is convergence the independent evolution of either the same trait or similar traits with the same function? A review of the concept of convergence in developmental biology and phylogenetic systematics may shed new light in studies of pollination syndromes., Scope: The aims of this article are (1) to explore the notion of convergence and other concepts (analogy, homoplasy and parallelism) within the theory and practice of developmental evolution and phylogenetic systematics; (2) to modify the definitions of syndromes in order to embrace the concepts of analogy and convergence; (3) to revisit the bat pollination syndrome in the context of angiosperm phylogeny, with focus on the showy 'petaloid' organs associated with the syndrome; (4) to revisit the genetic-developmental basis of flower colour; (5) to raise evolutionary hypotheses of floral evolution associated with the bat pollination syndrome; and (6) to highlight some of the current frontiers of research on the origin and evolution of flowers and its impact on pollination syndrome studies in the 21st century., Conclusions: The inclusion of the concepts of analogy and convergence within the concept of syndromes will constitute a new agenda of inquiry that integrates floral biology, phylogenetic systematics and developmental biology. Phyllostomid and pteropodid bat pollination syndrome traits in eudicots and monocots represent cases of analogous and convergent evolution. Pollination syndromes are a multivariate concept intrinsically related to the understanding of flower organogenesis and evolution. The formulation of hypotheses of pollination syndromes must consider the phylogenetic levels of universality for both plant and animal taxa, flower development, genetics, homology and evolution, and a clear definition of evolutionary concepts, including analogy, convergence, homoplasy and parallelism., (© The Author(s) 2023. Published by Oxford University Press on behalf of the Annals of Botany Company. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Interspecific hybridization events are on the rise in natural systems due to climate change disrupting species barriers. Across taxa, microsatellites have long been the molecular markers of choice to identify admixed individuals. However, with the advent of high-throughput sequencing easing the generation of genome-wide datasets, incorrect reports of hybridization resulting from microsatellite technical artefacts have been uncovered in a growing number of taxa. In the marine zooplankton genus Calanus (Copepoda), whose species are used as climate change indicators, microsatellite markers have suggested hybridization between C. finmarchicus and C. glacialis, while other nuclear markers (InDels) never detected any admixed individuals, leaving the scientific community divided. Here, for the first time, we investigated the potential for hybridization among C. finmarchicus, C. glacialis, C. helgolandicus and C. hyperboreus using two large and independent SNP datasets. These were derived firstly from a protocol of target-capture applied to 179 individuals collected from 17 sites across the North Atlantic and Arctic Oceans, including sympatric areas, and second from published RNA sequences. All SNP-based analyses were congruent in showing that Calanus species are distinct and do not appear to hybridize. We then thoroughly re-assessed the microsatellites showing hybrids, with the support of published transcriptomes, and identified technical issues plaguing eight out of 10 microsatellites, including size homoplasy, paralogy, potential for null alleles and even two primer pairs targeting the same locus. Our study illustrates how deceptive microsatellites can be when applied to the investigation of hybridization., (© 2023 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

Despite the obstacles facing marine colonists, most lineages of aquatic organisms have colonized and diversified in freshwaters repeatedly. These transitions can trigger rapid morphological or physiological change and, on longer timescales, lead to increased rates of speciation and extinction. Diatoms are a lineage of ancestrally marine microalgae that have diversified throughout freshwater habitats worldwide. We generated a phylogenomic data set of genomes and transcriptomes for 59 diatom taxa to resolve freshwater transitions in one lineage, the Thalassiosirales. Although most parts of the species tree were consistently resolved with strong support, we had difficulties resolving a Paleocene radiation, which affected the placement of one freshwater lineage. This and other parts of the tree were characterized by high levels of gene tree discordance caused by incomplete lineage sorting and low phylogenetic signal. Despite differences in species trees inferred from concatenation versus summary methods and codons versus amino acids, traditional methods of ancestral state reconstruction supported six transitions into freshwaters, two of which led to subsequent species diversification. Evidence from gene trees, protein alignments, and diatom life history together suggest that habitat transitions were largely the product of homoplasy rather than hemiplasy, a condition where transitions occur on branches in gene trees not shared with the species tree. Nevertheless, we identified a set of putatively hemiplasious genes, many of which have been associated with shifts to low salinity, indicating that hemiplasy played a small but potentially important role in freshwater adaptation. Accounting for differences in evolutionary outcomes, in which some taxa became locked into freshwaters while others were able to return to the ocean or become salinity generalists, might help further distinguish different sources of adaptive mutation in freshwater diatoms., (© The Author(s) 2023. Published by Oxford University Press, on behalf of the Society of Systematic Biologists. All rights reserved. For permissions, please email: journals.permissions@oup.com.)

Retrotransposon insertion patterns facilitate a virtually homoplasy-free picture of phylogenetic history. Still, a few most likely random parallel insertions or deletions result in rare cases of homoplasy in primates. The following question arises: how frequent is retrotransposon homoplasy in other phylogenetic clades? Here, we derived genome insertion data of toothed whales to evaluate the extension of homoplasy in a representative laurasiatherian group. Among more than a thousand extracted and aligned retrotransposon loci, we detected 37 cases of precise parallel insertions in species that are separated by over more than 10 million years, a time frame which minimizes the effects of incomplete lineage sorting. We compared the phylogenetic signal of insertions with the flanking sequences of these loci to further exclude potential polymorphic loci derived by incomplete lineage sorting. We found that the phylogenetic signals of retrotransposon insertion patterns exhibiting true homoplasy differ from the signals of their flanking sequences. In toothed whales, precise parallel insertions account for around 0.18-0.29% of insertion cases, which is about 12.5 times the frequency of such insertions among Alu s in primates. We also detected five specific deletions of retrotransposons on various lineages of toothed whale evolution, a frequency of 0.003%, which is slightly higher than such occurrences in primates. Overall, the level of retrotransposon homoplasy in toothed whales is still marginal compared to the phylogenetic diagnostic retrotransposon presence/absence signal.

Selaginella is the largest and most taxonomically complex genus in lycophytes. The fact that over 750 species are currently treated in a single genus makes Selaginellales/Selaginellaceae unique in pteridophytes. Here we assembled a dataset of six existing and newly sampled plastid and nuclear loci with a total of 684 accessions (74% increase of the earlier largest sampling) representing ca . 300 species to infer a new phylogeny. The evolution of 10 morphological characters is studied in the new phylogenetic context. Our major results include: (1) the nuclear and plastid phylogenies are congruent with each other and combined analysis well resolved and strongly supported the relationships of all but two major clades; (2) the Sinensis group is resolved as sister to S. subg. Pulviniella with strong support in two of the three analyses; (3) most morphological characters are highly homoplasious but some characters alone or combinations of characters well define the major clades in the family; and (4) an infrafamilial classification of Selaginellaceae is proposed and the currently defined Selaginella s.l. is split into seven subfamilies (corresponding to the current six subgenera + the Sinensis group) and 19 genera (the major diagnosable clades) with nine new species-poor genera. We support the conservation of Selaginella with a new type, S. flabellata , to minimize nomenclatural instability. We provide a key to subfamilies and genera, images illustrating their morphology, their morphological and geographical synopses, a list of constituent species, and necessary new combinations. This new classification will hopefully facilitate communication, promote further studies, and help conservation., Competing Interests: We declare no conflicts of interest., (© 2023 Kunming Institute of Botany, Chinese Academy of Sciences. Publishing services by Elsevier B.V. on behalf of KeAi Communications Co., Ltd.)

Several studies have focused on the phylogenetic relationships within the Geoplaninae land planarians (Tricladida). In those studies, ancient phylogenetic relationships remained obscure. In this work, the phylogeny of Geoplaninae is assessed through three different datasets, namely morphological, molecular, and both datasets combined, i.e, a total evidence approach (TE). The data matrix consisted of six DNA regions, including a newly developed marker (DOM5), and a morphological matrix with 37 characters. The study produced the best-resolved hypothesis so far for the phylogeny of Geoplaninae, although ancient clades still remain elusive. The effect of the morphological data on the TE tree topology and clade support is seemingly negligible. The phylogenetic tree also suggests that most of the diagnostic morphological characters of the genera are homoplastic, while unambiguous unique synapomorphies can characterize some supra-generic informal groupings., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023. Published by Elsevier Inc.)

The term "homology" is persistently polysemous, defying the expectation that extensive scientific research should yield semantic stability. A common response has been to seek a unification of various prominent definitions. This paper proposes an alternative strategy, based on the insight that scientific concepts function as tools for research: When analyzing various conceptualizations of homology, we should preserve those distinguishing features that support particular research goals. We illustrate the fruitfulness of our strategy by application to two cases. First, we revisit Lankester's celebrated evolutionary reappraisal of homology and argue that his analysis has been distorted by assimilation to modern agendas. His "homogeny" does not mean the same thing as modern evolutionary "homology," and his "homoplasy" is no mere antonym. Instead, Lankester uses both new terms to pose a question that remains strikingly relevant-how do mechanistic and historical causes of morphological resemblance interact? Second, we examine the puzzle of avian digit homology, which exemplifies disciplinary differences in homology conceptualization and assessment. Recent progress has been fueled by the development of new tools within the relevant disciplines (paleontology and developmental biology) and especially by increasing interdisciplinary cooperation. Conceptual unification has played very little role in this work, which instead seeks concrete evolutionary scenarios that integrate all the available evidence. Together these cases indicate the complex relationship between concepts and other tools in homology research., (© 2023 The Authors. Journal of Morphology published by Wiley Periodicals LLC.)

The hemiphractid frog genus Stefania is one of the many ancient (near-) endemic lineages of vertebrates inhabiting the biodiverse Pantepui biogeographical region in the Guiana Shield Highlands of northern South America-the famous "Lost World" of Arthur Conan Doyle. Previous molecular analyses of the genus Stefania have indicated that species boundaries and phylogenetic relationships are often incongruent with morphological traits in that clade. A substantial number of "taxonomically cryptic" species, often microendemic, remain to be described. This is notably the case for an isolated population from the summit of Wei-Assipu-tepui, a small table-top mountain at the border between Guyana and Brazil. That population was previously known as Stefania sp. 6 and belongs to the S. riveroi clade. The new species is phylogenetically distinct, but phenotypically extremely similar to S. riveroi, a taxon found only on the summit of Yuruaní-tepui in Venezuela and recovered as sister to all the other known species in the S. riveroi clade. The new taxon is described based on morphology and osteology. Data about genetic divergences within the S. riveroi clade are provided. A new synapomorphy for the genus Stefania is proposed: the presence of a distal process on the third metacarpal. Amended definitions are offered for the three other species in the S. riveroi clade (S. ayangannae, S. coxi, S. riveroi). The new species should be listed as Critically Endangered according to IUCN criteria., (© 2023. The Author(s).)

The phylogeny of sea spiders has been debated for more than a century. Despite several molecular studies in the last twenty years, interfamilial relationships remain uncertain. In the present study, relationships within Pycnogonida are examined in the light of a new dataset composed of 160 mitochondrial genomes (including 152 new sequences) and 130 18S rRNA gene sequences (including 120 new sequences), from 141 sea spider morphospecies representing 26 genera and 9 families. Node congruence between mitochondrial and nuclear markers was analysed to identify the most reliable relationships. We also reanalysed a multilocus dataset previously published and showed that the high percentages of missing data make phylogenetic conclusions difficult and uncertain. Our results support the monophyly of most families currently accepted, except Callipallenidae and Nymphonidae, the monophyly of the superfamilies Ammotheoidea (Ammotheidae + Pallenopsidae), Nymphonoidea (Nymphonidae + Callipallenidae), Phoxichilidioidea (Phoxichilidiidae + Endeidae) and Colossendeoidea (Colossendeidae + Pycnogonidae + Rhynchothoracidae), and the sister-group relationship between Ammotheoidea and Phoxichilidioidea. We discuss the morphological evolution of sea spiders, identifying homoplastic characters and possible synapomorphies. We also discuss the palaeontological and phylogenetic arguments supporting either a radiation of sea spiders prior to Jurassic or a progressive diversification from Ordovician or Cambrian., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 Elsevier Inc. All rights reserved.)

Hominoids diverged from cercopithecoids during the Oligocene in Afro-Arabia, initially radiating in that continent and subsequently dispersing into Eurasia. From the Late Miocene onward, the geographic range of hominoids progressively shrank, except for hominins, which dispersed out of Africa during the Pleistocene. Although the overall picture of hominoid evolution is clear based on available fossil evidence, many uncertainties persist regarding the phylogeny and paleobiogeography of Miocene apes (nonhominin hominoids), owing to their sparse record, pervasive homoplasy, and the decimated current diversity of this group. We review Miocene ape systematics and evolution by focusing on the most parsimonious cladograms published during the last decade. First, we provide a historical account of the progress made in Miocene ape phylogeny and paleobiogeography, report an updated classification of Miocene apes, and provide a list of Miocene ape species-locality occurrences together with an analysis of their paleobiodiversity dynamics. Second, we discuss various critical issues of Miocene ape phylogeny and paleobiogeography (hylobatid and crown hominid origins, plus the relationships of Oreopithecus) in the light of the highly divergent results obtained from cladistic analyses of craniodental and postcranial characters separately. We conclude that cladistic efforts to disentangle Miocene ape phylogeny are potentially biased by a long-branch attraction problem caused by the numerous postcranial similarities shared between hylobatids and hominids-despite the increasingly held view that they are likely homoplastic to a large extent, as illustrated by Sivapithecus and Pierolapithecus-and further aggravated by abundant missing data owing to incomplete preservation. Finally, we argue that-besides the recovery of additional fossils, the retrieval of paleoproteomic data, and a better integration between cladistics and geometric morphometrics-Miocene ape phylogenetics should take advantage of total-evidence (tip-dating) Bayesian methods of phylogenetic inference combining morphologic, molecular, and chronostratigraphic data. This would hopefully help ascertain whether hylobatid divergence was more basal than currently supported., Competing Interests: Declaration of competing interest The authors declare no competing interests., (Copyright © 2022 Elsevier Ltd. All rights reserved.)

The avian palaeognath phylogeny has been recently revised significantly due to the advancement of genome-wide comparative analyses and provides the opportunity to trace the evolution of the microstructure and crystallography of modern dinosaur eggshells. Here, eggshells of all major clades of Palaeognathae (including extinct taxa) and selected eggshells of Neognathae and non-avian dinosaurs are analysed with electron backscatter diffraction. Our results show the detailed microstructures and crystallographies of (previously) loosely categorized ostrich-, rhea-, and tinamou-style morphotypes of palaeognath eggshells. All rhea-style eggshell appears homologous, while respective ostrich-style and tinamou-style morphotypes are best interpreted as homoplastic morphologies (independently acquired). Ancestral state reconstruction and parsimony analysis additionally show that rhea-style eggshell represents the ancestral state of palaeognath eggshells both in microstructure and crystallography. The ornithological and palaeontological implications of the current study are not only helpful for the understanding of evolution of modern and extinct dinosaur eggshells, but also aid other disciplines where palaeognath eggshells provide useful archive for comparative contrasts (e.g. palaeoenvironmental reconstructions, geochronology, and zooarchaeology)., Competing Interests: SC, MH, LL, NK, YL, DV No competing interests declared, (© 2023, Choi et al.)

Background: Citrus species are among the most important and widely consumed fruit trees in the world and are subjected to increasing global cultivation. Sweet orange (Citrus sinensis L. Osbeck) is one of 30 species of citrus which is cultivated in different regions of Iran. In this study, 80 trees of 13 sweet orange cultivars of Mazandaran province were studied for genetic diversity and fingerprinting by five short simple repeat (SSR) marker., Results: The studied cultivars showed a high degree of genetic variability with an average genetic polymorphism of 98.46%. Behshahr and Jadeh Ghadim2 genotypes had the highest and lowest values in Nei genetic diversity, number of effective alleles, and Shannon index, respectively. Based on k-means clustering, the studied genotypes were divided into two main different groups. The high magnitude of genetic similarity between replicates of different cultivars indicated a potential case of homonymy or synonymy. DAPC analysis showed genetic admixture among some of the cultivars. The heatmap plot illustrated the alleles involved in cultivar differentiation. The CAPs analysis of monomorphic alleles of SSR loci indicated that these alleles differ in their sequences which add up to the genetic variability of citrus germplasm., Conclusion: In general, SSR markers, due to their codominant nature and abundance in genome, are a good indicator for cultivar fingerprinting and hybrid prediction in orange cultivars. The present results showed the high diversity of sweet orange trees in different cultivars in the north of the country., (© 2022. The Author(s).)

Fungi show a high degree of morphological convergence. Regarded for a long time as an obstacle for phylogenetic studies, homoplasy has also been proposed as a source of information about underlying morphogenetic patterning mechanisms. The "local-activation and long-range inhibition principle" (LALIP), underlying the famous reaction-diffusion model proposed by Alan Turing in 1952, appears to be one of the universal phenomena that can explain the ontogenetic origin of seriate patterns in living organisms. Reproductive structures of fungi in the class Agaricomycetes show a highly periodic structure resulting in, for example, poroid, odontoid, lamellate or labyrinthic hymenophores. In this paper, we claim that self-organized patterns might underlie the basic ontogenetic processes of these structures. Simulations based on LALIP-driven models and covering a wide range of parameters show an absolute mutual correspondence with the morphospace explored by extant agaricomycetes. This could not only explain geometric particularities but could also account for the limited possibilities displayed by hymenial configurations, thus making homoplasy a direct consequence of the limited morphospace resulting from the proposed patterning dynamics., (© 2022. The Author(s), under exclusive licence to Springer-Verlag GmbH Germany, part of Springer Nature.)

We have detected two mutations in the spike protein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) at amino acid positions 1163 and 1167 that appeared independently in multiple transmission clusters and different genetic backgrounds. Furthermore, both mutations appeared together in a cluster of 1,627 sequences belonging to clade 20E. This cluster is characterized by 12 additional single nucleotide polymorphisms but no deletions. The available structural information on the S protein in the pre- and postfusion conformations predicts that both mutations confer rigidity, which could potentially decrease viral fitness. Accordingly, we observed reduced infectivity of this spike genotype relative to the ancestral 20E sequence in vitro , and the levels of viral RNA in nasopharyngeal swabs were not significantly higher. Furthermore, the mutations did not impact thermal stability or antibody neutralization by sera from vaccinated individuals but moderately reduce neutralization by convalescent-phase sera from the early stages of the pandemic. Despite multiple successful appearances of the two spike mutations during the first year of SARS-CoV-2 evolution, the genotype with both mutations was displaced upon the expansion of the 20I (Alpha) variant. The midterm fate of the genotype investigated was consistent with the lack of advantage observed in the clinical and experimental data. IMPORTANCE We observed repeated, independent emergence of mutations in the SARS-CoV-2 spike involving amino acids 1163 and 1167, within the HR2 functional motif. Conclusions derived from evolutionary and genomic diversity analysis suggest that the co-occurrence of both mutations might pose an advantage for the virus and therefore a threat to effective control of the epidemic. However, biological characterization, including in vitro experiments and analysis of clinical data, indicated no clear benefit in terms of stability or infectivity. In agreement with this, continuous epidemiological surveillance conducted months after the first observations revealed that both mutations did not successfully outcompete other variants and stopped circulating 9 months after their initial detection. Additionally, we evaluated the potential of both mutations to escape neutralizing antibodies, finding that the presence of these two mutations on their own is not likely to confer antibody escape. Our results provide an example of how newly emerged spike mutations can be assessed to better understand the risk posed by new variants and indicate that some spike mutations confer no clear advantage to the virus despite independently emerging multiple times and are eventually displaced by fitter variants.

Homoplasic SNPs are considered important signatures of strong (positive) selective pressure, and hence of adaptive evolution for clinically relevant traits such as antibiotic resistance and virulence. Here we present a new tool, SNPPar, for efficient detection and analysis of homoplasic SNPs from large whole genome sequencing datasets (>1000 isolates and/or >100 000 SNPs). SNPPar takes as input an SNP alignment, tree and annotated reference genome, and uses a combination of simple monophyly tests and ancestral state reconstruction (ASR, via TreeTime) to assign mutation events to branches and identify homoplasies. Mutations are annotated at the level of codon and gene, to facilitate analysis of convergent evolution. Testing on simulated data (120 Mycobacterium tuberculosis alignments representing local and global samples) showed SNPPar can detect homoplasic SNPs with very high specificity (zero false-positives in all tests) and high sensitivity (zero false-negatives in 89 % of tests). SNPPar analysis of three empirically sampled datasets ( Elizabethkingia anophelis , Burkholderia dolosa and M. tuberculosis ) produced results that were in concordance with previous studies, in terms of both individual homoplasies and evidence of convergence at the codon and gene levels. SNPPar analysis of a simulated alignment of ~64 000 genome-wide SNPs from 2000 M. tuberculosis genomes took ~23 min and ~2.6 GB of RAM to generate complete annotated results on a laptop. This analysis required ASR be conducted for only 1.25 % of SNPs, and the ASR step took ~23 s and 0.4 GB of RAM. SNPPar automates the detection and annotation of homoplasic SNPs efficiently and accurately from large SNP alignments. As demonstrated by the examples included here, this information can be readily used to explore the role of homoplasy in parallel and/or convergent evolution at the level of nucleotide, codon and/or gene.

Premise: The genus Gynoxys and relatives form a species-rich lineage of Andean shrubs and trees with low genetic distances within the sunflower subtribe Tussilaginineae. Previous molecular phylogenetic investigations of the Tussilaginineae have included few, if any, representatives of this Gynoxoid group or reconstructed ambiguous patterns of relationships for it., Methods: We sequenced complete plastid genomes of 21 species of the Gynoxoid group and related Tussilaginineae and conducted detailed comparisons of the phylogenetic relationships supported by the gene, intron, and intergenic spacer partitions of these genomes. We also evaluated the impact of manual, motif-based adjustments of automatic DNA sequence alignments on phylogenetic tree inference., Results: Our results indicate that the inclusion of all plastid genome partitions is needed to infer well-supported phylogenetic trees of the Gynoxoid group. Whole plastome-based tree inference suggests that the genera Gynoxys and Nordenstamia are polyphyletic and form the core clade of the Gynoxoid group. This clade is sister to a clade of Aequatorium and Paragynoxys and also includes some but not all representatives of Paracalia., Conclusions: The concatenation and combined analysis of all plastid genome partitions and the construction of manually-curated, motif-based DNA sequence alignments are found to be instrumental in the recovery of well-supported relationships of the Gynoxoid group. We demonstrate that the correct assessment of homology in genome-level plastid sequence data sets is crucial for subsequent phylogeny reconstruction and that the manual post-processing of multiple sequence alignments improves the reliability of such reconstructions amid low genetic distances between taxa., (© 2021 The Authors. American Journal of Botany published by Wiley Periodicals LLC on behalf of Botanical Society of America.)

The majority of animal phyla have species that can regenerate. Comparing regeneration across animals can reconstruct the molecular and cellular evolutionary history of this process. Recent studies have revealed some similarity in regeneration mechanisms, but rigorous comparative methods are needed to assess whether these resemblances are ancestral pathways (homology) or are the result of convergent evolution (homoplasy). This review aims to provide a framework for comparing regeneration across animals, focusing on gene regulatory networks (GRNs), which are substrates for assessing process homology. The homology of the wound-induced activation of Wnt signaling and of adult stem cells provides examples of ongoing studies of regeneration that enable comparisons in a GRN framework. Expanding the study of regeneration GRNs in currently studied species and broadening taxonomic sampling for these approaches will identify processes that are unifying principles of regeneration biology across animals. These insights are important both for evolutionary studies of regeneration and for human regenerative medicine.

The process of non-allelic gene conversion acts on homologous sequences during recombination, replacing parts of one with the other to make them uniform. Such concerted evolution is best described as paralogous ribosomal RNA gene unification that serves to preserve the essential house-keeping functions of the converted genes. Transposed elements (TE), especially Alu short interspersed elements (SINE) that have more than a million copies in primate genomes, are a significant source of homologous units and a verified target of gene conversion. The consequences of such a recombination-based process are diverse, including multiplications of functional TE internal binding domains and, for evolutionists, confusing divergent annotations of orthologous transposable elements in related species. We systematically extracted and compared 68,097 Alu insertions in various primates looking for potential events of TE gene conversion and discovered 98 clear cases of Alu - Alu gene conversion, including 64 cases for which the direction of conversion was identified (e.g., Alu S conversion to Alu Y). Gene conversion also does not necessarily affect the entire homologous sequence, and we detected 69 cases of partial gene conversion that resulted in virtual hybrids of two elements. Phylogenetic screening of gene-converted Alu s revealed three clear hotspots of the process in the ancestors of Catarrhini, Hominoidea, and gibbons. In general, our systematic screening of orthologous primate loci for gene-converted TEs provides a new strategy and view of a post-integrative process that changes the identities of such elements.

Bitter taste facilitates the detection of potentially harmful substances and is perceived via bitter taste receptors (TAS2Rs) expressed on the tongue and oral cavity in vertebrates. In primates, TAS2R16 specifically recognizes β-glucosides, which are important in cyanogenic plants' use of cyanide as a feeding deterrent. In this study, we performed cell-based functional assays for investigating the sensitivity of TAS2R16 to β-glucosides in three species of bamboo lemurs ( Prolemur simus, Hapalemur aureus and H. griseus ), which primarily consume high-cyanide bamboo. TAS2R16 receptors from bamboo lemurs had lower sensitivity to β-glucosides, including cyanogenic glucosides, than that of the closely related ring-tailed lemur ( Lemur catta ). Ancestral reconstructions of TAS2R16 for the bamboo-lemur last common ancestor (LCA) and that of the Hapalemur LCA showed an intermediate sensitivity to β-glucosides between that of the ring-tailed lemurs and bamboo lemurs. Mutagenetic analyses revealed that P. simus and H. griseus had separate species - specific substitutions that led to reduced sensitivity. These results indicate that low sensitivity to β-glucosides at the cellular level-a potentially adaptive trait for feeding on cyanogenic bamboo-evolved independently after the Prolemur - Hapalemur split in each species.

The superfamily Alveolinoidea is a member of the Order Miliolida, and comprises three main families, the Alveolinidae, the Fabulariidae and the Rhapydioninidae. They are examples of Larger benthic foraminifera (LBF), which are single-celled organisms with specific characteristic endoskeletons. Alveolinoids are found globally from the Cretaceous to the present day, and are important biostratigraphic index fossils in shallow-marine carbonates. They are often associated with hydrocarbon reservoirs, and exhibit provincialism with characteristic genera often confined to one of the American, Tethyan or Indo-Pacific provinces. Previously, the systematic study of the global interrelationship between the various alveolinoid lineages has not been possible because of the absence of biostratigraphic correlation between the geographically scattered assemblages, and the scarcity of described material from the Indo-Pacific province. Here we use the literature and new material from the Americas, the French Alps, Iran, Tibet, India and South East Asia, coupled with the use of the planktonic foraminiferal zonal (PZ) correlation scheme to propose a comprehensive, global, systematic analysis of the biostratigraphic, phylogenetic and paleogeographic evolution of the alveolinoids. The alveolinoids originated in the Cretaceous in the Tethyan province. During a global sea-level low stand, a westward migration of some alveolinoids species to the Americas occurred, a behaviour previously reported in contemporaneous orbitolinid LBF. After the Cretaceous/Palaeogene (K-P) event, which saw the extinction of all Cretaceous alveolinoids, rare new forms of alveolinoids evolved again, first in the Americas and later independently in Tethys. As was found in previous studies of rotalid LBF, sea-level low stands in the Paleocene also allowed some alveolinoid forms to migrate, but this time in an eastward direction from the Americas to Tethys, and from Tethys on to the Indo-Pacific province. Alveolinoids still exist today ( Borelis and Alveolinella ), the former of which is cosmopolitan, while the latter is restricted to the Indo-Pacific province. Throughout their phylogenetic history, alveolinoids characteristically exhibit convergent evolution, with the repeated re-occurrence of certain morphological features. Understanding this propensity to homoplasy is essential in understanding and constructing the phylogenetic relationships within the alveolinoid superfamily., Competing Interests: The authors declare no conflicts of interest in connection to this article., (© 2021 The Authors.)

Homoplasy is a sort of noise in phylogenetic reconstructions, due to the accumulation of backmutations, convergent evolution and horizontal gene transfer (HGT), which is considered the major trigger of homoplasy in microorganism for its massive presence. It is also known that homoplasy increases with the complexity of the tree with both real and simulated data. In this paper, we analyzed the variation of homoplasy with the two widely used taxonomic markers ITS and LSU in four taxonomic models characterized by differences in the intra-specific distances. An algorithm (HomoDist) was developed to analyze the homoplasy index (HI) variation upon addition of a single element (strain or species) in increasing distance from a starting element. This algorithm allows to follow changes of the consistency index (CI), complementary to the HI, with the increase of the number of taxa and with the increase of the distance among elements. Results show that homoplasy increases-as expected-with the number of taxa, but also as a function of the overall distance among species, often with an almost linear relationship between distance and HI. No HI change was observed in trees with few taxa spanning through short distances, indicating that this noise is not prohibitive in this context, although the analysis of the ratio between HI and distance can be recommended as a criterion for tree acceptance. The absence of large changes of the HI within the species, and its increase when new species are added by HomoDist, suggest that homoplasy variation can be used as an auxiliary test in distance-based species delimitation with any type of marker.

Osteostraci and Galeaspida are stem-gnathostomes, occupying a key phylogenetic position for resolving the nature of the jawless ancestor from which jawed vertebrates evolved more than 400 million years ago. Both groups are characterized by the presence of rigid headshields that share a number of common morphological traits, in some cases hindering the resolution of their interrelationships and the exact nature of their affinities with jawed vertebrates. Here, we explore the morphological and functional diversity of osteostracan and galeaspid headshields using geometric morphometrics and computational fluid dynamics to constrain the factors that promoted the evolution of their similar morphologies and informing on the ecological scenario under which jawed vertebrates emerged. Phylomorphospace, Mantel analysis and Stayton metrics demonstrate a high degree of homoplasy. Computational fluid dynamics reveals similar hydrodynamic performance among morphologically convergent species, indicating the independent acquisition of the same morphofunctional traits and, potentially, equivalent lifestyles. These results confirm that a number of the characters typically used to infer the evolutionary relationships among galeaspids, osteostracans and jawed vertebrates are convergent in nature, potentially obscuring understanding of the assembly of the gnathostome bodyplan. Ultimately, our results reveal that while the jawless relatives of the earliest jawed vertebrates were ecologically diverse, widespread convergence on the same hydrodynamic adaptations suggests they had reached the limits of their potential ecological diversity-overcome by jawed vertebrates and their later innovations.

The fact that angiosperms are so species-rich and ecologically diverse - Darwin's second abominable mystery - could be explained by their ability to 'reinvent' themselves by evolving functional traits repeatedly over time, space and taxonomic clades. These trait innovations may facilitate adaptation and increase diversification rates. In this article, I quantify this 'trait flexibility' by reviewing the literature on trait transition rates and trait-dependent diversification rates in angiosperms and their extant sister clade, acrogymnosperms. I show that angiosperms indeed evolved elevated trait transition and trait-dependent diversification rates compared to gymnosperms, and rates are highest within species-rich angiosperm orders (e.g. Fabales, Lamiales). The (genetic) ability of certain angiosperm lineages to repeatedly evolve adaptive traits may have facilitated sustained high net diversification resulting from numerous episodic radiations., (© 2019 The Author New Phytologist © 2019 New Phytologist Foundation.)

More human deaths have been attributable to Mycobacterium tuberculosis than any other pathogen, and the epidemic is sustained by ongoing transmission. Various typing schemes have been developed to identify strain-specific differences and track transmission dynamics in affected communities, with recent introduction of whole genome sequencing providing the most accurate assessment. Mycobacterial interspersed repetitive unit (MIRU) typing is a family of variable number tandem repeat schemes that have been widely used to study the molecular epidemiology of M. tuberculosis . MIRU typing was used in most well-resourced settings to perform routine molecular epidemiology. Instances of MIRU homoplasy have been observed in comparison with sequence-based phylogenies, limiting its discriminatory value. A fundamental question is whether the observed homoplasy arises purely through stochastic processes, or whether there is evidence of natural selection. We compared repeat numbers at 24 MIRU loci with a whole genome sequence-based phylogeny of 245 isolates representing three modern M. tuberculosis lineages. This analysis demonstrated extensive homoplasy of repeat numbers, but did not detect any evidence of natural selection of repeat numbers, at least since the ancestral branching of the three modern lineages of M. tuberculosis . In addition, we observed good sensitivity but poor specificity and positive predictive values of MIRU-24 to detect clusters of recent transmission, as defined by whole-genome single nucleotide polymorphism analysis. These findings provide mechanistic insight, and support a transition away from VNTR-based typing toward sequence-based typing schemes for both research and public health purposes., (Copyright © 2020 Outhred, Gurjav, Jelfs, McCallum, Wang, Hill-Cawthorne, Marais and Sintchenko.)

The systematics of Hypotricha is one of the most puzzling problems in ciliate biology, having spanned numerous conflicting hypotheses with unstable relationships at various levels in molecular trees, for which the constant addition of newly discovered species has only increased the confusion. The hypotrichs comprise a remarkable morphologically diversified group of ciliates, and the phylogenetic potential of morphological traits is generally recognized. However, such characters were rarely used in phylogenetic reconstructions, and congruence with molecular data never assessed from simultaneous analyses. To properly reconciliate morphological and molecular information, maximum-likelihood and parsimony analyses of 79 morphological characters and 18S rDNA sequences were performed for 130 ingroup terminals, broadly sampled to represent the known hypotrich diversity. As result, well-supported and relatively stable clades were recovered, based on which the redefined Hypotricha comprises at least six higher taxa: The "arcuseriids", Holostichida, Parabirojimida, and the "amphisiellids", plus the two large clades Kentrurostylida nov. tax. (Hispidotergida nov. tax. and Simplicitergida nov. tax.) and Diatirostomata nov. tax. ("bistichellids", "kahliellids", Gonostomatida and Dorsomarginalia [Postoralida nov. tax. and Uroleptida]). Each taxon was circumscribed by synapomorphies, of which most were homoplastic, as the natural history of hypotrichs is portrayed by an outstanding quantity of convergences and reversions., (Copyright © 2020 Elsevier GmbH. All rights reserved.)

Many antibacterial drugs have multiple mechanisms of resistance, which are often represented simultaneously by a mixture of resistance mutations (some more frequent than others) in a clinical population. This presents a challenge for Genome-Wide Association Studies (GWAS) methods, making it difficult to detect less prevalent resistance mechanisms purely through (weak) statistical associations. Homoplasy, or the occurrence of multiple independent mutations at the same site, is often observed with drug resistance mutations and can be a strong indicator of positive selection. However, traditional GWAS methods, such as those based on allele counting or linear regression, are not designed to take homoplasy into account. In this article, we present a new method, called ECAT (for Evolutionary Cluster-based Association Test), that extends traditional regression-based GWAS methods with the ability to take advantage of homoplasy. This is achieved through a preprocessing step which identifies hypervariable regions in the genome exhibiting statistically significant clusters of distinct evolutionary changes, to which association testing by a linear mixed model (LMM) is applied using GEMMA (a well-established LMM-based GWAS tool). Thus, the approach can be viewed as extending GEMMA from the usual site- or gene-level analysis to focusing on clustered regions of mutations. This approach was evaluated on a large collection of more than 600 clinical isolates of multidrug-resistant (MDR) Mycobacterium tuberculosis from Lima, Peru. We show that ECAT does a better job of detecting known resistance mutations for several antitubercular drugs (including less prevalent mutations with weaker associations), compared with (site- or gene-based) GEMMA, as representative of existing GWAS methods. The power of the multiphase approach in ECAT comes from focusing association testing on the hypervariable regions of the genome, which reduces complexity in the model and increases statistical power., Competing Interests: Declaration of Conflicting Interests:The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article., (© The Author(s) 2020.)

Background: Changes in gene regulation are widely recognized as an important driver of adaptive phenotypic evolution. However, the specific molecular mechanisms that underpin such changes are still poorly understood. Chromatin state plays an essential role in gene regulation, by influencing the accessibility of coding loci to the transcriptional machinery. Changes in the function of chromatin remodellers are therefore strong candidates to drive changes in gene expression associated with phenotypic adaptation. Here, we identify amino acid homoplasies in the chromatin remodeller CHD9, shared between the extinct marsupial thylacine and eutherian wolf which show remarkable skull convergence. CHD9 is involved in osteogenesis, though its role in the process is still poorly understood. We examine whether CHD9 is able to regulate the expression of osteogenic target genes and examine the function of a key substitution in the CHD9 DNA binding domain., Results: We examined whether CHD9 was able to upregulate its osteogenic target genes, RUNX2, Osteocalcin (OC) and ALP in HEK293T cells. We found that overexpression of CHD9 upregulated RUNX2, the master regulator of osteoblast cell fate, but not the downstream genes OC or ALP, supporting the idea that CHD9 regulates osteogenic progenitors rather than terminal osteoblasts. We also found that the evolutionary substitution in the CHD9 DNA binding domain does not alter protein secondary structure, but was able to drive a small but insignificant increase in RUNX2 activation. Finally, CHD9 was unable to activate an episomal RUNX2 promoter-reporter construct, suggesting that CHD9 requires the full chromatin complement for its function., Conclusions: We provide new evidence to the role of CHD9 in osteogenic differentiation through its newly observed ability to upregulate the expression of RUNX2. Though we were unable to identify significant functional consequences of the evolutionary substitution in HEK293T cells, our study provides important steps forward in the functional investigation of protein homoplasy and its role in developmental processes. Mutations in coding genes may be a mechanism for driving adaptive changes in gene expression, and their validation is essential towards determining the functional consequences of evolutionary homoplasy.

Accurate reconstructions of phylogeny are essential for studying the evolution of a clade, and morphological characters are necessarily used for the reconstruction of the relationships of fossil organisms. However, variation in their evolutionary modes (for example rate variation and character non-independence) not accounted for in analyses may be leading to unreliable phylogenies. A recent study suggested that phylogenetic analyses of mammals may be suffering from a dominance of dental characters, which were shown to have lower phylogenetic signal than osteological characters and produced phylogenies less congruent with molecularly-derived benchmarks. Here we build on this previous work by testing five additional morphological partitions for phylogenetic signal and examining what aspects of dental and other character evolution may be affecting this, by fitting models of discrete character evolution to phylogenies inferred and time calibrated using molecular data. Results indicate that the phylogenetic signal of discrete characters correlate most strongly with rates of evolution, with increased rates driving increased homoplasy. In a dataset covering all Mammalia, dental characters have higher rates of evolution than other partitions. They do not, however, fit a model of independent character evolution any worse than other regions. Primates and marsupials show different patterns to other mammal clades, with dental characters evolving at slower rates and being more heavily integrated (less independent). While the dominance of dental characters in analyses of mammals could be leading to inaccurate phylogenies, the issue is not unique to dental characters and the results are not consistent across datasets. Molecular benchmarks (being entirely independent of the character data) provide a framework for examining each dataset individually to assess the evolution of the characters used., Competing Interests: The authors declare there are no competing interests., (©2020 Brocklehurst and Benevento.)

The dromaeosaurid theropod Halszkaraptor escuilliei is characterized by several unusual features absent in other paravians, part of which has been interpreted as diagnostic of a novel lineage adapted to a semiaquatic ecology. Recently, these evolutionary and ecological interpretations have been challenged, and Halszkaraptor has been claimed to be a transitional form between non-dromaeosaurid maniraptoriforms and other dromaeosaurids: following that reevaluation, its peculiar body plan would represent the retention of several maniraptoran plesiomorphies, lost among other dromaeosaurids, and not an adaptation to a novel ecology. This alternative scenario is here carefully investigated and tested. It is shown that most statements supporting this scenario are based on misinterpretation of anatomical traits and bibliography. Once these statements have been corrected, character state transition optimization over a well-supported phylogenetic framework indicates that the large majority of the peculiar features of the Halszkaraptor lineage are derived novelties acquired by the latter after its divergence from the last ancestor shared with eudromaeosaurs, and thus are not maniraptoriform plesiomorphies. At least seven novelties of the Halszkaraptor lineage are convergently acquired with spinosaurids, and are integrated in semiaquatic adaptations: one of these is reported here for the first time. The amount of morphological divergence of Halszkaraptorinae from the ancestral dromaeosaurid condition is comparable to those of Microraptorinae and Velociraptorinae. Among extant taxa, the sawbills (Mergini, Anseriformes) show the closest ecomorphological similarity with the peculiar body plan inferred for Halszkaraptor . The halszkaraptorine bauplan is thus confirmed as a derived amphibious specialization, and does not represent a "transitional" stage along the evolution of dromaeosaurids., Competing Interests: The author declares that he has no competing interests., (© 2020 Cau.)

Terrestrial artiodactyls (even-toed ungulates) inhabit some of the world's most extreme environments, including arid deserts and high elevations. As medium-to-large-bodied mammals, artiodactyls have a suite of specialized physiologies to facilitate occupation of regions unavailable to other large mammals. One such physiology is selective brain cooling, wherein reduction of brain temperature below core body temperature has been demonstrated to reduce evaporative water loss. This physiology is enabled by an arterial heat-exchanger called the carotid rete. The ubiquity of the carotid rete throughout the clade, as well as its evolutionary history, is currently uninvestigated. Here, I use osteological correlates to survey clade-wide presence and morphology of the carotid rete, prior to conducting a preliminary evolutionary analysis. Nearly all living artiodactyls possess a carotid rete and are capable of selective brain cooling; however, major arteries supplying the rete are derived from different embryonic aortic arches on a suborder-specific basis. Ancestral character estimation infers this pattern of variation to be the result of independent evolutionary processes, suggesting carotid rete homoplasy arising via parallelism. This is a surprising finding given the role this structure plays in driving a physiology that has been implicated in mitigating artiodactylan responses to extreme environmental conditions. Future studies should incorporate extinct species represented in the fossil record to better parse between parallel and convergent mechanisms, as well as to better understand the relationship between the carotid rete, selective brain cooling, and survivorship of climate perturbation. Anat Rec, 2018. © 2018 Wiley Periodicals, Inc. Anat Rec, 303:308-317, 2020. © 2018 American Association for Anatomy., (© 2018 Wiley Periodicals, Inc.)

Eugenia has a pantropical distribution and comprises ca. 1000 species found mostly in the Neotropics. Recent DNA based phylogenies show that unusual flower morphology of 'eugenioid' collections, e.g. fused calices that open by tearing, consistently emerged within Eugenia. These results emphasize a demand to revaluate flower morphology in a phylogenetic context within the genus. A reassessment of calyx fusion in Eugenia and traditionally related genera is here focused on clarification of the systematic relevance of this apparently recurrent characteristic. Twenty-four Eugenia species with some level of calyx fusion in the bud were newly used (one nuclear and four plastid markers) in conjunction with a representative sample of previously sequenced species to recover a time-calibrated Eugenia phylogeny of 86 accessions. Development of the fused calyx was analysed using scanning electron microscopy, differing patterns were re-coded and subsequently phylogenetic character reconstruction was performed. Eugenia was recovered as monophyletic including the traditionally segregated genera Calycorectes and Catinga. Ancestral character reconstruction uncovered free calyx lobes as the ancestral condition. Five development patterns leading to calyx fusion are reported in Eugenia including species with apparently six petals, which contrast with the standard tetramerous flowers. This condition is interpreted as the petaloid pattern, where two external fused calyx lobes cover the bud while two internal calyx lobes are free and petaloid. The fused calyx condition is homoplastic and evolved independently, several times in Eugenia, as did the different development patterns. Data presented here show that systematic incongruence resulting from multiple, independent origins of the fused calyx in Eugenia is further aggravated by an inability to distinguish parallelism and convergence within the recovered patterns., (Copyright © 2019 Elsevier Inc. All rights reserved.)

Understanding the evolution of biodiversity on Earth is a central aim in biology. Currently, various disciplines of science contribute to unravel evolution at all levels of life, from individual organisms to species and higher ranks, using different approaches and specific terminologies. The search for common origin, traditionally called homology, is a connecting paradigm of all studies related to evolution. However, it is not always sufficiently taken into account that defining homology depends on the hierarchical level studied (organism, population, and species), which can cause confusion. Therefore, we propose a framework to define homologies making use of existing terms, which refer to homology in different fields, but restricting them to an unambiguous meaning and a particular hierarchical level. We propose to use the overarching term "homology" only when "morphological homology," "vertical gene transfer," and "phylogenetic homology" are confirmed. Consequently, neither phylogenetic nor morphological homology is equal to homology. This article is intended for readers with different research backgrounds. We challenge their traditional approaches, inviting them to consider the proposed framework and offering them a new perspective for their own research., (© The Author(s) 2019. Published by Oxford University Press on behalf of the Society of Systematic Biologists.)

Accurate inferences of convergence require that the appropriate tree topology be used. If there is a mismatch between the tree a trait has evolved along and the tree used for analysis, then false inferences of convergence ('hemiplasy') can occur. To avoid problems of hemiplasy when there are high levels of gene tree discordance with the species tree, researchers have begun to construct tree topologies from individual loci. However, due to intralocus recombination, even locus-specific trees may contain multiple topologies within them. This implies that the use of individual tree topologies discordant with the species tree can still lead to incorrect inferences about molecular convergence. Here, we examine the frequency with which single exons and single protein-coding genes contain multiple underlying tree topologies, in primates and Drosophila, and quantify the effects of hemiplasy when using trees inferred from individual loci. In both clades, we find that there are most often multiple diagnosable topologies within single exons and whole genes, with 91% of Drosophila protein-coding genes containing multiple topologies. Because of this underlying topological heterogeneity, even using trees inferred from individual protein-coding genes results in 25% and 38% of substitutions falsely labelled as convergent in primates and Drosophila, respectively. While constructing local trees can reduce the problem of hemiplasy, our results suggest that it will be difficult to completely avoid false inferences of convergence. We conclude by suggesting several ways forward in the analysis of convergent evolution, for both molecular and morphological characters. This article is part of the theme issue 'Convergent evolution in the genomics era: new insights and directions'.

Bayesian analysis of morphological data is becoming increasingly popular mainly (but not only) because it allows for time-calibrated phylogenetic inference using relaxed morphological clocks and tip dating whenever fossils are available. As with molecular data, recent studies have shown that modeling among-character rate variation (ACRV) in morphological matrices greatly improves phylogenetic inference. In a likelihood framework this may be accomplished, for instance, by employing a hidden Markov model to assign characters to rate categories drawn from a (discretized) $\Gamma$ distribution and/or by partitioning data sets according to rate heterogeneity and estimating per-partition branch lengths, conditioned on a single topology. While the first approach is available in many phylogenetic analysis software, there is still no clear consensus on how to partition data, except perhaps in the simplest cases (e.g., "by codon" partitioning of coding sequences). Additionally, there is a trade-off between improvement in likelihood scores and the number of free parameters in the analysis, which rises quickly with the number of partitions. This trade-off may be dealt with by employing statistics that penalize overfitting of complex models, such as Akaike or Bayesian information criteria, or the more recently introduced stepping-stone method for marginal likelihood approximation. We applied the latter to three distinct matrices of discrete morphological data and demonstrated that sorting characters by homoplasy scores (obtained from implied weighting parsimony analysis) outperformed other partitioning strategies (anatomically-based and PartitionFinder2). The method was in fact so efficient in segregating characters by rates of evolution that no within-partition ACRV modeling was necessary, while among-partition rate variation was adequately accommodated by rate multipliers. We conclude that partitioning by homoplasy is a powerful and easy-to-implement strategy to address ACRV in complex data sets. We provide some guidelines focusing on morphological matrices, although this approach may be also applicable to molecular data sets., (© The Author(s) 2019. Published by Oxford University Press, on behalf of the Society of Systematic Biologists. All rights reserved. For permissions, please email: journals.permissions@oup.com.)

Background and Objectives: Interrogation of the genomic relations between Iranian Mycoplasma agalactiae vaccine strains of Taliqan, Lorestan and Shiraz., Materials and Methods: Two MLVA (covering VNTR loci of 5, 9, 17 and 19) and MLST (comprising dnaA, gltX, gyrB, C, tufA genes) genotyping systems plus nucleotide structure analysis of P80 gene, was conducted., Results: The shared MLVA pattern represented by the three strains differed to that of the Mag PG2 laboratory strain, only at locus VNTR19 where the PG2 genome hold a 3 bp longer stretch. In MLST analysis, at dnaA, gltX, gyrB, metS and tufA loci, the three strains displayed alleles 1, 21, 2, 2 and 1, respectively. At, gltX locus a new allele (21) was detected where a new sequence type (ST33) was identified. Besides, the trio strains hold an identical nucleotide structure in their ma-mp81 gene., Conclusion: In explanation, lack of efficient disease control measures, has possibly contributed in evolution of a clone or a few clones that gradually overwhelmed the population over the time. Besides, the similarity between the Iranian and the PG2 strains, might be due to homoplasy or farming exercises such as animal importation. Inclusion of further local isolates in next studies will help to assess these assumptions.

A homoplasy is a nucleotide identity resulting from a process other than inheritance from a common ancestor. Importantly, by distorting the ancestral relationships between nucleotide sequences, homoplasies can change the structure of the phylogeny. Homoplasies can emerge naturally, especially under high selection pressures and/or high mutation rates, or be created during the generation and processing of sequencing data. Identification of homoplasies is critical, both to understand their influence on the analyses of phylogenetic data and to allow an investigation into how they arose. Here we present HomoplasyFinder, a java application that can be used as a stand-a-lone tool or within the statistical programming environment R. Within R and Java, HomoplasyFinder is shown to be able to automatically, and quickly, identify any homoplasies present in simulated and real phylogenetic data. HomoplasyFinder can easily be incorporated into existing analysis pipelines, either within or outside of R, allowing the user to quickly identify homoplasies to inform downstream analyses and interpretation.

Premise of the Study: Simple sequence repeat (SSR) or microsatellite markers have been used in a broad range of studies mostly scoring alleles on the basis of amplicon size as a proxy for the number of repeat units of an SSR motif. However, additional sources of variation within the SSR or in the flanking regions have largely remained undetected., Methods: In this study, we implemented a next-generation sequencing-based genotyping approach in a newly characterized set of 18 nuclear SSR markers for the carob tree, Ceratonia siliqua . Our aim was to evaluate the effect of three different methods of scoring molecular variation present within microsatellite markers on the genetic diversity and structure results., Results: The analysis of the sequences of 77 multilocus genotypes from four populations revealed SSR variation and additional sources of polymorphism in 87% of the loci analyzed (42 single-nucleotide polymorphisms and five insertion/deletion polymorphisms), as well as divergent paralog copies in two loci. Ignoring sequence variation under standard amplicon size genotyping resulted in incorrect identification of 69% of the alleles, with important effects on the genetic diversity and structure estimates., Discussion: Next-generation sequencing allows the detection and scoring of SSRs, single-nucleotide polymorphisms, and insertion/deletion polymorphisms to increase the resolution of population genetic studies.

Convergent evolution-the appearance of the same character state in apparently unrelated organisms-is often inferred when a trait is incongruent with the species tree. However, trait incongruence can also arise from changes that occur on discordant gene trees, a process referred to as hemiplasy. Hemiplasy is rarely taken into account in studies of convergent evolution, despite the fact that phylogenomic studies have revealed rampant discordance. Here, we study the relative probabilities of homoplasy (including convergence and reversal) and hemiplasy for an incongruent trait. We derive expressions for the probabilities of the two events, showing that they depend on many of the same parameters. We find that hemiplasy is as likely-or more likely-than homoplasy for a wide range of conditions, even when levels of discordance are low. We also present a method to calculate the ratio of these two probabilities (the "hemiplasy risk factor") along the branches of a phylogeny of arbitrary length. Such calculations can be applied to any tree to identify when and where incongruent traits may be due to hemiplasy., Competing Interests: The authors declare no conflict of interest., (Copyright © 2018 the Author(s). Published by PNAS.)

Background: Members of the Bacillus genus have been isolated from a variety of environments. However, the relationship between potential metabolism and the niche from which bacteria of this genus have been isolated has not been extensively studied. The existence of a monophyletic aquatic Bacillus group, composed of members isolated from both marine and fresh water has been proposed. Here, we present a phylogenetic/phylogenomic analysis to investigate the potential relationship between the environment from which group members have been isolated and their evolutionary origin. We also carried out hierarchical clustering based on functional content to test for potential environmental effects on the genetic content of these bacteria., Results: The phylogenetic reconstruction showed that Bacillus strains classified as aquatic have evolutionary origins in different lineages. Although we observed the presence of a clade consisting exclusively of aquatic Bacillus, it is not comprised of the same strains previously reported. In contrast to phylogeny, clustering based on the functional categories of the encoded proteomes resulted in groups more compatible with the environments from which the organisms were isolated. This evidence suggests a detectable environmental influence on bacterial genetic content, despite their different evolutionary origins., Conclusion: Our results suggest that aquatic Bacillus species have polyphyletic origins, but exhibit convergence at the gene content level.

Introgressive hybridization between Micropterus dolomieu and Micropterus salmoides was assessed in their invaded South African range using nine microsatellite markers and two mtDNA gene regions. Although M. dolomieu and M. salmoides are distantly related, indicated by the large uncorrected pairwise distances observed between the two species, mitochondrial introgression and unidirectional admixture was detected., (© 2018 The Fisheries Society of the British Isles.)

Simultaneously analysing morphological, molecular and stratigraphic data suggests a potential resolution to a major remaining inconsistency in crocodylian evolution. The ancient, long-snouted thoracosaurs have always been placed near the Indian gharial Gavialis , but their antiquity ( ca 72 Ma) is highly incongruous with genomic evidence for the young age of the Gavialis lineage ( ca 40 Ma). We reconcile this contradiction with an updated morphological dataset and novel analysis, and demonstrate that thoracosaurs are an ancient iteration of long-snouted stem crocodylians unrelated to modern gharials. The extensive similarities between thoracosaurs and Gavialis are shown to be an almost 'perfect storm' of homoplasy, combining convergent adaptions to fish-eating, as well resemblances between genuinely primitive traits (thoracosaurs) and atavisms ( Gavialis ). Phylogenetic methods that ignore stratigraphy (parsimony and undated Bayesian methods) are unable to tease apart these similarities and invariably unite thoracosaurs and Gavialis. However, tip-dated Bayesian approaches additionally consider the large temporal gap separating ancient (thoracosaurs) and modern ( Gavialis ) iterations of similar long-snouted crocodyliforms. These analyses robustly favour a phylogeny which places thoracosaurs basal to crocodylians, far removed from modern gharials, which accordingly are a very young radiation. This phylogenetic uncoupling of ancient and modern gharial-like crocs is more consistent with molecular clock divergence estimates, and also the bulk of the crocodylian fossil record (e.g. all unequivocal gharial fossils are very young). Provided that the priors and models attribute appropriate relative weights to the morphological and stratigraphic signals-an issue that requires investigation-tip-dating approaches are potentially better able to detect homoplasy and improve inferences about phylogenetic relationships, character evolution and divergence dates., (© 2018 The Author(s).)

After decades of molecular phylogenetic studies, the deep phylogeny of gymnosperms has not been resolved, and the phylogenetic placement of Gnetales remains one of the most controversial issues in seed plant evolution. To resolve the deep phylogeny of seed plants and to address the sources of phylogenetic conflict, we conducted a phylotranscriptomic study with a sampling of all 13 families of gymnosperms and main lineages of angiosperms. Multiple datasets containing up to 1 296 042 sites across 1308 loci were analysed, using concatenation and coalescence approaches. Our study generated a consistent and well-resolved phylogeny of seed plants, which places Gnetales as sister to Pinaceae and thus supports the Gnepine hypothesis. Cycads plus Ginkgo is sister to the remaining gymnosperms. We also found that Gnetales and angiosperms have similar molecular evolutionary rates, which are much higher than those of other gymnosperms. This implies that Gnetales and angiosperms might have experienced similar selective pressures in evolutionary histories. Convergent molecular evolution or homoplasy is partially responsible for the phylogenetic conflicts in seed plants. Our study provides a robustly reconstructed backbone phylogeny that is important for future molecular and morphological studies of seed plants, in particular gymnosperms, in the light of evolution., (© 2018 The Author(s).)

This paper discusses the problem of whether creating a matrix with all the character state combinations that have a fixed number of steps (or extra steps) on a given tree T, produces the same tree T when analyzed with maximum parsimony or maximum likelihood. Exhaustive enumeration of cases up to 20 taxa for binary characters, and up to 12 taxa for 4-state characters, shows that the same tree is recovered (as unique most likely or most parsimonious tree) as long as the number of extra steps is within 1/4 of the number of taxa. This dependence, 1/4 of the number of taxa, is discussed with a general argumentation, in terms of the spread of the character changes on the tree used to select character state distributions. The present finding allows creating matrices which have as much homoplasy as possible for the most parsimonious or likely tree to be predictable, and examination of these matrices with hill-climbing search algorithms provides additional evidence on the (lack of a) necessary relationship between homoplasy and the ability of search methods to find optimal trees., (Copyright © 2018 Elsevier Inc. All rights reserved.)