Phylogenomic species delimitation and host‐symbiont coevolution in the fungus‐farming ant genus Sericomyrmex Mayr (Hymenoptera: Formicidae): ultraconserved elements (UCEs) resolve a recent radiation

Ants in the Neotropical genus Sericomyrmex Mayr cultivate fungi for food. Both ants and fungi are obligate, coevolved symbionts. The taxonomy of Sericomyrmex is problematic because the morphology of the worker caste is generally homogeneous across all of the species within the genus, species limits are vague, and the relationships between them are unknown. We used ultraconserved elements (UCEs) as genome‐scale markers to reconstruct evolutionary history and to infer species boundaries in Sericomyrmex. We recovered an average of ～990 UCE loci for 88 Sericomyrmex samples from across the geographical range of the genus as well as for five outgroup taxa. Using maximum likelihood and species‐tree approaches, we recovered nearly identical topologies across datasets with 50–95% matrix completeness. We identify nine species‐level lineages in Sericomyrmex, including two new species. This is less than the previously described 19 species, even accounting for two species for which we had no UCE samples, which brings the total number of Sericomyrmex species to 11. Divergence‐dating analyses recovered 4.3 Ma as the crown‐group age estimates for Sericomyrmex, indicating a recent, rapid radiation. We also sequenced mitochondrial cytochrome oxidase subunit I (COI) for 125 specimens. Resolution and support for clades in our COI phylogeny are weak, indicating that COI is not an appropriate species‐delimitation tool. However, taxa within species consistently cluster together, suggesting that COI is useful as a species identification (‘DNA barcoding’) tool. We also sequenced internal transcribed spacer (ITS) and large subunit (LSU) for 32 Sericomyrmex fungal cultivars. The fungal phylogeny confirms that Sericomyrmex fungi are generalized higher‐attine cultivars, interspersed with Trachymyrmex‐associated fungal species, indicating cultivar sharing and horizontal transfer between these two genera. Our results indicate that UCEs offer immense potential for delimiting and resolving relationships of problematic, recently diverged species.     

Advancing mite phylogenomics: Designing ultraconserved elements for Acari phylogeny

Mites (Acari) are one of the most diverse groups of life on Earth; yet, their evolutionary relationships are poorly understood. Also, the resolution of broader arachnid phylogeny has been hindered by an underrepresentation of mite diversity in phylogenomic analyses. To further our understanding of Acari evolution, we design targeted ultraconserved genomic elements (UCEs) probes, intended for resolving the complex relationships between mite lineages and closely related arachnids. We then test our Acari UCE baits in‐silico by constructing a phylogeny using 13 existing Acari genomes, as well as 6 additional taxa from a variety of genomic sources. Our Acari‐specific probe kit improves the recovery of loci within mites over an existing general arachnid UCE probe set. Our initial phylogeny recovers the major mite lineages, yet finds mites to be non‐monophyletic overall, with Opiliones (harvestmen) and Ricinuleidae (hooded tickspiders) rendering Parasitiformes paraphyletic.     

Phylogenomic analysis suggests Coreidae and Alydidae (Hemiptera: Heteroptera) are not monophyletic

Next‐generation sequencing technologies (NGS) allow systematists to amass a wealth of genomic data from non‐model species for phylogenetic resolution at various temporal scales. However, phylogenetic inference for many lineages dominated by non‐model species has not yet benefited from NGS, which can complement Sanger sequencing studies. One such lineage, whose phylogenetic relationships remain uncertain, is the diverse, agriculturally important and charismatic Coreoidea (Hemiptera: Heteroptera). Given the lack of consensus on higher‐level relationships and the importance of a robust phylogeny for evolutionary hypothesis testing, we use a large data set comprised of hundreds of ultraconserved element (UCE) loci to infer the phylogeny of Coreoidea (excluding Stenocephalidae and Hyocephalidae), with emphasis on the families Coreidae and Alydidae. We generated three data sets by including alignments that contained loci sampled for at least 50%, 60%, or 70% of the total taxa, and inferred phylogeny using maximum likelihood and summary coalescent methods. Twenty‐six external morphological features used in relatively comprehensive phylogenetic analyses of coreoids were also re‐evaluated within our molecular phylogenetic framework. We recovered 439–970 loci per species (16%–36% of loci targeted) and combined this with previously generated UCE data for 12 taxa. All data sets, regardless of analytical approach, yielded topologically similar and strongly supported trees, with the exception of outgroup relationships and the position of Hydarinae. We recovered a monophyletic Coreoidea, with Rhopalidae highly supported as the sister group to Alydidae + Coreidae. Neither Alydidae nor Coreidae were monophyletic; the coreid subfamilies Hydarinae and Pseudophloeinae were recovered as more closely related to Alydidae than to other coreid subfamilies. Coreinae were paraphyletic with respect to Meropachyinae. Most morphological traits were homoplastic with several clades defined by few, if any, synapomorphies. Our results demonstrate the utility of phylogenomic approaches in generating robust hypotheses for taxa with long‐standing phylogenetic problems and highlight that novel insights may come from such approaches.     

A multilocus phylogeny of archiheterodont bivalves (Mollusca,Bivalvia, Archiheterodonta)

The bivalve clade Heterodonta encompasses more than half of the extant bivalve species and is presently considered a derived group of the modern bivalves (Newell 1965 ; Waller 1998 ). Heterodonta is subdivided into two major lineages, the hyperdiverse Euheterodonta and Archiheterodonta. The latter comprises four relatively small extant families: Astartidae, Carditidae, Condylocardiidae and Crassatellidae, whose relationships and internal phylogeny are poorly understood. We assessed the phylogeny of archiheterodont bivalves using a multilocus data set comprised of molecular sequence data from six loci (18S rRNA, 28S rRNA, cytochrome c oxidase subunit I, cytochrome b, internal transcribed spacer 2 and histone H3). Resultant data sets of ~4 Kb of concatenated molecular sequence data were analysed using probabilistic approaches (maximum likelihood and Bayesian inference) and parsimony direct optimization. We recovered strong support for the monophyly of Archiheterodonta, within which Astartidae is the sister group of Crassatellidae, and these two constitute the sister clade of Carditidae, which is paraphyletic with respect to Condylocardiidae. The relationships among the constituent species groups were evaluated in the context of the archiheterodont fossil record through the estimation of divergence times. Diversification times of archiheterodont families were congruent with bounded estimates of origins based on palaeontological data: Archiheterodonta diversified during the Devonian, 373.1 Ma (95% highest posterior density interval [HPD] 325.8–428.2); Crassatelloidea around the Carboniferous, 330.1 Ma (95% HPD 291.0–372.7); Crassatellidae around the Triassic, 224.0 (95% HPD 140.6–320.2); Astartidae around the Permian, 288.2 Ma (95% HPD 269.2–307.3); and Carditoidea around the Jurassic, 178.8 Ma (95% HPD 120.9–228.3).     

Ultraconserved elements show utility in phylogenetic inference of Adephaga (Coleoptera) and suggest paraphyly of ‘Hydradephaga’

The beetle suborder Adephaga has been the subject of many phylogenetic reconstructions utilizing a variety of data sources and inference methods. However, no strong consensus has yet emerged on the relationships among major adephagan lineages. Ultraconserved elements (UCEs) have proved useful for inferring difficult or unresolved phylogenies at varying timescales in vertebrates, arachnids and Hymenoptera. Recently, a UCE bait set was developed for Coleoptera using polyphagan genomes and a member of the order Strepsiptera as an outgroup. Here, we examine the utility of UCEs for reconstructing the phylogeny of adephagan families, in the first in vitro application a UCE bait set in Coleoptera. Our final dataset included 305 UCE loci for 18 representatives of all adephagan families except Aspidytidae, and two polyphagan outgroups, with a total concatenated length of 83 547 bp. We inferred trees using maximum likelihood analyses of the concatenated UCE alignment and coalescent species tree methods (astral ii , ASTRID, svdquartets ). Although the coalescent species tree methods had poor resolution and weak support, concatenated analyses produced well‐resolved, highly supported trees. Hydradephaga was recovered as paraphyletic, with Gyrinidae sister to Geadephaga and all other adephagans. Haliplidae was recovered as sister to Dytiscoidea, with Hygrobiidae and Amphizoidae successive sisters to Dytiscidae. Finally, Noteridae was recovered as monophyletic and sister to Meruidae. Given the success of UCE data for resolving phylogenetic relationships within Adephaga, we suggest the potential for further resolution of relationships within Adephaga using UCEs with improved taxon sampling, and by developing Adephaga‐specific probes.     

Geographic mode of speciation in a mountain specialist Avian family endemic to the Palearctic

Mountains host greater avian diversity than lowlands at the same latitude due to their greater diversity of habitats stratified along an elevation gradient. Here we test whether this greater ecological heterogeneity promotes sympatric speciation. We selected accentors (Prunellidae), an avian family associated with mountains of the Palearctic, as a model system. Accentors differ in their habitat/elevation preferences and south‐central Siberia and Himalayan regions each host 6 of the 13 species in the family. We used sequences of the mtDNA ND2 gene and the intron 9 of the Z chromosome specific ACO1 gene to reconstruct a complete species‐level phylogeny of Prunellidae. The tree based on joint analysis of both loci was used to reconstruct the family's biogeographic history and to date the diversification events. We also analyzed the relationship between the node age and sympatry, to determine the geographic mode of speciation in Prunellidae. Our data suggest a Miocene origin of Prunellidae in the Himalayan region. The major division between alpine species (subgenus Laiscopus) and species associated with shrubs (subgenus Prunella) and initial diversification events within the latter happened within the Himalayan region in the Miocene and Pliocene. Accentors colonized other parts of the Palearctic during the Pliocene‐Pleistocene transition. This spread across the Palearctic resulted in rapid diversification of accentors. With only a single exception dating to 0.91 Ma, lineages younger than 1.5 Ma are allopatric. In contrast, sympatry values for older nodes are >0. There was no relationship between node age and range symmetry. Allopatric speciation (not to include peripatric) is the predominant geographic mode of speciation in Prunellidae despite the favorable conditions for ecological diversification in the mountains and range overlaps among species.     

Identification of CR1 retroposons in Arborophila rufipectus and their application to Phasianidae phylogeny

Chicken repeat 1 (CR1), a member of non‐LTR retroposon, is an important phylogenetic marker in avian systematics. In this study, we reported several characteristics of CR1 elements in a draft genome of Arborophila rufipectus (Sichuan partridge). According to the analyses of RepeatMasker, approximately 254 966 CR1 elements were identified in A. rufipectus, covering 6.7% of the genome. Subsequently, we selected eighteen novel CR1 elements by comparing the chicken genome, turkey genome and assembled A. rufipectus scaffolds. Here, a combined data set comprising of 22 CR1 loci, mitochondrial genomes and eight unlinked introns was analysed to infer the evolutionary relationships of twelve Phasianidae species. The applicability of CR1 sequences for inferring avian phylogeny relative to mtDNA and intron sequences was investigated as well. Our results elucidated the position of A. rufipectus in Phasianidae with robust supports that it presented a sister clade to Arborophila ardens/Arborophila brunneopectus, and implied that genus Arborophila was in a basal phylogenetic position within Phasianidae and a phylogenetic affinity between Meleagris gallopavo and Pucrasia macrolopha. Therefore, this work not only resolved some of the confounding relationships among Phasianidae, but also suggested CR1 sequences could provide powerful complementary data for phylogeny reconstruction.     

Molecular phylogeny of the <Emphasis Type="Italic">Drosophila obscura</Emphasis> species group,with emphasis on the Old World species

Background  

Species of the Drosophila obscura species group (e.g., D. pseudoobscura, D. subobscura) have served as favorable models in evolutionary studies since the 1930's. Despite numbers of studies conducted with varied types of data, the basal phylogeny in this group is still controversial, presumably owing to not only the hypothetical 'rapid radiation' history of this group, but also limited taxon sampling from the Old World (esp. the Oriental and Afrotropical regions). Here we reconstruct the phylogeny of this group by using sequence data from 6 loci of 21 species (including 16 Old World ones) covering all the 6 subgroups of this group, estimate the divergence times among lineages, and statistically test the 'rapid radiation' hypothesis.     

Historical biogeography of the ancient lycophyte genus Selaginella: early adaptation to xeric habitats on Pangea

The ancient and cosmopolitan lycophyte genus Selaginella has living representatives around the world, but their historical biogeography has not been assessed with modern methods. We estimated a time‐calibrated phylogeny using DNA marker regions rbcL and ITS1‐5.8S‐ITS2 from 200 species. Node density analyses revealed that Selaginellaceae has significantly older median and mean node ages than other putative “ancient” families. We used statistical model comparison to assess different biogeographical models on our dated tree, and to estimate ancestral ranges. These revealed that Selaginella originated on Euramerica around 383 Ma in the Devonian period, while its peak diversification began with the formation of Pangea. The divergence of the two main species‐rich Selaginella lineages occurred approximately 318 Ma on the supercontinent. The major divergences within these main lineages of Selaginella took place in the Late Permian and Early Triassic, along with lineages highly adapted for xeric habitats on Pangea.     

Ultraconserved elements anchor thousands of genetic markers spanning multiple evolutionary timescales

Although massively parallel sequencing has facilitated large-scale DNA sequencing, comparisons among distantly related species rely upon small portions of the genome that are easily aligned. Methods are needed to efficiently obtain comparable DNA fragments prior to massively parallel sequencing, particularly for biologists working with non-model organisms. We introduce a new class of molecular marker, anchored by ultraconserved genomic elements (UCEs), that universally enable target enrichment and sequencing of thousands of orthologous loci across species separated by hundreds of millions of years of evolution. Our analyses here focus on use of UCE markers in Amniota because UCEs and phylogenetic relationships are well-known in some amniotes. We perform an in silico experiment to demonstrate that sequence flanking 2030 UCEs contains information sufficient to enable unambiguous recovery of the established primate phylogeny. We extend this experiment by performing an in vitro enrichment of 2386 UCE-anchored loci from nine, non-model avian species. We then use alignments of 854 of these loci to unambiguously recover the established evolutionary relationships within and among three ancient bird lineages. Because many organismal lineages have UCEs, this type of genetic marker and the analytical framework we outline can be applied across the tree of life, potentially reshaping our understanding of phylogeny at many taxonomic levels.     

Phylogeny of the symphytan grade of Hymenoptera: new pieces into the old jigsaw(fly) puzzle

The Hymenoptera constitutes one of the largest, and ecologically and economically most important, insect orders. During the past decade, a number of hypotheses on the phylogenetic relationships among hymenopteran families and superfamilies have been presented, based on analyses of molecular and/or morphological data. Nevertheless, many questions still remain, particularly concerning relationships within the hyperdiverse suborder Apocrita, but also when it comes to the evolutionary history of the ancestrally herbivorous “sawfly” lineages that form the basal, paraphyletic grade Symphyta. Because a large part of the uncertainty appears to stem from limited molecular and taxonomic sampling, we set out to investigate the phylogeny of Hymenoptera using nine protein‐coding genes, of which five are new to analyses of the order. In addition, we more than tripled the taxon coverage across the symphytan grade, introducing representatives for many previously unsampled lineages. We recover a well supported phylogenetic structure for these early herbivorous hymenopteran clades, with new information regarding the monophyly of Xyelidae, the placement of the superfamily Pamphilioidea as sister to Tenthredinoidea + Unicalcarida, as well as the interrelationships among the tenthredinoid families Tenthredinidae, Cimbicidae, and Diprionidae. Based on the obtained phylogenies, and to prevent paraphyly of Tenthredinidae, we propose erection of the tribe Heptamelini to family status (Heptamelidae). In particular, our results give new insights into subfamilial relationships within the Tenthredinidae and other species‐rich sawfly families. The expanded gene set provides a useful toolbox for future detailed analyses of symphytan subgroups, especially within the diverse superfamily Tenthredinoidea.     

Occurrence of Can-SINEs and intron sequence evolution supports robust phylogeny of pinniped carnivores and their terrestrial relatives

Investigating the dog genome we found 178 965 introns with a moderate length of 200–1000 bp. A screening of these sequences against 23 different repeat libraries to find insertions of short interspersed elements (SINEs) detected 45 276 SINEs. Virtually all of these SINEs (98%) belong to the tRNA-derived Can-SINE family. Can-SINEs arose about 55 million years ago before Carnivora split into two basal groups, the Caniformia (dog-like carnivores) and the Feliformia (cat-like carnivores). Genome comparisons of dog and cat recovered 506 putatively informative SINE loci for caniformian phylogeny. In this study we show how to use such genome information of model organisms to research the phylogeny of related non-model species of interest. Investigating a dataset including representatives of all major caniformian lineages, we analysed 24 randomly chosen loci for 22 taxa. All loci were amplifiable and revealed 17 parsimony-informative SINE insertions. The screening for informative SINE insertions yields a large amount of sequence information, in particular of introns, which contain reliable phylogenetic information as well. A phylogenetic analysis of intron- and SINE sequence data provided a statistically robust phylogeny which is congruent with the absence/presence pattern of our SINE markers. This phylogeny strongly supports a sistergroup relationship of Musteloidea and Pinnipedia. Within Pinnipedia, we see strong support from bootstrapping and the presence of a SINE insertion for a sistergroup relationship of the walrus with the Otariidae.     

Squamate Conserved Loci (SqCL): A unified set of conserved loci for phylogenomics and population genetics of squamate reptiles

The identification of conserved loci across genomes, along with advances in target capture methods and high‐throughput sequencing, has helped spur a phylogenomics revolution by enabling researchers to gather large numbers of homologous loci across clades of interest with minimal upfront investment in locus design. Target capture for vertebrate animals is currently dominated by two approaches—anchored hybrid enrichment (AHE) and ultraconserved elements (UCE)—and both approaches have proven useful for addressing questions in phylogenomics, phylogeography and population genomics. However, these two sets of loci have minimal overlap with each other; moreover, they do not include many traditional loci that that have been used for phylogenetics. Here, we combine across UCE, AHE and traditional phylogenetic gene locus sets to generate the Squamate Conserved Loci set, a single integrated probe set that can generate high‐quality and highly complete data across all three loci types. We use these probes to generate data for 44 phylogenetically disparate taxa that collectively span approximately 33% of terrestrial vertebrate diversity. Our results generated an average of 4.29 Mb across 4709 loci per individual, of which an average of 2.99 Mb was sequenced to high enough coverage (≥10×) to use for population genetic analyses. We validate the utility of these loci for both phylogenomic and population genomic questions, provide a comparison among these locus sets of their relative usefulness and suggest areas for future improvement.