Rapid radiation, ancient incomplete lineage sorting and ancient hybridization in the endemic Lake Tanganyika cichlid tribe Tropheini

Reference Genome Choice and Filtering Thresholds Jointly Influence Phylogenomic Analyses

Molecular phylogenies are a cornerstone of modern comparative biology and are commonly employed to investigate a range of biological phenomena, such as diversification rates, patterns in trait evolution, biogeography, and community assembly. Recent work has demonstrated that significant biases may be introduced into downstream phylogenetic analyses from processing genomic data; however, it remains unclear whether there are interactions among bioinformatic parameters or biases introduced through the choice of reference genome for sequence alignment and variant calling. We address these knowledge gaps by employing a combination of simulated and empirical data sets to investigate the extent to which the choice of reference genome in upstream bioinformatic processing of genomic data influences phylogenetic inference, as well as the way that reference genome choice interacts with bioinformatic filtering choices and phylogenetic inference method. We demonstrate that more stringent minor allele filters bias inferred trees away from the true species tree topology, and that these biased trees tend to be more imbalanced and have a higher center of gravity than the true trees. We find the greatest topological accuracy when filtering sites for minor allele count (MAC) >3-4 in our 51-taxa data sets, while tree center of gravity was closest to the true value when filtering for sites with MAC >1-2. In contrast, filtering for missing data increased accuracy in the inferred topologies; however, this effect was small in comparison to the effect of minor allele filters and may be undesirable due to a subsequent mutation spectrum distortion. The bias introduced by these filters differs based on the reference genome used in short read alignment, providing further support that choosing a reference genome for alignment is an important bioinformatic decision with implications for downstream analyses. These results demonstrate that attributes of the study system and dataset (and their interaction) add important nuance for how best to assemble and filter short-read genomic data for phylogenetic inference.

Turnover of sex chromosomes in the Lake Tanganyika cichlid tribe Tropheini (Teleostei: Cichlidae)

Sex chromosome replacement is frequent in many vertebrate clades, including fish, frogs, and lizards. In order to understand the mechanisms responsible for sex chromosome turnover and the early stages of sex chromosome divergence, it is necessary to study lineages with recently evolved sex chromosomes. Here we examine sex chromosome evolution in a group of African cichlid fishes (tribe Tropheini) which began to diverge from one another less than 4 MYA. We have evidence for a previously unknown sex chromosome system, and preliminary indications of several additional systems not previously reported in this group. We find a high frequency of sex chromosome turnover and estimate a minimum of 14 turnovers in this tribe. We date the origin of the most common sex determining system in this tribe (XY-LG5/19) near the base of one of two major sub-clades of this tribe, about 3.4 MY ago. Finally, we observe variation in the size of one sex-determining region that suggests independent evolution of evolutionary strata in species with a shared sex-determination system. Our results illuminate the rapid rate of sex chromosome turnover in the tribe Tropheini and set the stage for further studies of the dynamics of sex chromosome evolution in this group.

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.

Evaluating a species phylogeny using ddRAD SNPs: Cyto-nuclear discordance and introgression in the salmonid genus Thymallus (Salmonidae)

Hybridization and introgression are very common among freshwater fishes due to the dynamic nature of hydrological landscapes. Cyclic patterns of allopatry and secondary contact provide numerous opportunities for interspecific gene flow, which can lead to discordant paths of evolution for mitochondrial and nuclear genomes. Here, we used double digest restriction-site associated DNA sequencing (ddRADseq) to obtain a genome-wide single nucleotide polymorphism (SNP) dataset comprehensive for allThymallus (Salmonidae)species to infer phylogenetic relationships and evaluate potential recent and historical gene flow among species. The newly obtained nuclear phylogeny was largely concordant with a previously published mitogenome-based topology but revealed a few cyto-nuclear discordances. These incongruencies primarily involved the placement of internal nodes rather than the resolution of species, except for one European species where anthropogenic stock transfers are thought to be responsible for the observed pattern. The analysis of four contact zones where multiple species are found revealed a few cases of mitochondrial capture and limited signals of nuclear introgression. Interestingly, the mechanisms restricting interspecific gene flow might be distinct; while in zones of secondary contact, small-scale physical habitat separation appeared as a limiting factor, biologically based reinforcement mechanisms are presumed to be operative in areas where species presumably evolved in sympatry. Signals of historical introgression were largely congruent with the routes of species dispersal previously inferred from mitogenome data. Overall, the ddRADseq dataset provided a robust phylogenetic reconstruction of the genus Thymallus including new insights into historical hybridization and introgression, opening up new questions concerning their evolutionary history.

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

African cichlids (subfamily: Pseudocrenilabrinae) are among the most diverse vertebrates, and their propensity for repeated rapid radiation has made them a celebrated model system in evolutionary research. Nonetheless, despite numerous studies, phylogenetic uncertainty persists, and riverine lineages remain comparatively underrepresented in higher-level phylogenetic studies. Heterogeneous gene histories resulting from incomplete lineage sorting (ILS) and hybridization are likely sources of uncertainty, especially during episodes of rapid speciation. We investigate relationships of Pseudocrenilabrinae and its close relatives while accounting for multiple sources of genetic discordance using species tree and hybrid network analyses with hundreds of single-copy exons. We improve sequence recovery for distant relatives, thereby extending the taxonomic reach of our probes, with a hybrid reference guided/de novo assembly approach. Our analyses provide robust hypotheses for most higher-level relationships and reveal widespread gene heterogeneity, including in riverine taxa. ILS and past hybridization are identified as sources of genetic discordance in different lineages. Sampling of various Blenniiformes (formerly Ovalentaria) adds strong phylogenomic support for convict blennies (Pholidichthyidae) as sister to Cichlidae, and points to other potentially useful protein-coding markers across the order. A reliable phylogeny with representatives from diverse environments will support ongoing taxonomic and comparative evolutionary research in the cichlid model system.

Phylogenomics of trophically diverse cichlids disentangles processes driving adaptive radiation and repeated trophic transitions

Cichlid fishes of the tribe Tropheini are a striking case of adaptive radiation, exemplifying multiple trophic transitions between herbivory and carnivory occurring in sympatry with other established cichlid lineages. Tropheini evolved highly specialized eco-morphologies to exploit similar trophic niches in different ways repeatedly and rapidly. To better understand the evolutionary history and trophic adaptations of this lineage, we generated a dataset of 532 targeted loci from 21 out of the 22 described Tropheini species. We resolved the Tropheini into seven monophyletic genera and discovered one to be polyphyletic. The polyphyletic genus, Petrochromis, represents three convergent origins of the algae grazing trophic specialization. This repeated evolution of grazing may have been facilitated by adaptive introgression as we found evidence for gene flow among algae grazing genera. We also found evidence of gene flow among algae browsing genera, but gene flow was restricted between herbivorous and carnivorous genera. Furthermore, we observed no evidence supporting a hybrid origin of this radiation. Our molecular evolutionary analyses suggest that opsin genes likely evolved in response to selection pressures associated with trophic ecology in the Tropheini. We found surprisingly little evidence of positive selection in coding regions of jaw-shaping genes in this trophically diverse lineage. This suggests low degrees of freedom for further change in these genes, and possibly a larger role for regulatory variation in driving jaw adaptations. Our study emphasizes Tropheini cichlids as an important model for studying the evolution of trophic specialization and its role in speciation.

Contradictory Phylogenetic Signals in the Laurasiatheria Anomaly Zone

Relationships among laurasiatherian clades represent one of the most highly disputed topics in mammalian phylogeny. In this study, we attempt to disentangle laurasiatherian interordinal relationships using two independent genome-level approaches: (1) quantifying retrotransposon presence/absence patterns, and (2) comparisons of exon datasets at the levels of nucleotides and amino acids. The two approaches revealed contradictory phylogenetic signals, possibly due to a high level of ancestral incomplete lineage sorting. The positions of Eulipotyphla and Chiroptera as the first and second earliest divergences were consistent across the approaches. However, the phylogenetic relationships of Perissodactyla, Cetartiodactyla, and Ferae, were contradictory. While retrotransposon insertion analyses suggest a clade with Cetartiodactyla and Ferae, the exon dataset favoured Cetartiodactyla and Perissodactyla. Future analyses of hitherto unsampled laurasiatherian lineages and synergistic analyses of retrotransposon insertions, exon and conserved intron/intergenic sequences might unravel the conflicting patterns of relationships in this major mammalian clade.

Phylogeographic relationships and the evolutionary history of the Carassius auratus complex with a newly born homodiploid raw fish (2nNCRC)

Background

An important aspect of studying evolution is to understand how new species are formed and their uniqueness is maintained. Hybridization can lead to the formation of new species through reorganization of the adaptive system and significant changes in phenotype. Interestingly, eight stable strains of 2nNCRC derived from interspecies hybridization have been established in our laboratory. To examine the phylogeographical pattern of the widely distributed genus Carassius across Eurasia and investigate the possible homoploid hybrid origin of the Carassius auratus complex lineage in light of past climatic events, the mitochondrial genome (mtDNA) and one nuclear DNA were used to reconstruct the phylogenetic relationship between the C. auratus complex and 2nNCRC and to assess how demographic history, dispersal and barriers to gene flow have led to the current distribution of the C. auratus complex.

Results

As expected, 2nNCRC had a very close relationship with the C. auratus complex and similar morphological characteristics to those of the C. auratus complex, which is genetically distinct from the other three species of Carassius. The estimation of divergence time and ancestral state demonstrated that the C. auratus complex possibly originated from the Yangtze River basin in China. There were seven sublineages of the C. auratus complex across Eurasia and at least four mtDNA lineages endemic to particular geographical regions in China. The primary colonization route from China to Mongolia and the Far East (Russia) occurred during the Late Pliocene, and the diversification of other sublineages of the C. auratus complex specifically coincided with the interglacial stage during the Early and Mid-Pleistocene in China.

Conclusion

Our results support the origin of the C. auratus complex in China, and its wide distribution across Eurasia was mainly due to natural Pleistocene dispersal and recent anthropogenic translocation. The sympatric distribution of the ancestral area for both parents of 2nNCRC and the C. auratus complex, as well as the significant changes in the structure of pharyngeal teeth and morphological characteristics between 2nNCRC and its parents, imply that homoploid hybrid speciation (HHS) for C. auratus could likely have occurred in nature. The diversification pattern indicated an independent evolutionary history of the C. auratus complex, which was not separated from the most recent common ancestor of C. carassius or C. cuvieri. Considering that the paleoclimate oscillation and the development of an eastward-flowing drainage system during the Pliocene and Pleistocene in China provided an opportunity for hybridization between divergent lineages, the formation of 2nNCRC in our laboratory could be a good candidate for explaining the HHS of C. auratus in nature.

OUP accepted manuscript

Whole genome sequences are beginning to revolutionize our understanding of phylogenetic relationships. Yet, even whole genome sequences can fail to resolve the evolutionary history of the most rapidly radiating lineages, where incomplete lineage sorting, standing genetic variation, introgression, and other factors obscure the phylogenetic history of the group. To overcome such challenges, one emerging strategy is to integrate results across different methods. Most such approaches have been implemented on reduced representation genomic data sets, but whole genomes should provide the maximum possible evidence approach. Here, we test the ability of single nucleotide polymorphisms extracted from whole genome resequencing data, implemented in an integrative genomic approach, to resolve key nodes in the phylogeny of the mbuna, rock-dwelling cichlid fishes of Lake Malaŵi, which epitomize the phylogenetic intractability that often accompanies explosive lineage diversification. This monophyletic radiation has diversified at an unparalleled rate into several hundred species in less than 2 million years. Using an array of phylogenomic methods, we consistently recovered four major clades of mbuna, but a large basal polytomy among them. Although introgression between clades apparently contributed to the challenge of phylogenetic reconstruction, reduction of the data set to nonintrogressed sites still did not help to resolve the basal polytomy. On the other hand, relationships among six congeneric species pairs were resolved without ambiguity, even in one case where existing data led us to predict that resolution would be difficult. We conclude that the bursts of diversification at the earliest stages of the mbuna radiation may be phylogenetically unresolvable, but other regions of the tree are phylogenetically clearly supported. Integration of multiple phylogenomic approaches will continue to increase confidence in relationships inferred from these and other whole-genome data sets. [Incomplete lineage sorting; introgression; linkage disequilibrium; multispecies coalescence; rapid radiation; soft polytomy.]

Ultraconserved Elements Improve the Resolution of Difficult Nodes within the Rapid Radiation of Neotropical Sigmodontine Rodents (Cricetidae: Sigmodontinae)

Sigmodontine rodents (Cricetidae, Sigmodontinae) represent the second largest muroid subfamily and the most species-rich group of New World mammals, encompassing above 410 living species and ca. 87 genera. Even with advances on the clarification of sigmodontine phylogenetic relationships that have been made recently, the phylogenetic relationships among the 12 main group of genera (i.e., tribes) remain poorly resolved, in particular among those forming the large clade Oryzomyalia. This pattern has been interpreted as consequence of a rapid radiation upon the group entrance into South America. Here, we attempted to resolve phylogenetic relationships within Sigmodontinae using target capture and high-throughput sequencing of ultraconserved elements (UCEs). We enriched and sequenced UCEs for 56 individuals and collected data from four already available genomes. Analyses of distinct data sets, based on the capture of 4,634 loci, resulted in a highly resolved phylogeny consistent across different methods. Coalescent species-tree based approaches, concatenated matrices, and Bayesian analyses recovered similar topologies that were congruent at the resolution of difficult nodes. We recovered good support for the intertribal relationships within Oryzomyalia; for instance, the tribe Oryzomyini appears as the sister taxa of the remaining oryzomyalid tribes. The estimates of divergence times agree with results of previous studies. We inferred the crown age of the sigmodontine rodents at the end of Middle Miocene, while the main lineages of Oryzomyalia appear to have radiated in a short interval during the Late Miocene. Thus, the collection of a genomic scale data set with a wide taxonomic sampling, provided resolution for the first time of the relationships among the main lineages of Sigmodontinae. We expect the phylogeny presented here will become the backbone for future systematic and evolutionary studies of the group.

The taxonomic diversity of the cichlid fish fauna of ancient Lake Tanganyika, East Africa

Ancient Lake Tanganyika in East Africa houses the world's ecologically and morphologically most diverse assemblage of cichlid fishes, and the third most species-rich after lakes Malawi and Victoria. Despite long-lasting scientific interest in the cichlid species flocks of the East African Great Lakes, for example in the context of adaptive radiation and explosive diversification, their taxonomy and systematics are only partially explored; and many cichlid species still await their formal description. Here, we provide a current inventory of the cichlid fish fauna of Lake Tanganyika, providing a complete list of all valid 208 Tanganyikan cichlid species, and discuss the taxonomic status of more than 50 undescribed taxa on the basis of the available literature as well as our own observations and collections around the lake. This leads us to conclude that there are at least 241 cichlid species present in Lake Tanganyika, all but two are endemic to the basin. We finally summarize some of the major taxonomic challenges regarding Lake Tanganyika's cichlid fauna. The taxonomic inventory of the cichlid fauna of Lake Tanganyika presented here will facilitate future research on the taxonomy and systematics and the ecology and evolution of the species flock, as well as its conservation.

Comparative transcriptomics reveals candidate carotenoid color genes in an East African cichlid fish

BACKGROUND

Carotenoids contribute significantly to animal body coloration, including the spectacular color pattern diversity among fishes. Fish, as other animals, derive carotenoids from their diet. Following uptake, transport and metabolic conversion, carotenoids allocated to body coloration are deposited in the chromatophore cells of the integument. The genes involved in these processes are largely unknown. Using RNA-Sequencing, we tested for differential gene expression between carotenoid-colored and white skin regions of a cichlid fish, Tropheus duboisi "Maswa", to identify genes associated with carotenoid-based integumentary coloration. To control for positional gene expression differences that were independent of the presence/absence of carotenoid coloration, we conducted the same analyses in a closely related population, in which both body regions are white.

RESULTS

A larger number of genes (n = 50) showed higher expression in the yellow compared to the white skin tissue than vice versa (n = 9). Of particular interest was the elevated expression level of bco2a in the white skin samples, as the enzyme encoded by this gene catalyzes the cleavage of carotenoids into colorless derivatives. The set of genes with higher expression levels in the yellow region included genes involved in xanthophore formation (e.g., pax7 and sox10), intracellular pigment mobilization (e.g., tubb, vim, kif5b), as well as uptake (e.g., scarb1) and storage (e.g., plin6) of carotenoids, and metabolic conversion of lipids and retinoids (e.g., dgat2, pnpla2, akr1b1, dhrs). Triglyceride concentrations were similar in the yellow and white skin regions. Extracts of integumentary carotenoids contained zeaxanthin, lutein and beta-cryptoxanthin as well as unidentified carotenoid structures.

CONCLUSION

Our results suggest a role of carotenoid cleavage by Bco2 in fish integumentary coloration, analogous to previous findings in birds. The elevated expression of genes in carotenoid-rich skin regions with functions in retinol and lipid metabolism supports hypotheses concerning analogies and shared mechanisms between these metabolic pathways. Overlaps in the sets of differentially expressed genes (including dgat2, bscl2, faxdc2 and retsatl) between the present study and previous, comparable studies in other fish species provide useful hints to potential carotenoid color candidate genes.

A molecular assessment of species diversity in Tympanopleura and Ageneiosus catfishes (Auchenipteridae: Siluriformes)

In order to test the congruence of genetic data to the morphologically defined Neotropical catfish genera Tympanopleura and Ageneiosus and explore species diversity, we generated 17 DNA barcodes from five of six species of Tympanopleura and 12 of 13 species of Ageneiosus. To discriminate limits between species, an automatic barcode gap discovery (ABGD), a generalised mixed yule-coalescent model (GYMC) and fixed distance thresholds Kimura two-parameter (K2P; 3%) were used to discriminate putative species limits from the DNA barcodes. The ABGD, GMYC and K2P methods agreed by each generating 13 clusters: six in Tympanopleura (five nominal plus one undescribed species) and seven in Ageneiosus. These clusters corresponded broadly to the described species, except in the case of the Ageneiosus ucayalensis group (A. akamai, A. dentatus, A. intrusus, A. ucayalensis, A. uranophthalmus and A. vittatus). Haplotype sharing and low divergences may have prevented molecular methods from distinguishing these species. We hypothesise that this is the result of a recent radiation of a sympatric species group distributed throughout the Amazon Basin. One putative new species of Tympanopleura was also supported by the molecular data. These results taken together highlight the utility of molecular methods such as DNA barcoding in understanding patterns of diversification across large geographic areas and in recognising overlooked diversity.

An examination of introgression and incomplete lineage sorting among three closely related species of chocolate-dipped damselfish (genus: Chromis)

AIM

To determine the impact of ecological and environmental histories on the evolution of coral reef damselfishes at two adjacent marine biogeographic suture zones.

LOCATION

Indo-West Pacific, notably including two suture zones: Socotra and Christmas and Cocos/Keeling Islands.

TAXON

Chromis dimidiata, Chromis margaritifer, and Chromis fieldi.

METHODS

We utilized a combination of nuclear and mitochondrial genetic markers in addition to visual abundance survey data of these fishes.

RESULTS

Despite genetic patterns consistent with incomplete lineage sorting and relatively low genetic differentiation among the three studied Chromis species, there is evidence of hybridization between C. margaritifer and C. fieldi at Christmas Island based on molecular and visual identification. Introgression appears to be spreading westwards to other C. fieldi populations based on COI haplotype comparison. Moreover, the genetic distance between C. margaritifer and C. fieldi suggests that Pleistocene sea-level fluctuations may have contributed to allopatric divergence and secondary contact between these two closely related species.

MAIN CONCLUSIONS

Our study highlights that evolutionary processes in coral reef fishes operate differently between suture zones, possibly due to different ecological and environmental predispositions regulating secondary contact of sister species. While secondary contact likely led to hybridization and introgression at Christmas and Cocos/Keeling Islands, none of those processes seem present at Socotra for the chocolate-dipped damselfish. This difference is likely due to an environmental barrier caused by hydrodynamic regimes in the Gulf of Aden.

East African cichlid lineages (Teleostei: Cichlidae) might be older than their ancient host lakes: new divergence estimates for the east African cichlid radiation

BACKGROUND

Cichlids are a prime model system in evolutionary research and several of the most prominent examples of adaptive radiations are found in the East African Lakes Tanganyika, Malawi and Victoria, all part of the East African cichlid radiation (EAR). In the past, great effort has been invested in reconstructing the evolutionary and biogeographic history of cichlids (Teleostei: Cichlidae). In this study, we present new divergence age estimates for the major cichlid lineages with the main focus on the EAR based on a dataset encompassing representative taxa of almost all recognized cichlid tribes and ten mitochondrial protein genes. We have thoroughly re-evaluated both fossil and geological calibration points, and we included the recently described fossil †Tugenchromis pickfordi in the cichlid divergence age estimates.

RESULTS

Our results estimate the origin of the EAR to Late Eocene/Early Oligocene (28.71 Ma; 95% HPD: 24.43-33.15 Ma). More importantly divergence ages of the most recent common ancestor (MRCA) of several Tanganyika cichlid tribes were estimated to be substantially older than the oldest estimated maximum age of the Lake Tanganyika: Trematocarini (16.13 Ma, 95% HPD: 11.89-20.46 Ma), Bathybatini (20.62 Ma, 95% HPD: 16.88-25.34 Ma), Lamprologini (15.27 Ma; 95% HPD: 12.23-18.49 Ma). The divergence age of the crown haplochromine H-lineage is estimated to 22.8 Ma (95% HPD: 14.40-26.32 Ma) and of the Lake Malawi radiation to 4.07 Ma (95% HDP: 2.93-5.26 Ma). In addition, we recovered a novel lineage within the Lamprologini tribe encompassing only Lamprologus of the lower and central Congo drainage with its divergence estimated to the Late Miocene or early Pliocene. Furthermore we recovered two novel mitochondrial haplotype lineages within the Haplochromini tribe: 'Orthochromis' indermauri and 'Haplochormis' vanheusdeni.

CONCLUSIONS

Divergence time estimates of the MRCA of several Tanganyika cichlid tribes predate the age of the extant Lake Tanganyika basin, and hence are in line with the recently formulated "Melting-Pot Tanganyika" hypothesis. The radiation of the 'Lower Congo Lamprologus clade' might be linked with the Pliocene origin of the modern lower Congo rapids as has been shown for other Lower Congo cichlid assemblages. Finally, the age of origin of the Lake Malawi cichlid flock agrees well with the oldest age estimate for lacustrine conditions in Lake Malawi.

Variations in the breeding behavior of cichlids and the evolution of the multi-functional seminal plasma protein, seminal plasma glycoprotein 120

Background

Seminal plasma proteins are associated with successful fertilization. However, their evolutionary correlation with fertilization mechanisms remains unclear. Cichlids from Lake Tanganyika show a variety-rich spawning behavior that is associated with the transfer of the sperm to the egg for fertilization. One of these behaviors, called “oral fertilization,” emerged during their speciation. In oral fertilization, females nuzzle the milt from male genitalia and pick up the released eggs in their mouths, which are then fertilized inside the oral cavity. Thus, the success of the fertilization is dependent on the retention of sperm in the oral cavity during spawning. Sperm aggregation and immobilization in viscous seminal plasma may help retain the sperm inside the oral cavity, which ultimately determines the success of the fertilization. Seminal plasma glycoprotein 120 (SPP120) is one of the major seminal plasma proteins present in cichlids. SPP120 has been implicated to immobilize sperm and increase the milt viscosity. However, the functional linkage between oral fertilization and seminal plasma proteins has not been investigated.

Results

During trials of simulated oral fertilization, it was observed that milt viscosity contributed to fertilization success by facilitating longer retention of the milt inside the mouth during spawning. Glycosylation of SPP120 was associated with high milt viscosity. Its glycosylation was specifically present in the milt of cichlid species exhibiting oral fertilization. Moreover, recombinant SPP120 from several the oral fertilization species strongly immobilized/aggregated sperm. Therefore, the functions of SPP120 (immobilization/aggregation and its glycosylation) may contribute to success of oral fertilization, and these functions of SPP120 are more prominent in oral fertilization species. In addition, comparative phylogenetic analyses showed a positive evolutionary correlation between SPP120 function and oral fertilization. Hence, these evolutions may have occurred to keep up with the transition in the mode of fertilization. In addition, rapid evolution in the molecular sequence might be associated with functional modifications of SPP120.

Conclusion

These results suggest that SPP120 might be associated with oral fertilization. In other words, reproductive traits that define the interaction between sperms and eggs could be the evolutionary selective force that cause the rapid functional modification of the fertilization-related reproductive protein, SPP120.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1292-0) contains supplementary material, which is available to authorized users.

Whole-genome sequences of Malawi cichlids reveal multiple radiations interconnected by gene flow

The hundreds of cichlid fish species in Lake Malawi constitute the most extensive recent vertebrate adaptive radiation. Here we characterize its genomic diversity by sequencing 134 individuals covering 73 species across all major lineages. The average sequence divergence between species pairs is only 0.1-0.25%. These divergence values overlap diversity within species, with 82% of heterozygosity shared between species. Phylogenetic analyses suggest that diversification initially proceeded by serial branching from a generalist Astatotilapia-like ancestor. However, no single species tree adequately represents all species relationships, with evidence for substantial gene flow at multiple times. Common signatures of selection on visual and oxygen transport genes shared by distantly related deep-water species point to both adaptive introgression and independent selection. These findings enhance our understanding of genomic processes underlying rapid species diversification, and provide a platform for future genetic analysis of the Malawi radiation.

Delineating species along shifting shorelines: Tropheus (Teleostei, Cichlidae) from the southern subbasin of Lake Tanganyika

BACKGROUND

Species delineation is particularly challenging in taxa with substantial intra-specific variation. In systematic studies of fishes, meristics and linear measurements that describe shape are often used to delineate species. Yet, little is known about the taxonomic value of these two types of morphological characteristics. Here, we used Tropheus (Teleostei, Cichlidae) from the southern subbasin of Lake Tanganyika to test which of these types of characters best matched genetic lineages that could represent species in this group of stenotypic rock-dwelling cichlids. We further investigated intra-population variation in morphology. By linking this to a proxy of a population's age, we could assess the evolutionary stability of different kinds of morphological markers.

RESULTS

Morphological data was collected from 570 specimens originating from 86 localities. An AFLP approach revealed the presence of five lineages in the southern subbasin: T. moorii, T. brichardi, T. sp. 'maculatus', T. sp. 'Mpimbwe' and T. sp. 'red', which we consider to represent distinct species. Although both types of morphological data supported this classification, a comparison of PST-values that describe inter-population morphological differentiation, revealed a better correspondence between the taxon delineation based on AFLP data and the patterns revealed by an analysis of meristics than between the AFLP-based taxon delineation and the patterns revealed by an analysis of shape. However, classifying southern populations of Tropheus was inherently difficult as they contained a large amount of clinal variation, both in genetic and in morphological data, and both within and among species. A scenario is put forward to explain the current-day distribution of the species and colour varieties and the observed clinal variation across the subbasin's shoreline. Additionally, we observed that variation in shape was larger in populations from shallow shores whereas populations from steep shores were more variable in meristics. This difference is explained in terms of the different timescales at which small and large scale lake level fluctuations affected populations of littoral cichlids at steep and shallow shores.

CONCLUSIONS

Our results showed meristics to be more evolutionary stable, and of higher taxonomic value for species delimitation in Tropheus, than linear measurements that describe shape. These results should be taken into account when interpreting morphological differences between populations of highly stenotypic species, such as littoral cichlids from the Great East African Lakes.

Phylogenomics uncovers early hybridization and adaptive loci shaping the radiation of Lake Tanganyika cichlid fishes

Lake Tanganyika is the oldest and phenotypically most diverse of the three East African cichlid fish adaptive radiations. It is also the cradle for the younger parallel haplochromine cichlid radiations in Lakes Malawi and Victoria. Despite its evolutionary significance, the relationships among the main Lake Tanganyika lineages remained unresolved, as did the general timescale of cichlid evolution. Here, we disentangle the deep phylogenetic structure of the Lake Tanganyika radiation using anchored phylogenomics and uncover hybridization at its base, as well as early in the haplochromine radiation. This suggests that hybridization might have facilitated these speciation bursts. Time-calibrated trees support that the radiation of Tanganyika cichlids coincided with lake formation and that Gondwanan vicariance concurred with the earliest splits in the cichlid family tree. Genes linked to key innovations show signals of introgression or positive selection following colonization of lake habitats and species' dietary adaptations are revealed as major drivers of colour vision evolution. These findings shed light onto the processes shaping the evolution of adaptive radiations.

Success of cuckoo catfish brood parasitism reflects coevolutionary history and individual experience of their cichlid hosts

Obligate brood parasites manipulate other species into raising their offspring. Avian and insect brood parasitic systems demonstrate how interacting species engage in reciprocal coevolutionary arms races through behavioral and morphological adaptations and counteradaptations. Mouthbrooding cichlid fishes are renowned for their remarkable evolutionary radiations and complex behaviors. In Lake Tanganyika, mouthbrooding cichlids are exploited by the only obligate nonavian vertebrate brood parasite, the cuckoo catfish Synodontis multipunctatus. We show that coevolutionary history and individual learning both have a major impact on the success of cuckoo catfish parasitism between coevolved sympatric and evolutionarily naïve allopatric cichlid species. The rate of cuckoo catfish parasitism in coevolved Tanganyikan hosts was 3 to 11 times lower than in evolutionarily naïve cichlids. Moreover, using experimental infections, we demonstrate that parasite egg rejection in sympatric hosts was much higher, leading to seven times greater parasite survival in evolutionarily naïve than sympatric hosts. However, a high rejection frequency of parasitic catfish eggs by coevolved sympatric hosts came at a cost of increased rejection of their own eggs. A significant cost of catfish parasitism was universal, except for coevolved sympatric cichlid species with previous experience of catfish parasitism, demonstrating that learning and individual experience both contribute to a successful host response.

Disentangling Incomplete Lineage Sorting and Introgression to Refine Species-Tree Estimates for Lake Tanganyika Cichlid Fishes

Adaptive radiation is thought to be responsible for the evolution of a great portion of the past and present diversity of life. Instances of adaptive radiation, characterized by the rapid emergence of an array of species as a consequence to their adaptation to distinct ecological niches, are important study systems in evolutionary biology. However, because of the rapid lineage formation in these groups, and occasional gene flow between the participating species, it is often difficult to reconstruct the phylogenetic history of species that underwent an adaptive radiation. In this study, we present a novel approach for species-tree estimation in rapidly diversifying lineages, where introgression is known to occur, and apply it to a multimarker data set containing up to 16 specimens per species for a set of 45 species of East African cichlid fishes (522 individuals in total), with a main focus on the cichlid species flock of Lake Tanganyika. We first identified, using age distributions of most recent common ancestors in individual gene trees, those lineages in our data set that show strong signatures of past introgression. This led us to formulate three hypotheses of introgression between different lineages of Tanganyika cichlids: the ancestor of Boulengerochromini (or of Boulengerochromini and Bathybatini) received genomic material from the derived H-lineage; the common ancestor of Cyprichromini and Perissodini experienced, in turn, introgression from Boulengerochromini and/or Bathybatini; and the Lake Tanganyika Haplochromini and closely related riverine lineages received genetic material from Cyphotilapiini. We then applied the multispecies coalescent model to estimate the species tree of Lake Tanganyika cichlids, but excluded the lineages involved in these introgression events, as the multispecies coalescent model does not incorporate introgression. This resulted in a robust species tree, in which the Lamprologini were placed as sister lineage to the H-lineage (including the Eretmodini), and we identify a series of rapid splitting events at the base of the H-lineage. Divergence ages estimated with the multispecies coalescent model were substantially younger than age estimates based on concatenation, and agree with the geological history of the Great Lakes of East Africa. Finally, we formally tested the three hypotheses of introgression using a likelihood framework, and find strong support for introgression between some of the cichlid tribes of Lake Tanganyika.

Phylogenomics of pike cichlids (Cichlidae: Crenicichla): the rapid ecological speciation of an incipient species flock

The rapid rise of phenotypic and ecological diversity in independent lake-dwelling groups of cichlids is emblematic of the East African Great Lakes. In this study, we show that similar ecologically based diversification has occurred in pike cichlids (Crenicichla) throughout the Uruguay River drainage of South America. We collected genomic data from nearly 500 ultraconserved element (UCEs) loci and >260 000 base pairs across 33 species, to obtain a phylogenetic hypothesis for the major species groups and to evaluate the relationships and genetic structure among five closely related, endemic, co-occurring species (the Uruguay River species flock; URSF). Additionally, we evaluated ecological divergence of the URSF based on body and lower pharyngeal jaw (LPJ) shape and gut contents. Across the genus, we recovered novel relationships among the species groups. We found strong support for the monophyly of the URSF; however, relationships among these species remain problematic, likely because of the rapid and recent evolution of this clade. Clustered co-ancestry analysis recovered most species as well delimited genetic groups. The URSF species exhibit species-specific body and LPJ shapes associated with specialized trophic roles. Collectively, our results suggest that the URSF consists of incipient species that arose via ecological speciation associated with the exploration of novel trophic roles.

The Role of Alternative Splicing and Differential Gene Expression in Cichlid Adaptive Radiation

Species diverge eco-morphologically through the continuous action of natural selection on functionally important structures, producing alternative adaptive morphologies. In cichlid fishes, the oral and pharyngeal jaws are such key structures. Adaptive variation in jaw morphology contributes to trophic specialization, which is hypothesized to fuel their rapid speciation in the East African Great Lakes. Much is known about the genes involved in cichlid jaw and craniofacial development. However, it is still unclear what salient sources of variation gave rise to trophic-niche specialization, facilitating adaptive radiation. Here, we explore two sources of transcriptional variation that may underlie species-specific disparities in jaw morphology. Using whole transcriptome RNA-sequencing, we analyze differences in gene expression and alternative splicing, at the end of postlarval development, in fully functional jaws of six species of cichlids from the Lake Tanganyika tribe Tropheini. Our data reveal a surprisingly high degree of alternative splicing events compared with gene expression differences among species and trophic types. This suggests that differential trophic adaptation of the jaw apparatus may have been shaped by transcriptional rewiring of splicing as well as gene expression variation during the rapid radiation of the Tropheini. Specifically, genes undergoing splicing across most species were found to be enriched for pharyngeal jaw gene ontology terms. Overall, jaw transcriptional patterns at postlarval developmental stage were highly dynamic and species-specific. In conclusion, this work indicates that shifts in alternative splicing could have played a more important role in cichlid adaptive radiation, and possibly adaptive radiation in general, than currently recognized.

Tempo and rates of diversification in the South American cichlid genus Apistogramma (Teleostei: Perciformes: Cichlidae)

Evaluating biodiversity and understanding the processes involved in diversification are noticeable conservation issues in fishes subject to large, sometimes illegal, ornamental trade purposes. Here, the diversity and evolutionary history of the Neotropical dwarf cichlid genus Apistogramma from several South American countries are investigated. Mitochondrial and nuclear markers are used to infer phylogenetic relationships between 31 genetically identified species. The monophyly of Apistogramma is suggested, and Apistogramma species are distributed into four clades, corresponding to three morphological lineages. Divergence times estimated with the Yule process and an uncorrelated lognormal clock dated the Apistogramma origin to the beginning of the Eocene (≈ 50 Myr) suggesting that diversification might be related to marine incursions. Our molecular dating also suggests that the Quaternary glacial cycles coincide with the phases leading to Apistogramma speciation. These past events did not influence diversification rates in the speciose genus Apistogramma, since diversification appeared low and constant through time. Further characterization of processes involved in recent Apistogramma diversity will be necessary.

Investigating Difficult Nodes in the Placental Mammal Tree with Expanded Taxon Sampling and Thousands of Ultraconserved Elements

The phylogeny of eutherian mammals contains some of the most recalcitrant nodes in the tetrapod tree of life. We combined comprehensive taxon and character sampling to explore three of the most debated interordinal relationships among placental mammals. We performed in silico extraction of ultraconserved element loci from 72 published genomes and invitro enrichment and sequencing of ultraconserved elements from 28 additional mammals, resulting in alignments of 3,787 loci. We analyzed these data using concatenated and multispecies coalescent phylogenetic approaches, topological tests, and exploration of support among individual loci to identify the root of Eutheria and the sister groups of tree shrews (Scandentia) and horses (Perissodactyla). Individual loci provided weak, but often consistent support for topological hypotheses. Although many gene trees lacked accepted species-tree relationships, summary coalescent topologies were largely consistent with inferences from concatenation. At the root of Eutheria, we identified consistent support for a sister relationship between Xenarthra and Afrotheria (i.e., Atlantogenata). At the other nodes of interest, support was less consistent. We suggest Scandentia is the sister of Primatomorpha (Euarchonta), but we failed to reject a sister relationship between Scandentia and Glires. Similarly, we suggest Perissodactyla is sister to Cetartiodactyla (Euungulata), but a sister relationship between Perissodactyla and Chiroptera remains plausible.

Investigating Difficult Nodes in the Placental Mammal Tree with Expanded Taxon Sampling and Thousands of Ultraconserved Elements

The phylogeny of eutherian mammals contains some of the most recalcitrant nodes in the tetrapod tree of life. We combined comprehensive taxon and character sampling to explore three of the most debated interordinal relationships among placental mammals. We performed in silico extraction of ultraconserved element loci from 72 published genomes and invitro enrichment and sequencing of ultraconserved elements from 28 additional mammals, resulting in alignments of 3,787 loci. We analyzed these data using concatenated and multispecies coalescent phylogenetic approaches, topological tests, and exploration of support among individual loci to identify the root of Eutheria and the sister groups of tree shrews (Scandentia) and horses (Perissodactyla). Individual loci provided weak, but often consistent support for topological hypotheses. Although many gene trees lacked accepted species-tree relationships, summary coalescent topologies were largely consistent with inferences from concatenation. At the root of Eutheria, we identified consistent support for a sister relationship between Xenarthra and Afrotheria (i.e., Atlantogenata). At the other nodes of interest, support was less consistent. We suggest Scandentia is the sister of Primatomorpha (Euarchonta), but we failed to reject a sister relationship between Scandentia and Glires. Similarly, we suggest Perissodactyla is sister to Cetartiodactyla (Euungulata), but a sister relationship between Perissodactyla and Chiroptera remains plausible.

Shifting barriers and phenotypic diversification by hybridisation

The establishment of hybrid taxa relies on reproductive isolation from the parental forms, typically achieved by ecological differentiation. Here, we present an alternative mechanism, in which shifts in the strength and location of dispersal barriers facilitate diversification by hybridisation. Our case study concerns the highly diverse, stenotopic rock-dwelling cichlids of the African Great Lakes, many of which display geographic colour pattern variation. The littoral habitat of these fish has repeatedly been restructured in the course of ancient lake level fluctuations. Genetic data and an experimental cross support the hybrid origin of a distinct yellow-coloured variant of Tropheus moorii from ancient admixture between two allopatric, red and bluish variants. Deficient assortative mating preferences imply that reproductive isolation continues to be contingent on geographic separation. Linking paleolimnological data with the establishment of the hybrid variant, we sketch a selectively neutral diversification process governed solely by rearrangements of dispersal barriers.

Genetic dissection of adaptive form and function in rapidly speciating cichlid fishes

Genes of major phenotypic effects and strong genetic correlations can facilitate adaptation, direct selective responses, and potentially lead to phenotypic convergence. However, the preponderance of this type of genetic architecture in repeatedly evolved adaptations remains unknown. Using hybrids between Haplochromis chilotes (thick-lipped) and Pundamilia nyererei (thin-lipped) we investigated the genetics underlying hypertrophied lips and elongated heads, traits that evolved repeatedly in cichlids. At least 25 loci of small-to-moderate and mainly additive effects were detected. Phenotypic variation in lip and head morphology was largely independent. Although several QTL overlapped for lip and head morphology traits, they were often of opposite effects. The distribution of effect signs suggests strong selection on lips. The fitness implications of several detected loci were demonstrated using a laboratory assay testing for the association between genotype and variation in foraging performance. The persistence of low fitness alleles in head morphology appears to be maintained through antagonistic pleiotropy/close linkage with positive-effect lip morphology alleles. Rather than being based on few major loci with strong positive genetic correlations, our results indicate that the evolution of the Lake Victoria thick-lipped ecomorph is the result of selection on numerous loci distributed throughout the genome.

Evidence for environmental and ecological selection in a microbe with no geographic limits to gene flow

The ability for organisms to disperse throughout their environment is thought to strongly influence population structure and thus evolution of diversity within species. A decades-long debate surrounds processes that generate and support high microbial diversity, particularly in the ocean. The debate concerns whether diversification occurs primarily through geographic partitioning (where distance limits gene flow) or through environmental selection, and remains unresolved due to lack of empirical data. Here we show that gene flow in a diatom, an ecologically important eukaryotic microbe, is not limited by global-scale geographic distance. Instead, environmental and ecological selection likely play a more significant role than dispersal in generating and maintaining diversity. We detected significantly diverged populations (F_ST > 0.130) and discovered temporal genetic variability at a single site that was on par with spatial genetic variability observed over distances of 15,000 km. Relatedness among populations was decoupled from geographic distance across the global ocean and instead, correlated significantly with water temperature and whole-community chlorophyll a Correlations with temperature point to the importance of environmental selection in structuring populations. Correlations with whole-community chlorophyll a, a proxy for autotrophic biomass, suggest that ecological selection via interactions with other plankton may generate and maintain population genetic structure in marine microbes despite global-scale dispersal. Here, we provide empirical evidence for global gene flow in a marine eukaryotic microbe, suggesting that everything holds the potential to be everywhere, with environmental and ecological selection rather than geography or dispersal dictating the structure and evolution of diversity over space and time.

Genomewide SNP data reveal cryptic phylogeographic structure and microallopatric divergence in a rapids-adapted clade of cichlids from the Congo River

The lower Congo River is a freshwater biodiversity hot spot in Africa characterized by some of the world's largest rapids. However, little is known about the evolutionary forces shaping this diversity, which include numerous endemic fishes. We investigated phylogeographic relationships in Teleogramma, a small clade of rheophilic cichlids, in the context of regional geography and hydrology. Previous studies have been unable to resolve phylogenetic relationships within Teleogramma due to lack of variation in nuclear genes and discrete morphological characters among putative species. To sample more broadly across the genome, we analysed double-digest restriction-associated sequencing (ddRAD) data from 53 individuals across all described species in the genus. We also assessed body shape and mitochondrial variation within and between taxa. Phylogenetic analyses reveal previously unrecognized lineages and instances of microallopatric divergence across as little as ~1.5 km. Species ranges appear to correspond to geographic regions broadly separated by major hydrological and topographic barriers, indicating these features are likely important drivers of diversification. Mitonuclear discordance indicates one or more introgressive hybridization events, but no clear evidence of admixture is present in nuclear genomes, suggesting these events were likely ancient. A survey of female fin patterns hints that previously undetected lineage-specific patterning may be acting to reinforce species cohesion. These analyses highlight the importance of hydrological complexity in generating diversity in certain freshwater systems, as well as the utility of ddRAD-Seq data in understanding diversification processes operating both below and above the species level.

Phylogenetic analyses of the subgenus Mollienesia (Poecilia, Poeciliidae, Teleostei) reveal taxonomic inconsistencies, cryptic biodiversity, and spatio-temporal aspects of diversification in Middle America

The subgenus Mollienesia is a diverse group of freshwater fishes, including species that have served as important models across multiple biological disciplines. Nonetheless, the taxonomic history of this group has been conflictive and convoluted, in part because the evolutionary relationships have not been rigorously resolved. We conducted a comprehensive molecular phylogenetic analysis of the subgenus Mollienesia to identify taxonomic discrepancies and potentially identify undescribed species, estimate ancestral areas of origin and estimate dates of divergence, as well as explore biogeographical patterns. Our findings confirm the presence of three main clades composed of the P. latipinna, P. sphenops, and P. mexicana species complexes. Unlike previously hypothesized morphology-based analyses, species found on the Caribbean Islands are not part of Mollienesia, but are more closely related to species of the subgenus Limia. Our study also revealed several taxonomic inconsistencies and distinct lineages in the P. mexicana species complex that may represent undescribed species. The diversity in the subgenus Mollienesia is a result of dynamic geologic activity leading to vicariant events, dispersal across geologic blocks, and ecological speciation.

A robust phylogeny among major lineages of the East African cichlids

The huge monophyletic group of the East African cichlid radiations (EAR) consists of thousands of species belonging to 12-14 tribes; the number of tribes differs among studies. Many studies have inferred phylogenies of EAR tribes using various genetic markers. However, these phylogenies partly contradict one another and can have weak statistic support. In this study, we conducted maximum-likelihood (ML) phylogenetic analyses using restriction site-associated DNA (RAD) sequences and propose a new robust phylogenetic hypothesis among Lake Tanganyika cichlid fishes, which cover most EAR tribes. Data matrices can vary in size and contents depending on the strategies used to process RAD sequences. Therefore, we prepared 23 data matrices with various processing strategies. The ML phylogenies inferred from 15 large matrices (2.0×10(6) to 1.1×10(7) base pairs) resolved every tribe as a monophyletic group with 100% bootstrap support and shared the same topology regarding relationships among the tribes. Most nodes among the tribes were supported by 100% bootstrap values, and the bootstrap support for the other node varied among the 15 ML trees from 70% to 100%. These robust ML trees differ partly in topology from those in earlier studies, and these phylogenetic relationships have important implications for the tribal classification of EAR.

First insights into the diversity of gill monogeneans of 'Gnathochromis' and Limnochromis (Teleostei, Cichlidae) in Burundi: do the parasites mirror host ecology and phylogenetic history?

Monogenea is one of the most species-rich groups of parasitic flatworms worldwide, with many species described only recently, which is particularly true for African monogeneans. For example, Cichlidogyrus, a genus mostly occurring on African cichlids, comprises more than 100 nominal species. Twenty-two of these have been described from Lake Tanganyika, a famous biodiversity hotspot in which many vertebrate and invertebrate taxa, including monogeneans, underwent unique and spectacular radiations. Given their often high degrees of host specificity, parasitic monogeneans were also used as a potential tool to uncover host species relationships. This study presents the first investigation of the monogenean fauna occurring on the gills of endemic ‘Gnathochromis’ species along the Burundese coastline of Lake Tanganyika. We test whether their monogenean fauna reflects the different phylogenetic position and ecological niche of ‘Gnathochromis’ pfefferi and Gnathochromis permaxillaris. Worms collected from specimens of Limnochromis auritus, a cichlid belonging to the same cichlid tribe as G. permaxillaris, were used for comparison. Morphological as well as genetic characterisation was used for parasite identification. In total, all 73 Cichlidogyrus individuals collected from ‘G.’ pfefferi were identified as C. irenae. This is the only representative of Cichlidogyrus previously described from ‘G.’ pfefferi, its type host. Gnathochromis permaxillaris is infected by a species of Cichlidogyrus morphologically very similar to C. gillardinae. The monogenean species collected from L. auritus is considered as new for science, but sample size was insufficient for a formal description. Our results confirm previous suggestions that ‘G.’ pfefferi as a good disperser is infected by a single monogenean species across the entire Lake Tanganyika. Although G. permaxillaris and L. auritus are placed in the same tribe, Cichlidogyrus sp. occurring on G. permaxillaris is morphologically more similar to C. irenae from ‘G.’ pfefferi, than to the Cichlidogyrus species found on L. auritus. Various evolutionary processes, such as host-switching or duplication events, might underlie the pattern observed in this particular parasite-host system. Additional samples for the Cichlidogyrus species occuring on G. permaxillaris and L. auritus are needed to unravel their evolutionary history by means of (co-)phylogenetic analyses.

HEXT, a software supporting tree-based screens for hybrid taxa in multilocus data sets, and an evaluation of the homoplasy excess test

The homoplasy excess test (HET) is a tree-based screen for hybrid taxa in multilocus nuclear phylogenies. Homoplasy between a hybrid taxon and the clades containing the parental taxa reduces bootstrap support in the tree. The HET is based on the expectation that excluding the hybrid taxon from the data set increases the bootstrap support for the parental clades, whereas excluding non-hybrid taxa has little effect on statistical node support. To carry out a HET, bootstrap trees are calculated with taxon-jackknife data sets, that is excluding one taxon (species, population) at a time. Excess increase in bootstrap support for certain nodes upon exclusion of a particular taxon indicates the hybrid (the excluded taxon) and its parents (the clades with increased support).We introduce a new software program, hext, which generates the taxon-jackknife data sets, runs the bootstrap tree calculations, and identifies excess bootstrap increases as outlier values in boxplot graphs. hext is written in r language and accepts binary data (0/1; e.g. AFLP) as well as co-dominant SNP and genotype data.We demonstrate the usefulness of hext in large SNP data sets containing putative hybrids and their parents. For instance, using published data of the genus Vitis (~6,000 SNP loci), hext output supports V. × champinii as a hybrid between V. rupestris and V. mustangensis.With simulated SNP and AFLP data sets, excess increases in bootstrap support were not always connected with the hybrid taxon (false positives), whereas the expected bootstrap signal failed to appear on several occasions (false negatives). Potential causes for both types of spurious results are discussed.With both empirical and simulated data sets, the taxon-jackknife output generated by hext provided additional signatures of hybrid taxa, including changes in tree topology across trees, consistent effects of exclusions of the hybrid and the parent taxa, and moderate (rather than excessive) increases in bootstrap support. hext significantly facilitates the taxon-jackknife approach to hybrid taxon detection, even though the simple test for excess bootstrap increase may not reliably identify hybrid taxa in all applications.

Hidden biodiversity in an ancient lake: phylogenetic congruence between Lake Tanganyika tropheine cichlids and their monogenean flatworm parasites

The stunning diversity of cichlid fishes has greatly enhanced our understanding of speciation and radiation. Little is known about the evolution of cichlid parasites. Parasites are abundant components of biodiversity, whose diversity typically exceeds that of their hosts. In the first comprehensive phylogenetic parasitological analysis of a vertebrate radiation, we study monogenean parasites infecting tropheine cichlids from Lake Tanganyika. Monogeneans are flatworms usually infecting the body surface and gills of fishes. In contrast to many other parasites, they depend only on a single host species to complete their lifecycle. Our spatially comprehensive combined nuclear-mitochondrial DNA dataset of the parasites covering almost all tropheine host species (N = 18), reveals species-rich parasite assemblages and shows consistent host-specificity. Statistical comparisons of host and parasite phylogenies based on distance and topology-based tests demonstrate significant congruence and suggest that host-switching is rare. Molecular rate evaluation indicates that species of Cichlidogyrus probably diverged synchronically with the initial radiation of the tropheines. They further diversified through within-host speciation into an overlooked species radiation. The unique life history and specialisation of certain parasite groups has profound evolutionary consequences. Hence, evolutionary parasitology adds a new dimension to the study of biodiversity hotspots like Lake Tanganyika.

Form, function and phylogeny: comparative morphometrics of Lake Tanganyika's cichlid tribe Tropheini

Lake Tanganyika's cichlid fishes represent one of the most diverse species assemblages of the world. In this study we focused on the tribe Tropheini which occupies several trophic niches, mostly in rocky habitats. We analysed morphological variation of seventeen closely related species by means of geometric morphometric methods and related these data to ecological characteristics and phylogeny of the study species. It turned out that morphology mostly correlated well with ecological parameters, but not always closely with the degree of the phylogenetic relatedness of the species. Overall, body shapes in the tribe Tropheini are of great evolutionary plasticity, but variation is restricted to particular body parts: the preorbital region once again emerged as a key factor that facilitated their impressive radiation.

Lake Tanganyika--a 'melting pot' of ancient and young cichlid lineages (Teleostei: Cichlidae)?

A long history of research focused on the East Africa cichlid radiations (EAR) revealed discrepancies between mtDNA and nuclear phylogenies, suggesting that interspecific hybridisation may have been significant during the radiation of these fishes. The approximately 250 cichlid species of Lake Tanganyika have their roots in a monophyletic African cichlid assemblage, but controversies remain about the precise phylogenetic origin and placement of different lineages and consequently about L. Tanganyika colonization scenarios. 3312 AFLP loci and the mitochondrial ND2 gene were genotyped for 91 species representing almost all major lacustrine and riverine haplotilapiine east African cichlid lineages with a focus on L. Tanganyika endemics. Explicitly testing for the possibility of ancient hybridisation events, a comprehensive phylogenetic network hypothesis is proposed for the origin and diversification of L. Tanganyika cichlids. Inference of discordant phylogenetic signal strongly suggests that the genomes of two endemic L. Tanganyika tribes, Eretmodini and Tropheini, are composed of an ancient mixture of riverine and lacustrine lineages. For the first time a strong monophyly signal of all non-haplochromine mouthbrooding species endemic to L. Tanganyika (“ancient mouthbrooders”) was detected. Further, in the genomes of early diverging L. Tanganyika endemics Trematocarini, Bathybatini, Hemibatini and Boulengerochromis genetic components of other lineages belonging to the East African Radiation appear to be present. In combination with recent palaeo-geological results showing that tectonic activity in the L. Tanganyika region resulted in highly dynamic and heterogeneous landscape evolution over the Neogene and Pleistocene, the novel phylogenetic data render a single lacustrine basin as the geographical cradle of the endemic L. Tanganyika cichlid lineages unlikely. Instead a scenario of a pre-rift origin of several independent L. Tanganyika precursor lineages which diversified in ancient rivers and precursor lakes and then amalgamated in the extant L. Tanganyika basin is put forward as an alternative: the 'melting pot Tanganyika' hypothesis.

Morphology, molecules, and monogenean parasites: an example of an integrative approach to cichlid biodiversity

The unparalleled biodiversity of Lake Tanganyika (Africa) has fascinated biologists for over a century; its unique cichlid communities are a preferred model for evolutionary research. Although species delineation is, in most cases, relatively straightforward, higher-order classifications were shown not to agree with monophyletic groups. Here, traditional morphological methods meet their limitations. A typical example are the tropheine cichlids currently belonging to Simochromis and Pseudosimochromis. The affiliations of these widespread and abundant cichlids are poorly understood. Molecular work suggested that genus and species boundaries should be revised. Moreover, previous morphological results indicated that intraspecific variation should be considered to delineate species in Lake Tanganyika cichlids. We review the genera Simochromis and Pseudosimochromis using an integrative approach. Besides a morphometric study and a barcoding approach, monogenean Cichlidogyrus (Platyhelminthes: Ancyrocephalidae) gill parasites, often highly species-specific, are used as complementary markers. Six new species are described. Cichlidogyrus raeymaekersi sp. nov., C. muterezii sp. nov. and C. banyankimbonai sp. nov. infect S. diagramma. Cichlidogyrus georgesmertensi sp. nov. was found on S. babaulti and S. pleurospilus, C. franswittei sp. nov. on both S. marginatus and P. curvifrons and C. frankwillemsi sp. nov. only on P. curvifrons. As relatedness between Cichlidogyrus species usually reflects relatedness between hosts, we considered Simochromis monotypic because the three Cichlidogyrus species found on S. diagramma belonged to a different morphotype than those found on the other Simochromis. The transfer of S. babaulti, S. marginatus, S. pleurospilus and S. margaretae to Pseudosimochromis was justified by the similarity of their Cichlidogyrus fauna and the intermediate morphology of S. margaretae. Finally parasite data also supported the synonymy between S. pleurospilus and S. babaulti, a species that contains a large amount of geographical morphological variation.

Asymmetric dominance and asymmetric mate choice oppose premating isolation after allopatric divergence

Assortative mating promotes reproductive isolation and allows allopatric speciation processes to continue in secondary contact. As mating patterns are determined by mate preferences and intrasexual competition, we investigated male–male competition and behavioral isolation in simulated secondary contact among allopatric populations. Three allopatric color morphs of the cichlid fish Tropheus were tested against each other. Dyadic male–male contests revealed dominance of red males over bluish and yellow-blotch males. Reproductive isolation in the presence of male–male competition was assessed from genetic parentage in experimental ponds and was highly asymmetric among pairs of color morphs. Red females mated only with red males, whereas the other females performed variable degrees of heteromorphic mating. Discrepancies between mating patterns in ponds and female preferences in a competition-free, two-way choice paradigm suggested that the dominance of red males interfered with positive assortative mating of females of the subordinate morphs and provoked asymmetric hybridization. Between the nonred morphs, a significant excess of negative assortative mating by yellow-blotch females with bluish males did not coincide with asymmetric dominance among males. Hence, both negative assortative mating preferences and interference of male–male competition with positive assortative preferences forestall premating isolation, the latter especially in environments unsupportive of competition-driven spatial segregation.