Target-capture phylogenomics provide insights on gene and species tree discordances in Old World treefrogs (Anura: Rhacophoridae)

Phylogenomic data exploration with increased sampling provides new insights into the higher-level relationships of butterflies and moths (Lepidoptera)

A robust and stable phylogenetic framework is a fundamental goal of evolutionary biology. As the third largest insect order in the world following Coleoptera and Diptera, Lepidoptera (butterflies and moths) play a central role in almost every terrestrial ecosystem as indicators of environmental change and serve as important models for biologists exploring questions related to ecology and evolutionary biology. However, for such a charismatic insect group, the higher-level phylogenetic relationships among its superfamilies are still poorly resolved. Compared to earlier phylogenomic studies, we increased taxon sampling among Lepidoptera (37 superfamilies and 68 families containing 263 taxa) and acquired a series of large amino-acid datasets from 69,680 to 400,330 for phylogenomic reconstructions. Using these datasets, we explored the effect of different taxon sampling with significant increases in the number of included genes on tree topology by considering a series of systematic errors using maximum-likelihood (ML) and Bayesian inference (BI) methods. Moreover, we also tested the effectiveness in topology robustness among the three ML-based models. The results showed that taxon sampling is an important determinant in tree robustness of accurate lepidopteran phylogenetic estimation. Long-branch attraction (LBA) caused by site-wise heterogeneity is a significant source of bias giving rise to unstable positions of ditrysian groups in phylogenomic reconstruction. Phylogenetic inference showed the most comprehensive framework to reveal the relationships among lepidopteran superfamilies, and presented some newly relationships with strong supports (Papilionoidea was sister to Gelechioidea and Immoidea was sister to Galacticoidea, respectively), but limited by taxon sampling, the relationships within the species-rich and relatively rapid radiation Ditrysia and especially Apoditrysia remain poorly resolved, which need to increase taxon sampling for further phylogenomic reconstruction. The present study demonstrates that taxon sampling is an important determinant for an accurate lepidopteran tree of life and provides some essential insights for future lepidopteran phylogenomic studies.

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.

Integrated phylogenetic analyses reveal the evolutionary, biogeographic, and diversification history of Asian warty treefrog genus Theloderma (Anura, Rhacophoridae)

Asian warty treefrogs, genus Theloderma, are morphologically variable arboreal frogs endemic to Southeast Asia and Southern China. However, integrated systematic studies are lacking, and knowledge of the genus in terms of diversity, origin, and historical diversification remains limited. To address these knowledge gaps, we used three mitochondrial and five nuclear gene fragments to reconstruct the Theloderma phylogeny, estimate divergence times, and examine the biogeography of the genus. Phylogenetic and species delimitation analyses suggest that the genus Theloderma comprises three major clades corresponding to two subgenera and seven species groups, and mPTP identified at least 12 putative cryptic species, suggesting that species diversity has been underestimated. Biogeographic analyses indicated that most recent common ancestor of Theloderma originated in the Indochina Peninsula during the Middle Oligocene (ca. 27.77 Ma) and the splitting of Clade A to C occurred in the Late Oligocene (ca. 23.55-25.57 Ma). Current biogeographic patterns result from two distinct processes: in situ diversification in the Indochina Peninsula and dispersal in multiple areas, namely southward dispersal to the Malay Peninsula and Borneo, northeastward dispersal to Southern China, northward dispersal to the Himalayas, and dispersal from Southern China to the Indochina Peninsula. Ancestral character reconstruction suggests that the ancestor of Theloderma may have possessed a small body size, rough dorsal skin, and absence of vomerine teeth and hand webbing, and that these four characters have undergone multiple evolutions. Principal component analysis based on eight bioclimatic variables did not clearly distinguish the three major clades of Theloderma, suggesting that species in these clades may occupy similar climatic ecological niches. Our research highlights the importance of orogeny and paleoclimatic changes, in shaping amphibian biodiversity in mountain ecosystems.

Their fates intertwined: diversification patterns of the Asian gliding vertebrates may have been forged by dipterocarp trees

The repeated evolution of gliding in diverse Asian vertebrate lineages is hypothesized to have been triggered by the dominance of tall dipterocarp trees in the tropical forests of Southeast Asia. These dipterocarp forests have acted as both centres of diversification and climatic refugia for gliding vertebrates, and support most of their extant diversity. We predict similarities in the diversification patterns of dipterocarp trees and gliding vertebrates, and specifically test whether episodic diversification events such as rate shifts and/or mass extinctions were temporally congruent in these groups. We analysed diversification patterns in reconstructed timetrees of Asian dipterocarps, the most speciose gliding vertebrates from different classes (Draco lizards, gliding frogs and Pteromyini squirrels) and compared them with similar-sized clades of non-gliding relatives (Diploderma lizards, Philautus frogs and Callosciurinae squirrels) from Southeast Asia. We found significant declines in net-diversification rates of dipterocarps and the gliding vertebrates during the Pliocene-Pleistocene, but not in the non-gliding groups. We conclude that the homogeneity and temporal coincidence of these rate declines point to a viable ecological correlation between dipterocarps and the gliding vertebrates. Further, we suggest that while the diversification decay in dipterocarps was precipitated by post-Miocene aridification of Asia, the crises in the gliding vertebrates were induced by both events concomitantly.

The genomic landscape, causes, and consequences of extensive phylogenomic discordance in Old World mice and rats

A species tree is a central concept in evolutionary biology whereby a single branching phylogeny reflects relationships among species. However, the phylogenies of different genomic regions often differ from the species tree. Although tree discordance is often widespread in phylogenomic studies, we still lack a clear understanding of how variation in phylogenetic patterns is shaped by genome biology or the extent to which discordance may compromise comparative studies. We characterized patterns of phylogenomic discordance across the murine rodents (Old World mice and rats) - a large and ecologically diverse group that gave rise to the mouse and rat model systems. Combining new linked-read genome assemblies for seven murine species with eleven published rodent genomes, we first used ultra-conserved elements (UCEs) to infer a robust species tree. We then used whole genomes to examine finer-scale patterns of discordance and found that phylogenies built from proximate chromosomal regions had similar phylogenies. However, there was no relationship between tree similarity and local recombination rates in house mice, suggesting that genetic linkage influences phylogenetic patterns over deeper timescales. This signal may be independent of contemporary recombination landscapes. We also detected a strong influence of linked selection whereby purifying selection at UCEs led to less discordance, while genes experiencing positive selection showed more discordant and variable phylogenetic signals. Finally, we show that assuming a single species tree can result in high error rates when testing for positive selection under different models. Collectively, our results highlight the complex relationship between phylogenetic inference and genome biology and underscore how failure to account for this complexity can mislead comparative genomic studies.

Low hybridization temperatures improve target capture success of invertebrate loci: a case study of leaf-footed bugs (Hemiptera: Coreoidea)

Target capture is widely used in phylogenomic, ecological and functional genomic studies. Bait sets that allow capture from a diversity of species can be advantageous, but high-sequence divergence from baits can limit yields. Currently, only four experimental comparisons of a critical target capture parameter, hybridization temperature, have been published. These have been in vertebrates, where bait divergences are typically low, and none include invertebrates where bait-target divergences may be higher. Most invertebrate capture studies use a fixed, high hybridization temperature to maximize the proportion of on-target data, but many report low locus recovery. Using leaf-footed bugs (Hemiptera: Coreoidea), we investigate the effect of hybridization temperature on capture success of ultraconserved elements targeted by (i) baits developed from divergent hemipteran genomes and (ii) baits developed from less divergent coreoid transcriptomes. Lower temperatures generally resulted in more contigs and improved recovery of targets despite a lower proportion of on-target reads, lower read depth and more putative paralogues. Hybridization temperatures had less of an effect when using transcriptome-derived baits, which is probably due to lower bait-target divergences and greater bait tiling density. Thus, accommodating low hybridization temperatures during target capture can provide a cost-effective, widely applicable solution to improve invertebrate locus recovery.

Description of a new Kurixalus species (Rhacophoridae, Anura) and a northwards range extension of the genus

Knowledge of biodiversity before species become extinct is paramount to conservation, especially when the relevant species are far from their expected distribution and, thus, likely overlooked. Here, we describe a new Kurixalus species corresponding to a range extension of Kurixalus on the Asian mainland, with the closest population in Taiwan. The species diverged from its closest relative during the Late Pliocene to Pleistocene, ca. 3.06 Mya (HPD 95%: 5.82-0.01), based on calibrations with a relaxed clock species tree of unlinked mtDNA 12S rRNA and nuclear DNA TYR. The status of the newly-described species is also supported by a divergence in call properties and morphometrics. We named the species described here as Kurixalusinexpectatus sp. nov. due to the nature of the discovery, as well as the adjunct distribution of the species relative to its closest congeners. The species was found in Zhejiang Province and it represents a range extension of 663 km for the Kurixalus genus.

The role of climate and islands in species diversification and reproductive-mode evolution of Old World tree frogs

Large diversifications of species are known to occur unevenly across space and evolutionary lineages, but the relative importance of their driving mechanisms, such as climate, ecological opportunity and key evolutionary innovations (KEI), remains poorly understood. Here, we explore the remarkable diversification of rhacophorid frogs, which represent six percent of global amphibian diversity, utilize four distinct reproductive modes, and span a climatically variable area across mainland Asia, associated continental islands, and Africa. Using a complete species-level phylogeny, we find near-constant diversification rates but a highly uneven distribution of species richness. Montane regions on islands and some mainland regions have higher phylogenetic diversity and unique assemblages of taxa; we identify these as cool-wet refugia. Starting from a centre of origin, rhacophorids reached these distant refugia by adapting to new climatic conditions (‘niche evolution’-dominant), especially following the origin of KEIs such as terrestrial reproduction (in the Late Eocene) or by dispersal during periods of favourable climate (‘niche conservatism’-dominant).

By examining climate, geographical and phylogenetic data, the diversification and evolution of rhacophorid frogs is explored

Phylogenomic assessment prompts recognition of the Serianthes clade and confirms the monophyly of Serianthes and its relationship with Falcataria and Wallaceodendron in the wider ingoid clade (Leguminosae, Caesalpinioideae)

The Indo-Pacific legume genus Serianthes was recently placed in the Archidendron clade (sensu Koenen et al. 2020), a subclade of the mimosoid clade in subfamily Caesalpinioideae, which also includes Acacia, Archidendron, Archidendropsis, Falcataria, Pararchidendron, Paraserianthes and Wallaceodendron. Serianthes comprises ca. 18 species, five subspecies and two varieties that are characterised by bipinnately compound leaves with alternate sessile leaflets, branched axillary corymbiform panicles and woody indehiscent pods. Generic relationships, as well as species relationships within genera in the Archidendron clade, remain uncertain. While the sister relationship between Serianthes and the genus Falcataria is strongly supported by molecular data, the distinction between Serianthes and the monotypic genus Wallaceodendron has been questioned, based on their similar flower and fruit morphologies. We combined three gene-enriched hybrid capture DNA sequence datasets (generated from the 964 mimobaits v1 probe set, the expanded 997 mimobaits v2 probe set and the GoFlag angiosperm 408 probe set) and used their overlapping markers (77 loci of the target exonic and flanking regions) to test the monophyly of Serianthes and to investigate generic relationships within the Archidendron clade using 55 ingoid plus two outgroup taxa. We show that Serianthes is monophyletic, confirm the Serianthes + Falcataria sister relationship to Wallaceodendron and recognise this combined clade as the Serianthes clade within the Archidendron clade. We also evaluated the use of overlapping loci across datasets in combination with concordance analyses to test generic relationships and further investigate previously unresolved relationships across the wider ingoid clade. Concordance analysis revealed limited gene tree conflicts near the tips of the Archidendron clade, but increased discordance at the base of the clade, which could be attributed to rapid lineage divergence (radiation) and/or incomplete lineage sorting.

FrogCap: A modular sequence capture probe-set for phylogenomics and population genetics for all frogs, assessed across multiple phylogenetic scales

Despite the prevalence of high-throughput sequencing in phylogenetics, many relationships remain difficult to resolve because of conflicting signal among genomic regions. Selection of different types of molecular markers from different genomic regions is required to overcome these challenges. For evolutionary studies in frogs, we introduce the publicly available FrogCap suite of genomic resources, which is a large collection of ~15,000 markers that unifies previous genetic sequencing efforts. FrogCap is designed to be modular, such that subsets of markers and SNPs can be selected based on the desired phylogenetic scale. FrogCap uses a variety of marker types that include exons and introns, ultraconserved elements, and previously sequenced Sanger markers, which span up to 10,000 bp in alignment lengths; in addition, we demonstrate potential for SNP-based analyses. We tested FrogCap using 121 samples distributed across five phylogenetic scales, comparing probes designed using a consensus- or exemplar genome-based approach. Using the consensus design is more resilient to issues with sensitivity, specificity, and missing data than picking an exemplar genome sequence. We also tested the impact of different bait kit sizes (20,020 vs. 40,040) on depth of coverage and found triple the depth for the 20,020 bait  kit. We observed sequence capture success (i.e., missing data, sequenced markers/bases, marker length, and informative sites) across phylogenetic scales. The incorporation of different marker types is effective for deep phylogenetic relationships and shallow population genetics studies. Having demonstrated FrogCap's utility and modularity, we conclude that these new resources are efficacious for high-throughput sequencing projects across variable timescales.

Seeing through the hedge: Phylogenomics of Thuja (Cupressaceae) reveals prominent incomplete lineage sorting and ancient introgression for Tertiary relict flora

The Eastern Asia (EA) - North America (NA) disjunction is a well-known biogeographic pattern of the Tertiary relict flora; however, few studies have investigated the evolutionary history of this disjunction using a phylogenomic approach. Here, we used 2369 single copy nuclear genes and nearly full plastomes to reconstruct the evolutionary history of the small Tertiary relict genus Thuja, which consists of five disjunctly distributed species. The nuclear species tree strongly supported an EA clade Thuja standishii-Thuja sutchuenensis and a "disjunct clade", where western NA species T. plicata is sister to an EA-eastern NA disjunct Thuja occidentalis-Thuja koraiensis group. Our results suggested that the observed topological discordance among the gene trees as well as the cytonuclear discordance is mainly due to incomplete lineage sorting, probably facilitated by the fast diversification of Thuja around the Early Miocene and the large effective population sizes of ancestral lineages. Furthermore, approximately 20% of the T. sutchuenensis nuclear genome is derived from an unknown ancestral lineage of Thuja, which might explain the close resemblance of its cone morphology to that of an ancient fossil species. Overall, our study demonstrates that single genes may not resolve interspecific relationships for disjunct taxa, and that more reliable results will come from hundreds or thousands of loci, revealing a more complex evolutionary history. This will steadily improve our understanding of their origin and evolution.

Target-capture phylogenomics provide insights on gene and species tree discordances in Old World treefrogs (Anura: Rhacophoridae)

Genome-scale data have greatly facilitated the resolution of recalcitrant nodes that Sanger-based datasets have been unable to resolve. However, phylogenomic studies continue to use traditional methods such as bootstrapping to estimate branch support; and high bootstrap values are still interpreted as providing strong support for the correct topology. Furthermore, relatively little attention has been given to assessing discordances between gene and species trees, and the underlying processes that produce phylogenetic conflict. We generated novel genomic datasets to characterize and determine the causes of discordance in Old World treefrogs (Family: Rhacophoridae)-a group that is fraught with conflicting and poorly supported topologies among major clades. Additionally, a suite of data filtering strategies and analytical methods were applied to assess their impact on phylogenetic inference. We showed that incomplete lineage sorting was detected at all nodes that exhibited high levels of discordance. Those nodes were also associated with extremely short internal branches. We also clearly demonstrate that bootstrap values do not reflect uncertainty or confidence for the correct topology and, hence, should not be used as a measure of branch support in phylogenomic datasets. Overall, we showed that phylogenetic discordances in Old World treefrogs resulted from incomplete lineage sorting and that species tree inference can be improved using a multi-faceted, total-evidence approach, which uses the most amount of data and considers results from different analytical methods and datasets.