Dense geographic and genomic sampling reveals paraphyly and a cryptic lineage in a classic sibling species complex

Incomplete or geographically biased sampling poses significant problems for research in phylogeography, population genetics, phylogenetics, and species delimitation. Despite the power of using genome-wide genetic markers in systematics and related fields, approaches such as the multispecies coalescent remain unable to easily account for unsampled lineages. The Empidonax difficilis/Empidonax occidentalis complex of small tyrannid flycatchers (Aves: Tyrannidae) is a classic example of widely distributed species with limited phenotypic geographic variation that was broken into two largely cryptic (or "sibling") lineages following extensive study. Though the group is well-characterized north of the US Mexico border, the evolutionary distinctiveness and phylogenetic relationships of southern populations remain obscure. In this article, we use dense genomic and geographic sampling across the majority of the range of the E. difficilis/E. occidentalis complex to assess whether current taxonomy and species limits reflect underlying evolutionary patterns, or whether they are an artifact of historically biased or incomplete sampling. We find that additional samples from Mexico render the widely recognized species-level lineage E. occidentalis paraphyletic, though it retains support in the best-fit species delimitation model from clustering analyses. We further identify a highly divergent unrecognized lineage in a previously unsampled portion of the group's range, which a cline analysis suggests is more reproductively isolated than the currently recognized species E. difficilis and E. occidentalis. Our phylogeny supports a southern origin of these taxa. Our results highlight the pervasive impacts of biased geographic sampling, even in well-studied vertebrate groups like birds, and illustrate what is a common problem when attempting to define species in the face of recent divergence and reticulate evolution

The biogeographic and evolutionary history of an endemic clade of Middle American sparrows: Melozone and Aimophila (Aves: Passerellidae)

The large number of endemic species in Middle America is frequently attributed to the interplay of geographical barriers and historical climatic changes in the region. This process promotes genetic divergence between populations, and given enough time, may yield new species. Animals that inhabit mid-elevation or highland habitats may be disproportionately affected in this way. Genetic analyses of animals in this region allow us to better understand how historical patterns of isolation have influenced the generation of new species in this biodiversity hotspot. We studied the biogeography and systematics of two closely related genera of sparrows (Passerellidae): Melozone and Aimophila. Collectively, this group is distributed from the southwestern United States and southward as far as central Costa Rica. We sampled 81 individuals of 8 Melozone and 2 Aimophila species, from 19 localities distributed throughout their ranges. We reconstructed phylogenetic relationships and time-calibrated species trees using multilocus sequence data comprised of one mitochondrial gene and five nuclear genes. We conducted an ancestral area reconstruction analysis to determine the probability of ancestral range at each divergent event. Despite analyzing six loci, we were unable to obtain a fully resolved phylogenetic tree. We recovered four main lineages: lineage 1 includes four Melozone species distributed north of Isthmus of Tehuantepec (M. albicollis, M. crissalis, M. aberti, M. fusca); lineage 2 includes three Melozone species distributed south of the Isthmus of Tehuantepec (M. biarcuata, M. cabanisi, M. leucotis); lineage 3 lineage consists of a single species endemic to the Pacific coast of Mexico (M. kieneri); and lineage 4 includes the more widely distributed sparrows in the genus Aimophila. Our analyses suggest that these genera probably originated during the late Miocene in the Madrean Highlands of southern Mexico. We identified dispersal as the prevalent cause of speciation in this clade with most lineages dispersing to their current distributions from southern Mexico either to the north following a developing and expanding Madro-Tertiary flora, or to the south across the Isthmus of Tehuantepec. A similar pattern of dispersal from this biogeographic region has been reported in other taxa including fishes, reptiles, and birds. Our results reveal that the four lineages identified represent geographically coherent and ecologically similar assemblages of taxa. Finally, when our genetic results are considered, along with apparent differences in morphology and song, the allopatric forms M. b. cabanisi and M. l. occipitalis warrant recognition as biological species. © 201

Pciris_reads_4of6

100bp SR sequence reads from an illumina Hiseq 2500 for Passerina ciris. Library prepped via ddRAD protocol with enzymes sbf1 and msp1. Reads were demultiplexed in pyrad v 3-0-65 (command "pyrad -p params.txt -s 1") but have had no other processing

Pciris_reads_2of6

100bp SR sequence reads from an illumina Hiseq 2500 for Passerina ciris. Library prepped via ddRAD protocol with enzymes sbf1 and msp1. Reads were demultiplexed in pyrad v 3-0-65 (command "pyrad -p params.txt -s 1") but have had no other processing

Pciris_reads_6of6

100bp SR sequence reads from an illumina Hiseq 2500 for Passerina ciris. Library prepped via ddRAD protocol with enzymes sbf1 and msp1. Reads were demultiplexed in pyrad v 3-0-65 (command "pyrad -p params.txt -s 1") but have had no other processing

uce_loci_50p

uce_loci_50

ND2_alignment

Nexus file for NADH dehydrogenase subunit 2 (ND2) sequence alignmen

Pciris_reads_3of6

100bp SR sequence reads from an illumina Hiseq 2500 for Passerina ciris. Library prepped via ddRAD protocol with enzymes sbf1 and msp1. Reads were demultiplexed in pyrad v 3-0-65 (command "pyrad -p params.txt -s 1") but have had no other processing

S2

Table S2: HZAR confidence interval

Data from: Dense geographic and genomic sampling reveals paraphyly and a cryptic lineage in a classic sibling species complex

Incomplete or geographically biased sampling poses significant problems for research in phylogeography, population genetics, phylogenetics, and species delimitation. Despite the power of using genome-wide genetic markers in systematics and related fields, approaches such as the multispecies coalescent remain unable to easily account for unsampled lineages. The Empidonax difficilis / E. occidentalis complex of small tyrannid flycatchers (Aves: Tyrannidae) is a classic example of widely-distributed species with limited phenotypic geographic variation that was broken into two largely cryptic (or “sibling”) lineages following extensive study. Though the group is well-characterized north of the U.S. Mexico border, the evolutionary distinctiveness and phylogenetic relationships of southern populations remain obscure. In this paper, we use dense genomic and geographic sampling across the majority of the range of the E. difficilis / E . occidentalis complex to assess whether current taxonomy and species limits reflect underlying evolutionary patterns, or whether they are an artifact of historically biased or incomplete sampling. We find that additional samples from Mexico render the widely recognized species-level lineage E. occidentalis paraphyletic, though it retains support in the best-fit species delimitation model from clustering analyses. We further identify a highly divergent unrecognized lineage in a previously unsampled portion of the group’s range, which a cline analysis suggests is more reproductively isolated than the currently recognized species E. difficilis and E. occidentalis. Our phylogeny supports a southern origin of these taxa. Our results highlight the pervasive impacts of biased geographic sampling, even in well-studied vertebrate groups like birds, and illustrate what is a common problem when attempting to define species in the face of recent divergence and reticulate evolution