Víctor Soria-Carrasco

Stick insect genomes reveal natural selection's role in parallel speciation.

Stick to the Bush Can the underlying genetic changes driving the divergence of populations into new species be predicted or repeated? Soria-Carrasco et al. (p. 738) investigated the genetic changes observed after one generation when stick insect (Timema cristinae) populations were transplanted from their preferred host plants to alternative hosts. Diverged genetic regions were relatively small, with most loci showing divergence in a single population pair. However, the number of loci showing parallel divergence was greater than expected by chance. Thus, selection can drive parallel phenotypic evolution via parallel genetic changes. Parallel speciation in insects shows both convergent and divergent selection after one generation. Natural selection can drive the repeated evolution of reproductive isolation, but the genomic basis of parallel speciation remains poorly understood. We analyzed whole-genome divergence between replicate pairs of stick insect populations that are adapted to different host plants and undergoing parallel speciation. We found thousands of modest-sized genomic regions of accentuated divergence between populations, most of which are unique to individual population pairs. We also detected parallel genomic divergence across population pairs involving an excess of coding genes with specific molecular functions. Regions of parallel genomic divergence in nature exhibited exceptional allele frequency changes between hosts in a field transplant experiment. The results advance understanding of biological diversification by providing convergent observational and experimental evidence for selection’s role in driving repeatable genomic divergence.

The K tree score

SUMMARY We introduce a new phylogenetic comparison method that measures overall differences in the relative branch length and topology of two phylogenetic trees. To do this, the algorithm first scales one of the trees to have a global divergence as similar as possible to the other tree. Then, the branch length distance, which takes differences in topology and branch lengths into account, is applied to the two trees. We thus obtain the minimum branch length distance or K tree score. Two trees with very different relative branch lengths get a high K score whereas two trees that follow a similar among-lineage rate variation get a low score, regardless of the overall rates in both trees. There are several applications of the K tree score, two of which are explained here in more detail. First, this score allows the evaluation of the performance of phylogenetic algorithms, not only with respect to their topological accuracy, but also with respect to the reproduction of a given branch length variation. In a second example, we show how the K score allows the selection of orthologous genes by choosing those that better follow the overall shape of a given reference tree. AVAILABILITY http://molevol.ibmb.csic.es/Ktreedist.html

Transitions between phases of genomic differentiation during stick-insect speciation

Speciation can involve a transition from a few genetic loci that are resistant to gene flow to genome-wide differentiation. However, only limited data exist concerning this transition and the factors promoting it. Here, we study phases of speciation using data from >100 populations of 11 species of Timema stick insects. Consistent with early phases of genic speciation, adaptive colour-pattern loci reside in localized genetic regions of accentuated differentiation between populations experiencing gene flow. Transitions to genome-wide differentiation are also observed with gene flow, in association with differentiation in polygenic chemical traits affecting mate choice. Thus, intermediate phases of speciation are associated with genome-wide differentiation and mate choice, but not growth of a few genomic islands. We also find a gap in genomic differentiation between sympatric taxa that still exchange genes and those that do not, highlighting the association between differentiation and complete reproductive isolation. Our results suggest that substantial progress towards speciation may involve the alignment of multi-faceted aspects of differentiation.

The evolution of novel host use is unlikely to be constrained by trade‐offs or a lack of genetic variation

The genetic and ecological factors that shape the evolution of animal diets remain poorly understood. For herbivorous insects, the expectation has been that trade‐offs exist, such that adaptation to one host plant reduces performance on other potential hosts. We investigated the genetic architecture of alternative host use by rearing individual Lycaeides melissa butterflies from two wild populations in a crossed design on two hosts (one native and one introduced) and analysing the genetic basis of differences in performance using genomic approaches. Survival during the experiment was highest when butterfly larvae were reared on their natal host plant, consistent with local adaptation. However, cross‐host correlations in performance among families (within populations) were not different from zero. We found that L. melissa populations possess genetic variation for larval performance and variation in performance had a polygenic basis. We documented very few genetic variants with trade‐offs that would inherently constrain diet breadth by preventing the optimization of performance across hosts. Instead, most genetic variants that affected performance on one host had little to no effect on the other host. In total, these results suggest that genetic trade‐offs are not the primary cause of dietary specialization in L. melissa butterflies.

Diversification rates and the latitudinal gradient of diversity in mammals

The latitudinal gradient of species richness has frequently been attributed to higher diversification rates of tropical groups. In order to test this hypothesis for mammals, we used a set of 232 genera taken from a mammalian supertree and, additionally, we reconstructed dated Bayesian phylogenetic trees of 100 genera. For each genus, diversification rate was estimated taking incomplete species sampling into account and latitude was assigned considering the heterogeneity in species distribution ranges. For both datasets, we found that the average diversification rate was similar among all latitudinal bands. Furthermore, when we used phylogenetically independent contrasts, we did not find any significant correlation between latitude and diversification parameters, including different estimates of speciation and extinction rates. Thus, other factors, such as the dynamics of dispersal through time, may be required to explain the latitudinal gradient of diversity in mammals.

Genome-Wide Association Mapping of Phenotypic Traits Subject to a Range of Intensities of Natural Selection in Timema cristinae*

The genetic architecture of adaptive traits can reflect the evolutionary history of populations and also shape divergence among populations. Despite this central role in evolution, relatively little is known regarding the genetic architecture of adaptive traits in nature, particularly for traits subject to known selection intensities. Here we quantitatively describe the genetic architecture of traits that are subject to known intensities of differential selection between host plant species in Timema cristinae stick insects. Specifically, we used phenotypic measurements of 10 traits and 211,004 single-nucleotide polymorphisms (SNPs) to conduct multilocus genome-wide association mapping. We identified a modest number of SNPs that were associated with traits and sometimes explained a large proportion of trait variation. These SNPs varied in their strength of association with traits, and both major and minor effect loci were discovered. However, we found no relationship between variation in levels of divergence among traits in nature and variation in parameters describing the genetic architecture of those same traits. Our results provide a first step toward identifying loci underlying adaptation in T. cristinae. Future studies will examine the genomic location, population differentiation, and response to selection of the trait-associated SNPs described here.

Estimation of phylogenetic inconsistencies in the three domains of life.

Discrepancies in phylogenetic trees of bacteria and archaea are often explained as lateral gene transfer events. However, such discrepancies may also be due to phylogenetic artifacts or orthology assignment problems. A first step that may help to resolve this dilemma is to estimate the extent of phylogenetic inconsistencies in trees of prokaryotes in comparison with those of higher eukaryotes, where no lateral gene transfer is expected. To test this, we used 21 proteomes each of eukaryotes (mainly opisthokonts), proteobacteria, and archaea that spanned equivalent levels of genetic divergence. In each domain of life, we defined a set of putative orthologous sequences using a phylogenetic-based orthology protocol and, as a reference topology, we used a tree constructed with concatenated genes of each domain. Our results show, for most of the tests performed, that the magnitude of topological inconsistencies with respect to the reference tree was very similar in the trees of proteobacteria and eukaryotes. When clade support was taken into account, prokaryotes showed some more inconsistencies, but then all values were very low. Discrepancies were only consistently higher in archaea but, as shown by simulation analysis, this is likely due to the particular tree of the archaeal species used here being more difficult to reconstruct, whereas the trees of proteobacteria and eukaryotes were of similar difficulty. Although these results are based on a relatively small number of genes, it seems that phylogenetic reconstruction problems, including orthology assignment problems, have a similar overall effect over prokaryotic and eukaryotic trees based on single genes. Consequently, lateral gene transfer between distant prokaryotic species may have been more rare than previously thought, which opens the way to obtain the tree of life of bacterial and archaeal species using genomic data and the concatenation of adequate genes, in the same way as it is usually done in eukaryotes.

Congruence test of molecular clock calibration hypotheses based on Bayes factor comparisons

This work was supported by the Spanish Ministry of Science and Innovation (CGL2006/06706 and CGL2009-10906 to C.A. and J.S. and CGL2008-00007 to J.G.-Z., the latter also with support of the European Regional Development Fund). CA received support from an FPU predoctoral studentship (Spanish Ministry of Education).

Natural selection and the predictability of evolution in Timema stick insects

Estimating the predictability of evolution Evolution results from expected effects, such as selection driving alleles toward fixation, and stochastic effects, such as unusual environmental variation and genetic drift. To determine the potential to predict evolutionary change, Nosil et al. examined three naturally occurring morphs of stick insects (see the Perspective by Reznick and Travis). They wanted to determine which selective parameters could be used to foresee changes, despite varying environmental conditions. One morph fit a model of negative frequency-dependent selection, likely owing to predation, but changes in other morph frequencies remained unpredictable. Thus, for specific cases, we can forecast short-term changes within populations, but evolution is more difficult to predict when it involves a balance between multiple selective factors and uncertainty in environmental conditions. Science, this issue p. 765; see also p. 738 Frequencies of stick insect coloration and patterning over time are used to determine the predictability of evolution. Predicting evolution remains difficult. We studied the evolution of cryptic body coloration and pattern in a stick insect using 25 years of field data, experiments, and genomics. We found that evolution is more difficult to predict when it involves a balance between multiple selective factors and uncertainty in environmental conditions than when it involves feedback loops that cause consistent back-and-forth fluctuations. Specifically, changes in color-morph frequencies are modestly predictable through time (r2 = 0.14) and driven by complex selective regimes and yearly fluctuations in climate. In contrast, temporal changes in pattern-morph frequencies are highly predictable due to negative frequency-dependent selection (r2 = 0.86). For both traits, however, natural selection drives evolution around a dynamic equilibrium, providing some predictability to the process.

Color phenotypes are under similar genetic control in two distantly related species of Timema stick insect

Ecology and genetics are both of general interest to evolutionary biologists as they can influence the phenotypic and genetic response to selection. The stick insects Timema podura and Timema cristinae exhibit a green/melanistic body color polymorphism that is subject to different ecologically based selective regimes in the two species. Here, we describe aspects of the genetics of this color polymorphism in T. podura, and compare this to previous results in T. cristinae. We first show that similar color phenotypes of the two species cluster in phenotypic space. We then use genome‐wide association mapping to show that in both species, color is controlled by few loci, dominance relationships between color alleles are the same, and SNPs associated with color phenotypes colocalize to the same linkage group. Regions within this linkage group that harbor genetic variants associated with color exhibit elevated linkage disequilibrium relative to genome wide expectations, but more strongly so in T. cristinae. We use these results to discuss predictions regarding how the genetics of color could influence levels of phenotypic and genetic variation that segregate within and between populations of T. podura and T. cristinae, drawing parallels with other organisms.