Evolution at two time frames: Polymorphisms from an ancient singular divergence event fuel contemporary parallel evolution.

When environments change, populations may adapt surprisingly fast, repeatedly and even at microgeographic scales. There is increasing evidence that such cases of rapid parallel evolution are fueled by standing genetic variation, but the source of this genetic variation remains poorly understood. In the saltmarsh beetle Pogonus chalceus, short-winged 'tidal' and long-winged 'seasonal' ecotypes have diverged in response to contrasting hydrological regimes and can be repeatedly found along the Atlantic European coast. By analyzing genomic variation across the beetles' distribution, we reveal that alleles selected in the tidal ecotype are spread across the genome and evolved during a singular and, likely, geographically isolated divergence event, within the last 190 Kya. Due to subsequent admixture, the ancient and differentially selected alleles are currently polymorphic in most populations across its range, which could potentially allow for the fast evolution of one ecotype from a small number of random individuals, as low as 5 to 15, from a population of the other ecotype. Our results suggest that cases of fast parallel ecological divergence can be the result of evolution at two different time frames: divergence in the past, followed by repeated selection on the same divergently evolved alleles after admixture. These findings highlight the importance of an ancient and, likely, allopatric divergence event for driving the rate and direction of contemporary fast evolution under gene flow. This mechanism is potentially driven by periods of geographic isolation imposed by large-scale environmental changes such as glacial cycles.

First genomic study on Lake Tanganyika sprat Stolothrissa tanganicae: a lack of population structure calls for integrated management of this important fisheries target species.

BACKGROUND: Clupeid fisheries in Lake Tanganyika (East Africa) provide food for millions of people in one of the world's poorest regions. Due to climate change and overfishing, the clupeid stocks of Lake Tanganyika are declining. We investigate the population structure of the Lake Tanganyika sprat Stolothrissa tanganicae, using for the first time a genomic approach on this species. This is an important step towards knowing if the species should be managed separately or as a single stock. Population structure is important for fisheries management, yet understudied for many African freshwater species. We hypothesize that distinct stocks of S. tanganicae could be present due to the large size of the lake (isolation by distance), limnological variation (adaptive evolution), or past separation of the lake (historical subdivision). On the other hand, high mobility of the species and lack of obvious migration barriers might have resulted in a homogenous population. RESULTS: We performed a population genetic study on wild-caught S. tanganicae through a combination of mitochondrial genotyping (96 individuals) and RAD sequencing (83 individuals). Samples were collected at five locations along a north-south axis of Lake Tanganyika. The mtDNA data had low global FST and, visualised in a haplotype network, did not show phylogeographic structure. RAD sequencing yielded a panel of 3504 SNPs, with low genetic differentiation (FST = 0.0054; 95% CI: 0.0046-0.0066). PCoA, fineRADstructure and global FST suggest a near-panmictic population. Two distinct groups are apparent in these analyses (FST = 0.1338 95% CI: 0.1239,0.1445), which do not correspond to sampling locations. Autocorrelation analysis showed a slight increase in genetic difference with increasing distance. No outlier loci were detected in the RADseq data. CONCLUSION: Our results show at most very weak geographical structuring of the stock and do not provide evidence for genetic adaptation to historical or environmental differences over a north-south axis. Based on these results, we advise to manage the stock as one population, integrating one management strategy over the four riparian countries. These results are a first comprehensive study on the population structure of these important fisheries target species, and can guide fisheries management.

Chromosomal inversions from an initial ecotypic divergence drive a gradual repeated radiation of Galápagos beetles.

Island faunas exhibit some of the most iconic examples where similar forms repeatedly evolve within different islands. Yet, whether these deterministic evolutionary trajectories within islands are driven by an initial, singular divergence and the subsequent exchange of individuals and adaptive genetic variation between islands remains unclear. Here, we study a gradual, repeated evolution of low-dispersive highland ecotypes from a dispersive lowland ecotype of Calosoma beetles along the island progression of the Galápagos. We show that repeated highland adaptation involved selection on multiple shared alleles within extensive chromosomal inversions that originated from an initial adaptation event on the oldest island. These highland inversions first spread through dispersal of highland individuals. Subsequent admixture with the lowland ecotype resulted in polymorphic dispersive populations from which the highland populations evolved on the youngest islands. Our findings emphasize the significance of an ancient divergence in driving repeated evolution and highlight how a mixed contribution of inter-island colonization and within-island evolution can shape parallel species communities.

A masculinizing supergene underlies an exaggerated male reproductive morph in a spider.

In many species, individuals can develop into strikingly different morphs, which are determined by a simple Mendelian locus. How selection shapes loci that control complex phenotypic differences remains poorly understood. In the spider Oedothorax gibbosus, males develop either into a 'hunched' morph with conspicuous head structures or as a fast-developing 'flat' morph with a female-like appearance. We show that the hunched-determining allele contains a unique genomic fragment of approximately 3 megabases that is absent in the flat-determining allele. This fragment comprises dozens of genes that duplicated from genes found at the same as well as different chromosomes. All functional duplicates, including a duplicate of the key sexual differentiation regulatory gene doublesex, show male-specific expression, which illustrates their integrated role as a masculinizing supergene. Our findings demonstrate how extensive indel polymorphisms and duplications of regulatory genes may contribute to the evolution of co-adapted gene clusters, sex-limited reproductive morphs and the enigmatic evolution of exaggerated sexual traits in general.

Large-scale and small-scale population genetic structure of the medically important gastropod species Bulinus truncatus (Gastropoda, Heterobranchia).

BACKGROUND: Gastropod snails remain strongly understudied, despite their important role in transmitting parasitic diseases. Knowledge of their distribution and population dynamics increases our understanding of the processes driving disease transmission. We report the first study to use high-throughput sequencing (HTS) to elucidate the population genetic structure of the hermaphroditic snail Bulinus truncatus (Gastropoda, Heterobranchia) on a regional (17-150 km) and inter-regional (1000-5400 km) scale. This snail species acts as an intermediate host of Schistosoma haematobium and Schistosoma bovis, which cause human and animal schistosomiasis respectively. METHODS: Bulinus truncatus snails were collected in Senegal, Cameroon, Egypt and France and identified through DNA barcoding. A single-end genotyping-by-sequencing (GBS) library, comprising 87 snail specimens from the respective countries, was built and sequenced on an Illumina HiSeq 2000 platform. Reads were mapped against S. bovis and S. haematobium reference genomes to identify schistosome infections, and single nucleotide polymorphisms (SNPs) were scored using the Stacks pipeline. These SNPs were used to estimate genetic diversity, assess population structure and construct phylogenetic trees of B. truncatus. RESULTS: A total of 10,750 SNPs were scored and used in downstream analyses. The phylogenetic analysis identified five clades, each consisting of snails from a single country but with two distinct clades within Senegal. Genetic diversity was low in all populations, reflecting high selfing rates, but varied between locations due to habitat variability. Significant genetic differentiation and isolation by distance patterns were observed at both spatial scales, indicating that gene flow is not strong enough to counteract the effects of population bottlenecks, high selfing rates and genetic drift. Remarkably, the population genetic differentiation on a regional scale (i.e. within Senegal) was as large as that between populations on an inter-regional scale. The blind GBS technique was able to pick up parasite DNA in snail tissue, demonstrating the potential of HTS techniques to further elucidate the role of snail species in parasite transmission. CONCLUSIONS: HTS techniques offer a valuable toolbox to further investigate the population genetic patterns of intermediate schistosome host snails and the role of snail species in parasite transmission.