Inferring phylogeny despite incomplete lineage sorting

Shared single copy genes are generally reliable for inferring phylogenetic relationships among polyploid taxa

Polyploidy, or whole-genome duplication, is expected to confound the inference of species trees with phylogenetic methods for two reasons. First, the presence of retained duplicated genes requires the reconciliation of the inferred gene trees to a proposed species tree. Second, even if the analyses are restricted to shared single copy genes, the occurrence of reciprocal gene loss, where the surviving genes in different species are paralogs from the polyploidy rather than orthologs, will mean that such genes will not have evolved under the corresponding species tree and may not produce gene trees that allow inference of that species tree. Here we analyze three different ancient polyploidy events, using synteny-based inferences of orthology and paralogy to infer gene trees from nearly 17,000 sets of homologous genes. We find that the simple use of single copy genes from polyploid organisms provides reasonably robust phylogenetic signals, despite the presence of reciprocal gene losses. Such gene trees are also most often in accord with the inferred species relationships inferred from maximum likelihood models of gene loss after polyploidy: a completely distinct phylogenetic signal present in these genomes. As seen in other studies, however, we find that methods for inferring phylogenetic confidence yield high support values even in cases where the underlying data suggest meaningful conflict in the phylogenetic signals.

Divergence and reticulation in the Mexican white oaks: ecological and phylogenomic evidence on species limits and phylogenetic networks in the Quercus laeta complex (Fagaceae)

BACKGROUND AND AIMS

Introgressive hybridization poses a challenge to taxonomic and phylogenetic understanding of taxa, particularly when there are high numbers of co-occurring, intercrossable species. The genus Quercus exemplifies this situation. Oaks are highly diverse in sympatry and cross freely, creating syngameons of interfertile species. Although a well-resolved, dated phylogeny is available for the American oak clade, evolutionary relationships within many of the more recently derived clades remain to be defined, particularly for the young and exceptionally diverse Mexican white oak clade. Here, we adopted an approach bridging micro- and macroevolutionary scales to resolve evolutionary relationships in a rapidly diversifying clade endemic to Mexico.

METHODS

Ecological data and sequences of 155 low-copy nuclear genes were used to identify distinct lineages within the Quercus laeta complex. Concatenated and coalescent approaches were used to assess the phylogenetic placement of these lineages relative to the Mexican white oak clade. Phylogenetic network methods were applied to evaluate the timing and genomic significance of recent or historical introgression among lineages.

KEY RESULTS

The Q. laeta complex comprises six well-supported lineages, each restricted geographically and with mostly divergent climatic niches. Species trees corroborated that the different lineages are more closely related to other species of Mexican white oaks than to each other, suggesting that this complex is polyphyletic. Phylogenetic networks estimated events of ancient introgression that involved the ancestors of three present-day Q. laeta lineages.

CONCLUSIONS

The Q. laeta complex is a morphologically and ecologically related group of species rather than a clade. Currently, oak phylogenetics is at a turning point, at which it is necessary to integrate phylogenetics and ecology in broad regional samples to figure out species boundaries. Our study illuminates one of the more complicated of the Mexican white oak groups and lays groundwork for further taxonomic study.

Homoplastic versus xenoplastic evolution: exploring the emergence of key intrinsic and extrinsic traits in the montane genus Soldanella (Primulaceae)

Specific ecological conditions in the high mountain environment exert a selective pressure that often leads to convergent trait evolution. Reticulations induced by incomplete lineage sorting and introgression can lead to discordant trait patterns among gene and species trees (hemiplasy/xenoplasy), providing a false illusion that the traits under study are homoplastic. Using phylogenetic species networks, we explored the effect of gene exchange on trait evolution in Soldanella, a genus profoundly influenced by historical introgression. At least three features evolved independently multiple times: the single-flowered dwarf phenotype, dysploid cytotype, and ecological generalism. The present analyses also indicated that the recurring occurrence of stoloniferous growth might have been prompted by an introgression event between an ancestral lineage and a still extant species, although its emergence via convergent evolution cannot be completely ruled out. Phylogenetic regression suggested that the independent evolution of larger genomes in snowbells is most likely a result of the interplay between hybridization events of dysploid and euploid taxa and hostile environments at the range margins of the genus. The emergence of key intrinsic and extrinsic traits in snowbells has been significantly impacted not only by convergent evolution but also by historical and recent introgression events.

A cyst-forming coccidian with large geographical range infecting forest and commensal rodents: Sarcocystis muricoelognathis sp. nov

BACKGROUND

The geographic distribution and host-parasite interaction networks of Sarcocystis spp. in small mammals in eastern Asia remain incompletely known.

METHODS

Experimental infections, morphological and molecular characterizations were used for discrimination of a new Sarcocystis species isolated from colubrid snakes and small mammals collected in Thailand, Borneo and China.

RESULTS

We identified a new species, Sarcocystis muricoelognathis sp. nov., that features a relatively wide geographic distribution and infects both commensal and forest-inhabiting intermediate hosts. Sarcocystis sporocysts collected from rat snakes (Coelognathus radiatus, C. flavolineatus) in Thailand induced development of sarcocysts in experimental SD rats showing a type 10a cyst wall ultrastructure that was identical with those found in Rattus norvegicus from China and the forest rat Maxomys whiteheadi in Borneo. Its cystozoites had equal sizes in all intermediate hosts and locations, while sporocysts and cystozoites were distinct from other Sarcocystis species. Partial 28S rRNA sequences of S. muricoelognathis from M. whiteheadi were largely identical to those from R. norvegicus in China but distinct from newly sequenced Sarcocystis zuoi. The phylogeny of the nuclear 18S rRNA gene placed S. muricoelognathis within the so-called S. zuoi complex, including Sarcocystis attenuati, S. kani, S. scandentiborneensis and S. zuoi, while the latter clustered with the new species. However, the phylogeny of the ITS1-region confirmed the distinction between S. muricoelognathis and S. zuoi. Moreover, all three gene trees suggested that an isolate previously addressed as S. zuoi from Thailand (KU341120) is conspecific with S. muricoelognathis. Partial mitochondrial cox1 sequences of S. muricoelognathis were almost identical with those from other members of the group suggesting a shared, recent ancestry. Additionally, we isolated two partial 28S rRNA Sarcocystis sequences from Low's squirrel Sundasciurus lowii that clustered with those of S. scandentiborneensis from treeshews.

CONCLUSIONS

Our results provide strong evidence of broad geographic distributions of rodent-associated Sarcocystis and host shifts between commensal and forest small mammal species, even if the known host associations remain likely only snapshots of the true associations.

Chloroplast Genome Provides Insights into Molecular Evolution and Species Relationship of Fleabanes (Erigeron: Tribe Astereae, Asteraceae) in the Juan Fernández Islands, Chile

Erigeron represents the third largest genus on the Juan Fernández Islands, with six endemic species, five of which occur exclusively on the younger Alejandro Selkirk Island with one species on both islands. While its continental sister species is unknown, Erigeron on the Juan Fernández Islands appears to be monophyletic and most likely evolved from South American progenitor species. We characterized the complete chloroplast genomes of five Erigeron species, including accessions of E. fernandezia and one each from Alejandro Selkirk and Robinson Crusoe Islands, with the purposes of elucidating molecular evolution and phylogenetic relationships. We found highly conserved chloroplast genomes in size, gene order and contents, and further identified several mutation hotspot regions. In addition, we found two positively selected chloroplast genes (ccsA and ndhF) among species in the islands. The complete plastome sequences confirmed the monophyly of Erigeron in the islands and corroborated previous phylogenetic relationships among species. New findings in the current study include (1) two major lineages, E. turricola-E. luteoviridis and E. fernandezia-E. ingae-E. rupicola, (2) the non-monophyly of E. fernandezia occurring on the two islands, and (3) the non-monophyly of the alpine species E. ingae complex.

A new feather mite species of the genus Mycterialges Gaud & Atyeo, 1981 (Acari, Xolalgidae) from the Oriental Stork, Ciconiaboyciana (Ciconiiformes, Ciconiidae) in Korea

A new feather mite species, Mycterialgesboycianaesp. nov. (Xolalgidae), was identified from the Oriental Stork, Ciconiaboyciana Swinhoe, 1873, in Korea. Males of M.boycianaesp. nov. are distinguished from Mycterialgesmesomorphus Gaud & Atyeo, 1981, in having a single triangular prodorsal shield, sinuous margins of the opisthosoma located between setae e2 and h2 on the hysteronotal shield, an oval-shaped epiandrum without posterior extensions, a shorter tibia + tarsus IV than femoragenu IV, and an absent ambulacral disc of leg IV. Females differ in having a prodorsal shield with a posterior margin that is blunt-angular, and a concave posterior margin of the hysteronotal shield with posterior extensions. This study presents the first record of the feather mite genus Mycterialges in birds of the genus Ciconia (Ciconiidae). Additionally, we determined the phylogenetic relationship among Ingrassiinae using the mitochondrial cytochrome c oxidase subunit (COI).

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

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

Morphological and molecular phylogenetic characterization of Sarcocystis kani sp. nov. and other novel, closely related Sarcocystis spp. infecting small mammals and colubrid snakes in Asia

We investigated the morphology and phylogenetic relationships of novel and previously recognized Sarcocystis spp. infecting small mammals and colubrid snakes in Asia. The nuclear 18S rRNA and mitochondrial cox1 of Sarcocystis sp.1 from mangrove snakes (Boiga dendrophila) in Thailand and Sarcocystis sp.2 from a ricefield rat (Rattus argentiventer) in Sumatra were partially sequenced. Sporocysts of Sarcocystis sp.1 induced development of sarcocysts in experimentally infected rats, which showed a unique ultrastructure that was observed previously by S.P. Kan in rats from Malaysia; therefore, we describe this species as Sarcocystis kani sp. nov. Its integration into the 18S rRNA phylogeny of Sarcocystis spp. cycling between small mammals and colubrid snakes helped clarify relationships among the so-called S. zuoi-complex of molecularly cryptic species: Sarcocystis kani sp. nov., S. sp.2, S. attenuati, S. scandentiborneensis, and S. zuoi were all included in this clade. Tree topology was resolved into dichotomies congruent with the morphological disparities between the taxa. However, cox1 gene sequencing (including newly sequenced S. singaporensis and S. zamani) revealed that Sarcocystis kani, S. attenuati, and S. scandentiborneensis were identical suggesting a recent, common ancestry. To identify other distinctive features, lineage-specific molecular patterns within both genes were examined revealing that all 18S rRNA sequences of the S. zuoi - complex possess a unique, 7-nt long motif in helix 38 of domain V7 that was different in S. clethrionomyelaphis which branched off basally from the complex. Three-dimensional homology modelling of COX1 protein structure identified amino acid substitutions within the barcode area specific for the S. zuoi-complex and substantial divergence in structurally important amino acids between Sarcocystis species of snakes as definitive hosts and other lineages of the Sarcocystidae. We discuss the utility of selected genes for species delimitation of the Sarcocystis spp. under investigation, which probably evolved during recent radiations of their intermediate and definitive hosts.

Incongruence in the phylogenomics era

Genome-scale data and the development of novel statistical phylogenetic approaches have greatly aided the reconstruction of a broad sketch of the tree of life and resolved many of its branches. However, incongruence - the inference of conflicting evolutionary histories - remains pervasive in phylogenomic data, hampering our ability to reconstruct and interpret the tree of life. Biological factors, such as incomplete lineage sorting, horizontal gene transfer, hybridization, introgression, recombination and convergent molecular evolution, can lead to gene phylogenies that differ from the species tree. In addition, analytical factors, including stochastic, systematic and treatment errors, can drive incongruence. Here, we review these factors, discuss methodological advances to identify and handle incongruence, and highlight avenues for future research.

Sampling affects population genetic inference: A case study of the Allen's (Selasphorus sasin) and rufous hummingbird (Selasphorus rufus)

Gene flow can affect evolutionary inference when species are undersampled. Here, we evaluate the effects of gene flow and geographic sampling on demographic inference of 2 hummingbirds that hybridize, Allen's hummingbird (Selasphorus sasin) and rufous hummingbird (Selasphorus rufus). Using whole-genome data and extensive geographic sampling, we find widespread connectivity, with introgression far beyond the Allen's × rufous hybrid zone, although the Z chromosome resists introgression beyond the hybrid zone. We test alternative hypotheses of speciation history of Allen's, rufous, and Calliope (S. calliope) hummingbird and find that rufous hummingbird is the sister taxon to Allen's hummingbird, and Calliope hummingbird is the outgroup. A model treating the 2 subspecies of Allen's hummingbird as a single panmictic population fit observed genetic data better than models treating the subspecies as distinct populations, in contrast to morphological and behavioral differences and analyses of spatial population structure. With additional sampling, our study builds upon recent studies that came to conflicting conclusions regarding the evolutionary histories of these 2 species. Our results stress the importance of thorough geographic sampling when assessing demographic history in the presence of gene flow.

Phylogenomics Reveals High Levels of Incomplete Lineage Sorting at the Ancestral Nodes of the Macaque Radiation

The genus Macaca includes 23 species assigned into 4 to 7 groups. It exhibits the largest geographic range and represents the most successful example of adaptive radiation of nonhuman primates. However, intrageneric phylogenetic relationships among species remain controversial and have not been resolved so far. In this study, we conducted a phylogenomic analysis on 16 newly generated and 8 published macaque genomes. We found strong evidence supporting the division of this genus into 7 species groups. Incomplete lineage sorting (ILS) was the primary factor contributing to the discordance observed among gene trees; however, we also found evidence of hybridization events, specifically between the ancestral arctoides/sinica and silenus/nigra lineages that resulted in the hybrid formation of the fascicularis/mulatta group. Combined with fossil data, our phylogenomic data were used to establish a scenario for macaque radiation. These findings provide insights into ILS and potential ancient introgression events that were involved in the radiation of macaques, which will lead to a better understanding of the rapid speciation occurring in nonhuman primates.

ChimeraTE: a pipeline to detect chimeric transcripts derived from genes and transposable elements

Transposable elements (TEs) produce structural variants and are considered an important source of genetic diversity. Notably, TE-gene fusion transcripts, i.e. chimeric transcripts, have been associated with adaptation in several species. However, the identification of these chimeras remains hindered due to the lack of detection tools at a transcriptome-wide scale, and to the reliance on a reference genome, even though different individuals/cells/strains have different TE insertions. Therefore, we developed ChimeraTE, a pipeline that uses paired-end RNA-seq reads to identify chimeric transcripts through two different modes. Mode 1 is the reference-guided approach that employs canonical genome alignment, and Mode 2 identifies chimeras derived from fixed or insertionally polymorphic TEs without any reference genome. We have validated both modes using RNA-seq data from four Drosophila melanogaster wild-type strains. We found ∼1.12% of all genes generating chimeric transcripts, most of them from TE-exonized sequences. Approximately ∼23% of all detected chimeras were absent from the reference genome, indicating that TEs belonging to chimeric transcripts may be recent, polymorphic insertions. ChimeraTE is the first pipeline able to automatically uncover chimeric transcripts without a reference genome, consisting of two running Modes that can be used as a tool to investigate the contribution of TEs to transcriptome plasticity.

Phylogenomics reveals extreme gene tree discordance in a lineage of dominant trees: hybridization, introgression, and incomplete lineage sorting blur deep evolutionary relationships despite clear species groupings in Eucalyptus subgenus Eudesmia

Eucalypts are a large and ecologically important group of plants on the Australian continent, and understanding their evolution is important in understanding evolution of the unique Australian flora. Previous phylogenies using plastome DNA, nuclear-ribosomal DNA, or random genome-wide SNPs, have been confounded by limited genetic sampling or by idiosyncratic biological features of the eucalypts, including widespread plastome introgression. Here we present phylogenetic analyses of Eucalyptus subgenus Eudesmia (22 species from western, northern, central and eastern Australia), in the first study to apply a target-capture sequencing approach using custom, eucalypt-specific baits (of 568 genes) to a lineage of Eucalyptus. Multiple accessions of all species were included, and target-capture data were supplemented by separate analyses of plastome genes (average of 63 genes per sample). Analyses revealed a complex evolutionary history likely shaped by incomplete lineage sorting and hybridization. Gene tree discordance generally increased with phylogenetic depth. Species, or groups of species, toward the tips of the tree are mostly supported, and three major clades are identified, but the branching order of these clades cannot be confirmed with confidence. Multiple approaches to filtering the nuclear dataset, by removing genes or samples, could not reduce gene tree conflict or resolve these relationships. Despite inherent complexities in eucalypt evolution, the custom bait kit devised for this research will be a powerful tool for investigating the evolutionary history of eucalypts more broadly.

Unveiling Cryptic Species Diversity and Genetic Variation of Lasiodiplodia (Botryosphaeriaceae, Botryosphaeriales) Infecting Fruit Crops in Taiwan

The genus Lasiodiplodia, a member of the family Botryosphaeriaceae, is an important fungal disease genus in agriculture. However, the Lasiodiplodia species survey and genetic diversity in Taiwan remain unclear. This study aimed to investigate the Lasiodiplodia species associated with various fruit species to explore the cryptic Lasiodiplodia species diversity, validate species delimitation, and unveil cryptic genetic diversity. Overall, six Lasiodiplodia species were identified, with several new records of infection identified. Additionally, phylogenetic analyses indicated that the relations of all isolates of L. theobromae might be paraphyletic. They were grouped with L. brasiliense based on Automatic Barcode Gap Discovery (ABGD), Automatic Partitioning (ASAP) and structure-based clustering analyses. These analyses did not provide conclusive evidence for L. brasiliensis as a stable species. It may be necessary to gather more information to clarify the species delineation. The multiple new records of Lasiodiplodia species with high genetic diversity and differentiation revealed that the diversity of Lasiodiplodia in Taiwan was underestimated in the past. We found that L. theobromae has the highest number of haplotypes but the lowest number of haplotype and nucleotide diversities, indicating a recent population expansion. This was supported by the significant negative Tajima's D and Fu and Li's D* tests. The high genetic diversity, low gene flow, and host-associated differentiation of Lasiodiplodia species indicate that they might harbour powerful evolutionary potential in Taiwan. This study provided critical insights into genetic variation, host-associated differentiation, and demography of Lasiodiplodia species, which would be helpful for disease management of related pathogens.

Mitonuclear discordance results from incomplete lineage sorting, with no detectable evidence for gene flow, in a rapid radiation of Todiramphus kingfishers

Many organisms possess multiple discrete genomes (i.e. nuclear and organellar), which are inherited separately and may have unique and even conflicting evolutionary histories. Phylogenetic reconstructions from these discrete genomes can yield different patterns of relatedness, a phenomenon known as cytonuclear discordance. In many animals, mitonuclear discordance (i.e. discordant evolutionary histories between the nuclear and mitochondrial genomes) has been widely documented, but its causes are often considered idiosyncratic and inscrutable. We show that a case of mitonuclear discordance in Todiramphus kingfishers can be explained by extensive genome-wide incomplete lineage sorting (ILS), likely a result of the explosive diversification history of this genus. For these kingfishers, quartet frequencies reveal that the nuclear genome is dominated by discordant topologies, with none of the internal branches in our consensus nuclear tree recovered in >50% of genome-wide gene trees. Meanwhile, a lack of inter-species shared ancestry, non-significant pairwise tests for gene flow, and little evidence for meaningful migration edges between species, leads to the conclusion that gene flow cannot explain the mitonuclear discordance we observe. This lack of evidence for gene flow combined with evidence for extensive genome-wide gene tree discordance, a hallmark of ILS, leads us to conclude that the mitonuclear discordance we observe likely results from ILS, specifically deep coalescence of the mitochondrial genome. Based on this case study, we hypothesize that similar demographic histories in other 'great speciator' taxa across the Indo-Pacific likely predispose these groups to high levels of ILS and high likelihoods of mitonuclear discordance.

Phylogenetic Analysis of Allotetraploid Species Using Polarized Genomic Sequences

Phylogenetic analysis of polyploid hybrid species has long posed a formidable challenge as it requires the ability to distinguish between alleles of different ancestral origins in order to disentangle their individual evolutionary history. This problem has been previously addressed by conceiving phylogenies as reticulate networks, using a two-step phasing strategy that first identifies and segregates homoeologous loci and then, during a second phasing step, assigns each gene copy to one of the subgenomes of an allopolyploid species. Here, we propose an alternative approach, one that preserves the core idea behind phasing-to produce separate nucleotide sequences that capture the reticulate evolutionary history of a polyploid-while vastly simplifying its implementation by reducing a complex multistage procedure to a single phasing step. While most current methods used for phylogenetic reconstruction of polyploid species require sequencing reads to be pre-phased using experimental or computational methods-usually an expensive, complex, and/or time-consuming endeavor-phasing executed using our algorithm is performed directly on the multiple-sequence alignment (MSA), a key change that allows for the simultaneous segregation and sorting of gene copies. We introduce the concept of genomic polarization that, when applied to an allopolyploid species, produces nucleotide sequences that capture the fraction of a polyploid genome that deviates from that of a reference sequence, usually one of the other species present in the MSA. We show that if the reference sequence is one of the parental species, the polarized polyploid sequence has a close resemblance (high pairwise sequence identity) to the second parental species. This knowledge is harnessed to build a new heuristic algorithm where, by replacing the allopolyploid genomic sequence in the MSA by its polarized version, it is possible to identify the phylogenetic position of the polyploid's ancestral parents in an iterative process. The proposed methodology can be used with long-read and short-read high-throughput sequencing data and requires only one representative individual for each species to be included in the phylogenetic analysis. In its current form, it can be used in the analysis of phylogenies containing tetraploid and diploid species. We test the newly developed method extensively using simulated data in order to evaluate its accuracy. We show empirically that the use of polarized genomic sequences allows for the correct identification of both parental species of an allotetraploid with up to 97% certainty in phylogenies with moderate levels of incomplete lineage sorting (ILS) and 87% in phylogenies containing high levels of ILS. We then apply the polarization protocol to reconstruct the reticulate histories of Arabidopsis kamchatica and Arabidopsis suecica, two allopolyploids whose ancestry has been well documented. [Allopolyploidy; Arabidopsis; genomic polarization; homoeologs; incomplete lineage sorting; phasing; polyploid phylogenetics; reticulate evolution.].

Identifying and addressing methodological incongruence in phylogenomics: A review

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

Population genomics indicate three different modes of divergence and speciation with gene flow in the green-winged teal duck complex

The processes leading to divergence and speciation can differ broadly among taxa with different life histories. We examine these processes in a small clade of ducks with historically uncertain relationships and species limits. The green-winged teal (Anas crecca) complex is a Holarctic species of dabbling duck currently categorized as three subspecies (Anas crecca crecca, A. c. nimia, and A. c. carolinensis) with a close relative, the yellow-billed teal (Anas flavirostris) from South America. A. c. crecca and A. c. carolinensis are seasonal migrants, while the other taxa are sedentary. We examined divergence and speciation patterns in this group, determining their phylogenetic relationships and the presence and levels of gene flow among lineages using both mitochondrial and genome-wide nuclear DNA obtained from 1,393 ultraconserved element (UCE) loci. Phylogenetic relationships using nuclear DNA among these taxa showed A. c. crecca, A. c. nimia, and A. c. carolinensis clustering together to form one polytomous clade, with A. flavirostris sister to this clade. This relationship can be summarized as (crecca, nimia, carolinensis)(flavirostris). However, whole mitogenomes revealed a different phylogeny: (crecca, nimia)(carolinensis, flavirostris). The best demographic model for key pairwise comparisons supported divergence with gene flow as the probable speciation mechanism in all three contrasts (crecca-nimia, crecca-carolinensis, and carolinensis-flavirostris). Given prior work, gene flow was expected among the Holarctic taxa, but gene flow between North American carolinensis and South American flavirostris (M ∼0.1-0.4 individuals/generation), albeit low, was not expected. Three geographically oriented modes of divergence are likely involved in the diversification of this complex: heteropatric (crecca-nimia), parapatric (crecca-carolinensis), and (mostly) allopatric (carolinensis-flavirostris). Our study shows that ultraconserved elements are a powerful tool for simultaneously studying systematics and population genomics in systems with historically uncertain relationships and species limits.

Taxonomic insights and evolutionary history in East Asian terrestrial slugs of the genus Meghimatium

East Asia, specifically the Japanese Archipelago, is a biodiversity hotspot of both vertebrates and invertebrates. Mollusks represent a burst of species diversity in this region due to the effects of biotic and abiotic factors on their morphological traits, such as shell shape and size. However, the evolutionary history of terrestrial slugs in East Asia remains unknown. In the present study, we investigated the molecular phylogeny of terrestrial slugs of the genus Meghimatium. This genus includes three described and eight undescribed species, and our study used all except for two. Based on phylogeny and the species delimitation tests, the genus Meghimatium was split into many putative species, suggesting higher species diversity than previously thought based on morphological and anatomical studies and that almost undescribed species may be inappropriate. Therefore, morphological traits, such as body size and colour, conventionally considered for classification may easily vary or be similar across geographic region. Moreover, the divergence time of this genus is almost concordant with the geographical time scale of the formation of the Japanese mainland. Our findings suggest that molecular phylogenetics helps classify Japanese Meghimatium slugs, but comprehensive taxonomic revisions using multi-locus analyses are needed.

Phylogenomic analyses of Camellia support reticulate evolution among major clades

Camellia (Theaceae) is a morphologically highly diverse genus of flowering plants and includes many famous species with high economic value, and the phylogeny of this genus is not fully resolved. We used 95 transcriptomes from 87 Camellia species and identified 1481 low-copy genes to conduct a detailed analysis of the phylogeny of this genus according to various data-screening criteria. The results show that, very different from the two existing classification systems of Camellia, 87 species are grouped into 8 main clades and two independent species, and that all 8 clades except Clade 8 were strongly supported by almost all the coalescent or concatenated trees using different gene subsets. However, the relationships among these clades were weakly supported and different from analyses using different gene subsets; furthermore, they do not agree with the phylogeny from chloroplast genomes of Camellia. Additional analyses support reticulate evolution (probably resulting from introgression or hybridization) among some major Camellia lineages, providing explanation for extensive gene tree conflicts. Furthermore, we inferred that together with the formation of East Asian subtropical evergreen broad-leaved forests, Camellia underwent a radiative divergence of major clades at 23 ∼ 19 Ma in the late Miocene then had a subsequent species burst at 10 ∼ 5 Ma. Principal component and cluster analyses provides new insights into morphological changes underlying the evolution of Camellia and a reference to further clarify subgenus and sections of this genus. The comprehensive study here including a nuclear phylogeny and other analyses reveal the rapid evolutionary history of Camellia.

Insights into cryptic speciation of quillworts in China

Cryptic species are commonly misidentified because of high morphological similarities to other species. One group of plants that may harbor large numbers of cryptic species is the quillworts (Isoëtes spp.), an ancient aquatic plant lineage. Although over 350 species of Isoëtes have been reported globally, only ten species have been recorded in China. The aim of this study is to better understand Isoëtes species diversity in China. For this purpose, we systematically explored the phylogeny and evolution of Isoëtes using complete chloroplast genome (plastome) data, spore morphology, chromosome number, genetic structure, and haplotypes of almost all Chinese Isoëtes populations. We identified three ploidy levels of Isoëtes in China-diploid (2n = 22), tetraploid (2n = 44), and hexaploid (2n = 66). We also found four megaspore and microspore ornamentation types in diploids, six in tetraploids, and three in hexaploids. Phylogenetic analyses confirmed that I. hypsophila as the ancestral group of the genus and revealed that Isoëtes diploids, tetraploids, and hexaploids do not form monophyletic clades. Most individual species possess a single genetic structure; however, several samples have conflicting positions on the phylogenetic tree based on SNPs and the tree based on plastome data. All 36 samples shared 22 haplotypes. Divergence time analysis showed that I. hypsophila diverged in the early Eocene (∼48.05 Ma), and most other Isoëtes species diverged 3-20 Ma. Additionally, different species of Isoëtes were found to inhabit different water systems and environments along the Yangtze River. These findings provide new insights into the relationships among Isoëtes species in China, where highly similar morphologic populations may harbor many cryptic species.

A New Assessment of Robust Capuchin Monkey (Sapajus) Evolutionary History Using Genome-Wide SNP Marker Data and a Bayesian Approach to Species Delimitation

Robust capuchin monkeys, Sapajus genus, are among the most phenotypically diverse and widespread groups of primates in South America, with one of the most confusing and often shifting taxonomies. We used a ddRADseq approach to generate genome-wide SNP markers for 171 individuals from all putative extant species of Sapajus to access their evolutionary history. Using maximum likelihood, multispecies coalescent phylogenetic inference, and a Bayes Factor method to test for alternative hypotheses of species delimitation, we inferred the phylogenetic history of the Sapajus radiation, evaluating the number of discrete species supported. Our results support the recognition of three species from the Atlantic Forest south of the São Francisco River, with these species being the first splits in the robust capuchin radiation. Our results were congruent in recovering the Pantanal and Amazonian Sapajus as structured into three monophyletic clades, though new morphological assessments are necessary, as the Amazonian clades do not agree with previous morphology-based taxonomic distributions. Phylogenetic reconstructions for Sapajus occurring in the Cerrado, Caatinga, and northeastern Atlantic Forest were less congruent with morphology-based phylogenetic reconstructions, as the bearded capuchin was recovered as a paraphyletic clade, with samples from the Caatinga biome being either a monophyletic clade or nested with the blond capuchin monkey.

The phylogenomics and evolutionary dynamics of the organellar genomes in carnivorous Utricularia and Genlisea species (Lentibulariaceae)

Utricularia and Genlisea are highly specialized carnivorous plants whose phylogenetic history has been poorly explored using phylogenomic methods. Additional sampling and genomic data are needed to advance our phylogenetic and taxonomic knowledge of this group of plants. Within a comparative framework, we present a characterization of plastome (PT) and mitochondrial (MT) genes of 26 Utricularia and six Genlisea species, with representatives of all subgenera and growth habits. All PT genomes maintain similar gene content, showing minor variation across the genes located between the PT junctions. One exception is a major variation related to different patterns in the presence and absence of ndh genes in the small single copy region, which appears to follow the phylogenetic history of the species rather than their lifestyle. All MT genomes exhibit similar gene content, with most differences related to a lineage-specific pseudogenes. We find evidence for episodic positive diversifying selection in PT and for most of the Utricularia MT genes that may be related to the current hypothesis that bladderworts' nuclear DNA is under constant ROS oxidative DNA damage and unusual DNA repair mechanisms, or even low fidelity polymerase that bypass lesions which could also be affecting the organellar genomes. Finally, both PT and MT phylogenetic trees were well resolved and highly supported, providing a congruent phylogenomic hypothesis for Utricularia and Genlisea clade given the study sampling.

Geographic Genetic Structure of Alectoris chukar in Türkiye: Post-LGM-Induced Hybridization and Human-Mediated Contaminations

Türkiye is considered an important evolutionary area for Chukar partridge (Alectoris chukar), since it is both a potential ancestral area and a diversification center for the species. Using 2 mitochondrial (Cty-b and D-loop) and 13 polymorphic microsatellite markers, we investigated the geographic genetic structure of A. chukar populations to determine how past climatic fluctuations and human activities have shaped the gene pool of this species in Türkiye. Our results indicate, firstly, that only A. chukar of the genus Alectoris is present in Türkiye (Anatolia and Thrace), with no natural or artificial gene flow from congenerics. Secondly, the geographic genetic structure of the species in Türkiye has been shaped by topographic heterogeneity, Pleistocene climatic fluctuations, and artificial transport by humans. Third, there appears to be three genetic clusters: Thracian, Eastern, and Western. Fourth, the post-LGM demographic expansion of the Eastern and Western populations has formed a hybrid zone in Central Anatolia (~8 kyBP). Fifth, the rate of China clade-B contamination in Türkiye is about 8% in mtDNA and about 12% in nuDNA, with the Southeastern Anatolian population having the highest contamination. Sixth, the Thracian population was the most genetically distinct, with the lowest genetic diversity and highest level of inbreeding and no China clad-B contamination. These results can contribute to the conservation regarding A. chukar populations, especially the Thracian population.

Testing Multispecies Coalescent Simulators Using Summary Statistics

As genomic-scale datasets motivate research on species tree inference, simulators of the multispecies coalescent (MSC) process have become essential for the testing and evaluation of new inference methods. However, the simulators themselves must be tested to ensure that they give valid samples. This work develops methods for checking whether a collection of gene trees is in accord with the MSC model on a given species tree. When applied to well-known simulators, we find that several give flawed samples. The tests presented are capable of validating both topological and metric properties of gene tree samples, and are implemented in a freely available R package MSCsimtester so that developers and users may easily apply them.

Not out of the Mediterranean: Atlantic populations of the gorgonian Paramuricea clavata are a separate sister species under further lineage diversification

The accurate delimitation of species boundaries in nonbilaterian marine taxa is notoriously difficult, with consequences for many studies in ecology and evolution. Anthozoans are a diverse group of key structural organisms worldwide, but the lack of reliable morphological characters and informative genetic markers hampers our ability to understand species diversification. We investigated population differentiation and species limits in Atlantic (Iberian Peninsula) and Mediterranean lineages of the octocoral genus Paramuricea previously identified as P. clavata. We used a diverse set of molecular markers (microsatellites, RNA-seq derived single-copy orthologues [SCO] and mt-mutS [mitochondrial barcode]) at 49 locations. Clear segregation of Atlantic and Mediterranean lineages was found with all markers. Species-tree estimations based on SCO strongly supported these two clades as distinct, recently diverged sister species with incomplete lineage sorting, P. cf. grayi and P. clavata, respectively. Furthermore, a second putative (or ongoing) speciation event was detected in the Atlantic between two P. cf. grayi color morphotypes (yellow and purple) using SCO and supported by microsatellites. While segregating P. cf. grayi lineages showed considerable geographic structure, dominating circalittoral communities in southern (yellow) and western (purple) Portugal, their occurrence in sympatry at some localities suggests a degree of reproductive isolation. Overall, our results show that previous molecular and morphological studies have underestimated species diversity in Paramuricea occurring in the Iberian Peninsula, which has important implications for conservation planning. Finally, our findings validate the usefulness of phylotranscriptomics for resolving evolutionary relationships in octocorals.

A turn in species conservation for hairpin banksias: demonstration of oversplitting leads to better management of diversity

PREMISE

Understanding evolutionary history and classifying discrete units of organisms remain overwhelming tasks, and lags in this workload concomitantly impede an accurate documentation of biodiversity and conservation management. Rapid advances and improved accessibility of sensitive high-throughput sequencing tools are fortunately quickening the resolution of morphological complexes and   thereby improving the estimation of species diversity. The recently described and critically endangered Banksia vincentia is morphologically similar to the hairpin banksia complex (B. spinulosa s.l.), a group of eastern Australian flowering shrubs whose continuum of morphological diversity has been responsible for taxonomic controversy and possibly questionable conservation initiatives.

METHODS

To assist conservation while testing the current taxonomy of this group, we used high-throughput sequencing to infer a population-scale evolutionary scenario for a sample set that is comprehensive in its representation of morphological diversity and a 2500-km distribution.

RESULTS

Banksia spinulosa s.l. represents two clades, each with an internal genetic structure shaped through historical separation by biogeographic barriers. This structure conflicts with the existing taxonomy for the group. Corroboration between phylogeny and population statistics aligns with the hypothesis that B. collina, B. neoanglica, and B. vincentia should not be classified as species.

CONCLUSIONS

The pattern here supports how morphological diversity can be indicative of a locally expressed suite of traits rather than relationship. Oversplitting in the hairpin banksias is atypical since genomic analyses often reveal that species diversity is underestimated. However, we show that erring on overestimation can yield negative consequences, such as the disproportionate prioritization of a geographically anomalous population.

Phylogenomic Analyses of 2,786 Genes in 158 Lineages Support a Root of the Eukaryotic Tree of Life between Opisthokonts and All Other Lineages

Advances in phylogenomics and high-throughput sequencing have allowed the reconstruction of deep phylogenetic relationships in the evolution of eukaryotes. Yet, the root of the eukaryotic tree of life remains elusive. The most popular hypothesis in textbooks and reviews is a root between Unikonta (Opisthokonta + Amoebozoa) and Bikonta (all other eukaryotes), which emerged from analyses of a single-gene fusion. Subsequent, highly cited studies based on concatenation of genes supported this hypothesis with some variations or proposed a root within Excavata. However, concatenation of genes does not consider phylogenetically-informative events like gene duplications and losses. A recent study using gene tree parsimony (GTP) suggested the root lies between Opisthokonta and all other eukaryotes, but only including 59 taxa and 20 genes. Here we use GTP with a duplication-loss model in a gene-rich and taxon-rich dataset (i.e., 2,786 gene families from two sets of 155 and 158 diverse eukaryotic lineages) to assess the root, and we iterate each analysis 100 times to quantify tree space uncertainty. We also contrasted our results and discarded alternative hypotheses from the literature using GTP and the likelihood-based method SpeciesRax. Our estimates suggest a root between Fungi or Opisthokonta and all other eukaryotes; but based on further analysis of genome size, we propose that the root between Opisthokonta and all other eukaryotes is the most likely.

Phylogenomics and Diversification of the Schistosomatidae Based on Targeted Sequence Capture of Ultra-Conserved Elements

Schistosomatidae Stiles and Hassall 1898 is a medically significant family of digenetic trematodes (Trematoda: Digenea), members of which infect mammals or birds as definitive hosts and aquatic or amphibious gastropods as intermediate hosts. Currently, there are 17 named genera, for many of which evolutionary interrelationships remain unresolved. The lack of a resolved phylogeny has encumbered our understanding of schistosomatid evolution, specifically patterns of host-use and the role of host-switching in diversification. Here, we used targeted sequence capture of ultra-conserved elements (UCEs) from representatives of 13 of the 17 named genera and 11 undescribed lineages that are presumed to represent either novel genera or species to generate a phylogenomic dataset for the estimation of schistosomatid interrelationships. This study represents the largest phylogenetic effort within the Schistosomatidae in both the number of loci and breadth of taxon sampling. We present a near-comprehensive family-level phylogeny providing resolution to several clades of long-standing uncertainty within Schistosomatidae, including resolution for the placement of the North American mammalian schistosomes, implying a second separate capture of mammalian hosts. Additionally, we present evidence for the placement of Macrobilharzia at the base of the Schistosoma + Bivitellobilharzia radiation. Patterns of definitive and intermediate host use and a strong role for intermediate host-switching are discussed relative to schistosomatid diversification.

Gene expression phylogenies and ancestral transcriptome reconstruction resolves major transitions in the origins of pregnancy

Structural and physiological changes in the female reproductive system underlie the origins of pregnancy in multiple vertebrate lineages. In mammals, the glandular portion of the lower reproductive tract has transformed into a structure specialized for supporting fetal development. These specializations range from relatively simple maternal nutrient provisioning in egg-laying monotremes to an elaborate suite of traits that support intimate maternal-fetal interactions in Eutherians. Among these traits are the maternal decidua and fetal component of the placenta, but there is considerable uncertainty about how these structures evolved. Previously, we showed that changes in uterine gene expression contributes to several evolutionary innovations during the origins of pregnancy (Mika et al., 2021b). Here, we reconstruct the evolution of entire transcriptomes ('ancestral transcriptome reconstruction') and show that maternal gene expression profiles are correlated with degree of placental invasion. These results indicate that an epitheliochorial-like placenta evolved early in the mammalian stem-lineage and that the ancestor of Eutherians had a hemochorial placenta, and suggest maternal control of placental invasiveness. These data resolve major transitions in the evolution of pregnancy and indicate that ancestral transcriptome reconstruction can be used to study the function of ancestral cell, tissue, and organ systems.

Admixture of evolutionary rates across a butterfly hybrid zone

Hybridization is a major evolutionary force that can erode genetic differentiation between species, whereas reproductive isolation maintains such differentiation. In studying a hybrid zone between the swallowtail butterflies Papilio syfanius and Papilio maackii (Lepidoptera: Papilionidae), we made the unexpected discovery that genomic substitution rates are unequal between the parental species. This phenomenon creates a novel process in hybridization, where genomic regions most affected by gene flow evolve at similar rates between species, while genomic regions with strong reproductive isolation evolve at species-specific rates. Thus, hybridization mixes evolutionary rates in a way similar to its effect on genetic ancestry. Using coalescent theory, we show that the rate-mixing process provides distinct information about levels of gene flow across different parts of genomes, and the degree of rate-mixing can be predicted quantitatively from relative sequence divergence ([Formula: see text]) between the hybridizing species at equilibrium. Overall, we demonstrate that reproductive isolation maintains not only genomic differentiation, but also the rate at which differentiation accumulates. Thus, asymmetric rates of evolution provide an additional signature of loci involved in reproductive isolation.

The bounded coalescent model: Conditioning a genealogy on a minimum root date

The coalescent model represents how individuals sampled from a population may have originated from a last common ancestor. The bounded coalescent model is obtained by conditioning the coalescent model such that the last common ancestor must have existed after a certain date. This conditioned model arises in a variety of applications, such as speciation, horizontal gene transfer or transmission analysis, and yet the bounded coalescent model has not been previously analysed in detail. Here we describe a new algorithm to simulate from this model directly, without resorting to rejection sampling. We show that this direct simulation algorithm is more computationally efficient than the rejection sampling approach. We also show how to calculate the probability of the last common ancestor occurring after a given date, which is required to compute the probability density of realisations under the bounded coalescent model. Our results are applicable in both the isochronous (when all samples have the same date) and heterochronous (where samples can have different dates) settings. We explore the effect of setting a bound on the date of the last common ancestor, and show that it affects a number of properties of the resulting phylogenies. All our methods are implemented in a new R package called BoundedCoalescent which is freely available online.

Publicly Available and Validated DNA Reference Sequences Are Critical to Fungal Identification and Global Plant Protection Efforts: A Use-Case in Colletotrichum

Publicly available and validated DNA reference sequences useful for phylogeny estimation and identification of fungal pathogens are an increasingly important resource in the efforts of plant protection organizations to facilitate safe international trade of agricultural commodities. Colletotrichum species are among the most frequently encountered and regulated plant pathogens at U.S. ports-of-entry. The RefSeq Targeted Loci (RTL) project at NCBI (BioProject no. PRJNA177353) contains a database of curated fungal internal transcribed spacer (ITS) sequences that interact extensively with NCBI Taxonomy, resulting in verified name-strain-sequence type associations for >12,000 species. We present a publicly available dataset of verified and curated name-type strain-sequence associations for all available Colletotrichum species. This includes an updated GenBank Taxonomy for 238 species associated with up to 11 protein coding loci and an updated RTL ITS dataset for 226 species. We demonstrate that several marker loci are well suited for phylogenetic inference and identification. We improve understanding of phylogenetic relationships among verified species, verify or improve phylogenetic circumscriptions of 14 species complexes, and reveal that determining relationships among these major clades will require additional data. We present detailed comparisons between phylogenetic and similarity-based approaches to species identification, revealing complex patterns among single marker loci that often lead to misidentification when based on single-locus similarity approaches. We also demonstrate that species-level identification is elusive for a subset of samples regardless of analytical approach, which may be explained by novel species diversity in our dataset and incomplete lineage sorting and lack of accumulated synapomorphies at these loci.

Phylogenetics in Space: How Continuous Spatial Structure Impacts Tree Inference

The tendency to discretize biology permeates taxonomy and systematics, leading to models that simplify the often continuous nature of populations. Even when the assumption of panmixia is relaxed, most models still assume some degree of discrete structure. The multispecies coalescent has emerged as a powerful model in phylogenetics, but in its common implementation is entirely space-independent - what we call the "missing z-axis". In this article, we review the many lines of evidence for how continuous spatial structure can impact phylogenetic inference. We illustrate and expand on these by using complex continuous-space demographic models that include distinct modes of speciation. We find that the impact of spatial structure permeates all aspects of phylogenetic inference, including gene tree stoichiometry, topological and branch-length variance, network estimation, and species delimitation. We conclude by utilizing our results to suggest how researchers can identify spatial structure in phylogenetic datasets.

Highly diversified mitochondrial genomes provide new evidence for interordinal relationships in the Arachnida

Arachnida is an exceptionally diverse class in the Arthropoda, consisting of 20 orders and playing crucial roles in the terrestrial ecosystems. However, their interordinal relationships have been debated for over a century. Rearranged or highly rearranged mitochondrial genomes (mitogenomes) were consistently found in this class, but their various extent in different lineages and efficiency for resolving arachnid phylogenies are unclear. Here, we reconstructed phylogenetic trees using mitogenome sequences of 290 arachnid species to decipher interordinal relationships as well as diversification through time. Our results recovered monophyly of ten orders (i.e. Amblypygi, Araneae, Ixodida, Mesostigmata, Opiliones, Pseudoscorpiones, Ricinulei, Sarcoptiformes, Scorpiones and Solifugae), while rejecting monophyly of the Trombidiformes due to the unstable position of the Eriophyoidea. The monophyly of Acari (subclass) was rejected, possibly due to the long-branch attraction of the Pseudoscorpiones. The monophyly of Arachnida was further rejected because the Xiphosura nested within arachnid orders with unstable positions. Mitogenomes that are highly rearranged in mites but less rearranged or conserved in the remaining lineages point to their exceptional diversification in mite orders; however, shared derived mitochondrial (mt) gene clusters were found within superfamilies rather than interorders, confusing phylogenetic signals in arachnid interordinal relationships. Molecular dating results show that arachnid orders have ancient origins, ranging from the Ordovician to the Carboniferous, yet have significantly diversified since the Cretaceous in orders Araneae, Mesostigmata, Sarcoptiformes, and Trombidiformes. By summarizing previously resolved key positions of some orders, we propose a plausible arachnid tree of life. Our results underline a more precise framework for interordinal phylogeny in the Arachnida and provide new insights into their ancient evolution.

Origin and diversification of the saguaro cactus (Carnegiea gigantea): a within-species phylogenomic analysis

Reconstructing accurate historical relationships within a species poses numerous challenges, not least in many plant groups in which gene flow is high enough to extend well beyond species boundaries. Nonetheless, the extent of tree-like history within a species is an empirical question on which it is now possible to bring large amounts of genome sequence to bear. We assess phylogenetic structure across the geographic range of the saguaro cactus, an emblematic member of Cactaceae, a clade known for extensive hybridization and porous species boundaries. Using 200 Gb of whole genome resequencing data from 20 individuals sampled from 10 localities, we assembled two data sets comprising 150,000 biallelic single nucleotide polymorphisms (SNPs) from protein coding sequences. From these we inferred within-species trees and evaluated their significance and robustness using five qualitatively different inference methods. Despite the low sequence diversity, large census population sizes, and presence of wide-ranging pollen and seed dispersal agents, phylogenetic trees were well resolved and highly consistent across both data sets and all methods. We inferred that the most likely root, based on marginal likelihood comparisons, is to the east and south of the region of highest genetic diversity, which lies along the coast of the Gulf of California in Sonora, Mexico. Together with striking decreases in marginal likelihood found to the north, this supports hypotheses that saguaro's current range reflects post-glacial expansion from the refugia in the south of its range. We conclude with observations about practical and theoretical issues raised by phylogenomic data sets within species, in which SNP-based methods must be used rather than gene tree methods that are widely used when sequence divergence is higher. These include computational scalability, inference of gene flow, and proper assessment of statistical support in the presence of linkage effects.

Biogeographic consequences of shifting climate for the western massasauga ( Sistrurus tergeminus )

The western massasauga (Sistrurus tergeminus) is a small pit viper with an extensive geographic range, yet observations of this species are relatively rare. They persist in patchy and isolated populations, threatened by habitat destruction and fragmentation, mortality from vehicle collisions, and deliberate extermination. Changing climates may pose an additional stressor on the survival of isolated populations. Here, we evaluate historic, modern, and future geographic projections of suitable climate for S. tergeminus to outline shifts in their potential geographic distribution and inform current and future management. We used maximum entropy modeling to build multiple models of the potential geographic distribution of S. tergeminus. We evaluated the influence of five key decisions made during the modeling process on the resulting geographic projections of the potential distribution, allowing us to identify areas of model robustness and uncertainty. We evaluated models with the area under the receiver operating curve and true skill statistic. We retained 16 models to project both in the past and future multiple general circulation models. At the last glacial maximum, the potential geographic distribution associated with S. tergeminus occurrences had a stronghold in the southern part of its current range and extended further south into Mexico, but by the mid‐Holocene, its modeled potential distribution was similar to its present‐day potential distribution. Under future model projections, the potential distribution of S. tergeminus moves north, with the strongest northward trends predicted under a climate scenario increase of 8.5 W/m². Some southern populations of S. tergeminus have likely already been extirpated and will continue to be threatened by shifting availability of suitable climate, as they are already under threat from desertification of grasslands. Land use and habitat loss at the northern edge of the species range are likely to make it challenging for this species to track suitable climates northward over time.

Molecular phylogeny of Asian Ardisia (Myrsinoideae, Primulaceae) and their leaf-nodulated endosymbionts, Burkholderia s.l. (Burkholderiaceae)

The genus Ardisia (Myrsinoideae, Primulaceae) has 16 subgenera and over 700 accepted names, mainly distributed in tropical Asia and America. The circumscription of Ardisia is not well-defined and sometimes confounded with the separation of some small genera. A taxonomic revision focusing on Ardisia and allies is necessary. In the Ardisia subgenus Crispardisia, symbiotic association with leaf-nodule bacteria is a unique character within the genus. The endosymbionts are vertically transmitted, highly specific and highly dependent on the hosts, suggesting strict cospeciation may have occurred in the evolutionary history. In the present study, we aimed to establish a phylogenetic framework for further taxonomic revision. We also aimed to test the cospeciation hypothesis of the leaf-nodulate Ardisia and their endosymbiotic bacteria. Nuclear ITS and two chloroplast intergenic spaces were used to reconstruct the phylogeny of Asian Ardisia and relatives in Myrsinoideae, Primulaceae. The 16S-23S rRNA were used to reconstruct the bacterial symbionts’ phylogeny. To understand the evolutionary association of the Ardisia and symbionts, topology tests and cophylogenetic analyses were conducted. The molecular phylogeny suggested Ardisia is not monophyletic, unless Sardiria, Hymenandra, Badula and Oncostemum are included. The results suggest the generic limit within Myrsinoideae (Primulaceae) needs to be further revised. The subgenera Crispardisia, Pimelandra, and Stylardisia were supported as monophyly, while the subgenus Bladhia was separated into two distant clades. We proposed to divide the subgenus Bladhia into subgenus Bladhia s.str. and subgenus Odontophylla. Both of the cophylogenetic analyses and topology tests rejected strict cospeciation hypothesis between Ardisia hosts and symbiotic Burkholderia. Cophylogenetic analyses showed general phylogenetic concordance of Ardisia and Burkholderia, and cospeciation events, host-switching events and loss events were all inferred.

Complex and reticulate origin of edible roses (Rosa, Rosaceae) in China

While roses are today among the most popular ornamental plants, the petals and fruits of some cultivars have flavored foods for millennia. The genetic origins of these edible cultivars remain poorly investigated. We collected the major varieties of edible roses available in China, assembled their plastome sequences, and phased the haplotypes for internal transcribed spacers (ITS1/ITS2) of the 18S-5.8S-26S nuclear ribosomal cistron. Our phylogenetic reconstruction using 88 plastid genomes, of primarily maternal origin, uncovered well-supported genetic relationships within Rosa, including all sections and all subgenera. We phased the ITS sequences to identify potential donor species ancestral to the development of known edible cultivars. The tri-parental Middle-Eastern origin of R. × damascena, the species most widely used in perfume products and food additives, was confirmed as a descendent of past hybridizations among R. moschata, R. gallica, and R. majalis/R. fedtschenkoana/R. davurica. In contrast, R. chinensis, R. rugosa, and R. gallica, in association with six other wild species, were the main donors for fifteen varieties of edible roses. The domesticated R. rugosa 'Plena' was shown to be a hybrid between R. rugosa and R. davurica, sharing a common origin with R. 'Fenghua'. Only R. 'Jinbian' and R. 'Crimson Glory' featured continuous flowering. All remaining cultivars of edible roses bloomed only once a year. Our study provides important resources for clarifying the origin of edible roses and suggests a future for breeding new cultivars with unique traits, such as continuous flowering.

Nucleotide Usage Biases Distort Inferences of the Species Tree

Despite the importance of natural selection in species’ evolutionary history, phylogenetic methods that take into account population-level processes typically ignore selection. The assumption of neutrality is often based on the idea that selection occurs at a minority of loci in the genome and is unlikely to compromise phylogenetic inferences significantly. However, genome-wide processes like GC-bias and some variation segregating at the coding regions are known to evolve in the nearly neutral range. As we are now using genome-wide data to estimate species trees, it is natural to ask whether weak but pervasive selection is likely to blur species tree inferences. We developed a polymorphism-aware phylogenetic model tailored for measuring signatures of nucleotide usage biases to test the impact of selection in the species tree. Our analyses indicate that although the inferred relationships among species are not significantly compromised, the genetic distances are systematically underestimated in a node-height-dependent manner: that is, the deeper nodes tend to be more underestimated than the shallow ones. Such biases have implications for molecular dating. We dated the evolutionary history of 30 worldwide fruit fly populations, and we found signatures of GC-bias considerably affecting the estimated divergence times (up to 23%) in the neutral model. Our findings call for the need to account for selection when quantifying divergence or dating species evolution.

OUP accepted manuscript

Most species are extinct, those that are not are often unknown. Sequenced and sampled species are often a minority of known ones. Past evolutionary events involving horizontal gene flow, such as horizontal gene transfer, hybridization, introgression, and admixture, are therefore likely to involve “ghosts,” that is extinct, unknown, or unsampled lineages. The existence of these ghost lineages is widely acknowledged, but their possible impact on the detection of gene flow and on the identification of the species involved is largely overlooked. It is generally considered as a possible source of error that, with reasonable approximation, can be ignored. We explore the possible influence of absent species on an evolutionary study by quantifying the effect of ghost lineages on introgression as detected by the popular D-statistic method. We show from simulated data that under certain frequently encountered conditions, the donors and recipients of horizontal gene flow can be wrongly identified if ghost lineages are not taken into account. In particular, having a distant outgroup, which is usually recommended, leads to an increase in the error probability and to false interpretations in most cases. We conclude that introgression from ghost lineages should be systematically considered as an alternative possible, even probable, scenario. [ABBA–BABA; D-statistic; gene flow; ghost lineage; introgression; simulation.]

A customised target capture sequencing tool for molecular identification of Aloe vera and relatives

Plant molecular identification studies have, until recently, been limited to the use of highly conserved markers from plastid and other organellar genomes, compromising resolution in highly diverse plant clades. Due to their higher evolutionary rates and reduced paralogy, low-copy nuclear genes overcome this limitation but are difficult to sequence with conventional methods and require high-quality input DNA. Aloe vera and its relatives in the Alooideae clade (Asphodelaceae, subfamily Asphodeloideae) are of economic interest for food and health products and have horticultural value. However, pressing conservation issues are increasing the need for a molecular identification tool to regulate the trade. With > 600 species and an origin of ± 15 million years ago, this predominantly African succulent plant clade is a diverse and taxonomically complex group for which low-copy nuclear genes would be desirable for accurate species discrimination. Unfortunately, with an average genome size of 16.76 pg, obtaining high coverage sequencing data for these genes would be prohibitively costly and computationally demanding. We used newly generated transcriptome data to design a customised RNA-bait panel targeting 189 low-copy nuclear genes in Alooideae. We demonstrate its efficacy in obtaining high-coverage sequence data for the target loci on Illumina sequencing platforms, including degraded DNA samples from museum specimens, with considerably improved phylogenetic resolution. This customised target capture sequencing protocol has the potential to confidently indicate phylogenetic relationships of Aloe vera and related species, as well as aid molecular identification applications.

The dynamics of introgression across an avian radiation

Hybridization and resulting introgression can play both a destructive and a creative role in the evolution of diversity. Thus, characterizing when and where introgression is most likely to occur can help us understand the causes of diversification dynamics. Here, we examine the prevalence of and variation in introgression using phylogenomic data from a large (1300+ species), geographically widespread avian group, the suboscine birds. We first examine patterns of gene tree discordance across the geographic distribution of the entire clade. We then evaluate the signal of introgression in a subset of 206 species triads using Patterson's D‐statistic and test for associations between introgression signal and evolutionary, geographic, and environmental variables. We find that gene tree discordance varies across lineages and geographic regions. The signal of introgression is highest in cases where species occur in close geographic proximity and in regions with more dynamic climates since the Pleistocene. Our results highlight the potential of phylogenomic datasets for examining broad patterns of hybridization and suggest that the degree of introgression between diverging lineages might be predictable based on the setting in which they occur.

Development and comparison of loop-mediated isothermal amplification with quantitative PCR for the specific detection of Saprolegnia spp

Saprolegniasis is an important disease in freshwater aquaculture, and is associated with oomycete pathogens in the genus Saprolegnia. Early detection of significant levels of Saprolegnia spp. pathogens would allow informed decisions for treatment which could significantly reduce losses. This study is the first to report the development of loop-mediated isothermal amplification (LAMP) for the detection of Saprolegnia spp. and compares it with quantitative PCR (qPCR). The developed protocols targeted the internal transcribed spacer (ITS) region of ribosomal DNA and the cytochrome C oxidase subunit 1 (CoxI) gene and was shown to be specific only to Saprolegnia genus. This LAMP method can detect as low as 10 fg of S. salmonis DNA while the qPCR method has a detection limit of 2 pg of S. salmonis DNA, indicating the superior sensitivity of LAMP compared to qPCR. When applied to detect the pathogen in water samples, both methods could detect the pathogen when only one zoospore of Saprolegnia was present. We propose LAMP as a quick (about 20–60 minutes) and sensitive molecular diagnostic tool for the detection of Saprolegnia spp. suitable for on-site applications.