A vertebrate reproductive system involving three ploidy levels : hybrid origin of triploids in a contact zone of diploid and tetraploid green toads (Bufo viridis subgroup)

The rise and consequences of polyploidy in vertebrates, whose origin was associated with genome duplications, may be best studied in natural diploid and polyploid populations. In a diploid/tetraploid (2n/4n) geographic contact zone of Palearctic green toads in northern Kyrgyzstan, we examine 4ns and triploids (3n) of unknown genetic composition and origins. Using mitochondrial and nuclear sequence, and nuclear microsatellite markers in 84 individuals, we show that 4n (Bufo pewzowi) are allopolyploids, with a geographically proximate 2n species (B. turanensis) being their maternal ancestor and their paternal ancestor as yet unidentified. Local 3n forms arise through hybridization. Adult 3n mature males (B. turanensis mtDNA) have 2n mothers and 4n fathers, but seem distinguishable by nuclear profiles from partly aneuploid 3n tadpoles (with B. pewzowi mtDNA). These observations suggest multiple pathways to the formation of triploids in the contact zone, involving both reciprocal origins. To explain the phenomena in the system, we favor a hypothesis where 3n males (with B. turanensis mtDNA) backcross with 4n and 2n females. Together with previous studies of a separately evolved, sexually reproducing 3n lineage, these observations reveal complex reproductive interactions among toads of different ploidy levels and multiple pathways to the evolution of polyploid lineages.

Genome size rather than content might affect call properties in toads of three ploidy levels (Anura: Bufonidae: Bufo viridis subgroup).

In vertebrates, genome size has been shown to correlate with nuclear and cell sizes, and influences phenotypic features, such as brain complexity. In three different anuran families, advertisement calls of polyploids exhibit longer notes and intervals than diploids, and difference in cellular dimensions have been hypothesized to cause these modifications. We investigated this phenomenon in green toads (Bufo viridis subgroup) of three ploidy levels, in a different call type (release calls) that may evolve independently from advertisement calls, examining 1205 calls, from ten species, subspecies, and hybrid forms. Significant differences between pulse rates of six diploid and four polyploid (3n, 4n) green toad forms across a range of temperatures from 7 to 27 °C were found. Laboratory data supported differences in pulse rates of triploids vs. tetraploids, but failed to reach significance when including field recordings. This study supports the idea that genome size, irrespective of call type, phylogenetic context, and geographical background, might affect call properties in anurans and suggests a common principle governing this relationship. The nuclear-cell size ratio, affected by genome size, seems the most plausible explanation. However, we cannot rule out hypotheses under which call-influencing genes from an unexamined diploid ancestral species might also affect call properties in the hybrid-origin polyploids.

Low rates of X-Y recombination, not turnovers, account for homomorphic sex chromosomes in several diploid species of Palearctic green toads (Bufo viridis subgroup).

Contrasting with birds and mammals, most ectothermic vertebrates present homomorphic sex chromosomes, which might be due either to a high turnover rate or to occasional X-Y recombination. We tested these two hypotheses in a group of Palearctic green toads that diverged some 3.3 million years ago. Using sibship analyses of sex-linked markers, we show that all four species investigated share the same pair of sex chromosomes and a pattern of male heterogamety with drastically reduced X-Y recombination in males. Phylogenetic analyses of sex-linked sequences show that X and Y alleles cluster by species, not by gametolog. We conclude that X-Y homomorphy and fine-scale sequence similarity in these species do not stem from recent sex-chromosome turnovers, but from occasional X-Y recombination.

Strong reproductive barriers in a narrow hybrid zone of West-Mediterranean green toads (Bufo viridis subgroup) with Plio-Pleistocene divergence

Background One key question in evolutionary biology deals with the mode and rate at which reproductive isolation accumulates during allopatric speciation. Little is known about secondary contacts of recently diverged anuran species. Here we conduct a multi-locus field study to investigate a contact zone between two lineages of green toads with an estimated divergence time of 2.7 My, and report results from preliminary experimental crosses. Results The Sicilian endemic Bufo siculus and the Italian mainland-origin B. balearicus form a narrow hybrid zone east of Mt. Etna. Despite bidirectional mtDNA introgression over a ca. 40 km North-South cline, no F1 hybrids could be found, and nuclear genomes display almost no admixture. Populations from each side of the contact zone showed depressed genetic diversity and very strong differentiation (FST = 0.52). Preliminary experimental crosses point to a slightly reduced fitness in F1 hybrids, a strong hybrid breakdown in backcrossed offspring (F1 x parental, with very few reaching metamorphosis) and a complete and early mortality in F2 (F1 x F1). Conclusion Genetic patterns at the contact zone are molded by drift and selection. Local effective sizes are reduced by the geography and history of the contact zone, B. balearicus populations being at the front wave of a recent expansion (late Pleistocene). Selection against hybrids likely results from intrinsic genomic causes (disruption of coadapted sets of genes in backcrosses and F2-hybrids), possibly reinforced by local adaptation (the ranges of the two taxa roughly coincide with the borders of semiarid and arid climates). The absence of F1 in the field might be due to premating isolation mechanisms. Our results, show that these lineages have evolved almost complete reproductive isolation after some 2.7 My of divergence, contrasting sharply with evidence from laboratory experiments that some anuran species may still produce viable F1 offspring after > 20 My of divergence.

Inferring the degree of incipient speciation in secondary contact zones of closely related lineages of Palearctic green toads (Bufo viridis subgroup).

Reproductive isolation between lineages is expected to accumulate with divergence time, but the time taken to speciate may strongly vary between different groups of organisms. In anuran amphibians, laboratory crosses can still produce viable hybrid offspring >20 My after separation, but the speed of speciation in closely related anuran lineages under natural conditions is poorly studied. Palearctic green toads (Bufo viridis subgroup) offer an excellent system to address this question, comprising several lineages that arose at different times and form secondary contact zones. Using mitochondrial and nuclear markers, we previously demonstrated that in Sicily, B. siculus and B. balearicus developed advanced reproductive isolation after Plio-Pleistocene divergence (2.6 My, 3.3-1.9), with limited historic mtDNA introgression, scarce nuclear admixture, but low, if any, current gene flow. Here, we study genetic interactions between younger lineages of early Pleistocene divergence (1.9 My, 2.5-1.3) in northeastern Italy (B. balearicus, B. viridis). We find significantly more, asymmetric nuclear and wider, differential mtDNA introgression. The population structure seems to be molded by geographic distance and barriers (rivers), much more than by intrinsic genomic incompatibilities. These differences of hybridization between zones may be partly explained by differences in the duration of previous isolation. Scattered research on other anurans suggests that wide hybrid zones with strong introgression may develop when secondary contacts occur <2 My after divergence, whereas narrower zones with restricted gene flow form when divergence exceeds 3 My. Our study strengthens support for this rule of thumb by comparing lineages with different divergence times within the same radiation.

New polymorphic microsatellites markers and development of mitotyping primers for West-Mediterranean green toad species (Bufo viridis subgroup)

We report new polymorphic microsatellites for three species of Palearctic green toads (Bufo viridis subgroup): 10 in B. balearicus and seven each in B. siculus and B. boulengeri. Diversity at these loci, measured for 27 B. balearicus, 23 B. siculus and 11 B. boulengeri, ranged from low to high (two to 10 alleles). Mitotyping primers, specific to the control region, which allow fast screening of parapatric Sicilian endemic B. siculus and Italian mainland-origin B. balearicus, were developed.

Thirteen polymorphic microsatellite markers for the European green toad Bufo viridis viridis, a declining amphibian species.

We report 13 new polymorphic microsatellite markers for the European green toad Bufo viridis viridis (B. viridis subgroup), a declining amphibian from Central, Southeastern and Eastern Europe. Diversity at these loci estimated for 19 individuals ranged from two to ten alleles. Most of these primers also cross-amplify in related West-Mediterranean green toad species (Bufo balearicus, B. siculus and B. boulengeri). These microsatellites will be useful for conservation genetics of threatened Bufo viridis viridis populations and evolutionary studies of green toad taxa in secondary contact to examine hybridization.

A genetic isolation gradient of populations of the Balearic green toad (Bufo balearicus) follows rising eastward fragmentation of the rural landscapes on the island of Menorca.

We present the first approach to the genetic diversity and structure of the Balearic toad (Bufo balearicus Boettger, 1880) for the island of Menorca. Forty-one individ- uals from 21 localities were analyzed for ten microsatellite loci. We used geo-refer- enced individual multilocus genotypes and a model-based clustering method for the inference of the number of populations and of the spatial location of genetic dis- continuities between those populations.¦Only six of the microsatellites analyzed were polymorphic. We revealed a northwest- ern area inhabited by a single population with several well-connected localities and another set of populations in the southeast that includes a few unconnected small units with genetically significant differences among them as well as with the individ- uals from the northwest of the island. The observed fragmentation may be explained by shifts from agricultural to tourism practices that have been taking place on the island of Menorca since the 1960s. The abandonment of rural activities in favor of urbanization and concomitant service areas has mostly affected the southeast of the island and is currently threatening the overall geographic connectivity between the different farming areas of the island that are inhabited by the Balearic toad.

Development and cross-amplification of thirty microsatellite loci in five diploid and polyploid Central Asian species of Palearctic green toads (Bufo viridis subgroup)

We report 30 polymorphic microsatellite markers for five species of Palearctic green toads (Bufo viridis subgroup): 23 in the diploid B. latastii, 19 in diploid B. turanensis, 20 in diploid B. shaartusiensis, 27 in tetraploid B. pewzowi, and 30 in triploid B. baturae. Genetic diversity at these loci, measured for 10-18 individuals per species, ranged from 2 to 19 alleles. These microsatellite loci will be useful for conservation plans (genetic diversity, population structure, evolutionary units), inheritance patterns, and evolution of green toads.

The fourth transmembrane segment of the Na,K-ATPase alpha subunit: a systematic mutagenesis study.

The Na,K-ATPase is a major ion-motive ATPase of the P-type family responsible for many aspects of cellular homeostasis. To determine the structure of the pathway for cations across the transmembrane portion of the Na,K-ATPase, we mutated 24 residues of the fourth transmembrane segment into cysteine and studied their function and accessibility by exposure to the sulfhydryl reagent 2-aminoethyl-methanethiosulfonate. Accessibility was also examined after treatment with palytoxin, which transforms the Na,K-pump into a cation channel. Of the 24 tested cysteine mutants, seven had no or a much reduced transport function. In particular cysteine mutants of the highly conserved "PEG" motif had a strongly reduced activity. However, most of the non-functional mutants could still be transformed by palytoxin as well as all of the functional mutants. Accessibility, determined as a 2-aminoethyl-methanethiosulfonate-induced reduction of the transport activity or as inhibition of the membrane conductance after palytoxin treatment, was observed for the following positions: Phe(323), Ile(322), Gly(326), Ala(330), Pro(333), Glu(334), and Gly(335). In accordance with a structural model of the Na,K-ATPase obtained by homology modeling with the two published structures of sarcoplasmic and endoplasmic reticulum calcium ATPase (Protein Data Bank codes 1EUL and 1IWO), the results suggest the presence of a cation pathway along the side of the fourth transmembrane segment that faces the space between transmembrane segments 5 and 6. The phenylalanine residue in position 323 has a critical position at the outer mouth of the cation pathway. The residues thought to form the cation binding site II ((333)PEGL) are also part of the accessible wall of the cation pathway opened by palytoxin through the Na,K-pump.

Functional effects of Na+,K+-ATPase gene mutations linked to familial hemiplegic migraine.

Familial hemiplegic migraine type 2, an autosomal dominant form of migraine with aura, has been associated with four distinct mutations in the alpha2-subunit of the Na+,K+-ATPase. We have introduced these mutations in the alpha2-subunit of the human Na+,K+-ATPase and the corresponding mutations in the Bufo marinus alpha1-subunit and studied these mutants by expression in Xenopus oocyte. Metabolic labeling studies showed that the mutants were synthesized and associated with the beta-subunit, except for the alpha2HW887R mutant, which was poorly synthesized, and the alpha1BW890R, which was partially retained in the endoplasmic reticulum. [3H]ouabain binding showed the presence of the alpha2HR689Q and alpha2HM731T at the membrane, whereas the alpha2HL764P and alpha2HW887R could not be detected. Functional studies with the mutants of the B. marinus Na+,K+-ATPase showed a reduced or abolished electrogenic activity and a low K+ affinity for the alpha1BW890R mutant. Through different mechanisms, all these mutations result in a strong decrease of the functional expression of the Na+,K+-pump. The decreased activity in alpha2 isoform of the Na+,K+-pump expressed in astrocytes seems an essential component of hemiplegic migraine pathogenesis and may be responsible for the cortical spreading depression, which is one of the first events in migraine attacks.

Adrenergic, dopaminergic, and muscarinic receptor stimulation leads to PKA phosphorylation of Na-K-ATPase.

Na-K-adenosinetriphosphatase (Na-K-ATPase) is a potential target for phosphorylation by protein kinase A (PKA) and C (PKC). We have investigated whether the Na-K-ATPase alpha-subunit becomes phosphorylated at its PKA or PKC phosphorylation sites upon stimulation of G protein-coupled receptors primarily linked either to the PKA or the PKC pathway. COS-7 cells, transiently or stably expressing Bufo marinus Na-K-ATPase wild-type alpha- or mutant alpha-subunits affected in its PKA or PKC phosphorylation site, were transfected with recombinant DNA encoding beta 2- or alpha 1-adrenergic (AR), dopaminergic (D1A-R), or muscarinic cholinergic (M1-AChR) receptor subspecies. Agonist stimulation of beta 2-AR or D1A-R led to phosphorylation of the wild-type alpha-subunit, as well as the PKC mutant, but not of the PKA mutant, indicating that these receptors can phosphorylate the Na-K-ATPase via PKA activation. Surprisingly, stimulation of the alpha 1B-AR, alpha 1C-AR, and M1-AChR also increased the phosphorylation of the wild-type alpha-subunit and its PKC mutant but not of its PKA mutant. Thus the phosphorylation induced by these primarily phospholipase C-linked receptors seems mainly mediated by PKA activation. These data indicate that the Na-K-ATPase alpha-subunit can act as an ultimate target for PKA phosphorylation in a cascade starting with agonist-receptor interaction and leading finally to a phosphorylation-mediated regulation of the enzyme.

Hydrophobic labeling of (Na+,K+)-ATPase: further evidence that the beta subunit is embedded in the membrane bilayer.

O-Hexanoyl-3,5-diiodo-N-(4-azido-2-nitro-phenyl)tyramine has been used after photochemical conversion into the reactive nitrene to label (Na+,K+)-ATPase from Bufo marinus toad kidney. Immunochemical evidence indicates that the reagent labels both subunits of the enzyme in partially purified form as well as in microsomal membranes. These results support the view that the glycoprotein subunit, like the catalytic subunit, possesses hydrophobic domains by which it is integrated into the plasma membrane.

Influence of environmental conditions on the distribution of Central Asian green toads with three ploidy levels.

We studied the distribution of Palearctic green toads (Bufo viridis subgroup), an anuran species group with three ploidy levels, inhabiting the Central Asian Amudarya River drainage. Various approaches (one-way, multivariate, components variance analyses and maximum entropy modelling) were used to estimate the effect of altitude, precipitation, temperature and land vegetation covers on the distribution of toads. It is usually assumed that polyploid species occur in regions with harsher climatic conditions (higher latitudes, elevations, etc.), but for the green toads complex, we revealed a more intricate situation. The diploid species (Bufo shaartusiensis and Bufo turanensis) inhabit the arid lowlands (from 44 to 789 m a.s.l.), while tetraploid Bufo pewzowi were recorded in mountainous regions (340-3492 m a.s.l.) with usually lower temperatures and higher precipitation rates than in the region inhabited by diploid species. The triploid species Bufo baturae was found in the Pamirs (Tajikistan) at the highest altitudes (2503-3859 m a.s.l.) under the harshest climatic conditions.

Phosphorylation of the Na,K-ATPase alpha-subunit by protein kinase A and C in vitro and in intact cells. Identification of a novel motif for PKC-mediated phosphorylation.

Na,K-ATPase is a potential target for regulatory phosphorylation by protein kinase A and C (PKA and PKC). To identify the phosphorylation sites, we have mutated the alpha 1-subunit of Bufo marinus in a highly conservative PKA and in 20 different PKC consensus sequences. The mutants were expressed in Xenopus oocytes and their phosphorylation capacity tested in homogenates upon stimulation of PKA or PKC. While serine 943 (Ser-943) was identified as a unique target site for PKA, none of the PKC consensus serine or threonine residues are implicated in PKC phosphorylation. Controlled trypsinolysis of phosphorylated alpha-subunits of various purified enzyme preparations and of alpha/beta complexes from oocyte homogenates revealed that PKC phosphorylation was exclusively associated with the N terminus. A fusion protein containing the first 32 amino acids of the Bufo alpha-subunit was phosphorylated in vitro and serine and threonine residues (Thr-15 and Ser-16) in this region were identified by site-directed mutagenesis as the PKC phosphorylation sites. Finally, the Bufo alpha-subunit was phosphorylated by protein kinases in transfected COS-7 cells. In intact cells, PKA stimulation induced phosphorylation exclusively on Ser-943 and PKC stimulation mainly on Thr-15 and Ser-16, which are contained in a novel PKC phosphorylation motif.