Sympatric and independently evolving lineages in the Thoropa miliaris - T. taophora species complex (Anura: Cycloramphidae)

Species delimitation can be challenging and affected by subjectivity. Sibling lineages that occur in sympatry constitute good candidates for species delimitation regardless of the adopted species concept. The Thoropa miliaris + T. taophora species complex exhibits high genetic diversity distributed in several lineages that occur sympatrically in the southeastern Atlantic Forest of Brazil. We used 414 loci obtained by anchored hybrid enrichment to characterize genetic variation in the Thoropa miliaris species group (T. saxatilis, T megatympanum, T. miliaris, and T. taophora), combining assignment analyses with traditional and coalescent phylogeny reconstruction. We also investigated evolutionary independence in co-occurring lineages by estimating gene flow, and validated lineages under the multispecies coalescent. We recovered most previously described lineages as unique populations in assignment analyses; exceptions include two lineages within T. miliaris that are further substructured, and the merging of all T. taophora lineages. We found very low probabilities of gene flow between sympatric lineages, suggesting independent evolution. Species tree inferences and species delimitation yielded resolved relationships and indicate that all lineages constitute putative species that diverged during the Pliocene and Pleistocene, later than previously estimated.

Phylogenomics of African radiation of Praomyini (Muridae: Murinae) rodents: First fully resolved phylogeny, evolutionary history and delimitation of extant genera

The tribe Praomyini is a diversified group including 64 species and eight extant rodent genera. They live in a broad spectrum of habitats across whole sub-Saharan Africa. Members of this tribe are often very abundant, they have a key ecological role in ecosystems, they are hosts of many potentially pathogenic microorganisms and comprise numerous agricultural pests. Although this tribe is well supported by both molecular and morphological data, its intergeneric relationships and the species contents of several genera are not yet fully resolved. Recent molecular data suggest that at least three genera in current sense are paraphyletic. However, in these studies the species sampling was sparse and the resolution of relationships among genera was poor, probably due to a fast radiation of the tribe dated to the Miocene and insufficient amount of genetic data. Here we used genomic scale data (395 nuclear loci = 610,965 bp long alignment and mitogenomes = 14,745 bp) and produced the first fully resolved species tree containing most major lineages of the Praomyini tribe (i.e. all but one currently delimited genera and major intrageneric clades). Results of a fossil-based divergence dating analysis suggest that the radiation started during the Messinian stage (ca. 7 Ma) and was likely linked to a fragmentation of the pan-African Miocene forest. Some lineages remained in the rain forests, while many others adapted to a broad spectrum of new open lowland and montane habitats that appeared at the beginning of Pliocene. Our analyses clearly confirmed the presence of three polyphyletic genera (Praomys, Myomyscus and Mastomys). We review current knowledge of these three genera and suggest corresponding taxonomic changes. To keep genera monophyletic, we propose taxonomic re-arrangements and delimit four new genera. Furthermore, we discovered a new highly divergent genetic lineage of Praomyini in southwestern Ethiopia, which is described as a new species and genus.

Gene flow in phylogenomics: Sequence capture resolves species limits and biogeography of Afromontane forest endemic frogs from the Cameroon Highlands

Puddle frogs of the Phrynobatrachus steindachneri species complex are a useful group for investigating speciation and phylogeography in Afromontane forests of the Cameroon Volcanic Line, western Central Africa. The species complex is represented by six morphologically relatively cryptic mitochondrial DNA lineages, only two of which are distinguished at the species level - southern P. jimzimkusi and Lake Oku endemic P. njiomock, leaving the remaining four lineages identified as 'P. steindachneri'. In this study, the six mtDNA lineages are subjected to genomic sequence capture analyses and morphological examination to delimit species and to study biogeography. The nuclear DNA data (387 loci; 571,936 aligned base pairs) distinguished all six mtDNA lineages, but the topological pattern and divergence depths supported only four main clades: P. jimzimkusi, P. njiomock, and only two divergent evolutionary lineages within the four 'P. steindachneri' mtDNA lineages. One of the two lineages is herein described as a new species, P. amieti sp. nov. Reticulate evolution (hybridization) was detected within the species complex with morphologically intermediate hybrid individuals placed between the parental species in phylogenomic analyses, forming a ladder-like phylogenetic pattern. The presence of hybrids is undesirable in standard phylogenetic analyses but is essential and beneficial in the network multispecies coalescent. This latter approach provided insight into the reticulate evolutionary history of these endemic frogs. Introgressions likely occurred during the Middle and Late Pleistocene climatic oscillations, due to the cyclic connections (likely dominating during cold glacials) and separations (during warm interglacials) of montane forests. The genomic phylogeographic pattern supports the separation of the southern (Mt. Manengouba to Mt. Oku) and northern mountains at the onset of the Pleistocene. Further subdivisions occurred in the Early Pleistocene, separating populations from the northernmost (Tchabal Mbabo, Gotel Mts.) and middle mountains (Mt. Mbam, Mt. Oku, Mambilla Plateau), as well as the microendemic lineage restricted to Lake Oku (Mt. Oku). This unique model system is highly threatened as all the species within the complex have exhibited severe population declines in the past decade, placing them on the brink of extinction. In addition, Mount Oku is identified to be of particular conservation importance because it harbors three species of this complex. We, therefore, urge for conservation actions in the Cameroon Highlands to preserve their diversity before it is too late.

Tempo and mode of evolution of oryzomyine rodents (Rodentia, Cricetidae, Sigmodontinae): A phylogenomic approach

The tribe Oryzomyini is an impressive group of rodents, comprising 30 extant genera and an estimated 147 species. Recent remarkable advances in the understanding of the diversity, taxonomy and systematics of the tribe have mostly derived from analyses of single or few genetic markers. However, the evolutionary history and biogeography of Oryzomyini, its origin and diversification across the Neotropics, remain unrevealed. Here we use a multi-locus dataset (over 400 loci) obtained through anchored phylogenomics to provide a genome-wide phylogenetic hypothesis for Oryzomyini and to investigate the tempo and mode of its evolution. Species tree and supermatrix analyses produced topologies with strong support for most branches, with all genera confirmed as monophyletic, a result that previous studies failed to obtain. Our analyses also corroborated the monophyly and phylogenetic relationship of three main clades of Oryzomyini (B, C and D). The origin of the tribe is estimated to be in the Miocene (8.93-5.38 million years ago). The cladogenetic events leading to the four main clades occurred during the late Miocene and early Pliocene and most speciation events in the Pleistocene. Geographic range estimates suggested an east of Andes origin for the ancestor of oryzomyines, most likely in the Boreal Brazilian region, which includes the north bank of Rio Amazonas and the Guiana Shield. Oryzomyini rodents are an autochthonous South America radiation, that colonized areas and dominions of this continent mainly by dispersal events. The evolutionary history of the tribe is deeply associated with the Andean cordillera and the landscape history of Amazon basin.

Nuclear phylogenomics, but not mitogenomics, resolves the most successful Late Miocene radiation of African mammals (Rodentia: Muridae: Arvicanthini)

The tribe Arvicanthini (Muridae: Murinae) is a highly diversified group of rodents (ca. 100 species) and with 18 African genera (plus one Asiatic) represents probably the most successful adaptive radiation of extant mammals in Africa. They colonized a broad spectrum of habitats (from rainforests to semi-deserts) in whole sub-Saharan Africa and their members often belong to most abundant parts of mammal communities. Despite intensive efforts, the phylogenetic relationships among major lineages (i.e. genera) remained obscured, which was likely caused by the intensive radiation of the group, dated to the Late Miocene. Here we used genomic scale data (377 nuclear loci; 581,030 bp) and produced the first fully resolved species tree containing all currently delimited genera of the tribe. Mitogenomes were also extracted, and while the results were largely congruent, there was less resolution at basal nodes of the mitochondrial phylogeny. Results of a fossil-based divergence dating analysis suggest that the African radiation started early after the colonization of Africa by a single arvicanthine ancestor from Asia during the Messinian stage (ca. 7 Ma), and was likely linked with a fragmentation of the pan-African Miocene forest. Some lineages remained in the rain forest, while many others successfully colonized broad spectrum of new open habitats (e.g. savannas, wetlands or montane moorlands) that appeared at the beginning of Pliocene. One lineage even evolved partially arboricolous life style in savanna woodlands, which allowed them to re-colonize equatorial forests. We also discuss delimitation of genera in Arvicanthini and propose corresponding taxonomic changes.

Anchored hybrid enrichment phylogenomics resolves the backbone of erebine moths

The subfamily Erebinae (Lepidoptera, Erebidae) includes approximately 10,000 species with many still undescribed. It is one of the most diverse clades within the moth superfamily Noctuoidea and encompasses a diversity of ecological habits. Erebine caterpillars feed on a broad range of host plants including several economically important crops. Adults possess a unique array of adaptations for predator defense, including some of the most sensitive hearing organs (tympana) across the Lepidoptera and striking wing coloration to startle visual predators. Despite the relevance of these moths to agriculture and ecological research, a robust phylogenetic framework is lacking. Here we used anchored hybrid enrichment, a relatively new approach in phylogenomics, to resolve relationships among the subfamily. Using the recently developed Lep1 anchored hybrid enrichment probe set, 658 gene fragments with an average length of 320 bp were captured from an exemplar set of 75 erebine species, representing 73 genera and 23 tribes. While the total number of erebine tribes is not firmly established, this represents at least 75% of known tribal level diversity. Anchored hybrid enrichment data were partitioned by locus and by codon position for maximum likelihood phylogenetic analysis and coalescent-based species-tree approaches. Results from our study provided strong nodal support (BP ≥ 95) for nearly all nodes in the partitioned ML tree, solidifying many relationships that were previously uncertain or moderately supported based on morphology or a smaller number of gene fragments. Likelihood analyses confidently resolved the placement of Acantholipini as a sister tribe to Sypnini and all other Erebinae. The remaining tribes were placed in a single, strongly supported clade split into two major subclades. Additionally, 25 tropical species that did not have previous tribal assignments are confidently placed on the phylogeny. Statistical comparisons with Shimodaira-Hasegawa (SH) tests found that our maximum likelihood trees were significantly more likely than alternative hypotheses. This study demonstrates the utility of anchored phylogenomics for resolving relationships within subfamilies of Lepidoptera.

New insights on the sister lineage of percomorph fishes with an anchored hybrid enrichment dataset

Percomorph fishes represent over 17,100 species, including several model organisms and species of economic importance. Despite continuous advances in the resolution of the percomorph Tree of Life, resolution of the sister lineage to Percomorpha remains inconsistent but restricted to a small number of candidate lineages. Here we use an anchored hybrid enrichment (AHE) dataset of 132 loci with over 99,000 base pairs to identify the sister lineage of percomorph fishes. Initial analyses of this dataset failed to recover a strongly supported sister clade to Percomorpha, however, scrutiny of the AHE dataset revealed a bias towards high GC content at fast-evolving codon partitions (GC bias). By combining several existing approaches aimed at mitigating the impacts of convergence in GC bias, including RY coding and analyses of amino acids, we consistently recovered a strongly supported clade comprised of Holocentridae (squirrelfishes), Berycidae (Alfonsinos), Melamphaidae (bigscale fishes), Cetomimidae (flabby whalefishes), and Rondeletiidae (redmouth whalefishes) as the sister lineage to Percomorpha. Additionally, implementing phylogenetic informativeness (PI) based metrics as a filtration method yielded this same topology, suggesting PI based approaches will preferentially filter these fast-evolving regions and act in a manner consistent with other phylogenetic approaches aimed at mitigating GC bias. Our results provide a new perspective on a key issue for studies investigating the evolutionary history of more than one quarter of all living species of vertebrates.

Expanding anchored hybrid enrichment to resolve both deep and shallow relationships within the spider tree of life

BACKGROUND

Despite considerable effort, progress in spider molecular systematics has lagged behind many other comparable arthropod groups, thereby hindering family-level resolution, classification, and testing of important macroevolutionary hypotheses. Recently, alternative targeted sequence capture techniques have provided molecular systematics a powerful tool for resolving relationships across the Tree of Life. One of these approaches, Anchored Hybrid Enrichment (AHE), is designed to recover hundreds of unique orthologous loci from across the genome, for resolving both shallow and deep-scale evolutionary relationships within non-model systems. Herein we present a modification of the AHE approach that expands its use for application in spiders, with a particular emphasis on the infraorder Mygalomorphae.

RESULTS

Our aim was to design a set of probes that effectively capture loci informative at a diversity of phylogenetic timescales. Following identification of putative arthropod-wide loci, we utilized homologous transcriptome sequences from 17 species across all spiders to identify exon boundaries. Conserved regions with variable flanking regions were then sought across the tick genome, three published araneomorph spider genomes, and raw genomic reads of two mygalomorph taxa. Following development of the 585 target loci in the Spider Probe Kit, we applied AHE across three taxonomic depths to evaluate performance: deep-level spider family relationships (33 taxa, 327 loci); family and generic relationships within the mygalomorph family Euctenizidae (25 taxa, 403 loci); and species relationships in the North American tarantula genus Aphonopelma (83 taxa, 581 loci). At the deepest level, all three major spider lineages (the Mesothelae, Mygalomorphae, and Araneomorphae) were supported with high bootstrap support. Strong support was also found throughout the Euctenizidae, including generic relationships within the family and species relationships within the genus Aptostichus. As in the Euctenizidae, virtually identical topologies were inferred with high support throughout Aphonopelma.

CONCLUSIONS

The Spider Probe Kit, the first implementation of AHE methodology in Class Arachnida, holds great promise for gathering the types and quantities of molecular data needed to accelerate an understanding of the spider Tree of Life by providing a mechanism whereby different researchers can confidently and effectively use the same loci for independent projects, yet allowing synthesis of data across independent research groups.

Methodological congruence in phylogenomic analyses with morphological support for teiid lizards (Sauria: Teiidae)

A well-known issue in phylogenetics is discordance among gene trees, species trees, morphology, and other data types. Gene-tree discordance is often caused by incomplete lineage sorting, lateral gene transfer, and gene duplication. Multispecies-coalescent methods can account for incomplete lineage sorting and are believed by many to be more accurate than concatenation. However, simulation studies and empirical data have demonstrated that concatenation and species tree methods often recover similar topologies. We use three popular methods of phylogenetic reconstruction (one concatenation, two species tree) to evaluate relationships within Teiidae. These lizards are distributed across the United States to Argentina and the West Indies, and their classification has been controversial due to incomplete sampling and the discordance among various character types (chromosomes, DNA, musculature, osteology, etc.) used to reconstruct phylogenetic relationships. Recent morphological and molecular analyses of the group resurrected three genera and created five new genera to resolve non-monophyly in three historically ill-defined genera: Ameiva, Cnemidophorus, and Tupinambis. Here, we assess the phylogenetic relationships of the Teiidae using "next-generation" anchored-phylogenomics sequencing. Our final alignment includes 316 loci (488,656bp DNA) for 244 individuals (56 species of teiids, representing all currently recognized genera) and all three methods (ExaML, MP-EST, and ASTRAL-II) recovered essentially identical topologies. Our results are basically in agreement with recent results from morphology and smaller molecular datasets, showing support for monophyly of the eight new genera. Interestingly, even with hundreds of loci, the relationships among some genera in Tupinambinae remain ambiguous (i.e. low nodal support for the position of Salvator and Dracaena).

A hybrid phylogenetic-phylogenomic approach for species tree estimation in African Agama lizards with applications to biogeography, character evolution, and diversification

Africa is renowned for its biodiversity and endemicity, yet little is known about the factors shaping them across the continent. African Agama lizards (45 species) have a pan-continental distribution, making them an ideal model for investigating biogeography. Many species have evolved conspicuous sexually dimorphic traits, including extravagant breeding coloration in adult males, large adult male body sizes, and variability in social systems among colorful versus drab species. We present a comprehensive time-calibrated species tree for Agama, and their close relatives, using a hybrid phylogenetic-phylogenomic approach that combines traditional Sanger sequence data from five loci for 57 species (146 samples) with anchored phylogenomic data from 215 nuclear genes for 23 species. The Sanger data are analyzed using coalescent-based species tree inference using (*)BEAST, and the resulting posterior distribution of species trees is attenuated using the phylogenomic tree as a backbone constraint. The result is a time-calibrated species tree for Agama that includes 95% of all species, multiple samples for most species, strong support for the major clades, and strong support for most of the initial divergence events. Diversification within Agama began approximately 23 million years ago (Ma), and separate radiations in Southern, East, West, and Northern Africa have been diversifying for >10Myr. A suite of traits (morphological, coloration, and sociality) are tightly correlated and show a strong signal of high morphological disparity within clades, whereby the subsequent evolution of convergent phenotypes has accompanied diversification into new biogeographic areas.