Host-microbe interactions are intimately linked to eukaryotic evolution, particularly in sap-sucking insects that often rely on obligate microbial symbionts for nutrient provisioning. Cicadas (Cicadidae: Auchenorrhyncha) specialize on xylem fluid and derive many essential amino acids and vitamins from intracellular bacteria or fungi (Hodgkinia, Sulcia, and Ophiocordyceps) that are propagated via transmission from mothers to offspring. Despite the beneficial role of these non-gut symbionts in nutrient provisioning, the role of beneficial microbiota within the gut remains unclear. Here, we investigate the relative abundance and impact of host phylogeny and ecology on gut microbial diversity in cicadas using 16S ribosomal RNA gene amplicon sequencing data from 197 wild-collected cicadas and new mitochondrial genomes across 38 New Zealand cicada species, including natural hybrids between one pair of two species. We find low abundance and a lack of phylogenetic structure and hybrid effects but a significant role of elevation in explaining variation in gut microbiota.

Background

Despite considerable effort, progress in spider molecular systematics has lagged behind many other comparable arthropod groups, thereby hindering family-level resolution, classification, and testing of important macroevolutionary hypotheses. Recently, alternative targeted sequence capture techniques have provided molecular systematics a powerful tool for resolving relationships across the Tree of Life. One of these approaches, Anchored Hybrid Enrichment (AHE), is designed to recover hundreds of unique orthologous loci from across the genome, for resolving both shallow and deep-scale evolutionary relationships within non-model systems. Herein we present a modification of the AHE approach that expands its use for application in spiders, with a particular emphasis on the infraorder Mygalomorphae.

Results

Our aim was to design a set of probes that effectively capture loci informative at a diversity of phylogenetic timescales. Following identification of putative arthropod-wide loci, we utilized homologous transcriptome sequences from 17 species across all spiders to identify exon boundaries. Conserved regions with variable flanking regions were then sought across the tick genome, three published araneomorph spider genomes, and raw genomic reads of two mygalomorph taxa. Following development of the 585 target loci in the Spider Probe Kit, we applied AHE across three taxonomic depths to evaluate performance: deep-level spider family relationships (33 taxa, 327 loci); family and generic relationships within the mygalomorph family Euctenizidae (25 taxa, 403 loci); and species relationships in the North American tarantula genus Aphonopelma (83 taxa, 581 loci). At the deepest level, all three major spider lineages (the Mesothelae, Mygalomorphae, and Araneomorphae) were supported with high bootstrap support. Strong support was also found throughout the Euctenizidae, including generic relationships within the family and species relationships within the genus Aptostichus. As in the Euctenizidae, virtually identical topologies were inferred with high support throughout Aphonopelma.

Conclusions

The Spider Probe Kit, the first implementation of AHE methodology in Class Arachnida, holds great promise for gathering the types and quantities of molecular data needed to accelerate an understanding of the spider Tree of Life by providing a mechanism whereby different researchers can confidently and effectively use the same loci for independent projects, yet allowing synthesis of data across independent research groups.

We present the first genome-wide molecular phylogeny of jumping spiders (Araneae: Salticidae), inferred from Anchored Hybrid Enrichment (AHE) sequence data. From 12 outgroups plus 34 salticid taxa representing all but one subfamily and most major groups recognized in previous work, we obtained 447 loci totalling 96,946 aligned nucleotide sites. Our analyses using concatenated likelihood, parsimony, and coalescent methods (ASTRAL and SVDQuartets) strongly confirm most previous results, resolving as monophyletic the Spartaeinae, Salticinae (with the hisponines sister), Salticoida, Amycoida, Saltafresia, and Simonida. The agoriines, previously difficult to place beyond subfamily, are finally placed confidently within the saltafresians as relatives of the chrysillines and hasariines. Relationships among the baviines, astioids, marpissoids, and saltafresians remain uncertain, though our analyses tentatively conclude the first three form a clade together. Deep relationships, among the seven subfamilies, appear to be largely resolved, with spartaeines, lyssomanines, and asemoneines forming a clade. In most analyses, Onomastus (representing the onomastines) is strongly supported as sister to the hisponines plus salticines. Overall, the much-improved resolution of many deep relationships despite a relatively sparse taxon sample suggests AHE is a promising technique for salticid phylogenetics.

Background

The past decade has witnessed remarkable progress towards resolution of the Tree of Life. However, despite the increased use of genomic scale datasets, some phylogenetic relationships remain difficult to resolve. Here we employ anchored phylogenomics to capture 107 nuclear loci in 29 species of acanthomorph teleost fishes, with 25 of these species sampled from the recently delimited clade Ovalentaria. Previous studies employing multilocus nuclear exon datasets have not been able to resolve the nodes at the base of the Ovalentaria tree with confidence. Here we test whether a phylogenomic approach will provide better support for these nodes, and if not, why this may be.

Results

After using a novel method to account for paralogous loci, we estimated phylogenies with maximum likelihood and species tree methods using DNA sequence alignments of over 80,000 base pairs. Several key relationships within Ovalentaria are well resolved, including 1) the sister taxon relationship between Cichlidae and Pholidichthys, 2) a clade containing blennies, grammas, clingfishes, and jawfishes, and 3) monophyly of Atherinomorpha (topminnows, flyingfishes, and silversides). However, many nodes in the phylogeny associated with the early diversification of Ovalentaria are poorly resolved in several analyses. Through the use of rarefaction curves we show that limited phylogenetic resolution among the earliest nodes in the Ovalentaria phylogeny does not appear to be due to a deficiency of data, as average global node support ceases to increase when only 1/3rd of the sampled loci are used in analyses. Instead this lack of resolution may be driven by model misspecification as a Bayesian mixed model analysis of the amino acid dataset provided good support for parts of the base of the Ovalentaria tree.

Conclusions

Although it does not appear that the limited phylogenetic resolution among the earliest nodes in the Ovalentaria phylogeny is due to a deficiency of data, it may be that both stochastic and systematic error resulting from model misspecification play a role in the poor resolution at the base of the Ovalentaria tree as a Bayesian approach was able to resolve some of the deeper nodes, where the other methods failed.

AIM: Rump plumage coloration of the Eastern Yellow Robin (Eopsaltria australis), a widespread Australian songbird, varies from bright yellow in the tropical north to oliveâgreen in the temperate south. Here, we test whether colour variation: (1) correlates most strongly with neutral genetic variation and so is best explained by historical processes, (2) reflects selection associated with different visual environments (dense versus open habitats) and/or (3) reflects selection associated with climatic variation. LOCATION: Eastern Australia. METHODS: We quantified colour variation using reflectance spectrometry and visual models. We performed geographical cline analysis of colour and neutral genetic variation (genomeâwide single nucleotide polymorphisms). We tested for correlations of colour variation with climate, vegetation density, geographical location and genetic variation. We accounted for covariation and spatial autocorrelation, and conducted analyses at continental and regional spatial scales. RESULTS: Clinal variation of colour traits and neutral genetic markers were largely concordant. At the continental scale, colour variation was strongly associated with neutral genetic structure and geography, and to a lesser extent with environment. At the regional scale, environmental variation was a better predictor of colour variation than it was at the larger scale. MAIN CONCLUSION: At the continental scale, colour variation is strongly associated with largeâscale population history. In contrast, at the regional scale, where the influence of history and geography is weaker, environmental variation has a role in facilitating the maintenance of colour variation. Our results highlight the need to assess selective and neutral alternatives at multiple spatial scales when studying geographical variation.