A trans-Amazonian screening of mtDNA reveals deep intraspecific divergence in forest birds and suggests a vast underestimation of species diversity

The Amazonian avifauna remains severely understudied relative to that of the temperate zone, and its species richness is thought to be underestimated by current taxonomy. Recent molecular systematic studies using mtDNA sequence reveal that traditionally accepted species-level taxa often conceal genetically divergent subspecific lineages found to represent new species upon close taxonomic scrutiny, suggesting that intraspecific mtDNA variation could be useful in species discovery. Surveys of mtDNA variation in Holarctic species have revealed patterns of variation that are largely congruent with species boundaries. However, little information exists on intraspecific divergence in most Amazonian species. Here we screen intraspecific mtDNA genetic variation in 41 Amazonian forest understory species belonging to 36 genera and 17 families in 6 orders, using 758 individual samples from Ecuador and French Guiana. For 13 of these species, we also analyzed trans-Andean populations from the Ecuadorian Chocó. A consistent pattern of deep intraspecific divergence among trans-Amazonian haplogroups was found for 33 of the 41 taxa, and genetic differentiation and genetic diversity among them was highly variable, suggesting a complex range of evolutionary histories. Mean sequence divergence within families was the same as that found in North American birds (13%), yet mean intraspecific divergence in Neotropical species was an order of magnitude larger (2.13% vs. 0.23%), with mean distance between intraspecific lineages reaching 3.56%. We found no clear relationship between genetic distances and differentiation in plumage color. Our results identify numerous genetically and phenotypically divergent lineages which may result in new species-level designations upon closer taxonomic scrutiny and thorough sampling, although lineages in the tropical region could be older than those in the temperate zone without necessarily representing separate species. In-depth phylogeographic surveys are urgently needed to avoid underestimating tropical diversity, and the use of mtDNA markers can be instrumental in identifying and prioritizing taxa for species discovery.     

Genetic variation and geographical structure of five mouse-sized opossums (Marsupialia, Didelphidae) throughout the Guiana Region

Aim and location To study the patterns of genetic variation and geographical structure of five forest‐dwelling didelphid species distributed throughout the Guiana Region (Marmosops parvidens, M. pinheiroi, Monodelphis brevicaudata, Marmosa murina, Micoureus demerarae) and place the results for wide‐ranging species (M. murina, M. demerarae) in the broader geographical context of the rest of the taxon ranges based on published data. Methods Variation in the mitochondrial cytochrome b gene were analysed in order to assess intra‐specific sequence divergence as well as haplotypic and nucleotide diversity among populations. Relationships between haplotypes were inferred by cladistic (maximum parsimony) and probabilistic (maximum likelihood) approaches, allowing comparisons between tree topologies and patterns of populations geographical distribution. Results Phylogenetic analysis of the Guiana Region populations suggest the lack of common patterns of geographical structure among mouse‐sized opossums except for both Marmosops species. Marmosops parvidens and M. pinheiroi showed individuals from Surinam strongly related with those from Guyana, and haplotypes from French Guiana having a basal position. Concerning Micoureus demerarae, haplotypes from Venezuela were positioned as ancestral, in contrast with Monodelphis brevicaudata, whose French Guiana haplotypes were clustered in a basal position. No evidence of geographical structure was observed for Marmosa murina. Genetic variation within Guiana Region populations fluctuated between 1% for M. murina and 7% in the case of M. brevicaudata. Regarding the widely distributed species, phylogeographical structure at the scale of the Amazon Basin suggests that Guiana Region populations are the sister group of south‐eastern Amazonian and Atlantic Forest populations, thus resulting in an eastern clade well separated from the western Amazonian regions. Main conclusions Monophyly of mouse‐sized opossums haplotypes over the Guiana Region confirms the presence of a well‐differentiated zoogeographical area relative to other Amazonia bioregions. Low genetic diversity among Guiana Region samples suggests a relatively recent origin of populations dating from later Miocene and Pleistocene periods. Moreover, nucleotide and haplotypic diversity values suggest an evolutionary scenario of rapid population growth and dispersal over the Guiana territory, from an ancestral population with small effective size. Dispersal events – rather than vicariance – seem to be responsible for the present genetic and phylogeographical patterns observed in the Guiana Region.     

Population genetics of the hispid cotton rat (Sigmodon hispidus): patterns of genetic diversity at the major histocompatibility complex

The hispid cotton rat, Sigmodon hispidus, is a common rodent widely distributed across the southern USA and south into South America. To characterize major histocompatibility complex (MHC) diversity in this species and to elucidate large-scale patterns of genetic partitioning, we examined MHC genetic variability within and among 13 localities, including a disjunct population in Arizona and a population from Costa Rica that may represent an undescribed species. We also tested the hypothesis that populations within the USA are at equilibrium with regard to gene flow and genetic drift, resulting in isolation-by-distance. Using single-strand conformation polymorphism (SSCP) analysis we identified 25 alleles from 246 individuals. Gene diversity within populations ranged from 0.000 to 0.908. Analysis of molecular variance (AMOVA) revealed that 83.7% of observed variation was accounted for by within-population diversity and 16.3% was accounted for by among-population divergence. The disjunct population in Arizona was fixed for a single allele. The Costa Rican population was quite divergent based on allelic composition and was the only population with unique alleles. Within the main portion of the geographical distribution of S. hispidus in the USA there was considerable divergence among some populations; however, there was no significant pattern of isolation-by-distance overall (P = 0.090). Based on the significant divergence of the only sampled population to its east, the Mississippi River appears to represent a substantial barrier to gene flow.     

Vector soup: high‐throughput identification of Neotropical phlebotomine sand flies using metabarcoding

Phlebotomine sand flies are haematophagous dipterans of primary medical importance. They represent the only proven vectors of leishmaniasis worldwide and are involved in the transmission of various other pathogens. Studying the ecology of sand flies is crucial to understand the epidemiology of leishmaniasis and further control this disease. A major limitation in this regard is that traditional morphological‐based methods for sand fly species identifications are time‐consuming and require taxonomic expertise. DNA metabarcoding holds great promise in overcoming this issue by allowing the identification of multiple species from a single bulk sample. Here, we assessed the reliability of a short insect metabarcode located in the mitochondrial 16S rRNA for the identification of Neotropical sand flies, and constructed a reference database for 40 species found in French Guiana. Then, we conducted a metabarcoding experiment on sand flies mixtures of known content and showed that the method allows an accurate identification of specimens in pools. Finally, we applied metabarcoding to field samples caught in a 1‐ha forest plot in French Guiana. Besides providing reliable molecular data for species‐level assignations of phlebotomine sand flies, our study proves the efficiency of metabarcoding based on the mitochondrial 16S rRNA for studying sand fly diversity from bulk samples. The application of this high‐throughput identification procedure to field samples can provide great opportunities for vector monitoring and eco‐epidemiological studies.     

Molecular systematics and phylogeography of Amazonian poison frogs of the genus Dendrobates

The study of Amazonian biodiversity requires detailed knowledge of the phylogenetic relationships of closely related taxa distributed across Amazonia. The Amazonian poison frogs of the genus Dendrobates have undergone many taxonomic revisions, but the phylogenetic relationships within this group remain poorly understood. Most previous classifications were based on morphology and skin toxin analyses, with limited use of DNA sequence data. Using mtDNA sequence data from four gene regions (cytochrome b, cytochrome oxidase I, 16S rRNA, and 12S rRNA), we present a molecular phylogenetic analysis of the evolutionary relationships within a representative group of Amazonian Dendrobates. We use the resulting phylogenetic hypothesis to investigate different biogeographic hypotheses concerning genetic divergence and species diversity in Amazonia. The results of the analysis support the presence of ancient paleogeographic barriers to gene flow between eastern and western Amazonia, and indicate substantial genetic divergence between species found in the northern and southern regions of western Amazonia.     

Variable phylogeographic histories of five forest birds with populations in Upper and Lower Guinea: implications for taxonomy and evolutionary conservation

This study evaluates phylogeographic structure in selected forest bird species distributed on either side of the Dahomey Gap in West Africa to extend analyses of avian diversification across the Guineo-Congolean forest. Mitochondrial DNA sequence data were collected for five species, including four with a break in distribution across the Dahomey Gap and one with a more continuous distribution. Data for the five focal species were included in a single analysis along with appropriate outgroup taxa, allowing estimates of divergence time on a common timescale. Both phylogeographic patterns and estimated divergence times were highly variable. In two cases (Green Hylia Hylia prasina and Yellow-whiskered Greenbul Eurillas latirostris), samples from east and west of the Dahomey Gap belong to divergent mitochondrial lineages, but with substantially older divergence in Eurillas (～4.1 MY). In Western Bluebill Spermophaga haematina, relatively divergent eastern and western lineages (～1.1 MY) were sampled in the same locality in eastern Ghana, suggesting east to west gene flow across the gap following earlier vicariance. In Forest Robin Stiphrornis erythrothorax, birds sampled within Ghana were from two highly divergent lineages (～1.9 MY), one of which is more closely related to birds from Cameroon and Equatorial Guinea, such that divergence across the Dahomey Gap was comparatively recent (～0.36 MY). For both Hylia and Stiphrornis, we corroborate recent studies documenting multiple divergent lineages within each genus, some of which may represent distinct species. Consistent with previous studies, our analysis shows that the Dahomey Gap can be a significant barrier to gene flow for species restricted to forest habitats. At the same time, individual species exhibit a broad range of divergence times and phylogeographic histories of variable complexity. We discuss the taxonomic implications of our results for each focal species.     

Phylogenetic analysis,taxonomic revision,and dental ontogeny of the Cretaceous Zhelestidae (Mammalia: Eutheria)

Iphisa elegans Gray, 1851 is a ground‐dwelling lizard widespread over Amazonia that displays a broadly conserved external morphology over its range. This wide geographical distribution and conservation of body form contrasts with the expected poor dispersal ability of the species, the tumultuous past of Amazonia, and the previously documented prevalence of cryptic species in widespread terrestrial organisms in this region. Here we investigate this homogeneity by examining hemipenial morphology and conducting phylogenetic analyses of mitochondrial (CYTB) and nuclear (C‐MOS) DNA sequence data from 49 individuals sampled across Amazonia. We detected remarkable variation in hemipenial morphology within this species, with multiple cases of sympatric occurrence of distinct hemipenial morphotypes. Phylogenetic analyses revealed highly divergent lineages corroborating the patterns suggested by the hemipenial morphotypes, including co‐occurrence of different lineages. The degrees of genetic and morphological distinctness, as well as instances of sympatry among mtDNA lineages/morphotypes without nuDNA allele sharing, suggest that I. elegans is a complex of cryptic species. An extensive and integrative taxonomic revision of the I. elegans complex throughout its wide geographical range is needed. © 2012 The Linnean Society of London, Zoological Journal of the Linnean Society, 2012, 166 , 361–376.     

Systematic review of Atractus schach (Serpentes: Dipsadidae) species complex from the Guiana Shield with description of three new species

The Guiana Shield harbours one of the best preserved and largest extents of tropical forest on Earth and an immense biodiversity. The herpetofauna of this region remains poorly known. The species-rich snake genus Atractus contains ～140 species, many with complicated taxonomic histories, including A. schach. Examination of specimens in museums and newly collected material from French Guiana has allowed the illustration of hemipenial morphology for the first time and an expanded diagnosis. Concatenated molecular phylogenetic (mitochondrial and nuclear genes) and phenotypic (morphometrics, external and hemipenial morphology) analyses confirm non-monophyly of the A. flammigerus group and indicate that A. schach is a species complex with three new species described here. The geographic distribution of A. schach sensu stricto is restricted to Guiana, Surinam, and French Guiana north of Tumucumaque massif. Populations tentatively assigned to A. schach from the east from French Guiana in the Roura lowlands to Almeirim, and from central Amazonia between the Negro and Trombetas rivers in Brazil are also recognized as new species. Our results suggest that populations from south of the Amazon River are not conspecific with those from the Guiana Shield.

http://www.zoobank.org/urn:lsid:zoobank.org:pub:A7AE40BC-4716-4302-B3BE-1F43600B0A72     

Variation in tropical bird survival across longitude and guilds: a case study from the Amazon

Survival estimates of tropical birds have been used to examine life‐history variation across latitudes and dietary guilds. Here, we provide apparent survival estimates of 31 rainforest species from central Amazonian Brazil and compare our results with estimates from Ecuador, Peru (western Amazonia) and French Guiana (eastern Amazonia) to examine geographic variation in demography. Our averaged annual survival estimate for central Amazonian species (?= 0.59, SE = 0.10) was concordant with previously published estimates from the western Amazon (?= 0.58, SE = 0.02), and eastern Amazon (?= 0.63, SE = 0.06). Two predominate demographic patterns were detected across the study areas: within species, survival was highest or lowest in eastern or western Amazonia, but rarely in central Amazonia. The most striking demographic variation was exhibited by Pithys albifrons, for which apparent survival estimates were nearly twice as high in eastern Amazonia (?= 0.80, SE = 0.06) than in western Amazonia (?= 0.42, SE = 0.06) but intermediate in central Amazonia (?= 0.54, SE = 0.04). Although variation in survival may be associated with differences in life history characteristics, our analysis of flocking guild, body size, and nest architecture revealed only moderate differences in survival associated with nest architecture. These results suggest that geographic variation in population processes may be significant for widespread Amazonian species.     

Reinforcing and expanding the predictions of the disturbance vicariance hypothesis in Amazonian harlequin frogs: a molecular phylogenetic and climate envelope modelling approach

The disturbance vicariance hypothesis (DV) has been proposed to explain speciation in Amazonia, especially its edge regions, e.g. in eastern Guiana Shield harlequin frogs (Atelopus) which are suggested to have derived from a cool-adapted Andean ancestor. In concordance with DV predictions we studied that (i) these amphibians display a natural distribution gap in central Amazonia; (ii) east of this gap they constitute a monophyletic lineage which is nested within a pre-Andean/western clade; (iii) climate envelopes of Atelopus west and east of the distribution gap show some macroclimatic divergence due to a regional climate envelope shift; (iv) geographic distributions of climate envelopes of western and eastern Atelopus range into central Amazonia but with limited spatial overlap. We tested if presence and apparent absence data points of Atelopus were homogenously distributed with Ripley’s K function. A molecular phylogeny (mitochondrial 16S rRNA gene) was reconstructed using Maximum Likelihood and Bayesian Inference to study if Guianan Atelopus constitute a clade nested within a larger genus phylogeny. We focused on climate envelope divergence and geographic distribution by computing climatic envelope models with MaxEnt based on macroscale bioclimatic parameters and testing them by using Schoener’s index and modified Hellinger distance. We corroborated existing DV predictions and, for the first time, formulated new DV predictions aiming on species’ climate envelope change. Our results suggest that cool-adapted Andean Atelopus ancestors had dispersed into the Amazon basin and further onto the eastern Guiana Shield where, under warm conditions, they were forced to change climate envelopes.     

Phylogenetic relationships and biodiversity in Hylids (Anura: Hylidae) from French Guiana

We evaluated two biodiversity criteria, higher taxonomic diversity and phylogenetic diversity in French Guiana. For this, we used a recent assessment of the knowledge accumulated since 30 years of study on the amphibian species currently known in French Guiana. We focused on two well-represented genera, Hyla and Scinax, belonging to the subfamily Hylinae. We used partial sequences of two mitochondrial genes (16S rDNA and 12S rDNA, 813 bp) and two nuclear genes (tyrosinase and 18S rRNA, 1590 bp) covering a total of 2403 bp. According to the high bootstrap support in phylogenetic analysis of the complete dataset, the genus Scinax is a homophyletic clade formed by two species groups (rubra and rostrata) in French Guiana. The genus Hyla was confirmed to be a paraphyletic group formed by two species groups as well (30 chromosomes and the 'gladiator frogs'). We confirmed that these genera should be taxonomically reconsidered. Moreover, at the genus, subfamily and family levels, the use of only morphological characters or only molecular DNA markers would hamper estimations of biodiversity. Thus, we strongly advise the combined use of both morphology and molecular data (nuclear and mitochondrial markers).     

Multiple quaternary refugia in the eastern Guiana shield revealed by comparative phylogeography of 12 frog species

The Guiana Shield (GS) is one of the most pristine regions of Amazonia and biologically one of the richest areas on Earth. How and when this massive diversity arose remains the subject of considerable debate. The prevailing hypothesis of Quaternary glacial refugia suggests that a part of the eastern GS, among other areas in Amazonia, served as stable forested refugia during periods of aridity. However, the recently proposed disturbance-vicariance hypothesis proposes that fluctuations in temperature on orbital timescales, with some associated aridity, have driven Neotropical diversification. The expectations of the temporal and spatial organization of biodiversity differ between these two hypotheses. Here, we compare the genetic structure of 12 leaf-litter inhabiting frog species from the GS lowlands using a combination of mitochondrial and nuclear sequences in an integrative analytical approach that includes phylogenetic reconstructions, molecular dating, and Geographic Information System methods. This comparative and integrated approach overcomes the well-known limitations of phylogeographic inference based on single species and single loci. All of the focal species exhibit distinct phylogeographic patterns highlighting taxon-specific historical distributions, ecological tolerances to climatic disturbance, and dispersal abilities. Nevertheless, all but one species exhibit a history of fragmentation/isolation within the eastern GS during the Quaternary with spatial and temporal concordance among species. The signature of isolation in northern French Guiana (FG) during the early Pleistocene is particularly clear. Approximate Bayesian Computation supports the synchrony of the divergence between northern FG and other GS lineages. Substructure observed throughout the GS suggests further Quaternary fragmentation and a role for rivers. Our findings support fragmentation of moist tropical forest in the eastern GS during this period when the refuge hypothesis would have the region serving as a contiguous wet-forest refuge.     

Black flies (Diptera: Simuliidae) of French Guiana: cytotaxonomy and a preliminary list of species

The objectives of the present study were to broaden the survey of simuliid species in French Guiana and to cytologically analyze the species in the Simulium perflavum species group. Twelve species of Simulium were collected from which S. goeldii, S. quadrifidum, S. trombetense, S. near incrustatum, S. metallicum (s.l.) sp1, S. metallicum (s.l.) sp2 and S. ochraceum (s.l.) are reported for the first time for this region. The only species collected in the S. perflavum group was S. rorotaense; 34 larvae of this species were cytologically analyzed, all of which had the standard sequence. S. metallicum (s.l.), S. ochraceum (s.l.), S. guianense (s.l.) and S. oyapockense (s.l.) are involved with transmission of onchocerciasis in Central and South America, however, in French Guiana these species were not found biting humans during the sampling period. With the few collections made during this study, we increase the number of simulid species known in French Guiana from 6 to 13. It is clear that more simuliid species can be expected to be found when more sampling is done, including collections in other ecoregions in French Guiana.     

Two deep evolutionary lineages in the circumtropical glasseye Heteropriacanthus cruentatus (Teleostei,Priacanthidae) with admixture in the south‐western Indian Ocean

A phylogeographic study of the circumtropical glasseye Heteropriacanthus cruentatus was conducted. Molecular analyses indicate two mitochondrial cytochrome c oxidase subunit I (coI) lineages that are 10·4% divergent: one in the western Atlantic (Caribbean) and another that was detected across the Indo‐Pacific. A fixed single nucleotide polymorphism (SNP) was detected at a nuclear locus (S7 ribosomal protein) and is consistent with this finding. There is evidence of recent dispersal from the Atlantic to the Indian Ocean with individuals of mixed lineages detected in South Africa and the Mozambique Channel. Using coalescent analyses of the mitochondrial dataset, time of divergence between lineages was estimated to be c. 15·3 million years. The deep divergence between these two lineages indicates distinct evolutionary units, however, due to the lack of morphological differences and evidence of hybridization between lineages, taxonomic revision is not suggested at this time.     

Phylogeography of the chestnut‐tailed antbird (Myrmeciza hemimelaena) clarifies the role of rivers in Amazonian biogeography

Aim We examined patterns of spatial and temporal diversification of the Amazonian endemic chestnut‐tailed antbird, Mymeciza hemimelaena (Thamnophilidae), to evaluate the diversification of a widespread avian taxon across rivers that potentially represent major natural barriers. Location Lowland Amazonia. Methods Sequences of the mitochondrial ND2 and cytochrome b genes were investigated from 65 individuals distributed throughout the entire range of M. hemimelaena, and including the two currently valid subspecies M. h. hemimelaena and M. h. pallens. Based on a combination of phylogeographic tools, molecular dating, and population genetic methods, we reconstructed a spatio‐temporal scenario of diversification of M. hemimelaena in the Amazon. Results The data revealed three genetically divergent and monophyletic groups in M. hemimelaena, which can also be distinguished by a combination of morphological and vocal characters. Two of these clades correspond to the previously described taxa M. h. hemimelaena and M. h. pallens, which are separated by the upper Madeira River, a main Amazonian tributary. The third clade is distributed between the middle reaches of the Madeira River and the much smaller tributaries Jiparaná and Aripuanã, and, although currently treated as M. h. pallens, clearly constitutes an independent evolutionary lineage probably deserving separate species status. Molecular clock and population genetic analyses indicate that diversification in this group occurred throughout the Pleistocene, with demographic fluctuations assumed for M. h. hemimelaena and M. h. pallens. Main conclusions The findings implicate rivers as barriers driving diversification in the M. hemimelaena complex. Levels of mitochondrial DNA divergence and associated morphological and vocal traits support its division into at least three separate species with comparatively small ranges. The existence of a previously unrecognized lineage in the M. hemimelaena complex, and the high degree of population structuring found in M. h. hemimelaena underscore the pervasiveness of cryptic endemism throughout Amazonia and the importance of DNA‐based taxonomic and phylogeographic studies in providing the accurate estimates of diversity that are essential for conservation planning.     

Molecular phylogeny of French Guiana Hylinae: implications for the systematic and biodiversity of the Neotropical frogs

In this study we used nucleotide sequences from a segment of mitochondrial 16S ribosomal DNA gene to investigate the evolutionary relationships of some French Guiana Hylinae. New sequences, representing the members of different French Guiana frogs-five specimens of the Scinax genus, two Hyla, one Osteocephalus, one Hyalinobatrachium and two Rana as out-group-were examined. In addition, 26 sequences available from GenBank database representing the other subfamilies of the Hylidae were added to our study. This work allowed us to clarify relationships within the four hylids subfamilies (Pelodryadinae, Phyllomedusinae, Hemiphractinae and Hylinae) and the phylogenetic placement of the enigmatic Scinax genus within the Hylidae. We found that: (1) the Scinax genus displays a high level of differentiation in comparison to two other genera (Litoria and Hyla) belonging to 'Hylidae' family; (2) the Hylinae are paraphyletic given the position of the Litoria, which was the sister-group of the Hyla and the Osteocephalus genera; (3) the anterior works and our results (based on two different data sets) showed the paraphyly of the Hylidae questioning the validity of this family; (4) the reassessment of these different taxonomic groups will induce a huge implication on the estimation (past, present and future) of the biodiversity (in Neotropical frogs).     

Spatial expansion and population structure of the neotropical malaria vector, Anopheles darlingi (Diptera: Culicidae)

Extensive population structuring is known to occur in Anopheles darlingi , the primary malaria vector of the Neotropics. We analysed the phylogeographic structure of the species using the mitochondrial cytochrome oxidase I marker. Diversity is divided into six main population groups in South America: Colombia, central Amazonia, southern Brazil, south-eastern Brazil, and two groups in north-east Brazil. The ancestral distribution of the taxon is hypothesized to be central Amazonia, and there is evidence of expansion from this region during the late Pleistocene. The expansion was not a homogeneous front, however, with at least four subgroups being formed due to geographic barriers. As the species spread, populations became isolated from each other by the Amazon River and the coastal mountain ranges of south-eastern Brazil and the Andes. Analyses incorporating distances around these barriers suggest that the entire South American range of An. darlingi is at mutation–dispersal–drift equilibrium. Because the species is distributed throughout such a broad area, the limited dispersal across some landscape types promotes differentiation between otherwise proximate populations. Moreover, samples from the An. darlingi holotype location in Rio de Janeiro State are substantially derived from all other populations, implying that there may be additional genetic differences of epidemiological relevance. The results obtained contribute to our understanding of gene flow in this species and allow the formulation of human mosquito health protocols in light of the potential population differences in vector capacity or tolerance to control strategies.  © 2009 The Linnean Society of London, Biological Journal of the Linnean Society , 2009, 97 , 854–866.     

High-throughput identification of informative nuclear Loci for shallow-scale phylogenetics and phylogeography

One of the major challenges for researchers studying phylogeography and shallow-scale phylogenetics is the identification of highly variable and informative nuclear loci for the question of interest. Previous approaches to locus identification have generally required extensive testing of anonymous nuclear loci developed from genomic libraries of the target taxon, testing of loci of unknown utility from other systems, or identification of loci from the nearest model organism with genomic resources. Here, we present a fast and economical approach to generating thousands of variable, single-copy nuclear loci for any system using next-generation sequencing. We performed Illumina paired-end sequencing of three reduced-representation libraries (RRLs) in chorus frogs (Pseudacris) to identify orthologous, single-copy loci across libraries and to estimate sequence divergence at multiple taxonomic levels. We also conducted PCR testing of these loci across the genus Pseudacris and outgroups to determine whether loci developed for phylogeography can be extended to deeper phylogenetic levels. Prior to sequencing, we conducted in silico digestion of the most closely related reference genome (Xenopus tropicalis) to generate expectations for the number of loci and degree of coverage for a particular experimental design. Using the RRL approach, we: (i) identified more than 100,000 single-copy nuclear loci, 6339 of which were obtained for divergent conspecifics and 904 of which were obtained for heterospecifics; (ii) estimated average nuclear sequence divergence at 0.1% between alleles within an individual, 1.1% between conspecific individuals that represent two different clades, and 1.8% between species; and (iii) determined from PCR testing that 53% of the loci successfully amplify within-species and also many amplify to the genus-level and deeper in the phylogeny (16%). Our study effectively identified nuclear loci present in the genome that have levels of sequence divergence on par with mitochondrial loci commonly used in phylogeography. Specifically, we estimated that ～7% of loci in the chorus frog genome are >3% divergent within species; this translates to a prediction of approximately 50,000 single-copy loci in the genome with >3% divergence. Moreover, successful amplification of many loci at deeper phylogenetic levels indicates that the RRL approach represents an efficient method for rapid identification of informative loci for both phylogenetics and phylogeography. We conclude by making recommendations for minimizing the cost and maximizing the efficiency of locus identification for future studies in this field.     

Phylogenetic relationships among Southern Balkan Rutilus species inferred from cytochrome b sequence analysis: Micro-geographic resolution and taxonomic implications

The phylogenetic relationships and geographic distribution of Greek roaches (Rutilus spp.) were investigated by analyzing the complete mitochondrial cytochrome b gene sequence of 84 specimens collected from 15 southern Balkan lakes. Phylogenetic analysis revealed the existence of five highly divergent haplotype groups with mean pairwise sequence divergence between them ranging from 4.1 to 9%, namely at the level of values reported for different species. These five groups correspond to four species namely Rutilus rutilus, Rutilus prespensis, Rutilus panosi and Rutilus ylikiensis. On the other hand, the existence of two highly divergent haplotype groups, which are currently attributed to R. rutilus, suggests the re-examination of the current taxonomic status of this species.