Introduction
The red seaweed genus Porphyra, sister to the genus Pyropia (Sutherland et al., Reference Sutherland, Lindstrom, Nelson, Brodie, Lynch, Hwang, Choi, Miyata, Kikuchi, Oliveira, Farr, Neefus, Mols-Mortensen, Milstein and Müller2011), forms part of very valuable (Blouin et al., Reference Blouin, Brodie, Grossman, Xu and Brawley2011) marine food products commonly known as ‘Nori’ (Sahoo et al., Reference Sahoo, Tang and Yarish2002). Porphyra umbilicalis Kützing is common in Europe and it is traditionally harvested in Great Britain (‘purple laver’) for human consumption (Blouin et al., Reference Blouin, Brodie, Grossman, Xu and Brawley2011). Farming of these valuable crops was revolutionized by the discovery of their life cycle by Drew (Reference Drew1949). Porphyra species typically have a heteromorphic life cycle with a gametophytic blade phase (haplophasic phase) and a filamentous sporophyte called the conchocelis (diplophasic phase). The sporophytic conchocelis phase is difficult to find in nature with only a few in situ records (e.g. Miura, Reference Miura1961; Martinez, Reference Martinez1990). To get further insight into the biology of P. umbilicalis, seven novel polymorphic microsatellite markers, isolated from enriched DNA libraries are reported in this paper.
Experimental
For library preparation, ten gametophytic fronds of P. umbilicalis were collected from Punta Herminia, Coruña, Spain, in April 2010 (43°23′23″N, 8°24′01″W). Specimens from this population were placed at the herbarium of the University of Santiago de Compostela, Spain (SANT-Algae No: 24072). Genomic DNA was isolated using a nuclei isolation protocol (Varela-Álvarez et al., Reference Varela-Álvarez, Andreakis, Lago-Lestón, Pearson, Serrão, Procaccini, Duarte and Marbá2006) followed by a LiCl extraction protocol described by Hong et al. (Reference Hong, Coury, Polne-Fuller and Gibor1992) as modified by van Oppen et al. (Reference van Oppen, Olsen and Stam1995). Four microsatellite enriched libraries were constructed following Varela-Álvarez et al. (Reference Varela-Álvarez, Paulino and Serrão2017), based on Glenn and Schable (Reference Glenn, Schable, Zimmer and Roalson2005) and Billote et al. (Reference Billote, Lagoda, Risterucci and Baurens1999). We also used nucleotide repeats that were complementary to those in other Porphyra spp. microsatellites: P. haitanensis (Zuo et al., Reference Zuo, Wang, Cao, Su, Liao and Chen2007) and P. yezoensis (Kong et al., Reference Kong, Mao, Yang, Qu, Yan and Wang2009) (recently transferred to the genus Pyropia (Sutherland et al., Reference Sutherland, Lindstrom, Nelson, Brodie, Lynch, Hwang, Choi, Miyata, Kikuchi, Oliveira, Farr, Neefus, Mols-Mortensen, Milstein and Müller2011)). From the libraries (development details in supplementary materials), 146 sequences were recovered, of which, 97 contained a total number of 241 repeats identified using MSATCOMMANDER (Faircloth, Reference Faircloth2008). Using these, 68 primer pairs were designed with PRIMER3 (Rozen and Skaletsky, Reference Rozen, Skaletsky, Krawetz and Misener2000). Primers were tested under two different strategies. A subset of 38 primers pairs were tested with a labelled tail, where one primer in each pair was modified on the 5′ end with a CAG sequence tail (5′-CAGTCGGGCGTCATCA-3′) to allow the use of a fluorescently labelled third primer in the PCR. The CAG tail was removed after primer testing and optimization. The other 30 primers were designed without any tail. For primer optimization we used first six individuals from the same population used for the library (Punta Herminia, Coruña, Spain) and two positive clones (the plasmid isolated from the library and the DNA template used to build the library). We selected seven primers to test for polymorphism in a panel of 44 individuals from four populations (11 individuals per population). These populations covered were: St Ives and Polzeath, UK (50°12′56″N, 5°28′56″W; 50°34′34″N, 4°55′12″W), Pta. Herminia, Coruña (43°23′23″N, 8°24′01″W) and Sanxenxo, Pontevedra, Spain (42°23′58″N, 8°48′54″W). Cross amplification on other species was tested, using on five individuals of Porphyra dioica Brodie and L.M. Irvine from Oeiras, Portugal (38°41′08″N, 9°18′33″W) and four individuals of Porphyra linearis Greville from São Bartolomeu do Mar, Portugal (41°34′27″N, 8°48′00″W).
Polymerase chain reactions (PCR) for all markers were performed separately in a 20 µl reaction volume containing 5–50 ng genomic DNA. For PCR using non-tailed primers, PCR reactions contained 0.25 mM of each primer, 0.2 mM of dNTPs, 2–3 mM of MgCl2 (see Table 1 for loci optimizations), 1× GoTaq Flexi Buffer and 1 U GoTaq DNA polymerase (Promega). Forward primers were labelled either with FAM or HEX dye. For PCR including a third labelled primer, we used 0.25 µM unlabelled primer, 0.025 µM CAG tail-primer, and 0.225 µM universal dye-labelled primer. Amplifications using an Applied Biosystems thermal cycler (GeneAmp 2720) were conducted with the following three programs: (a) Program 1: 95°C for 5 min followed by 35 cycles at 95°C for 20 s, annealing at a specific temperature for each locus (see Table 1 for details) for 20 s (except for Por_021 and Por_042, only 15 s annealing), 72°C for 30 s, and a final extension step at 72°C for 10 min (b) Program 2: Program 1 but with 30 cycles (c) Program 3: Touchdown thermal cycling program consisting of an initial denaturing step of 3 min at 95°C, followed by 14 cycles of ‘touchdown’ PCR consisting of 15 s at 95°C, 15 s at 62°C (reduced by 0.5°C each subsequent cycle), and 15 s at 72°C, followed by 24 cycles consisting of 15 s at 95°C, 15 s at 55°C and 15 s at 72°C; with a final elongation step at 72°C for 10 min. PCR products were run on an ABI PRISM 3130xl sequencer (Applied Biosystems) and sized with GeneScan-350ROX size standard. Results were analysed using Genemapper version 4.0 (Applied Biosystems) and Genodive (Meirmans and Van Tienderen, Reference Meirmans and Van Tienderen2004). Binning and allele rounding were also checked with TANDEM (Matschiner and Salzburger, Reference Matschiner and Salzburger2009).
Null alleles were detected in three loci (below in the table).
Discussion
Seven working primer pairs that showed polymorphism (amplification conditions are given in Table 1) with 2–12 alleles per locus. Null alleles were detected in three loci, and cross amplification with P. dioica and P. linearis was detected in five loci. This study found multiple alleles per locus (up to six alleles per individual), in agreement with Niwa and Sakamoto (Reference Niwa and Sakamoto2010) where polyploids (allopolyploids) with multiple alleles were found in some Japanese species of Porphyra sensu latto (now Pyropia spp. in Sutherland et al., Reference Sutherland, Lindstrom, Nelson, Brodie, Lynch, Hwang, Choi, Miyata, Kikuchi, Oliveira, Farr, Neefus, Mols-Mortensen, Milstein and Müller2011). Other possible causes for multiple alleles per locus could be: (a) blades can be genetic chimeras (Mitman and van der Meer, Reference Mitman and van der Meer1994), or (b) microsatellites amplify across as different members of gene families (Brawley et al., Reference Brawley, Blouin, Ficko-Blean, Wheeler, Lohr, Goodson, Jenkins, Blaby-Haas, Helliwell, Chan, Marriage, Bhattacharya, Klein, Badis, Brodie, Cao, Collén, Dittami, Gachon, Green, Karpowicz, Kim, Kudahl, Lin, Michel, Mittag, Olson, Pangilinan, Peng, Qiu, Shu, Singer, Smith, Sprecher, Wagner, Wang, Wang, Yan, Yarish, Zäuner-Riek, Zhuang, Zou, Lindquist, Grimwood, Barry, Rokhsar, Schmutz, Stiller, Grossman and Prochnik2017). These polymorphic microsatellite markers should be useful for investigating these different hypotheses and for inferring genetic differentiation between populations of P. umbilicalis in the North East Atlantic.
Supplementary material
The supplementary material for this article can be found at https://doi.org/10.1017/S147926211700034X.
Acknowledgements
Funded by FCT (Fundação para a Ciência e a Tecnologia, Portugal) via postdoctoral fellowship SFRH/BPD/109452/2015 to EVÁ, and the following FCT projects: (a) NORIGENOMICS – PTDC/MAR/099698/2008 to EVÁ, (b) GENEKELP-PTDC/MAR-EST/6053/2014 to EAS and (c) UID/Multi/04326/2013 to CCMAR.
We would like to thank Travis Glenn for the help/comments with the design of the probes for the enrichments, Mafalda Afonso and Luci Jesus for the lab work, and the people who helped the sample collection: Liam Cronin, Ignacio Bárbara and Tânia Pereira.