Variation in repeat length and heteroplasmy of the mitochondrial DNA control region along a core-edge gradient in the eastern spadefoot toad (Pelobates syriacus)

Peripheral populations are those situated at the distribution margins of a species and are often subjected to more extreme abiotic and biotic conditions than those at the core. Here, we hypothesized that shorter repeat length and fewer heteroplasmic mitochondrial DNA (mtDNA) copies, which are associated with more efficient mitochondrial function, may be related to improved survival under extreme environmental conditions. We sampled eastern spadefoot toads (mostly as tadpoles) from 43 rain pools distributed along a 300-km gradient from core to edge of the species' distribution. We show that mean pool tandem repeat length and heteroplasmy increase from edge to core, even after controlling for body size. We evaluate several alternative hypotheses and propose the Fisher hypothesis as the most likely explanation. However, additional sequential sampling and experimental studies are required to determine whether selection under extreme conditions, or alternative mechanisms, could account for the gradient in heteroplasmy and repeat length in the mtDNA control region.     

Polystoma nacialtuneli n. sp. (Monogenea: Polystomatidae) from the eastern spadefoot, Pelobates syriacus (Pelobatidae) in Turkey

Polystoma nacialtuneli n. sp. is described from the urinary bladder of the eastern spadefoot, Pelobates syriacus from Turkey. This is the fifth polystome species known from Turkey and the third species in Pelobates. We show that this new parasite species can be distinguished from other polystome species in the area by a combination of characteristics, including parasite size and the shape and size of the hamuli. Polystoma pelobatis from Pelobates cultripes has a pair of well-developed hamuli, while P. fuscus from Pe. fuscus characteristically has a pair of underdeveloped hamuli barely larger than the marginal hooklets. Polystoma nacialtuneli n. sp. has well-developed hamuli that vary significantly in shape. Phylogenetic relationships of P. nacialtuneli n. sp. within Polystoma, supplemented with molecular divergences estimated from internal transcribed spacer 1 (ITS1) sequences, indicate that they are well separated from their closest relatives, i.e. P. fuscus and P. pelobatis from Pe. fuscus and Pe. cultripes, respectively.     

Mapping the intramolecular signal transduction of G‐protein coupled receptors

G‐protein coupled receptors (GPCRs), a major gatekeeper of extracellular signals on plasma membrane, are unarguably one of the most important therapeutic targets. Given the recent discoveries of allosteric modulations, an allosteric wiring diagram of intramolecular signal transductions would be of great use to glean the mechanism of receptor regulation. Here, by evaluating betweenness centrality (C_B) of each residue, we calculate maps of information flow in GPCRs and identify key residues for signal transductions and their pathways. Compared with preexisting approaches, the allosteric hotspots that our C_B‐based analysis detects for A_2A adenosine receptor (A_2AAR) and bovine rhodopsin are better correlated with biochemical data. In particular, our analysis outperforms other methods in locating the rotameric microswitches, which are generally deemed critical for mediating orthosteric signaling in class A GPCRs. For A_2AAR, the inter‐residue cross‐correlation map, calculated using equilibrium structural ensemble from molecular dynamics simulations, reveals that strong signals of long‐range transmembrane communications exist only in the agonist‐bound state. A seemingly subtle variation in structure, found in different GPCR subtypes or imparted by agonist bindings or a point mutation at an allosteric site, can lead to a drastic difference in the map of signaling pathways and protein activity. The signaling map of GPCRs provides valuable insights into allosteric modulations as well as reliable identifications of orthosteric signaling pathways. Proteins 2014; 82:727–743. © 2013 Wiley Periodicals, Inc.     

Phylogeographic patterns of genetic diversity in the common spadefoot toad, Pelobates fuscus (Anura: Pelobatidae), reveals evolutionary history, postglacial range expansion and secondary contact

Based on allozyme variation of 410 newly collected individuals from 52 localities, we reconstructed range-wide phylogeography of the widespread Western Palearctic anuran, Pelobates fuscus. To study genetic diversity, evolutionary history, postglacial range expansion and secondary contact zones, we used a multidisciplinary approach combining information from various genetic analyses and ecological niche modeling. We confirmed the presence of two main groups in P. fuscus, initially revealed by genome size variation. Pelobates f. vespertinus presents a monomorphic group, but two main groups can be identified in P. f. fuscus: an East European and a West European group. We suggest the existence of at least four different Last Glacial refugia for P. fuscus: (1) the area between the Caspian and Azov Seas as the origin for the expansion of P. f. vespertinus; (2) the northwestern part of the Black Sea area for the East European P. f. fuscus; (3) the southwestern part of the Pannonian Plain and (4) the Po Plain for the West European P. f. fuscus. The routes of postglacial range expansions from the refugia are discussed. We newly identified a hybrid zone between P. f. fuscus and P. f. vespertinus. The width of this zone is about 12.5 km. In light of these findings, the two subspecies of P. fuscus constitute distinct evolutionary lineages and merit recognition as separate species. Our data do not provide support for the validity of P.f. insubricus. We therefore propose to synonymize this subspecies with P. f. fuscus.     

Genetic variation in two forms of the common spadefoot toad Pelobates fuscus (Pelobatidae, Anura, Amphibia) distinguished by genome size

Genetic differences (presumed 23 loci) between two cryptic forms of Pelobates fuscus, differing in genome size, were examined by means of polyacrylamide-gel electrophoresis. This method allowed to discriminate between these forms. Average genetic distance (DN) between both the forms was equal to 0.311, ranging from 0.055 to 0.563. As a rule, differences within these groups were smaller (0.021-0.388). The data show obvious genetic differentiation between these two cryptic forms of P. fuscus. Differences between these forms and P. syriacus were significantly higher (in average 0.943). Genetic distances in relation to speciation and species concepts are discussed.     

River islands,refugia and genetic structuring in the endemic brown frog Rana kukunoris (Anura,Ranidae) of the Qinghai‐Tibetan Plateau

Frequently, Pleistocene climatic cycling has been found to be the diver of genetic structuring in populations, even in areas that did not have continental ice sheets, such as on the Qinghai‐Tibetan Plateau (QTP). Typically, species distributed on the plateau have been hypothesized to re‐treat to south‐eastern refugia, especially during the Last Glacial Maximum (LGM). We evaluated sequence variation in the mitochondrial DNA gene Cytb and the nuclear DNA gene RAG‐1 in Rana kukunoris, a species endemic to the QTP. Two major lineages, N and S, were identified, and lineage N was further subdivided into N1 and N2. The geographical distribution and genealogical divergences supported the hypothesis of multiple refugia. However, major lineages and sublineages diverged prior to the LGM. Demographical expansion was detected only in lineage S and sublineage N2. Sublineage N1 might have survived several glacial cycles in situ and did not expand after the LGM because of the absence of suitable habitat; it survived in river islands. Genetic analysis and environment modelling suggested that the north‐eastern edge of QTP contained a major refugium for R. kukunoris. From here, lineage S dispersed southwards after the LGM. Two microrefugia in northern Qilian Mountains greatly contributed to current level of intraspecific genetic diversity. These results were found to have important implications for the habitat conservation in Northwest China.     

Fossorial but widespread: the phylogeography of the common spadefoot toad (Pelobates fuscus), and the role of the Po Valley as a major source of genetic variability

Pelobates fuscus is a fossorial amphibian that inhabits much of the European plain areas. To unveil traces of expansion and contraction events of the species' range, we sequenced 702 bp of the mitochondrial cytochrome b gene. To infer the population history we applied phylogeographical methods, such as nested clade phylogeographical analysis (NCPA), and used summary statistics to analyse population structure under a neutral model of evolution. Populations were assigned to different drainage systems and we tested hypotheses of explicit refugial models using information from analysis of molecular variance, nucleotide diversity, effective population size estimation, NCPA, mismatch distribution and Bayesian dating. Coalescent simulations were used as post hoc tests for plausibility of derived or a priori assumed biogeographical hypotheses. Our combination of all approaches enabled the reconstruction of the colonization history and phylogeography of P. fuscus and confirmed a previous assumption of the existence of two major genetic lineages within P. fuscus. Using the Afro-European vicariance of Pelobates cultripes and Pelobates varaldii and applying Bayesian dating we estimated the divergence of these phylogeographical lineages to the Pliocene. We suggest the existence of three different glacial refugia: (i) the area between the Caspian and Black Seas as the origin for the expansion of the 'eastern lineage'; (ii) the Danube system as a centre of diversity for part of the 'western lineage'; (iii) the Po Valley, the largest centre of genetic variability. This fits the hypothesis that climatic fluctuation was a key event for differentiation processes in P. fuscus.     

On glacial refugia,genetic diversity,and microevolutionary processes: deep phylogeographical structure in the endemic newt Lissotriton italicus1

We investigated the phylogeographical patterns of Lissotriton italcus, a newt endemic to the Italian peninsula, aiming to determine why hotspots of intraspecific diversity so ‘hot’. We found two main mitochindrial DNA lineages (net sequence divergence of 6.8% at two fragments of total length of 1897 bp): one restricted to part of the Calabrian peninsula (i.e the southernmost portion of the species range) and the other widespread throughout the rest of the species range. Both lineages, which had a parapatric distribution, showed evidence of further subdivisions, with an overall number of eight terminal haplogroups, most of whose times to the most recent common ancestors were estimated at the Late Pleistocene. Analysis of molecular variance suggested that partitioning populations according to the geographical distribution of these haplogroups can explain 97% of the observed genetic variation. These results suggest that L. italicus underwent repeated cycles of allopatric fragmentation throughout the Pleistocene, and that it likely survived the Late Pleistocene paleoenvironmental changes within eight separate refugia. Thus, the current hotspot of intraspecific diversity of L. italicus (within the Calabrian peninsula) has not been moulded by long‐term stability of large populations but rather by multiple events of allopatric fragmentation and divergence. When compared with the patterns recently identified in other species, these results suggest that the occurrence of phases of allopatric divergence (eventually followed by secondary admixture) could be a common, albeit probably underrated feature in the history of formation of hotspots of intraspecific diversity. © 2011 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, 105 , 42–55.     

Diversification of African tree frogs (genus Leptopelis) in the highlands of Ethiopia

The frog genus Leptopelis is composed of ~50 species that occur across sub‐Saharan Africa. The majority of these frogs are typically arboreal; however, a few species have evolved a fossorial lifestyle. Most species inhabit lowland forests, but a few species have adapted to high elevations. Five species of Leptopelis occupy the Ethiopian highlands and provide a good opportunity to study the evolutionary transition from an arboreal to a fossorial lifestyle, as well as the diversification in this biodiversity hot spot. We sequenced 14 nuclear and three mitochondrial genes, and generated thousands of SNPs from ddRAD sequencing to study the evolutionary relationships of Ethiopian Leptopelis. The five species of highland Leptopelis form a monophyletic group, which diversified during the late Miocene and Pliocene. We found strong population structure in the fossorial species L. gramineus, with levels of genetic differentiation between populations similar to those found between arboreal species. This could indicate that L. gramineus is a complex of cryptic species. We propose that after the original colonization of the Ethiopian highlands by the ancestor of the L. gramineus group, episodes of vicariance fragmented the ancestral populations of this group. We also report the re‐evolution of arboreality in L. susanae, which evolved from a fossorial ancestor, a rare ecological switch in frogs that had previously been reported only once.     

Landscape models for nuclear genetic diversity and genetic structure in white-footed mice (Peromyscus leucopus)

Dramatic changes in the North American landscape over the last 12 000 years haveshaped the genomes of the small mammals, such as the white-footed mouse (Peromyscusleucopus), which currently inhabit the region. However, very recent interactions ofpopulations with each other and the environment are expected to leave the most pronouncedsignature on rapidly evolving nuclear microsatellite loci. We analyzed landscapecharacteristics and microsatellite markers of P. leucopus populations along atransect from southern Ohio to northern Michigan, in order to evaluate hypotheses aboutthe spatial distribution of genetic heterogeneity. Genetic diversity increased to thenorth and was best approximated by a single-variable model based on habitat availabilitywithin a 0.5-km radius of trapping sites. Interpopulation differentiation measured byclustering analysis was highly variable and not significantly related to latitude orhabitat availability. Interpopulation differentiation measured as F_STvalues and chord distance was correlated with the proportion of habitat intervening, butwas best explained by agricultural distance and by latitude. The observed gradients indiversity and interpopulation differentiation were consistent with recent habitatavailability being the major constraint on effective population size in this system, andcontradicted the predictions of both the postglacial expansion and core-peripheryhypotheses.     

The shaping of genetic variation in edge‐of‐range populations under past and future climate change

With rates of climate change exceeding the rate at which many species are able to shift their range or adapt, it is important to understand how future changes are likely to affect biodiversity at all levels of organisation. Understanding past responses and extent of niche conservatism in climatic tolerance can help predict future consequences. We use an integrated approach to determine the genetic consequences of past and future climate changes on a bat species, Plecotus austriacus. Glacial refugia predicted by palaeo‐modelling match those identified from analyses of extant genetic diversity and model‐based inference of demographic history. Former refugial populations currently contain disproportionately high genetic diversity, but niche conservatism, shifts in suitable areas and barriers to migration mean that these hotspots of genetic diversity are under threat from future climate change. Evidence of population decline despite recent northward migration highlights the need to conserve leading‐edge populations for spearheading future range shifts.     

Low genetic diversity in the masked palm civet Paguma larvata (Viverridae)

The masked palm civet is distributed through south-east Asia, China and the Himalayas. Because of its potential role in the severe acute respiratory syndrome (SARS) epidemic, it has become important to gather information on this species, and notably to provide a tool to determine the origin of farm and market animals. For this purpose, we studied the genetic variability and the phylogeographic pattern of the masked palm civet Paguma larvata . First, two portions of mitochondrial genes, cytochrome b and the control region, were sequenced for a total of 76 individuals sampled from China, the Indochinese region and the Sundaic region. Results indicated a low genetic variability and suggested a lack of a phylogeographic structure in this species, which do not allow inferring the geographic origin of samples of unknown origin, although it is possible to distinguish individuals from China and the Sundaic region. This low variation is in contrast to the well-marked morphological differentiation between the populations in the Sundaic and Chinese–Indochinese regions. We also used five microsatellite loci to genotype 149 samples from two wild and four farmed populations in China, where the masked palm civet is farmed and where the SARS coronavirus was isolated. These analyses also showed a reduced variability in Chinese civets and showed that farmed populations did not exhibit a lower genetic diversity than wild populations, suggesting frequent introductions of wild individuals into farms.     

The first find of a spadefoot toad (Anura,Pelobatidae) in the Miocene of Moldova

A skeleton fragment of a spadefoot toad (family Pelobatidae Bonaparte, 1850) from the Middle Miocene (Early Sarmatian) of the Karpov Yar locality (northern Moldova) is described. Since the material is incomplete, it is only determined as Pelobatidae gen. indet. This is the first spadefoot toad from the Miocene of Moldova.     

Genetic footprints of late Quaternary climate change in the diversity of Patagonian‐Fueguian rodents

Species are impacted by climate change at both ecological and evolutionary time scales. Studies in northern continents have provided abundant evidence of dramatic shifts in distributions of species subsequent to the last glacial maximum (LGM), particularly at high latitudes. However, little is known about the history of southern continents, especially at high latitudes. South America is the only continent, other than Antarctica, that extends beyond 40 °S. Genetic studies of a few Patagonian species have provided seemingly conflicting results, indicating either postglacial colonization from restricted glacial refugia or persistence through glacial cycles and in situ differentiation. Using mitochondrial DNA sequences of 14 species of sigmodontine rodents, a major faunal ensemble of Patagonia and Tierra del Fuego, we show that at least nine of these species bear genetic footprints of demographic expansion from single restricted sources. However, timing of demographic expansion precedes the LGM in most of these species. Four species are fragmented phylogeographically within the region. Our results indicate that (i) demographic instability in response to historical climate change has been widespread in the Patagonian‐Fueguian region, and is generally more pronounced at high latitudes in both southern and northern continents; (ii) colonization from lower latitudes is an important component of current Patagonian‐Fueguian diversity; but (iii) in situ differentiation has also contributed to species diversity.     

Plio‐Pleistocene diversification and genetic population structure of an endangered lizard (the Blue Mountains water skink,Eulamprus leuraensis) in south‐eastern Australia

Aim Although climatic fluctuations occurred world‐wide during the Pleistocene, the severity of glacial and drought events – and hence their influence on animal and plant biogeography – differed among regions. Many Holarctic species were forced to warmer‐climate refugia during glacial periods, leaving the genetic signature of recent expansion and gene flow among modern‐day populations. Montane south‐eastern Australia experienced less extreme glaciation, but the effects of drier and colder climatic conditions over this period on biotic distributions, and hence on the present‐day genetic structure of animal and plant populations, are poorly known. Location South‐eastern Australia. Methods The endangered Blue Mountains water skink (Eulamprus leuraensis) is a viviparous lizard known from fewer than 40 isolated small swamps at 560–1060 m elevation in south‐eastern Australia. We conducted molecular phylogenetic, dating and population genetics analyses using the mitochondrial NADH dehydrogenase 4 (ND4) of 224 individuals of E. leuraensis sampled across the species’ distribution. Results Ancient divergences in haplotype groups between lizards from the Blue Mountains and the Newnes Plateau, and strong genetic differences, even between swamps separated by only a few kilometres, suggest that the species has persisted as a series of relatively isolated populations within its current distribution for about a million years. Presumably, habitat patches similar to current‐day swamps persisted throughout glacial–interglacial cycles in this region, allowing the development of high levels of genetic structuring within and among present‐day populations. Main conclusions Our results suggest that less extreme glacial conditions occurred in the Southern Hemisphere compared with the Northern Hemisphere, allowing cold‐adapted species (such as E. leuraensis) to persist in montane areas. However, additional studies are needed before we can assemble a comprehensive view of the impact of Pleistocene climatic variation on the phylogeography of Southern Hemisphere taxa.     

Intraspecific genetic variation in the common midwife toad (Alytes obstetricans): subspecies assignment using mitochondrial and microsatellite markers

The common midwife toad (Alytes obstetricans), widely distributed in the northern half of the Iberian Peninsula and part of Western Europe, is currently subdivided into four subspecies: A. o. obstetricans, A. o. boscai, A. o. pertinax and A. o. almogavarii. However, the delimitation of these subspecies and their ranges are still under discussion because strong discordances have been found between morphological and molecular data, and especially among different genetic markers. Here, we screen a set of novel microsatellite loci and mtDNA sequences of A. obstetricans populations representative of all currently recognized subspecies to investigate the correspondence between genetic groupings inferred from clustering analysis of microsatellite genotypes and the described subspecies and test whether patterns of mtDNA variation are concordant with those genetic clusters. Our results confirm previous expectations of extremely high intraspecific diversity in A. obstetricans in Iberia. Analyses of microsatellite and mtDNA data were concordant in recovering five well‐defined groups, of which three correspond to previously defined subspecies, while the two additional clusters correspond to populations of subspecies A. o. boscai separated by the Douro River. Our results suggest the occurrence of two distinct genetic units within A. o. boscai that likely result from a long independent evolutionary history, thus deserving special attention from a conservation point of view.