Demography and genome divergence of lake and stream populations of an East African cichlid fish

Disentangling the processes and mechanisms underlying adaptive diversification is facilitated by the comparative study of replicate population pairs that have diverged along a similar environmental gradient. Such a setting is realized in a cichlid fish from southern Lake Tanganyika, Astatotilapia burtoni, which occurs within the lake proper as well as in various affluent rivers. Previously, we demonstrated that independent lake and stream populations show similar adaptations to the two habitat regimes. However, little is known about the evolutionary and demographic history of the A. burtoni populations in question and the patterns of genome divergence among them. Here, we apply restriction site‐associated DNA sequencing (RADseq) to examine the evolutionary history, the population structure and genomic differentiation of lake and stream populations in A. burtoni. A phylogenetic reconstruction based on genome‐wide molecular data largely resolved the evolutionary relationships among populations, allowing us to re‐evaluate the independence of replicate lake–stream population clusters. Further, we detected a strong pattern of isolation by distance, with baseline genomic divergence increasing with geographic distance and decreasing with the level of gene flow between lake and stream populations. Genome divergence patterns were heterogeneous and inconsistent among lake‐stream population clusters, which is explained by differences in divergence times, levels of gene flow and local selection regimes. In line with the latter, we only detected consistent outlier loci when the most divergent lake–stream population pair was excluded. Several of the thus identified candidate genes have inferred functions in immune and neuronal systems and show differences in gene expression between lake and stream populations.     

Change in sexual signalling traits outruns morphological divergence across an ecological gradient in the post‐glacial radiation of the songbird genus Junco

The relative roles of natural and sexual selection in promoting evolutionary lineage divergence remains controversial and difficult to assess in natural systems. Local adaptation through natural selection is known to play a central role in promoting evolutionary divergence, yet secondary sexual traits can vary widely among species in recent radiations, suggesting that sexual selection may also be important in the early stages of speciation. Here, we compare rates of divergence in ecologically relevant traits (morphology) and sexually selected signalling traits (coloration) relative to neutral structure in genome‐wide molecular markers and examine patterns of variation in sexual dichromatism to explore the roles of natural and sexual selection in the diversification of the songbird genus Junco (Aves: Passerellidae). Juncos include divergent lineages in Central America and several dark‐eyed junco (J. hyemalis) lineages that diversified recently as the group recolonized North America following the last glacial maximum (ca. 18,000 years ago). We found an accelerated rate of divergence in sexually selected characters relative to ecologically relevant traits. Moreover, sexual dichromatism measurements suggested a positive relationship between the degree of colour divergence and the strength of sexual selection when controlling for neutral genetic distance. We also found a positive correlation between dichromatism and latitude, which coincides with the geographic axis of decreasing lineage age in juncos but also with a steep ecological gradient. Finally, we found significant associations between genome‐wide variants linked to functional genes and proxies of both sexual and natural selection. These results suggest that the joint effects of sexual and ecological selection have played a prominent role in the junco radiation.     

Inbreeding and outbreeding depression in Stylidium hispidum: implications for mixing seed sources for ecological restoration

The benefits of composite rather than local seed provenances for ecological restoration have recently been argued, largely on the basis of maximizing evolutionary potential. However, these arguments have downplayed the potentially negative consequences of outbreeding depression once mixed provenances interbreed. In this study, we compared intraspecific F1 hybrid performance and molecular marker differentiation among four populations of Stylidium hispidum, a species endemic to Southwestern Australia. Multivariate ordination of 134 AFLP markers analyzed genetic structure and detected two clusters of paired sites that diverged significantly for marker variation along a latitudinal boundary. To test for outbreeding depression and to determine the consequences of molecular population divergence for hybrid fitness, we conducted controlled pollinations and studied germination and survival for three cross categories (within‐population crosses, short‐ and long‐distance F1 hybrids) for paired sites distributed within and between the two genetically differentiated regions. We found evidence of outbreeding depression in long‐distance hybrids (111–124 km), and inbreeding depression among progeny of within‐population crosses, relative to short‐distance (3–10 km) hybrids, suggesting an intermediate optimal outcrossing distance in this species. These results are discussed in light of the evolutionary consequences of mixing seed sources for biodiversity restoration.     

Divergent selection and drift shape the genomes of two avian sister species spanning a saline–freshwater ecotone

The role of species divergence due to ecologically based divergent selection—or ecological speciation—in generating and maintaining biodiversity is a central question in evolutionary biology. Comparison of the genomes of phylogenetically related taxa spanning a selective habitat gradient enables discovery of divergent signatures of selection and thereby provides valuable insight into the role of divergent ecological selection in speciation. Tidal marsh ecosystems provide tractable opportunities for studying organisms' adaptations to selective pressures that underlie ecological divergence. Sharp environmental gradients across the saline–freshwater ecotone within tidal marshes present extreme adaptive challenges to terrestrial vertebrates. Here, we sequence 20 whole genomes of two avian sister species endemic to tidal marshes—the saltmarsh sparrow (Ammospiza caudacutus) and Nelson's sparrow (A. nelsoni)—to evaluate the influence of selective and demographic processes in shaping genome‐wide patterns of divergence. Genome‐wide divergence between these two recently diverged sister species was notably high (genome‐wide F_ST = 0.32). Against a background of high genome‐wide divergence, regions of elevated divergence were widespread throughout the genome, as opposed to focused within islands of differentiation. These patterns may be the result of genetic drift resulting from past tidal march colonization events in conjunction with divergent selection to different environments. We identified several candidate genes that exhibited elevated divergence between saltmarsh and Nelson's sparrows, including genes linked to osmotic regulation, circadian rhythm, and plumage melanism—all putative candidates linked to adaptation to tidal marsh environments. These findings provide new insights into the roles of divergent selection and genetic drift in generating and maintaining biodiversity.     

Ecological Predictors of Range Areas and Use of Burrow Systems in the Diurnal Rodent, Octodon degus

Variation in animal space use patterns may be linked to numerous ecological factors affecting survival and reproduction. We examined the relationship between ecology and above‐ and below‐ground components of space use by Octodon degus, a semi‐fossorial rodent in Chile. We monitored the daytime minimum convex polygon and adaptive kernel range areas of 26 individuals and determined the number of burrow systems used by degus during night‐time radiotelemetry and trapping of burrow systems on two study grids at Rinconada de Maipú, a semi‐arid Matorral in central Chile. We quantified food biomass, soil hardness, distance to overhead vegetative cover, and density of burrow openings at putative nest burrows. Degus living on the grid with more shrub cover had larger range areas than degus living on the grid with less cover. The range areas of degus decreased with increasing distance from overhead vegetative cover. There was a weak (but statistically significant) negative relationship between the number of burrow systems used by degus and the distance to vegetative cover and density of burrow openings at burrow systems. Male and female degus had similar range areas. Our results suggest that overhead cover decreases the risk of predation to male and female degus. Degus probably balance the benefits of numerous burrow openings (reduced predation risk) with time and energy requirements of burrow construction and maintenance. Models of space use that consider the effect of multiple ecological variables should measure different dimensions of space use.     

Cyto‐nuclear discordance suggests complex evolutionary history in the cave‐dwelling salamander,Eurycea lucifuga

Our understanding of the evolutionary history and ecology of cave‐associated species has been driven historically by studies of morphologically adapted cave‐restricted species. Our understanding of the evolutionary history and ecology of nonrestricted cave species, troglophiles, is limited to a few studies, which present differing accounts of troglophiles’ relationship with the cave habitat, and its impact on population dynamics. Here, we used phylogenetics, demographic statistics, and population genetic methods to study lineage divergence, dates of divergence, and population structure in the Cave Salamander, Eurycea lucifuga, across its range. In order to perform these analyses, we sampled 233 individuals from 49 populations, using sequence data from three gene loci as well as genotyping data from 19 newly designed microsatellite markers. We find, as in many other species studied in a phylogeographic context, discordance between patterns inferred from mitochondrial relationships and those inferred by nuclear markers indicating a complicated evolutionary history in this species. Our results suggest Pleistocene‐based divergence among three main lineages within E. lucifuga corresponding to the western, central, and eastern regions of the range, similar to patterns seen in species separated in multiple refugia during climatic shifts. The conflict between mitochondrial and nuclear patterns is consistent with what we would expect from secondary contact between regional populations following expansion from multiple refugia.     

The stochastic model and deviations from randomness in eukaryotictRNAs: Comparison with the PAM approach

Summary The stochastic model of molecular evolution was used to makea priori predictions for the total number of one-step nucleotide changes required to account for a given observed number of nucleotide substitutions between two homologous nucleic acids. The experimental deviations from randomness found for eukaryotic transfer RNAs and summarized by Dayhoff and McLaughlin (1972) are shown to affect only slightly the quantitative predictions of the model. This is true for both short and long periods of evolutionary divergence. The model can thus be used with some confidence for quantitatively correcting the branch lengths of phylogenetic trees derived from either nucleic acid sequence or hybridization data.     

Mitochondrial genes and mammalian phylogenies: increasing the reliability of branch length estimation

Since branch lengths provide important information about the timing and the extent of evolutionary divergence among taxa, accurate resolution of evolutionary history depends as much on branch length estimates as on recovery of the correct topology. However, the empirical relationship between the choice of genes to sequence and the quality of branch length estimation remains ill defined. To address this issue, we evaluated the accuracy of branch lengths estimated from subsets of the mitochondrial genome for a mammalian phylogeny with known subordinal relationships. Using maximum-likelihood methods, we estimated branch lengths from an 11-kb sequence of all 13 protein-coding genes and compared them with estimates from single genes (0.2-1.8 kb) and from 7 different combinations of genes (2-3.5 kb). For each sequence, we separated the component of the log-likelihood deviation due to branch length differences associated with alternative topologies from that due to those that are independent of the topology. Even among the sequences that recovered the same tree topology, some produced significantly better branch length estimates than others did. The combination of correct topology and significantly better branch length estimation suggests that these gene combinations may prove useful in estimating phylogenetic relationships for mammalian divergences below the ordinal level. Thus, the proper choice of genes to sequence is a critical factor for reliable estimation of evolutionary history from molecular data.     

Environmental selection is a main driver of divergence in house sparrows (Passer domesticus) in Romania and Bulgaria

Both neutral and adaptive evolutionary processes can cause population divergence, but their relative contributions remain unclear. We investigated the roles of these processes in population divergence in house sparrows (Passer domesticus) from Romania and Bulgaria, regions characterized by high landscape heterogeneity compared to Western Europe. We asked whether morphological divergence, complemented with genetic data in this human commensal species, was best explained by environmental variation, geographic distance, or landscape resistance—the effort it takes for an individual to disperse from one location to the other—caused by either natural or anthropogenic barriers. Using generalized dissimilarity modeling, a matrix regression technique that fits biotic beta diversity to both environmental predictors and geographic distance, we found that a small set of climate and vegetation variables explained up to ~30% of the observed divergence, whereas geographic and resistance distances played much lesser roles. Our results are consistent with signals of selection on morphological traits and of isolation by adaptation in genetic markers, suggesting that selection by natural environmental conditions shapes population divergence in house sparrows. Our study thus contributes to a growing body of evidence that adaptive evolution may be a major driver of diversification.     

<Emphasis Type="Italic">Wolbachia</Emphasis> Strains Typing in Different Geographic Population Spider, <Emphasis Type="Italic">Hylyphantes Graminicola</Emphasis> (Linyphiidae)

The Wolbachia endosymbiont of spiders has not been extensively examined. In order to investigate the distribution, evolutionary history, and reproductive phenotype of Wolbachia in spiders in China, we tested 11 geographic populations of Hylyphantes graminicola. Wolbachia infection has been detected in each population. 10 Wolbachia strains have been characterized by multilocus sequence typing (MLST). Phylogenetic analyses indicated that eight Wolbachia strains in H. graminicola belonged to supergroup B, and two belonged to supergroup A. No correlation existed between Wolbachia diversity and host’s geographic distance. The significant correlation was observed between pairwise distance of H. graminicola COI and genetic divergence of associated Wolbachia strains. We also found that Wolbachia infection frequencies in hosts varied over geographic space.     

Foraging at the thermal limit: burrowing spiders (Seothyra,Eresidae) in the Namib desert dunes

Summary In the Namib Desert dunes, the web of Seothyra sp. (Eresidae) comprises sticky silk lining the edges of a horizontal mat on the sand surface. The spider sits in a silk-lined burrow attached to the mat. Arthropods become entangled in the sticky silk of the mat and are attacked and pulled into the burrow by the spider. We investigated the influence of sand surface temperature on the activity of spiders during the summer. We determined the range of thermal conditions encountered by spiders, their temperature tolerance and the influence of temperature on foraging activity and prey handling behavior. The environmental temperatures available to Seothyra vary from 17–33° C at the coolest time of day to 33–73° C at the hottest. When prevented from retreating into burrows, spiders showed signs of thermal stress at about 49° C, whereas unrestrained spiders continued to forage at web temperatures above 65° C by moving between the hot surface mat and the cooler burrow. Spiders responded quicker to prey stimuli during the hot hours of the day and completed prey capture sequences in significantly less time at surface temperatures above 49° C than below. Furthermore, captured arthropods succumbed more quickly at high surface temperatures. Our study supports the hypothesis that web design and thermoregulatory behaviors enable Seothyra to hunt under extreme thermal conditions.     

Evolution of sociality in spiders leads to depleted genomic diversity at both population and species levels

Across several animal taxa, the evolution of sociality involves a suite of characteristics, a “social syndrome,” that includes cooperative breeding, reproductive skew, primary female‐biased sex ratio, and the transition from outcrossing to inbreeding mating system, factors that are expected to reduce effective population size (Ne). This social syndrome may be favoured by short‐term benefits but come with long‐term costs, because the reduction in Ne amplifies loss of genetic diversity by genetic drift, ultimately restricting the potential of populations to respond to environmental change. To investigate the consequences of this social life form on genetic diversity, we used a comparative RAD‐sequencing approach to estimate genomewide diversity in spider species that differ in level of sociality, reproductive skew and mating system. We analysed multiple populations of three independent sister‐species pairs of social inbreeding and subsocial outcrossing Stegodyphus spiders, and a subsocial outgroup. Heterozygosity and within‐population diversity were sixfold to 10‐fold lower in social compared to subsocial species, and demographic modelling revealed a tenfold reduction in Ne of social populations. Species‐wide genetic diversity depends on population divergence and the viability of genetic lineages. Population genomic patterns were consistent with high lineage turnover, which homogenizes the genetic structure that builds up between inbreeding populations, ultimately depleting genetic diversity at the species level. Indeed, species‐wide genetic diversity of social species was 5–8 times lower than that of subsocial species. The repeated evolution of species with this social syndrome is associated with severe loss of genomewide diversity, likely to limit their evolutionary potential.     

Among‐lake reciprocal transplants induce convergent expression of immune genes in threespine stickleback

Geographic variation in parasite communities can drive evolutionary divergence in host immune genes. However, biotic and abiotic environmental variation can also induce plastic differences in immune function among populations. At present, there is little information concerning the relative magnitudes of heritable vs. induced immune divergence in natural populations. We examined immune gene expression profiles of threespine stickleback (Gasterosteus aculeatus) from six lakes on Vancouver Island, British Columbia. Parasite community composition differs between lake types (large or small, containing limnetic‐ or benthic‐like stickleback) and between watersheds. We observed corresponding differences in immune gene expression profiles among wild‐caught stickleback, using a set of seven immune genes representing distinct branches of the immune system. To evaluate the role of environmental effects on this differentiation, we experimentally transplanted wild‐caught fish into cages in their native lake, or into a nearby foreign lake. Transplanted individuals' immune gene expression converged on patterns typical of their destination lake, deviating from their native expression profile. Transplant individuals' source population had a much smaller effect, suggesting relatively weak genetic underpinning of population differences in immunity, as viewed through gene expression. This strong environmental regulation of immune gene expression provides a counterpoint to the large emerging literature documenting microevolution and genetic diversification of immune function. Our findings illustrate the value of studying immunity in natural environmental settings where the immune system has evolved and actively functions.     

Natural selection contributes to geographic patterns of thermal plasticity in Plantago lanceolata

A long‐standing debate in evolutionary biology concerns the relative importance of different evolutionary forces in explaining phenotypic diversification at large geographic scales. For example, natural selection is typically assumed to underlie divergence along environmental gradients. However, neutral evolutionary processes can produce similar patterns. We collected molecular genetic data from 14 European populations of Plantago lanceolata to test the contributions of natural selection versus neutral evolution to population divergence in temperature‐sensitive phenotypic plasticity of floral reflectance. In P. lanceolata, reflectance plasticity is positively correlated with latitude/altitude. We used population pairwise comparisons between neutral genetic differentiation (F_ST and Jost's D) and phenotypic differentiation (P_ST) to assess the contributions of geographic distance and environmental parameters of the reproductive season in driving population divergence. Data are consistent with selection having shaped large‐scale geographic patterns in thermal plasticity. The aggregate pattern of P_ST versus F_ST was consistent with divergent selection. F_ST explained thermal plasticity differences only when geographic distance was not included in the model. Differences in the extent of cool reproductive season temperatures, and not overall temperature variation, explained plasticity differences independent of distance. Results are consistent with the hypothesis that thermal plasticity is adaptive where growing seasons are shorter and cooler, that is, at high latitude/altitude.     

Spatial evolutionary and ecological vicariance analysis (SEEVA), a novel approach to biogeography and speciation research,with an example from Brazilian Gentianaceae

Aim Spatial evolutionary and ecological vicariance analysis (SEEVA) is a simple analytical method that evaluates environmental or ecological divergence associated with evolutionary splits. It integrates evolutionary hypotheses, phylogenetic data, and spatial, temporal, environmental and geographical information to elucidate patterns. Using a phylogeny of Prepusa Mart. and Senaea Taub. (Angiospermae: Gentianaceae), SEEVA is used to describe the radiation and ecological patterns of this basal gentian group across south‐eastern Brazil. Location Latin America, global. Methods Environmental data for 151 geolocated botanical collections, associated with specimens from seven species, were compiled with Arc GIS, and were matched with geolocated base layers of eight climatological variables, as well as one each of geological, soil type, elevational and vegetation variables. Sister groups were defined on the basis of the six nested nodes that defined the phylogenetic tree of these two genera. A (0, 1)‐scaled divergence index (D) was defined and tested for each of 12 environmental and for each of the six phylogenetic nodes, by means of contingency analyses. We contrast divergence indices of nested clades, allopatric and sympatric sister clades. Results The level of ecological divergence between sister clades/species, defined in terms of D measures, was substantial for five of six nodes, with 21 of 72 environmental comparisons having D > 0.75. Soil types and geological age of bedrock were strongly divergent only for basal nodes in the phylogeny, by contrast with temperature and precipitation, which exhibited strong divergence at all nodes. There has been strong divergence and progressive occupation of wetter and colder habitats throughout the history of Prepusa. Nodes separating allopatric sister clades exhibited larger niche divergence than did those separating sympatric sister clades. Main conclusions SEEVA provides a multi‐source, direct analysis method for correlating field collections, phylogenetic hypotheses, species distributions and georeferenced environmental data. Using SEEVA, it was possible to quantify and test the divergence between sister lineages, illustrating both niche conservatism and ecological specialization. SEEVA permits elucidation of historical and ecological vicariance for evolutionary lineages, and is amenable to wide application, taxonomically, geographically and ecologically.