Duplication and divergence: the evolution of nematode globins

In common with many other groups, nematodes express globins with unknown functions. Nematode globin-like genes can be divided into class 1 globins, similar to vertebrate myoglobins, and a wide range of additional classes. Here we show that class 1 nematode globins possess a huge amount of diversity in gene sequence and structure. There is evidence for multiple events of gene duplication, intron insertion and loss between species, and for allelic variation effecting both synonymous and non-synonymous sites within species. We have also examined gene expression patterns in class I globins from a variety of species. The results show variation in the degree of gene expression, but the tissue specificity and temporal specificity of expression may be more conserved in the phylum. Because the structure-function relationships for the binding and transport of oxygen by globins are well understood, the consequences of genetic variation causing amino acid changes are explored. The gene family shows great promise for discovering unique insights into both structure-function relationships of globins and their physiologial roles.     

Molecular evolution of the carbonic anhydrase genes: Calculation of divergence time for mouse carbonic anhydrase I and II

Summary A cDNA clone in pBR322 that cross-hybridizes with a mouse carbonic anhydrase form II (CAII) probe has been sequenced and identified as mouse carbonic anhydrase form I (CAI). The 1224-base-pair clone encodes the entire 260-amino-acid protein and appears to contain an Alu-like element in the 3 untranslated region. The deduced amino acid sequence exhibits 77% homology to human CAI and contains 17 of the 20 residues that are considered unique to and invariant for all mammalian CAI isozymes. The results of a detailed comparison of the nucleic acid sequences spanning the coding regions of mouse CAI and rabbit CAI have been used to calibrate an evolutionary clock for the carbonic anhydrases (CAs). These data have been applied to a comparison of the mouse CAI and CAII nucleic acid sequences to calculate the divergence time between the two genes. The divergence-time calculation provides the first estimation of the evolutionary relationship between CAs based entirely on nucleotide sequence comparison.     

Presence and partial characterisation of a major proteolytic enzyme in the larval gut of Callosobruchus maculatus

Protease activity in the gut of larvae of the bruchid beetle, Callosobruchus maculatus (a storage pest of cowpea seeds), has been investigated to help clarify nutritional mechanisms in view of reports that these insects carry out little or no proteolysis (Applebaum, 1964). Larval gut homogenates showed protease activity against a variety of different protein substrates, but did not hydrolyse a synthetic trypsin substrate. The proteolytic activity had a pH optimum of 5.4. It was not inhibited by serine protease inhibitors, but was inhibited by reagents reactive against-SH groups. Protein trypsin inhibitors from legume seeds which are not hosts to C. maculatus (soybean, limabean) were not effective inhibitors of the larval proteolytic activity but a cowpea protease inhibitor preparation and aprotinin partially inhibited proteolysis. The latter two inhibitors also inhibited the plant thiol protease papain. It is suggested that C. maculatus has replaced the normal insect proteases with an enzyme similar to plant proteases to evade the antimetabolic effects of trypsin/chymotrypsin inhibitors in seeds. Besides trypsin/chymotrypsin inhibitors, cowpea seeds also contain proteins which inhibit papain; these inhibitors were purified and were shown to be effective inhibitors of C. maculatus larval protease.

---

Mise en évidence et caractérisation partielle d'une enzyme protéolytique importante du tube digestif des larves de Callosobruchus maculatus

Résumé Callosobruchus maculatus attaque les graines stockées particulièrement de Vigna unguiculata. L'activité protéolytique du tube digestif des larves a été examinée. Aucune hydrolyse n'a été observée contre la N-benzoyl arginine p-nitroanilide (substrat de la trypsine synthétique), mais la protéolyse a été mise en évidence en utilisant la viciline de V. unguiculata comme substrat naturel ou la myoglobine comme substrat artificiel, le p.H optimum est de 5,4. Les inhibiteurs chimiques des protéases de la sérine n'ont pas altéré l'activité enzymatique, mais les réactifs du groupe SH l'ont inhibée. Nous suggérons que cette protéase est une thiol protéase, c'est à dire avec un groupe actif semblable à celui de la papaïne, protéase végétale. En accord avec cette hypothèse, les inhibiteurs de protéine efficaces seulement contre la trypsine, c'est à dire les inhibiteurs de la trypsine de Glycine max et de la trypsine de Phaseolus lunatus n'affectent pas cette protéase larvaire, tandis que les protéines avec une faible action inhibitrice contre la papaïne (préparation inhibitrice de la protéase de V. unguiculata, aprotinine) inhibent partiellement la protéolyse par des extraits de tube digestif larvaire.Des inhibiteurs protéiques de papaïne extraits de graines de V. unguiculata sont des inhibiteurs très efficaces de la protéase larvaire. On peut penser que C. maculatus contient une protéase semblable aux protéases végétales, plutôt qu'aux protéases classiques, d'insectes, ce qui permet d'éviter les effets antimétaboliques directs des inhibiteurs de protéase (inhibiteurs de trypsine et de chymotrypsine) trouvés en quantité relativement importantes dans les graines de légumineuses.

     

A Novel Neutral Protease from <Emphasis Type="Italic">Thermoactinomyces</Emphasis> Species 27a: Sequencing of the Gene,Purification, and Characterization of the Enzyme

The nucleotide sequence of the previously cloned (Zabolotskaya, M. V., Nosovskaya, E. A., Kaplun, M. A., and Akimkina, T. V. (2001). Mol. Gen. Mikrobiol. Virusol. No 1, 32–34) DNA fragment from Thermoactinomyces sp. 27a (GenBank Accession No. AY280367) containing the metalloproteinase gene was determined. A continuous open reading frame encoding a polypeptide of 673 aa was revealed. Analysis of this sequence demonstrated that the metalloproteinase from Thermoactinomyces sp. 27a is synthesized as a preproprotein and includes a leader peptide (26 aa), N-terminal propeptide (215 aa), mature region (317 aa), and additional C-terminal domain (115 aa). The recombinant enzyme from Thermoactinomyces sp. 27a was expressed in Bacillus subtilis AJ73 cells and purified by anion exchange chromatography to an electrophoretically homogeneous state. The determined N-terminal amino acid sequence of the mature protein was identical to that deduced from the gene. The obtained data suggest that the mature protein should include 432 aa and have a calculated molecular weight of 46,262 Da. However, the molecular weight of the mature protein determined by mass spectrometry was 34,190 ± 70 Da indicating a C-terminal processing. Theproteinase was not inhibited by phenylmethyl sulfonyl fluoride but was inhibited by o-phenanthroline and ethylenediaminetetraacetic acid. The enzyme had maximum activity by azocasein hydrolysis at 55°C and pH 6.5–7.5; it was stable at pH 7.5–8.5 and remained stable at 50°C for several hours. The k_cat/K_m for 3-(2-furyl)acryloyl-glycyl-L-leucine amide hydrolysis was (2.8 ± 0.1) ×10³ M^–1×s^–1.     

Molecular phylogeny of the Canary Island lacertids (Gallotia): mitochondrial DNA restriction fragment divergence in relation to sequence divergence and geological time

Restriction fragment length polymorphisms of 6 base pair recognising endonucleases are used to reconstruct the phylogeny of the endemic Canary Island lacertid, Gallotia. The division into conventional species is upheld by this molecular analysis and the western Canary Island lizard (G. galloti) and eastern Canary Island lizard (G. atlantica) are hypothesized to be sister species. A more comprehensive study of the intraspecific relationships of G. galloti, based on nineteen restriction enzymes, indicates that there are distinct southern and northern lineages within this species. The phylogenetic analysis does not uphold the conventional subspecies, but suggests an alternative arrangement with one northern (La Palma, Tenerife) and one southern (Gomera, Hierro) subspecies. The inferred timing of molecular divergence of populations of G. galloti, based on RFLP analysis, is compatible with the geological timing for island origin and fossil data. Mantel tests show that mitochondrial RFLP divergence is correlated with mitrochondrial 12s rRNA and cytochrome oxidase I sequence divergence and highly correlated with mitochondrial cytochrome b sequence divergence.     

Phenotypic distribution models corroborate species distribution models: A shift in the role and prevalence of a dominant prairie grass in response to climate change

Phenotypic distribution within species can vary widely across environmental gradients but forecasts of species’ responses to environmental change often assume species respond homogenously across their ranges. We compared predictions from species and phenotype distribution models under future climate scenarios for Andropogon gerardii, a widely distributed, dominant grass found throughout the central United States. Phenotype data on aboveground biomass, height, leaf width, and chlorophyll content were obtained from 33 populations spanning a ~1000 km gradient that encompassed the majority of the species’ environmental range. Species and phenotype distribution models were trained using current climate conditions and projected to future climate scenarios. We used permutation procedures to infer the most important variable for each model. The species‐level response to climate was most sensitive to maximum temperature of the hottest month, but phenotypic variables were most sensitive to mean annual precipitation. The phenotype distribution models predict that A. gerardii could be largely functionally eliminated from where this species currently dominates, with biomass and height declining by up to ~60% and leaf width by ~20%. By the 2070s, the core area of highest suitability for A. gerardii is projected to shift up to ~700 km northeastward. Further, short‐statured phenotypes found in the present‐day short grass prairies on the western periphery of the species’ range will become favored in the current core ~800 km eastward of their current location. Combined, species and phenotype models predict this currently dominant prairie grass will decline in prevalence and stature. Thus, sourcing plant material for grassland restoration and forage should consider changes in the phenotype that will be favored under future climate conditions. Phenotype distribution models account for the role of intraspecific variation in determining responses to anticipated climate change and thereby complement predictions from species distributions models in guiding climate adaptation strategies.     

Relationships of the chromosomal species in the eurasian mole rats of theSpalax ehrenbergi group as determined by DNA-DNA hybridization,and an estimate of the spalacid-murid divergence time

Summary DNA-DNA hybridization was used to measure the average genomic divergence among the four chromosomal species of the Eurasian mole rats belonging to theSpalax ehrenbergi complex (Rodentia: Spalacidae). The percent nucleotide substitutions in the single-copy nuclear DNA among the species ranged from 0 to 5%, suggesting that speciation has occurred with minor genomic changes in these animals. The youngest chromosomal species appear to differ by 0.2–0.6% base pair mismatch, which is only between one and three base differences in a 500-bp fragment. The interspecific values of percent nucleotide differences permit the recognition of two well-separated speciation events in theS. ehrenbergi complex, the older (of Lower Pleistocene age) having isolated the chromosomal species 2n=54 before the divergence of the three other species.DNA-DNA hybridization was also used to compare the Spalacinae (Eurasian mole rats), Murinae (Old World rats and mice), and Arvicolinae (voles and lemmings). These data enabled us to estimate the time of divergence of the spalacids at ca. 19 million years ago. The dates of divergence among the other rodent lineages, as predicted by DNA hybridization results, agree well with paleontological data. These dates of divergence are obtained by the relation between geological time and single-copy nuclear DNA change, a relation that was calibrated by Catzeflis et al. (1987) through the use of fossil Arvicolinae and Murinae data.     

Effects of vicariant barriers, habitat stability, population isolation and environmental features on species divergence in the south-western Australian coastal reptile community

Identifying explicit hypotheses regarding the factors determining genetic structuring within species can be difficult, especially in species distributed in historically dynamic regions. To contend with these challenges, we use a framework that combines species distribution models, environmental data and multi-locus genetic data to generate and explore phylogeographic hypotheses for reptile species occupying the coastal sand-dune and sand-plain habitats of the south-western Australian biodiversity hotspot, a community which has both a high diversity of endemics and has varied dramatically in spatial extent over time. We use hierarchical amova, summary statistic and distance-based analyses to explicitly test specific phylogeographic hypotheses. Namely, we test if biogeographic vicariance across barriers, habitat stability, population isolation along a linear habitat or fragmentation across different environments can explain genetic divergence within five co-distributed squamate reptile species. Our results show that patterns of genetic variation reflect complex and species-specific interactions related to the spatial distribution of habitats present currently and during repeated glacial minima, as opposed to being associated with historical factors such as habitat stability between glacial and inter-glacial periods or vicariant barriers. We suggest that the large impact of habitat characteristics over time (i.e. relative levels of habitat connectivity, climatic gradients and spatial heterogeneity of soil types) reflects the ecological restrictions of the sand-dune and sand-plain reptile communities and may explain the lack of concordance across taxa. The study demonstrates the general utility of the approach for assemblage-level, as well as single species, phylogeographic study, including its usefulness for exploring biologically informed hypotheses about what factors have influenced patterns of genetic variation.     

High molecular and phenotypic diversity in the Merodon avidus complex (Diptera,Syrphidae): cryptic speciation in a diverse insect taxon

This paper examines molecular and phenotypic variability in the widely spread European hoverfly species complex Merodon avidus. Herein, based on the mitochondrial DNA (mtDNA) sequences of the cytochrome c oxidase subunit I (COI) and morphometric wing parameters, M. avidus is shown to comprise a complex of cryptic species, and one variety is redefined as a valid species: M. bicolor Gil Collado, 1930 (as var. of spinipes). The species M. bicolor, M. avidus A, and M. avidus B were clearly delimited based on their wing size. A total of 29 M. avidus and M. bicolor individuals presented 20 mtDNA haplotypes, four of which were shared by M. avidus A and M. avidus B, three were confined to M. bicolor, seven to M. avidus A, and six to M. avidus B. Sequence divergences between lineages occurring in the Balkan and in Spain ranged from 4.93 to 6.0 (uncorrected p in %) whereas divergences between M. avidus A and M. avidus B were 0.26 to 1.56. Divergence among morphologically identified individuals of M. avidus A and M. avidus B species ranged from 0.13 to 1.58, and from 0.13 to 0.52, respectively. The phenotypic substructuring and observed genetic uniqueness of populations in spatially and temporally fragmented M. avidus taxa were used to identify genetic units. The early split of two allopatric lineages, Spanish M. bicolor and Balkan M. avidus, was followed by diversification in each lineage. Present‐day morphological uniformity masks much of the genetic complexity of lineages within the M. avidus complex. © 2009 The Linnean Society of London, Zoological Journal of the Linnean Society, 2009, 155 , 819–833.     

Amino acids in transmembrane helix 1 confer major functional differences between human and mouse orthologs of the polyspecific membrane transporter OCT1

Organic cation transporter 1 (OCT1) is a membrane transporter that affects hepatic uptake of cationic and weakly basic drugs. OCT1 transports structurally highly diverse substrates. The mechanisms conferring this polyspecificity are unknown. Here, we analyzed differences in transport kinetics between human and mouse OCT1 orthologs to identify amino acids that contribute to the polyspecificity of OCT1. Following stable transfection of HEK293 cells, we observed more than twofold differences in the transport kinetics of 22 out of 28 tested substrates. We found that the β2-adrenergic drug fenoterol was transported with eightfold higher affinity but at ninefold lower capacity by human OCT1. In contrast, the anticholinergic drug trospium was transported with 11-fold higher affinity but at ninefold lower capacity by mouse Oct1. Using human–mouse chimeric constructs and site-directed mutagenesis, we identified nonconserved amino acids Cys36 and Phe32 as responsible for the species-specific differences in fenoterol and trospium uptake. Substitution of Cys36 (human) to Tyr36 (mouse) caused a reversal of the affinity and capacity of fenoterol but not trospium uptake. Substitution of Phe32 to Leu32 caused reversal of trospium but not fenoterol uptake kinetics. Comparison of the uptake of structurally similar β2-adrenergics and molecular docking analyses indicated the second phenol ring, 3.3 to 4.8 Å from the protonated amino group, as essential for the affinity for fenoterol conferred by Cys36. This is the first study to report single amino acids as determinants of OCT1 polyspecificity. Our findings suggest that structure–function data of OCT1 is not directly transferrable between substrates or species.     

Predicting the spread of Aedes albopictus in Australia under current and future climates: Multiple approaches and datasets to incorporate potential evolutionary divergence

When predicting the potential and future invasive range of a species, there is a growing appreciation that insights about factors limiting distributions can be provided by using multiple modelling approaches and by incorporating information from different parts of a species range. Here we apply this strategy to build on previous CLIMEX models to predict the invasion potential of Aedes albopictus, the Asian tiger mosquito, in mainland Australia. A combination of CLIMEX and MAXENT modelling indicated that the mosquito was expected to become widespread along the eastern seaboard and extend into northern Tasmania, but to remain restricted to the coastal fringe, a pattern which is not expected to shift much under climate change. However, a recent expansion of A. albopictus in North America points to evolutionary changes affecting the distribution of this species; when the North American range is included in models, A. albopictus is predicted to become much more widespread and extend inland and into Western Australia. These patterns highlight the potential impact of evolution on species distributions arising from multiple introductions or in situ evolution. By considering future climate scenarios, we demonstrate that there is likely to be a persistent public health threat associated with invasion by this species.     

Redox atlas of the mouse : Immunohistochemical detection of glutaredoxin-, peroxiredoxin-, and thioredoxin-family proteins in various tissues of the laboratory mouse

Background

Oxidoreductases of the thioredoxin family of proteins have been thoroughly studied in numerous cellular and animal models mimicking human diseases. Despite of their well documented role in various disease conditions, no systematic information on the presence of these proteins is available.

Methods

Here, we have systematically analyzed the presence of some of the major constituents of the glutaredoxin (Grx)-, peroxiredoxin (Prx)-, and thioredoxin (Trx)-systems, i.e. Grx1, Grx2, Grx3 (TXNL-2/PICOT), Grx5, nucleoredoxin (Nrx), Prx1, Prx2, Prx3, Prx4, Prx5, Prx6, Trx1, thioredoxin reductase 1 (TrxR1), Trx2, TrxR2, and γ-glutamyl cysteine synthetase (γ-GCS) in various tissues of the mouse using immunohistochemistry.

Results

The identification of the Trx family proteins in the central nervous system, sensory organs, digestive system, lymphatic system, reproductive system, urinary system, respiratory system, endocrine system, skin, heart, and muscle revealed a number of significant differences between these proteins with respect to their distribution in these tissues.

Conclusion

Our results imply more specific functions and interactions between the proteins of this family than previously assumed.

General significance

Crucial functions of Trx family proteins have been demonstrated in various disease conditions. A detailed overview on their distribution in various tissues will be helpful to fully comprehend their potential role and the interactions of these proteins in the most thoroughly studied model for human diseases—the laboratory mouse.This article is part of a Special Issue entitled Human and Murine Redox Protein Atlases.     

Identification,expression and enzyme activity of the group III sPLA2s in Cyprinus carpio L

Five group III secreted phospholipase (pla2g3s) homologous genes located on different linkage groups were identified from common carp (Cyprinus carpio), which we named Ccpla2g3a1, Ccpla2g3a2, Ccpla2g3b, Ccpla2g3c1 and Ccpla2g3c2. The five genes encode 530, 525, 461, 752 and 753 amino acids, respectively. Sequence analysis showed that the Ccpla2g3as contain seven exons and the others contain four exons. Synteny analysis of fish pla2g3s indicated that pla2g3a and pla2g3b were from the same ancestor gene, and Ccpla2g3a1, Ccpla2g3a2, Ccpla2g3c1 and Ccpla2g3c2 were from the specific genome duplication of common carp. Due to the significant variation of the pla2g3bs from common carp and zebrafish (Danio rerio), they formed a separate group in the phylogenetic tree. The tissue distributions of Ccpla2g3s coincided with their expression profiles during the embryo stages. The expression levels of Ccpla2g3as and Ccpla2g3cs were low at the embryo stages, and they were abundant in the liver and brain, respectively, whereas the expression of Ccpla2g3b was high at 0.5 h after fertilization and in the ovary. We obtained three soluble recombinant proteins of the bee venom-like PLA2 (BVLP) from Ccpla2g3 and evaluated their PLA₂ enzyme properties. The optimum pHs of MBP-a1-BVLP, MBP-b-BVLP and MBP-c1-BVLP were 7.5, 7.0 and 8.0, respectively, and specific activities were 7.68 ± 0.66, 4.155 ± 0.158 and 1.93 ± 0.05 U μmol^–1, respectively. The K_d for Ca²⁺ of MBP-b-BVLP was the lowest (2.6 μM), whereas the values for both MBP-a1-BVLP and MBP-c1-BVLP were about 15 μM. The K_m values of three proteins ranged from 31.9 to 41.91 μM.     

Spatial, seasonal and species variations of harmful algal blooms in the South Yellow Sea and East China Sea

The occurrences of harmful algal blooms (HABs), in terms of frequency and area in the Chinese coastal waters, have been increasing since 1980s and caused considerable economic losses. In the present study, we have analyzed spatial and seasonal characteristics of HAB events in the southern Yellow Sea and East China Sea along Chinese coast from 1933 to 2004. With a total 435 HAB records, the most frequent HAB occurrence area (FHA) is off the Yangtze River mouth and another two FHA areas are located south of the Yangtze River estuary along about isobaths of 30–60 m coastal water in the East China Sea. The time of HAB occurrence shifted during our study period: from autumn (August–October) before 1980s to July–August in 1980s, during May–July in 1990s, and May–June for the period of 2000–2004. Causative species were found to be different: Noctiluca scintillans and Skeletonema costatum were dominant causative species prior to 2000; and Prorocentrum donghaiense Lu was dominant from 2000 to 2004 and also caused large blooms in May. Trichodesmium sp. caused many HABs in autumn (August–October) prior to 1980s with only one HAB between 1980 and 2004. The changes of the dominant HAB species may have affected the timings of HAB occurrence, as well as the increasing HAB-affected areas in recent years.     

Variation in life-history characteristics among populations of North American wood turtles: a view from the north

Life-history traits such as age at maturity, body size and clutch size tend to vary across a species' distribution. The purpose of our study was to describe the demography of a newly discovered population of North American wood turtles Glyptemys insculpta at the species' northern range limit, and to compare our findings to those of other studies to test hypotheses about adaptive life-history variation. Turtles were hand-captured from May to October 2005 and 2006 along a 4.5 km stretch of river located in the Sudbury District, ON, Canada (46°N). Fifty-five captured individuals provided a population density estimate of 1.3 turtles/100 m of river. Juveniles comprised 35% of wood turtles captured, and growth ring counts (i.e. age estimates) indicated recruitment in each of the past 11 years. Among populations, we found a nonlinear pattern in body size variation with the largest turtles in the north, smallest turtles in the centre of the range, and intermediate-sized turtles in the south. This nonlinear pattern in body size was reflected in clutch size variation. Selective pressures to overcome years of low recruitment may have resulted in larger body sizes and hence large clutch sizes at northern latitudes while conspecifics at southern latitudes can achieve larger body sizes because they live in a more productive environment. Population density decreased with latitude, likely as a result of a gradient in habitat productivity.     

The role of niche divergence and phenotypic adaptation in promoting lineage diversification in the Sage Sparrow (Artemisiospiza belli,Aves: Emberizidae)

Niche divergence or conservatism and phenotypic adaptation are important in lineage diversification. We used mitochondrial DNA (mtDNA), morphology and ecological niche models to examine these processes in three subspecies of Sage Sparrow (Artemisiospiza belli belli, A. b. canescens and A. b. nevadensis) that breed in bioclimatically diverse ecoregions in western North America. Overall, mtDNA and morphology are congruent with subspecies, ecoregion and bioclimatic niche. Niche divergence, rather than niche conservatism, accompanied by phenotypic adaptation, is associated with lineage diversification between subspecies. This diversification has occurred with and without physical barriers or accompanying genetic divergence. Populations of A. b. canescens are divided by a montane barrier into two bioclimatic regions (San Joaquin Valley, Mojave Desert), where they are indistinguishable phenotypically, but show distinctive genetic patterns. Although there is no physical barrier between A. b. canescens in the San Joaquin Valley and A. b. belli in the Coast Ranges, these populations occupy different bioclimatic niches and are phenotypically, but not genetically, diagnosable. Niche overlap is greatest between A. b. canescens from the Mojave Desert and A. b. nevadensis from the Great Basin, yet these subspecies maintain distinctive phenotypes and mtDNA, even in local secondary contact and sympatry. Palaeoclimatic niche models for the Last Glacial Maximum (c. 21 000 bp ) and the Last Interglacial (c. 120 000 bp ) suggest that ecoregionally distinct populations of Artemisiospiza belli experienced different Pleistocene range fluctuations and glacial refugia, with temporal niche conservatism. Populations probably reached their current distributions as favourable climates and habitats expanded after the last glaciation. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, ??, ??–??.     

Phylogeographic structure in long‐tailed voles (Rodentia: Arvicolinae) belies the complex Pleistocene history of isolation,divergence, and recolonization of Northwest North America's fauna

Quaternary climate fluctuations restructured biodiversity across North American high latitudes through repeated episodes of range contraction, population isolation and divergence, and subsequent expansion. Identifying how species responded to changing environmental conditions not only allows us to explore the mode and tempo of evolution in northern taxa, but also provides a basis for forecasting future biotic response across the highly variable topography of western North America. Using a multilocus approach under a Bayesian coalescent framework, we investigated the phylogeography of a wide‐ranging mammal, the long‐tailed vole, Microtus longicaudus. We focused on populations along the North Pacific Coast to refine our understanding of diversification by exploring the potentially compounding roles of multiple glacial refugia and more recent fragmentation of an extensive coastal archipelago. Through a combination of genetic data and species distribution models (SDMs), we found that historical climate variability influenced contemporary genetic structure, with multiple isolated locations of persistence (refugia) producing multiple divergent lineages (Beringian or northern, southeast Alaska or coastal, and southern or continental) during glacial advances. These vole lineages all occur along the North Pacific Coast where the confluence of numerous independent lineages in other species has produced overlapping zones of secondary contact, collectively a suture zone. Finally, we detected high levels of neoendemism due to complex island geography that developed in the last 10,000 years with the rising sea levels of the Holocene.