Phylogenomic analyses reveals gene flow between populations of the freshwater shrimp Potimirim brasiliana (Caridea, Atyidae) along its wide distribution

Abstract

Potimirim is one of the 40 genera of Atyidae restricted to America, which occurs in coastal freshwater habitats and questions about population status and variability have been emerging. Potimirim brasiliana occurs in Brazil from the northeastern to southeastern region. In order to evaluate the hypothesis of genetic structure among populations, we performed molecular analyses with specimens from all known limit of distribution of the species. The molecular markers used were COI and 28S. Phylogenetic trees were obtained by maximum likelihood and Bayesian analyses, and a haplotype network was obtained based only on COI. We found clear separation between P. brasiliana, P. potimirim, P. glabra and Potimirim sp 2. No pattern of structuration was found among P. brasiliana, but the haplotype network showed geographic pattern of structuration for the congener P. potimirim. The lack of genetic structuration among P. brasiliana can be explained by its life cycle that requires brackish water to complete their larval development. The larvae and juvenile in contact with these habitats can spread through oceanic currents, especially in higher rainfall seasons, maintaining the gene flow. The explanation for the geographical pattern found among P. potimirim is still missing and aspects about its lifecycle and larval development should be investigated.

Key words
Brazilian drainages; COI; connectivity; Decapoda; 28S

INTRODUCTION

The freshwater shrimp’s diversity found in the tropical and subtropical regions of South America is represented by 5 families, 11 genera, 88 species (Magalhães et al. 2016). Among them, Potimirim Holthuis, 1954, is one of the 40 genera of the speciose family Atyidae De Haan, 1849 (De Grave & Fransen 2011) restricted to the American continent. Members of the genus show short rostrum, with teeth only on the ventral edge, antennal and pterigostomial spines, no supraorbital spines, cornea not wider than the ocular peduncles, basis of the pereiopods without exopod, basis of the first three or four pairs of pereiopods with epipod, slender chelae and endopod of the first pair of pereiopods of males different from the females (Holthuis 1954).

The controversial taxonomic background of the genus Potimirim received an important contribution by Torati & Mantelatto (2012) that used combination of alpha morphology and molecular analyses to elucidate some aspects of its evolution and help the identification of species. Many of these previous difficulties on species identification arise from the wide distribution in similar habitats and high intraspecific variability (Müller 1892, Villalobos 1959, Smalley 1963, Page et al. 2008).

According to the most current literature (Torati & Mantelatto 2012, WoRMS 2020), there are five recognized freshwater species in the genus Potimirim: Potimirim americana Guérin-Méneville, 1855, Potimirim mexicana De Saussure, 1857, Potimirim potimirim (Müller, 1881), Potimirim glabra Kingsley, 1878 and Potimirim brasiliana Villalobos, 1959. Other than these, there are two new species pending description (Torati & Mantelatto 2012): Potimirim sp 1 and Potimirim sp 2. All these species are restricted to America and occur along both Pacific and Atlantic coast in Neotropical region (Torati & Mantelatto 2012).

Potimirim brasiliana is a tiny freshwater shrimp, endemic of coastal rivers and streams of Brazil, distributed from the states of Bahia (northeastern region), Espírito Santo, Rio de Janeiro, São Paulo (southeastern region), to Paraná and Santa Catarina in southern region (Villalobos 1959, Barros & Fontoura 1996a, Torati 2009, Torati & Mantelatto 2012). It is found beneath submerged marginal vegetation, hidden under rocks and gravel (Chace & Hobbs 1969, Ramos-Porto & Palácios 1981, Lima et al. 2006, Torati 2009), or in fast waters (Chace & Hobbs 1969, Ramos-Porto & Palácios 1981). This shrimp genus depends on brackish water to complete its lifecycle (Molina 1987) and some abiotic factors can influence the population structure, such as temperature affecting the time of larval development (Hoffmann & Negreiros-Fransozo 2010), and the reproductive patterns (Barros & Fontoura 1996b, Rocha et al. 2013).

Potimirim brasiliana is a gonochoric species (Grilli et al. 2014), and shows sexual dimorphism, in which the females are commonly larger than males, besides the presence of appendix masculina in males (Parker 1992). The seasons of year and geographical location affect the reproductive activity and fertility of this species (e.g. populations of southern Brazil show lower reproductive activity) (Rocha et al. 2013). According to the IUCN (International Union for Conservation of Nature and Natural Resources) criteria, the conservation status of this species has previously been suggested as Least Concern (LC), but additional studies in different areas are necessary to make inferences more precise on this status (Mantelatto et al. 2016).

Despite of the knowledge accumulated during the past years, and the recent molecular phylogenetic picture between species of Potimirim, there is no conclusive study about the genetic relationships between populations of P. brasiliana. Variable molecular markers can provide substantial background on the genetic structure of populations (Palsbøll et al. 2007). COI is a mitochondrial marker with high mutation rate compared with nuclear genes, and it is used to determine intra and interspecific phylogenetic relationships. The gene 28S is a nuclear marker and was chosen for its low mutation rate, being efficient for the study of interspecific relationships (see Timm & Bracken-Grissom 2015 for review on genetic markers).

Due its wide geographical distribution along the Brazilian drainages, we performed more robust genetic analyses in order to elucidate the relationships and to evaluate the hypothesis of genetic structure among the populations of P. brasiliana. In addition to the comparative analysis, genetic information on the congener P. potimirim was furnished.

MATERIALS AND METHODS

Sampling

We analyzed specimens of the genus Potimirim obtained from field collections in Brazil and from the Crustacean Collection of the Department of Biology (CCDB) at the Faculty of Philosophy, Science and Letters at Ribeirão Preto (FFCLRP), University of São Paulo (USP). Specimens of P. brasiliana were obtained from all known limit of its distribution. The collections complied with current applicable state and federal Brazilian laws (SISBIO license to FLM for the collection and genetic analysis of decapods No. 11777-1, 16/09/2007).

We checked the identification of the material based on diagnostic morphological characteristics of the species (Holthuis 1954, Villalobos 1959, Chace & Hobbs 1969, Chace 1972, Melo 2003, Torati & Mantelatto 2012), especially the shape and size of the appendix masculina, rostrum and dactyl of the third and fourth pereiopods.

Molecular Data

Nucleotide sequences were obtained based on protocols described by Schubart et al. (2000), modified by Pileggi & Mantelatto (2010) and Torati & Mantelatto (2012).

We obtained new sequences of Potimirim for two molecular markers: COI and 28S; complementary sequences were retrieved from GenBank [COI: 7 nominated in GenBank as Potimirim glabra but identified as Potimitim sp 2 by Torati & Mantelatto (2012), and 1 of P. glabra. 28S: 1 of P. potimirim and 1 of Potimirim glabra, identified as Potimirim sp 2 by Torati & Mantelatto (2012)]. COI sequences of Atya scabra Leach, 1815 and Micratya poeyi Guérin-Méneville, 1855 and 28S sequences of Jonga serrei Bouvier, 1909, Micratya poeyi, Atya ortmannioides Villalobos, 1956, Atya scabra and Atyoida bisulcata Randall, 1840 (all from GenBank) were used as outgroups (Table I).

Table I
Data of the sequences obtained from GenBank. *Species named as P. glabra in GenBank, but posteriorly identified as Potimirim sp 2.

Species	Locality	GenBank	Gene	Reference
P. glabra	Panama	JF811005.1	COI	Page et al. 2011
P. glabra*	Trinidad and Tobago	JF811012.1	COI	Page et al. 2011
P. glabra*	Trinidad and Tobago	JF811011.1	COI	Page et al. 2011
P. glabra*	Trinidad and Tobago	JF811010.1	COI	Page et al. 2011
P. glabra*	Trinidad and Tobago	JF811009.1	COI	Page et al. 2011
P. glabra*	Trinidad and Tobago	JF811008.1	COI	Page et al. 2011
P. glabra*	Trinidad and Tobago	JF811007.1	COI	Page et al. 2011
P. glabra*	Trinidad and Tobago	JF811006.1	COI	Page et al. 2011
P. potimirim	Panama	FN995614.1	28S	Von Rintelen et al. 2012
P. glabra*	Trinidad and Tobago	FN995613.1	28S	Von Rintelen et al. 2012
Atya scabra	Trinidad and Tobago	JF810990.1	COI	Page et al. 2013
Micratya poeyi	Caribbean	FJ348827.1	COI	Cook et al. 2010
Atya scabra	----	FN995550.1	28S	Von Rintelen et al. 2012
Micratya poeyi	----	FN995595.1	28S	Von Rintelen et al. 2012
Jonga serrei	----	FN995587.1	28S	Von Rintelen et al. 2012
Atya ortmannioides	----	FN995549.1	28S	Von Rintelen et al. 2012
Atyoida bisulcata	----	FN995553.1	28S	Von Rintelen et al. 2012

DNA Extraction

Adult males were selected whenever possible, which helps the identification of the species. Genetic vouchers were deposited at CCDB.

Genomic DNA was extracted from muscle tissue of the abdomen of selected individuals. The tissue was incubated for 24 h in 600 μL of lysis buffer with 200 μL of proteinase K (500 μg/mL) at 55°C, followed by addition of 200 μL of ammonium acetate (7.5M). The sample was centrifugated at 14,000 rpm for 10 minutes at 18°C. The supernatant was collected and transferred to 600 μL of isopropanol. The resulting sample was centrifuged at 14,000 rpm for 10 minutes at 18°C and kept in freezer at −20°C. After 24 to 48 h the sample was centrifuged again at 14,000 rpm for 10 minutes at 18°C and the supernatant discarded. The resulting residue (pellet) was washed with ethanol 70%, dried, and resuspended in 30 μL of TE 1× buffer. The extracted DNA concentration was measured by spectrophotometer (NanoDrop® 2000/2000c).

DNA Amplification

The regions of interest [fragments of the mitochondrial gene cytochrome c oxidase subunit I (COI) – primers LCO1490 (5´ - GGT CAA CAA ATC ATA AAG ATA TTG G - 3´) and HCO2198 (5´ - TAA ACT TCA GGG TGA CCA AAA AAT CA - 3´ (Folmer et al. 1994) – and nuclear gene 28S – primes 28S Rd1a (5’ - CCC SCG TAA YTT AAG CAT AT - 3’) and 28S Rd4b (5’ - CCT TGG TCC GTG TTT CAA GAC - 3’) (Edgecombe & Giribet 2006)] – were amplified by PCR (Polymerase Chain Reaction). The PCR products were obtained in reactions of 25 μL containing 6.5 μL of deionized water, 5 μL of 5 M betaine, 3 μL of 10× PCR buffer, 3 μL of 25 mM MgCl2, 4 μL of 0.25 μM each DNTP, 2 μL of primer (1 μL of forward and 1 μL of reverse), 0.5 μL of Taq DNA recombinant polymerase (5 U/μL), and 1 μL of DNA at 10 ng/μL. The amplification was performed in an Applied Biosystems Veriti 96 Well Thermal Cycler® (thermal cycles: initial denaturing for 2 min at 94°C; pairing for 40 cycles: 30 s at 94°C, 30 s at 48°C, 1 min at 72°C; final extension 10 min at 72°C for COI; and initial denaturing for 5 min at 95°C; pairing for 40 cycles: 45 s at 95°C, 45 s at 52°C, 1 min at 72°C; final extension 3 min at 72°C for 28S). Electrophoresis was performed using 1% agarose gel with PCR product. We photographed the results with a digital camera C-7070 Olympus® in a transilluminator UV M20 UVP®.

Purification and Sequencing

We purified the PCR products using the SureClean Plus® kit, following the manufacturer’s instructions. The sequencing of the samples was performed in automated sequencer (ABI 3730 XL DNA Analyzer®, Applied Biosystems) at Department of Technology of the Faculty of Agrarian and Veterinary Sciences at Jaboticabal, Paulista State University “Júlio de Mesquita Filho,” with the reaction kit ABI Big Dye® Terminator Mix (Applied Biosystems).

Editing Sequences

Both DNA strands (sense and antisense) were sequenced for higher reliability. Sequences were edited to generate a consensus sequence using the computational program BioEdit v7.0.5 (Hall 2005). The sequence alignment was conducted in the software MAFFT v7.158 (Katoh & Standley 2013), using the default setting. Identical sequences of the same locality were excluded from analyses.

Phylogenetic Analysis

Statistical selection of nucleotide substitution models was conducted in jModelTest v2.1.4 (Darriba et al. 2012). Bayesian analyses were conducted in MrBayes v3.2.2 program (Ronquist et al. 2012), using parameters (frequency of nucleotides and transition/transversion ratio, gamma distribution and variable sites proportion, when it was present in the model) chosen from among the models available in jModelTest. Analyses were performed with 10⁹ generations in two independent parallel simulations with four chains each and were stopped when they reached stationarity (average standard deviation between two simulations less than 0.01). The parameter values were saved once every 1000 rounds. The first quarter of the parameters and trees were discarded. Only nodes with a posteriori probability higher or equal to 50% were shown (see Ronquist et al. 2012). Initial analyses were performed using COI and 28S matrices separately. Analysis of the concatenate matrix of both genes of interest was also conducted in MrBayes program. Data were divided into two independent partitions, each following the model that best applied to each gene. Maximum likelihood analyses were conducted in the RAxML program (7.6.3) (Stamatakis 2006), implemented in Cyberinfrastructure for Phylogenetic Research (CIPRES; available on line at http://www.phylo.org). The assumed nucleotide substitution model was GTR+Γ+I. The consistency of the topologies was measured by bootstrap method (1000 pseudoreplicates), and only confidence values higher than 50% were reported. The semi-strict consensus tree between Bayesian inference method and maximum likelihood was obtained using the software Mesquite v3.01 (Maddison 2003).

Haplotype Network

The number of haplotypes was obtained in DnaSP v.5.10.1 (Rozas & Rozas 1999). The haplotype network was obtained using the method median joining network in the software PopArt (Population Analysis with Reticulate Trees) (Leigh & Bryant 2015). A GeoTags block was used to establish the geographical coordinates were each specimen was collected and to color the network according to them. Only COI sequences were used in this analysis.

RESULTS

We analyzed the greatest possible diversity of specimens throughout their distribution, as detailed below. In the list of analyzed material there are the specimens that served as genetic vouchers, of which sequences were generated by us, deposited in GenBank and used in the phylogenetic analyzes. Cols. = collectors; ♂ = male; ♀ = female; ♀♀ = ovigerous female.

Potimirim brasiliana

Analyzed material. Uruçuca, Bahia, Brazil, 08/xi/2011, Cols. Carvalho, F.L. et al., 8 ♀, 1 ♀♀, COI – GenBank Access KP202789-93, 28S – KP202727-31, KP202753-56 (CCDB 1585); Ilhéus, Bahia, Brazil, 31/iii/2009, Cols. Mantelatto, F.L. and Almeida, A.O., 7 ♀, COI – KP202794-97, KP202821, 28S – KP202732, KP202757-58 (CCDB 2633); Mangaratiba, Rio de Janeiro, Brazil, 28/ii/2008, Col. Lima, G.V., 6 ♂, 2 ♀, 1 ♀♀, COI – KP202811-14, KP202819-20, 28S – KP202736, KP202743-45, KP202759-61 (CCDB 2142); Paraty, Rio de Janeiro, Brazil, 15/viii/2007, Cols. Mantelatto, F.L. et al., 3 ♂, 3 ♀, COI – KP202786, KP202802-05, 28S – KP202725, KP202737-39, KP202746 (CCDB 2005); Paraty, Rio de Janeiro, Brazil, 12/xii/1996, Cols. unknown, 3 ♂, 1 ♀, 2 ♀♀ (CCDB 2693); Ubatuba, São Paulo, Brazil, 16/viii/2007, Cols. Mantelatto, F.L. et al., 3 ♂, 5 ♀, COI – KP202781, KP202787, KP202806, KP202815-17, 28S – KP202724, KP202740-41, KP202747, KP202749-51, KP202784 (CCDB 2006); Caraguatatuba, São Paulo, Brazil, 09/xi/2007, Cols. Mantelatto, F.L., 2 ♀ (CCDB 0049); São Sebastião, São Paulo, Brazil, 11/vii/2006, Cols. Mantelatto, F.L. et al., 3 ♂, 7 ♀ (CCDB 3105); Ilhabela, São Paulo, Brazil, 13/vii/2006, Col. Mossolin, E.C., 2 ♂, COI – KP202785, 28S – KP202724 (CCDB 1991); Ilhabela, São Paulo, Brazil, 14/vii/2006, Col. Mossolin, E.C., 5 ♂, 2 ♀ (CCDB 6404); Iguape, São Paulo, Brazil, 13/v/2006, Cols. Mantelatto, F.L. et al., 1 ♂, 1 ♀, COI – KP202808, 28S – KP202752 (CCDB 2386); Cananéia, São Paulo, Brazil, 18/iv/2011, Col. Mantelatto, F.L., 4 ♀, 2 ♀♀ (CCDB 3213); Matinhos, Paraná, Brazil, 20/02/2008, Cols. Mantelatto, F.L. & Mossolin, E.C., 3 ♂, COI – KP202788, KP202807, KP202818 (CCDB 2115); BR376, Santa Catarina, Brazil, 19/ii/2008, Col. Mantelatto, F.L. and Mossolin, E.C., 1 ♂, 3 ♀, COI – KP202798-99, KP202809-10, 28S – KP202733-35, KP202742 (CCDB 2116).

Potimirim potimirim

Analyzed material. Isla Colon, Bocas del Toro, Panama, 13/viii/2008, Cols. Torati, L.S. and Page, T., 3 ♂, 3 ♀, COI – KP202783, KP202827-29, 28S – KP202764, KP202777-80 (CCDB 2399); Iconha, Espírito Santo, Brazil, 19/vi/2012, Cols. Carvalho, F.L. et al., 5 ♂, 3 ♀, 1 3 ♀♀, COI – KP202801, KP202822-24, 28S – KP202766-71 (CCDB 4066); Piúma, Espírito Santo, Brazil, 15/vi/2012, Col. Carvalho, F.L. et al., 2 ♀♀, COI – KP202825-26 (CCDB 4069); Ilhéus, Bahia, Brazil, 31/iii/2009, Cols. Mantelatto, F.L. and Almeida, A.O., 1 ♂, 28S – KP202765 (CCDB 2634); Canavieiras, Bahia, Brazil, 17/viii/2010, Col. Carvalho, F.L., 10 ♂, 3 ♀, 28S – KP202772 (CCDB 3037); Ilhabela, São Paulo, Brazil, 07/vii/2011, Cols. Mantelatto, F.L. et al., 5 ♀, 28S – KP202775 (CCDB 3942); Cananéia, São Paulo, Brazil, 29/viii/2011, Cols. Mantelatto, F.L. et al., 15 ♂, 35 ♀, 28S – KP202773 (CCDB 3720); Cananéia, São Paulo, Brazil, 29/viii/2011, Cols. Mantelatto, F.L. et al., 1 ♀, 28S – KP202763 (CCDB 3721); Cananéia, São Paulo, Brazil, 09/xi/2011, Cols. Carvalho, F.L. et al., 1 ♀, COI – KP202784, 28S – KP202762 (CCDB 3743); Ariri, São Paulo, Brazil, 10/xi/2011, Cols. Carvalho, F.L. et al., 15 ♂, 6 ♀, COI – KP202800, 28S – KP202774 (CCDB 3741); Caieira do Norte, Santa Catarina, Brazil, 18/iv/2007, Cols. Mantelatto, F.L. et al., 3 ♂, 3 ♀, 28S – KP202776 (CCDB 1893).

Potimirim glabra

Analyzed material. Pacific Ocean, Costa Rica, 10/ii/2009, Cols. Mantelatto, F.L. et al., 7 ♂, 1 ♀, COI – KP202782 (CCDB 3341).

We obtained 56 sequences of ~750 base pairs (bp) of Potimirim for 28S (37 of P. brasiliana and 19 of P. potimirim). For COI, we obtained 49 sequences of ~530 bp (36 of P. brasiliana, 1 of P. glabra and 12 of P. potimirim). The number of haplotypes of each species was: 13 of P. brasiliana, 1 of P. glabra, 5 of Potimirim sp 2 and 7 of P. potimirim, including the new sequences and the ones retrieved from GenBank.

The analysis conducted in the software jModelTest showed that the substitution models that best fit our data were HKY+I for 28S and HKY+ Γ for COI, considering both Akaike Information Criterion (AIC) and Bayesian Information Criterion (BIC).

The consensus tree between the Bayesian inference and maximum likelihood, based on the concatenate matrix of both genes (fig. 1), showed a clear separation between Potimirim potimirim, Potimirim glabra and the other species. Potimirim sp 2 is also a well-supported clade, but among a polytomy of Potimirim brasiliana.

Figure 1
Consensus tree between the maximum likelihood and Bayesian inference based on the concatenate matrix of COI and 28S. Node numbers represent the bootstrap values and a posteriori probability. *Sequence from GenBank. Brackets show Potimirim species that were retrieved as monophyletic groups by the analysis. BA – Bahia, ES – Espírito Santo, PR – Paraná, RJ – Rio de Janeiro, SC – Santa Catarina, SP – São Paulo.

Separately, the nuclear gene 28S was inefficient for the study of phylogenetic relationships among Potimirim. The trees obtained based on this molecular marker showed very poor resolution of the clades and only showed the separation of P. potimirim from other species. On the other hand, the mitochondrial gene COI showed clear delimitations between the studied species. The consensus tree between the Bayesian inference and maximum likelihood, based only on COI (fig. 2), showed that each one of the studied species form a distinct clade, including Potimirim brasiliana.

Figure 2
Consensus tree between the maximum likelihood and Bayesian inference based on COI. Node numbers represent the bootstrap values and a posteriori probability. *Sequence from GenBank. BA – Bahia, ES – Espírito Santo, PR – Paraná, RJ – Rio de Janeiro, SC – Santa Catarina, SP – São Paulo.

Although Potimirim brasiliana has a wide distribution along the Brazilian drainage, there is no geographical pattern of structuration among the populations, what is shown by the tree and the haplotype network (figs. 2 and 3) based on COI. Nevertheless, the haplotype network showed that there is a geographical pattern of distribution of the haplotypes of Potimirim potimirim (i.e. Panama – São Paulo/Parana – Espirito Santo). P. brasiliana’s COI intraspecific distances ranged from 0.0 to 8.6%. Interspecific distance between Potimirim species ranged from 12.8 to 16.9%.

Figure 3
Haplotype network based on COI (530 bp) for populations of Potimirim and two outgroups (Atya scabra and Micratya poeyi).

DISCUSSION

Potimirim brasiliana’s COI intraspecific distances (0.0 to 8.6%) were lower than the mean interspecific gap of Atyidae shrimps (13.2 to 13.8%) calculated by Silva et al. (2011), and lower than the interspecific distance between Potimirim species (12.8 to 16.9%) calculated in the present study. In addition, we found 13 haplotypes of COI among P. brasiliana, in which the topology showed no geographical pattern, which means that the populations present gene flow between them. The explanation about this intriguing genetic connection even within smaller scales is not totally conclusive since it was expected that species with wide geographic distribution in association with the differences in the environmental conditions, presents expressive genetic divergence. In this sense and no exclusive factor, we found some characteristics of larval cycle of atyids to support our hypothesis. Although P. brasiliana is an amphidromic shrimp, in which adult individuals live and reproduce in freshwater, the larvae of Potimirim species require brackish water to successfully complete their larval development and lifecycle (Villalobos 1959, Molina 1987, Hoffmann & Negreiros-Fransozo 2010). In most its area of occurrence, P. brasiliana shows reproductive activity throughout the whole year, with more intensity during the summer, which is coincident with the higher rainfall period in southern Atlantic region and that are favorable to the passive transport of larvae to the estuary and their development in different salinity (see Molina 1987, Hoffmann & Negreiros-Fransozo 2010, Rocha et al. 2013 for review). Therefore, it is likely that the gene flow occurs when these potential colonizers (larvae) are in contact with these habitats in conditions which the dispersion of larvae or juvenile by ocean currents in facilitated. This is considered the main explanation of the widespread distribution and genetic connectivity seen in amphidromous species (Page et al. 2008, Fujita et al. 2016).

As observed in some other members of Atyidae (Molina 1987, Hoffmann & Negreiros-Fransozo 2010, Rocha et al. 2013), Potimirim brasiliana presents many favorable and convergent characteristics that fits with this model of dispersion: 1) produces many small elliptic eggs per brood (mean fecundity = 515 eggs/female; mean diameter = 0.57 mm x 0.37 mm; Molina 1987); 2) high reproductive activity throughout the year (Hoffmann & Negreiros-Fransozo 2010, Rocha et al. 2013); 3) a moderate incubation period (extends from 14 to 27 days according to the temperature) (Molina 1987), and 4) extended larval developmental period (11 zoea and 1 decapodid stages reached in about 45 days (14-16 salinities range) under laboratory conditions (Molina 1987). These population and larval cycle profiles are key features to the success of this dispersion in low scales and gene flow.

In contrast, the haplotype network and the tree based only on COI showed geographical pattern of separation for populations of Potimirim potimirim, divided into populations from Panama, São Paulo and Espírito Santo. Unfortunately, as far as we know, nothing is known about larval cycle of P. potimirim which do not permit us make a conclusive inference regarding the capacity of dispersion and the absence of connection among the studied populations. Considering the phylogenetic relationship between both species (see below more details), studies on larval development and experiments on survival in different salinities of P. potimirim could help to explain why there is gene flow between populations of P. brasiliana. They would also explain why this does not occur in P. potimirim and allow to test the idea that those populations might be experiencing different selective pressures during the evolution in America. In addition, due the phylogenetic separation in two different clades (see Torati & Mantelatto 2012 and discussion below), we can hypothesize that P. potimirim much probably present a different pattern of larval development that should be investigated.

The mitochondrial marker COI showed a clear separation of the four studied species of Potimirim, corroborated Torati & Mantelatto (2012) and reinforced its utility, as a highly variable mitochondrial molecular marker, to delimit species and even populations (Harrison 2004, Toon et al. 2009). On the other hand, the nuclear gene 28S was not efficient to show the relationships between populations of P. brasiliana but showed a clear separation between this species and P. potimirim. In general, nuclear markers are well conserved and are normally used to delimit genera (Toon et al. 2009).

The evolution of the genus Potimirim is strongly related to the closure of the Isthmus of Panama (Torati 2009). This biogeographic barrier was responsible for separating Potimirim populations from the Pacific Ocean from those in the Atlantic Ocean, which made possible the speciation of these populations. According to Torati (2009), the cladogenesis of P. brasiliana and Potimirim sp 2 occurred more recently than the cladogenesis of P. glabra and the ancestral population of P. brasiliana and Potimirim sp 2, and our results corroborate this once the trees shows that Potimirim sp 2 is more closely related to P. brasiliana than to P. glabra.

Although P. glabra and Potimirim sp 2 are morphologically similar (Torati & Mantelatto 2012), our molecular analyses showed that Potimirim sp 2 is more closely related to P. brasiliana. These three species share an important morphological character: the shape of the appendix masculina, which distances them from P. potimirim. Torati & Mantelatto (2012; Fig. 3, page 635) obtained two groups of Potimirim supported by both the 16Smt topology and morphology of the appendix masculina: one formed by P. potimirim, P. mexicana and Potimirim sp 1, and the other containing P. brasiliana, P. glabra and Potimirim sp 2. This separation, however, was not supported by Grilli et al. (2014), considering the sexual system of the species, once P. mexicana is protrandric hermaphrodite and P. potimirim is gonochoric, a difference that should be investigated.

Our findings provide a new perspective about the evolution of the genus Potimirim, especially regarding the life cycle and gene flow particularities observed in the studied species. This scenario serve as appropriate theoretical basis for the development of management strategies for conservation proposes since P. brasiliana is endemic to Brazilian drainages and the genetic connectivity found within the geographical regions is important to maintain the genetic variability and prevent local extinction, as observed in some invertebrates (Bohonak & Jenkins 2003), including mangrove crab (Buranelli et al. 2019). Additionally, our findings highlighted the need of studies on larval development of Potimirim potimirim to fully understand its genetic structuration and the taxonomic status.

ACKNOWLEDGMENTS

The present study is part of a long-term project to evaluate the taxonomy of freshwater/estuarine decapods in Brazil, and was supported by scientific grants provided to FLM by the Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP, 2002/08178-9; PhD to LGP 2005/50651-1; Temáticos Biota 2010/50188-8 and INTERCRUSTA 2018/13685-5; Coleções Científicas 2009/54931-0; PROTAX 2016/50376-5) and the Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq – Brazil, 473050/2007-2, 471011/2011-8; 490314/2011; PQ 302748/2010-5, 304968/2014-5 and 302253/2019-0). BMP received support from University of São Paulo (Bolsa Institucional RUSP – Procs. 2012/2741 and 2013/1803). LGP was supported by post-doctoral fellowship (PNPD Proc. 02630/09-5) and ALVS by Master Science fellowship (Proc. 88882.158755/2017-01), both from Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES - Código de Financiamento 001, Ciências do Mar II Proc. 2005/2014 − 23038.004308/2014-14). We are grateful to many colleagues and friends (Alexandre Almeida, Edvanda Souza-Carvalho, Emerson Mossolin, Fabricio Carvalho, Lucas Torati, Ingo Wehrtmann) for their help in collections, for making available some essential fresh specimens, for lending material from collections used in our research, and for critical discussion during the preparation of this manuscript. Thanks are due to all members of LBSC for their assistance during the development of this study. We are thankful to the anonymous reviewers for their suggestions that helped improve the quality of the manuscript.

REFERENCES

• BARROS MP & FONTOURA NF. 1996a. Biologia reprodutiva de Potimirim glabra (Kingsley, 1878) (Crustacea, Decapoda, Atyidae), na praia da Vigia, Garopaba, Santa Catarina, Brasil. Nauplius 4: 1-10.
• Barros MP & Fontoura NF. 1996b. Crescimento de Potimirim glabra (Kingsley, 1878) (Crustacea, Decapoda, Atyidae), na Praia da Vigia, Garopaba, Santa Catarina, Brasil. Nauplius 4: 11-28.
• Bohonak AJ & Jenkins DG. 2003. Ecological and evolutionary significance of dispersal by freshwater invertebrates. Ecol Letters 6: 783-796.
• Bouvier MEL. 1909. Les Crevettes d’eau Douce de la Famille des Atyidés qui se Trouvent dans l’île de Cuba. Bull Mus Nat D’Hist Nat 15(6): 329-336.
• BURANELLI RC, FELDER DL & MANTELATTO FL. 2019. Genetic diversity among populations of the Western Atlantic mangrove crab Ucides cordatus (Linnaeus, 1763) (Decapoda: Brachyura: Ocypodidae): evidence for panmixia and useful data for future management and conservation. J Crust Biol 39(4): 386-395.
• Chace FA & Hobbs HH. 1969. The freshwater and terrestrial Decapod Crustaceans of the West Indies with special reference to Dominica. U S Nat Mus Bull 292: 1-258.
• Chace FA. 1972. The shrimps of the Smithsonian-Bredin Caribbean Expeditions with a summary of the West Indian. Shallow-water species (Crustacea: Decapoda: Natantia). Smith Contrib Zool 98: 1-179.
• Cook BD, Pringle CM & Hughes JM. 2010. Immigration history of amphidromous species on a Greater Antillean Island. J Biogeogr 37(2): 270-277.
• Darriba D, Taboada GL, Doallo R & Posada D. 2012. jModeltest 2: more models, new heuristics and parallel computing. Nat Meth 9(8): 772.
• De Grave S & Fransen CHJM. 2011. Carideorum catalogus: the recent species of the dendrobranchiate, stenopodidean, procaridean and caridean shrimps (Crustacea: Decapoda). Zool Med 89: 195-589.
• De Haan W. 1849. Crustacea elaborate. In: Von Siebold PF (Ed), Fauna Japonica. Vol. 6, xxi + 230, p. 184-186.
• De Saussure MH. 1857. Memoire sur Divers Crustaces Nouveaux des Antilles et du Mexique. Mem Soc Phys d’Hist Nat Gen 14(2): 417-496.
• Edgecombe GD & Giribet G. 2006. A century later – a total evidence re-evaluation of the phylogeny of scutigeromorph centipedes (Myriapoda: Chilopoda). Invert Syst 20: 503-525.
• Folmer M, Black W, Hoeh RL & Vrijenhoek R. 1994. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol Mar Biol Biotech 3(5): 294-299.
• Fujita S, Zenimoto K, Iguchi A, Kai Y, Ueno M & Yamashita S. 2016. Comparative phylogeography to test for predictions of marine larval dispersal in three amphidromous shrimps. Mar Ecol Progr Ser 560: 105-120.
• Grilli NM, Terossi M & Mantelatto FL. 2014. Sexual system of the freshwater shrimps of the genus Potimirim Holthuis (Decapoda: Caridea: Atyidae): is there a pattern in this genus? Mar Fresh Res 65: 759-765.
• Guérin-Méneville FE. 1855. Animales Articulados con piés articulados. Crustaceos, Parte Descriptiva. In: De La Sagra DR (Ed), Historia fisica, politica y natural de la Isla de Cuba. Librero de la Sociedad de Geografia, Paris.
• Hall T. 2005. BioEdit 7.0.5. North Carolina State University, Department of Microbiology.
• Harrison JS. 2004. Evolution, biogeography, and the utility of mitochondrial 16s and COI genes in phylogenetic analysis of the crab genus Austinixa (Decapoda: Pinnotheridae). Mol Phylog Evol 30: 743-754.
• Hoffmann P & Negreiros-Fransozo ML. 2010. Reproductive cycle and fecundity of Potimirim glabra (Kingsley, 1954) (Caridea, Atyidae) from a littoral stream. Invert Reprod Develop 54(3): 133-141.
• Holthuis LB. 1954. On a collection of decapod Crustacea from the Republic of El Salvador (Central America). Zool Verh 23: 1-43.
• Katoh K & Standley DM. 2013. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol 30(4): 772-780.
• Kingsley JS. 1878. Notes on the North American Caridea in the Museum of the Peabody Academy of Science at Salem. Proc Acad Nat Sci Philad 1878: 89-98.
• Leach WE. 1815. A tabular view of the external characters of four classes of animals, which Linné arrange under Insecta, with the descriptions of the genera composing three of these classes into orders, and c. and descriptions of several new genera and species. Trans Linn Soc Lond 11: 306-400.
• Leigh JW & Bryant D. 2015. PopART: Full-feature software for haplotype network construction. Meth Ecol Evol 6(9): 1110-1116.
• Lima GV, Silveira CM & Oshiro LMY. 2006. Estrutura populacional dos camarões simpátricos Potimirim glabra e Potimirim potimirim (Crustacea, Decapoda, Atyidae) no rio Sahy, Rio de Janeiro, Brasil. Iheringia, Ser Zool 96(1): 81-87.
• Maddison DR. 2003. Are strepsipterans related to flies? Exploring longbranch attraction. Study 2. In: Mesquite: a modular system for evolutionary analysis, version 1.0. Available at: http://mesquiteproject.org Last access: 21 October 2018.
  » http://mesquiteproject.org
• Magalhães C, Campos MR, Collins PA & Mantelatto FL. 2016. Chapter 10: Diversity, distribution and conservation of freshwater crabs and shrimps in South America, p. 303-322. In: Kawai T & Cumberlidge N (Eds), A Global Overview of the Conservation of Freshwater Decapod Crustaceans. Springer, Cham, Switzerland.
• Mantelatto FL, Torati LS, Pileggi LG, Mossolin EC, Terossi M, Carvalho FL, Rocha SS & Magalhães C. 2016. Avaliação dos Camarões Atiídeos (Decapoda: Atyidae), Cap. 5: p. 93-102. In: Pinheiro M & Boos H (Eds), Livro Vermelho dos Crustáceos do Brasil: Avaliação 2010-2014. Porto Alegre, RS, Sociedade Brasileira de Carcinologia - SBC, 466 p.
• Melo GAS. 2003. Famílias Atyidae, Palaemonidae e Sergestidae. In: Melo GAS (Ed), Manual de Identificação dos Crustacea Decapoda de Água Doce do Brasil, p. 289-415.
• Molina FML. 1987. Biologia de Potimirim brasiliana Villalobos, 1950 (Crustacea, Decapoda, Atyidae). PhD Thesis, Universidade de São Paulo. São Paulo, Brasil.
• Müller F. 1881. Atyoida potimirim, eine schlamm fressende Süsswassergarneele. Kosmos 9: 117-124.
• Müller F. 1892. O camarão miúdo do Itajahy, Atyoida potimirim. Arch Mus Nac Rio de Janeiro 8: 155-178.
• Page TJ, Cook BD, Rintelen T, Rintelen K & Hughes JM. 2008. Evolutionary relationships of atyid shrimps imply both ancient Caribbean radiations and common marine dispersals. J North Amer Benthol Soc 27(1): 68-83.
• Page TJ, Cook BD, Torati LS, Binderup A, Pringle CM & Hughes JM. 2011. Direct Submission. Australian Rivers Institute, Griffith University, Molecular Lab, Nathan, QLD 4111, Australia.
• Page TJ, Torati LS, Cook BD, Binderup A, Pringle CM, Reuschel S, Schubart CD & Hughes JM. 2013. Invertebres Sans Frontieres: Large Scales of Connectivity of Selected Freshwater Species among Caribbean Islands. Biotropica 45: 236-244.
• Palsbøll PJ, Berube M & Allendorf FW. 2007. Identification of management units using population genetic data. Trends Ecol Evol 22: 11–16.
• Parker GA. 1992. The evolution of sexual dimorphism in fish. J Fish Biol 41(B): 1-20.
• Pileggi LG & Mantelatto FL. 2010. Molecular phylogeny of the freshwater prawn genus Macrobrachium (Decapoda, Palaemonidae), with emphasis on the relationships among selected American species. Invert Syst 24(2): 194-208.
• Ramos-Porto M & Palácios JAP. 1981. Estudos ecológicos do rio Capibaribe-Mirim, PE. IV. Crustáceos decápodas natantes. Trab Oceanogr UFPE 16: 265-296.
• Randall JW. 1840. Catalogue of the Crustacea brought by Thomas Nuttall and J.K. Townsend, from the west coast of North America and the Sandwich Islands, with 88 descriptions of such species as are apparently new, among which are included several species of different localities, previously existing in the collection of the Academy. J Acad Nat Sci Philad 8: 106-147.
• Rocha SS, Bueno SLS, Shimizu RM & Mantelatto FL. 2013. Reproductive biology and population structure of Potimirim brasiliana Villalobos, 1959 (Decapoda, Atyidae) from a littoral fast-flowing stream, São Paulo state, Brazil. Crustaceana 86(1): 67-83.
• Ronquist F, Teslenko M, Van Der Mark P, Ayres D, Darling A, Höhna S, Larget B, Lui L, Suchard MA & Huelsenbeck JP. 2012. MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst Biol 61: 539-542.
• Rozas J & Rozas R. 1999. DnaSP version 3.0: an integrated program for molecular population genetics and molecular evolution analysis. Bioinformatics 15(2): 174-175.
• Schubart CD, Neigel JE & Felder DL. 2000. Use of the mitochondrial 16S rRNA gene for phylogenetic and population studies of Crustacea. Crust Issues 12: 817-830.
• Silva MJ, Creer S, dos Santos A, Costa AC, Cunha MR, Costa FO & Carvalho GR. 2011. Systematic and evolutionary insights derived from mtDNA COI barcode diversity in the Decapoda (Crustacea: Malacostraca). PLoS ONE 6(5): e19449.
• Smalley AE. 1963. The genus Potimirim in Central America (Crustacea, Atyidae). Rev Biol Trop 11(2): 177-183.
• Stamatakis A. 2006. RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics 22(21): 2688-2690.
• Timm L & Bracken-Grissom HD. 2015. The forest for the trees: evaluating molecular phylogenies with an emphasis on higher-level Decapoda. J Crust Biol 35: 577-592.
• Toon A, Finley M, Staples J & Crandall KA. 2009. Decapod Phylogenetics and Molecular Evolution. Crust Issues 18: 15-29.
• Torati LS. 2009. Revisão taxonômica das espécies brasileiras de Potimirim Holthuis, 1954 e filogenia do grupo baseado em dados moleculares. Dissertação de Mestrado, Universidade de São Paulo. São Paulo, Brasil. (Unpublished).
• Torati LS & Mantelatto FL. 2012. Ontogenetic and evolutionary change of external morphology of the neotropical shrimp Potimirim (Holthuis, 1954) explained by a molecular phylogeny of the genus. J Crust Biol 32(4): 625-640.
• Villalobos FA. 1956. Contribución al conocimiento de los Atyidae de Mexico, I: Una nueva especie de Atya de la vertiente del Pacifico del Estado de Michoacan. Ann Inst Biol UNAM 26(2): 459-475.
• Villalobos FA. 1959. Contribución al conocimiento de los Atyidae de México. II (Crustacea, Decapoda). Estudio de algunas especies del género Potimirim (= Ortmannia), com descripción de una especie nueva en Brasil. Ann Inst Biol UNAM 30: 269-330.
• Von Rintelen K, Page TJ, Cai Y, Roe K, Stelbrink B, Kuhajda BR, Iliffe TM, Hughes J & Von Rintelen T. 2012. Drawn to the dark side: A molecular phylogeny of freshwater shrimps (Crustacea: Decapoda: Caridea: Atyidae) reveals frequent cave invasions and challenges current taxonomic hypotheses. Mol Phylog Evol 36: 82-96.
• WORMS. 2020. Potimirim Holthuis, 1954. Accessed at: http://www.marinespecies.org/aphia.php?p=taxdetails&id=586694 on 2020-12-20.
  » http://www.marinespecies.org/aphia.php?p=taxdetails&id=586694

Publication Dates

• Publication in this collection
  15 Mar 2021
• Date of issue
  2021

History

• Received
  30 Mar 2019
• Accepted
  19 Sept 2019