Whole-genome data reveal the complex history of a diverse ecological community

External insect gall morphology influences the functional guilds of natural enemy communities

The evolution of diverse and novel morphological traits is poorly understood, especially how symbiotic interactions can drive these adaptations. The extreme diversity of external traits in insect-induced galls is currently explained by the Enemy Hypothesis, in which these traits have selective advantage in deterring parasitism. While previous tests of this hypothesis used only taxonomic identity, we argue that ecologically functional traits of enemies (i.e. mode of parasitism, larval development strategy) are a crucial addition. Here, we characterize parasitoid guild composition across four disparate gall systems and find consistent patterns of association between enemy guild and gall morphology. Specifically, galls with a longer average larva-to-surface distance host a significantly higher proportion of enemies with a distinct combination of functional traits (i.e. ectoparasitic, idiobiont, elongate ovipositor). Our results support the Enemy Hypothesis and highlight the importance of species ecology in examining insect communities and the evolution of novel defensive characters.

Mapping restricted introgression across the genomes of admixed indigenous African cattle breeds

BACKGROUND

The genomes of indigenous African cattle are composed of components with Middle Eastern (taurine) and South Asian (indicine) origins, providing a valuable model to study hybridization and to identify genetic barriers to gene flow. In this study, we analysed indigenous African cattle breeds as models of hybrid zones, considering taurine and indicine samples as ancestors. In a genomic cline analysis of whole-genome sequence data, we considered over 8 million variants from 144 animals, which allows for fine-mapping of potential genomic incompatibilities at high resolution across the genome.

RESULTS

We identified several thousand variants that had significantly steep clines ('SCV') across the whole genome, indicating restricted introgression. Some of the SCV were clustered into extended regions, with the longest on chromosome 7, spanning 725 kb and including 27 genes. We found that variants with a high phenotypic impact (e.g. indels, intra-genic and missense variants) likely represent greater genetic barriers to gene flow. Furthermore, our findings provide evidence that a large proportion of breed differentiation in African cattle could be linked to genomic incompatibilities and reproductive isolation. Functional evaluation of genes with SCV suggest that mitonuclear incompatibilities and genes associated with fitness (e.g. resistance to paratuberculosis) could account for restricted gene flow in indigenous African cattle.

CONCLUSIONS

To our knowledge, this is the first time genomic cline analysis has been applied to identify restricted introgression in the genomes of indigenous African cattle and the results provide extended insights into mechanisms (e.g. genomic incompatibilities) contributing to hybrid differentiation. These results have important implications for our understanding of genetic incompatibilities and reproductive isolation and provide important insights into the impact of cross-breeding cattle with the aim of producing offspring that are both hardy and productive.

Inferring the ecological and evolutionary determinants of community genetic diversity

Understanding the relative contributions of ecological and evolutionary processes to the structuring of ecological communities is needed to improve our ability to predict how communities may respond to future changes in an increasingly human-modified world. Metabarcoding methods make it possible to gather population genetic data for all species within a community, unlocking a new axis of data to potentially unveil the origins and maintenance of biodiversity at local scales. Here, we present a new eco-evolutionary simulation model for investigating community assembly dynamics using metabarcoding data. The model makes joint predictions of species abundance, genetic variation, trait distributions and phylogenetic relationships under a wide range of parameter settings (e.g. high speciation/low dispersal or vice versa) and across a range of community states, from pristine and unmodified to heavily disturbed. We first demonstrate that parameters governing metacommunity and local community processes leave detectable signatures in simulated biodiversity data axes. Next, using a simulation-based machine learning approach we show that neutral and non-neutral models are distinguishable and that reasonable estimates of several model parameters within the local community can be obtained using only community-scale genetic data, while phylogenetic information is required to estimate those describing metacommunity dynamics. Finally, we apply the model to soil microarthropod metabarcoding data from the Troodos mountains of Cyprus, where we find that communities in widespread forest habitats are structured by neutral processes, while high-elevation and isolated habitats act as an abiotic filter generating non-neutral community structure. We implement our model within the ibiogen R package, a package dedicated to the investigation of island, and more generally community-scale, biodiversity using community-scale genetic data.

Genetic diversity varies with species traits and latitude in predatory soil arthropods (Myriapoda: Chilopoda)

Aim

To investigate the drivers of intra-specific genetic diversity in centipedes, a group of ancient predatory soil arthropods.

Location

Asia, Australasia and Europe.

Time Period

Present.

Major Taxa Studied

Centipedes (Class: Chilopoda).

Methods

We assembled a database of 1245 mitochondrial cytochrome c oxidase subunit I sequences representing 128 centipede species from all five orders of Chilopoda. This sequence dataset was used to estimate genetic diversity for centipede species and compare its distribution with estimates from other arthropod groups. We studied the variation in centipede genetic diversity with species traits and biogeography using a beta regression framework, controlling for the effect of shared evolutionary history within a family.

Results

A wide variation in genetic diversity across centipede species (0-0.1713) falls towards the higher end of values among arthropods. Overall, 27.57% of the variation in mitochondrial COI genetic diversity in centipedes was explained by a combination of predictors related to life history and biogeography. Genetic diversity decreased with body size and latitudinal position of sampled localities, was greater in species showing maternal care and increased with geographic distance among conspecifics.

Main Conclusions

Centipedes fall towards the higher end of genetic diversity among arthropods, which may be related to their long evolutionary history and low dispersal ability. In centipedes, the negative association of body size with genetic diversity may be mediated by its influence on local abundance or the influence of ecological strategy on long-term population history. Species with maternal care had higher genetic diversity, which goes against expectations and needs further scrutiny. Hemispheric differences in genetic diversity can be due to historic climatic stability and lower seasonality in the southern hemisphere. Overall, we find that despite the differences in mean genetic diversity among animals, similar processes related to life-history strategy and biogeography are associated with the variation within them.

Delimiting the cryptic diversity and host preferences of Sycophila parasitoid wasps associated with oak galls using phylogenomic data

Cryptic species diversity is a major challenge for the species-rich community of parasitoids attacking oak gall wasps due to a high degree of sexual dimorphism, morphological plasticity, small size, and poorly known biology. As such, we know very little about the number of species present, nor the evolutionary forces responsible for generating this diversity. One hypothesis is that trait diversity in the gall wasps, including the morphology of the galls they induce, has evolved in response to selection imposed by the parasitoid community, with reciprocal selection driving diversification of the parasitoids. Using a rare, continental-scale data set of Sycophila parasitoid wasps reared from 44 species of cynipid galls from 18 species of oak across the US, we combined mitochondrial DNA barcodes, Ultraconserved Elements (UCEs), morphological, and natural history data to delimit putative species. Using these results, we generate the first large-scale assessment of ecological specialization and host association in this species-rich group, with implications for evolutionary ecology and biocontrol. We find most Sycophila target specific subsets of available cynipid host galls with similar morphologies, and generally attack larger galls. Our results suggest that parasitoid wasps such as Sycophila have adaptations allowing them to exploit particular host trait combinations, while hosts with contrasting traits are resistant to attack. These findings support the tritrophic niche concept for the structuring of plant-herbivore-parasitoid communities.

Why sequence all eukaryotes?

Life on Earth has evolved from initial simplicity to the astounding complexity we experience today. Bacteria and archaea have largely excelled in metabolic diversification, but eukaryotes additionally display abundant morphological innovation. How have these innovations come about and what constraints are there on the origins of novelty and the continuing maintenance of biodiversity on Earth? The history of life and the code for the working parts of cells and systems are written in the genome. The Earth BioGenome Project has proposed that the genomes of all extant, named eukaryotes-about 2 million species-should be sequenced to high quality to produce a digital library of life on Earth, beginning with strategic phylogenetic, ecological, and high-impact priorities. Here we discuss why we should sequence all eukaryotic species, not just a representative few scattered across the many branches of the tree of life. We suggest that many questions of evolutionary and ecological significance will only be addressable when whole-genome data representing divergences at all of the branchings in the tree of life or all species in natural ecosystems are available. We envisage that a genomic tree of life will foster understanding of the ongoing processes of speciation, adaptation, and organismal dependencies within entire ecosystems. These explorations will resolve long-standing problems in phylogenetics, evolution, ecology, conservation, agriculture, bioindustry, and medicine.

Concordant phylogeographic responses to large-scale coastal disturbance in intertidal macroalgae and their epibiota

Major ecological disturbance events can provide opportunities to assess multispecies responses to upheaval. In particular, catastrophic disturbances that regionally extirpate habitat-forming species can potentially influence the genetic diversity of large numbers of codistributed taxa. However, due to the rarity of such disturbance events over ecological timeframes, the genetic dynamics of multispecies recolonization processes have remained little understood. Here, we use single nucleotide polymorphism (SNP) data from multiple coastal species to track the dynamics of cocolonization events in response to ancient earthquake disturbance in southern New Zealand. Specifically, we use a comparative phylogeographic approach to understand the extent to which epifauna (with varying ecological associations with their macroalgal hosts) share comparable spatial and temporal recolonization patterns. Our study reveals concordant disturbance-related phylogeographic breaks in two intertidal macroalgal species along with two associated epibiotic species (a chiton and an isopod). By contrast, two codistributed species, one of which is an epibiotic amphipod and the other a subtidal macroalga, show few, if any, genetic effects of palaeoseismic coastal uplift. Phylogeographic model selection reveals similar post-uplift recolonization routes for the epibiotic chiton and isopod and their macroalgal hosts. Additionally, codemographic analyses support synchronous population expansions of these four phylogeographically similar taxa. Our findings indicate that coastal paleoseismic activity has driven concordant impacts on multiple codistributed species, with concerted recolonization events probably facilitated by macroalgal rafting. These results highlight that high-resolution comparative genomic data can help reconstruct concerted multispecies responses to recent ecological disturbance.

The evolution of comparative phylogeography: putting the geography (and more) into comparative population genomics

Comparative population genomics is an ascendant field using genomic comparisons between species to draw inferences about forces regulating genetic variation. Comparative phylogeography, by contrast, focuses on the shared lineage histories of species codistributed geographically and is decidedly organismal in perspective. Comparative phylogeography is approximately 35 years old, and, by some metrics, is showing signs of reduced growth. Here, we contrast the goals and methods of comparative population genomics and comparative phylogeography and argue that comparative phylogeography offers an important perspective on evolutionary history that succeeds in integrating genomics with landscape evolution in ways that complement the suprageographic perspective of comparative population genomics. Focusing primarily on terrestrial vertebrates, we review the history of comparative phylogeography, its milestones and ongoing conceptual innovations, its increasingly global focus, and its status as a bridge between landscape genomics and the process of speciation. We also argue that, as a science with a strong “sense of place,” comparative phylogeography offers abundant “place-based” educational opportunities with its focus on geography and natural history, as well as opportunities for collaboration with local communities and indigenous peoples. Although comparative phylogeography does not yet require whole-genome sequencing for many of its goals, we conclude that it nonetheless plays an important role in grounding our interpretation of genetic variation in the fundamentals of geography and Earth history.

Discordant Pleistocene population size histories in a guild of hymenopteran parasitoids

Signatures of past changes in population size have been detected in genome-wide variation in many species. However, the causes of such demographic changes and the extent to which they are shared across co-distributed species remain poorly understood. During Pleistocene glacial maxima, many temperate European species were confined to southern refugia. While vicariance and range expansion processes associated with glacial cycles have been widely documented, it is unclear whether refugial populations of co-distributed species have experienced shared histories of population size change. We analyse whole-genome sequence data to reconstruct and compare demographic histories during the Quaternary for Iberian refuge populations in a single ecological guild (seven species of chalcid parasitoid wasps associated with oak cynipid galls). For four of these species, we find support for large changes in effective population size (N_e ) through the Pleistocene that coincide with major climate events. However, there is little evidence that the timing, direction and magnitude of demographic change are shared across species, suggesting that demographic histories in this guild are largely idiosyncratic, even at the scale of a single glacial refugium.

The Pleistocene species pump past its prime: Evidence from European butterfly sister species

The Pleistocene glacial cycles had a profound impact on the ranges and genetic make-up of organisms. While it is clear that the contact zones that have been described for many sister taxa are secondary and have formed in the current interglacial, it is unclear when the taxa involved began to diverge. Previous estimates based on small numbers of loci are unreliable given the stochasticity of genetic drift and the contrasting effects of incomplete lineage sorting and gene flow on gene divergence. Here, we use genome-wide transcriptome data to estimate divergence for 18 sister species pairs of European butterflies showing either sympatric or contact zone distributions. We find that in most cases, species divergence predates the mid-Pleistocene transition or even the entire Pleistocene period. We also show that although post-divergence gene flow is restricted to contact zone pairs, they are not systematically younger than sympatric pairs. This suggests that contact zones are not limited to the initial stages of the speciation process, but can involve notably old taxa. Finally, we show that mitochondrial divergence and nuclear divergence are only weakly correlated and mitochondrial divergence is higher for contact zone pairs.

Moving apart together: co-movement of a symbiont community and their ant host, and its importance for community assembly

BACKGROUND

Species interactions may affect spatial dynamics when the movement of one species is determined by the presence of another one. The most direct species-dependence of dispersal is vectored, usually cross-kingdom, movement of immobile parasites, diseases or seeds by mobile animals. Joint movements of species should, however, not be vectored by definition, as even mobile species are predicted to move together when they are tightly connected in symbiont communities.

METHODS

We studied concerted movements in a diverse and heterogeneous community of arthropods (myrmecophiles) associated with red wood ants. We questioned whether joint-movement strategies eventually determine and speed-up community succession.

RESULTS

We recorded an astonishingly high number of obligate myrmecophiles outside red wood ant nests. They preferentially co-moved with the host ants as the highest densities were found in locations with the highest density of foraging red wood ants, such as along the network of ant trails. These observations suggest that myrmecophiles resort to the host to move away from the nest, and this to a much higher extent than hitherto anticipated. Interestingly, functional groups of symbionts displayed different dispersal kernels, with predatory myrmecophiles moving more frequently and further from the nest than detritivorous myrmecophiles. We discovered that myrmecophile diversity was lower in newly founded nests than in mature red wood ant nests. Most myrmecophiles, however, were able to colonize new nests fast suggesting that the heterogeneity in mobility does not affect community assembly.

CONCLUSIONS

We show that co-movement is not restricted to tight parasitic, or cross-kingdom interactions. Movement in social insect symbiont communities may be heterogeneous and functional group-dependent, but clearly affected by host movement. Ultimately, this co-movement leads to directional movement and allows a fast colonisation of new patches, but not in a predictable way. This study highlights the importance of spatial dynamics of local and regional networks in symbiont metacommunities, of which those of symbionts of social insects are prime examples.

Why are there not more herbivorous insect species?

Insect species richness is estimated to exceed three million species, of which roughly half is herbivorous. Despite the vast number of species and varied life histories, the proportion of herbivorous species among plant-consuming organisms is lower than it could be due to constraints that impose limits to their diversification. These include ecological factors, such as vague interspecific competition; anatomical and physiological limits, such as neural limits and inability of handling a wide range of plant allelochemicals; phylogenetic constraints, like niche conservatism; and most importantly, a low level of concerted genetic variation necessary to a phyletic conversion. It is suggested that diversification ultimately depends on what we call the intrinsic trend of diversification of the insect genome. In support of the above, we survey the major types of host-specificity, the mechanisms and constraints of host specialization, possible pathways of speciation, and hypotheses concerning insect diversification.

Low host specificity and broad geographical ranges in a community of parasitic non-pollinating fig wasps (Sycoryctinae; Chalcidoidea)

Plants, phytophagous insects and their parasitoids form the most diverse assemblages of macroscopic organisms on earth. Enclosed assemblages in particular represent a tractable system for studying community assembly and diversification. Communities associated with widespread plant species are especially suitable as they facilitate a comparative approach. Pantropical fig-wasp communities represent a remarkably well-replicated system, ideal for studying these historical processes. We expect high dispersal ability in non-pollinating fig wasps to result in lower geographical turnover in comparison to pollinating fig wasps. The ability of non-pollinating wasps to utilise a number of hosts (low host specificity) is a key determinant of overall geographical range, with intraspecific competition becoming a constraining factor should diet breadth overlap among species. Finally, we expect conserved community structure throughout the host range. We aim to test these expectations, derived from population genetic and community studies, using the multi-trophic insect community associated with Ficus hirta throughout its 3,500 km range across continental and insular Asia. We collect molecular evidence from one coding mitochondrial gene, one non-coding nuclear gene and multiple microsatellites across 25 geographical sites. Using these data, we establish species boundaries, determine levels of host specificity among non-pollinating fig wasps and quantify geographical variation in community composition. We find low host specificity in two genera of non-pollinating fig wasps. Functional community structure is largely conserved across the range of the host fig, despite limited correspondence between the ranges of non-pollinator and pollinator species. While nine pollinators are associated with Ficus hirta, the two non-pollinator tribes developing in its figs each contained only four species. Contrary to predictions, we find stronger isolation by distance in non-pollinators than pollinators. Long-lived non-pollinators may disperse more gradually and be less reliant on infrequent long-distance dispersal by wind currents. Segregation among non-pollinating species across their range is suggestive of competitive exclusion and we propose that this may be a result of increased levels of local adaptation and moderate, but regular, rates of dispersal. Our findings provide one more example of lack of strict codiversification in the geographical diversification of plant-associated insect communities.

DILS: Demographic inferences with linked selection by using ABC

We present DILS, a deployable statistical analysis platform for conducting demographic inferences with linked selection from population genomic data using an Approximate Bayesian Computation framework. DILS takes as input single-population or two-population data sets (multilocus fasta sequences) and performs three types of analyses in a hierarchical manner, identifying: (a) the best demographic model to study the importance of gene flow and population size change on the genetic patterns of polymorphism and divergence, (b) the best genomic model to determine whether the effective size Ne and migration rate N, m are heterogeneously distributed along the genome (implying linked selection) and (c) loci in genomic regions most associated with barriers to gene flow. Also available via a Web interface, an objective of DILS is to facilitate collaborative research in speciation genomics. Here, we show the performance and limitations of DILS by using simulations and finally apply the method to published data on a divergence continuum composed by 28 pairs of Mytilus mussel populations/species.

From Inquilines to Gall Inducers: Genomic Signature of a Life-Style Transition in Synergus Gall Wasps

Gall wasps (Hymenoptera: Cynipidae) induce complex galls on oaks, roses, and other plants, but the mechanism of gall induction is still unknown. Here, we take a comparative genomic approach to revealing the genetic basis of gall induction. We focus on Synergus itoensis, a species that induces galls inside oak acorns. Previous studies suggested that this species evolved the ability to initiate gall formation recently, as it is deeply nested within the genus Synergus, whose members are mostly inquilines that develop inside the galls of other species. We compared the genome of S. itoensis with that of three related Synergus inquilines to identify genomic changes associated with the origin of gall induction. We used a novel Bayesian selection analysis, which accounts for branch-specific and gene-specific selection effects, to search for signatures of selection in 7,600 single-copy orthologous genes shared by the four Synergus species. We found that the terminal branch leading to S. itoensis had more genes with a significantly elevated dN/dS ratio (positive signature genes) than the other terminal branches in the tree; the S. itoensis branch also had more genes with a significantly decreased dN/dS ratio. Gene set enrichment analysis showed that the positive signature gene set of S. itoensis, unlike those of the inquiline species, is enriched in several biological process Gene Ontology terms, the most prominent of which is “Ovarian Follicle Cell Development.” Our results indicate that the origin of gall induction is associated with distinct genomic changes, and provide a good starting point for further characterization of the genes involved.

Molecular markers reveal diversity in composition of Megastigmus (Hymenoptera: Megastigmidae) from eucalypt galls

Since outbreaks of the invasive blue gum chalcids Leptocybe spp. began, the genus Megastigmus (Hymenoptera: Megastigmidae) has been increasingly studied as containing potential biocontrol agents against these pests. Megastigmus species have been collected and described from Australia, the presumed origin of Leptocybe spp., with M. zvimendeli and M. lawsoni reported as Leptocybe spp. parasitoids established outside of Australia. Parasitic Megastigmus have been reported to occur locally in the Neotropics, Afrotropic, Palearctic, and Indomalaya biogeographic realms, and in many cases described as new to science. However, molecular tools have not been used in studying parasitic Megastigmus, and difficulties in morphological taxonomy have compromised further understanding of eucalypt-associated Megastigmus as well as the Megastigmus-Leptocybe association. In this study, we used molecular markers to study the species composition and phylogeny of Megastigmus collected from eucalypt galls in Australia and from Leptocybe spp. galls from South Africa, Kenya, Israel, China, and Vietnam. We record thirteen discrete species and a species complex associated with eucalypt galls. A summary of morphological characters is provided to assist morphological delimitation of the studied group. A phylogeny based on 28S rDNA identified species groups of importance to Leptocybe spp. biocontrol agents from four clades with nine species. Relationships between Megastigmus from eucalypt galls and their phytophagous congeners were unresolved. Further molecular work is needed to clarify the identity of many species.

Dispersal barriers and opportunities drive multiple levels of phylogeographic concordance in the Southern Alps of New Zealand

Phylogeographic concordance, or the sharing of phylogeographic patterns among codistributed species, suggests similar responses to topography or climatic history. While the orientation and timing of breaks between lineages are routinely compared, spatial dynamics within regions occupied by individual lineages provide a second opportunity for comparing responses to past events. In environments with complex topography and glacial history, such as New Zealand's South Island, geographically nested comparisons can identify the processes leading to phylogeographic concordance between and within regional genomic clusters. Here, we used single nucleotide polymorphisms (obtained via ddRADseq) for two codistributed forest beetle species, Agyrtodes labralis (Leiodidae) and Brachynopus scutellaris (Staphylinidae), to evaluate the role of climate change and topography in shaping phylogeographic concordance at two, nested spatial scales: do species diverge over the same geographic barriers, with similar divergence times? And within regions delimited by these breaks, do species share similar spatial dynamics of directional expansion or isolation-by-distance? We found greater congruence of phylogeographic breaks between regions divided by the strongest dispersal barriers (i.e., the Southern Alps). However, these shared breaks were not indicative of shared spatial dynamics within the regions they delimit, and the most similar spatial dynamics between species occurred within regions with the strongest gradients in historical climatic stability. Our results indicate that lack of concordance as traditionally detected by lineage turnover does not rule out the possibility of shared histories, and variation in the presence and type of concordance may provide insights into the different processes shaping phylogeographic patterns across geologically dynamic regions.

Low-coverage genomic data resolve the population divergence and gene flow history of an Australian rain forest fig wasp

Population divergence and gene flow are key processes in evolution and ecology. Model-based analysis of genome-wide data sets allows discrimination between alternative scenarios for these processes even in nonmodel taxa. We used two complementary approaches (one based on the blockwise site frequency spectrum [bSFS], the second on the pairwise sequentially Markovian coalescent [PSMC]) to infer the divergence history of a fig wasp, Pleistodontes nigriventris. Pleistodontes nigriventris and its fig tree mutualist Ficus watkinsiana are restricted to rain forest patches along the eastern coast of Australia and are separated into The Northern population is to the north of the Southern populations by two dry forest corridors (the Burdekin and St. Lawrence Gaps). We generated whole genome sequence data for two haploid males per population and used the bSFS approach to infer the timing of divergence between northern and southern populations of P. nigriventris, and to discriminate between alternative isolation with migration (IM) and instantaneous admixture (ADM) models of postdivergence gene flow. Pleistodontes nigriventris has low genetic diversity (π = 0.0008), to our knowledge one of the lowest estimates reported for a sexually reproducing arthropod. We find strongest support for an ADM model in which the two populations diverged ca. 196 kya in the late Pleistocene, with almost 25% of northern lineages introduced from the south during an admixture event ca. 57 kya. This divergence history is highly concordant with individual population demographies inferred from each pair of haploid males using PSMC. Our analysis illustrates the inferences possible with genome-level data for small population samples of tiny, nonmodel organisms and adds to a growing body of knowledge on the population structure of Australian rain forest taxa.

The Role of Evolution in Shaping Ecological Networks

The structure of ecological networks reflects the evolutionary history of their biotic components, and their dynamics are strongly driven by ecoevolutionary processes. Here, we present an appraisal of recent relevant research, in which the pervasive role of evolution within ecological networks is manifest. Although evolutionary processes are most evident at macroevolutionary scales, they are also important drivers of local network structure and dynamics. We propose components of a blueprint for further research, emphasising process-based models, experimental evolution, and phenotypic variation, across a range of distinct spatial and temporal scales. Evolutionary dimensions are required to advance our understanding of foundational properties of community assembly and to enhance our capability of predicting how networks will respond to impending changes.

Genomic dissection of an extended phenotype: Oak galling by a cynipid gall wasp

Galls are plant tissues whose development is induced by another organism for the inducer's benefit. 30,000 arthropod species induce galls, and in most cases the inducing effectors and target plant systems are unknown. Cynipid gall wasps are a speciose monophyletic radiation that induce structurally complex galls on oaks and other plants. We used a model system comprising the gall wasp Biorhiza pallida and the oak Quercus robur to characterise inducer and host plant gene expression at defined stages through the development of galled and ungalled plant tissues, and tested alternative hypotheses for the origin and type of galling effectors and plant metabolic pathways involved. Oak gene expression patterns diverged markedly during development of galled and normal buds. Young galls showed elevated expression of oak genes similar to legume root nodule Nod factor-induced early nodulin (ENOD) genes and developmental parallels with oak buds. In contrast, mature galls showed substantially different patterns of gene expression to mature leaves. While most oak transcripts could be functionally annotated, many gall wasp transcripts of interest were novel. We found no evidence in the gall wasp for involvement of third-party symbionts in gall induction, for effector delivery using virus-like-particles, or for gallwasp expression of genes coding for plant hormones. Many differentially and highly expressed genes in young larvae encoded secretory peptides, which we hypothesise are effector proteins exported to plant tissues. Specifically, we propose that host arabinogalactan proteins and gall wasp chitinases interact in young galls to generate a somatic embryogenesis-like process in oak tissues surrounding the gall wasp larvae. Gall wasp larvae also expressed genes encoding multiple plant cell wall degrading enzymes (PCWDEs). These have functional orthologues in other gall inducing cynipids but not in figitid parasitoid sister groups, suggesting that they may be evolutionary innovations associated with cynipid gall induction.

Plant galls are induced by organisms that manipulate host plant development to produce novel structures. The organisms involved range from mutualistic (such as nitrogen fixing bacteria) to parasitic. In the case of parasites, the gall benefits only the gall-inducing partner. A wide range of organisms can induce galls, but the processes involved are understood only for some bacterial and fungal galls. Cynipid gall wasps induce diverse and structurally complex galls, particularly on oaks (Quercus). We used transcriptome and genome sequencing for one gall wasp and its host oak to identify genes active in gall development. On the plant side, when compared to normally developing bud tissues, young gall tissues showed elevated expression of loci similar to those found in nitrogen-fixing root nodules of leguminous plants. On the wasp side, we found no evidence for involvement of viruses or microorganisms carried by the insects in gall induction or delivery of inducing stimuli. We found that gall wasps express many genes whose products may be secreted to the host, including enzymes that degrade plant cell walls. Genome comparisons between galling and non-galling relatives showed cell wall-degrading enzymes are restricted to gall inducers, and hence potentially key components of a gall inducing lifestyle.

The determinants of genetic diversity in butterflies

Under the neutral theory, genetic diversity is expected to increase with population size. While comparative analyses have consistently failed to find strong relationships between census population size and genetic diversity, a recent study across animals identified a strong correlation between propagule size and genetic diversity, suggesting that r-strategists that produce many small offspring, have greater long-term population sizes. Here we compare genome-wide genetic diversity across 38 species of European butterflies (Papilionoidea), a group that shows little variation in reproductive strategy. We show that genetic diversity across butterflies varies over an order of magnitude and that this variation cannot be explained by differences in current abundance, propagule size, host or geographic range. Instead, neutral genetic diversity is negatively correlated with body size and positively with the length of the genetic map. This suggests that genetic diversity is determined both by differences in long-term population size and the effect of selection on linked sites.

ABLE: blockwise site frequency spectra for inferring complex population histories and recombination

We introduce ABLE (Approximate Blockwise Likelihood Estimation), a novel simulation-based composite likelihood method that uses the blockwise site frequency spectrum to jointly infer past demography and recombination. ABLE is explicitly designed for a wide variety of data from unphased diploid genomes to genome-wide multi-locus data (for example, RADSeq) and can also accommodate arbitrarily large samples. We use simulations to demonstrate the accuracy of this method to infer complex histories of divergence and gene flow and reanalyze whole genome data from two species of orangutan. ABLE is available for download at https://github.com/champost/ABLE.