Inferring the history of speciation in house mice from autosomal,X‐linked,Y‐linked and mitochondrial genes

Patterns of genetic differentiation among taxa at early stages of divergence provide an opportunity to make inferences about the history of speciation. Here, we conduct a survey of DNA‐sequence polymorphism and divergence at loci on the autosomes, X chromosome, Y chromosome and mitochondrial DNA in samples of Mus domesticus, M. musculus and M. castaneus. We analyzed our data under a divergence with gene flow model and estimate that the effective population size of M. castaneus is 200 000–400 000, of M. domesticus is 100 000–200 000 and of M. musculus is 60 000–120 000. These data also suggest that these species started to diverge approximately 500 000 years ago. Consistent with this recent divergence, we observed considerable variation in the genealogical patterns among loci. For some loci, all alleles within each species formed a monophyletic group, while at other loci, species were intermingled on the phylogeny of alleles. This intermingling probably reflects both incomplete lineage sorting and gene flow after divergence. Likelihood ratio tests rejected a strict allopatric model with no gene flow in comparisons between each pair of species. Gene flow was asymmetric: no gene flow was detected into M. domesticus, while significant gene flow was detected into both M. castaneus and M. musculus. Finally, most of the gene flow occurred at autosomal loci, resulting in a significantly higher ratio of fixed differences to polymorphisms at the X and Y chromosomes relative to autosomes in some comparisons, or just the X chromosome in others, emphasizing the important role of the sex chromosomes in general and the X chromosome in particular in speciation.     

Population genetics in nonmodel organisms: II. natural selection in marginal habitats revealed by deep sequencing on dual platforms

Population genetics of species living in marginal habitats could be particularly informative about the genetics of adaptation, but such analyses have not been readily feasible until recently. Sonneratia alba, a mangrove species widely distributed in the Indo-West Pacific, provides a very suitable system for the study of local adaptation. In this study, we analyzed DNA variation by pooling 71 genes from 85-100 individuals for DNA sequencing. For each of the two nearby S. alba populations, we obtained ~2,500 × coverage on the Illumina GA platform and for the Sanya population, an additional 5,400 × coverage on the AB SOLiD platform. For the Sanya sample, although each sequencing method called many putative single nucleotide polymorphisms, the two sets of calls did not overlap, suggesting platform-dependent errors. Conventional sequencing corroborated that each population is monomorphic. The two populations differ by 54 bp of 79,000 sites, but 90% of the variants are found in 10% of the genes. Strong local adaptation and high migration may help to explain the extensive monomorphism shared by the two populations in the presence of a small number of highly differentiated loci.     

Ancestral polymorphism in exon 2 of bluethroat (Luscinia svecica) MHC class II B genes

The genes of the major histocompatibility complex (MHC) are important model genes for understanding selective forces in evolution. Here, we document, using a cloning and sequencing approach, high polymorphism at the exon 2 of the MHC class II B (MHCIIB) genes in the bluethroat (Luscinia svecica); a minimum of 61 unique alleles were detected in 20 individuals, and at least 11 functional loci. In addition, several pseudogenes were revealed. The specimens originated from three different bluethroat subspecies (azuricollis, cyanecula and svecica), and we also analysed four specimens of the closely related thrush nightingale (L. luscinia) for comparison. Phylogenetic analyses of the functional alleles revealed 258 equally parsimonious trees with poor statistical support for the majority of nodes. The distribution of the sequences in the trees point to an ancestral origin of the polymorphism in MHC class II B genes, a portion of which predated the phylogenetic split between the bluethroat and the thrush nightingale. Strong signatures of balancing selection were uncovered for the codons coding for the peptide‐binding residues of the functional MHCIIB exon 2 alleles. Our results highlight the importance of duplication and recombination events for shaping passerine MHC and give insights in the evolutionary dynamics of MHC variation among closely related taxa.     

The population genetics of the apple maggot fly,Rhagoletis pomonella and the snowberry maggot,R. zephyria: implications for models of sympatric speciation

True fruit flies belonging to the Rhagoletis pomonella (Diptera:Tephritidae) sibling species complex possess several attributes consistent with a history of sympatric divergence via host plant shifts. Here, we investigate whether hybridization and genetic introgression is occurring between two members of the group, Rhagoletis pomonella (Walsh), whose primary hosts are domestic apples (Malus pumila) and hawthorns (Crataegus spp., and R. zephyria (Snow) whose host is snowberries (Symphoricarpos spp.). These flies are important because they appear to be at a transition stage between taxa reproductively isolated solely on the basis of host plant-related adaptations and those that have evolved additional non-host dependent sterility and inviability. Observing extensive hybridization and introgression between R. pomonella and R. zephyria would have major repercussions for current models of sympatric speciation.In a survey of allozyme and mitochondrial DNA (mtDNA) variation for 1105 flies collected throughout the northern United States we found two results suggesting that low level hybridization occurs between R. pomonella and R. zephyria. (1) Two flies reared from snowberries and one fly reared from hawthoms had genotypes indicative of them being R. pomonella and R. zephyria, respectively. Rhagoletis pomonella and R. zephyria adults may therefore occasionally frequent each others host plant, providing the opportunity for hybridization. (2) Four flies collected from hawthorns and one from snowberries had genotypes that made them likely to be F1 hybrids.Likelihood analysis revealed the data were also consistent with an hypothesis of shared ancestral alleles (i.e., the pattern of genetic variation could also be explained by R. pomonella and R. zephyria sharing alleles/haplotypes whose origins date to a common ancestor). We estimated that, in the absence of interspecific mating, random assortment of genes within R. pomonella and R. zephyria populations would produce an average of 5.4 flies with genotypes suggesting they were F1 hybrids – a number equivalent to the 5 putative F1 hybrids observed in the study. Our results therefore underscore the difficulty in distinguishing between hypotheses of low level introgression and shared ancestral polymorphism. But even if hybridization is occurring, the data suggest that it is happening at a very low and probably evolutionarily insignificant level (perhaps 0.09% per generation), consistent with sympatric speciation theory. Future tests are discussed that could help resolve the hybridization issue for R. pomonella and R. zephyria.     

A young clade repeating an old pattern: diversity in Nothonotus darters (Teleostei: Percidae) endemic to the Cumberland River

Hypotheses of diversification in eastern North American freshwater fishes have focused primarily on allopatric distributions of species between disjunct highland areas and major river systems. However, these hypotheses do not fully explain the rich diversity of species within highland regions and river systems. Relatively old diversification events at small geographic scales have been observed in the Barcheek Darter subclade that occurs in the Cumberland River drainage (CRD) in Kentucky and Tennessee, United States of America, but it is unknown if this pattern is consistent in other darter subclades. We explored phylogeographic diversity in two species of Nothonotus darters, N. microlepidus and N. sanguifluus, endemic to the CRD to compare phylogenetic patterns between Barcheek Darters and species of Nothonotus. We collected sequence data for a mitochondrial gene (cytb) and three nuclear genes (MLL, S7 and RAG1) from 19 N. microlepidus and 35 N. sanguifluus specimens. Gene trees were estimated using maximum likelihood and Bayesian methods, and a 'species tree' was inferred using a Bayesian method. These trees indicate that species diversity in Nothonotus is underestimated. Five distinct lineages were evident, despite retained ancestral polymorphism and unsampled extirpated populations. Comparison of chronograms for Barcheek Darters and Nothonotus revealed that microendemism resulting from species diversification at small geographic scales in the CRD is a consistent pattern in both old and young darter subclades. Our analyses reveal that geographic isolating mechanisms that result in similar phylogeographic patterns in the CRD are persistent through long expanses of evolutionary time.     

Population genetics of glyoxalase I (GLO) polymorphism in India

Phenotype and gene frequency data are presented on the glyoxalase I (GLO) polymorphism in seven endogamous caste groups: Jat Sikh, Ramdasia Sikh, Ramgarhia Sikh, Khatri, Brahmin and Bania of Patiala district, and Jat Sikh of Faridkot district of Punjab, North-West India. Apparently, there is considerable heterogeneity in the frequency distribution of the GLO1 gene that varies from 0.168 in Bania to 0.287 in Brahmin. However, these differences are not statistically significant, and the overall GLO1 frequency in Punjab is well within the North Indian range.     

Meeting review: American Genetics Association Symposium on the genetics of speciation

Yeast can be engineered to carry human chromosomes; highly diverged ducks can produce viable, fertile offspring; and mitochondrial genes can move between widely divergent groups of plants. Some sunflower or oak species have porous genomes; mice, crickets, birds, and butterflies form hybrid zones; and bacterial lineages have been exchanging genes for several billion years. Even so, nature is discrete and full of species. Here, we discuss some of the ingredients that make nature discrete and can lead to clustering even in the presence of gene flow. Many of these results have been recently published, in this issue and elsewhere, and were discussed at the Genetics of Speciation Symposium held at the annual meeting of the American Genetics Association, Vancouver, Canada, in 2006.     

A gene's eye view of epistasis,selection and speciation

In this mini‐review, I discuss the effects of gene interaction or epistasis from a `gene's eye view.' By a `gene's eye view' of epistasis, I mean that I will consider a single, bi‐allelic locus, A , whose effects on fitness result only from its interactions with alleles of another, unknown locus, X . I will show how changes in the frequencies of alleles at the background locus affect the relationship of alleles at the A ‐locus to fitness. Changing the genetic background changes the fundamental characteristics of the A ‐locus, such as the magnitude and sign of allelic effects on fitness, and, consequently, it changes the strength and pattern of selection. I consider each of the four kinds of pair–wise interactions between two loci and show that some kinds of epistasis are more sensitive than others to population genetic subdivision. Lastly, I show that some kinds of epistasis are more likely than others to affect the process of speciation and contribute to or be responsible for general genetic features of interspecific hybrids, such as Haldane's rule.     

Population structure and speciation in the genus Tursiops based on microsatellite and mitochondrial DNA analyses

Bottlenose dolphins (Tursiops truncatus) have a world-wide distribution, and show morphotypic variation among regions. Distinctions between coastal and pelagic populations have been documented; however, regional patterns of differentiation had not been previously investigated in a wider geographic context. We analysed up to nine different populations from seven different areas of the world by mitochondrial DNA and microsatellite DNA markers, and found differentiation among all putative regional populations. Both mtDNA and microsatellite DNA data show significant differentiation, suggesting restricted gene flow for both males and females. Dolphins in coastal habitat showed less variability and were in most cases differentiated from a pelagic lineage, which could suggest local founder events in some cases. Two coastal populations recently classified as belonging to a new species, T. aduncus, were each highly differentiated from populations of the truncatus morphotype, and from each other, suggesting a possible third species represented by the South African aduncus type.     

Gradual speciation in a global hotspot of plant diversity

The speciation process that underlies recent, rapid radiations of plants is controversial, and suggested mechanisms range from pollinator or ecological niche differentiation to allopatry and nonadaptive divergence. Phylogenetic approaches to locating the most appropriate speciation models have been constrained by the low levels of molecular divergence between recently diverged species, which are typical of recent, rapid radiations. In this issue of Molecular Ecology, Rymer et al. (2010) used coalescence analyses of sequence data and genome scans of Amplified Fragment Length Polymorphism (AFLP) loci to demonstrate that in a species complex in the irid genus Gladiolus, a member of the hyper diverse Cape flora of southern Africa, speciation is a gradual process. Older divergences are genetically more differentiated, and show a greater difference in flowering time and floral morphology, than taxa that diverged more recently. There is no evidence of any abrupt events. Gene flow is limited by shifts in flowering time and floral morphology; thus, by pre-zygotic rather than by post-zygotic mechanisms, these evolved together with the occupation of somewhat different habitats. This research gives the first critical insight into how the remarkable diversity in a diversity hotspot could have arisen. More importantly, it demonstrates that the speciation process in recent, rapid radiations is tractable and can be investigated with suitable genetic tools.     

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.     

Microsatellite polymorphism and population subdivision in natural populations of European sea bass Dicentrarchus labrax (Linnaeus, 1758)

Polymorphism of microsatellite markers was used to study the genetic variability and structure in natural populations of European sea bass Dicentrarchus labrax. The data consisted of six microsatellite loci analysed for 172 individuals from three samples collected in the Golfe-du-Lion (France) and one sample collected in the Golfo-de-Valencia (Spain). Our goals were (i) to assess the level of genetic variability as revealed by these markers, (ii) to estimate the genetic differentiation among natural populations within a restricted area, and (iii) to evaluate how microsatellite loci fit the predictions of the two most widely used mutation models (the infinite allele model and the stepwise mutation model). As expected, our results indicate that the genetic polymorphism is very high when compared with previously used genetic markers, the mean expected heterozygosity per locus ranging between 0.69 and 0.93. We also found that all loci but one fitted the infinite allele model better. Using this model as a lower limit, we could extrapolate from the observed diversity effective population sizes on the order of 35 000 individuals. Our results also suggest that there may be a slight genetic differentiation between the two gulfs (F_ST= 0.007, P < 0.05), indicating that the corresponding populations are likely to be dynamically independent. This finding for a species with high dispersal abilities, if confirmed, has important beatings on fish-stock assessment.     

Meiosis and speciation: A study in a speciatingMus terricolor complex

The three chromosomal species of theMus terricolor complex possess 2n = 40 chromosomes. We show that their karyotypes differ in stable heterochromatin variations fixed in homozygous condition as prominent short arms in autosomes 1, 3 and 6. The three chromosomal species exhibit a high incidence of polymorphisms for Robertsonian fusions and pericentric inversions. Breeding experiments and histological analysis of testis show that heterozygosity for pericentric inversions and Robertsonian fusions had no effect on fertility. Meiotic analysis shows normal overall progression of meiosis in the heterozygotes, which is consistent with their normal gametogenesis. Nevertheless, both the inversion and fusion heterozygotes had undergone some alterations in the regular process of homologous synapsis, and it appeared that certain features of the meiotic system circumvented the potential negative effects of these polymorphic chromosomal rearrangements. The results indicate that the attributes of the meiotic system in a given organism could modulate the potential of a chromosomal rearrangement as reproductive barrier. The meiotic modulation hypothesis offers an explanation for the contradictory effects of the similar kinds of chromosomal mutations reported in different species.     

Ancient lakes revisited: from the ecology to the genetics of speciation

Explosive speciation in ancient lakes has fascinated biologists for centuries and has inspired classical work on the tempo and modes of speciation. Considerable attention has been directed towards the extrinsic forces of speciation—the geological, geographical and ecological peculiarities of ancient lakes. Recently, there has been a resurgence of interest in the intrinsic nature of these radiations, the biological characteristics conducive to speciation. While new species are thought to arise mainly by the gradual enhancement of reproductive isolation among geographically isolated populations, ancient lakes provide little evidence for a predominant role of geography in speciation. Recent phylogenetic work provides strong evidence that multiple colonization waves were followed by parallel intralacustrine radiations that proceeded at relatively rapid rates despite long‐term gene flow through hybridization and introgression. Several studies suggest that hybridization itself might act as a key evolutionary mechanism by triggering major genomic reorganization/revolution and enabling the colonization of new ecological niches in ancient lakes. These studies propose that hybridization is not only of little impediment to diversification but could act as an important force in facilitating habitat transitions, promoting postcolonization adaptations and accelerating diversification. Emerging ecological genomic approaches are beginning to shed light on the long‐standing evolutionary dilemma of speciation in the face of gene flow. We propose an integrative programme for future studies on speciation in ancient lakes.     

Micro-evolution in the weevil genus Larinus: the formation of host biotypes and speciation

Data are presented on allozyme variation between 15 populations of the stenophagous capitulum weevil, Larinus cynarae , and three populations of its congener, L. latus , that had been collected throughout the northern mediterranean range of these species. A phenetic analysis of these data revealed no direct relationship between genetic variation and host-plant association within L. cynarae , but there was a strong geographical structuring of allozyme patterns. Most of the genetic variation was due to differences between geographical regions and variation within these was small. Wright's F _ST values showed that Italian and Greek populations of L. cynarae were most distinct from L. latus , with southern Iberian, northern Spanish and French populations increasingly less so. This pattern was associated with a cline in the frequencies of certain alleles along this geographical arc from France to Greece. A phenogram of Nei's genetic distances indicated the close genetic relationship between the two species of Larinus and separated the populations of L. cynarae into three allopatric groups. These groups have different host-plant spectra - dominated by Cynara cardunculus in Italy and Greece, Cynara humilis/Onopordum in southern Iberia and Onopordum spp. in France/Northern Spain - and can be considered to be host biotypes of L. cynarae. L. latus , which occurs in Greece and further east is also an Onopordum specialist. An analysis of the phylogeny of this group of Larinus indicates a primary separation into eastern ( L. latus ) and western ( L. cynarae ) taxa, with further branching of the L. cynarae lineage into the putative host-biotypes. An hypothesis for the evolution of these taxa is given, based on the evolutionary history of host-plant taxa and geographical constraints.     

Evidence of microsatellite hitch-hiking selection in Atlantic cod (Gadus morhua L.): implications for inferring population structure in nonmodel organisms

Microsatellites have gained wide application for elucidating population structure in nonmodel organisms. Since they are generally noncoding, neutrality is assumed but rarely tested. In Atlantic cod (Gadus morhua L.), microsatellite studies have revealed highly heterogeneous estimates of genetic differentiation among loci. In particular one locus, Gmo 132, has demonstrated elevated genetic differentiation. We investigated possible hitch-hiking selection at this and other microsatellite loci in Atlantic cod. We employed 11 loci for analysing samples from the Baltic Sea, North Sea, Barents Sea and Newfoundland covering a large part of the species' distributional range. The 'classical' Lewontin-Krakauer test for selection based on variance in estimates of F(ST) and (standardized genetic differentiation) revealed only one significant pairwise test (North Sea-Barents Sea), and the source of the elevated variance could not be ascribed exclusively to Gmo 132. In contrast, different variants of the recently developed ln Rtheta test for selective sweeps at microsatellite loci revealed a high number of significant outcomes of pair-wise tests for Gmo 132. Further, the presence of selection was indicated in at least one other locus. The results suggest that many previous estimates of genetic differentiation in cod based on microsatellites are inflated, and in some cases relationships among populations are obscured by one or more loci being the subject to hitch-hiking selection. Likewise, temporal estimates of effective population sizes in Atlantic cod may be flawed. We recommend, generally, to use a higher number of microsatellite loci to elucidate population structure in marine fishes and other nonmodel species to allow for identification of outlier loci that are subject to selection.     

Proceedings of the SMBE Tri-National Young Investigators' Workshop 2005. Molecular genetics of natural populations

Significant progress in evolutionary genetics has been made by studying, on the one hand, patterns of DNA sequence polymorphism and, on the other, genetic architecture of complex adaptive traits. However, connections between nucleotide variants under selection and adaptively relevant phenotypes are missing. Such connections can be established using precise gene replacement. We review the recent successful introduction of this technique to the analysis of two evolutionarily interesting loci--Odysseus and desaturase2. Both genes have subtle phenotypes that nevertheless could be identified using gene replacement, demonstrating that effects of naturally occurring alleles can be measured in the laboratory. This is an important first step in connecting statistical signatures of selection with adaptation in nature. More candidate genes involved in adaptation, for example, through cloning of genes responsible for reproductive isolation, now need to be identified. Molecular genetic manipulation, DNA polymorphism analysis, and field studies then have to be integrated to provide fresh insights into the mechanisms of evolutionary change.     

Tracing the first step to speciation: ecological and genetic differentiation of a salamander population in a small forest

Mechanisms and processes of ecologically driven adaptive speciation are best studied in natural situations where the splitting process is still occurring, i.e. before complete reproductive isolation is achieved. Here, we present a case of an early stage of adaptive differentiation under sympatric conditions in the fire salamander, Salamandra salamandra, that allows inferring the underlying processes for the split. Larvae of S. salamandra normally mature in small streams until metamorphosis, but in an old, continuous forest area near Bonn (the Kottenforst), we found salamander larvae not only in small streams but also in shallow ponds, which are ecologically very different from small streams. Common-environment experiments with larvae from both habitat types reveal specific adaptations to these different ecological conditions. Mitochondrial and microsatellite analyses show that the two ecologically differentiated groups also show signs of genetic differentiation. A parallel analysis of animals from a neighbouring much larger forest area (the Eifel), in which larvae mature only in streams, shows no signs of genetic differentiation, indicating that gene flow between ecologically similar types can occur over large distances. Hence, geographical factors cannot explain the differential larval habitat adaptations in the Kottenforst, in particular since adult life and mating of S. salamandra is strictly terrestrial and not associated with larval habitats. We propose therefore that the evolution of these adaptations was coupled with the evolution of cues for assortative mating which would be in line with models of sympatric speciation that suggest a co-evolution of habitat adaptations and associated mating signals.     

Insertion polymorphism of transposable elements and population structure of Anopheles gambiae M and S molecular forms in Cameroon

The insertion polymorphism of five transposable element (TE) families was studied by Southern blots in several populations of the M and S molecular forms of the mosquito Anopheles gambiae sensu stricto from southern Cameroon. We showed that the mean TE insertion site number and the within-population insertion site polymorphism globally differed between the M and S molecular forms. The comparison of the TE insertion profiles of the populations revealed a significant differentiation between these two molecular forms (0.163 < Phi(ST) < 0.371). We cloned several insertions of a non-LTR retrotransposon (Aara8) that were fixed in one form and absent in the other one. The only insertion that could be clearly located on a chromosome arm mapped to cytological division 6 of chromosome X, confirming the importance of this region in the ongoing speciation between the M and S molecular forms.     

Population genetics of Drosophila ananassae: inversion polymorphism and body size in Indian geographical populations

Twelve Indian natural populations of Drosophila ananassae, a cosmopolitan and domestic species, were sampled and laboratory populations (mass cultures) were established from naturally impregnated females. These populations were maintained in the laboratory for some generations and were analysed chromosomally to know the frequency of different inversions. The chromosomal analysis revealed the presence of three cosmopolitan inversions. The data on the whole show that there are significant differences in the frequencies of different chromosome arrangements in these populations. Body size (wing length and thorax length) was measured in both sexes (50 females and 50 males), in all the 12 geographical populations of D. ananassae. There are statistically significant differences in wing length as well as in thorax length of both sexes among different geographical populations. Five geographical strains were crossed reciprocally and body size (wing length and thorax length) was measured in F₁ and F₂ progeny. The comparison of body size (both traits) between mid‐parent, F₁ and F₂ shows that there is an increase in body size in F₁ and F₂ progeny as compared with parents. Thus, there is no break down of heterosis in F₂, which suggests absence of coadaptation in geographical populations of D. ananassae. Scaling test statistical analysis showed additive, dominance and epistatic effects in certain crosses involving geographical strains of D. ananassae. Correlation between chromosome arrangement frequency and body size has also been tested and significant negative correlation has been found between 2L – ST chromosome arrangement and male thorax.