Phylogenetics of Coenonymphina (Nymphalidae: Satyrinae) and the problem of rooting rapid radiations

Chromosome-Level Assembly and Annotation of the Pearly Heath Coenonympha arcania Butterfly Genome

We present the first chromosome-level genome assembly and annotation of the pearly heath Coenonympha arcania, generated with a PacBio HiFi sequencing approach and complemented with Hi-C data. We additionally compare synteny, gene, and repeat content between C. arcania and other Lepidopteran genomes. This reference genome will enable future population genomics studies with Coenonympha butterflies, a species-rich genus that encompasses some of the most highly endangered butterfly taxa in Europe.

Phylogeny of North Americas largest cicada radiation redefines Tibicinoides and Okanagana (Hemiptera: Auchenorrhyncha: Cicadidae: Tibicininae)

Tibicinoides, with three small endemic California cicada species, has a confusing, intertwined systematic history with Okanagana that we unravel here. An ingroup including all species of Tibicinoides and the majority (84.7%) of Okanagana species were sampled for six gene regions, polarized with Clidophleps, Okanagodes, Subpsaltria, and Tibicina outgroups, and subjected to Bayesian phylogenetic analysis. Although the ingroup was monophyletic from all outgroups including Tibicina, Tibicinoides rendered Okanagana paraphyletic among two major ingroup clades. To bring classification into agreement with phylogeny, we redescribe and redefine Tibicinoides to include all Okanagana species with a hooked uncus in the male genitalia, all of which grouped with the type T. cupreosparsa (Uhler, 1889) in the first of these clades: T. boweni (Chatfield-Taylor & Cole, 2020) comb. n., T. catalina (Davis, 1936) comb. n., T. hesperia (Uhler, 1876) comb. n., T. mercedita (Davis, 1915), T. minuta (Davis, 1915), T. pallidula (Davis, 1917a) comb. n., T. pernix (Bliven, 1964) comb. n., T. rubrovenosa (Davis, 1915) comb. n., T. simulata (Davis, 1921) comb. n., T. striatipes (Haldeman, 1852) comb. n., T. uncinata (Van Duzee, 1915) comb. n., T. utahensis (Davis, 1919) comb. n., and T. vanduzeei (Distant, 1914) comb. n. Okanagana is redescribed and restricted to the species of the second major clade which contained the type O. rimosa (Say, 1830). We describe two new genera for morphologically distinct orphan lineages: Chlorocanta gen. nov. for C. viridis (Davis, 1918) comb. n. and Hewlettia gen. nov. for H. nigriviridis (Davis, 1921) comb. n. We recognize O. rubrobasalis Davis, 1926 stat. rev. as a species and relegate two former species to junior subjective synonyms: O. noveboracensis (Emmons, 1854) = O. canadensis (Provancher, 1889) and O. occidentalis (Walker in Lord, 1866) = O. lurida Davis, 1919. Tibicinoides and Okanagana together represent a rapid radiation that presents challenges to phylogenetic analysis including suboptimal outgroups and short internodes.

A review of the genus Triphysa Zeller, 1858 (Lepidoptera, Satyridae)*

A review of the genus Triphysa Zeller, 1858 is presented. One new species Triphysa issykkulica sp.n. (type locality: Kazakhstan, W of Almaty, 800 m) and 8 new subspecies are described: Triphysa phryne kasikoporana ssp. n. (type locality: Kasikoporan [NE Turkey, Agri prov.]), Triphysa striatula urumtchiensis ssp. n. (type locality: Urumtchi), Triphysa issykkulica pljustchi ssp. n. (type locality: W. Kirgiziya, Talasskii Mts., Manas), Triphysa nervosa tuvinica ssp. n. (type locality: N. Tuva, near Kyzyl, Tuge Mt.), Triphysa nervosa arturi ssp. n. (type locality: S. Tuva, 15 km WSW Erzin), Triphysa nervosa kobdoensis ssp. n. (type locality: W. Mongolia, Hovd aimak, 15 km S Khara-Us-Nuur lake, 1300 m), Triphysa nervosa mongolaltaica ssp. n. (type locality: Mongolia, Hovd aimak, Bulgan-Gol basin, middle stream of Ulyasutai-Gol river, 2500−3000 m) and Triphysa nervosa brinikhi ssp. n. (type locality: Russia, Chita Reg., Onon distr., 18 km WSW Nizhniy Zasuchey vill., Butyvken lake, Pinus forest, steppe) are described. New status for Triphysa striatula Elwes, 1899, stat. n. is established. The lectotypes of Triphysa nervosa gartoki O.Bang-Haas, 1927, Triphysa phryne kintschouensis O. Bang-Haas, 1939, Triphysa phryne biocellata Staudinger, 1901, Triphysa nervosa tscherski Grum-Grshimailo, 1889 [1890], Triphysa nervosa glacialis A. Bang-Haas, 1912 are designated, the neotype of T. dohrnii Zeller, 1850 (type locality: [Russia], Sarepta) is designated.

Miocene Climate and Habitat Change Drove Diversification in Bicyclus, Africa's Largest Radiation of Satyrine Butterflies

Compared to other regions, the drivers of diversification in Africa are poorly understood. We studied a radiation of insects with over 100 species occurring in a wide range of habitats across the Afrotropics to investigate the fundamental evolutionary processes and geological events that generate and maintain patterns of species richness on the continent. By investigating the evolutionary history of Bicyclus butterflies within a phylogenetic framework, we inferred the group's origin at the Oligo-Miocene boundary from ancestors in the Congolian rainforests of central Africa. Abrupt climatic fluctuations during the Miocene (ca. 19-17 Ma) likely fragmented ancestral populations, resulting in at least eight early-divergent lineages. Only one of these lineages appears to have diversified during the drastic climate and biome changes of the early Miocene, radiating into the largest group of extant species. The other seven lineages diversified in forest ecosystems during the late Miocene and Pleistocene when climatic conditions were more favourable-warmer and wetter. Our results suggest changing Neogene climate, uplift of eastern African orogens, and biotic interactions might have had different effects on the various subclades of Bicyclus, producing one of the most spectacular butterfly radiations in Africa.

mPartition: A Model-Based Method for Partitioning Alignments

Maximum likelihood (ML) analysis of nucleotide or amino-acid alignments is widely used to infer evolutionary relationships among species. Computing the likelihood of a phylogenetic tree from such alignments is a complicated task because the evolutionary processes typically vary across sites. A number of studies have shown that partitioning alignments into sub-alignments of sites, where each sub-alignment is analyzed using a different model of evolution (e.g., GTR + I + G), is a sensible strategy. Current partitioning methods group sites into subsets based on the inferred rates of evolution at the sites. However, these do not provide sufficient information to adequately reflect the substitution processes of characters at the sites. Moreover, the site rate-based methods group all invariant sites into one subset, potentially resulting in wrong phylogenetic trees. In this study, we propose a partitioning method, called mPartition, that combines not only the evolutionary rates but also substitution models at sites to partition alignments. Analyses of different partitioning methods on both real and simulated datasets showed that mPartition was better than the other partitioning methods tested. Notably, mPartition overcame the pitfall of grouping all invariant sites into one subset. Using mPartition may lead to increased accuracy of ML-based phylogenetic inference, especially for multiple loci or whole genome datasets.

A simple method for data partitioning based on relative evolutionary rates

BACKGROUND

Multiple studies have demonstrated that partitioning of molecular datasets is important in model-based phylogenetic analyses. Commonly, partitioning is done a priori based on some known properties of sequence evolution, e.g. differences in rate of evolution among codon positions of a protein-coding gene. Here we propose a new method for data partitioning based on relative evolutionary rates of the sites in the alignment of the dataset being analysed. The rates are inferred using the previously published Tree Independent Generation of Evolutionary Rates (TIGER), and the partitioning is conducted using our novel python script RatePartitions. We conducted simulations to assess the performance of our new method, and we applied it to eight published multi-locus phylogenetic datasets, representing different taxonomic ranks within the insect order Lepidoptera (butterflies and moths) and one phylogenomic dataset, which included ultra-conserved elements as well as introns.

METHODS

We used TIGER-rates to generate relative evolutionary rates for all sites in the alignments. Then, using RatePartitions, we partitioned the data into partitions based on their relative evolutionary rate. RatePartitions applies a simple formula that ensures a distribution of sites into partitions following the distribution of rates of the characters from the full dataset. This ensures that the invariable sites are placed in a partition with slowly evolving sites, avoiding the pitfalls of previously used methods, such as k-means. Different partitioning strategies were evaluated using BIC scores as calculated by PartitionFinder.

RESULTS

Simulations did not highlight any misbehaviour of our partitioning approach, even under difficult parameter conditions or missing data. In all eight phylogenetic datasets, partitioning using TIGER-rates and RatePartitions was significantly better as measured by the BIC scores than other partitioning strategies, such as the commonly used partitioning by gene and codon position. We compared the resulting topologies and node support for these eight datasets as well as for the phylogenomic dataset.

DISCUSSION

We developed a new method of partitioning phylogenetic datasets without using any prior knowledge (e.g. DNA sequence evolution). This method is entirely based on the properties of the data being analysed and can be applied to DNA sequences (protein-coding, introns, ultra-conserved elements), protein sequences, as well as morphological characters. A likely explanation for why our method performs better than other tested partitioning strategies is that it accounts for the heterogeneity in the data to a much greater extent than when data are simply subdivided based on prior knowledge.

Paleoclimate determines diversification patterns in the fossorial snake family Uropeltidae Cuvier, 1829

Understanding how and why diversification rates vary across evolutionary time is central to understanding how biodiversity is generated and maintained. Recent mathematical models that allow estimation of diversification rates across time from reconstructed phylogenies have enabled us to make inferences on how biodiversity copes with environmental change. Here, we explore patterns of temporal diversification in Uropeltidae, a diverse fossorial snake family. We generate a time-calibrated phylogenetic hypothesis for Uropeltidae and show a significant correlation between diversification rate and paleotemperature during the Cenozoic. We show that the temporal diversification pattern of this group is punctuated by one rate shift event with a decrease in diversification and turnover rate between ca. 11Ma to present, but there is no strong support for mass extinction events. The analysis indicates higher turnover during periods of drastic climatic fluctuations and reduced diversification rates associated with contraction and fragmentation of forest habitats during the late Miocene. Our study highlights the influence of environmental fluctuations on diversification rates in fossorial taxa such as uropeltids, and raises conservation concerns related to present rate of climate change.

Ten genes and two topologies: an exploration of higher relationships in skipper butterflies (Hesperiidae)

Despite multiple attempts to infer the higher-level phylogenetic relationships of skipper butterflies (Family Hesperiidae), uncertainties in the deep clade relationships persist. The most recent phylogenetic analysis included fewer than 30% of known genera and data from three gene markers. Here we reconstruct the higher-level relationships with a rich sampling of ten nuclear and mitochondrial markers (7,726 bp) from 270 genera and find two distinct but equally plausible topologies among subfamilies at the base of the tree. In one set of analyses, the nuclear markers suggest two contrasting topologies, one of which is supported by the mitochondrial dataset. However, another set of analyses suggests mito-nuclear conflict as the reason for topological incongruence. Neither topology is strongly supported, and we conclude that there is insufficient phylogenetic evidence in the molecular dataset to resolve these relationships. Nevertheless, taking morphological characters into consideration, we suggest that one of the topologies is more likely.

Anatomy of the Enigmatic Reptile Elachistosuchus huenei Janensch, 1949 (Reptilia: Diapsida) from the Upper Triassic of Germany and Its Relevance for the Origin of Sauria

The holotype and only known specimen of the enigmatic small reptile Elachistosuchus huenei Janensch, 1949 from the Upper Triassic (Norian) Arnstadt Formation of Saxony-Anhalt (Germany) is redescribed using μCT scans of the material. This re-examination revealed new information on the morphology of this taxon, including previously unknown parts of the skeleton such as the palate, braincase, and shoulder girdle. Elachistosuchus is diagnosed especially by the presence of the posterolateral process of the frontal, the extension of the maxillary tooth row to the posterior margin of the orbit, the free posterior process of the jugal, and the notched anterior margin of the interclavicle. Phylogenetic analyses using two recently published character-taxon matrices recovered conflicting results for the phylogenetic position of Elachistosuchus-either as an archosauromorph, as a lepidosauromorph or as a more basal, non-saurian diapsid. These different placements highlight the need of a thorough revision of critical taxa and new character sets used for inferring neodiapsid relationships.

Diversity dynamics in Nymphalidae butterflies: effect of phylogenetic uncertainty on diversification rate shift estimates

The species rich butterfly family Nymphalidae has been used to study evolutionary interactions between plants and insects. Theories of insect-hostplant dynamics predict accelerated diversification due to key innovations. In evolutionary biology, analysis of maximum credibility trees in the software MEDUSA (modelling evolutionary diversity using stepwise AIC) is a popular method for estimation of shifts in diversification rates. We investigated whether phylogenetic uncertainty can produce different results by extending the method across a random sample of trees from the posterior distribution of a Bayesian run. Using the MultiMEDUSA approach, we found that phylogenetic uncertainty greatly affects diversification rate estimates. Different trees produced diversification rates ranging from high values to almost zero for the same clade, and both significant rate increase and decrease in some clades. Only four out of 18 significant shifts found on the maximum clade credibility tree were consistent across most of the sampled trees. Among these, we found accelerated diversification for Ithomiini butterflies. We used the binary speciation and extinction model (BiSSE) and found that a hostplant shift to Solanaceae is correlated with increased net diversification rates in Ithomiini, congruent with the diffuse cospeciation hypothesis. Our results show that taking phylogenetic uncertainty into account when estimating net diversification rate shifts is of great importance, as very different results can be obtained when using the maximum clade credibility tree and other trees from the posterior distribution.

Host plant utilization, host range oscillations and diversification in nymphalid butterflies: a phylogenetic investigation

It has been suggested that phenotypic plasticity is a major factor in the diversification of life, and that variation in host range in phytophagous insects is a good model for investigating this claim. We explore the use of angiosperm plants as hosts for nymphalid butterflies, and in particular the evidence for past oscillations in host range and how they are linked to host shifts and to diversification. At the level of orders of plants, a relatively simple pattern of host use and host shifts emerges, despite the 100 million years of history of the family Nymphalidae. We review the evidence that these host shifts and the accompanying diversifications were associated with transient polyphagous stages, as suggested by the "oscillation hypothesis." In addition, we investigate all currently polyphagous nymphalid species and demonstrate that the state of polyphagy is rare, has a weak phylogenetic signal, and a very apical distribution in the phylogeny; we argue that these are signs of its transient nature. We contrast our results with data from the bark beetles Dendroctonus, in which a more specialized host use is instead the apical state. We conclude that plasticity in host use is likely to have contributed to diversification in nymphalid butterflies.

Unraveling reticulate evolution in North American Dryopteris (Dryopteridaceae)

BACKGROUND

The thirteen species of Dryopteris in North America have long been suspected of having undergone a complicated history of reticulate evolution via allopolyploid hybridization. Various explanations for the origins of the allopolyploid taxa have been suggested, and though most lines of evidence have supported the so-called "semicristata" hypothesis, contention over the group's history has continued in several recent, conflicting studies.

RESULTS

Sequence data from nine plastid and two nuclear markers were collected from 73 accessions representing 35 species of Dryopteris. Sequences from each of the allopolyploids are most closely related to their progenitor species as predicted by the "semicristata" hypothesis. Allotetraploid D. campyloptera appears to be derived from a hybrid between diploid D. expansa and D. intermedia; D. celsa, from diploid D. ludoviciana x D. goldiana; and D. carthusiana and D. cristata, from diploid "D. semicristata" x D. intermedia and D. ludoviciana, respectively. Allohexaploid D. clintoniana appears to be derived from D. cristata x D.goldiana. The earliest estimated dates of formation of the allopolyploids, based on divergence time analyses, were within the last 6 Ma. We found no evidence for recurrent formation of any of the allopolyploids. The sexual allopolyploid taxa are derived from crosses between parents that show intermediate levels of genetic divergence relative to all pairs of potential progenitors. In addition, the four allotetraploids are transgressive with respect to geographic range relative to one or both of their parents (their ranges extend beyond those of the parents), suggesting that ecological advantages in novel habitats or regions may promote long-term regional coexistence of the hybrid taxa with their progenitors.

CONCLUSIONS

This study provides the first thorough evaluation of the North American complex of woodferns using extensive sampling of taxa and genetic markers. Phylogenies produced from each of three datasets (one plastid and two nuclear) support the "semicristata" hypothesis, including the existence of a missing diploid progenitor, and allow us to reject all competing hypotheses. This study demonstrates the value of using multiple, biparentally inherited markers to evaluate reticulate complexes, assess the frequency of recurrent polyploidization, and determine the relative importance of introgression vs. hybridization in shaping the histories of such groups.

Phylogenetic structure in African HIV-1 subtype C revealed by selective sequential pruning

Subtype C is the most prevalent clade in the HIV-1 pandemic. Previous studies suggested that African HIV-1 subtype C lacks a well-defined phylogenetic structure. Here we show that, by sequential pruning of ambiguously positioned taxa, a well-defined intrasubtype C phylogenetic structure becomes apparent, with 52% African HIV-1 subtype C isolates analyzed in envelope sequences branching within 11 clusters, also supported in a tree of full-length genomes, and all with geographical associations. Among 46 viruses recently transmitted in South Africa, 70% branched within 7 clusters (41% in the largest one) and 15% additional isolates were intercluster recombinants. Choice of the outgroup sequence and inclusion of intrasubtype recombinant viruses in the analyses could greatly affect support of clusters. The identification of clusters comprising a large proportion of African HIV-1 subtype C viruses may have implications for the design of vaccines intended for use in areas where subtype C is prevalent.

Phylogenetics and biogeography of a spectacular Old World radiation of butterflies: the subtribe Mycalesina (Lepidoptera: Nymphalidae: Satyrini)

Background

Butterflies of the subtribe Mycalesina (Nymphalidae: Satyrinae) are important model organisms in ecology and evolution. This group has radiated spectacularly in the Old World tropics and presents an exciting opportunity to better understand processes of invertebrate rapid radiations. However, the generic-level taxonomy of the subtribe has been in a constant state of flux, and relationships among genera are unknown. There are six currently recognized genera in the group. Mycalesis, Lohora and Nirvanopsis are found in the Oriental region, the first of which is the most speciose genus among mycalesines, and extends into the Australasian region. Hallelesis and Bicyclus are found in mainland Africa, while Heteropsis is primarily Madagascan, with a few species in Africa. We infer the phylogeny of the group with data from three genes (total of 3139 bp) and use these data to reconstruct events in the biogeographic history of the group.

Results

The results indicate that the group Mycalesina radiated rapidly around the Oligocene-Miocene boundary. Basal relationships are unresolved, but we recover six well-supported clades. Some species of Mycalesis are nested within a primarily Madagascan clade of Heteropsis, while Nirvanopsis is nested within Lohora. The phylogeny suggests that the group had its origin either in Asia or Africa, and diversified through dispersals between the two regions, during the late Oligocene and early Miocene. The current dataset tentatively suggests that the Madagascan fauna comprises two independent radiations. The Australasian radiation shares a common ancestor derived from Asia. We discuss factors that are likely to have played a key role in the diversification of the group.

Conclusions

We propose a significantly revised classification scheme for Mycalesina. We conclude that the group originated and radiated from an ancestor that was found either in Asia or Africa, with dispersals between the two regions and to Australasia. Our phylogeny paves the way for further comparative studies on this group that will help us understand the processes underlying diversification in rapid radiations of invertebrates.