Erratum : A new Felimare ( Mollusca : Heterobranchia : Nudibranchia ) of the Atlantic blue chromodorid chromatic group from Cape Verde

1 Departamento de Biología, Facultad de Ciencias del Mar y Ambientales, Campus de Excelencia Internacional del Mar (CEI·MAR), Universidad de Cádiz, Av. República Saharahui, s/n, Apdo. 40, CP 11510 Puerto Real, Cádiz, Spain. (DO) (Corresponding author) E-mail: jazmin.ortigosa@gmail.com. ORCID iD: https://orcid.org/0000-0002-1857-4630 (JLC) E-mail: lucas.cervera@uca.es. ORCID iD: http://orcid.org/0000-0002-8337-2867 2 Present address: Calle 13#180-202 Ignacio Zaragoza; CP 15000 Ciudad de México, México. 3 Departamento de Biología, Edificio de Biología, Campus de Excelencia Internacional UAM+CSIC, Universidad Autónoma de Madrid, Madrid, Spain. (MP) E-mail: marta.pola@uam.es. ORCID iD: https://orcid.org/0000-0003-0518-346X 4 Instituto Universitario de Investigación Marina (INMAR), Campus de Excelencia Internacional del Mar (CEI·MAR), Universidad de Cádiz, Avenida República Saharaui, s/n, Ap. 40, 11510 Puerto Real, Cádiz, Spain.


INTRODUCTION
defined the Atlantic blue chromodorid chromatic group as species of blue colour (from pale blue to navy) with white, blue or orange spots, blotches, or lines, and with a distribution range in the Atlantic Ocean (Mediterranean and Caribbean included).This group comprised species of five formerly different genera: Hypselodoris Stimpson, 1855and Mexichromis Bertsch, 1977, nowadays transferred to Felimare Ev. Marcus and Er. Marcus, 1967;Glossodoris Ehrenberg, 1831 andChromodoris Alder andHancock, 1855, both as part of the provisional genus 'Felimida' Ev. Marcus, 1971;and Risbecia Odhner, 1934, synonymized by Johnson and Gosliner (2012) with Hypselodoris.Ortea's et al. (1996) review, however, did not include in this group other Atlantic species clearly belonging to the blue chromatic group, such as Felimare zebra (Heilprin, 1889) from Bermuda.
To date, in the Atlantic coast of Africa, there is evidence of nine species of Felimare that fit into the blue chromodorid chromatic group (Table 1).In this paper, we describe a new species of Felimare from specimens collected from Cape Verde using morphological characters as well as molecular analyses based on two mitochondrial genes, cytochrome c oxidase subunit I (COI) and 16S rRNA (16S), and one nuclear gene, histone-3 (H3).Thus, the number of the blue chromatic chromodorids from the western coast of Africa has risen to ten.

Phylogenetic analyses
Molecular analysis included a total of 21 specimens including eight species of Felimare, four specimens of three other genera of Chromodorididae Bergh, 1891 and one specimen of Prodoris clavigera (Thiele, 1912), originally ascribed to Bathydoris Bergh, 1884, as an outgroup (Table 2).
DNA sequences were assembled, edited and aligned using Geneious 8.1.2(Kearse et al. 2012).The alignments were checked by eye.All the sequences were checked for contamination with BLAST (Altschul et al. 1990) implemented in the GenBank database.Protein-coding sequences were translated into amino acids for confirmation of alignment.Pairwise uncorrected p-distance values between each taxon, uncorrected p-distances between all taxa, and level of saturation for first, second, and third codon positions (p-distances against transitions plus transversions) were calculated in MEGA 5.0 (Tamura et al. 2011) for the COI and H3

MB28-004390
Tarrafal, Cape Verde genes.The most variable regions from the 16S rRNA alignment were removed in the first analyses, using both the default settings and the standard options for stringent and less stringent selection in Gblocks (Talavera and Castresana 2007).When these regions were excluded from the analyses, the combined phylogenetic tree was poorly resolved and with low node support.Therefore, final analyses were performed including all bases.The best-fit models of evolution for each gene were determined using the Akaike information criterion (Akaike 1974) implemented in MrModeltest v. 2.3 (Nylander 2004), resulting in the GTR+I+G model for COI, 16S and H3.Maximum likelihood (ML) analyses were performed using the RA×ML software v7.0.4 (Stamatakis 2006) and node support was assessed with non-parametric bootstrap (BS) with 5000 replicates, random starting trees, and parameters estimated from each dataset under the model selected for the original dataset.Bayesian Inference (BI) analyses were conducted using MrBayes version 3.1.2b(Ronquist and Huelsenbeck 2003) for five million generations with two independent runs and a sampling frequency of 1000.The models implemented were those estimated with MrModeltest v. 2.3.Convergence was diagnosed graphically by plotting for each run the likelihood against the number of generations using the software Tracer version 1.4.1 (Drummond and Rambaut 2007).For each analysis, the first 1250 trees were discarded ('burn-in' period).Node support was assessed with posterior probability (PP).BI and RA×ML phylogenetic analyses were performed on the 280-core "Phylo-Cluster" hosted at the Center for Comparative Genomics, California Academy of Sciences.Only nodes supported by PP≥0.95 and BS≥75 were considered as resolved.The combined tree provided better resolution than COI (658 pb), H3 (328 pb) or 16S (up to 476 pb) separately.The combined dataset yielded a sequence alignment of 1462 positions.The ABGD method (Puillandre et al. 2012) was performed for the COI alignment using the online version of the software (available at http://wwwabi.snv.jussieu.fr/public/abgd/)(18 Nov 2016).ABGD was run by selecting Jukes-Cantor and Kimura parameters distance, Pmin=0.001,Pmax=0.1,Steps=10, and relative gap width=1.

Samples for morphological studies
Two specimens of Felimare were obtained in two different surveys in May 2009 and July 2011 through scuba diving at Tarrafal, Cape Verde, and preserved in 96% ethanol.One specimen was dissected by dorsal incision.The internal features were examined using a dissecting microscope and drawn with a camera lucida.The buccal mass was removed and dissolved in 10% sodium hydroxide until the radula and the labial cuticle were isolated from the surrounding tissue.Both were then rinsed in water, dried and mounted for examination under a Quanta 200 scanning electron microscope.
Morphological and anatomical comparison between the new species and congeners from the eastern Atlantic blue chromodorid chromatic group was based on published information and personal observations.

Molecular analyses
We successfully obtained 64 new sequences; 12 additional sequences were obtained from GenBank (Table 2).The combined tree of COI, 16S and H3 provided better resolution than individual genes trees (Fig. 1).The individual genes trees can be seen in the supplementary material Figure S1. Figure 1 shows the phylogenetic hypothesis based on the combined dataset constructed by BI.The topology of the ML tree was almost identical (not shown).
The maximum intraspecific pairwise uncorrected p-distance for COI between specimens of F. bilineata was 4.23%.The minimum interspecific p-distance within Felimare was between Felimare n. sp. and F. bilineata and was 8.63% (Table 3).Both ABGD species delimitation method analyses recovered eight partitions with prior maximal distance P=0.001 for 1 group and P=0.035 for 8 groups.(Fig. 1).External morphology.Body high and elongate, with dark blue colour and smooth surface.Small round orange spots all over the body, including the foot (Fig. 2).Mantle edge narrow.Dorsum with a series of large and unclustered MDFs on the edge of the mantle, totally absent in the anterior area facing the rhinophores, 15 on the 50 mm specimen (MB28-004391), and 9 on the 30 mm specimen (MB28-004390).MDFs easily seen through the mantle.Posterior end of the foot not covered by the notum (Fig. 2).Rhinophores and branchial leaves dark blue, slightly lighter in colour.Rhinophores of two specimens with 22 lamellae, arranged nearly horizontal.Ten unipinnate branchial leaves in the 50 mm specimen (MB28-004391).Anatomy (Figs 3,4).Jaw composed of two pieces covered by unicuspid rodlets (Fig. 3A, B).The radular formula of the 50-mm specimen (MB28-004391) is 69×114.0.114.(Fig. 3C).Without rachidian radular tooth.Innermost lateral teeth sharp bifid, inner lateral teeth with two small denticles on their outer edge (Fig. 3D).Middle lateral teeth curved and bicuspid (Fig. 3E).Outermost teeth slightly bifid, broad tips, with up to four denticles on outer edge (Fig. 3F).

Taxonomy
Reproductive system (Fig. 4A, B) with vestibular gland elongated.Vagina very wide and muscular.Pyriform seminal receptacle with a short and thin conduct that joins along the middle length of the vagina.Bursa copulatrix spherical.Uterine duct short and narrow, entering the female gland near the entrance of the oviduct.Ampulla elongated and wide.Penis unarmed.Prostate elongated, located anteriorly to female gland, narrowing slightly to a convoluted deferent duct, all over the reproductive system.Female gland half the size of the entire reproductive system.Distribution.To date only known from Tarrafal, Santiago Island, Cape Verde.
Etymology.The name refers to the Latin words aurantiacus, meaning orange, and maculatus, meaning spot, which refers to the chromatic pattern of the body, scattered with orange polka dots.
Comparative diagnosis.Though all the Felimare species distributed along the Atlantic coast of Africa and belonging to the blue Atlantic chromatic group have some chromatic yellow pattern over light to dark blue, Felimare aurantimaculata n. sp. is the only one having small orange dots instead of a white middle line with a distinct and unique pattern as in Felimare pinna; parallel white line or yellow stripes as in Felimare bilineata, Felimare ciminoi (Ortea, Valdés and García-Gómez, 1996), Felimare francoisae (Bouchet in Bouchet and Ortea, 1980), Felimare garciagomezi (Ortea and Valdés, 1996), Felimare tema, Felimare tricolor and Felimare xicoi (Ortea and Valdés, 1996); or irregular patterns as in Felimare gofasi (Ortea, Valdés and García-Gómez, 1996).The morphotype of Felimare picta from the Azores Islands has small yellow spots in a lower density all over the body that make it resemble Felimare aurantimaculata n. sp., but our results show that they are different species nested in different subclades (Fig. 1).The size of the   MDFs of Felimare aurantimaculata sp.nov. is larger than that of the remaining known species of the genus.Moreover, the arrangement of the MDFs in Felimare aurantimaculata n. sp. is also unique among the above blue chromatic group, since the remaining species lack MDFs in the middle region of their mantle edge or they are limited to the posterior end (Table 4).The species F. francoisae and F. garciagomezi (previously named as Mexichromis) differ from Felimare aurantimaculata n. sp., besides the chromatic pattern, in having jaw rodlets with 3-4 cuspids, instead of unicuspid rodlets as in F. ciminoi and F. gofasi (Table 4).The reproductive system is like that of other species of Felimare, formerly ascribed to Hypselodoris, such as the size and shape of the vestibular gland; the muscularized wide vagina; the spherical shape of the bursa copulatrix; and the pyriform seminal receptacle, generally attached in some part of the vagina, with the exception of Felimare molloi, that joins just at the entrance of the bursa copulatrix (Ortea et al. 1996, Gosliner andJohnson 1999).
Recently, using molecular techniques, Almada et al. (2016) and Furfaro et al. (2016) ended the controversy related to the Felimare picta complex.They both stated that Felimare tema (Edmunds, 1981) is a valid species distributed from Cape Verde and Senegal to Angola, in the southern hemisphere, while Felimare picta is restricted to the northern hemisphere, but with a broader distribution, from the type locality in the Mediterranean sea up to the Gulf of Mexico.Although there is now no doubt that F. picta is not distributed in the southeastern Atlantic, we decided to include F. picta in our phylogenetic analyses to discard the possibility that Felimare aurantimaculata n. sp.could be a chromatic morphotype of this species.IB, ML, p-distance and ABGD strongly support the hypothesis that Felimare aurantimaculata n. sp. is a new valid species.

Fig. 1 .
Fig. 1. -Phylogenetic hypothesis based on BI of the combined dataset (H3+COI+16S).Numbers above branches represent PP.Numbers below branches represent BS.Unsupported branches not labelled.Dashed rectangles are ABGD groups (Jukes-Cantor parameter).p-uncorrected distances for COI data set for each species.The new species is within the grey rectangle.

Table 1 .
-Species of Felimare of the blue Atlantic chromatic group distributed along the western African coast.

Table 2 .
-Specimens used for molecular analyses, sample codes, vouchers, collection sites, GenBank accession numbers and references.

Table 3 .
-Minimum and maximum and COI gene pairwise uncorrected p-distances between F. aurantimaculatata n. sp. and the remaining species of the genus Felimare.

Table 4 .
-Comparative table of the main features for Felimare species belonging to the blue Atlantic chromatic group distributed along the western African coast.