Metabarcoding by capture using a single COI probe (MCSP) to identify and quantify fish species in ichthyoplankton swarms.

The ability to determine the composition and relative frequencies of fish species in large ichthyoplankton swarms could have extremely important ecological applications However, this task is currently hampered by methodological limitations. We proposed a new method for Amazonian species based on hybridization capture of the COI gene DNA from a distant species (Danio rerio), absent from our study area (the Amazon basin). The COI sequence of this species is approximately equidistant from all COI of Amazonian species available. By using this sequence as probe we successfully facilitated the simultaneous identification of fish larvae belonging to the order Siluriformes and to the Characiformes represented in our ichthyoplankton samples. Species relative frequencies, estimated by the number of reads, showed almost perfect correlations with true frequencies estimated by a Sanger approach, allowing the development of a quantitative approach. We also proposed a further improvement to a previous protocol, which enables lowering the sequencing effort by 40 times. This new Metabarcoding by Capture using a Single Probe (MCSP) methodology could have important implications for ecology, fisheries management and conservation in fish biodiversity hotspots worldwide. Our approach could easily be extended to other plant and animal taxa.

DNA Metabarcoding of Amazonian Ichthyoplankton Swarms.

Tropical rainforests harbor extraordinary biodiversity. The Amazon basin is thought to hold 30% of all river fish species in the world. Information about the ecology, reproduction, and recruitment of most species is still lacking, thus hampering fisheries management and successful conservation strategies. One of the key understudied issues in the study of population dynamics is recruitment. Fish larval ecology in tropical biomes is still in its infancy owing to identification difficulties. Molecular techniques are very promising tools for the identification of larvae at the species level. However, one of their limits is obtaining individual sequences with large samples of larvae. To facilitate this task, we developed a new method based on the massive parallel sequencing capability of next generation sequencing (NGS) coupled with hybridization capture. We focused on the mitochondrial marker cytochrome oxidase I (COI). The results obtained using the new method were compared with individual larval sequencing. We validated the ability of the method to identify Amazonian catfish larvae at the species level and to estimate the relative abundance of species in batches of larvae. Finally, we applied the method and provided evidence for strong temporal variation in reproductive activity of catfish species in the Ucayalí River in the Peruvian Amazon. This new time and cost effective method enables the acquisition of large datasets, paving the way for a finer understanding of reproductive dynamics and recruitment patterns of tropical fish species, with major implications for fisheries management and conservation.

Temporal and spatial distribution of young Brachyplatystoma spp. (Siluriformes: Pimelodidae) along the rapids stretch of the Madeira River (Brazil) before the construction of two hydroelectric dams.

Monthly (April 2009 to May 2010) bottom-trawl sampling for Brachyplatystoma species along the rapids stretch of the Madeira River in Brazil revealed that Brachyplatystoma rousseauxii larvae and juveniles were present in low abundances in all areas and during all hydrological periods. The presence of larvae and juveniles throughout the hydrological cycle suggests asynchronous spawning in the headwaters of the Madeira River.

Histological development of the digestive system of the Amazonian pimelodid catfish Pseudoplatystoma punctifer.

The organogenesis of the digestive system was described in the Amazonian pimelodid catfish species Pseudoplatystoma punctifer from hatching (3.5 mm total length, TL) to 41 days post-fertilization (dpf) (58.1 mm TL) reared at 28°C. Newly hatched larvae showed a simple digestive tract, which appeared as a straight undifferentiated and unfolded tube lined by a single layer of columnar epithelial cells (future enterocytes). During the endogenous feeding period, comprised between 20 and 96 h post-fertilization (3.5 to 6.1 mm TL), the larval digestive system experienced a fast transformation with the almost complete development and differentiation of most of digestive organs (buccopahrynx, oesophagus, intestine, liver and exocrine pancreas). Yolk reserves were not completely depleted at the onset of exogenous feeding (4 dpf, 6.1 mm TL), and a period of mixed nutrition was observed up to 6 to 7 dpf (6.8 to 7.3 mm TL) when yolk was definitively exhausted. The stomach was the organ that latest achieved its complete differentiation, characterized by the development of abundant gastric glands in the fundic stomach between 10 and 15 dpf (10.9 to 15.8 mm TL) and the formation of the pyloric sphincter at the junction of the pyloric stomach and the anterior intestine at 15 dpf (15.8 mm TL). The above-mentioned morphological and histological features observed suggested the achievement of a digestive system characteristic of P. punctifer juveniles and adults. The ontogeny of the digestive system in P. punctifer followed the same general pattern as in most Siluriform species so far, although some species-specific differences in the timing of differentiation of several digestive structures were noted, which might be related to different reproductive guilds, egg and larval size or even different larval rearing practices. According to present findings on the histological development of the digestive system in P. punctifer, some recommendations regarding the rearing practices of this species are also provided in order to improve the actual larval rearing techniques of this fast-growing Neotropical catfish species.

Genetic structure in the Amazonian catfish Brachyplatystoma rousseauxii: influence of life history strategies.

The Dorado or Plateado (Gilded catfish) Brachyplatystoma rousseauxii (Pimelodidae, Siluriformes) is a commercially valuable migratory catfish performing the largest migration in freshwaters: from the Amazonian headwaters in the Andean foothills (breeding area) to the Amazon estuary (nursery area). In spite of its importance to inform management and conservation efforts, the genetic variability of this species has only recently begun to be studied. The aim of the present work was to determine the population genetic structure of B. rousseauxii in two regions: the Upper Madera Basin (five locations in the Bolivian Amazon) and the Western Amazon Basin (one regional sample from the Uyucalí-Napo-Marañon-Amazon basin, Peru). Length polymorphism at nine microsatellite loci (284 individuals) was used to determine genetic variability and to identify the most probable panmictic units (using a Bayesian approach), after a significant departure from Hardy-Weinberg equilibrium was observed in the overall dataset (Western Amazon + Upper Madera). Bayesian analyses revealed at least three clusters in admixture in the five locations sampled in the Bolivian Amazon, whereas only two of these clusters were observed in the Western Amazon. Considering the migratory behaviour of B. rousseauxii, different life history strategies, including homing, are proposed to explain the cluster distribution. Our results are discussed in the light of the numerous threats to the species survival in the Madera basin, in particular dam and reservoir construction.

Population genetic structure of Cichla pleiozona (Perciformes: Cichlidae) in the Upper Madera basin (Bolivian Amazon): sex-biased dispersal?

This study investigates the population structure of the Tucunaré (Cichla pleiozona) in the Bolivian Amazon (Upper Madera) by using nuclear (EPIC-PCR, 67 individuals) and mitochondrial (Control Region, 41 published and 76 new sequences) DNA analyses, in relation with ecological (water quality: muddy, clear and mix) and geographic factors. Our analyses of both markers showed the highest diversity in clear waters (Yata, Middle and Upper Iténez), and the existence of two populations in muddy waters (Sécure and Ichilo) and one in mix waters (Manuripi). On the other hand, mitochondrial analyses identified three populations in clear waters where nuclear analyses identified a panmictic population. The highest diversity observed in the Yata-Iténez system suggests that an aquatic refuge occurred during the past in this area. The possible explanations for the observed discrepancy between nuclear and mitochondrial markers are discussed, and a sex-biased dispersal seems to be the most plausible hypothesis in the light of the available information and field observations.

Gender determination in the Paiche or Pirarucu (Arapaima gigas) using plasma vitellogenin, 17beta-estradiol, and 11-ketotestosterone levels.

Arapaima gigas is an air-breathing giant fish of Amazonian rivers. Given its great economic and cultural importance, the aquaculture development of this species represents an evident solution to face the decline of wild populations. In captivity, reproduction occurs generally in large earthen ponds where stocks of a few tens of brooders are maintained together at the beginning of the rainy season (December-March in the Peruvian Amazon). Fry production relies on the spontaneous formation of male and female pairs, which build a nest, delimit a territory and guard the offspring for at least 20 days from other congeners and predators. However, as sex determination of A. gigas is not possible by morphological criteria, it is very difficult to optimize reproduction conditions and fry production in each pond, which seriously hampers the culture of this species. This situation prompted us to develop sexing methodologies based on (1) the detection of female specific plasma Vitellogenin (Vtg) using an enzyme immuno assay (EIA), and (2) the determination of plasma 17beta-estradiol and 11-ketotestosterone levels for immature specimens. The Vtg purification was performed by electro-elution after polyacrilamide gel electrophoresis (PAGE) from plasma of 17beta-estradiol treated A. gigas juveniles. Two different Vtg molecules were isolated, (Vtg(1) and Vtg(2)) with 184 and 112 kDa apparent molecular masses, respectively, and two antibodies were raised in rabbits for each Vtg molecule. Adult fish were 100% accurately sexed by Vtg EIA, while 100% of immature fish and 95% of adults were accurately sexed by 17beta-Estradiol and 11-Ketestosterone ratios. We also observed different color pattern development in male and female adult fish (6-year-olds) around the reproductive period.

Life-history characteristics of the large Amazonian migratory catfish Brachyplatystoma rousseauxii in the Iquitos region, Peru.

The main life-history traits of the dorado Brachyplatystoma rousseauxii, a large Amazonian catfish undertaking the largest migration known for a freshwater fish species (from the nursery area in the estuary of the Amazon to the breeding zones in the head waters of the western Amazon basin close to the Andes), were determined from a 5 year sampling of >15,000 specimens in the Peruvian Amazon. The breeding season occurred during the falling and low-water periods, which is hypothesized to be an adaptation to maximize the chances of young stages to reach the estuary. The size at first sexual maturity was slightly larger for females than males, c. 91 and 83 cm standard length (L(S)), respectively. Both males and females reproduce for the first time at >3 years old. The fecundity per spawning event ranged from 481,734 to 1,045,284 oocytes for females weighing 25 and 34 kg, respectively. Seasonal variations of body condition were similar among sexes, but differed between immature specimens that had a higher condition during the low-water period and lower condition during rising waters, and mature individuals that showed the opposite pattern. The growth characteristics were estimated by L(S) frequency analysis. For females, the best fitting models gave a mean birth date in August, during the height of the breeding cycle, with the following von Bertalanffy growth function parameters: L(Sinfinity) = 153.3, K = 0.29 and t(0) =- 0.37 years. For males, the best fitting model gave a mean birth date in July, also during the height of the breeding period, with L(Sinfinity) = 142, K = 0.30 and t(0) =- 0.36 years. At a given age, females were systematically larger than males and the size difference increased with age. The largest females sampled (148 cm L(S)) was 11 years old and the largest male (134 cm L(S)) was 9 years old. The mortality estimates were higher for males total (Z) = 1.34, natural (M) = 0.52 and fishing (F) = 0.82 than for females (Z = 0.98, M = 0.50, F = 0.48). The life-history patterns of B. rousseauxii are discussed in light of the available knowledge about this species and the understanding of its complex life cycle.