Exon-based phylogenomics and the relationships of African cichlid fishes: tackling the challenges of reconstructing phylogenies with repeated rapid radiations

Introgression Underlies Phylogenetic Uncertainty But Not Parallel Plumage Evolution in a Recent Songbird Radiation

Instances of parallel phenotypic evolution offer great opportunities to understand the evolutionary processes underlying phenotypic changes. However, confirming parallel phenotypic evolution and studying its causes requires a robust phylogenetic framework. One such example is the "black-and-white wagtails," a group of 5 species in the songbird genus Motacilla: 1 species, Motacilla alba, shows wide intra-specific plumage variation, while the 4r others form 2 pairs of very similar-looking species (M. aguimp + M. samveasnae and M. grandis + M. maderaspatensis, respectively). However, the 2 species in each of these pairs were not recovered as sisters in previous phylogenetic inferences. Their relationships varied depending on the markers used, suggesting that gene tree heterogeneity might have hampered accurate phylogenetic inference. Here, we use whole genome resequencing data to explore the phylogenetic relationships within this group, with a special emphasis on characterizing the extent of gene tree heterogeneity and its underlying causes. We first used multispecies coalescent methods to generate a "complete evidence" phylogenetic hypothesis based on genome-wide variants, while accounting for incomplete lineage sorting (ILS) and introgression. We then investigated the variation in phylogenetic signal across the genome to quantify the extent of discordance across genomic regions and test its underlying causes. We found that wagtail genomes are mosaics of regions supporting variable genealogies, because of ILS and inter-specific introgression. The most common topology across the genome, supporting M. alba and M. aguimp as sister species, appears to be influenced by ancient introgression. Additionally, we inferred another ancient introgression event, between M. alba and M. grandis. By combining results from multiple analyses, we propose a phylogenetic network for the black-and-white wagtails that confirms that similar phenotypes evolved in non-sister lineages, supporting parallel plumage evolution. Furthermore, the inferred reticulations do not connect species with similar plumage coloration, suggesting that introgression does not underlie parallel plumage evolution in this group. Our results demonstrate the importance of investing genome-wide patterns of gene tree heterogeneity to help understand the mechanisms underlying phenotypic evolution. [Gene tree heterogeneity; incomplete lineage sorting; introgression; parallel evolution; phylogenomics; plumage evolution; wagtails.].

Disentangling historical relationships within Poeciliidae (Teleostei: Cyprinodontiformes) using ultraconserved elements

Poeciliids (Cyprinodontiformes: Poeciliidae), commonly known as livebearers, are popular fishes in the aquarium trade (e.g., guppies, mollies, swordtails) that are widely distributed in the Americas, with 274 valid species in 27 genera. This group has undergone various taxonomic changes recently, spurred by investigations using traditional genetic markers. Here we used over 1,000 ultraconserved loci to infer the relationships within Poeciliidae in the first attempt at understanding their diversification based on genome-scale data. We explore gene tree discordance and investigate potential incongruence between concatenation and coalescent inference methods. Our aim is to examine the influence of incomplete lineage sorting and reticulate evolution on the poeciliids' evolutionary history and how these factors contribute to the observed gene tree discordace. Our concatenated and coalescent phylogenomic inferences recovered four major clades within Poeciliidae. Most supra-generic level relationships we inferred were congruent with previous molecular studies, but we found some disagreements; the Middle American taxa Phallichthys and Poecilia (Mollienesia) were recovered as non-monophyletic, and unlike other recent molecular studies, we recovered Brachyrhaphis as monophyletic. Our study is the first to provide signatures of reticulate evolution in Poeciliidae at the family level; however, continued finer-scale investigations are needed to understand the complex evolutionary history of the family along with a much-needed taxonomic re-evaluation.

Comparison of seven complete mitochondrial genomes from Lamprologus and Neolamprologus (Chordata, Teleostei, Perciformes) and the phylogenetic implications for Cichlidae

In this study, mitochondrial genomes (mitogenomes) of seven cichlid species (Lamprologuskungweensis, L.meleagris, L.ornatipinnis, Neolamprologusbrevis, N.caudopunctatus, N.leleupi, and N.similis) are characterized for the first time. The newly sequenced mitogenomes contained 37 typical genes [13 protein-coding genes (PCGs), two ribosomal RNA genes (rRNAs) and 22 transfer RNA genes (tRNAs)]. The mitogenomes were 16,562 ~ 16,587 bp in length with an A + T composition of 52.1~58.8%. The cichlid mitogenomes had a comparable nucleotide composition, A + T content was higher than the G + C content. The AT-skews of most mitogenomes were inconspicuously positive and the GC-skews were negative, indicating higher occurrences of C than G. Most PCGs started with the conventional start codon, ATN. There was no essential difference in the codon usage patterns of these seven species. Using Ka/Ks, we found the fastest-evolving gene were atp8. But the results of p-distance indicated that the fastest-evolving gene was nad6. Phylogenetic analysis revealed that L.meleagris did not cluster with Lamprologus species, but with species from the genus Neolamprologus. The novel information obtained about these mitogenomes will contribute to elucidating the complex relationships among cichlid species.

Ancient standing genetic variation facilitated the adaptive radiation of Lake Victoria cichlids

Cichlid fishes are textbook examples of explosive speciation and adaptive radiation, providing a great opportunity to understand how the genomic substrate yields extraordinary species diversity. Recently, we performed comparative genomic analyses of three Lake Victoria cichlids to reveal the genomic substrates underlying their rapid speciation and adaptation. We found that long divergent haplotypes derived from large-scale standing genetic variation, which originated before the adaptive radiation of Lake Victoria cichlids, may have contributed to their rapid diversification. In addition, the present study on genomic data from other East African cichlids suggested the reuse of alleles that may have originated in the ancestral lineages of Lake Tanganyika cichlids during cichlid evolution. Therefore, our results highlight that the primary factor that could drive repeated adaptive radiation across East African cichlids was allelic reuse from standing genetic variation to adapt to their own specific environment. In this report, we summarize the main results and discuss the evolutionary mechanisms of cichlids, based on our latest findings.

Mitogenomic Characterization of Cameroonian Endemic Coptodon camerunensis (Cichliformes: Cichlidae) and Matrilineal Phylogeny of Old-World Cichlids

The mitogenomic evolution of old-world cichlids is still largely incomplete in Western Africa. In this present study, the complete mitogenome of the Cameroon endemic cichlid, Coptodon camerunensis, was determined by next-generation sequencing. The mitogenome was 16,557 bp long and encoded with 37 genes (13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes, and a control region). The C. camerunensis mitogenome is AT-biased (52.63%), as exhibited in its congener, Coptodon zillii (52.76% and 53.04%). The majority of PCGs start with an ATG initiation codon, except COI, which starts with a GTG codon and five PCGs and ends with the TAA termination codon and except seven PCGs with an incomplete termination codon. In C. camerunensis mitogenome, most tRNAs showed classical cloverleaf secondary structures, except tRNA-serine with a lack of DHU stem. Comparative analyses of the conserved blocks of two Coptodonini species control regions revealed that the CSB-II block was longer than other blocks and contained highly variable sites. Using 13 concatenated PCGs, the mitogenome-based Bayesian phylogeny easily distinguished all the examined old-world cichlids. Except for Oreochromini and Coptodinini tribe members, the majority of the taxa exhibited monophyletic clustering within their respective lineages. C. camerunensis clustered closely with Heterotilapia buttikoferi (tribe Heterotilapiini) and had paraphyletic clustering with its congener, C. zillii. The Oreochromini species also displayed paraphyletic grouping, and the genus Oreochromis showed a close relationship with Coptodinini and Heterotilapiini species. In addition, illustrating the known distribution patterns of old-world cichlids, the present study is congruent with the previous hypothesis and proclaims that prehistoric geological evolution plays a key role in the hydroclimate of the African continent during Mesozoic, which simultaneously disperses and/or colonizes cichlids in different ichthyological provinces and Rift Lake systems in Africa. The present study suggests that further mitogenomes of cichlid species are required, especially from western Africa, to understand their unique evolution and adaptation.