Genome Dosage Effect and Hybrid Morphology — The Case of the Hybridogenetic Water Frogs of thePelophylax esculentusComplex

Karyotypes of water frogs from the Pelophylax esculentus complex: results of cross-species chromosomal painting

Amphibian species have the largest genome size enriched with repetitive sequences and relatively similar karyotypes. Moreover, many amphibian species frequently hybridize causing nuclear and mitochondrial genome introgressions. In addition, hybridization in some amphibian species may lead to clonality and polyploidization. All such events were found in water frogs from the genus Pelophylax. Among the species within the genus Pelophylax, P. esculentus complex is the most widely distributed and well-studied. This complex includes two parental species, P. ridibundus and P. lessonae, and their hybrids, P. esculentus, reproducing hemiclonally. Parental species and their hybrids have similar but slightly polymorphic karyotypes, so their precise identification is still required. Here, we have developed a complete set of 13 chromosome painting probes for two parental species allowing the precise identification of all chromosomes. Applying chromosomal painting, we identified homologous chromosomes in both parental species and orthologous chromosomes in their diploid hemiclonal hybrids. Comparative painting did not reveal interchromosomal exchanges between the studied water frog species and their hybrids. Using cross-specific chromosome painting, we detected unequal distribution of the signals along chromosomes suggesting the presence of species-specific tandem repeats. Application of chromosomal paints to the karyotypes of hybrids revealed differences in the intensity of staining for P. ridibundus and P. lessonae chromosomes. Thus, both parental genomes have a divergence in unique sequences. Obtained chromosome probes may serve as a powerful tool to unravel chromosomal evolution in phylogenetically related species, identify individual chromosomes in different cell types, and investigate the elimination of chromosomes in hybrid water frogs.

Genome elimination from the germline cells in diploid and triploid male water frogs Pelophylax esculentus

Hybridogenesis is a hemiclonal reproductive strategy in diploid and triploid hybrids. Our study model is a frog P. esculentus (diploid RL and triploids RLL and RRL), a natural hybrid between P. lessonae (LL) and P. ridibundus (RR). Hybridogenesis relies on elimination of one genome (L or R) from gonocytes (G) in tadpole gonads during prespermatogenesis, but not from spermatogonial stem cells (SSCs) in adults. Here we provide the first comprehensive study of testis morphology combined with chromosome composition in the full spectrum of spermatogenic cells. Using genomic in situ hybridization (GISH) and FISH we determined genomes in metaphase plates and interphase nuclei in Gs and SSCs. We traced genomic composition of SSCs, spermatocytes and spermatozoa in individual adult males that were crossed with females of the parental species and gave progeny. Degenerating gonocytes (24%–39%) and SSCs (18%–20%) led to partial sterility of juvenile and adult gonads. We conclude that elimination and endoreplication not properly completed during prespermatogenesis may be halted when gonocytes become dormant in juveniles. After resumption of mitotic divisions by SSCs in adults, these 20% of cells with successful genome elimination and endoreplication continue spermatogenesis, while in about 80% spermatogenesis is deficient. Majority of abnormal cells are eliminated by cell death, however some of them give rise to aneuploid spermatocytes and spermatozoa which shows that hybridogenesis is a wasteful process.

Maintenance of pure hybridogenetic water frog populations: Genotypic variability in progeny of diploid and triploid parents

An intriguing outcome of hybridisation is the emergence of clonally and hemiclonally reproducing hybrids, that can sustain, reproduce, and lead to the emergence of polyploid forms. However, the maintenance of diploid and polyploid hybrid complexes in natural populations remains unresolved. We selected water frogs from the Pelophylax esculentus complex to study how diploid and triploid hybrids, which reproduce hemiclonally via hybridogenesis, are maintained in natural populations. During gametogenesis in diploid hybrids, one of the parental genomes is eliminated, and the remaining genome is endoreplicated. In triploid hybrids, the single-copy genome is typically eliminated, while genome endoreplication does not occur. To investigate how diploid and triploid hybrid frogs reproduce in populations without parental species, we crossed these hybrid animals from two separate pure hybrid populations located in Poland. Using cytogenetic analysis of tadpoles that emerged from the crosses, we established which gametes were produced by parental hybrids. The majority of hybrid females and hybrid males produced one type of gamete with the P. ridibundus genome. However, in both studied populations, approximately half of the diploid and triploid hybrids simultaneously produced gametes with different genome compositions and ploidy levels, specifically, the P. ridibundus and P. lessonae genomes, as well as diploid gametes with genomes of both parental species. Triploid hybrid males and females mostly produced haploid gametes with the P. lessonae genome; however, gametes with the P. ridibundus genome have also been observed. These results suggest that not all hybrids follow the classical hybridogenetic reproduction program and reveal a significant level of alterations in the gametogenesis pathways. In addition, we found a variable survival rate of particular progeny genotypes when we crossed hybrid females with different males suggesting the important role of postzygotic barriers on the maintenance of pure hybrid systems. We suggest that the observed variability in produced gametes and the different survival rate of the progeny with certain genotypes is crucial for the existence of pure hybrid systems.

Genomic basis for an informed conservation management of Pelophylax water frogs in Luxembourg

Genetic identification methods have become increasingly important for species that are difficult to identify in the field. A case in point is Pelophylax water frogs. While their morphological determination is highly complex, they include species protected under EU law and some that are classified as invasive. Additionally, genetic data can provide insights into their complex breeding systems, which may or may not involve the reproductive dependency of one species on another. Here, we generate baseline data for water frog monitoring in Luxembourg. We applied a countrywide sampling approach and used SNPs generated by ddRAD sequencing to identify individuals and infer the breeding systems present in the country. We found Pelophylax lessonae and P. kl. esculentus throughout Luxembourg, mostly living in syntopy. In general, a reproductive dependency of P. kl. esculentus on P. lessonae (L‐E system) was revealed. Besides this general system, we detected triploid P. kl. esculentus in six ponds. This indicates a modified L‐E system with reproductive dependency of the triploids on the diploid P. kl. esculentus. The invasive P. cf. bedriagae was detected in three ponds in southern Luxembourg, with evidence for hybridization with native water frogs. In addition to the ddRAD data, we tested a simple genetic method for future monitoring based on the MND1 marker. It showed in almost all cases, an identical species identification as the ddRAD data and was successfully applied to DNA extracts from mouth swabs. Combining this method with our baseline data will enable informed choices for the protection of native water frog species in Luxembourg.

Morphometry of two cryptic tree frog species at their hybrid zone reveals neither intermediate nor transgressive morphotypes

Under incomplete reproductive isolation, secondary contact of diverged allopatric lineages may lead to the formation of hybrid zones that allow to study recombinants over several generations as excellent systems of genomic interactions resulting from the evolutionary forces acting on certain genes and phenotypes. Hybrid phenotypes may either exhibit intermediacy or, alternatively, transgressive traits, which exceed the extremes of their parents due to epistasis and segregation of complementary alleles. While transgressive morphotypes have been examined in fish, reptiles, birds, and mammals, studies in amphibians are rare. Here, we associate microsatellite-based genotypes with morphometrics-based morphotypes of two tree frog species of the Hyla arborea group, sampled across a hybrid zone in Poland, to understand whether the genetically differentiated parental species also differ in morphology between each other and their hybrids and whether secondary contact leads to the evolution of intermediate or transgressive morphotypes. Using univariate approaches, explorative multivariate methods (principal component analyses) as well as techniques with prior grouping (discriminant function analyses), we find that morphotypes of both parental species and hybrids differ from each other. Importantly, hybrid morphotypes are neither intermediate nor transgressive but found to be more similar to H. orientalis than to H. arborea.

Gene flow in phylogenomics: Sequence capture resolves species limits and biogeography of Afromontane forest endemic frogs from the Cameroon Highlands

Puddle frogs of the Phrynobatrachus steindachneri species complex are a useful group for investigating speciation and phylogeography in Afromontane forests of the Cameroon Volcanic Line, western Central Africa. The species complex is represented by six morphologically relatively cryptic mitochondrial DNA lineages, only two of which are distinguished at the species level - southern P. jimzimkusi and Lake Oku endemic P. njiomock, leaving the remaining four lineages identified as 'P. steindachneri'. In this study, the six mtDNA lineages are subjected to genomic sequence capture analyses and morphological examination to delimit species and to study biogeography. The nuclear DNA data (387 loci; 571,936 aligned base pairs) distinguished all six mtDNA lineages, but the topological pattern and divergence depths supported only four main clades: P. jimzimkusi, P. njiomock, and only two divergent evolutionary lineages within the four 'P. steindachneri' mtDNA lineages. One of the two lineages is herein described as a new species, P. amieti sp. nov. Reticulate evolution (hybridization) was detected within the species complex with morphologically intermediate hybrid individuals placed between the parental species in phylogenomic analyses, forming a ladder-like phylogenetic pattern. The presence of hybrids is undesirable in standard phylogenetic analyses but is essential and beneficial in the network multispecies coalescent. This latter approach provided insight into the reticulate evolutionary history of these endemic frogs. Introgressions likely occurred during the Middle and Late Pleistocene climatic oscillations, due to the cyclic connections (likely dominating during cold glacials) and separations (during warm interglacials) of montane forests. The genomic phylogeographic pattern supports the separation of the southern (Mt. Manengouba to Mt. Oku) and northern mountains at the onset of the Pleistocene. Further subdivisions occurred in the Early Pleistocene, separating populations from the northernmost (Tchabal Mbabo, Gotel Mts.) and middle mountains (Mt. Mbam, Mt. Oku, Mambilla Plateau), as well as the microendemic lineage restricted to Lake Oku (Mt. Oku). This unique model system is highly threatened as all the species within the complex have exhibited severe population declines in the past decade, placing them on the brink of extinction. In addition, Mount Oku is identified to be of particular conservation importance because it harbors three species of this complex. We, therefore, urge for conservation actions in the Cameroon Highlands to preserve their diversity before it is too late.

The Legacy of Sexual Ancestors in Phenotypic Variability, Gene Expression, and Homoeolog Regulation of Asexual Hybrids and Polyploids

Hybridization and polyploidization are important evolutionary processes whose impacts range from the alteration of gene expression and phenotypic variation to the triggering of asexual reproduction. We investigated fishes of the Cobitis taenia-elongatoides hybrid complex, which allowed us to disentangle the direct effects of both processes, due to the co-occurrence of parental species with their diploid and triploid hybrids. Employing morphological, ecological, and RNAseq approaches, we investigated the molecular determinants of hybrid and polyploid forms. In contrast with other studies, hybridization and polyploidy induced relatively very little transgressivity. Instead, Cobitis hybrids appeared intermediate with a clear effect of genomic dosing when triploids expressed higher similarity to the parent contributing two genome sets. This dosage effect was symmetric in the germline (oocyte gene expression), interestingly though, we observed an overall bias toward C. taenia in somatic tissues and traits. At the level of individual genes, expression-level dominance vastly prevailed over additivity or transgressivity. Also, trans-regulation of gene expression was less efficient in diploid hybrids than in triploids, where the expression modulation of homoeologs derived from the "haploid" parent was stronger than those derived from the "diploid" parent. Our findings suggest that the apparent intermediacy of hybrid phenotypes results from the combination of individual genes with dominant expression rather than from simple additivity. The efficiency of cross-talk between trans-regulatory elements further appears dosage dependent. Important effects of polyploidization may thus stem from changes in relative concentrations of trans-regulatory elements and their binding sites between hybridizing genomes. Links between gene regulation and asexuality are discussed.

Variation in hybridogenetic hybrid emergence between populations of water frogs from the Pelophylax esculentus complex

Many closely related species are capable of mating to produce hybrid offspring, which are usually sterile. Nevertheless, altering the gametogenesis of hybrid offspring can rescue hybrids from sterility by enabling asexual reproduction. Hybridogenesis is one of the most complicated asexual reproductive modes, and it includes drastic genome reorganization only in the germline; this is achieved through elimination of one parental genome and duplication of the remaining one to restore diploid chromosomal set and overcome blocks in meiotic progression. We investigated a model of hybridogenesis, namely, water frogs from the Pelophylax esculentus complex, for the emergence of asexual reproduction. Further, we assessed the impact of its asexual reproduction on the maintenance of interspecies hybrids from two populations on the western edge of the P. esculentus range, in which hybrids coexist with either both parental species or with only one parental species. After analysing tadpole karyotypes, we conclude that in both studied populations, the majority of diploid hybrid males produced haploid gametes with the P. ridibundus genome after elimination of the P. lessonae genome. Hybrid females exhibited problems with genome elimination and duplication; they usually produced oocytes with univalents, but there were observations of individual oocytes with 13 bivalents and even 26 bivalents. In some hybrid tadpoles, especially F1 crosses, we observed failed germ cell development, while in tadpoles from backcrosses, germ cells were normally distributed and contained micronuclei. By identifying chromosomes present in micronuclei, we estimated that the majority of tadpoles from all crosses were able to selectively eliminate the P. lessonae chromosomes. According to our results, hybridogenesis in hybrids can appear both from crosses of parental species and crosses between sexual species with hybrid individuals. The ability to eliminate a genome and perform endoreplication to ensure gamete formation differed between male and female hybrids from the studied populations. Some diploid hybrid females can rarely produce not only haploid gametes but also diploid gametes, which is a crucial step in the formation of triploid hybrids.

Cerebrovascular β-dystroglycan immunoreactivity in vertebrates: not detected in anurans and in the teleosts Ostariophysi and Euteleostei

The aim of the present paper was to check for the presence of cerebrovascular dystroglycan in vertebrates, because dystroglycan, which is localized in the vascular astroglial end-feet, has a pivotal function in glio-vascular connections. In mammalian brains, the immunoreactivity of β-dystroglycan subunit delineates the vessels. The results of the present study demonstrate similar patterns in other vertebrates, except for anurans and the teleost groups Ostariophysi and Euteleostei. In this study, we investigated 1 or 2 representative species of the main groups of Chondrichthyes, teleost and non-teleost ray-finned fishes, urodeles, anurans, and reptiles. We also investigated 5 mammalian and 3 bird species. Animals were obtained from breeders or fishermen. The presence of β-dystroglycan was investigated immunohistochemically in free-floating sections. Pre-embedding electron microscopical immunohistochemistry on Heterodontus japonicus shark brains demonstrated that in Elasmobranchii, β-dystroglycan is also localized in the perivascular glial end-feet despite the different construction of their blood-brain barrier. The results indicated that the cerebrovascular β-dystroglycan immunoreactivity disappeared separately in anurans, and in teleosts, in the latter group before its division to Ostariophysi and Euteleostei. Immunohistochemistry in muscles and western blots from brain homogenates, however, detected the presence of β-dystroglycan, even in anurans and all teleosts. A possible explanation is that in the glial end-feet, β-dystroglycan is masked in these animals, or disappeared during adaptation to the freshwater habitat.

The programmed DNA elimination and formation of micronuclei in germ line cells of the natural hybridogenetic water frog Pelophylax esculentus

DNA elimination is a radical form of gene silencing and occurs both in somatic and germ cells. The programmed DNA elimination occurs during gametogenesis in interspecies hybrids that reproduce by hybridogenesis (stick insects, fishes, and amphibians) and concerns removal of whole genomes of one of the parental species and production of clonal gametes propagating the genome of the other species. The cellular mechanisms differ considerably in hybridogenetic insects and fishes but remains unknown in edible frogs Pelophylax esculentus, natural hybrids between Pelophylax lessonae and Pelophylax ridibundus. Here we report DNA elimination mechanism in early developing gonads of diploid and triploid hybrid frogs, studied by TEM, immunofluorescence, and cytochemistry. In gonocytes of both sexes (primary oogonia and prespermatogonia), micronuclei emerge as detached nuclear buds formed during interphase. We found depletion of nuclear pore complexes in micronuclear membrane and chromatin inactivation via heterochromatinization followed by degradation of micronuclei by autophagy. Micronuclei formation does not lead to apoptotic cell death showing that genome elimination is a physiological process. Chromatin elimination via micronuclei in P. esculentus is unique among hybridogenetic animals and contributes to broadening the knowledge about reproductive modes in animals.

hybriddetective: A workflow and package to facilitate the detection of hybridization using genomic data in r

The ability to detect and characterize hybridization in nature has long been of interest to many fields of biology and often has direct implications for wildlife management and conservation. The capacity to identify the presence of hybridization, and quantify the numbers of individuals belonging to different hybrid classes, permits inference on the magnitude of, and timescale over which, hybridization has been or is occurring. Here, we present an r package and associated workflow developed for the detection, with estimates of efficiency and accuracy, of multigenerational hybrid individuals using genetic or genomic data in conjunction with the program newhybrids. This package includes functions for the identification and testing of diagnostic panels of markers, the simulation of multigenerational hybrids, and the quantification and visualization of the efficiency and accuracy with which hybrids can be detected. Overall, this package delivers a streamlined hybrid analysis platform, providing improvements in speed, ease of use and repeatability over current ad hoc approaches. The latest version of the package and associated documentation are available on GitHub (https://github.com/bwringe/hybriddetective).

Variability of NOR patterns in European water frogs of different genome composition and ploidy level

We studied water frogs from a complex composed of two species: Pelophylax lessonae (Camerano, 1882) (genome LL, 2n = 26) and P. ridibundus (Pallas, 1771) (RR, 2 = 26), and their natural hybrid P. esculentus (Fitzinger, 1843) of various ploidy and genome composition (RL, 2n = 26, and RRL or RLL, 3n = 39). Tetraploids RRLL were found (4n = 52) in juveniles. We applied cytogenetic techniques: AgNO3, chromomycin A3, PI and fluorescent in situ hybridization with a 28S rDNA probe. Results obtained by silver staining corresponded well with those stained with CMA3, PI and FISH. As a rule, NORs are situated on chromosomes 10. The number of Ag-NORs visible on metaphase plates was the same as the number of Ag-nucleoli present in interphase nuclei of the same individual. In all analyzed metaphases, NORs exhibited variations in size after AgNO3 and CMA3 stainings. Sixty-six individuals (out of 407 analyzed) were polymorphic for the localization and number of NORs. Fifty-one diploids had NORs only on one chromosome of pair 10. Three triploids (LLR and RRL) displayed two NORs, and two other triploid RRL individuals displayed one, instead of expected three NORs. In ten individuals extra NORs were detected on chromosomes other than 10 (chromosomes 2 and 9).

A general model of distant hybridization reveals the conditions for extinction in Atlantic salmon and brown trout

Interspecific hybridization is common in nature but can be increased in frequency or even originated by human actions, such as species introduction or habitat modification, which may threaten species persistence. When hybridization occurs between distantly related species, referred to as "distant hybridization," the resulting hybrids are generally infertile or fertile but do not undergo chromosomal recombination during gametogenesis. Here, we present a model describing this frequent but poorly studied interspecific hybridization to assess its consequences on parental species and to anticipate the conditions under which they can reach extinction. Our general model fully incorporates three important processes: density-dependent competition, dominance/recessivity inheritance of traits and assortative mating. We demonstrate its use and flexibility by assessing population extinction risk between Atlantic salmon and brown trout in Norway, whose interbreeding has recently increased due to farmed fish releases into the wild. We identified the set of conditions under which hybridization may threaten salmonid species. Thanks to the flexibility of our model, we evaluated the effect of an additional risk factor, a parasitic disease, and showed that the cumulative effects dramatically increase the extinction risk. The consequences of distant hybridization are not genetically, but demographically mediated. Our general model is useful to better comprehend the evolution of such hybrid systems and we demonstrated its importance in the field of conservation biology to set up management recommendations when this increasingly frequent type of hybridization is in action.

Genomic effects on advertisement call structure in diploid and triploid hybrid waterfrogs (Anura, Pelophylax esculentus)

BACKGROUND

In anurans, differences in male mating calls have intensively been studied with respect to taxonomic classification, phylogeographic comparisons among different populations and sexual selection. Although overall successful, there is often much unexplained variation in these studies. Potential causes for such variation include differences among genotypes and breeding systems, as well as differences between populations. We investigated how these three factors affect call properties in male water frogs of Pelophylax lessonae (genotype LL), P. ridibundus (RR) and their interspecific hybrid P. esculentus which comes in diploid (LR) and triploid types (LLR, LRR).

RESULTS

We investigated five call parameters that all showed a genomic dosage effect, i.e. they either decreased or increased with the L/R ratio in the order LL-LLR-LR-LRR-RR. Not all parameters differentiated equally well between the five genotypes, but combined they provided a good separation. Two of the five call parameters were also affected by the breeding system. Calls of diploid LR males varied, depending on whether these males mated with one or both of the parental species (diploid systems) or triploid hybrids (mixed ploidy systems). With the exception of the northernmost mixed-ploidy population, call differences were not related to the geographic location of the population and they were not correlated with genetic distances in the R and L genomes.

CONCLUSIONS

We found an influence of all three tested factors on call parameters, with the effect size decreasing from genotype through breeding system to geographic location of the population. Overall, results were in line with predictions from a dosage effect in L/R ratios, but in three call parameters all three hybrid types were more similar to one or the other parental species. Also calls of diploid hybrids varied between breeding systems in agreement with the sexual host required for successful reproduction. The lack of hybrid call differences in a mixed-ploidy population at the northern edge of the water frog distribution is likely to be associated with genetic particularities, including a) low genetic variability and/or b) a local loss of genes coding for genotype-dependent call differentiation under conditions where female discrimination between diploid and triploid males is not beneficial.

A simplified molecular method for distinguishing among species and ploidy levels in European water frogs (Pelophylax)

Western Palearctic water frogs in the genus Pelophylax are a set of morphologically similar anuran species that form hybridogenetic complexes. Fully reliable identification of species and especially of hybrid ploidy depends on karyological and molecular methods. In central Europe, native water frog populations consist of the Pelophylax esculentus complex, that is, P. lessonae (LL), P. ridibundus (RR) and the hybrid form P. esculentus that can have different karyotypes (RL, LLR and RRL). We developed existing molecular methods further and propose a simple PCR method based on size-differences in the length of the serum albumin intron-1 and the RanaCR1, a non-LTR retrotransposon of the chicken repeat (CR) family. This PCR yields taxon-specific banding patterns that can easily be screened by standard agarose gel electrophoresis and correctly identify species in all of the 160 samples that had been identified to karyotype with other methods. To distinguish ploidy levels in LR, LLR and RRL specimens, we used the ratio of the peak heights of the larger (ridibundus specific) to the smaller (lessonae specific) bands of fluorescently labelled PCR products resolved on a capillary DNA sequencer and obtained a correct assignment of the karyotype in 93% of cases. Our new method will cut down time and expenses drastically for a reliable identification of water frogs of the P. esculentus complex and potentially for identification of other hybridogenetic complexes and/or taxa, and it even serves as a good indicator of the ploidy status of hybrid individuals.