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.

Sex-specific clines support incipient speciation in a common European mammal

Hybrid zones provide excellent opportunities to study processes and mechanisms underlying reproductive isolation and speciation. Here we investigated sex-specific clines of molecular markers in hybrid zones of morphologically cryptic yet genetically highly-diverged evolutionary lineages of the European common vole (Microtus arvalis). We analyzed the position and width of four secondary contact zones along three independent transects in the region of the Alps using maternally (mitochondrial DNA) and paternally (Y-chromosome) inherited genetic markers. Given male-biased dispersal in the common vole, a selectively neutral secondary contact would show broader paternal marker clines than maternal ones. In a selective case, for example, involving a form of Haldane’s rule, Y-chromosomal clines would not be expected to be broader than maternal markers because they are transmitted by the heterogametic sex and thus gene flow would be restricted. Consistent with the selective case, paternal clines were significantly narrower or at most equal in width to maternal clines in all contact zones. In addition, analyses using maximum likelihood cline-fitting detected a shift of paternal relative to maternal clines in three of four contact zones. These patterns suggest that processes at the contact zones in the common vole are not selectively neutral, and that partial reproductive isolation is already established between these evolutionary lineages. We conclude that hybrid zone movement, sexual selection and/or genetic incompatibilities are likely associated with an unusual unidirectional manifestation of Haldane’s rule in this common European mammal.

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.

Thermal adaptation and ecological speciation

Ecological speciation is defined as the emergence of reproductive isolation as a direct or indirect consequence of divergent ecological adaptation. Several empirical examples of ecological speciation have been reported in the literature which very often involve adaptation to biotic resources. In this review, we investigate whether adaptation to different thermal habitats could also promote speciation and try to assess the importance of such processes in nature. Our survey of the literature identified 16 animal and plant systems where divergent thermal adaptation may underlie (partial) reproductive isolation between populations or may allow the stable coexistence of sibling taxa. In many of the systems, the differentially adapted populations have a parapatric distribution along an environmental gradient. Isolation often involves extrinsic selection against locally maladapted parental or hybrid genotypes, and additional pre- or postzygotic barriers may be important. Together, the identified examples strongly suggest that divergent selection between thermal environments is often strong enough to maintain a bimodal genotype distribution upon secondary contact. What is less clear from the available data is whether it can also be strong enough to allow ecological speciation in the face of gene flow through reinforcement-like processes. It is possible that intrinsic features of thermal gradients or the genetic basis of thermal adaptation make such reinforcement-like processes unlikely but it is equally possible that pertinent systems are understudied. Overall, our literature survey highlights (once again) the dearth of studies that investigate similar incipient species along the continuum from initial divergence to full reproductive isolation and studies that investigate all possible reproductive barriers in a given system.

A genetically explicit model of speciation by sensory drive within a continuous population in aquatic environments

BACKGROUND: The sensory drive hypothesis predicts that divergent sensory adaptation in different habitats may lead to premating isolation upon secondary contact of populations. Speciation by sensory drive has traditionally been treated as a special case of speciation as a byproduct of adaptation to divergent environments in geographically isolated populations. However, if habitats are heterogeneous, local adaptation in the sensory systems may cause the emergence of reproductively isolated species from a single unstructured population. In polychromatic fishes, visual sensitivity might become adapted to local ambient light regimes and the sensitivity might influence female preferences for male nuptial color. In this paper, we investigate the possibility of speciation by sensory drive as a byproduct of divergent visual adaptation within a single initially unstructured population. We use models based on explicit genetic mechanisms for color vision and nuptial coloration. RESULTS: We show that in simulations in which the adaptive evolution of visual pigments and color perception are explicitly modeled, sensory drive can promote speciation along a short selection gradient within a continuous habitat and population. We assumed that color perception evolves to adapt to the modal light environment that individuals experience and that females prefer to mate with males whose nuptial color they are most sensitive to. In our simulations color perception depends on the absorption spectra of an individual's visual pigments. Speciation occurred most frequently when the steepness of the environmental light gradient was intermediate and dispersal distance of offspring was relatively small. In addition, our results predict that mutations that cause large shifts in the wavelength of peak absorption promote speciation, whereas we did not observe speciation when peak absorption evolved by stepwise mutations with small effect. CONCLUSION: The results suggest that speciation can occur where environmental gradients create divergent selection on sensory modalities that are used in mate choice. Evidence for such gradients exists from several animal groups, and from freshwater and marine fishes in particular. The probability of speciation in a continuous population under such conditions may then critically depend on the genetic architecture of perceptual adaptation and female mate choice.

Phenotypic divergence but not genetic distance predicts assortative mating among species of a cichlid fish radiation

The hypothesis of ecological divergence giving rise to premating isolation in the face of gene flow is controversial. However, this may be an important mechanism to explain the rapid multiplication of species during adaptive radiation following the colonization of a new environment when geographical barriers to gene flow are largely absent but underutilized niche space is abundant. Using cichlid fish, we tested the prediction of ecological speciation that the strength of premating isolation among species is predicted by phenotypic rather than genetic distance. We conducted mate choice experiments between three closely related, sympatric species of a recent radiation in Lake Mweru (Zambia/DRC) that differ in habitat use and phenotype, and a distantly related population from Lake Bangweulu that resembles one of the species in Lake Mweru. We found significant assortative mating among all closely related, sympatric species that differed phenotypically, but none between the distantly related allopatric populations of more similar phenotype. Phenotypic distance between species was a good predictor of the strength of premating isolation, suggesting that assortative mating can evolve rapidly in association with ecological divergence during adaptive radiation. Our data also reveals that distantly related allopatric populations that have not diverged phenotypically, may hybridize when coming into secondary contact, e.g. upon river capture because of diversion of drainage systems.

Genomic signatures of relaxed disruptive selection associated with speciation reversal in whitefish

Background: Speciation reversal: the erosion of species differentiation via an increase in introgressive hybridization due to the weakening of previously divergent selection regimes, is thought to be an important, yet poorly understood, driver of biodiversity loss. Our study system, the Alpine whitefish (Coregonus spp.) species complex is a classic example of a recent postglacial adaptive radiation: forming an array of endemic lake flocks, with the independent origination of similar ecotypes among flocks. However, many of the lakes of the Alpine radiation have been seriously impacted by anthropogenic nutrient enrichment, resulting in a collapse in neutral genetic and phenotypic differentiation within the most polluted lakes. Here we investigate the effects of eutrophication on the selective forces that have shaped this radiation, using population genomics. We studied eight sympatric species assemblages belonging to five independent parallel adaptive radiations, and one species pair in secondary contact. We used AFLP markers, and applied FST outlier (BAYESCAN, DFDIST) and logistic regression analyses (MATSAM), to identify candidate regions for disruptive selection in the genome and their associations with adaptive traits within each lake flock. The number of outlier and adaptive trait associated loci identified per lake were then regressed against two variables (historical phosphorus concentration and contemporary oxygen concentration) representing the strength of eutrophication. Results: Whilst we identify disruptive selection candidate regions in all lake flocks, we find similar trends, across analysis methods, towards fewer disruptive selection candidate regions and fewer adaptive trait/candidate loci associations in the more polluted lakes. Conclusions: Weakened disruptive selection and a concomitant breakdown in reproductive isolating mechanisms in more polluted lakes has lead to increased gene flow between coexisting Alpine whitefish species. We hypothesize that the resulting higher rates of interspecific recombination reduce either the number or extent of genomic islands of divergence surrounding loci evolving under disruptive natural selection. This produces the negative trend seen in the number of selection candidate loci recovered during genome scans of whitefish species flocks, with increasing levels of anthropogenic eutrophication: as the likelihood decreases that AFLP restriction sites will fall within regions of heightened genomic divergence and therefore be classified as FST outlier loci. This study explores for the first time the potential effects of human-mediated relaxation of disruptive selection on heterogeneous genomic divergence between coexisting species.

Putative functions of extracellular matrix glycoproteins in secondary palate morphogenesis

Cleft palate is a common birth defect in humans. Elevation and fusion of paired palatal shelves are coordinated by growth and transcription factors, and mutations in these can cause malformations. Among the effector genes for growth factor signaling are extracellular matrix (ECM) glycoproteins. These provide substrates for cell adhesion (e.g., fibronectin, tenascins), but also regulate growth factor availability (e.g., fibrillins). Cleft palate in Bmp7 null mouse embryos is caused by a delay in palatal shelf elevation. In contrast, palatal shelves of Tgf-β3 knockout mice elevate normally, but a cleft develops due to their failure to fuse. However, nothing is known about a possible functional interaction between specific ECM proteins and Tgf-β/Bmp family members in palatogenesis. To start addressing this question, we studied the mRNA and protein distribution of relevant ECM components during secondary palate development, and compared it to growth factor expression in wildtypewild type and mutant mice. We found that fibrillin-2 (but not fibrillin-1) mRNA appeared in the mesenchyme of elevated palatal shelves adjacent to the midline epithelial cells, which were positive for Tgf-β3 mRNA. Moreover, midline epithelial cells started expressing fibronectin upon contact of the two palatal shelves. These findings support the hypothesis that fibrillin-2 and fibronectin are involved in regulating the activity of Tgf-β3 at the fusing midline. In addition, we observed that tenascin-W (but not tenascin-C) was misexpressed in palatal shelves of Bmp7-deficient mouse embryos. In contrast to tenascin-C, tenascin-W secretion was strongly induced by Bmp7 in embryonic cranial fibroblasts in vitro. These results are consistent with a putative function for tenascin-W as a target of Bmp7 signaling during palate elevation. Our results indicate that distinct ECM proteins are important for morphogenesis of the secondary palate, both as downstream effectors and as regulators of Tgf-β/Bmp activity.

Tuberculosis in a Swiss army training camp: contact investigation using an Interferon gamma release assay

BACKGROUND: In tuberculosis (TB), the risk of exposure is determined mainly by the proximity to and the hours of direct contact with an infectious patient. We describe the contact investigation after detection of an infectious form of TB in a military camp using an Interferon-g-Release-Assay (IGRA, QuantiFERON-TB Gold In Tube [QTF-GIT]) eight weeks after detection of the index case. INDEX PATIENT: The index patient presented with fever, cough and weight loss in the military hospital six weeks after entering the camp. TB was suspected and anti-tuberculous therapy given immediately. Subsequently, TB was microbiologically confirmed. METHODS: Four exposure groups were formed a priori based on the proximity and the hours of direct contact to the index case. 168 (95.5%) agreed to be investigated: - Group A: sharing the same dormitory (15 persons) - Group B: same platoon, but not sharing the dormitory (20 persons) - Group C: staff and patients of the military hospital (22 persons) - Group D: other three platoons and senior military staff (111 persons). RESULTS: 34 (20.2%) out of 168 contacts tested positive in the QFT-GIT assay. For the exposure groups, the respective QFT-GIT testing results were: group A, 14/15 (93%); group B, 4/20 (20%); group C, 5/22 (22.7%); and group D, 11/111 (9.9%). No secondary TB cases were identified. CONCLUSIONS: In our study, test results show a correlation with the risk of exposure, suggesting that IGRA may be useful for the assessment of TB infection in TB contacts. The high mobility of recruits reduced traceability of contacts. In this context, QFT-GIT allowed for an efficient screening of contacts at a single time point.

The effects of impingement and dysplasia on stress distributions in the hip joint during sitting and walking: a finite element analysis

Soft tissue damage has been observed in hip joints with pathological geometries. Our primary goal was to study the relationship between morphological variations of the bony components of the hip and resultant stresses within the soft tissues of the joint during routine daily activities. The secondary goal was to find the range of morphological parameters in which stresses are minimized. Computational models of normal and pathological joints were developed based on variations of morphological parameters of the femoral head (Alpha angle) and acetabulum (CE angle). The Alpha angle was varied between 40 degrees (normal joint) and 80 degrees (cam joint). The CE angle was varied between 0 degrees (dysplastic joint) and 40 degrees (pincer joint). Dynamic loads and motions for walking and standing to sitting were applied to all joint configurations. Contact pressures and stresses were calculated and crosscompared to evaluate the influence of morphology. The stresses in the soft tissues depended strongly on the head and acetabular geometry. For the dysplastic joint, walking produced high acetabular rim stresses. Conversely, for impinging joints, standing-to-sitting activities that involved extensive motion were critical, inducing excessive distortion and shearing of the tissue-bone interface. Zones with high von Mises stresses corresponded with clinically observed damage zones in the acetabular cartilage and labrum. Hip joint morphological parameters that minimized were 20 degrees

Seismic anisotropy in the Morcles nappe shear zone: Implications for seismic imaging of crustal scale shear zones

Microstructures and textures of calcite mylonites from the Morcles nappe large-scale shear zone in southwestern Switzerland develop principally as a function of 1) extrinsic physical parameters including temperature, stress, strain, strain rate and 2) intrinsic parameters, such as mineral composition. We collected rock samples at a single location from this shear zone, on which laboratory ultrasonic velocities, texture and microstructures were investigated and quantified. The samples had different concentration of secondary mineral phases (< 5 up to 40 vol.%). Measured seismic P wave anisotropy ranges from 6.5% for polyphase mylonites (~ 40 vol.%) to 18.4% in mylonites with < 5 vol.% secondary phases. Texture strength of calcite is the main factor governing the seismic P wave anisotropy. Measured S wave splitting is generally highest in the foliation plane, but its origin is more difficult to explain solely by calcite texture. Additional texture measurements were made on calcite mylonites with low concentration of secondary phases (≤ 10 vol.%) along the metamorphic gradient of the shear zone (15 km distance). A systematic increase in texture strength is observed moving from the frontal part of the shear zone (anchimetamorphism; 280 °C) to the higher temperature, basal part (greenschist facies; 350–400 °C). Calculated P wave velocities become increasingly anisotropic towards the high-strain part of the nappe, from an average of 5.8% in the frontal part to 13.2% in the root of the basal part. Secondary phases raise an additional complexity, and may act either to increase or decrease seismic anisotropy of shear zone mylonites. In light of our findings we reinterpret the origin of some seismically reflective layers in the Grône–Zweisimmen line in southwestern Switzerland (PNR20 Swiss National Research Program). We hypothesize that reflections originate in part from the lateral variation in textural and microstructural arrangement of calcite mylonites in shear zones.

Grain coarsening in contact metamorphic carbonates: effects of second-phase particles, fluid flow and thermal perturbations

Under contact metamorphic conditions, carbonate rocks in the direct vicinity of the Adamello pluton reflect a temperature-induced grain coarsening. Despite this large-scale trend, a considerable grain size scatter occurs on the outcrop-scale indicating local influence of second-order effects such as thermal perturbations, fluid flow and second-phase particles. Second-phase particles, whose sizes range from nano- to the micron-scale, induce the most pronounced data scatter resulting in grain sizes too small by up to a factor of 10, compared with theoretical grain growth in a pure system. Such values are restricted to relatively impure samples consisting of up to 10 vol.% micron-scale second-phase particles, or to samples containing a large number of nano-scale particles. The obtained data set suggests that the second phases induce a temperature-controlled reduction on calcite grain growth. The mean calcite grain size can therefore be expressed in the form D 1⁄4 C2 eQ*/RT(dp/fp)m*, where C2 is a constant, Q* is an activation energy, T the temperature and m* the exponent of the ratio dp/fp, i.e. of the average size of the second phases divided by their volume fraction. However, more data are needed to obtain reliable values for C2 and Q*. Besides variations in the average grain size, the presence of second-phase particles generates crystal size distribution (CSD) shapes characterized by lognormal distributions, which differ from the Gaussian-type distributions of the pure samples. In contrast, fluid-enhanced grain growth does not change the shape of the CSDs, but due to enhanced transport properties, the average grain sizes increase by a factor of 2 and the variance of the distribution increases. Stable d18O and d13C isotope ratios in fluid-affected zones only deviate slightly from the host rock values, suggesting low fluid/rock ratios. Grain growth modelling indicates that the fluid-induced grain size variations can develop within several ka. As inferred from a combination of thermal and grain growth modelling, dykes with widths of up to 1 m have only a restricted influence on grain size deviations smaller than a factor of 1.1.To summarize, considerable grain size variations of up to one order of magnitude can locally result from second-order effects. Such effects require special attention when comparing experimentally derived grain growth kinetics with field studies.