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.

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.

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.

Host-parasite relationships during a biologic invasion: 75 years postinvasion, cane toads and sympatric Australian frogs retain separate lungworm faunas.

Invasive species may carry with them parasites from their native range, differing from parasite taxa found in the invaded range. Host switching by parasites (either from the invader to native fauna or from native fauna to the invader) may have important consequences for the viability of either type of host (e.g., their survivorship, fecundity, dispersal ability, or geographic distribution). Rhabdias pseudosphaerocephala (Nematoda) is a common parasite of cane toads (Rhinella marina) in the toad's native range (South and Central America) and also in its introduced Australian range. This lungworm can depress host viability and is capable of infecting Australian frogs in laboratory trials. Despite syntopy between toads and frogs for up to 75 yr, our analyses, based on DNA sequence data of lungworms from 80 frogs and 56 toads, collected from 2008 to 2011, did not reveal any cases of host switching in nature: toads and native frogs retain entirely different lungworm faunas. All lungworms in cane toads were the South and Central American species Rhabdias pseudosphaerocephala, whereas Australian frogs contained at least four taxa (mostly undescribed and currently lumped under the name Rhabdias cf. hylae). General patterns of prevalence and intensity, based on the dissection of 1,315 frogs collected between 1989 and 2011 across the toads' Australian range, show that these Australian endemic Rhabdias spp. are widely distributed geographically and across host taxa but are more common in some frog species (especially, large-bodied species) than they are in others.

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.

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.

Invader immunology: invasion history alters immune system function in cane toads (Rhinella marina) in tropical Australia.

Because an individual's investment into the immune system may modify its dispersal rate, immune function may evolve rapidly in an invader. We collected cane toads (Rhinella marina) from sites spanning their 75-year invasion history in Australia, bred them, and raised their progeny in standard conditions. Evolved shifts in immune function should manifest as differences in immune responses among the progeny of parents collected in different locations. Parental location did not affect the offspring's cell-mediated immune response or stress response, but blood from the offspring of invasion-front toads had more neutrophils, and was more effective at phagocytosis and killing bacteria. These latter measures of immune function are negatively correlated with rate of dispersal in free-ranging toads. Our results suggest that the invasion of tropical Australia by cane toads has resulted in rapid genetically based compensatory shifts in the aspects of immune responses that are most compromised by the rigours of long-distance dispersal.

A cryptic heterogametic transition revealed by sex-linked DNA markers in Palearctic green toads.

In sharp contrast to birds and mammals, most cold-blooded vertebrates have homomorphic (morphologically undifferentiated) sex chromosomes. This might result either from recurrent X-Y recombination (occurring e.g. during occasional events of sex reversal) or from frequent turnovers (during which sex-determining genes are overthrown by new autosomal mutations). Evidence for turnovers is indeed mounting in fish, but very few have so far been documented in amphibians, possibly because of practical difficulties in identifying sex chromosomes. Female heterogamety (ZW) has long been established in Bufo bufo, based on sex reversal and crossing experiments. Here, we investigate a sex-linked marker identified from a laboratory cross between Palearctic green toads (Bufo viridis subgroup). The F(1) offspring produced by a female Bufo balearicus and a male Bufo siculus were phenotypically sexed, displaying an even sex ratio. A sex-specific marker detected in highly reproducible AFLP genotypes was cloned. Sequencing revealed a noncoding, microsatellite-containing fragment. Reamplification and genotyping of families of this and a reciprocal cross showed B. siculus to be male heterogametic (XY) and suggested the same system for B. balearicus. Our results thus reveal a cryptic heterogametic transition within bufonid frogs and help explain patterns of hybrid fitness within the B. viridis subgroup. Turnovers of genetic sex-determination systems may be more frequent in amphibians than previously thought and thus contribute to the prevalence of homomorphic sex chromosomes in this group.

Origin and genome evolution of polyploid green toads in Central Asia: evidence from microsatellite markers.

Polyploidization, which is expected to trigger major genomic reorganizations, occurs much less commonly in animals than in plants, possibly because of constraints imposed by sex-determination systems. We investigated the origins and consequences of allopolyploidization in Palearctic green toads (Bufo viridis subgroup) from Central Asia, with three ploidy levels and different modes of genome transmission (sexual versus clonal), to (i) establish a topology for the reticulate phylogeny in a species-rich radiation involving several closely related lineages and (ii) explore processes of genomic reorganization that may follow polyploidization. Sibship analyses based on 30 cross-amplifying microsatellite markers substantiated the maternal origins and revealed the paternal origins and relationships of subgenomes in allopolyploids. Analyses of the synteny of linkage groups identified three markers affected by translocation events, which occurred only within the paternally inherited subgenomes of allopolyploid toads and exclusively affected the linkage group that determines sex in several diploid species of the green toad radiation. Recombination rates did not differ between diploid and polyploid toad species, and were overall much reduced in males, independent of linkage group and ploidy levels. Clonally transmitted subgenomes in allotriploid toads provided support for strong genetic drift, presumably resulting from recombination arrest. The Palearctic green toad radiation seems to offer unique opportunities to investigate the consequences of polyploidization and clonal transmission on the dynamics of genomes in vertebrates.

Associated anomalies in multi-malformed infants with cleft lip and palate: An epidemiologic study of nearly 6 million births in 23 EUROCAT registries.

We studied 5,449 cases of cleft lip (CL) with or without cleft palate (CL/P) identified between 1980 and 2000 from the EUROCAT network of 23 registers (nearly 6 million births) in 14 European countries. We investigated specific types of defects associated with clefts. Among CL/P cases (prevalence = 9.1 per 10,000), 1,996 (36.6%) affected only the lip (CL) and 3,453 (63.4%) involved CL and palate (CLP). A total of 3,860 CL/P cases (70.8%) occurred as isolated anomalies and 1,589 (29.2%) were associated with other defects such as multiple congenital anomalies of unknown origin (970), chromosomal (455) and recognized syndromes (164). Associated malformations were more frequent in infants who had CLP (34.0%) than in infants with CL only (20.8%). Among multi-malformed infants, 2 unrelated anomalies were found in 351 cases, 3 in 242 cases, and 4 or more in 377 cases. Among 5,449 CL/P cases, 4,719 were live births (LB) (86.6%), 203 stillbirths (SB) (3.7%), while 508 (9.3%) were terminations of pregnancy (ToP). CL/P occurred significantly more frequently in males (M/F = 1.70), especially among total isolated cases (M/F = 1.87) and CLP isolated cases (M/F = 1.92). The study confirmed that musculoskeletal, cardiovascular, and central nervous system defects are frequently associated with CL/P. An association with reduction anomalies of the brain was found. This association suggests that clinicians should seek to identify structural brain anomalies in these patients with CL/P as the potential functional consequences may be important for rehabilitation and clinical management.

Use of asthma medication during pregnancy and risk of specific congenital anomalies: A European case-malformed control study.

BACKGROUND: Pregnant women with asthma need to take medication during pregnancy. OBJECTIVE: We sought to identify whether there is an increased risk of specific congenital anomalies after exposure to antiasthma medication in the first trimester of pregnancy. METHODS: We performed a population-based case-malformed control study testing signals identified in a literature review. Odds ratios (ORs) of exposure to the main groups of asthma medication were calculated for each of the 10 signal anomalies compared with registrations with nonchromosomal, nonsignal anomalies as control registrations. In addition, exploratory analyses were done for each nonsignal anomaly. The data set included 76,249 registrations of congenital anomalies from 13 EUROmediCAT registries. RESULTS: Cleft palate (OR, 1.63; 95% CI, 1.05-2.52) and gastroschisis (OR, 1.89; 95% CI, 1.12-3.20) had significantly increased odds of exposure to first-trimester use of inhaled β2-agonists compared with nonchromosomal control registrations. Odds of exposure to salbutamol were similar. Nonsignificant ORs of exposure to inhaled β2-agonists were found for spina bifida, cleft lip, anal atresia, severe congenital heart defects in general, or tetralogy of Fallot. None of the 4 literature signals of exposure to inhaled steroids were confirmed (cleft palate, cleft lip, anal atresia, and hypospadias). Exploratory analyses found an association between renal dysplasia and exposure to the combination of long-acting β2-agonists and inhaled corticosteroids (OR, 3.95; 95% CI, 1.99-7.85). CONCLUSIONS: The study confirmed increased odds of first-trimester exposure to inhaled β2-agonists for cleft palate and gastroschisis and found a potential new signal for renal dysplasia associated with combined long-acting β2-agonists and inhaled corticosteroids. Use of inhaled corticosteroids during the first trimester of pregnancy seems to be safe in relation to the risk for a range of specific major congenital anomalies.