Mitochondrial DNA diversity in an apomictic Daphnia complex from the Canadian High Arctic

Cyclic parthenogenesis is the ancestral mode of reproduction in the cladoceran crustacean, Daphnia pulex, but some populations have made the transition to obligate parthenogenesis and this is the only mode of reproduction known to occur in arctic populations. Melanism and polyploidy are also common in arctic populations of this species. Prior allozyme studies of arctic D. pulex revealed substantial levels of clonal diversity on a regional scale. Clonal groupings based on cluster analysis of allozyme genotypes do not conform to groupings based on the presence/absence of melanin or on ploidy level. In order to further elucidate genetic relationships among arctic D. pulex clones, mitochondrial DNA (mtDNA) variation was examined in 31 populations from two Canadian high-arctic sites. The data were also compared to a previous study of mtDNA variation in populations from a Canadian low-arctic site. Cladistic analysis of restriction site variation of the entire mitochondrial genome and nucleotide sequence variation of the mitochondrial control region was used to construct genetic relationships among mitochondrial genotypes. Three distinct mitochondrial lineages were detected. One lineage was associated with diploid, nonmelanic clones and is the same as the lineage that is found in temperate populations of D. pulex. The other two lineages (A & B) were associated with polyploid, melanic clones. Sequence divergence between the A and B lineages was 2.4%. Sequence divergence between D. pulex and either of these two lineages exceeded 3%. It is suggested that the melanic, polyploid clones are hybrids between males of D. pulex (and/or a closely related congener, D. pulicaria) and females of either of two ancestral melanic species that have mitochondrial lineages A and B. Geographic patterns of mitochondrial diversity in ‘melanic’ lineage B support the hypothesis of an high-arctic refuge for the ancestral species during the last glacial period.     

THE EXTENT OF INTROGRESSION OUTSIDE THE CONTACT ZONE BETWEEN NOTROPIS CORNUTUS AND NOTROPIS CHRYSOCEPHALUS (TELEOSTEI: CYPRINIDAE)

The cyprinid fishes, Notropis cornutus and N. chrysocephalus, hybridize in a long, narrow zone in the midwestern United States. To quantify the extent of introgression of genetic markers outside of this zone, samples were collected along transects starting near the region of contact (as defined by morphological characters), followed by samples progressively more distant. Diagnostic allozymic and mitochondrial DNA (mtDNA) restriction site markers were used to estimate the extent of introgression outside of the zone, while polymorphic allozyme and mtDNA markers were used to evaluate the potential for gene flow among populations within transects. Analysis of populations from the northern transect provided evidence for differentiation of populations for some of the markers; however, on average, enough gene flow has occurred to overcome substantial differentiation. Introgressed mtDNA and allozyme haplotypes were rare and found only in the population closest to the contact zone. The rarity of introgressed alleles in the more northern populations is consistent with the recent origin of these populations after the Wisconsin glaciation (less than 12,000 years bp) and/or selection maintaining the northern boundary of the contact zone. Analysis of populations from the southern transect revealed evidence for population subdivision but no evidence for introgression at the diagnostic allozyme loci; however, nearly all individuals from this transect possessed introgressed mtDNA haplotypes, with samples furthest from the contact zone exhibiting the highest frequencies of introgression. Patterns of variation for one of the polymorphic allozyme markers (Est-A) and introgressed mtDNAs were highly correlated, suggesting that allozymic heterogeneity at this locus is also the result of introgression. The most likely explanation for these data is that these introgressed haplotypes are indicators of a more southern position of the contact zone during the Pleistocene, with the contact zone shifting northward with the recession of the glacial front. Such movement implicates selection in the maintenance of distributional limits of these species, and hence, the width and position of the contact zone.     

Molecular genetic analysis of Daphnia in the English Lake District: species identity, hybridisation and resting egg banks

• 1 Taxonomic issues in many Daphnia species complexes are often confused by a high degree of phenotypic plasticity and by interspecific hybridisation. Here, we employ molecular genetic tools to confirm the species composition and incidence of hybridisation in extant and resting egg populations of Daphnia from Windermere and Esthwaite Water in the English Lake District.
• 2 A combination of species‐diagnostic allozymes and mtDNA, confirms that contemporary populations are dominated by a single species, D. galeata.
• 3 We present the first account of genetic characterisation of dated ‘resting’ eggs using microsatellites and mtDNA, employing PCR‐based DNA recovery, thus providing a temporal dimension to taxonomic patterns. Thirty years ago, two species were present in Esthwaite, D. galeata and D. hyalina, but Windermere populations were dominated by D. galeata only.
• 4 The use of PCR‐based mtDNA RFLP analysis as a species‐diagnostic tool, and microsatellites to monitor clonal diversity, provide a valuable approach for long‐term studies, especially in populations free from the complicating effect of frequent hybridisation. The detailed limnological records available for many large lakes, and associated changes in land‐use, pollutants and climate, combined with long‐term ephippial molecular genetic data, provide opportunities for exploring natural and anthropogenic impacts on genetic and community structure.

     

Molecular and morphological divergence in the butterfly genus Lycaeides (Lepidoptera: Lycaenidae) in North America: evidence of recent speciation

Male genital morphology, allozyme allele frequencies and mtDNA sequence variation were surveyed in the butterfly species Lycaeides idas and L. melissa from across much of their range in North America. Despite clear differences in male genital morphology, wing colour patterns and habitat characteristics, genetic variation was not taxonomically or geographically structured and the species were not identifiable by either genetic data set. Genetic distances (Nei's D=0.002–0.078, calculated from allozyme data) between all populations of both species were within the range commonly observed for conspecific populations of other butterflies. The most frequent mtDNA haplotype was present in individuals of both species in populations from southern California to Wisconsin. We conclude that speciation has probably happened recently and the lack of genetic differentiation between the species is the product of either (1) recent or ongoing gene flow at neutral loci, and/or (2) an insufficiency of time for lineage sorting. The evolution of male genital morphology, wing colour patterns and ecological characteristics has proceeded more rapidly than allozyme or mtDNA evolution.     

World phylogeography and male‐mediated gene flow in the sandbar shark,Carcharhinus plumbeus

The sandbar shark, Carcharhinus plumbeus, is a large, cosmopolitan, coastal species. Females are thought to show philopatry to nursery grounds while males potentially migrate long distances, creating an opportunity for male‐mediated gene flow that may lead to discordance in patterns revealed by mitochondrial DNA (mtDNA) and nuclear markers. While this dynamic has been investigated in elasmobranchs over small spatial scales, it has not been examined at a global level. We examined patterns of historical phylogeography and contemporary gene flow by genotyping 329 individuals from nine locations throughout the species’ range at eight nuclear microsatellite markers and sequencing the complete mtDNA control region. Pairwise comparisons often resulted in fixation indices and divergence estimates of greater magnitude using mtDNA sequence data than microsatellite data. In addition, multiple methods of estimation suggested fewer populations based on microsatellite loci than on mtDNA sequence data. Coalescent analyses suggest divergence and restricted migration among Hawaii, Taiwan, eastern and western Australia using mtDNA sequence data and no divergence and high migration rates, between Taiwan and both Australian sites using microsatellite data. Evidence of secondary contact was detected between several localities and appears to be discreet in time rather than continuous. Collectively, these data suggest complex spatial/temporal relationships between shark populations that may feature pulses of female dispersal and more continuous male‐mediated gene flow.     

EVIDENCE FOR A MOSAIC HYBRID ZONE IN THE GRASS SHRIMP PALAEMONETES KADIAKENSIS (PALAEMONIDAE) AS REVEALED BY MULTIPLE GENETIC MARKERS

Molecular techniques provide powerful tools for studying the geographic structure of hybrid zones and the dynamics of gene exchange between incipient species. We examined allozyme variation at five loci (PGM, GPI, MDH-1, MDH-2, and LDH) for 27 populations of Palaemonetes kadiakensis from the central, coastal, and eastern regions of Texas. Central Texas populations of P. kadiakensis exhibited highly significant linkage disequilibrium and departures from Hardy-Weinberg genotype proportions. In populations with linkage disequilibrium, allelic differences at GPI defined two types of P. kadiakensis, designated A and B. Both types existed in central Texas with little or no evidence of interbreeding, whereas the populations from all other localities showed complete introgression of type B alleles into the type A gene pool. We also examined ribosomal DNA (rDNA) and mitochondrial DNA (mtDNA) variation in a subset of populations, chosen to cover a range of geographic locations and levels of linkage disequilibrium. Two groups of mtDNA haplotypes and two restriction fragment patterns for the rDNA corresponded to allozyme type A and B individuals in populations exhibiting linkage disequilibrium. In populations with ongoing hybridization, all hybrid animals (N= 15) exhibited type A mtDNA. Exhibition of type A mtDNA indicated that type A females had mated successfully with type B males, but type B females had not mated successfully with type A males. Genotype distributions suggest reduced reproduction by hybrid offspring in central Texas populations. These patterns are consistent with a mosaic model of hybrid zone dynamics.     

Intra- and inter-specific relationships in the six Italian species of the fairy shrimp genus Chirocephalus: combining allozyme and mtDNA data

This study analysed the levels of genetic differentiation within and among the six Italian species of the fairy shrimp genus Chirocephalus by analysing electrophoretic polymorphisms at 22 enzymatic loci and by sequencing a 665‐bp fragment of the mitochondrial gene encoding for subunit I of cytochrome oxidase. The allozyme data revealed different levels of intra‐specific differentiation; mean θ estimates were low in Chirocephalus salinus, higher and comparable in C. diaphanus and C. kerkyrensis, while C. ruffoi was the most genetically structured species. At the inter‐specific level, C. marchesonii was the most differentiated species, both for allozymes and mtDNA. Phylogenetic relationships deduced from allozymes and mtDNA were not always consistent with each other. This highlights the differences in performance of the two classes of molecular markers and the need of different independent strategies of data analysis to search for possible incongruence. Neither allozymes nor mtDNA supported monophyly of the diaphanus‐group, previously recognized on the basis of the morphology of appendages (antennae and penes). In contrast the molecular results were consistent with the great heterogeneity in resting egg morphology among representatives of the diaphanus‐group.     

DEVELOPMENT OF THE HYBRID SWARM BETWEEN PECOS PUPFISH (CYPRINODONTIDAE: CYPRINODON PECOSENSIS) AND SHEEPSHEAD MINNOW (CYPRINODON VARIEGATUS): A PERSPECTIVE FROM ALLOZYMES AND mtDNA

A comparison of allozyme and mtDNA frequencies was used for insight into a situation in the Pecos River, Texas where contact between the endemic pupfish (Cyprinodon pecosensis) and an introduced congener (C. variegatus) has resulted in rapid, geographically extensive genetic introgression. Temporal changes in mean frequencies of diagnostic allozyme markers indicate that the clinal pattern of introduced genetic material (Echelle and Connor 1989) is slowly decreasing in amplitude. Significant rank concordance in diagnostic allele frequencies among sites and across sampling years indicates directional influences upon temporal allele frequency change. These observations are consistent with the theory of gene flow in neutral clines. Levels of introgression indicated by each of four allozyme loci and mtDNA were roughly equivalent. The early history of the hybrid swarm is explained by genetic swamping, possibly mediated by selection for C. variegatus or C. variegatus × C. pecosensis, at a time when the normally abundant endemic species had been catastrophically depleted. High frequencies of an introduced GPI-A allele in all samples of intergrades suggests that the introduced genome originated with a single founding event.     

Isolation‐driven divergence: speciation in a widespread North American songbird (Aves: Certhiidae)

Lineage, or true ‘species’, trees may differ from gene trees because of stochastic processes in molecular evolution leading to gene‐tree heterogeneity. Problems with inferring species trees because of excessive incomplete lineage sorting may be exacerbated in lineages with rapid diversification or recent divergences necessitating the use of multiple loci and individuals. Many recent multilocus studies that investigate divergence times identify lineage splitting to be more recent than single‐locus studies, forcing the revision of biogeographic scenarios driving divergence. Here, we use 21 nuclear loci from regional populations to re‐evaluate hypotheses identified in an mtDNA phylogeographic study of the Brown Creeper (Certhia americana), as well as identify processes driving divergence. Nuclear phylogeographic analyses identified hierarchical genetic structure, supporting a basal split at approximately 32°N latitude, splitting northern and southern populations, with mixed patterns of genealogical concordance and discordance between data sets within the major lineages. Coalescent‐based analyses identify isolation, with little to no gene flow, as the primary driver of divergence between lineages. Recent isolation appears to have caused genetic bottlenecks in populations in the Sierra Madre Oriental and coastal mountain ranges of California, which may be targets for conservation concerns.     

PHYLOGEOGRAPHY OF THE SPOTTED SAND BASS, PARALABRAX MACULATOFASCIATUS: DIVERGENCE OF GULF OF CALIFORNIA AND PACIFIC COAST POPULATIONS

Abstract Have the warm tropical waters and currents of the southern Gulf of California, Mexico (also known as the Sea of Cortez), formed a barrier to gene flow, resulting in disjunct populations in the upper gulf that are isolated from the outer Pacific Coast? Phylogeographic and genetic divergences of the spotted sand bass, Paralabrax maculatofasciatus, from three Gulf of California and two outer Pacific coastal locations were tested using mitochondrial DNA (mtDNA) control region sequences. Sequence data from two congeners that are sympatrically distributed along the outer Pacific Coast, the barred sand bass, P. nebulifer, and the kelp bass, P. clathratus, were used to gauge the levels of genetic divergences. Differences among the three species and between the northern gulf and outer Pacific coastal populations of P. maculatofasciatus also were analyzed using 40 allozymic presumptive gene loci. Allozyme and mtDNA analyses each revealed many fixed differences among the species. Three significant allozymic frequency differences and two fixed mtDNA substitutions differentiated the gulf and outer Pacific coastal populations of P. maculatofasciatus. Three unique mtDNA haplotypes and three unique allozyme alleles were identified from the outer Pacific coastal population. The gulf sites contained four unique mtDNA haplotypes and six unique allozyme alleles. Partitioning of the mtDNA variation revealed that 72% of the variance occurred between the gulf and outer Pacific Coast, 20% between sampling sites in the two regions, and 8% within the sites. There appears to be little gene flow across the waters of the southern Baja Penninsula, producing divergence estimated as 120,000 to 600,000 years between the outer Pacific coastal and the Gulf of California populations. This separation level may date to a hypothesized seaway closure near La Paz, Mexico, during the mid‐Pleistocene, and characterizes other fish populations. A second pattern of deeper allopatric species‐level divergences in some other fishes may date to a Pliocene closure of a mid‐Baja Penninsular seaway. Significant differences also were discerned in P. maculatofasciatus between the San Diego and central Baja California coastal sites and between the upper/central and the lower gulf locations. Variation between locations in the two regions may be indicative of larval retention and low adult migration, which needs to be tested further.     

Population and Species Divergence of Chemical Cues that Influence Male Recognition of Females in Desmognathine Salamanders

Growing evidence indicates that males may be more discriminating of mating partners than often has been assumed. In the North American Ocoee dusky salamander, Desmognathus ocoee (Plethodontidae: Desmognathinae), sexual incompatibility among conspecific populations is high in encounters staged in the laboratory, at least in part because males fail to recognize ‘other’ females as appropriate targets for courtship. I used Y-mazes to test the hypothesis that males of D. ocoee discriminate between substrate-borne chemical cues produced by ‘own’ (homotypic) and ‘other’ (heterotypic) females. Males of four populations discriminated in favor of substrates soiled by homotypic females over clean (control) substrates (expt 1), suggesting that females produce chemical cues of sociosexual significance to males. Furthermore, males from these populations discriminated in favor of substrates soiled by homotypic females vs. substrates soiled by heterotypic females (expt 2), both conspecific and heterospecific (D. carolinensis and D. orestes). Thus, differences among populations and species in female chemical cues appear to affect the chemotactic responses of males. I suggest that, together with differences in behavioral signals and responses exhibited during courtship, differences in female chemical cues likely contribute to sexual incompatibility among populations and taxa of desmognathine salamanders.     

SIGNIFICANT ROLE FOR HISTORICAL EFFECTS IN THE EVOLUTION OF REPRODUCTIVE ISOLATION: EVIDENCE FROM PATTERNS OF INTROGRESSION BETWEEN THE CYPRINID FISHES,LUXILUS CORNUTUS AND LUXILUS CHRYSOCEPHALUS

Samples of Luxilus cornutus, Luxilus chrysocephalus, and their hybrids were collected along hypothesized routes of dispersal from Pleistocene refugia to examine the significance of geographic variation in patterns of introgression between these species. Patterns of allozyme and mitochondrial DNA (mtDNA) variation were generally consistent with those from previous studies. Tests of Hardy-Weinberg equilibrium revealed significant deficiencies of heterozygotes in all samples, indicating some form of reproductive isolation. Mitochondrial DNAs of each species were not equally represented in F₁ hybrids; however, this bias was eliminated when the two largest samples were excluded from the analysis. Backcross hybrids exhibited biased mtDNA introgression, as samples from Lake Erie (eastern) and Lake Michigan (western) drainages showed significant excesses of mtDNAs from L. chrysocephalus and L. cornutus, respectively, relative to frequencies of diagnostic allozyme markers. The extent and direction of allozyme and mtDNA introgression was quantified by calculating isolation index values from morphologically “pure” individuals of each species from each locality. Analysis of variance of these measures identified limited introgression of allozyme variants with no geographic pattern, but significant differences in direction of mtDNA introgression between drainages (i.e., postglacial dispersal route). Association between patterns of mtDNA introgression and dispersal route across the latitudinal width of the contact zone is best explained by genetic divergence during past isolation of ancestral populations from these drainages. These results identify a significant role for historical effects in the evolution of reproductive isolation and the process of speciation.     

Lack of genetic differentiation in the fat snook Centropomus parallelus (Teleostei: Centropomidae) along the Brazilian coast

The patterns of genetic variation in a fragment of the cytochrome b mtDNA gene in 130 individuals of the fat snook (Centropomus parallelus) did not show evidence of genetic differentiation along the Brazilian coast as shown in F‐statistics and AMOVA analyses, indicating the existence of high gene flow among the studied populations. In addition, the study of the evolutionary demography of the species using mismatch distribution analysis and statistical tests such as Tajima’s D and Fu’s F indicate that it experienced a pattern of population expansion during its recent evolutionary history.     

One Species or Several? Discordant Patterns of Geographic Variation between Allozymes and mtDNA Sequences among Spiders in the Genus Metepeira (Araneae: Araneidae)

Paradoxically, an allozyme study of Metepeira “spinipes” (sensu lato) demonstrated extensive gene flow among four populations whose members are nevertheless morphologically and behaviorally distinct. Initially, the authors tentatively concluded that the populations exhibited panmixis and suggested that local environmental effects accounted for the apparent morphological and behavioral differences. However, they later concluded that such differences were too great to be accounted for by the environment alone and that the four populations actually represented three different species. To confirm that the allozyme results were, in fact, artifactual, we reexamined the relationships among these populations by sequencing a portion of the 12S mtDNA ribosomal subunit. In contrast to the allozyme result, our results demonstrate good agreement between patterns of genetic and morphological/behavioral variation. We suggest (1) that the allozyme allele frequencies are homogenized by balancing selection, not gene flow as was previously concluded, and therefore (2) that this study provides another instance in which inferences about population structures from allozyme data are misleading.     

Genetic Variation and Population Differentiation in Natural Populations of Cryptomeria japonica

Abstract Sugi (Cryptomeria japonica D. Don) is a valuable tree species in Japan. The present natural distribution is limited to small scattered areas in temperate moist regions, and most of these areas are surrounded by vast artificial plantations. We studied natural populations of C. japonica in an effort to determine the amount and distribution of genetic diversity using 12 allozyme markers. The amount of genetic variation within the species is high (H_T=0.196) but most is found within populations with little among populations (G_ST=0.034) despite their isolated distribution. This pattern of genetic diversity is inferred to be the consequence of the following: (1) the distribution of this species in the past was wider and more continuous than it is now; (2) a high rate of gene flow occurs, perhaps including gene flow between natural populations and plantations; and (3) the long lifespan. However, the distribution of allele frequencies at the 6Pg-1 in northern populations on the side near the Sea of Japan is clearly different from those in other populations. This observation is inferred to result from founding events.     

Ring species as demonstrations of the continuum of species formation

In the mid‐20th century, Ernst Mayr (1942) and Theodosius Dobzhansky (1958) championed the significance of ‘circular overlaps’ or ‘ring species’ as the perfect demonstration of the gradual nature of species formation. As an ancestral species expands its range, wrapping around a geographic barrier, derived taxa within the ring display interactions typical of populations, such as genetic and morphological intergradation, while overlapping taxa at the terminus of the ring behave largely as sympatric, reproductively isolated species. Are ring species extremely rare or are they just difficult to detect? What conditions favour their formation? Modelling studies have attempted to address these knowledge gaps by estimating the biological parameters that result in stable ring species (Martins et al. 2013), and determining the necessary topographic parameters of the barriers encircled (Monahan et al. 2012). However, any generalization is undermined by a major limitation: only a handful of ring species are known to exist in nature. In addition, many of them have been broken into multiple species presumed to be evolving independently, usually obscuring the evolutionary dynamics that generate diversity. A paper in this issue of Molecular Ecology by Fuchs et al. (2015), focused on the entire genealogy of a bulbul (Alophoixus) species complex, offers key insights into the evolutionary processes underlying diversification of this Indo‐Malayan bird. Their findings fulfil most of the criteria that can be expected for ring species (Fig.  1 ): an ancestor has colonized the mainland from Sundaland, expanded along the forested habitat wrapping around Thailand's lowlands, adjacent taxa intergrade around the ring distribution, and terminal taxa overlap at the ring closure. Although it remains unclear whether ring divergence has resulted in restrictive gene flow relative to that observed around the ring, their results suggest that circular overlaps might be more common in nature than currently recognized in the literature. Most importantly, this work shows that the continuum of species formation that Mayr and Dobzhansky praised in circular overlaps is found in biological systems currently described as ‘rings of species’, in addition to the idealized ‘ring species’.     

THE EVOLUTIONARY HISTORY OF PARTHENOGENETIC CNEMIDOPHORUS LEMNISCATUS (SAURIA: TEIIDAE). II. MATERNAL ORIGIN AND AGE INFERRED FROM MITOCHONDRIAL DNA ANALYSES

Restriction endonuclease analyses were performed on mitochondrial DNAs (mtDNAs) representing unisexual parthenogenetic (cytotypes A, B, and C) and bisexual (cytotypes D and E) populations of Amazonian lizards presently regarded as Cnemidophorus lemniscatus. The results of mtDNA cleavage map comparisons among these C. lemniscatus indicated that (1) there was no cleavage site variation among the unisexuals, (2) mtDNAs from the bisexual cytotypes D and E differed in sequence from one another by about 13%, and (3) mtDNAs from cytotypes A–C differed from those of cytotype D by about 5% and from those of cytotype E by about 13%. Higher resolution restriction fragment size comparisons confirmed the high degree of similarity among the unisexual mtDNAs, but identified 12 cleavage site variants among the 13 cytotype D mtDNAs examined. Both cladistic and phenetic (UPGMA) analyses of the data indicate that the unisexual and cytotype D mtDNAs form a single clade, suggesting that a female of cytotype D was the maternal progenitor of the unisexuals. The similarity among the unisexual mtDNAs and the variability among those of cytotype D suggest that the three unisexual cytotypes arose recently from a common maternal lineage. The mtDNA variability observed among cytotype D individuals has a strong geographic component, suggesting that the unisexuals arose from one or a few geographically proximal populations. The mtDNA comparisons also support the conclusion, based on allozyme comparisons (Sites et al., 1990, this issue), that cytotypes D and E, although presently allocated to C. lemniscatus, are separate species.     

Population structure and genetic diversity of metamorphic and paedomorphic populations of the tiger salamander,Ambystoma tigrinum

Genetic relationships, population subdivision and genetic diversity were estimated from mtDNA and allozyme data for two subspecies of tiger salamander, one of which is obligately metamorphic and the other polymorphic for paedomorphosis (larval reproduction). Far greater genetic differentiation exists between subspecies than within subspecies, suggesting that the subspecies have evolved in allopatry. Values of F_st calculated from both mtDNA and allozymes were greater than 0.400 for each subspecies. Significant population subdivision was detected even on a microgeographic scale. This extensive population subdivision indicates that populations can respond to extremely localized selection pressures. In the case of paedomorphosis, populations in permanent water should evolve paedomorphosis as long as the appropriate genes exist. For both mtDNA and allozymes, comparisons of population structure within the polymorphic subspecies and between polymorphic and metamorphic subspecies reveal no discernible effects of paedomorphosis. However, a comparison of paedomorphic and metamorphic populations of the polymorphic subspecies showed significantly higher mtDNA diversity in paedomorphic populations. The discrepancy between the allozyme and mtDNA results may be due to the lower effective population size of mtDNA compared to autosomal genes.     

Geography disentangles introgression from ancestral polymorphism in Lake Malawi cichlids

Phenotypically diverse Lake Malawi cichlids exhibit similar genomes. The extensive sharing of genetic polymorphism among forms has both intrigued and frustrated biologists trying to understand the nature of diversity in this and other rapidly evolving systems. Shared polymorphism might result from hybridization and/or the retention of ancestrally polymorphic alleles. To examine these alternatives, we used new genomic tools to characterize genetic differentiation in widespread, geographically structured populations of Labeotropheus fuelleborni and Metriaclima zebra. These phenotypically distinct species share mitochondrial DNA (mtDNA) haplotypes and show greater mtDNA differentiation among localities than between species. However, Bayesian analysis of nuclear single nucleotide polymorphism (SNP) data revealed two distinct genetic clusters corresponding perfectly to morphologically diagnosed L. fuelleborni and M. zebra. This result is a function of the resolving power of the multi‐locus dataset, not a conflict between nuclear and mitochondrial partitions. Locus‐by‐locus analysis showed that mtDNA differentiation between species (F_CT) was nearly identical to the median single‐locus SNP F_CT. Finally, we asked whether there is evidence for gene flow at sites of co‐occurrence. We used simulations to generate a null distribution for the level of differentiation between co‐occurring populations of L. fuelleborni and M. zebra expected if there was no hybridization. The null hypothesis was rejected for the SNP data; populations that co‐occur at rock reef sites were slightly more similar than expected by chance, suggesting recent gene flow. The coupling of numerous independent markers with extensive geographic sampling and simulations utilized here provides a framework for assessing the prevalence of gene flow in recently diverged species.     

Polytomies,signal and noise: revisiting the mitochondrial phylogeny and phylogeography of the Eurasian blindsnake species complex (Typhlopidae,Squamata)

Polytomies are multifurcating nodes on a phylogenetic tree that represent unresolved relationships. Contrary to ‘hard’ polytomies that hide simultaneous splitting events, ‘soft’ polytomies can theoretically be resolved with the addition of phylogenetic signal. This is not always successful, especially when a radiation is old and rapid, because adding more signal will inevitably increase the noise. The Xerotyphlops vermicularis species complex is an example of a largely unresolved old and rapid radiation. This mtDNA phylogeny is revisited by analysing samples representative of all lineages and generating a 3700‐bp final data set, after preliminary tests found two of the analysed markers responsible for long‐branch attraction. The new enhanced data set increased phylogenetic resolution, resulting in a robust time‐calibrated phylogeny and phylogeographic conclusions. Jordan/south Syria populations separated during the Middle Miocene, and the Messinian salinity crisis (MSC) appears responsible for the second diversification wave. A third radiation occurred during the Early Pliocene, resulting in four mtDNA groups west, east and south of the Amanos mountains in east Turkey/north Syria. Finally, three hypotheses regarding the number of potential ‘species’ within the complex are made, with the most moderate ‘four species’ hypothesis fitting very well with the diversity patterns and established systematics of east Mediterranean reptiles.