AFLP reveals low genetic diversity of the bryozoan Pectinatella magnifica (Leidy, 1851) in the Czech Republic
© The Author(s) 2017
Received: 3 March 2017
Accepted: 14 November 2017
Published: 25 November 2017
Non-native species have aroused scientific interest because of their ability to successfully colonise areas to which they have been introduced, despite their sometimes limited genetic variation compared to their native range. These species establish themselves with the aid of some pre-existing features favouring them in the new environment. Pectinatella magnifica (Leidy, 1851), the freshwater magnificent bryozoan, is non-native in Europe and Asia. This study was designed to determine the genetic diversity and population structure of P. magnifica colonies collected from the Protected Landscape Area (PLA) and UNESCO Biosphere Reserve Třeboňsko (the Czech Republic) in the 2009 and 2011–2014 periods using Amplified Fragment Length Polymorphism (AFLP).
The vast majority of the examined non-native colonies, except three colonies sampled in 2012, expressed very low levels of genetic variation, not differentiating from the USA native colony. The Bayesian clustering approach grouped the 28 accessions into two genetically different populations.
The data suggest relatively low gene diversity within all colonies, which might reflect the recent expansion of P. magnifica in the Czech Republic.
Pectinatella magnifica (Leidy, 1851), a freshwater bryozoan species that naturally inhabits lakes and rivers in eastern North America [1, 2], has been recently spreading in Europe and western Asia [3–5]. Its colonies begin life from a germinated statoblast (a product of asexual reproduction) or settled larva (the result of sexual reproduction). Once formed, colonies exude a protective gelatin-like matrix of exopolymers, enabling attachment to different submerged natural or artificial substrates  as flat sheets (in the form of thin films), as sessile globular colonies, or as gelatinous balls floating in water. The colonies are covered on their surface with rosettes, each containing 12–18 genetically identical zooids and can lead to considerable dimension and weight, as in the case of nearly 1 m diameter and 70 kg colony described by Balounová et al. . Recent studies have focused on the determination of evolutionary relationships within P. magnifica, incorporating morphological traits (such as the morphology of the statoblasts, white spots on lophophora, pigmented mouth), molecular approaches, or a combination of both morphological traits and genetic data . Most systematic (phylogenetic) studies including P. magnifica have used nuclear ribosomal (ssrDNA, lsrDNA) and mitochondrial (rrnL, rrnS, cox1, cox3, cytb) genes [9–13]. However, to determine intraspecific relatedness, ribosomal and mitochondrial genes are not sufficient, as demonstrated for instance in the case of using mitochondrial cytochrome c oxidase subunit I gene (cox1) for intraspecific analysis of geographically distinct colonies of P. magnifica in Korea . AFLP markers are considered as one of the most popular tools for genetic analysis in the fields of evolutionary genetics, ecology and conservation of genetic resources . The technique combines a high-information content and fidelity with the possibility of carrying out genome wide scans. However, a potential problem with this technique is the lack of homology of bands with the same electrophoretic mobility, what is known as fragment-size homoplasy . This study is the first of its kind, where the random whole genome AFLP analysis was applied to determine relatedness of geographically and time-distant colonies of P. magnifica. The results provide the first comprehensive summary of the methodology used.
Total genomic DNA was extracted from statoblasts by using DNeasy® Blood & Tissue kit (Qiagen, Valencia, CA, USA) according to the manufacturer’s protocol. The quality and quantity of obtained DNA was verified by using a spectrophotometer BioSpec Nano (Shimadzu, Tokyo, Japan). The DNA was tested by using specific primers (Pm28SF: 5′-CTTTTCAGCCAAGCACATGA-3′ and Pm28SR: 5′-AAACCTCGCGTAGAGCAAAA-3′), designed for this study with Primer3 software , amplifying a fragment of 398 bp. Primers were designed based on the unique P. magnifica 28S rRNA gene sequence available in GenBank (FJ409576.1), and were chosen on the basis of their performance and specificity against alternative primers. The specificity was confirmed by sequencing. PCR was performed with 3 ng of DNA template, mixed with 1 µl (10 µM) of each primer, 12.5 µl PPP MM (TopBio, Prague, Czech Republic) and distilled water to make a final volume of 20 µl. PCR was performed using the following programme: pre-denaturation for 5 min at 95 °C, 35 cycles of 30 s at 95 °C, 1 min 60 °C and 40 s at 72 °C, finally, 10 min at 72 °C. PCR products were subjected to electrophoresis on 1.5% agarose gel in 1 × TBE buffer, and detected through ethidium bromide staining.
AFLP analysis was carried out as described by Vos et al.  with four primer combinations: fluorescent dye-labelled primer EcoRI + ACG/MseI + AGT, EcoRI + ACG/MseI + ACC, EcoRI + ACG/MseI + AAC and EcoRI + ACG/MseI + ATT.
One microliter of selective PCR amplicons, 0.4 μl of the GeneScan™ 500 LIZ® (Life Technologies, USA) and 11 μl of formamide, were denatured at 95 °C for 5 min, cooled on ice immediately and directly forwarded to an capillary electrophoresis on an ABI 3500 genetic analyzer (Life Technologies, USA) and data were analysed by GeneMapper 5.0 software (Life Technologies, USA). Signals obtained from fragment analyzer were transformed into a binary character matrix with 1 for the presence or 0 for the absence of a band at a particular position.
Cluster and principal component analysis (PCA) were performed. These analyses were calculated by using R software (cran.r-project.org, version 3.2.3, accessed in 2015). Analysis of the genetic structure of the population was carried out assigning samples to K populations according to allele frequencies of each locus by using Bayesian model-based clustering method with Structure 2.3.4 . The program was run 10 times at each of 10 different K values (1–10) with a period of 5000 burn-in steps followed by 50,000 iterations. The admixture model with independent allele frequencies was used in all analyses. The number of populations was inferred by plotting the log probability of the data [Ln P(D)] for each K value. Also, in order to determine the most appropriate number of K, ΔK was calculated as described by Evanno et al. .
Restriction enzyme primer combinations from Amplified Fragment Polymorphism
AFLP selective primer extensions
No. polymorphic bands
The STRUCTURE results divided samples into two groups (K = 2) with no consistent trend differentiating individual localities or years of collection by using these groups (Additional file 2: Table S2).
Recent studies show that the level of expansion of the non-native Pectinatella magnifica is increasing not only in Europe and Asia, but also within the country of its native origin [2, 4, 5, 8, 20]. The AFLP technique is very popular in population and evolutionary genetic studies of organisms with no prior sequence knowledge, aiming to assess the degree of variability and genetic structure among tested individuals. The advantage of this technique is genome wide screening, which provides a high information content. However, every methodology can have challenges and in the case of AFLP, the evaluation of fragments can be misleading due to size migration, which may lead to assigning fragment lengths to loci of different origin in the genome .
In this study, the AFLP technique has been used to characterize the population structure of P. magnifica, with the aim to better understand the invasion pattern of this species. The results showed a low level of genetic diversity within all colonies sampled, with low differentiation between sampling locations. Also, STRUCTURE results support the hypothesis that a single metapopulation exists without any significant genetic differentiation. In the study of Kang and An , the comparison of 20 sequences of the cox1 gene from P. magnifica colonies collected from three different watersheds in Korea did not show any variation between the colonies. Low genetic diversity was also observed in other plant and animal invaders, including species which utilize both, sexual and asexual reproduction [21–24]. The main reason for low genetic diversity within populations of non-native organisms in secondary areas of occurrence may be attributed to their recent establishment. Most likely, only a small fraction of the source germplasm may have been the initial contaminant, generating an extreme diversity bottleneck in the newly established population with weak or low population structure. The results indicate a recent invasion, where colonies likely originate from single or closely related sources. This fact clearly demonstrates that a high amount of genetic diversity is not necessarily a condition for global spread over contemporary timescales. Given the fact that P. magnifica was first reported in Europe in 1883 , a low diversity across Europe is expected, although it could also be possible that P. magnifica was previously undetected. The results of genetic analyses performed in this study on P. magnifica colonies sampled within 5 years across different localities correspond with the so called genetic paradox in invasive species, where invasive populations are able to overcome low genetic diversity and inbreeding to thrive in the invaded environment. During optimal conditions, colonies increase in size through rapid repeated asexual budding of new zooids, but sexual events also occur regularly and play an important role in creating genetic variation in progeny. However, based on the assumption of a low initial degree of total allelic variation in all colonizing individuals, rapid increase of genetic diversity cannot be expected through recombination events and it is not clear how important the amounts of genetic diversity are for persistence and spread in the future.
The genetic diversity between different transects and over the years was not very different, but showed a clustering pattern in term of years, where samples taken in the same year grouped more closely than the samples from the same locality collected within years. The USA sample representing the native genetic resource did not differ from other samples representing invasive genomes. This phenomenon could be explained as being a result of low general diversity between these genotypes but also as erroneous estimation of co-migrating bands in AFLP analysis.
In conclusion, our results showed the presence of a single metapopulation in all recently known localities across the Czech Republic with low or slight differentiation of P. magnifica colonies from the different locations or over the years respectively, and relatively low gene diversity. We suggest that the genetic pattern observed in this study might reflect the recent expansion of P. magnifica in PLA (during the last decades), and the colonization of new water streams, dams and ponds. In spite of its geographic distance, the sample from the native origin has shown no particular AFLP band pattern.
JR, VC, LH conceived and planned the project. LH, VM, ZB, carried out the experiments. LH, JM and VM performed data analysis. LH wrote the manuscript and all authors reviewed it. All authors read and approved the final manuscript.
The authors thank Dr. Roxana I. Teodor and Philip Hughes for critically editing the manuscript and Dr. Barbora Kubátová for computational assistance during experiment evaluation.
The authors declare that they have no competing interests.
Availability of data and materiáls
All supporting data are available as Additional data.
Consent for publication
All the authors have read the paper and have agreed to having their names listed as authors. All authors authorize publication in the Journal of Biological Research-Thessaloniki.
Ethics approval and consent to participate
This study was supported by the Czech Science Foundation (Grant No. P503/12/0337) and by Grant Agency of the University of South Bohemia (Grant No. GA JU 120/2016/Z and GA JU 081/2016/Z).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Michelena TM, Ostman C, Boylen CW, Nierzwicki-Bauer SA. Records of bryozoans in the freshwater reach of the Hudson River Estuary. Northeast Nat. 2014;21:369–79.View ArticleGoogle Scholar
- Balounová Z, Pechoušková E, Rajchard J, Joza V, Šinko J. World-wide distribution of the bryozoan Pectinatella magnifica (Leidy 1851). Eur J Environ Sci. 2013;3:96–100.Google Scholar
- Jo H, Joo G-J, Byeon M, Hong D-G, Gim J-S, Kim J-Y, Choi J-Y. Distribution pattern of Pectinatella magnifica (Leidy, 1851), an invasive species, in the Geum river and the Nakdong River, South Korea. J Ecol Environ. 2014;37:217–23.View ArticleGoogle Scholar
- Aleksandrov B, Voloshkevich O, Kurakin A, Rybalko A, Gontar V. The first finding of bryozoan Pectinatella magnifica (Lophopodidae) in lower Danube. Vestnik Zool. 2014;48:307–12.Google Scholar
- Zorić K, Szekeres J, Csányi B, Kolarevic S, Marković V, Paunović M. Distribution of the non-native bryozoan Pectinatella magnifica (Leidy, 1851) in the Danube River. Acta Zool Bulg. 2015;67:241–7.Google Scholar
- Choi J-Y, Joo G-J, Kim S-K, Hong D-G, Jo H. Importance of substrate material for sustaining the bryozoan Pectinatella magnifica following summer rainfall in lotic freshwater ecosystem, South Korea. J Ecol Environ. 2015;38:375–81.View ArticleGoogle Scholar
- Balounová Z, Brezina V, Susterova K, Rajchard J. Primary screening of potentially bio-active substances in the lyophilisate of Pectinatella magnifica biomass. Vet Med-Czech. 2015;60:141–6.View ArticleGoogle Scholar
- Kang N, An K-G. Statoblast ultrastructure and genetic identity of Pectinatella magnifica population based on COI gene, from three different watersheds in Korea. Anim Cells Syst. 2015;19:78–84.View ArticleGoogle Scholar
- Wood TS, Lore M. The higher phylogeny of Phylactolaemate bryozoans inferred from 18S ribosomal DNA sequences. In: Jackson PNW, Cancino JM, Moyano HI, editors. Bryozoan Studies. London: Taylor and Francis Group; 2004. p. 361–8.Google Scholar
- Okuyama M, Wada H, Ishii T. Phylogenetic relationships of freshwater bryozoans (Ectoprocta, Phylactolaemata) inferred from mitochondrial ribosomal DNA sequences. Zool Scr. 2006;35:243–9.View ArticleGoogle Scholar
- Fuchs J, Obst M, Sundberg P. The first comprehensive molecular phylogeny of Bryozoa (Ectoprocta) based on combined analyses of nuclear and mitochondrial genes. Mol Phylogenet Evol. 2009;52:225–33.View ArticlePubMedGoogle Scholar
- Waesenbach A, Cox CJ, Littlewood DTJ, Porter JS, Taylor PD. First molecular estimate of cyclostome bryozoan phylogeny confirms extensive homoplasy among skeletal characters used in traditional taxonomy. Mol Phylogenet Evol. 2009;52:241–51.View ArticleGoogle Scholar
- Waesenbach A, Taylor PD, Littlewood DTJ. A molecular phylogeny of bryozoans. Mol Phylogenet Evol. 2012;62:718–35.View ArticleGoogle Scholar
- Caballero A, Quesada H, Rolán-Alvarez E. Impact of amplified fragment length polymorphism size homoplasy on the estimation of population genetic diversity and the detection of selective loci. Genetics. 2008;179:539–54.View ArticlePubMedPubMed CentralGoogle Scholar
- Caballero A, Quesada H. Homoplasy and distribution of AFLP fragments: an analysis in silico of the genome of different species. Mol Biol Evol. 2010;27:1139–51.View ArticlePubMedGoogle Scholar
- Untergasser A, Cutcutache I, Koressaar T, Ye J, Faircloth BC, Remm M, et al. Primer3—new capabilities and interfaces. Nucleic Acids Res. 2012;40:e115.View ArticlePubMedPubMed CentralGoogle Scholar
- Vos P, Hogers R, Bleeker M, Reijans M, van de Lee T, Hornes M, et al. AFLP: a new technique for DNA fingerprinting. Nucleic Acids Res. 1995;23:4407–14.View ArticlePubMedPubMed CentralGoogle Scholar
- Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000;155:945–59.PubMedPubMed CentralGoogle Scholar
- Evanno G, Regnaut S, Goudet J. Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol. 2005;14:2611–20.View ArticlePubMedGoogle Scholar
- Taticchi MI, Elia AC, Battoe L, Havens KE. First report about freshwater Bryozoa in Florida (Lake Apopka). Ital J Zool. 2009;76:194–200.View ArticleGoogle Scholar
- Maron JL, Vilá M, Bommarco R, Elmendorf S, Beardsley P. Rapid evolution of an invasive plant. Ecol Monographs. 2004;74:261–80.View ArticleGoogle Scholar
- Dlugosch KM, Parker M. Founding events in species invasions: genetic variation, adaptive evolution, and the role of multiple introductions. Mol Ecol. 2008;17:431–49.View ArticlePubMedGoogle Scholar
- Loomis ES, Fishman L. A continent-wide clone: population genetic variation of the invasive plant Hieracium aurantiacum (Orange hawkweed; Asteraceae) in North America. Int J Plant Sci. 2009;170:759–65.View ArticleGoogle Scholar
- Lebret K, Kritzberg ES, Rengefors K. Population genetic structure of a microalgal species under expansion. PLoS ONE. 2013;8:e82510.View ArticlePubMedPubMed CentralGoogle Scholar