Hierarchical Hybrid Enrichment: Multitiered Genomic Data Collection Across Evolutionary Scales, With Application to Chorus Frogs (Pseudacris).

Determining the optimal targets of genomic subsampling for phylogenomics, phylogeography, and population genomics remains a challenge for evolutionary biologists. Of the available methods for subsampling the genome, hybrid enrichment (sequence capture) has become one of the primary means of data collection for systematics, due to the flexibility and cost efficiency of this approach. Despite the utility of this method, information is lacking as to what genomic targets are most appropriate for addressing questions at different evolutionary scales. In this study, first, we compare the benefits of target loci developed for deep- and shallow scales by comparing these loci at each of three taxonomic levels: within a genus (phylogenetics), within a species (phylogeography), and within a hybrid zone (population genomics). Specifically, we target evolutionarily conserved loci that are appropriate for deeper phylogenetic scales and more rapidly evolving loci that are informative for phylogeographic and population genomic scales. Second, we assess the efficacy of targeting multiple-locus sets for different taxonomic levels in the same hybrid enrichment reaction, an approach we term hierarchical hybrid enrichment. Third, we apply this approach to the North American chorus frog genus Pseudacris to answer key evolutionary questions across taxonomic and temporal scales. We demonstrate that in this system the type of genomic target that produces the most resolved gene trees differs depending on the taxonomic level, although the potential for error is substantially lower for the deep-scale loci at all levels. We successfully recover data for the two different locus sets with high efficiency. Using hierarchical data targeting deep and shallow levels: we 1) resolve the phylogeny of the genus Pseudacris and introduce a novel visual and hypothesis testing method that uses nodal heat maps to examine the robustness of branch support values to the removal of sites and loci; 2) estimate the phylogeographic history of Pseudacris feriarum, which reveals up to five independent invasions leading to sympatry with congener Pseudacris nigrita to form replicated reinforcement contact zones with ongoing gene flow into sympatry; and 3) quantify with high confidence the frequency of hybridization in one of these zones between P. feriarum and P. nigrita, which is lower than microsatellite-based estimates. We find that the hierarchical hybrid enrichment approach offers an efficient, multitiered data collection method for simultaneously addressing questions spanning multiple evolutionary scales. [Anchored hybrid enrichment; heat map; hybridization; phylogenetics; phylogeography; population genomics; reinforcement; reproductive character displacement.].

New insights into the identities of the elasmobranch fauna of Sri Lanka.

Identities of elasmobranchs from Sri Lanka encountered during collections conducted in an intensive nine-day survey of fish markets and landing sites at 11 localities in the North Western, Northern, and Eastern Provinces in March of 2018 were assessed. In total, 111 specimens representing 34 elasmobranch species were examined. Sequence data for the NADH2 gene were generated for all specimens. Independent Neighbor-Joining analyses, which included data for related taxa, were conducted for 25 subgroups of elasmobranchs to help confirm specific identifications. Five of the 34 species encountered are likely new to science. These consist of one species each of the batoid genera Brevitrygon, Narcine, and Torpedo, and the selachian genera Centrophorus, and Chiloscyllium. The specific identities of 12 species previously known to occur in Sri Lanka are updated to conform to current taxonomy; four of these (Gymnura cf. poecilura 2, Carcharhinus cf. limbatus, Echinorhinus sp. 1, and Iago cf. omanensis 1) represent what appear to be undescribed species reported previously from other localities. Three species (Maculabatis arabica, Acroteriobatus variegatus, and Centroscymnus owstonii) are reported from Sri Lanka for the first time; the latter species also represents the first documented record of this genus and family for the island nation. One of the two specimens on which the recent description of the new species of Planonasus indicus was based was also collected as part of this survey. Although some of the species confirmed to occur in Sri Lanka have also been found in India, others were previously known only from the Persian Gulf, Gulf of Oman, or localities in Southeast Asia. The high amount of novelty discovered as a result of a survey of such short duration emphasizes the importance of more intensive survey efforts in this region now that the civil unrest that precluded such work for nearly three decades has come to an end.

Assessing the links among environmental contaminants, endocrinology, and parasites to understand amphibian declines in montane regions of Costa Rica.

Amphibians inhabiting montane riparian zones in the Neotropics are particularly vulnerable to decline, but the reasons are poorly understood. Because environmental contaminants, endocrine disruption, and pathogens often figure prominently in amphibian declines it is imperative that we understand how these factors are potentially interrelated to affect montane populations. One possibility is that increased precipitation associated with global warming promotes the deposition of contaminants in montane regions. Increased exposure to contaminants, in turn, potentially elicits chronic elevations in circulating stress hormones that could contribute to montane population declines by compromising resistance to pathogens and/or production of sex steroids regulating reproduction. Here, we test this hypothesis by examining contaminant levels, stress and sex steroid levels, and nematode abundances in male drab treefrogs, Smilisca sordida, from lowland and montane populations in Costa Rica. We found no evidence that montane populations were more likely to possess contaminants (i.e., organochlorine, organophosphate and carbamate pesticides or benzidine and chlorophenoxy herbicides) than lowland populations. We also found no evidence of elevational differences in circulating levels of the stress hormone corticosterone, estradiol or progesterone. However, montane populations possessed lower androgen levels, hosted more nematode species, and had higher nematode abundances than lowland populations. Although these results suggested that nematodes contributed to lower androgens in montane populations, we were unable to detect a significant inverse relationship between nematode abundance and androgen level. Our results suggest that montane populations of this species are not at greater risk of exposure to contaminants or chronic stress, but implicate nematodes and compromised sex steroid levels as potential threats to montane populations.

Species tree estimation of North American chorus frogs (Hylidae: Pseudacris) with parallel tagged amplicon sequencing.

The field of phylogenetics is changing rapidly with the application of high-throughput sequencing to non-model organisms. Cost-effective use of this technology for phylogenetic studies, which often include a relatively small portion of the genome but several taxa, requires strategies for genome partitioning and sequencing multiple individuals in parallel. In this study we estimated a multilocus phylogeny for the North American chorus frog genus Pseudacris using anonymous nuclear loci that were recently developed using a reduced representation library approach. We sequenced 27 nuclear loci and three mitochondrial loci for 44 individuals on 1/3 of an Illumina MiSeq run, obtaining 96.5% of the targeted amplicons at less than 20% of the cost of traditional Sanger sequencing. We found heterogeneity among gene trees, although four major clades (Trilling Frog, Fat Frog, crucifer, and West Coast) were consistently supported, and we resolved the relationships among these clades for the first time with strong support. We also found discordance between the mitochondrial and nuclear datasets that we attribute to mitochondrial introgression and a possible selective sweep. Bayesian concordance analysis in BUCKy and species tree analysis in (*)BEAST produced largely similar topologies, although we identify taxa that require additional investigation in order to clarify taxonomic and geographic range boundaries. Overall, we demonstrate the utility of a reduced representation library approach for marker development and parallel tagged sequencing on an Illumina MiSeq for phylogenetic studies of non-model organisms.