Discovery Of Genetically Distinct Sympatric Lineages in the Freshwater Mussel Cyprogenia Aberti (Bivalvia: Unionidae)

The western fanshell, Cyprogenia aberti, occurs in the CentralHighlands region of North America and is the only congener ofthe federally endangered fanshell, C. stegaria. Due to a reductionin range size and local extirpations, extant populations ofC. aberti have become increasingly isolated. Consequently, C.aberti has been the focus of numerous conservation efforts,yet no work has examined the degree of genetic variation amongthe species' disjunct populations. Phylogenetic analyses ofnucleotide sequences from two mitochondrial gene portions (CO1and ND1) revealed results important to the conservation managementof C. aberti. First, phylogenetic analyses do not support amonophyletic C. aberti. Second, C. aberti is comprised of asmany as five independent lineages, one of which includes thefederally endangered C. stegaria. Third, analysis recoveredtwo genetically distinct sympatric lineages of C. aberti inthe Ouachita River Drainage. These phylogenetic results suggestthat the five distinct evolutionary lineages of C. aberti shouldbe managed as separate conservation units. This study illustrateswhy it is critical to evaluate genetic variability in endangeredand threatened ‘species’ prior to implementing arecovery program. These data also reveal the value of assessingbiological diversity of nonimperiled taxa before populationsare extirpated and valuable genetic data necessary for reconstructionof evolutionary relationships is lost. (Received 24 October 2005; accepted 8 June 2006)     

Muskrat predation on a diverse unionid mussel community: impacts of prey species composition,size and shape

1. The muskrat, Ondatra zibethicus, is a predator known to affect species composition, population size and age structure of freshwater unionid mussel communities. Muskrats leave large piles of dead shells (middens) on the edges of rivers, lakes and streams. We compared the species composition and size structure of shells collected from muskrat middens to the nearby live unionid community in the lower Licking River (Kentucky, USA). 2. Like previous studies, our results suggest that muskrats are both size‐selective and species‐specific predators; however, our results indicate that mussel shape is also an important factor. 3. We generated a shape metric (‘cubocity’) sensitive to the overall shape of the mussel. Species with relatively lower cubocity values (around 0.85) are plate‐like or spike‐like, while mussels with more cuboidal shells have higher cubocity values (near 1.0). 4. Our results suggest muskrats prefer cuboidal mussels and generally avoid spike‐shaped mussels. The endangered fanshell, Cyprogenia stegaria, was the most favoured prey; the fanshell’s relative size and shape appear to make it particularly vulnerable to muskrats. 5. We believe the predictive capabilities of this shape metric will be of benefit to those who monitor and manage threatened mussel populations.     

Phylogeography and post-glacial colonization patterns of the rainbow darter, Etheostoma caeruleum (Teleostei: Percidae)

Aim To examine the effects of historical climate change and drainage isolation on the distribution of mitochondrial DNA cytochrome b genetic variation within the rainbow darter, Etheostoma caeruleum (Percidae: Etheostomatinae). Location Eastern North American streams including tributaries to the Mississippi River, Great Lakes, Potomac River and Hudson Bay drainages. Methods Parsimony analyses, Bayesian analyses and haplotype networks of mitochondrial DNA sequences. Results Four major clades were recovered from sampled populations of E. caeruleum. Three of four clades are distributed in the western portion of the species’ range (primarily west of the Mississippi River). Samples from this region do not form a monophyletic group, and sequences often vary greatly between samples from adjacent stream systems (up to 7.2% divergence). A basal clade includes samples from the White River system in the Ozark Highlands. The northern Ozarks–upper Midwest clade includes samples from Missouri River tributaries and the upper Midwest (Hudson Bay, upper Mississippi River, and western Lake Michigan drainage). The eastern clade is composed of individuals from the Ohio River, Great Lakes and Potomac River. The Mississippi River corridor clade includes samples from middle and lower Mississippi River tributaries. Main conclusions The four major clades of E. caeruleum are deep allopatric lineages with well‐defined boundaries and have additional phylogeographical structure within each clade. The Ozark Highlands have the greatest levels of diversity relative to distributional area, with marked cytochrome b subdivisions between adjacent stream systems. Samples from previously glaciated areas do not have a subset of the cytochrome b diversity found in unglaciated areas, but four separate source areas are identified based on phylogenetic analyses. Dispersal into previously glaciated areas followed several known glacial outlets and, based on sequence divergence between populations, may have occurred during different glacial or interglacial stages. The disjunct distribution and cytochrome b pattern of E. caeruleum in the Mississippi River corridor clade is consistent with late Pleistocene and Recent changes in the course and characteristics of the middle and lower Mississippi River. Phylogeographical boundaries between clades of E. caeruleum correspond to independent sources of biogeographical information and provide insight into historical stream drainage relationships, post‐glacial colonization and drainage isolation patterns.     

Positive selection on the mitochondrial ATP synthase 6 and the NADH dehydrogenase 2 genes across 22 hare species (genus Lepus)

Amino acid changes in mitochondrial (mt) oxidative phosphorylation (OXPHOS) genes have been suggested as a key adaptation to environmental variation. Here, we analyzed 416 sequences of ATPase synthase 6 (MT‐ATP6) and NADH dehydrogenase 2 (MT‐ND2) in 22 different hare (Lepus) species from across a wide range of habitats and climates. We used site‐ and branch‐based methods to test for positive selection on specific codons and lineages. We found four codons in MT‐ATP6 and five in MT‐ND2 under positive selection, affecting several species lineages. We investigated the association of protein variants at each locus with climate zone, using multinomial generalized linear models (glm), including species, regions, historical introgression events, and the co‐occurring protein variant at the other locus as additional explanatory variables. A significant climate effect as based on the “Köppen climate classification” was observed for MT‐ND2 protein variants as translated from our nucleotide sequences. Moreover, MT‐ND2 protein variants were significantly affected by the co‐occurring MT‐ATP6 protein variant in the same mtDNA molecule. Contrary to the expectation for non‐recombining mitochondrial DNA molecules, the presence of an evolutionarily relatively ancestral protein variant at one locus was associated with a relatively derived protein at the other locus in the same mitochondrial molecule, respectively. The relative evolutionary status of a protein variant was evaluated according to its positions relative to the respective out‐group protein variant in a network analysis of nucleotide sequences. All our results suggest a complex effect of various climatic parameters acting on multiple mtOXPHOS genes in a co‐adaptive way, favoring combinations of ancestral and derived variants.     

Homogeneity at nuclear microsatellite loci masks mitochondrial haplotype diversity in the endangered fanshell pearlymussel (Cyprogenia stegaria)

We report on multiple patterns of differentiation and connectivity in the fanshell pearlymussel (Cyprogenia stegaria), based on different markers. Knowledge of genetic variation and genetic connectivity among remaining populations of this federally endangered species is needed to initiate implementation of the species recovery plan. We collected tissue samples from 96 specimens from the Green, Rolling Fork, and Licking Rivers, tributaries to the Ohio River, and the Clinch River, a tributary to the Tennessee River, providing broad coverage of the current distributional range of the species. Results from 7 nuclear DNA microsatellite markers suggested minimal population-level differentiation, whereas a mitochondrial DNA (mtDNA) marker (ND1) exhibited significant differentiation between C. stegaria in the Clinch River and the Ohio River populations. The ND1 data also confirm the existence of 2 distinct mtDNA lineages in the genus that transcends species boundaries. Further analyses suggest that the disproportionally strong signal from 2 very divergent ND1 lineages possibly masks finer-grained structure in the Ohio River population, based on one of the mtDNA lineages only. We recommend further sampling to confirm the absence of one lineage from the upper Clinch River drainage and suggest that provisional management guidelines should limit reciprocal exchanges among C. stegaria populations from the Clinch River and those in the Ohio River system.     

Evidence for cryptic hybridization between different evolutionary lineages of the invasive clam genus Corbicula (Veneroida,Bivalvia)

We studied two Corbicula morphotypes in a syntopic population in the Rhine River in order to reveal their taxonomic, reproductive and phylogenetic relationship, using morphometrics, DAF‐fingerprinting, mitochondrial COI and nuclear ITS1 sequence variation. Morphometric analysis showed that two statistically distinguishable morphotypes with few intermediates were present.Mitochondrial sequence analysis detected two divergent clades. DAF‐fingerprinting revealed three highly distinctive multilocus genotypes. Two of the multilocus genotypes were significantly associated with different morphotypes and mitochondrial lineages. The third genotype B, however, was found in both morphotypes, intermediates and mitochondrial lineages. Conclusive evidence for hybridization came from RFLP analysis of the nuclear ITS1 locus. We interpret the hybrids as F1 hybrids between different evolutionary lineages. Integration of Corbicula sequences from all over the world into Maximum Parsimony analysis suggested a simultaneous radiation resulting in several evolutionary lineages whose species status remained doubtful. An unequivocal taxonomic assignment of the two evolutionary lineages in the Rhine population was therefore not possible.     

Phylogeography of the bigeye chub Hybopsis amblops (Teleostei: Cypriniformes): early Pleistocene diversification and post‐glacial range expansion

The bigeye chub, Hybopsis amblops, is a member of the Central Highlands ichthyofauna of eastern North America. Phylogenetic analyses of the H. amblops species group based on a 1059 bp fragment of the mitochondrial DNA cytochrome b gene did not recover a monophyletic group. The inclusion of Hybopsis hypsinotus in the species complex is questionable. Within H. amblops, five strongly supported clades were identified; two clades containing haplotypes from the Ozark Highlands and three clades containing haplotypes from the Eastern Highlands and previously glaciated regions of the Ohio and Wabash River drainages. Estimates of the timing of divergence indicated that prior to the onset of glaciation, vicariant events separated populations east and west of the Mississippi River. East of the Mississippi River glacial cycles associated with the blocking and rerouting of the Teays River system caused populations to be pushed southward into refugia of the upper Ohio River. Following the most recent Wisconsinan glaciation, populations expanded northward into previously glaciated regions and southward into the Cumberland River drainage. In the Ozarks, west of the Mississippi River, isolation of clades appears to be maintained by the lack of stream capture events between the upper Arkansas and the White River systems and a barrier formed by the Arkansas River.     

Positive and purifying selection in mitochondrial genomes of a bird with mitonuclear discordance

Diversifying selection on metabolic pathways can reduce intraspecific gene flow and promote population divergence. An opportunity to explore this arises from mitonuclear discordance observed in an Australian bird Eopsaltria australis. Across >1500 km, nuclear differentiation is low and latitudinally structured by isolation by distance, whereas two highly divergent, parapatric mitochondrial lineages (>6.6% in ND2) show a discordant longitudinal geographic pattern and experience different climates. Vicariance, incomplete lineage sorting and sex‐biased dispersal were shown earlier to be unlikely drivers of the mitonuclear discordance; instead, natural selection on a female‐linked trait was the preferred hypothesis. Accordingly, here we tested for signals of positive, divergent selection on mitochondrial genes in E. australis. We used codon models and physicochemical profiles of amino acid replacements to analyse complete mitochondrial genomes of the two mitochondrial lineages in E. australis, its sister species Eopsaltria griseogularis, and outgroups. We found evidence of positive selection on at least five amino acids, encoded by genes of two oxidative phosphorylation pathway complexes NADH dehydrogenase (ND4 and ND4L) and cytochrome bc₁ (cyt‐b) against a background of widespread purifying selection on all mitochondrial genes. Three of these amino acid replacements were fixed in ND4 of the geographically most widespread E. australis lineage. The other two replacements were fixed in ND4L and cyt‐b of the geographically more restricted E. australis lineage. We discuss whether this selection may reflect local environmental adaptation, a by‐product of other selective processes, or genetic incompatibilities, and propose how these hypotheses can be tested in future.     

Molecular patterns of differentiation in canyon treefrogs (Hyla arenicolor): evidence for introgressive hybridization with the Arizona treefrog (H. wrightorum) and correlations with advertisement call differences

Detection of genetic and behavioural diversity within morphologically similar species has led to the discovery of cryptic species complexes. We tested the hypothesis that US populations of the canyon treefrog (Hyla arenicolor) may consist of cryptic species by examining mate‐attraction signals among three divergent clades defined by mtDNA. Using a multi‐locus approach, we re‐analysed phylogenetic relationships among the three clades and a closely related, but morphologically and behaviourally dissimilar species, the Arizona treefrog (H. wrightorum). We found evidence for introgression of H. wrightorum’s mitochondrial genome into H. arenicolor. Additionally, the two‐clade topology based on nuclear data is more congruent with patterns of call variation than the three‐clade topology from the mitochondrial dataset. The magnitude of the call divergence is probably insufficient to promote isolation of the nuclear DNA‐defined clades should they become sympatric, but further divergence in call properties significant in species identification could promote speciation in the future.     

Population genetic history of the dreissenid mussel invasions: expansion patterns across North America

This study tests population genetic patterns across the Eurasian dreissenid mussel invasions of North America—encompassing the zebra mussel Dreissena polymorpha (1986 detection) and the quagga mussel D. rostriformis bugensis (detected in 1990, which now has largely displaced the former in the Great Lakes). We evaluate their source-spread relationships and invasion genetics using 9–11 nuclear microsatellite loci for 583 zebra mussels (21 sites) and 269 quagga mussels (12 sites) from Eurasian and North American range locations, with the latter including the Great Lakes, Mississippi River basin, Atlantic coastal waterways, Colorado River system, and California reservoirs. Additionally, mtDNA cytochrome b gene sequences are used to verify species identity. Our results indicate that North American zebra mussels originate from multiple non-native northern European populations, whereas North American quagga mussels trace to native estuaries in the Southern Bug and Dnieper Rivers. Invasive populations of both species show considerable genetic diversity and structure (zebra F _ST = 0.006–0.263, quagga F _ST = 0.008–0.267), without founder effects. Most newer zebra mussel populations have appreciable genetic diversity, whereas quagga mussel populations from the Colorado River and California show some founder effects. The population genetic composition of both species changed over time at given sites; with some adding alleles from adjacent populations, some losing them, and all retaining closest similarity to their original composition. Zebra mussels from Kansas and California appear genetically similar and assign to a possible origin from the St. Lawrence River, whereas quagga mussels from Nevada and California assign to a possible origin from Lake Ontario. These assignments suggest that overland colonization pathways via recreational boats do not necessarily reflect the most proximate connections. In conclusion, our microsatellite results comprise a valuable baseline for resolving present and future dreissenid mussel invasion pathways.     

Deceptive single‐locus taxonomy and phylogeography: Wolbachia‐associated divergence in mitochondrial DNA is not reflected in morphology and nuclear markers in a butterfly species

The satyrine butterfly Coenonympha tullia (Nymphalidae: Satyrinae) displays a deep split between two mitochondrial clades, one restricted to northern Alberta, Canada, and the other found throughout Alberta and across North America. We confirm this deep divide and test hypotheses explaining its phylogeographic structure. Neither genitalia morphology nor nuclear gene sequence supports cryptic species as an explanation, instead indicating differences between nuclear and mitochondrial genome histories. Sex‐biased dispersal is unlikely to cause such mito‐nuclear differences; however, selective sweeps by reproductive parasites could have led to this conflict. About half of the tested samples were infected by Wolbachia bacteria. Using multilocus strain typing for three Wolbachia genes, we show that the divergent mitochondrial clades are associated with two different Wolbachia strains, supporting the hypothesis that the mito‐nuclear differences resulted from selection on the mitochondrial genome due to selective sweeps by Wolbachia strains.     

Phylogeny of slug species of the genus Arion: evidence of monophyly of Iberian endemics and of the existence of relict species in Pyrenean refuges

The Iberian Peninsula contains the majority of the Paleartic land slug species of the genus Arion, which exhibits diverse taxonomic problems. The present study investigated Arion taxonomy on the basis of analyses of the mitochondrial ND1 gene and nuclear internal transcribed spacer 1 (ITS1) sequences. The Iberian endemic species were monophyletically clustered in two divergent sister clades. The topotype specimens of Arion lusitanicus and the closely related species Arion nobrei and Arion fuligineus, as well as Arion hispanicus and Arion flagellus, were grouped into an ‘Atlantic’ clade, whereas Arion baeticus, Arion gilvus, Arion anguloi, Arion wiktori and Arion paularensis were included in a ‘Continental–Mediterranean’ clade. Calibration of mutation rate in the ND1 gene suggested that the divergence of these two clades occurred around the Pliocene–Pleistocene boundary, with subsequent speciation events during the Pleistocene. A group of ancestral and divergent endemic species with distribution centred in the Pyrenean mountain range (Arion molinae, Arion lizarrusti, Arion antrhacius and Arion iratii) arose in the Pliocene and survived through the Pleistocene in geographically confined small populations. Arion lusitanicus showed up to be polyphyletic: specimens, sampled outside the geographic range of the topotype in the north‐western Iberian Peninsula, were included in a non‐monophyletic clade together with the widely distributed species Arion ater and Arion rufus. The divergent species with a wide European distribution (Arion subfuscus, Arion hortensis, Arion fagophilus and Arion intermedius) were located in basal positions in all topologies. The evolutionary history of these slug species (highly sensitive to climatic factors, with capacity for both outcrossing and selfing, and with low dispersal ability) appears to have been moulded by Pliocene–Pleistocene climate events and by the rugged topography of southern Europe, giving rise to repeated cycles of population isolation during periods of glaciation alternating with interglacial expansions limited by geographic barriers.     

Genetic diversity and evolutionary history of the Schizothorax species complex in the Lancang River (upper Mekong)

The genus Schizothorax (Cyprinidae), one of the most diverse genera of ichthyofauna of the Qinghai‐Tibetan Plateau (QTP), is a good candidate for investigating patterns of genetic variation and evolutionary mechanisms. In this study, sequences from the mitochondrial control region, the cytochrome b gene, and two nuclear genes were used to re‐examine the genetic diversity and investigate the evolutionary history of the Schizothorax species complex inhabiting the Lancang River. Three maternal clades were detected in the Schizothorax species complex, but frequent nuclear allele sharing also occurred among the three maternal clades. A discrepancy between topologies of mitochondrial and nuclear loci might result from introgression or/and incomplete lineage sorting. The divergence of the clades of the Schizothorax species complex was closely related to the Late Pliocene and Early Pleistocene orogenesis of the QTP and Southwest Mountains of China. Demographic analyses indicated that the species complex subsequently persisted in situ with stable populations during Pleistocene glacial cycling, which suggested that Pleistocene climate changes did not exert a remarkable influence on the species complex. Our study provides a comprehensive analysis of the genetic diversity and evolutionary history of the Schizothorax species complex in the Lancang River.     

Phylogeography of the northern hogsucker, Hypentelium nigricans (Teleostei: Cypriniformes): genetic evidence for the existence of the ancient Teays River

Aim To assess the roles of dispersal and vicariance in shaping the present distribution and diversity within Hypentelium nigricans, the northern hogsucker (Teleostei: Cypriniformes). Location Eastern United States. Methods Parsimony analyses, Bayesian analyses, pairwise genetic divergence and mismatch plots are used to examine patterns of genetic variation across H. nigricans. Results Species relationships within the genus Hypentelium were consistent with previous hypotheses. However, relationships between haplotypes within H. nigricans revealed two deeply divergent groups, a clade containing haplotypes from the New and Roanoke rivers (Atlantic Slope) plus Interior Highlands and upper Mississippi River and a clade containing haplotypes from the Eastern Highlands, previously glaciated regions of the Ohio and Wabash rivers, and the Amite and Homochitto rivers of south‐western Mississippi. Main conclusions The phylogenetic history of Hypentelium was shaped by old vicariant events associated with erosion of the Blue Ridge and separation of the Mobile and Mississippi river basins. Within H. nigricans two clades existed prior to the Pleistocene; a widespread clade in the pre‐glacial Teays‐Mississippi River system and a clade in Cumberland and Tennessee rivers. Pleistocene events fragmented the Teays‐Mississippi fauna. Following the retreat of the glaciers H. nigricans dispersed northward into previously glaciated regions. These patterns are replicated in other clades of fishes and are consistent with some of the predictions of Mayden's (Systematic Zoology, 37, 329, 1988) pre‐Pleistocene vicariance hypothesis.     

Genome‐wide mosaicism in divergence between zoonotic malaria parasite subpopulations with separate sympatric transmission cycles

Plasmodium knowlesi is a significant cause of human malaria transmitted as a zoonosis from macaque reservoir hosts in South‐East Asia. Microsatellite genotyping has indicated that human infections in Malaysian Borneo are an admixture of two highly divergent sympatric parasite subpopulations that are, respectively, associated with long‐tailed macaques (Cluster 1) and pig‐tailed macaques (Cluster 2). Whole‐genome sequences of clinical isolates subsequently confirmed the separate clusters, although fewer of the less common Cluster 2 type were sequenced. Here, to analyse population structure and genomic divergence in subpopulation samples of comparable depth, genome sequences were generated from 21 new clinical infections identified as Cluster 2 by microsatellite analysis, yielding a cumulative sample size for this subpopulation similar to that for Cluster 1. Profound heterogeneity in the level of intercluster divergence was distributed across the genome, with long contiguous chromosomal blocks having high or low divergence. Different mitochondrial genome clades were associated with the two major subpopulations, but limited exchange of haplotypes from one to the other was evident, as was also the case for the maternally inherited apicoplast genome. These findings indicate deep divergence of the two sympatric P. knowlesi subpopulations, with introgression likely to have occurred recently. There is no evidence yet of specific adaptation at any introgressed locus, but the recombinant mosaic types offer enhanced diversity on which selection may operate in a currently changing landscape and human environment. Loci responsible for maintaining genetic isolation of the sympatric subpopulations need to be identified in the chromosomal regions showing fixed differences.     

Mitochondrial DNA sequences suggest unexpected phylogenetic position of Corso-Sardinian grass snakes (<Emphasis Type="Italic">Natrix cetti</Emphasis>) and do not support their species status,with notes on phylogeography and subspecies delineation of grass snakes

We supplement a previously published mitochondrial DNA data set of grass snake sequences (ND1, ND2, ND4, cyt b, in total 3,806 bp) with sequences of Corso-Sardinian and Tuscan specimens and infer their phylogeny using Bayesian, maximum likelihood and maximum parsimony methods. In addition, we estimate divergence times of grass snake clades using a relaxed molecular clock calibrated with fossil evidence, and, in a second approach, the post-Messinian reopening of the Strait of Gibraltar. Recently it was suggested that Corso-Sardinian grass snakes represent a distinct species: Natrix cetti. All tree-building methods revealed well-supported branching patterns and deep divergences among grass snakes. However, sequences of N. natrix were consistently paraphyletic with respect to Corso-Sardinian sequences. The sister group of Corso-Sardinian grass snakes is a clade embracing N. n. helvetica and N. n. lanzai. Extensive gene flow between N. n. helvetica and a more distantly related subspecies (N. n. natrix) is well known, which is why we conclude that the status of Corso-Sardinian grass snakes as subspecies of N. natrix should be reinstated. Many currently recognized grass snake subspecies conflict with mitochondrial clades, suggestive of inappropriate morphological taxon delineation and mitochondrial introgression. Divergences among grass snakes are old, and the results of the two independent dating approaches are largely congruent. Accordingly, the Alpine orogenesis seems to have caused the origin of the oldest clade, corresponding to Iberian N. n. astreptophora. The formation of Corso-Sardinian grass snakes was dated to the Early Pliocene and could result from post-Messinian flooding of the Mediterranean Basin. Another deeply divergent clade of approximately the same age, endemic in central and northern Europe, suggests the Pleistocene survival of grass snakes north of the Alps. At least one glacial refuge in which old lineages survived Pleistocene cold periods was located on each of the three major southern European peninsulas and in Anatolia. Due to pronounced sequence divergences among Italian and southern Swiss grass snakes, we hypothesize multiple refugia south of the Alps and in the Apennine Peninsula, and there is evidence for two refuges on the Balkan Peninsula.     

Phylogeny of the green lacewing Chrysoperla nipponensis species‐complex (Neuroptera: Chrysopidae) in China,based on mitochondrial sequences and AFLP data

Abstract Several studies have indicated that the green lacewing, Chrysoperla nipponensis (Neuroptera: Chrysopidae) may include more than one valid species. We investigated the phylogenetic status of Chrysoperla nipponensis s.l. in China and Japan using mitochondrial sequences and AFLP data. The molecular phylogenetic analyses based on mitochondrial genes showed that the C. nipponensis species‐complex comprises four clades, each having high support values. In addition, the phylogenetic tree based on AFLP data indicates that the species‐complex comprises three groups. These results confirm that C. nipponensis s.l. comprises at least three genetically distinct clades and suggests that two of these clades may be closely related to populations of C. nipponensis in Japan. However, these clades cannot be recognized as species until analysis of their courtship songs has been completed.     

Dual influence of terrestrial and marine historical processes on the phylogeography of the Brazilian intertidal red alga Gracilaria caudata

In this study, we explored how past terrestrial and marine climate changes have interacted to shape the phylogeographic patterns of the intertidal red seaweed Gracilaria caudata, an economically important species exploited for agar production in the Brazilian north‐east. Seven sites were sampled along the north‐east tropical and south‐east sub‐tropical Brazilian coast. The genetic diversity and structure of G. caudata was inferred using a combination of mitochondrial (COI and cox2‐3), chloroplast (rbcL) and 15 nuclear microsatellite markers. A remarkable congruence between nuclear, mitochondrial and chloroplast data revealed clear separation between the north‐east (from 03° S to 08° S) and the south‐east (from 20° S to 23° S) coast of Brazil. These two clades differ in their demographic histories, with signatures of recent demographic expansions in the north‐east and divergent populations in the south‐east, suggesting the maintenance of several refugia during the last glacial maximum due to sea‐level rise and fall. The Bahia region (around 12° S) occupies an intermediate position between both clades. Microsatellites and mtDNA markers showed additional levels of genetic structure within each sampled site located south of Bahia. The separation between the two main groups in G. caudata is likely recent, probably occurring during the Quaternary glacial cycles. The genetic breaks are concordant with (i) those separating terrestrial refugia, (ii) major river outflows and (iii) frontiers between tropical and subtropical regions. Taken together with previously published eco‐physiological studies that showed differences in the physiological performance of the strains from distinct locations, these results suggest that the divergent clades in G. caudata correspond to distinct ecotypes in the process of incipient speciation and thus should be considered for the management policy of this commercially important species.     

Development of microsatellite loci from a reference genome for the Neotropical butterfly Heliconius numata and its close relatives

The Neotropical butterfly Heliconius numata (Lepidoptera: Nymphalidae: Heliconiinae) is known for its striking diversity of wing color patterns driven by the Müllerian mimicry of multiple local models and controlled by a single supergene locus. Such fine‐scale variation of traits under strong selection offers a unique opportunity for the study of the ecology and genetics of adaptation. However, little is still known of the population processes driving geographical variation in wing‐pattern phenotypes. We report the characterization of 26 microsatellite markers for the butterfly H. numata, including six located inside the wing color‐pattern supergene region. All markers are polymorphic, with allele numbers ranging from 2 to 21 per locus, an observed heterozygosity of 0.111 to 0.848 and an expected heterozygosity of 0.126 to 0.942. A subset of 18 of these markers was tested on five closely related sympatric Heliconius species with an amplification success ranging from 88% to 94%. The obtained set of microsatellite markers provides a new and useful set of tools to investigate patterns of differentiation and selection in populations of mimetic Heliconius butterflies. Moreover, markers developed within the color‐pattern supergene will facilitate characterization of the association between the genetic architecture and the functional diversity of wing patterns. Finally, the cross‐species amplification success of the described markers extends their utility to also encompass comparative population genetic studies of closely related species within a clade of rapidly diversifying species.     

Population Structure of the Giant Garter Snake, Thamnophis gigas

The giant garter snake, Thamnophis gigas, is a threatened species endemic to California’s Central Valley. We tested the hypothesis that current watershed boundaries have caused genetic differentiation among populations of T. gigas. We sampled 14 populations throughout the current geographic range of T. gigas and amplified 859 bp from the mitochondrial gene ND4 and one nuclear microsatellite locus. DNA sequence variation from the mitochondrial gene indicates there is some genetic structuring of the populations, with high F_ST values and unique haplotypes occurring at high frequency in several populations. We found that clustering populations by watershed boundary results in significant between-region genetic variance for mtDNA. However, analysis of allele frequencies at the microsatellite locus NSU3 reveals very low F_ST values and little between-region variation in allele frequencies. The discordance found between mitochondrial and microsatellite data may be explained by aspects of molecular evolution and/or T. gigas life history characteristics. Differences in effective population size between mitochondrial and nuclear DNA, or male-biased gene flow, result in a lower migration rate of mitochondrial haplotypes relative to nuclear alleles. However, we cannot exclude homoplasy as one explanation for homogeneity found for the single microsatellite locus. The mitochondrial nucleotide sequence data supports conservation practices that identify separate management units for T. gigas.