Evidence of introgression, ecological divergence and adaptation in Asterias sea stars.

Hybrid zones are important windows into the evolutionary dynamics of populations, revealing how processes like introgression and adaptation structure population genomic variation. Importantly, they are useful for understanding speciation and how species respond to their environments. Here, we investigate two closely related sea star species, Asterias rubens and A. forbesi, distributed along rocky European and North American coastlines of the North Atlantic, and use genome-wide molecular markers to infer the distribution of genomic variation within and between species in this group. Using genomic data and environmental niche modelling, we document hybridization occurring between northern New England and the southern Canadian Maritimes. We investigate the factors that maintain this hybrid zone, as well as the environmental variables that putatively drive selection within and between species. We find that the two species differ in their environmental niche breadth; Asterias forbesi displays a relatively narrow environmental niche while conversely, A. rubens has a wider niche breadth. Species distribution models accurately predict hybrids to occur within environmental niche overlap, thereby suggesting environmental selection plays an important role in the maintenance of the hybrid zone. Our results imply that the distribution of genomic variation in North Atlantic sea stars is influenced by the environment, which will be crucial to consider as the climate changes.

Intensity of infection with intracellular Eimeria spp. and pinworms is reduced in hybrid mice compared to parental subspecies.

Genetic diversity in animal immune systems is usually beneficial. In hybrid recombinants, this is less clear, as the immune system could also be impacted by genetic conflicts. In the European house mouse hybrid zone, the long-standing impression that hybrid mice are more highly parasitized and less fit than parentals persists despite the findings of recent studies. Working across a novel transect, we assessed infections by intracellular protozoans, Eimeria spp., and infections by extracellular macroparasites, pinworms. For Eimeria, we found lower intensities in hybrid hosts than in parental mice but no evidence of lowered probability of infection or increased mortality in the centre of the hybrid zone. This means ecological factors are very unlikely to be responsible for the reduced load of infected hybrids. Focusing on parasite intensity (load in infected hosts), we also corroborated reduced pinworm loads reported for hybrid mice in previous studies. We conclude that intensity of diverse parasites, including the previously unstudied Eimeria, is reduced in hybrid mice compared to parental subspecies. We suggest caution in extrapolating this to differences in hybrid host fitness in the absence of, for example, evidence for a link between parasitemia and health.

Bat population recoveries give insight into clustering strategies during hibernation.

BACKGROUND: Behaviour during hibernation contributes to energy conservation in winter. Hibernating bats select roosts with respect to physiological and environmental stressors, available local microclimate and species-specific requirements. RESULTS: We found that, in the period between 1977 and 2018, hibernating Myotis myotis and Rhinolophus hipposideros bats showed exponential population growth. The growth rates, corrected for local winter seasonal severity and winter duration, were equal to 10 and 13%, respectively. While R. hipposideros only utilised the thermally stable and, at survey time, warmer corridors in the hibernaculum, an increasing proportion of M. myotis roosted in the thermally stable corridors as their abundance increased. About 14% of all hibernating M. myotis displayed solitary roosting, irrespective of other covariates. Those bats that clustered together formed progressively larger clusters with increasing abundance, particularly in cold corridors. We found no statistically significant relationship for clustering behaviour or cluster size with winter severity or winter duration. CONCLUSIONS: Abundance of hibernating bats is increasing in Central Europe. As the number of M. myotis bats increases, thermally unstable corridors become saturated with large clusters and the animals begin to roost deeper underground.

When Viruses Don't Go Viral: The Importance of Host Phylogeographic Structure in the Spatial Spread of Arenaviruses.

Many emerging infections are RNA virus spillovers from animal reservoirs. Reservoir identification is necessary for predicting the geographic extent of infection risk, but rarely are taxonomic levels below the animal species considered as reservoir, and only key circumstances in nature and methodology allow intrinsic virus-host associations to be distinguished from simple geographic (co-)isolation. We sampled and genetically characterized in detail a contact zone of two subtaxa of the rodent Mastomys natalensis in Tanzania. We find two distinct arenaviruses, Gairo and Morogoro virus, each spatially confined to a single M. natalensis subtaxon, only co-occurring at the contact zone's centre. Inter-subtaxon hybridization at this centre and a continuum of quality habitat for M. natalensis show that both viruses have the ecological opportunity to spread into the other substaxon's range, but do not, strongly suggesting host-intrinsic barriers. Such barriers could explain why human cases of another M. natalensis-borne arenavirus, Lassa virus, are limited to West Africa.

Holobiont suture zones: Parasite evidence across the European house mouse hybrid zone.

Parasite hybrid zones resulting from host secondary contact have never been described in nature although parasite hybridization is well known and secondary contact should affect them similarly to free-living organisms. When host populations are isolated, diverge and recontact, intimate parasites (host specific, direct life cycle) carried during isolation will also meet and so may form parasite hybrid zones. If so, we hypothesize these should be narrower than the host's hybrid zone as shorter parasite generation time allows potentially higher divergence. We investigate multilocus genetics of two parasites across the European house mouse hybrid zone. We find each host taxon harbours its own parasite taxa. These also hybridize: Parasite hybrid zones are significantly narrower than the host's. Here, we show a host hybrid zone is a suture zone for a subset of its parasite community and highlight the potential of such systems as windows on the evolutionary processes of host-parasite interactions and recombinant pathogen emergence.

Murine cytomegalovirus is not restricted to the house mouse Mus musculus domesticus: prevalence and genetic diversity in the European house mouse hybrid zone.

UNLABELLED: Murine cytomegalovirus (MCMV) is a betaherpesvirus of the house mouse, Mus musculus domesticus. It is a common infectious agent of wild mice and a highly studied pathogen of the laboratory mouse. Betaherpesviruses are specific to their hosts, and it is not known if other Mus taxa carry MCMV or if it is restricted to M. m. domesticus. We sampled mice over a 145-km transect of Bavaria-Bohemia crossing a hybrid zone between M. m. domesticus and Mus musculus musculus in order to investigate the occurrence of MCMV in two Mus subspecies and to test the limits of the specificity of the virus for its host. We hypothesized that if the two subspecies carry MCMV and if the virus is highly specific to its host, divergent MCMV lineages would have codiverged with their hosts and would have a geographical distribution constrained by the host genetic background. A total of 520 mice were tested by enzyme-linked immunosorbent assay (ELISA) and/or nested PCR targeting the M94 gene. Seropositive and PCR-positive individuals were found in both Mus subspecies. Seroprevalence was high, at 79.4%, but viral DNA was detected in only 41.7% of mice. Sequencing revealed 20 haplotypes clustering in 3 clades that match the host genetic structure in the hybrid zone, showing 1 and 2 MCMV lineages in M. m. domesticus and M. m. musculus, respectively. The estimated time to the most recent common ancestor (1.1 million years ago [Mya]) of the MCMVs matches that of their hosts. In conclusion, MCMV has coevolved with these hosts, suggesting that its diversity in nature may be underappreciated, since other members of the subgenus Mus likely carry different MCMVs. IMPORTANCE: Murine cytomegalovirus (MCMV) is a betaherpesvirus of the house mouse, Mus musculus domesticus, an important lab model for human cytomegalovirus (HCMV) infection. The majority of lab studies are based on only two strains of MCMVs isolated from M. m. domesticus, Smith and K181, the latter derived from repeated passage of Smith in mouse submaxillary glands. The presence of MCMV in other members of the Mus subgenus had not even been investigated. By screening mouse samples collected in the European house mouse hybrid zone between M. m. domesticus and M. m. musculus, we show that MCMV is not restricted to the M. m. domesticus subspecies and that MCMVs likely codiverged with their Mus hosts. Thus, the diversity of MCMV in nature may be seriously underappreciated, since other members of the subgenus Mus likely carry their own MCMV lineages.

Steep clines within a highly permeable genome across a hybrid zone between two subspecies of the European rabbit.

Maintenance of genetic distinction in the face of gene flow is an important aspect of the speciation process. Here, we provide a detailed spatial and genetic characterization of a hybrid zone between two subspecies of the European rabbit. We examined patterns of allele frequency change for 22 markers located on the autosomes, X-chromosome, Y-chromosome and mtDNA in 1078 individuals sampled across the hybrid zone. While some loci revealed extremely wide clines (w ≥ 300 km) relative to an estimated dispersal of 1.95-4.22 km/generation, others showed abrupt transitions (w ≈ 10 km), indicating localized genomic regions of strong selection against introgression. The subset of loci showing steep clines had largely coincident centers and stepped changes in allele frequency that did not co-localize with any physical barrier or ecotone, suggesting that the rabbit hybrid zone is a tension zone. The steepest clines were for X- and Y-chromosome markers. Our results are consistent with previous inference based on DNA sequence variation of individuals sampled in allopatry in suggesting that a large proportion of each genome has escaped the overall barrier to gene flow in the middle of the hybrid zone. These results imply an old history of hybridization and high effective gene flow and anticipate that isolation factors should often localize to small genomic regions.

A gene copy number arms race in action: X,Y-chromosome transmission distortion across a species barrier.

A remarkable gene copy number (CN) arms race system has recently been described in laboratory mice, where Slx;Slxl1 and Sly genes compete over transmission by altering the fertilization success of X and Y chromosome-bearing sperm, respectively. Here, we focus on this system in nature, where natural selection can counter CN/gene product escalation. Our model is house mouse subspecies hybridizing in Europe. In some regions, Y chromosomes of the Eastern subspecies have introgressed onto Western genomic backgrounds, accompanied by sex ratio distortion in favor of males, consistent with the inbred lines suggested mechanism: Overabundance of SLY protein expressed by invading Y chromosomes. We take Slx as representative of the X side of this arms race and measure Slx|Sly CN and expression across an "Invasion" transect where Ys introgress and a "Control" transect with negligible introgression. Since we found similar Slx|Sly ratios in both transects, SLY overabundance is unlikely to explain the introgression. However, Slx CN is relatively low in the introgression area, suggesting that Slx is less able to combat Sly effects here. Furthermore, deterministic changes in Slx;Sly expression proportions versus CN proportions suggest standing variation for trans regulation of Slx|Sly is being co-opted in nature where their arms race reduces population fitness.

Subspecific rodent taxa as the relevant host taxonomic level for mammarenavirus host specificity.

Mastomys natalensis-borne mammarenaviruses appear specific to subspecific M. natalensis taxa rather than to the whole species. Yet mammarenaviruses carried by M. natalensis are known to spill over and jump hosts in northern sub-Saharan Africa. Phylogeographic studies increasingly show that, like M. natalensis, small mammals in sub-Saharan Africa are often genetically structured into several subspecific taxa. Other mammarenaviruses may thus also form virus-subspecific host taxon associations. To investigate this, and if mammarenaviruses carried by M. natalensis in southern Africa are less prone to spill-over, we screened 1225 non-M. natalensis samples from Tanzania where many small mammal taxa meet. We found mammarenavirus RNA in 6 samples. Genetic/genomic characterisation confirmed they were not spill-over from M. natalensis. We detected host jumps among rodent tribe members and an association between mammarenaviruses and subspecific taxa of Mus minutoides and Grammomys surdaster, indicating host genetic structure may be crucial to understand virus distribution and host specificity.

Sperm-related phenotypes implicated in both maintenance and breakdown of a natural species barrier in the house mouse.

The house mouse hybrid zone (HMHZ) is a species barrier thought to be maintained by a balance between dispersal and natural selection against hybrids. While the HMHZ is characterized by frequency discontinuities for some sex chromosome markers, there is an unexpected large-scale regional introgression of a Y chromosome across the barrier, in defiance of Haldane's rule. Recent work suggests that a major force maintaining the species barrier acts through sperm traits. Here, we test whether the Y chromosome penetration of the species barrier acts through sperm traits by assessing sperm characteristics of wild-caught males directly in a field laboratory set up in a Y introgression region of the HMHZ, later calculating the hybrid index of each male using 1401 diagnostic single nucleotide polymorphisms (SNPs). We found that both sperm count (SC) and sperm velocity were significantly reduced across the natural spectrum of hybrids. However, SC was more than rescued in the presence of the invading Y. Our results imply an asymmetric advantage for Y chromosome introgression consistent with the observed large-scale introgression. We suggest that selection on sperm-related traits probably explains a large component of patterns observed in the natural hybrid zone, including the Y chromosome penetration.

Genomic evidence of paternal genome elimination in the globular springtail Allacma fusca.

Paternal genome elimination-a type of reproduction in which males inherit but fail to pass on their father's genome-evolved independently in 6-8 arthropod clades. Thousands of species, including several important for agriculture, reproduce via this mode of reproduction. While paternal genome elimination is well established in some of the clades, the evidence in globular springtails (Symphypleona) remains elusive, even though they represent the oldest and most species-rich clade putatively reproducing via paternal genome elimination. We sequenced genomic DNA from whole bodies of Allacma fusca males with high fractions (>27.5%) of sperm to conclusively confirm that all the sperm carry 1 parental haplotype only. Although it is suggestive that the single haplotype present in sperm is maternally inherited, definitive genetic proof of the parent of origin is still needed. The genomic approach we developed allows for the detection of genotypic differences between germline and soma in all species with sufficiently high fraction of germline in their bodies. This opens new opportunities for scans of reproductive modes in small organisms.

Asymmetric reproductive isolation between terminal forms of the salamander ring species Ensatina eschscholtzii revealed by fine-scale genetic analysis of a hybrid zone.

BACKGROUND: Ring species, exemplified by salamanders of the Ensatina eschscholtzii complex, represent a special window into the speciation process because they allow the history of species formation to be traced back in time through the geographically differentiated forms connecting the two terminal forms of the ring. Of particular interest is the nature and extent of reproductive isolation between the geographically terminal forms, in this case E. e. eschscholtzii and E. e. klauberi. Previous studies have documented infrequent hybridization at the end of the ring. Here, we report the first fine-scale genetic analysis of a hybrid zone between the terminal forms in southern California using individual-based Bayesian analyses of multilocus genetic data to estimate levels and direction of hybridization and maximum-likelihood analysis of linkage disequilibrium and cline shape to make inferences about migration and selection in the hybrid zone. RESULTS: The center of the hybrid zone has a high proportion of hybrids, about half of which were classified as F1s. Clines are narrow with respect to dispersal, and there are significant deviations from Hardy-Weinberg equilibrium as well as nonrandom associations (linkage disequilibria) between alleles characteristic of each parental type. There is cytonuclear discordance, both in terms of introgression and the geographic position of mitochondrial versus nuclear clines. Genetic disequilibrium is concentrated on the eschscholtzii side of the zone. Nearly all hybrids possess klauberi mtDNA, indicating that most hybrids are formed from female klauberi mating with male eschscholtzii or male hybrids (but not vice versa). CONCLUSIONS: Our results are consistent with a tension zone trapped at an ecotone, with gene combinations characteristic of klauberi showing up on the eschscholtzii side of the zone due to asymmetric hybridization. We suggest that the observed asymmetry is best explained by increased discriminatory power of eschscholtzii females, or asymmetric postzygotic isolation. The relatively high frequency of hybrids, particularly F1s, contrasts with other contacts between the terminal forms, and with other contacts between other divergent Ensatina lineages, highlighting the diverse outcomes of secondary contact within a single species complex.

Reinforcement selection acting on the European house mouse hybrid zone.

Behavioural isolation may lead to complete speciation when partial postzygotic isolation acts in the presence of divergent-specific mate-recognition systems. These conditions exist where Mus musculus musculus and M. m. domesticus come into contact and hybridize. We studied two mate-recognition signal systems, based on urinary and salivary proteins, across a Central European portion of the mouse hybrid zone. Introgression of the genomic regions responsible for these signals: the major urinary proteins (MUPs) and androgen binding proteins (ABPs), respectively, was compared to introgression at loci assumed to be nearly neutral and those under selection against hybridization. The preference of individuals taken from across the zone regarding these signals was measured in Y mazes, and we develop a model for the analysis of the transition of such traits under reinforcement selection. The strongest assortative preferences were found in males for urine and females for ABP. Clinal analyses confirm nearly neutral introgression of an Abp locus and two loci closely linked to the Abp gene cluster, whereas two markers flanking the Mup gene region reveal unexpected introgression. Geographic change in the preference traits matches our reinforcement selection model significantly better than standard cline models. Our study confirms that behavioural barriers are important components of reproductive isolation between the house mouse subspecies.

A dense linkage map for a large repetitive genome: discovery of the sex-determining region in hybridizing fire-bellied toads (Bombina bombina and Bombina variegata).

Genomic analysis of hybrid zones offers unique insights into emerging reproductive isolation and the dynamics of introgression. Because hybrid genomes consist of blocks inherited from one or the other parental taxon, linkage information is essential. In most cases, the spectrum of local ancestry tracts can be efficiently uncovered from dense linkage maps. Here, we report the development of such a map for the hybridizing toads, Bombina bombina and Bombina variegata (Anura: Bombinatoridae). Faced with the challenge of a large (7-10 Gb), repetitive genome, we set out to identify a large number of Mendelian markers in the nonrepetitive portion of the genome that report B. bombina vs B. variegata ancestry with appropriately quantified statistical support. Bait sequences for targeted enrichment were selected from a draft genome assembly, after filtering highly repetitive sequences. We developed a novel approach to infer the most likely diplotype per sample and locus from the raw read mapping data, which is robust to over-merging and obviates arbitrary filtering thresholds. Validation of the resulting map with 4755 markers underscored the large-scale synteny between Bombina and Xenopus tropicalis. By assessing the sex of late-stage F2 tadpoles from histological sections, we identified the sex-determining region in the Bombina genome to 7 cM on LG5, which is homologous to X. tropicalis chromosome 5, and inferred male heterogamety. Interestingly, chromosome 5 has been repeatedly recruited as a sex chromosome in anurans with XY sex determination.

Sympatric occurrence of 3 arenaviruses, Tanzania.

To determine the specificity of Morogoro virus for its reservoir host, we studied its host range and genetic diversity in Tanzania. We found that 2 rodent species other than Mastomys natalensis mice carry arenaviruses. Analysis of 340 nt of the viral RNA polymerase gene showed sympatric occurrence of 3 distinct arenaviruses.

Genetic and historic evidence for climate-driven population fragmentation in a top cetacean predator: the harbour porpoises in European water.

Recent climate change has triggered profound reorganization in northeast Atlantic ecosystems, with substantial impact on the distribution of marine assemblages from plankton to fishes. However, assessing the repercussions on apex marine predators remains a challenging issue, especially for pelagic species. In this study, we use Bayesian coalescent modelling of microsatellite variation to track the population demographic history of one of the smallest temperate cetaceans, the harbour porpoise (Phocoena phocoena) in European waters. Combining genetic inferences with palaeo-oceanographic and historical records provides strong evidence that populations of harbour porpoises have responded markedly to the recent climate-driven reorganization in the eastern North Atlantic food web. This response includes the isolation of porpoises in Iberian waters from those further north only approximately 300 years ago with a predominant northward migration, contemporaneous with the warming trend underway since the 'Little Ice Age' period and with the ongoing retreat of cold-water fishes from the Bay of Biscay. The extinction or exodus of harbour porpoises from the Mediterranean Sea (leaving an isolated relict population in the Black Sea) has lacked a coherent explanation. The present results suggest that the fragmentation of harbour distribution range in the Mediterranean Sea was triggered during the warm 'Mid-Holocene Optimum' period (approx. 5000 years ago), by the end of the post-glacial nutrient-rich 'Sapropel' conditions that prevailed before that time.

The impact of global selection on local adaptation and reproductive isolation.

Despite the homogenizing effect of strong gene flow between two populations, adaptation under symmetric divergent selection pressures results in partial reproductive isolation: adaptive substitutions act as local barriers to gene flow, and if divergent selection continues unimpeded, this will result in complete reproductive isolation of the two populations, i.e. speciation. However, a key issue in framing the process of speciation as a tension between local adaptation and the homogenizing force of gene flow is that the mutation process is blind to changes in the environment and therefore tends to limit adaptation. Here we investigate how globally beneficial mutations (GBMs) affect divergent local adaptation and reproductive isolation. When phenotypic divergence is finite, we show that the presence of GBMs limits local adaptation, generating a persistent genetic load at the loci that contribute to the trait under divergent selection and reducing genome-wide divergence. Furthermore, we show that while GBMs cannot prohibit the process of continuous differentiation, they induce a substantial delay in the genome-wide shutdown of gene flow. This article is part of the theme issue 'Towards the completion of speciation: the evolution of reproductive isolation beyond the first barriers'.

Three arenaviruses in three subspecific natal multimammate mouse taxa in Tanzania: same host specificity, but different spatial genetic structure?

Mastomys natalensis is widespread in sub-Saharan Africa and hosts several arenavirus species, including the pathogenic zoonotic Lassa virus in West Africa. Mitochondrial lineages sub-divide the range of M. natalensis and have been associated with cryptic structure within the species. To test specificity of arenaviruses to hosts carrying these lineages, we screened 1772 M. natalensis in a large area of Tanzania where three mitochondrial lineages meet. We detected fifty-two individuals that were positive for one of three arenaviruses: Gairo, Morogoro, and Luna virus. This is the first record of Luna virus in Tanzania. We confirmed the specificity of each arenavirus to a distinct host mitochondrial lineage except for three cases in one locality at the centre of a host hybrid zone. No arenaviruses were detected in a large part of the study area. Morogoro and Gairo virus showed differences in prevalence (Morogoro virus lower than Gairo virus) and in genetic structure (Morogoro virus more structured than Gairo virus). However, both viruses have genetic neighbourhood size estimates of the same order of magnitude as Lassa virus. While differences in arenavirus and/or host evolutionary and ecological dynamics may exist, Tanzanian arenaviruses could be suited to model Lassa virus dynamics in M. natalensis.

Evolutionary history of an MHC gene in two leporid species: characterisation of Mhc-DQA in the European brown hare and comparison with the European rabbit.

We surveyed the genetic diversity of the expressed major histocompatibility complex class II DQA locus in natural populations of European brown hares, Lepus europaeus, from Austria and Belgium (267 individuals in total). Based on cDNA sequences, we designed hare-specific primers to amplify the highly variable second exon of the DQA gene. Using cloning-sequencing methodology and capillary electrophoresis single-strand conformation polymorphism, we found ten alleles of the DQA exon 2 locus across these two European regions, of which eight are described for the first time. To search for signals of selection and recombination in the evolution of the DQA gene within the leporids, we augmented our sample with orthologous DQA alleles from the European rabbit, Oryctolagus cuniculus, in order to carry out a species level, species pairwise comparison. We found evidence of recombination in the history of the DQA sequences in leporids with some recombinant alleles bridging the species divide. In both species, selection on peptide binding site codons can be detected, though stronger for the rabbit. This result suggests that there may be a differential selection pressure in the deeper evolutionary history of these two species due to differences in several demographic and ecological traits likely subjecting them to differential selection by parasites. Finally, evolutionary relationships show a widespread and statistically significant intermingling of alleles from the two species. The many macroparasites shared between hares and rabbits may explain this pattern of trans-species polymorphism.