Divergent Macroparasite Infections in Parapatric Swiss Lake-Stream Pairs of Threespine Stickleback (Gasterosteus aculeatus)

Spatial heterogeneity in diversity and intensity of parasitism is a typical feature of most host-parasite interactions, but understanding of the evolutionary implications of such variation is limited. One possible outcome of infection heterogeneities is parasite-mediated divergent selection between host populations, ecotypes or species which may facilitate the process of ecological speciation. However, very few studies have described infections in population-pairs along the speciation continuum from low to moderate or high degree of genetic differentiation that would address the possibility of parasite-mediated divergent selection in the early stages of the speciation process. Here we provide an example of divergent parasitism in freshwater fish ecotypes by examining macroparasite infections in threespine stickleback (Gasterosteus aculeatus) of four Swiss lake systems each harbouring parapatric lake-stream ecotype pairs. We demonstrate significant differences in infections within and between the pairs that are driven particularly by the parasite taxa transmitted to fish from benthic invertebrates. The magnitude of the differences tended to correlate positively with the extent of neutral genetic differentiation between the parapatric lake and stream populations of stickleback, whereas no such correlation was found among allopatric populations from similar or contrasting habitats. This suggests that genetic differentiation is unrelated to the magnitude of parasite infection contrasts when gene flow is constrained by geographical barriers while in the absence of physical barriers, genetic differentiation and the magnitude of differences in infections tend to be positively correlated.

Structure and dynamics of hybrid zones at different stages of speciation in the common vole (Microtus arvalis)

The genetic structure and dynamics of hybrid zones provide crucial information for understanding the processes and mechanisms of evolutionary divergence and speciation. In general, higher levels of evolutionary divergence between taxa are more likely to be associated with reproductive isolation and may result in suppressed or strongly restricted hybridization. In this study, we examined two secondary contact zones between three deep evolutionary lineages in the common vole (Microtus arvalis). Differences in divergence times between the lineages can shed light on different stages of reproductive isolation and thus provide information on the ongoing speciation process in M. arvalis. We examined more than 800 individuals for mitochondrial (mtDNA), Y-chromosome and autosomal markers and used assignment and cline analysis methods to characterize the extent and direction of gene flow in the contact zones. Introgression of both autosomal and mtDNA markers in a relatively broad area of admixture indicates selectively neutral hybridization between the least-divergent lineages (Central and Eastern) without evidence for partial reproductive isolation. In contrast, a very narrow area of hybridization, shifts in marker clines and the quasi-absence of Y-chromosome introgression support a moving hybrid zone and unidirectional selection against male hybrids between the lineages with older divergence (Central and Western). Data from a replicate transect further support non-neutral processes in this hybrid zone and also suggest a role for landscape history in the movement and shaping of geneflow profiles.

Color polymorphism and sex ratio distortion in a cichlid fish as an incipient stage in sympatric speciation by sexual selection

We investigated a Lake Victoria cichlid with a complex colour polymorphism that apparently represents one original species and two incipient species, all of which are sympatric. In laboratory breeding experiments we observed sex ratio distortion in certain matings between original and incipient species. Mate choice experiments show that males of the incipient species exhibit mating preferences against the original species, and males and females of the original species exhibit strong mating preferences against the incipient species. Mating preferences might evolve by sex ratio selection to avoid matings with distorted progeny sex ratios. Phenotype frequencies in nature suggest that mating preferences translate into mating frequencies, thus restricting gene flow and exerting disruptive sexual selection between the original and incipient species. The incipient species do not differ in morphology or ecology from the original species, implying that colour polymorphism, associated with sex ratio distortion, can be an incipient stage in sympatric speciation, and that disruption of gene flow can precede ecological differentiation

Direct mate choice maintains diversity among sympatric cichlids in Lake Victoria

Mate choice may play an important role in animal speciation. The haplochromine cichlids of Lake Victoria are suitable to test this hypothesis. Diversity in ecology, coloration and anatomy evolved in these fish faster than postzygotic barriers to gene flow, and little is known about how this diversity is maintained. It was tested whether recognizable forms are selection-maintained morphs or reproductively isolated species by investigating in the field reproductive timing, location of spawning sites, and mate choice behaviour. There was a large interspecific overlap in timing of breeding and location of spawning sites, which was largest in members of the same genus. Behavioural mate choice of such closely related taxa was highly assortative, such that it is likely that they are sexually isolated species and that direct mate choice is the major force that directs gene flow and maintains form diversity. The results differ from what is known about recent radiations of other lacustrine fish groups where speciation seems to be driven by diverging microhabitat preferences or diverging timing of reproduction, but are in agreement with predictions from models of speciation by diverging mate preferences.

Evolutionary Genomics and Conservation of the Endangered Przewalski's Horse.

Przewalski's horses (PHs, Equus ferus ssp. przewalskii) were discovered in the Asian steppes in the 1870s and represent the last remaining true wild horses. PHs became extinct in the wild in the 1960s but survived in captivity, thanks to major conservation efforts. The current population is still endangered, with just 2,109 individuals, one-quarter of which are in Chinese and Mongolian reintroduction reserves [1]. These horses descend from a founding population of 12 wild-caught PHs and possibly up to four domesticated individuals [2-4]. With a stocky build, an erect mane, and stripped and short legs, they are phenotypically and behaviorally distinct from domesticated horses (DHs, Equus caballus). Here, we sequenced the complete genomes of 11 PHs, representing all founding lineages, and five historical specimens dated to 1878-1929 CE, including the Holotype. These were compared to the hitherto-most-extensive genome dataset characterized for horses, comprising 21 new genomes. We found that loci showing the most genetic differentiation with DHs were enriched in genes involved in metabolism, cardiac disorders, muscle contraction, reproduction, behavior, and signaling pathways. We also show that DH and PH populations split ∼45,000 years ago and have remained connected by gene-flow thereafter. Finally, we monitor the genomic impact of ∼110 years of captivity, revealing reduced heterozygosity, increased inbreeding, and variable introgression of domestic alleles, ranging from non-detectable to as much as 31.1%. This, together with the identification of ancestry informative markers and corrections to the International Studbook, establishes a framework for evaluating the persistence of genetic variation in future reintroduced populations.

Tracing reinforcement through asymmetrical partner preference in the European common vole Microtus arvalis

Background The mechanistic basis of speciation and in particular the contribution of behaviour to the completion of the speciation process is often contentious. Contact zones between related taxa provide a situation where selection against hybridization might reinforce separation by behavioural mechanisms, which could ultimately fully isolate the taxa. One of the most abundant European mammals, the common vole Microtus arvalis, forms multiple natural hybrid zones where rapidly diverging evolutionary lineages meet in secondary contact. Very narrow zones of hybridization spanning only a few kilometres and sex-specific gene flow patterns indicate reduced fitness of natural hybrids and incipient speciation between some of the evolutionary lineages. In this study, we examined the contribution of behavioural mechanisms to the speciation process in these rodents by fine-mapping allopatric and parapatric populations in the hybrid zone between the Western and Central lineages and experimental testing of the partner preferences of wild, pure-bred and hybrid female common voles. Results Genetic analysis based on microsatellite markers revealed the presence of multiple parapatric and largely non-admixed populations at distances of about 10 km at the edge of the area of natural hybridization between the Western and Central lineages. Wild females from Western parapatric populations and lab-born F1 hybrids preferred males from the Western lineage whereas wild females of Central parapatric origin showed no measurable preference. Furthermore, wild and lab-born females from allopatric populations of the Western or Central lineages showed no detectable preference for males from either lineage. Conclusions The detected partner preferences are consistent with asymmetrical reinforcement of pre-mating reproductive isolation mechanisms in the European common vole and with earlier results suggesting that hybridization is more detrimental to the Western lineage. As a consequence, these differences in behaviour might contribute to a further geographical stabilization of this moving hybrid zone. Such behavioural processes could also provide a mechanistic perspective for frequently-detected asymmetrical introgression patterns in the largely allopatrically diversifying Microtus genus and other rapidly speciating rodents.

Low but contrasting neutral genetic differentiation shaped by winter temperature in European great tits

Gene flow is usually thought to reduce genetic divergence and impede local adaptation by homogenising gene pools between populations. However, evidence for local adaptation and phenotypic differentiation in highly mobile species, experiencing high levels of gene flow, is emerging. Assessing population genetic structure at different spatial scales is thus a crucial step towards understanding mechanisms underlying intraspecific differentiation and diversification. Here, we studied the population genetic structure of a highly mobile species – the great tit Parus major – at different spatial scales. We analysed 884 individuals from 30 sites across Europe including 10 close-by sites (< 50 km), using 22 microsatellite markers. Overall we found a low but significant genetic differentiation among sites (FST = 0.008). Genetic differentiation was higher, and genetic diversity lower, in south-western Europe. These regional differences were statistically best explained by winter temperature. Overall, our results suggest that great tits form a single patchy metapopulation across Europe, in which genetic differentiation is independent of geographical distance and gene flow may be regulated by environmental factors via movements related to winter severity. This might have important implications for the evolutionary trajectories of sub-populations, especially in the context of climate change, and calls for future investigations of local differences in costs and benefits of philopatry at large scales.

Genomics of Rapid Incipient Speciation in Sympatric Threespine Stickleback

Ecological speciation is the process by which reproductively isolated populations emerge as a consequence of divergent natural or ecologically-mediated sexual selection. Most genomic studies of ecological speciation have investigated allopatric populations, making it difficult to infer reproductive isolation. The few studies on sympatric ecotypes have focused on advanced stages of the speciation process after thousands of generations of divergence. As a consequence, we still do not know what genomic signatures of the early onset of ecological speciation look like. Here, we examined genomic differentiation among migratory lake and resident stream ecotypes of threespine stickleback reproducing in sympatry in one stream, and in parapatry in another stream. Importantly, these ecotypes started diverging less than 150 years ago. We obtained 34,756 SNPs with restriction-site associated DNA sequencing and identified genomic islands of differentiation using a Hidden Markov Model approach. Consistent with incipient ecological speciation, we found significant genomic differentiation between ecotypes both in sympatry and parapatry. Of 19 islands of differentiation resisting gene flow in sympatry, all were also differentiated in parapatry and were thus likely driven by divergent selection among habitats. These islands clustered in quantitative trait loci controlling divergent traits among the ecotypes, many of them concentrated in one region with low to intermediate recombination. Our findings suggest that adaptive genomic differentiation at many genetic loci can arise and persist in sympatry at the very early stage of ecotype divergence, and that the genomic architecture of adaptation may facilitate this.

Cryptic invasion drives phenotypic changes in central European threespine stickleback

Cryptic invasions are commonly associated with genetic changes of the native species or genetic lineage that the invaders replace. Phenotypic shifts resulting from cryptic invasions are less commonly reported given the relative paucity of historical specimens that document such phenotypic changes. Here, I study such a case in two populations of threespine stickleback from central Europe, comparing contemporary patterns of gene flow with phenotypic changes between historical and contemporary population samples. I find gene flow from an invasive lineage to be associated with significant phenotypic changes, where the degree of phenotypic change corresponds with the level of gene flow that a population receives. These findings underline the utility of combining genetic approaches with phenotypic data to estimate the impact of gene flow in systems where anthropogenic alterations have removed former geographic barriers promoting cryptic invasions.

Genetic diversity and differentiation follow secondary succession in a multi-species study on woody plants from subtropical China

Aims: Species diversity and genetic diversity may be affected in parallel by similar environmental drivers. However, genetic diversity may also be affected independently by habitat characteristics. We aim at disentangling relationships between genetic diversity, species diversity and habitat characteristics of woody species in subtropical forest. Methods: We studied 11 dominant tree and shrub species in 27 plots in Gutianshan, China, and assessed their genetic diversity (Ar) and population differentiation (F’ST) with microsatellite markers. We tested if Ar and population specific F’ST were correlated to local species diversity and plot characteristics. Multi-model inference and model averaging were used to determine the relative importance of each predictor. Additionally we tested for isolation-by-distance and isolation-by-elevation by regressing pairwise F’ST against pairwise spatial and elevational distances. Important findings: Genetic diversity was not related to species diversity for any of the study species. Thus, our results do not support joint effects of habitat characteristics on these two levels of biodiversity. Instead, genetic diversity in two understory shrubs, Rhododendron simsii and Vaccinium carlesii, was affected by plot age with decreasing genetic diversity in successionally older plots. Population differentiation increased with plot age in Rhododendron simsii and Lithocarpus glaber. This shows that succession can reduce genetic diversity within, and increase genetic diversity between populations. Furthermore, we found four cases of isolation-by-distance and two cases of isolation-by-elevation. The former indicates inefficient pollen and seed dispersal by animals whereas the latter might be due to phenological asynchronies. These patterns indicate that succession can affect genetic diversity without parallel effects on species diversity and that gene flow in a continuous subtropical forest can be restricted even at a local scale.

Ecosystem size matters: The dimensionality of intralacustrine diversification in Icelandic stickleback is predicted by lake size

Cases of evolutionary diversification can be characterized along a continuum from weak to strong genetic and phenotypic differentiation. Several factors may facilitate or constrain the differentiation process. Comparative analyses of replicates of the same taxon at different stages of differentiation can be useful to identify these factors. We estimated the number of distinct phenotypic groups in threespine stickleback populations from nine lakes in Iceland and in one marine population. Using the inferred number of phenotypic groups in each lake, genetic divergence from the marine population, and physical lake and landscape variables, we tested if ecosystem size, approximated by lake size and depth, or isolation from the ancestral marine gene pool predict the occurrence and the extent of phenotypic and genetic diversification within lakes. We find intralacustrine phenotypic diversification to be the rule rather than the exception, occurring in all but the youngest lake population and being manifest in ecologically important phenotypic traits. Neutral genetic data further indicates non-random mating in four out of nine studied lakes, and restricted gene flow between sympatric phenotypic groups in two. Although neither the phenotypic variation nor the number of intralacustrine phenotypic groups were associated with any of our environmental variables, the number of phenotypic traits that were differentiated was significantly positively related to lake size, and evidence for restricted gene flow between sympatric phenotypic groups was only found in the largest lakes where trait specific phenotypic differentiation was highest.

Using multiple landscape genetic approaches to test the validity of genetic clusters in a species characterized by an isolation-by-distance pattern

Bayesian clustering methods are typically used to identify barriers to gene flow, but they are prone to deduce artificial subdivisions in a study population characterized by an isolation-by-distance pattern (IbD). Here we analysed the landscape genetic structure of a population of wild boars (Sus scrofa) from south-western Germany. Two clustering methods inferred the presence of the same genetic discontinuity. However, the population in question was characterized by a strong IbD pattern. While landscape-resistance modelling failed to identify landscape features that influenced wild boar movement, partial Mantel tests and multiple regression of distance matrices (MRDMs) suggested that the empirically inferred clusters were separated by a genuine barrier. When simulating random lines bisecting the study area, 60% of the unique barriers represented, according to partial Mantel tests and MRDMs, significant obstacles to gene flow. By contrast, the random-lines simulation showed that the boundaries of the inferred empirical clusters corresponded to the most important genetic discontinuity in the study area. Given the degree of habitat fragmentation separating the two empirical partitions, it is likely that the clustering programs correctly identified a barrier to gene flow. The differing results between the work published here and other studies suggest that it will be very difficult to draw general conclusions about habitat permeability in wild boar from individual studies.

Pulsatile shear and Gja5 modulate arterial identity and remodeling events during flow-driven arteriogenesis

In the developing chicken embryo yolk sac vasculature, the expression of arterial identity genes requires arterial hemodynamic conditions. We hypothesize that arterial flow must provide a unique signal that is relevant for supporting arterial identity gene expression and is absent in veins. We analyzed factors related to flow, pressure and oxygenation in the chicken embryo vitelline vasculature in vivo. The best discrimination between arteries and veins was obtained by calculating the maximal pulsatile increase in shear rate relative to the time-averaged shear rate in the same vessel: the relative pulse slope index (RPSI). RPSI was significantly higher in arteries than veins. Arterial endothelial cells exposed to pulsatile shear in vitro augmented arterial marker expression as compared with exposure to constant shear. The expression of Gja5 correlated with arterial flow patterns: the redistribution of arterial flow provoked by vitelline artery ligation resulted in flow-driven collateral arterial network formation and was associated with increased expression of Gja5. In situ hybridization in normal and ligation embryos confirmed that Gja5 expression is confined to arteries and regulated by flow. In mice, Gja5 (connexin 40) was also expressed in arteries. In the adult, increased flow drives arteriogenesis and the formation of collateral arterial networks in peripheral occlusive diseases. Genetic ablation of Gja5 function in mice resulted in reduced arteriogenesis in two occlusion models. We conclude that pulsatile shear patterns may be central for supporting arterial identity, and that arterial Gja5 expression plays a functional role in flow-driven arteriogenesis.

Intussusceptive angiogenesis: pillars against the blood flow

Adaptation of vascular networks to functional demands needs vessel growth, vessel regression and vascular remodelling. Biomechanical forces resulting from blood flow play a key role in these processes. It is well-known that metabolic stimuli, mechanical forces and flow patterns can affect gene expression and remodelling of vascular networks in different ways. For instance, in the sprouting type of angiogenesis related to hypoxia, there is no blood flow in the rising capillary sprout. In contrast, it has been shown that an increase of wall shear stress initiates the splitting type of angiogenesis in skeletal muscle. Otherwise, during development, both sprouting and intussusception act in parallel in building the vascular network, although with differences in spatiotemporal distribution. Thereby, in addition to regulatory molecules, flow dynamics support the patterning and remodelling of the rising vascular tree. Herewith, we present an overview of angiogenic processes with respect to intussusceptive angiogenesis as related to local haemodynamics.

A KEL gene encoding serine at position 193 of the Kell glycoprotein results in expression of KEL1 antigen

BACKGROUND: The KEL2/KEL1 (k/K) blood group polymorphism represents 578C>T in the KEL gene and Thr193Met in the Kell glycoprotein. Anti-KEL1 can cause severe hemolytic disease of the fetus and newborn. Molecular genotyping for KEL*1 is routinely used for assessing whether a fetus is at risk. Red blood cells (RBCs) from a KEL:1 blood donor (D1) were found to have abnormal KEL1 expression during evaluation of anti-KEL1 reagents. STUDY DESIGN AND METHODS: Kell genotyping methods, including KEL exon 6 direct sequencing, were applied. KEL cDNA from D1 was sequenced. Flow cytometry was used to assess KEL1 and KEL2 RBC expression. RESULTS: RBCs from the donor, her mother, and an unrelated donor gave weak or negative reactions with some anti-KEL1 reagents. Other Kell-system antigens appeared normal. The three individuals were homozygous for KEL C578 (KEL*2) but heterozygous for a 577A>T transversion, encoding Ser193. They appeared to be KEL*2 homozygotes by routine genotyping methods. Flow cytometry revealed weak KEL1 expression and normal KEL2, similar to that of KEL*2 homozygotes. CONCLUSION: Ser193 in the Kell glycoprotein appears to result in expression of abnormal KEL1, in addition to KEL2. The mutation is not detected by routine Kell genotyping methods and, because of unpredicted KEL1 expression, could lead to a misdiagnosis.