key: cord-0429069-00icjjip authors: De León, Luis F.; Cornejo, Aydeé; Gavilán, Ronnie G.; Aguilar, Celestino title: Hidden biodiversity in Neotropical streams: DNA barcoding uncovers high endemicity of freshwater macroinvertebrates at small spatial scales date: 2020-03-31 journal: bioRxiv DOI: 10.1101/2020.03.31.018457 sha: 204020507c227cf6f02ed17eaa9916099ff42c5e doc_id: 429069 cord_uid: 00icjjip Aquatic macroinvertebrates play a crucial role in freshwater ecosystems, but their diversity remains poorly known, particularly in the tropics. This “taxonomic void” represents limits our understanding of biodiversity patterns and processes in freshwater ecosystems, and the scale at which they operate. We used DNA barcoding to estimate lineage diversity (and the diversity of unique haplotypes) in 224 specimens of freshwater macroinvertebrates at a small spatial scale within the Panama Canal Watershed (PCW). In addition, we compiled available barcoding data to assess macroinvertebrate diversity at a broader spatial scale spanning the Isthmus of Panama. Consistently across two species delimitation algorithms (i.e., ABGD and GMYC), we found high lineage diversity within the PCW, with ~ 100-106 molecular operational taxonomic units (MOTUs at 2% sequence divergence) across 168 unique haplotypes. We also found a high lineage diversity along the Isthmus of Panama, but this diversity peaked within the PCW. However, our rarefaction/extrapolation approach showed that this diversity remains under sampled. As expected, these results indicate that the diversity of Neotropical freshwater macroinvertebrates is higher than previously thought, with the possibility of high endemicity even at narrow spatial scales. Geographic isolation is likely a main factor shaping these patterns of diversity. However, local disturbances such as the rupture of the continental divide due to the construction of the Panama Canal might be reshaping these patterns of diversity at a local scale. Although further work is needed to better understand the processes driving diversification in freshwater macroinvertebrates, we suggest that Neotropical streams represent continental islands of diversity. Understanding these islands of diversity is crucial in the face of increasing human disturbance. 6 115 agarose gel and purified with the enzymes Exonuclease I (EXO) and Shrimp Alkaline 116 Phosphatase (SAP) [43] , before sequencing using an Applied Biosystems Genetic Analyzer (ABI 117 3130xl, Applied Biosystems, Carlsbad, California). Sequences were aligned using Geneious 118 V7.03 [44] . Sequence alignments were also inspected by eye to confirm overall sequence quality. 119 Project sequences are available in GenBank (accession numbers: KX039451-KX039650, 120 KU980966-KU981004). 121 122 Data analysis 123 To confirm morphological identification for our sequenced specimens, we performed 168 were used as input. All analyses were run using a relative barcoding gap width (X value) set to 169 1.0. Only the recursive results were used because they allowed for different gap thresholds , and 4 from the present study). One site (Frijolitos) was sampled during both studies, but 183 we analyzed them separately to preserve independence between the two studies. We then applied 184 the same rarefaction/extrapolation approach described above to generate rarefaction curves as 185 function of the number of individuals sampled. Although there might be important 186 methodological differences between the two studies, our objective here is to provide a general 9 187 overview molecular diversity across the sites, rather than a precise estimate within each site. 189 Results 190 We collected approximately 300 specimens across the four sites; however, our analysis focused 191 on the 224 individuals that were successfully barcoded (Table 1 ; S1 Table) . We were able to 192 identify nearly 70% of individuals to genus level using the morphological approach, but species-193 level identification was only possible for 56 individuals (25%). Some of the most numerous taxa 194 across sites included Leptophlebiidae (11.2 % of individuals), Libellulidae (11.2 %), Naucoridae 195 (6.7 %), Notonectidae (6.3 %), Chironomidae (5.4%), Gerridae (4.9 %), Hydropsychidae (4.0 196 %), Perlidae (3.6 %), and Baetidae (3.6 %). 197 The 224 COI sequences revealed a total of 168 haplotypes (S1 Table) . 225 Both ML and BI inference trees for all specimens showed well-defined clades at the level of 226 order and family, with some differences in the topology, but overall support was higher for the 227 BI tree, which we used to represent the number of molecular species (Fig 2; S1 Fig) . (Fig 4) , and this pattern was consistent across the first three Hill numbers 247 (Fig S2) . Similar results were found when looking at diversity Hill across the Isthmus of Panama 248 using the compiled barcoding data set. In particular, we observed substantial diversity of both 249 MOTUs and unique haplotypes across sites (Fig 4) , but the accumulation curves did not reach 250 saturation (Fig 4) . In addition, both MOTUs and haplotype diversity tended to increase at sites 251 within the PCW, in contrast to sites located at the eastern and western portion of the country (Fig 252 4) . Data on assignment and diversity of MOTUS across study sites are available in the 253 supplementary material (S1 Table) . (Table 1 ; Fig 2) , and a large portion of these lineages appear to be unique to 272 each site (Fig 3) . In addition, our rarefaction/extrapolation approach showed that this diversity is 273 still undersampled across sites both within the PCW and along the Isthmus of Panama (Fig 4) . These findings confirm that the diversity of freshwater macroinvertebrates in Neotropical 275 environments is largely understudied [10,34,68], and could be much higher than previously 276 thought. In addition, the fact that only a small number of specimens from these taxa matched 277 available sequences in public databases further highlights the potential for biodiversity discovery 278 in Neotropical freshwater environments. This seems particularly relevant for taxa such as 279 Hydropsychidae, Gerridae, Chironomidae, Leptophlebiidae, Libellulidae and Notonectidae, 280 which showed high lineage/haplotype diversity across sites (Table 1 ; Fig 2) . Some of these taxa 281 also showed high haplotype diversity in a previous molecular study across Panama [26] , and are 13 283 [6,10]. They also span several functional feeding groups such as filterers (Hydropsychidae), 284 scrapers (Chironomidae) and predators (Notonectidae), which play a crucial role in freshwater The role of macroinvertebrates in stream ecosystem function The Role of Benthic Invertebrate Species in Freshwater 335 Ecosystems: Zoobenthic species influence energy flows and nutrient cycling The Role 338 of Biodiversity in the Functioning of Freshwater and Marine Benthic Ecosystems Importance of regional 341 diversity and environmental conditions on local species richness of aquatic macro-342 invertebrates in tropical forested streams Guía para el estudio de los macroinvertebrados acuáticos del Departamento de 345 Fondo para la Protección del Medio Ambiente Odonata of Costa Rica: Diversity and checklist of 347 species Clave ilustrada y diversidad de las larvas de Trichoptera 350 en el Nor-Occidente de Colombia / Illustrated key for the larvae of Trichoptera in the 351 northwest of Colombia Clave taxonómica para larvas de las familias del orden Trichoptera (Insecta) de 354 Aquatic insect diversity of Costa Rica: state of knowledge 11. Flowers RW. Diversity of stream-living insects in northwestern Panama Insects of Panama and Mesoamerica: Selected studies Macroinvertebrates in neotropical streams: richness patterns along 366 a catchment and assemblage structure between 2 seasons Organization of macroinvertebrate communities at a hierarchy of 369 spatial scales in a tropical stream Spatial scale and the diversity of macroinvertebrates in a 372 Neotropical catchment An introduction to the aquatic insects of North America 4 -Macroinvertebrates: Composition Life Histories and Production Latitudinal variations in organic diversity Species Diversity in Space and Time Structure and diversity of stream invertebrate 384 assemblages: the influence of temperature with altitude and latitude Short Communication: Insect biodiversity in freshwater ecosystems: is there any 387 latitudinal gradient? New vistas in Neotropical stream ecology-389 Are tropical streams really different? Why are there so many insect species? Perspectives from fossils and 18 393 phylogenies Freshwater biodiversity and aquatic insect 396 diversification Long-term isolation and endemicity of 399 Neotropical aquatic insects limit the community responses to recent amphibian decline. 400 Diversity and Distributions The biogeography of lower Mesoamerican freshwater fishes Tempo and mode of 404 allopatric divergence in the weakly electric fish Sternopygus dariensis in the Accelerated 407 species inventory on Madagascar using coalescent-based models of species delineation Beta diversity at multiple 410 hierarchical levels: explaining the high diversity of scarab beetles in tropical montane 411 forests 416 Freshwater biodiversity: importance, threats, status and conservation challenges Global 419 patterns of stream detritivore distribution: implications for biodiversity loss in changing 420 climates Global 426 patterns of freshwater species diversity, threat and endemism Effects of multiple 429 stressors associated with agriculture on stream macroinvertebrate communities in a tropical 430 catchment Low predictability of colour 432 polymorphism in introduced Guppy (Poecilia reticulata) Populations in Panama. PLOS 433 ONE Tropical fish community does not 435 recover 45 years after predator introduction Extinction rates of North American freshwater fauna Exploring the effects of 440 salinization on trophic diversity in freshwater ecosystems: a quantitative review Guía para el estudio de los macroinvertebrados acuáticos del Departamento 443 de Antioquia An inexpensive, automation-445 friendly protocol for recovering high-quality DNA DNA primers for amplification of 448 mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates Ten species in one: DNA 452 barcoding reveals cryptic species in the neotropical skipper butterfly Astraptes fulgerator Convenient single-step, one tube 455 purification of PCR products for direct sequencing Geneious 21 458 Basic: An integrated and extendable desktop software platform for the organization and 459 analysis of sequence data MAFFT multiple sequence alignment software version 7: 462 improvements in performance and usability New Tools for Microbial Genomics, Phylogenetics, and Molecular 465 Evolution IQ-TREE: a fast and effective 467 stochastic algorithm for estimating maximum-likelihood phylogenies BEAST 2: a software 470 platform for Bayesian evolutionary analysis Creating the CIPRES Science Gateway for inference of 473 large phylogenetic trees Ultrafast approximation for phylogenetic 477 bootstrap New algorithms 22 479 and methods to estimate maximum-likelihood phylogenies: assessing the performance of 480 PhyML 3.0 LogCombiner v1. 8.2. 2013. 482 54 Bayesian phylogenetics with BEAUti and 484 the BEAST 1.7 Turning DNA Barcode 487 Sequences into Annotated Operational Taxonomic Units MEGA7: Molecular Evolutionary Genetics Analysis 490 Version 7.0 for Bigger Datasets Genetic 493 diversity of armored scales (Hemiptera: Diaspididae) and soft scales (Hemiptera: Coccidae) 494 in Chile 496 From water striders to water bugs: the molecular diversity of aquatic Heteroptera 497 (Gerromorpha, Nepomorpha) of Germany based on DNA barcodes Incorporation 23 500 of DNA barcoding into a large-scale biomonitoring program: opportunities and pitfalls Can DNA barcodes of stream 503 macroinvertebrates improve descriptions of community structure and water quality Cryptic 506 biodiversity in streams: a comparison of macroinvertebrate communities based on 507 morphological and DNA barcode identifications. Freshwater Science Delimiting species using single-locus data and the Generalized 510 Mixed Yule Coalescent approach: a revised method and evaluation on simulated data sets Automatic Barcode Gap Discovery 513 for primary species delimitation Models and 516 estimators linking individual-based and sample-based rarefaction, extrapolation and 517 comparison of assemblages Rarefaction and 519 extrapolation with Hill numbers: a framework for sampling and estimation in species 520 diversity studies iNEXT: an R package for rarefaction and extrapolation of 522 species diversity (Hill numbers) Structure of macroinvertebrate 525 communities in riffles of a Neotropical karst stream in the wet and dry seasons Cross-Cordillera exchange mediated by the Panama 528 Canal increased the species richness of local freshwater fish assemblages Biodiversity beyond trees: Panama's Canal provides limited 532 conservation lessons for Nicaragua Threats to the running water ecosystems of the world Effects of deforestation on macroinvertebrate diversity and 537 assemblage structure in Ecuadorian Amazon streams Hybridization and 540 contemporary evolution in an introduced cichlid fish from Lake Malawi National Park Rarefaction and extrapolation curves for molecular diversity (MOTUs) at each site Number at the top represent fit curves for the first three Hill numbers: species richness (q = 0), 546 the exponential of Shannon entropy Capira (CAP), Frijolitos (FRI) Frijolitos 550 (FTO), Frijoles (FES) and Indio (IND) For each specimen we show taxonomic group (i.e., order, family 558 and genus), molecular species identity (MOTUs: based on ABGD and GMYC)