Polymorphism in hybrid male sterility in wild-derived Mus musculus musculus strains on proximal chromosome 17

The hybrid sterility-1 (Hst1) locus at Chr 17 causes male sterility in crosses between the house mouse subspecies Mus musculus domesticus (Mmd) and M. m. musculus (Mmm). This locus has been defined by its polymorphic variants in two laboratory strains (Mmd genome) when mated to PWD/Ph mice (Mmm genome): C57BL/10 (carrying the sterile allele) and C3H (fertile allele). The occurrence of sterile and/or fertile (wild Mmm × C57BL)F₁ males is evidence that polymorphism for this trait also exists in natural populations of Mmm; however, the nature of this polymorphism remains unclear. Therefore, we derived two wild-origin Mmm strains, STUS and STUF, that produce sterile and fertile males, respectively, in crosses with C57BL mice. To determine the genetic basis underlying male fertility, the (STUS × STUF)F₁ females were mated to C57BL/10 J males. About one-third of resulting hybrid males (33.8%) had a significantly smaller epididymis and testes than parental animals and lacked spermatozoa due to meiotic arrest. A further one-fifth of males (20.3%) also had anomalous reproductive traits but produced some spermatozoa. The remaining fertile males (45.9%) displayed no deviation from values found in parental individuals. QTL analysis of the progeny revealed strong associations of male fitness components with the proximal end of Chr 17, and a significant effect of the central section of Chr X on testes mass. The data suggest that genetic incompatibilities associated with male sterility have evolved independently at the proximal end of Chr 17 and are polymorphic within both Mmd and Mmm genomes.     

Coevolution of Cryptosporidium tyzzeri and the house mouse (Mus musculus)

Two house mouse subspecies occur in Europe, eastern and northern Mus musculus musculus (Mmm) and western and southern Mus musculus domesticus (Mmd). A secondary hybrid zone occurs where their ranges meet, running from Scandinavia to the Black Sea. In this paper, we tested a hypothesis that the apicomplexan protozoan species Cryptosporidium tyzzeri has coevolved with the house mouse. More specifically, we assessed to what extent the evolution of this parasite mirrors divergence of the two subspecies. In order to test this hypothesis, we analysed sequence variation at five genes (ssrRNA, Cryptosporidium oocyst wall protein (COWP), thrombospondin-related adhesive protein of Cryptosporidium 1 (TRAP-C1), actin and gp60) in C. tyzzeri isolates from Mmd and Mmm sampled along a transect across the hybrid zone from the Czech Republic to Germany. Mmd samples were supplemented with mice from New Zealand. We found two distinct isolates of C. tyzzeri, each occurring exclusively in one of the mouse subspecies (C. tyzzeri-Mmm and C. tyzzeri-Mmd). In addition to genetic differentiation, oocysts of the C. tyzzeri-Mmd subtype (mean: 4.24 × 3.69 μm) were significantly smaller than oocysts of C. tyzzeri-Mmm (mean: 4.49 × 3.90 μm). Mmm and Mmd were susceptible to experimental infection with both C. tyzzeri subtypes; however, the subtypes were not infective for the rodent species Meriones unguiculatus, Mastomys coucha, Apodemus flavicollis or Cavia porcellus. Overall, our results support the hypothesis that C. tyzzeri is coevolving with Mmm and Mmd.     

Phylogeography of Chinese house mice (Mus musculus musculus/castaneus): distribution,routes of colonization and geographic regions of hybridization

House mice (Mus musculus) are human commensals and have served as a primary model in biomedical, ecological and evolutionary research. Although there is detailed knowledge of the biogeography of house mice in Europe, little is known of the history of house mice in China, despite the fact that China encompasses an enormous portion of their range. In the present study, 535 house mice caught from 29 localities in China were studied by sequencing the mitochondrial D‐loop and genotyping 10 nuclear microsatellite markers distributed on 10 chromosomes. Phylogenetic analyses revealed two evolutionary lineages corresponding to Mus musculus castaneus and Mus musculus musculus in the south and north, respectively, with the Yangtze River approximately representing the boundary. More detailed analyses combining published sequence data from mice sampled in neighbouring countries revealed the migration routes of the two subspecies into China: M. m. castaneus appeared to have migrated through a southern route (Yunnan and Guangxi), whereas M. m. musculus entered China from Kazakhstan through the north‐west border (Xinjiang). Bayesian analysis of mitochondrial sequences indicated rapid population expansions in both subspecies, approximately 4650–9300 and 7150–14 300 years ago for M. m. castaneus and M. m. musculus, respectively. Interestingly, the migration routes of Chinese house mice coincide with the colonization routes of modern humans into China, and the expansion times of house mice are consistent with the development of agriculture in southern and northern China, respectively. Finally, our study confirmed the existence of a hybrid zone between M. m. castaneus and M. m. musculus in China. Further study of this hybrid zone will provide a useful counterpart to the well‐studied hybrid zone between M. m. musculus and Mus musculus domesticus in central Europe.     

The evolution of bat nucleic acid‐sensing Toll‐like receptors

We characterized the nucleic acid‐sensing Toll‐like receptors (TLR) of a New World bat species, the common vampire bat (Desmodus rotundus), and through a comparative molecular evolutionary approach searched for general adaptation patterns among the nucleic acid‐sensing TLRs of eight different bats species belonging to three families (Pteropodidae, Vespertilionidae and Phyllostomidae). We found that the bat TLRs are evolving slowly and mostly under purifying selection and that the divergence pattern of such receptors is overall congruent with the species tree, consistent with the evolution of many other mammalian nuclear genes. However, the chiropteran TLRs exhibited unique mutations fixed in ligand‐binding sites, some of which involved nonconservative amino acid changes and/or targets of positive selection. Such changes could potentially modify protein function and ligand‐binding properties, as some changes were predicted to alter nucleic acid binding motifs in TLR 9. Moreover, evidence for episodic diversifying selection acting specifically upon the bat lineage and sublineages was detected. Thus, the long‐term adaptation of chiropterans to a wide variety of environments and ecological niches with different pathogen profiles is likely to have shaped the evolution of the bat TLRs in an order‐specific manner. The observed evolutionary patterns provide evidence for potential functional differences between bat and other mammalian TLRs in terms of resistance to specific pathogens or recognition of nucleic acids in general.     

Gastrointestinal microbiota of wild and inbred individuals of two house mouse subspecies assessed using high‐throughput parallel pyrosequencing

The effects of gastrointestinal tract microbiota (GTM) on host physiology and health have been the subject of considerable interest in recent years. While a variety of captive bred species have been used in experiments, the extent to which GTM of captive and/or inbred individuals resembles natural composition and variation in wild populations is poorly understood. Using 454 pyrosequencing, we performed 16S rDNA GTM barcoding for 30 wild house mice (Mus musculus) and wild‐derived inbred strain mice belonging to two subspecies (M. m. musculus and M. m. domesticus). Sequenced individuals were selected according to a 2 × 2 experimental design: wild (14) vs. inbred origin (16) and M. m. musculus (15) vs. M. m. domesticus (15). We compared alpha diversity (i.e. number of operational taxonomic units – OTUs), beta diversity (i.e. interindividual variability) and microbiota composition across the four groups. We found no difference between M. m. musculus and M. m. domesticus subspecies, suggesting low effect of genetic differentiation between these two subspecies on GTM structure. Both inbred and wild populations showed the same level of microbial alpha and beta diversity; however, we found strong differentiation in microbiota composition between wild and inbred populations. Relative abundance of ~ 16% of OTUs differed significantly between wild and inbred individuals. As laboratory mice represent the most abundant model for studying the effects of gut microbiota on host metabolism, immunity and neurology, we suggest that the distinctness of laboratory‐kept mouse microbiota, which differs from wild mouse microbiota, needs to be considered in future biomedical research.     

Female House Mice do not Differ in Their Exploratory Behaviour from Males

Dispersal is an important ecological trait affecting genetic variation and dynamics of populations. Hence, the exploratory behaviour prior to actual dispersal may be crucial for potentially dispersing individuals. In mammals, females are traditionally seen as the more philopatric sex and dispersal as male‐biased behaviour, and so behavioural strategies related to the exploration of novel resources should be differentially expressed in males and females. In addition, due to sexual selection exploratory strategies may be expected to vary according to females’ reproductive phase. We employed a standard open‐field test as an approximation of the first phase of dispersal, using adult house mice representing two subspecies, M. m. musculus and M. m. domesticus. We tested the prediction that exploration of neutral area varies in females during different phases of the oestrus cycle and is different between both sexes and subspecies. We expected to find higher exploration in males, as the more dispersing sex and less pronounced subspecies‐specific differences in females than in males. We found no significant effect of the oestrous phase on any of the parameters of the exploratory behaviour measured. Sexual dimorphism was found only in latency to enter the arena in M. m. domesticus where females hesitated longer to enter a new area than males. Significant subspecies‐specific differences were found in three of four tested exploration parameters, so we conclude that females of both subspecies follow similar strategies to those displayed by males. Musculus mice show shorter latency to enter a new area, but once inside, domesticus mice explore the arena significantly longer, with less frequent retreats to a shelter. Our results thus highlight that the role of female dispersal in interdemic gene flow should not be neglected.     

Immunogenetic heterogeneity in a widespread ungulate: the European roe deer (Capreolus capreolus)

Understanding how immune genetic variation is shaped by selective and neutral processes in wild populations is of prime importance in both evolutionary biology and epidemiology. The European roe deer (Capreolus capreolus) has considerably expanded its distribution range these last decades, notably by colonizing agricultural landscapes. This range shift is likely to have led to bottlenecks and increased roe deer exposure to a new range of pathogens that until recently predominantly infected humans and domestic fauna. We therefore investigated the historical and contemporary forces that have shaped variability in a panel of genes involved in innate and acquired immunity in roe deer, including Mhc‐Drb and genes encoding cytokines or toll‐like receptors (TLRs). Together, our results suggest that genetic drift is the main contemporary evolutionary force shaping immunogenetic variation within populations. However, in contrast to the classical view, we found that some innate immune genes involved in micropathogen recognition (e.g. Tlrs) continue to evolve dynamically in roe deer in response to pathogen‐mediated positive selection. Most studied Tlrs (Tlr2, Tlr4 and Tlr5) had similarly high levels of amino acid diversity in the three studied populations including one recently established in southwestern France that showed a clear signature of genetic bottleneck. Tlr2 implicated in the recognition of Gram‐positive bacteria in domestic ungulates, showed strong evidence of balancing selection. The high immunogenetic variation revealed here implies that roe deer are able to cope with a wide spectrum of pathogens and to respond rapidly to emerging infectious diseases.     

ICAM‐1 and MKL‐1 polymorphisms impose considerable impacts on coronary heart disease occurrence

This study was aimed to explore the correlation of intercellular adhesion molecule‐1 (ICAM‐1) K469E and megakaryoblastic leukaemia factor‐1 (MKL‐1) ?184C/T polymorphisms with the susceptibility to coronary heart disease (CHD) in the Chinese Han population. 100 CHD patients and 91 healthy people that had no blood connection with each other were enrolled in this case‐control study. ICAM‐1 and MKL‐1 polymorphisms were genotyped by polymerase chain reaction‐restriction fragment length polymorphism (PCR‐RFLP) approach. Multiple logistic regression was used to analyse the correlation between polymorphisms of ICAM‐1 and MKL‐1 and CHD susceptibility. Differences of genotype and allele frequencies of the two SNPs between case and control groups were analysed by chi‐square test. Odds ratios (ORs) and 95% confidence intervals (CIs) were indicated relative susceptibility of CHD. The distributions of ICAM‐1 and MKL‐1 polymorphisms in each group conformed to Hardy‐Weinberg equilibrium (HWE). After adjusting for traditional risk factors, the TT genotype frequency of MKL‐1 ?184C/T polymorphism was found significantly higher in case group than in control group (P < .05). Meanwhile, T allele frequency increased in case group compared with control group, and the differences had statistical significance (P = .04, OR = 2.34, 95% CI = 1.34‐5.26). Logistic regression analysis in this study proved that smoking, hypertension, diabetes and triglyceride (TG) were all risk factors for CHD ICAM‐1 K469E polymorphism has no association with the onset of CHD. But MKL‐1 ?184C/T polymorphism is associated with the risk of CHD and T allele might be a susceptibility factor for CHD.     

A genome scan of diversifying selection in Ophiocordyceps zombie‐ant fungi suggests a role for enterotoxins in co‐evolution and host specificity

Identification of the genes underlying adaptation sheds light on the biological functions targeted by natural selection. Searches for footprints of positive selection, in the form of rapid amino acid substitutions, and the identification of species‐specific genes have proved to be powerful approaches to identifying the genes involved in host specialization in plant‐pathogenic fungi. We used an evolutionary comparative genomic approach to identify genes underlying host adaptation in the ant‐infecting genus Ophiocordyceps, which manipulates ant behaviour. A comparison of the predicted genes in the genomes of species from three species complexes—O. unilateralis, O. australis and O. subramanianii—revealed an enrichment in pathogenesis‐associated functions, including heat‐labile enterotoxins, among species‐specific genes. Furthermore, these genes were overrepresented among those displaying significant footprints of positive selection. Other categories of genes suspected to be important for virulence and pathogenicity in entomopathogenic fungi (e.g., chitinases, lipases, proteases, core secondary metabolism genes) were much less represented, although a few candidate genes were found to evolve under positive selection. An analysis including orthologs from other entomopathogenic fungi in a broader context showed that positive selection on enterotoxins was specific to the ant‐infecting genus Ophiocordyceps. Together with previous studies reporting the overexpression of an enterotoxin during behavioural manipulation in diseased ants, our findings suggest that heat‐labile enterotoxins are important effectors in host adaptation and co‐evolution in the Ophiocordyceps entomopathogenic fungi.     

Selection on intra‐individual variation in stigma–anther distance in the tropical tree Ipomoea wolcottiana (Convolvulaceae)

• It is well known that animals can exert strong selective pressures on plant traits. However, studies on the evolutionary consequences of plant–animal interactions have mainly focused on understanding how these interactions shape trait means, while overlooking its potential direct effect on the variability among structures within a plant (e.g. flowers and fruits). The degree of within‐plant variability can have strong fitness effects but few studies have evaluated its role as a potential target of selection.
• Here we reanalysed data on Ipomoea wolcottiana stigma–anther distance to test alternate mechanisms driving selection on the mean as well as on intra‐individual variance in 2 years. We found strong negative selection acting on intra‐individual variation but not on mean stigma–anther distance, suggesting independent direct selection on the latter.
• Our result suggests that intra‐individual variance has the potential to be an important target of selection in nature, and that ignoring it could lead to the wrong characterisation of the selection regime.
• We highlight the need for future studies to consider patterns of selection on the mean as well as on intra‐individual variance if we want to understand the full extent of plant–animal interactions as an evolutionary force in nature.

     

Toll‐like receptor variation in the bottlenecked population of the Seychelles warbler: computer simulations see the ‘ghost of selection past’ and quantify the ‘drift debt’

Balancing selection can maintain immunogenetic variation within host populations, but detecting its signal in a postbottlenecked population is challenging due to the potentially overriding effects of drift. Toll‐like receptor genes (TLRs) play a fundamental role in vertebrate immune defence and are predicted to be under balancing selection. We previously characterized variation at TLR loci in the Seychelles warbler (Acrocephalus sechellensis), an endemic passerine that has undergone a historical bottleneck. Five of seven TLR loci were polymorphic, which is in sharp contrast to the low genomewide variation observed. However, standard population genetic statistical methods failed to detect a contemporary signature of selection at any TLR locus. We examined whether the observed TLR polymorphism could be explained by neutral evolution, simulating the population's demography in the software DIYABC. This showed that the posterior distributions of mutation rates had to be unrealistically high to explain the observed genetic variation. We then conducted simulations with an agent‐based model using typical values for the mutation rate, which indicated that weak balancing selection has acted on the three TLR genes. The model was able to detect evidence of past selection elevating TLR polymorphism in the prebottleneck populations, but was unable to discern any effects of balancing selection in the contemporary population. Our results show drift is the overriding evolutionary force that has shaped TLR variation in the contemporary Seychelles warbler population, and the observed TLR polymorphisms might be merely the ‘ghost of selection past’. Forecast models predict immunogenetic variation in this species will continue to be eroded in the absence of contemporary balancing selection. Such ‘drift debt’ occurs when a gene pool has not yet reached its new equilibrium level of polymorphism, and this loss could be an important threat to many recently bottlenecked populations.     

Genetic variation in Tertiary relics: The case of eastern‐Mediterranean Abies (Pinaceae)

The eastern‐Mediterranean Abies taxa, which include both widely distributed species and taxa with minuscule ranges, represent a good model to study the impacts of range size and fragmentation on the levels of genetic diversity and differentiation. To assess the patterns of genetic diversity and phylogenetic relationships among eastern‐Mediterranean Abies taxa, genetic variation was assessed by eight nuclear microsatellite loci in 52 populations of Abies taxa with a focus on those distributed in Turkey and the Caucasus. Both at the population and the taxon level, the subspecies or regional populations of Abies nordmanniana s.l. exhibited generally higher allelic richness, private allelic richness, and expected heterozygosity compared with Abies cilicica s.l. Results of both the Structure analysis and distance‐based approaches showed a strong differentiation of the two A. cilicica subspecies from the rest as well as from each other, whereas the subspecies of A. nordmanniana were distinct but less differentiated. ABC simulations were run for a set of scenarios of phylogeny and past demographic changes. For A. ×olcayana, the simulation gave a poor support for the hypothesis of being a taxon resulting from a past hybridization, the same is true for Abies equi‐trojani: both they represent evolutionary branches of Abies bornmuelleriana.     

Polymorphisms in a desaturase 2 ortholog associate with cuticular hydrocarbon and male mating success variation in a natural population of Drosophila serrata

Elucidating the nature of genetic variation underlying both sexually selected traits and the fitness components of sexual selection is essential to understanding the broader consequences of sexual selection as an evolutionary process. To date, there have been relatively few attempts to connect the genetic variance in sexually selected traits with segregating DNA sequence polymorphisms. We set out to address this in a well‐characterized sexual selection system – the cuticular hydrocarbons (CHCs) of Drosophila serrata – using an indirect association study design that allowed simultaneous estimation of the genetic variance in CHCs, sexual fitness and single nucleotide polymorphism (SNP) effects in an outbred population. We cloned and sequenced an ortholog of the D. melanogaster desaturase 2 gene, previously shown to affect CHC biosynthesis in D. melanogaster, and associated 36 SNPs with minor allele frequencies > 0.02 with variance in CHCs and sexual fitness. Three SNPs had significant multivariate associations with CHC phenotype (q‐value < 0.05). At these loci, minor alleles had multivariate effects on CHCs that were weakly associated with the multivariate direction of sexual selection operating on these traits. Two of these SNPs had pleiotropic associations with male mating success, suggesting these variants may underlie responses to sexual selection due to this locus. There were 15 significant male mating success associations (q‐value < 0.1), and interestingly, we detected a nonrandom pattern in the relationship between allele frequency and direction of effect on male mating success. The minor‐frequency allele usually reduced male mating success, suggesting a positive association between male mating success and total fitness at this locus.     

Do changes in gene expression contribute to sexual isolation and reinforcement in the house mouse?

Expression divergence, rather than sequence divergence, has been shown to be important in speciation, particularly in the early stages of divergence of traits involved in reproductive isolation. In the two European subspecies of house mice, Mus musculus musculus and Mus musculus domesticus, earlier studies have demonstrated olfactory‐based assortative mate preference in populations close to their hybrid zone. It has been suggested that this behaviour evolved following the recent secondary contact between the two taxa (~3,000 years ago) in response to selection against hybridization. To test for a role of changes in gene expression in the observed behavioural shift, we conducted a RNA sequencing experiment on mouse vomeronasal organs. Key candidate genes for pheromone‐based subspecies recognition, the vomeronasal receptors, are expressed in these organs. Overall patterns of gene expression varied significantly between samples from the two subspecies, with a large number of differentially expressed genes between the two taxa. In contrast, only ~200 genes were found repeatedly differentially expressed between populations within M. m. musculus that did or did not display assortative mate preferences (close to or more distant from the hybrid zone, respectively), with an overrepresentation of genes belonging to vomeronasal receptor family 2. These receptors are known to play a key role in recognition of chemical cues that handle information about genetic identity. Interestingly, four of five of these differentially expressed receptors belong to the same phylogenetic cluster, suggesting specialization of a group of closely related receptors in the recognition of odorant signals that may allow subspecies recognition and assortative mating.     

Experimental evidence of genome‐wide impact of ecological selection during early stages of speciation‐with‐gene‐flow

Theory predicts that speciation‐with‐gene‐flow is more likely when the consequences of selection for population divergence transitions from mainly direct effects of selection acting on individual genes to a collective property of all selected genes in the genome. Thus, understanding the direct impacts of ecologically based selection, as well as the indirect effects due to correlations among loci, is critical to understanding speciation. Here, we measure the genome‐wide impacts of host‐associated selection between hawthorn and apple host races of Rhagoletis pomonella (Diptera: Tephritidae), a model for contemporary speciation‐with‐gene‐flow. Allele frequency shifts of 32 455 SNPs induced in a selection experiment based on host phenology were genome wide and highly concordant with genetic divergence between co‐occurring apple and hawthorn flies in nature. This striking genome‐wide similarity between experimental and natural populations of R. pomonella underscores the importance of ecological selection at early stages of divergence and calls for further integration of studies of eco‐evolutionary dynamics and genome divergence.     

Independent evolution of sexual dimorphism and female‐limited mimicry in swallowtail butterflies (Papilio dardanus and Papilio phorcas)

Several species of swallowtail butterflies (genus Papilio) are Batesian mimics that express multiple mimetic female forms, while the males are monomorphic and nonmimetic. The evolution of such sex‐limited mimicry may involve sexual dimorphism arising first and mimicry subsequently. Such a stepwise scenario through a nonmimetic, sexually dimorphic stage has been proposed for two closely related sexually dimorphic species: Papilio phorcas, a nonmimetic species with two female forms, and Papilio dardanus, a female‐limited polymorphic mimetic species. Their close relationship indicates that female‐limited polymorphism could be a shared derived character of the two species. Here, we present a phylogenomic analysis of the dardanus group using 3964 nuclear loci and whole mitochondrial genomes, showing that they are not sister species and thus that the sexually dimorphic state has arisen independently in the two species. Nonhomology of the female polymorphism in both species is supported by population genetic analysis of engrailed, the presumed mimicry switch locus in P. dardanus. McDonald–Kreitman tests performed on SNPs in engrailed showed the signature of balancing selection in a polymorphic population of P. dardanus, but not in monomorphic populations, nor in the nonmimetic P. phorcas. Hence, the wing polymorphism does not balance polymorphisms in engrailed in P. phorcas. Equally, unlike in P. dardanus, none of the SNPs in P. phorcas engrailed were associated with either female morph. We conclude that sexual dimorphism due to female polymorphism evolved independently in both species from monomorphic, nonmimetic states. While sexual selection may drive male–female dimorphism in nonmimetic species, in mimetic Papilios, natural selection for protection from predators in females is an alternative route to sexual dimorphism.     

Mito‐nuclear discordance with evidence of shared ancestral polymorphism and selection in cactophilic species of Drosophila

The Drosophila serido haplogroup is a monophyletic group composed of the following four cryptic and cactophilic species that are endemic to eastern Brazil: D. borborema, D. gouveai, D. seriema and D. serido. Here, we investigate the mito‐nuclear discordance in these species found among the cytochrome c oxidase subunit I (COI) mitochondrial gene, the autosomal alpha‐Esterase‐5 (α‐Est5) and the X‐linked period gene (per). Our analysis indicates that shared polymorphisms in these three molecular markers may be explained by the maintenance of ancestral polymorphisms rather than introgressive hybridization. The primary structures of COI, per and α‐Est5 genes evolve primarily under purifying selection, but we detected some sites that evolved under positive selection in α‐Est5. Considering the high variability of cacti species in eastern Brazil and the role attributed to Drosophila esterases in digestion metabolism and/or the detoxification of several compounds found in cactus tissues, we conjecture about the role of natural selection triggered by host shifts as an important factor in the intraspecific diversification of the D. serido haplogroup.