Lessons from a natural experiment: Allopatric morphological divergence and sympatric diversification in the Midas cichlid species complex are largely influenced by ecology in a deterministic way

Meristic co-evolution and genomic co-localization of lateral line scales and vertebrae in Central American cichlid fishes

Meristic traits are often treated as distinct phenotypes that can be used to differentiate and delineate recently diverged species. For instance, the number of lateral line scales and vertebrae, two traits that vary substantially among Neotropical Heroine cichlid species, have been previously suggested to co-evolve. These meristic traits could co-evolve due to shared adaptive, developmental, or genetic factors. If they were found to be genetically or developmentally non-independent, this might require a more general re-evaluation of their role in evolutionary or taxonomic studies. We expanded a previous analysis of correlated evolution of meristic traits (lateral line scales and vertebrae counts) in these fishes to include a range of phylogenetic reconstructions as well as the analyses of 13 Nicaraguan Midas cichlid species (Amphilophus spp.). Additionally, we performed qualitative traits locus (QTL) mapping in a F2 laboratory-reared hybrid population from two ecologically divergent Midas cichlid fish species to discover and evaluate whether genomic co-segregation might explain the observed patterns of meristic co-evolution. Meristic values for these traits were found to morphologically differentiate some species of the Midas cichlid adaptive radiation. Our QTL analysis pinpointed several genomic regions associated with divergence in these traits and highlighted the potential for genomic co-segregation of the lateral line and vertebrae numbers on two chromosomes. Further, our phylogenetic comparative analyses consistently recovered a significant positive evolutionary correlation between the counts of lateral line scale and vertebrae numbers in Neotropical cichlids. Hence, the findings of genomic co-segregation could partially explain the co-evolution of these two meristic traits in these species. Continuing to unravel the genetic architecture governing meristic divergence helps to better understand both trait correlations and the utility of meristic traits in taxonomic diagnoses and how traits in phenotypes might be expected to co-evolve.

Repeated Divergence in Opsin Gene Expression Mirrors Photic Habitat Changes in Rapidly Evolving Crater Lake Cichlid Fishes

AbstractSelection pressures differ along environmental gradients, and traits tightly linked to fitness (e.g., the visual system) are expected to track such variation. Along gradients, adaptation to local conditions might be due to heritable and nonheritable environmentally induced variation. Disentangling these sources of phenotypic variation requires studying closely related populations in nature and in the laboratory. The Nicaraguan lakes represent an environmental gradient in photic conditions from clear crater lakes to very turbid great lakes. From two old, turbid great lakes, Midas cichlid fish (Amphilophus cf. citrinellus) independently colonized seven isolated crater lakes of varying light conditions, resulting in a small adaptive radiation. We estimated variation in visual sensitivities along this photic gradient by measuring cone opsin gene expression among lake populations. Visual sensitivities observed in all seven derived crater lake populations shifted predictably in direction and magnitude, repeatedly mirroring changes in photic conditions. Comparing wild-caught and laboratory-reared fish revealed that 48% of this phenotypic variation is genetically determined and evolved rapidly. Decreasing intrapopulation variation as environments become spectrally narrower suggests that different selective landscapes operate along the gradient. We conclude that the power to predict phenotypic evolution along gradients depends on both the magnitude of environmental change and the selective landscape shape.

Camera Traps Uncover the Behavioral Ecology of an Endemic, Cryptic Monkey Species in the Congo Basin

Guenons are the most diverse clade of African primates, and many species living within the core of the Congo Basin rainforest are still understudied. The recently described guenon species, Cercopithecus lomamiensis, known as lesula, is a cryptic, semi-terrestrial species endemic to the central Congo Basin in the Democratic Republic of the Congo. The recent IUCN Red List Assessment recognizes lesula's risk of extinction in the wild as Vulnerable. The objective of our study was to use camera traps to expand knowledge on the behavioral ecology of lesula. We conducted three systematic, terrestrial camera trap (CT) surveys within Lomami National Park and buffer zone (Okulu: 2013; Losekola: 2014; E15: 2015). We accumulated 598 independent events of lesula over 5960 CT days from 92 CTs. Typical of Cercopithecus species, camera trap videos reveal that lesula has a diurnal activity pattern, birth seasonality, a group size of up to 32 individuals, and social organization with female philopatry and male dispersal. Results also suggest that lesula are highly terrestrial, distinguishing them from other Cercopithecus species, which are mostly arboreal. Our study provides new information about the behavioral ecology of this little-studied primate, generating species-specific knowledge of a threatened species for successful conservation planning.

Divergent and non-parallel evolution of MHC IIB in the Neotropical Midas cichlid species complex

Background

Ecological diversification is the result of divergent natural selection by contrasting habitat characteristics that favours the evolution of distinct phenotypes. This process can happen in sympatry and in allopatry. Habitat-specific parasite communities have the potential to drive diversification among host populations by imposing selective pressures on their host's immune system. In particular, the hyperdiverse genes of the major histocompatibility complex (MHC) are implicated in parasite-mediated host divergence. Here, we studied the extent of divergence at MHC, and discuss how it may have contributed to the Nicaraguan Midas cichlid species complex diversification, one of the most convincing examples of rapid sympatric parallel speciation.

Results

We genotyped the MHC IIB for individuals from six sympatric Midas cichlid assemblages, each containing species that have adapted to exploit similar habitats. We recovered large allelic and functional diversity within the species complex. While most alleles were rare, functional groups of alleles (supertypes) were common, suggesting that they are key to survival and that they were maintained during colonization and subsequent radiations. We identified lake-specific and habitat-specific signatures for both allelic and functional diversity, but no clear pattern of parallel divergence among ecomorphologically similar phenotypes.

Conclusions

Colonization and demographic effects of the fish could have contributed to MHC evolution in the Midas cichlid in conjunction with habitat-specific selective pressures, such as parasites associated to alternative preys or environmental features. Additional ecological data will help evaluating the role of host–parasite interactions in the Midas cichlid radiations and aid in elucidating the potential role of non-parallel features differentiating crater lake species assemblages.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-01997-9.

Thyroid hormone tinkering elicits integrated phenotypic changes potentially explaining rapid adaptation of color vision in cichlid fish

Vision is critical for most vertebrates, including fish. One challenge that aquatic habitats pose is the high variability in spectral properties depending on depth and the inherent optical properties of the water. By altering opsin gene expression and chromophore usage, cichlid fish modulate visual sensitivities to maximize sensory input from the available light in their respective habitat. Thyroid hormone (TH) has been proposed to play a role in governing adaptive diversification in visual sensitivity in Nicaraguan Midas cichlids, which evolved in less than 4000 generations. As suggested by indirect measurements of TH levels (i.e., expression of deiodinases), populations adapted to short wavelength light in clear lakes have lower TH levels than ones inhabiting turbid lakes enriched in long-wavelength light. We experimentally manipulated TH levels by exposing 2-week-old Midas cichlids to exogenous TH or a TH inhibitor and measured opsin gene expression and chromophore usage (via cyp27c1 expression). Although exogenous TH induces long-wavelength sensitivity by changing opsin gene expression and chromophore usage in a concerted manner, TH-inhibited fish exhibit a visual phenotype with sensitivities shifted to shorter wavelengths. Tinkering with TH levels in eyes results in concerted phenotypic changes that can provide a rapid mechanism of adaptation to novel light environments.

The macroparasite fauna of cichlid fish from Nicaraguan lakes, a model system for understanding host-parasite diversification and speciation

The Nicaraguan lakes represent an ideal continent-island-like setting to study the colonization patterns of both fish and their parasites. The dominant fish fauna are cichlids, particularly the Midas cichlid species complex Amphilophus spp., a well-studied model for recent sympatric speciation. Here, we characterized the Midas cichlid macroparasite diversity in Nicaraguan lakes. We evaluated patterns of parasite diversity across host populations. Morphological and molecular analyses were conducted, revealing a macroparasite fauna composed by 37 taxa, including platyhelminths, nematodes, copepods, branchiurans, hirudineans and oribatids. Three invasive species are reported for the first time. The Midas cichlid was infected by 22 parasite taxa, 18 shared with other cichlids. Eight taxa conformed the core parasite fauna of the Midas cichlid. The large lakes had higher parasite diversity than the smaller and isolated crater lakes, although parasite infracommunity diversity was lower. Environmental factors along with the differential distribution of intermediate hosts, the potential resistance gained by their hosts after colonization of new lakes, competitive exclusion among parasites, or the introduction of exotic fish, may determine the observed pattern of parasite heterogeneous distribution. Our study provides a ground to explore the evolutionary history of both, hosts and parasites within the context of speciation and diversification processes.

Benefits and limitations of a new genome-based PCR-RFLP genotyping assay (GB-RFLP): A SNP-based detection method for identification of species in extremely young adaptive radiations

High-throughput DNA sequencing technologies make it possible now to sequence entire genomes relatively easily. Complete genomic information obtained by whole-genome resequencing (WGS) can aid in identifying and delineating species even if they are extremely young, cryptic, or morphologically difficult to discern and closely related. Yet, for taxonomic or conservation biology purposes, WGS can remain cost-prohibitive, too time-consuming, and often constitute a "data overkill." Rapid and reliable identification of species (and populations) that is also cost-effective is made possible by species-specific markers that can be discovered by WGS. Based on WGS data, we designed a PCR restriction fragment length polymorphism (PCR-RFLP) assay for 19 Neotropical Midas cichlid populations (Amphilophus cf. citrinellus), that includes all 13 described species of this species complex. Our work illustrates that identification of species and populations (i.e., fish from different lakes) can be greatly improved by designing genetic markers using available "high resolution" genomic information. Yet, our work also shows that even in the best-case scenario, when whole-genome resequencing information is available, unequivocal assignments remain challenging when species or populations diverged very recently, or gene flow persists. In summary, we provide a comprehensive workflow on how to design RFPL markers based on genome resequencing data, how to test and evaluate their reliability, and discuss the benefits and pitfalls of our approach.

Dual function and associated costs of a highly exaggerated trait in a cichlid fish

Exaggerated secondary sexual characteristics are apparently costly and seem to defy natural selection. This conundrum promoted the theory of sexual selection. Accordingly, exaggerated secondary sexual characteristics might be ornaments on which female choice is based and/or armaments used during male-male competition. Males of many cichlid fish species, including the adaptive radiation of Nicaraguan Midas cichlids, develop a highly exaggerated nuchal hump, which is thought to be a sexually selected trait. To test this hypothesis, we conducted a series of behavioral assays in F2 hybrids obtained from crossing a species with a relatively small hump and one with an exaggerated hump. Mate-choice experiments showed a clear female preference for males with large humps. In an open-choice experiment with limited territories, couples including large humped males were more successful in acquiring these territories. Therefore, nuchal humps appear to serve dual functions as an ornament for attracting mates and as an armament for direct contest with rivals. Although being beneficial in terms of sexual selection, this trait also imposes fitness costs on males possessing disproportionally large nuchal humps since they exhibit decreased endurance and increased energetic costs when swimming. We conclude that these costs illustrate trade-offs associated with large hump size between sexual and natural selection, which causes the latter to limit further exaggeration of this spectacular male trait.

Rapid adaptive radiation in a hillstream cyprinid fish in the East African White Nile River basin

Adaptive radiation of freshwater fishes was long thought to be possible only in lacustrine environments. Recently, several studies have shown that riverine and stream environments also provide the ecological opportunity for adaptive radiation. In this study, we report on a riverine adaptive radiation of six ecomorphs of cyprinid hillstream fishes of the genus Garra in a river located in the Ethiopian Highlands in East Africa. Garra are predominantly highly specialized algae-scrapers with a wide distribution ranging from Southeast Asia to West Africa. However, adaptive phenotypic diversification in mouth type, sucking disc morphology, gut length and body shape have probably been found among these ecomorphs in a single Ethiopian river. Moreover, we found two novel phenotypes of Garra ("thick-lipped" and "predatory") that had not been discovered before in this species-rich genus (>160 species). Mitochondrial and genome-wide data suggest monophyletic, intrabasin evolution of Garra phenotypic diversity with signatures of gene flow from other local populations. Although sympatric ecomorphs are genetically distinct and can be considered to being young species as suggested by genome-wide single nucleotide polymorphism data, mitochondrial DNA was unable to identify any genetic structure suggesting recent and rapid speciation events. Some data suggest a hybrid origin of the novel "thick-lipped" ecomorph. Here we highlight how, driven by ecological opportunity, an ancestral trophically highly specialized lineage is likely to have rapidly radiated in a riverine environment promoted by the evolution of novel feeding strategies.

The comparative genomic landscape of adaptive radiation in crater lake cichlid fishes

Factors ranging from ecological opportunity to genome composition might explain why only some lineages form adaptive radiations. While being rare, particular systems can provide natural experiments within an identical ecological setting where species numbers and phenotypic divergence in two closely related lineages are notably different. We investigated one such natural experiment using two de novo assembled and 40 resequenced genomes and asked why two closely related Neotropical cichlid fish lineages, the Amphilophus citrinellus species complex (Midas cichlids; radiating) and Archocentrus centrarchus (Flyer cichlid; nonradiating), have resulted in such disparate evolutionary outcomes. Although both lineages inhabit many of the same Nicaraguan lakes, whole-genome inferred demography suggests that priority effects are not likely to be the cause of the dissimilarities. Also, genome-wide levels of selection, transposable element dynamics, gene family expansion, major chromosomal rearrangements and the number of genes under positive selection were not markedly different between the two lineages. To more finely investigate particular subsets of the genome that have undergone adaptive divergence in Midas cichlids, we also examined if there was evidence for 'molecular pre-adaptation' in regions identified by QTL mapping of repeatedly diverging adaptive traits. Although most of our analyses failed to pinpoint substantial genomic differences, we did identify functional categories containing many genes under positive selection that provide candidates for future studies on the propensity of Midas cichlids to radiate. Our results point to a disproportionate role of local, rather than genome-wide factors underlying the propensity for these cichlid fishes to adaptively radiate.

Contrasting signatures of genomic divergence during sympatric speciation

The transition from 'well-marked varieties' of a single species into 'well-defined species'-especially in the absence of geographic barriers to gene flow (sympatric speciation)-has puzzled evolutionary biologists ever since Darwin1,2. Gene flow counteracts the buildup of genome-wide differentiation, which is a hallmark of speciation and increases the likelihood of the evolution of irreversible reproductive barriers (incompatibilities) that complete the speciation process3. Theory predicts that the genetic architecture of divergently selected traits can influence whether sympatric speciation occurs4, but empirical tests of this theory are scant because comprehensive data are difficult to collect and synthesize across species, owing to their unique biologies and evolutionary histories5. Here, within a young species complex of neotropical cichlid fishes (Amphilophus spp.), we analysed genomic divergence among populations and species. By generating a new genome assembly and re-sequencing 453 genomes, we uncovered the genetic architecture of traits that have been suggested to be important for divergence. Species that differ in monogenic or oligogenic traits that affect ecological performance and/or mate choice show remarkably localized genomic differentiation. By contrast, differentiation among species that have diverged in polygenic traits is genomically widespread and much higher overall, consistent with the evolution of effective and stable genome-wide barriers to gene flow. Thus, we conclude that simple trait architectures are not always as conducive to speciation with gene flow as previously suggested, whereas polygenic architectures can promote rapid and stable speciation in sympatry.

Genetic Variation and Hybridization in Evolutionary Radiations of Cichlid Fishes

Evolutionary radiations are responsible for much of the variation in biodiversity across taxa. Cichlid fishes are well known for spectacular evolutionary radiations, as they have repeatedly evolved into large and phenotypically diverse arrays of species. Cichlid genomes carry signatures of past events and, at the same time, are the substrate for ongoing evolution. We survey genome-wide data and the available literature covering 438 cichlid populations (412 species) across multiple radiations to synthesize information about patterns and sharing of genetic variation. Nucleotide diversity within species is low in cichlids, with 92% of surveyed populations having less diversity than the median value found in other vertebrates. Divergence within radiations is also low, and a large proportion of variation is shared among species due to incomplete lineage sorting and widespread hybridization. Population genetics therefore provides a suitable conceptual framework for evolutionary genomic studies of cichlid radiations. We focus in detail on the roles of hybridization in shaping the patterns of genetic variation and in promoting cichlid diversification.

Parallel and non-parallel changes of the gut microbiota during trophic diversification in repeated young adaptive radiations of sympatric cichlid fish

BACKGROUND

Recent increases in understanding the ecological and evolutionary roles of microbial communities have underscored the importance of their hosts' biology. Yet, little is known about gut microbiota dynamics during the early stages of ecological diversification and speciation. We sequenced the V4 region of the 16s rRNA gene to study the gut microbiota of Nicaraguan Midas cichlid fish (Amphilophus cf. citrinellus). Specifically, we tested the hypothesis that parallel divergence in trophic ecology in extremely young adaptive radiations from two crater lakes is associated with parallel changes of their gut microbiota.

RESULTS

Bacterial communities of fish guts and lake water were highly distinct, indicating that the gut microbiota is shaped by host-specific factors. Among individuals of the same crater lake, differentiation in trophic ecology was weakly associated with gut microbiota differentiation, suggesting that diet, to some extent, affects the gut microbiota. However, differences in trophic ecology were much more pronounced across than within species whereas similar patterns were not observed for taxonomic and functional differences of the gut microbiota. Across the two crater lakes, we could not detect conclusive evidence for parallel changes of the gut microbiota associated with trophic ecology.

CONCLUSIONS

A lack of clearly differentiated niches during the early stages of ecological diversification might result in non-parallel changes of gut microbial communities, as observed in our study system as well as in other recently diverged fish species. Video Abstract.

Microbial co-occurrence networks of gut microbiota reveal community conservation and diet-associated shifts in cichlid fishes

Background

The extent to which deterministic rather than stochastic processes guide gut bacteria co-existence and ultimately their assembling into a community remains largely unknown. Co-occurrence networks of bacterial associations offer a powerful approach to begin exploring gut microbial community structure, maintenance and dynamics, beyond compositional aspects alone. Here we used an iconic model system, the cichlid fishes, with their multiple lake assemblages and extraordinary ecological diversity, to investigate a) patterns of microbial associations that were robust to major phylogeographical variables, and b) changes in microbial network structure along dietary shifts. We tackled these objectives using the large gut microbiota sequencing dataset available (nine lakes from Africa and America), building geographical and diet-specific networks and performing comparative network analyses.

Results

Major findings indicated that lake and continental microbial networks were highly resembling in global topology and node taxonomic composition, despite the heterogeneity of the samples. A small fraction of the observed co-occurrences among operational taxonomic units (OTUs) was conserved across all lake assemblages. These were all positive associations and involved OTUs within the genera Cetobacterium and Turicibacter and several OTUs belonging to the families of Peptostreptococcaceae and Clostridiaceae (order Clostridiales). Mapping of diet contribution on the African Lake Tanganyika network (therefore excluding the geographic variable) revealed a clear community change from carnivores (C) to omnivores (O) to herbivores (H). Node abundances and effect size for pairwise comparisons between diets supported a strong contrasting pattern between C and H. Moreover, diet-associated nodes in H formed complex modules of positive interactions among taxonomically diverse bacteria (mostly Verrucomicrobia and Proteobacteria).

Conclusions

Conservation of microbial network topologies and specific bacterial associations across distinct lake assemblages point to a major host-associated effect and potential deterministic processes shaping the cichlid gut microbiota. While the origin and biological relevance of these common associations remain unclear, their persistence suggests an important functional role in the cichlid gut. Among the very diverse cichlids of L. Tanganyika, diet nonetheless represents a major driver of microbial community changes. By intersecting results from predictive network inferences and experimental trials, future studies will be directed to explore the strength of these associations, predict the outcome of community alterations driven by diet and ultimately help understanding the role of gut microbiota in cichlid trophic diversification.

Evolutionary dynamics of pre- and postzygotic reproductive isolation in cichlid fishes

Cichlid fishes are exceptionally species-rich, speciated at explosive rates and, hence, are a model system in speciation research. Yet, their reproductive isolating barriers have, so far, not been comprehensively studied. Here, we review current knowledge on pre- and postzygotic mechanisms in cichlids. While premating isolation is the norm in cichlids, its strength varies across lineages and with the geographical setting. Moreover, manipulations of ambient conditions tended to reduce assortative mating among closely related species, suggesting that premating isolation in cichlids is often fragile and context dependent. The observed lack of complete reproductive isolation is supported by past and present hybridization events that have contributed to diversity by creating novel allelic combinations. On the other hand, our meta-analysis highlights that intrinsic postzygotic isolation might accumulate faster than assumed. Mild forms of genetic incompatibilities, such as sex ratio distortion, can already be observed among closely related species. Therefore, cessation of gene flow by strong reproductive isolation in cichlids requires a combination of premating prezygotic isolation supplemented with intrinsic and extrinsic postzygotic barriers. Further, we suggest crucial next steps to improve our knowledge about reproductive barriers in cichlids to understand the evolutionary dynamics of pre- and postzygotic isolation mechanisms during adaptive radiations. This article is part of the theme issue 'Towards the completion of speciation: the evolution of reproductive isolation beyond the first barriers'.

Parallelism in eco-morphology and gene expression despite variable evolutionary and genomic backgrounds in a Holarctic fish

Understanding the extent to which ecological divergence is repeatable is essential for predicting responses of biodiversity to environmental change. Here we test the predictability of evolution, from genotype to phenotype, by studying parallel evolution in a salmonid fish, Arctic charr (Salvelinus alpinus), across eleven replicate sympatric ecotype pairs (benthivorous-planktivorous and planktivorous-piscivorous) and two evolutionary lineages. We found considerable variability in eco-morphological divergence, with several traits related to foraging (eye diameter, pectoral fin length) being highly parallel even across lineages. This suggests repeated and predictable adaptation to environment. Consistent with ancestral genetic variation, hundreds of loci were associated with ecotype divergence within lineages of which eight were shared across lineages. This shared genetic variation was maintained despite variation in evolutionary histories, ranging from postglacial divergence in sympatry (ca. 10-15kya) to pre-glacial divergence (ca. 20-40kya) with postglacial secondary contact. Transcriptome-wide gene expression (44,102 genes) was highly parallel across replicates, involved biological processes characteristic of ecotype morphology and physiology, and revealed parallelism at the level of regulatory networks. This expression divergence was not only plastic but in part genetically controlled by parallel cis-eQTL. Lastly, we found that the magnitude of phenotypic divergence was largely correlated with the genetic differentiation and gene expression divergence. In contrast, the direction of phenotypic change was mostly determined by the interplay of adaptive genetic variation, gene expression, and ecosystem size. Ecosystem size further explained variation in putatively adaptive, ecotype-associated genomic patterns within and across lineages, highlighting the role of environmental variation and stochasticity in parallel evolution. Together, our findings demonstrate the parallel evolution of eco-morphology and gene expression within and across evolutionary lineages, which is controlled by the interplay of environmental stochasticity and evolutionary contingencies, largely overcoming variable evolutionary histories and genomic backgrounds.

The physiology of rapid ecological specialization: A look at the Midas cichlids

Understanding the process of speciation is a primary goal of evolutionary biology, yet the question of whether speciation can reach completion in the presence of gene flow remains controversial. For more than 50 years, the cichlids of Africa, and more recently those in South and Central America, have served as model systems for the study of speciation in nature. Cichlids are distinguished by their enormous species richness, their diversity of behavioural and trophic adaptations, and their rapid rate of divergence. In both Africa and South and Central America, the repeated interaction of geology, new founder events and adaptive evolution has created a series of natural experiments with speciation occurring both within and between waterbodies of differing ages. In the "From the Cover" paper in this issue of the Journal of Molecular Ecology, Raffini, Schneider, Franchini, Kautt and Meyer move beyond the question of which mechanisms drive speciation, and instead show that divergent morphologies and physiologies translate into adaptive traits. They investigate differences in physiology and gene expression profiles in a benthic/limnetic species pair of Midas cichlidsin a 24,000-year-old Nicaraguan crater lake. While recently diverged, these two species demonstrate significant ecological, but limited genetic differentiation. The authors find that the distinct morphotypes translate into relevant differences in swimming performance and metabolic rates that correspond to differential gene expression profiles. Hence, the authors take an integrative approach examining the impacts of morphological differences on performance and niche partitioning: an approach that can advance our understanding of the drivers of morphological and physiological divergence during speciation.

Implementing Large Genomic Single Nucleotide Polymorphism Data Sets in Phylogenetic Network Reconstructions: A Case Study of Particularly Rapid Radiations of Cichlid Fish

The Midas cichlids of the Amphilophus citrinellus spp. species complex from Nicaragua (13 species) are an extraordinary example of adaptive and rapid radiation ($<$24,000 years old). These cichlids are a very challenging group to infer its evolutionary history in phylogenetic analyses, due to the apparent prevalence of incomplete lineage sorting (ILS), as well as past and current gene flow. Assuming solely a vertical transfer of genetic material from an ancestral lineage to new lineages is not appropriate in many cases of genes transferred horizontally in nature. Recently developed methods to infer phylogenetic networks under such circumstances might be able to circumvent these problems. These models accommodate not just ILS, but also gene flow, under the multispecies network coalescent (MSNC) model, processes that are at work in young, hybridizing, and/or rapidly diversifying lineages. There are currently only a few programs available that implement MSNC for estimating phylogenetic networks. Here, we present a novel way to incorporate single nucleotide polymorphism (SNP) data into the currently available PhyloNetworks program. Based on simulations, we demonstrate that SNPs can provide enough power to recover the true phylogenetic network. We also show that it can accurately infer the true network more often than other similar SNP-based programs (PhyloNet and HyDe). Moreover, our approach results in a faster algorithm compared to the original pipeline in PhyloNetworks, without losing power. We also applied our new approach to infer the phylogenetic network of Midas cichlid radiation. We implemented the most comprehensive genomic data set to date (RADseq data set of 679 individuals and $>$37K SNPs from 19 ingroup lineages) and present estimated phylogenetic networks for this extremely young and fast-evolving radiation of cichlid fish. We demonstrate that the MSNC is more appropriate than the multispecies coalescent alone for the analysis of this rapid radiation. [Genomics; multispecies network coalescent; phylogenetic networks; phylogenomics; RADseq; SNPs.]

Phylogeography and Ecological Niche Shape the Cichlid Fish Gut Microbiota in Central American and African Lakes

Cichlid fishes, with their repeated colonization of lakes and subsequent radiations at different scales of phylogenetic and ecological diversification, offer an excellent model system to understand the factors shaping the host-gut microbiota association in nature. Here, we characterized the gut microbiota of the Amphilophus species complex from Central America (known as the Midas cichlid complex), encompassing 158 wild specimens (13 species) collected from seven Nicaraguan lakes, and combined these data with previously published data from two African lakes (spanning 29 species). Our aim was to comprehensively explore trends in microbiota variation and persistence along the large spatial and temporal scales of cichlid diversification (from the oldest radiation in L. Tanganyika, 9-12 My old, to young ones in Nicaraguan crater lakes, <0.5 My old), in allopatry and sympatry (within and across lakes), and across the range of dietary niches (from highly specialized to generalist feeders). Despite their extraordinary diversity, cichlids shared a remarkably conserved microbial taxonomic profile, which argues for a primary role of the host genetics in the assembly and maintenance of these microbial communities. Within this partly constrained microbiota profile, geographic isolation (continent and lake) represented the first level of discrimination. For the Midas cichlid, a partial congruency was found between host microbiota and genetic distances, suggesting that microbial communities have partly diversified along their cichlid phylogeographic history of crater lake colonization. In sympatry (within lakes), the young and poorly ecologically diversified cichlid assemblages of Central American lakes display largely unresolved gut microbiotas (in terms of both alpha and beta diversities), whereas the phylogenetically and ecologically diverse species found in African lakes showed greater microbial interspecific diversity. This pattern largely points to the level of habitat segregation, trophic niche overlap, and reproductive barriers as major modulators of the gut microbiota connectivity among sympatric species.

John ME. How to Investigate the Origins of Novelty: Insights Gained from Genetic, Behavioral, and Fitness Perspectives

Biologists are drawn to the most extraordinary adaptations in the natural world, often referred to as evolutionary novelties, yet rarely do we understand the microevolutionary context underlying the origins of novel traits, behaviors, or ecological niches. Here we discuss insights gained into the origins of novelty from a research program spanning biological levels of organization from genotype to fitness in Caribbean pupfishes. We focus on a case study of the origins of novel trophic specialists on San Salvador Island, Bahamas and place this radiation in the context of other rapid radiations. We highlight questions that can be addressed about the origins of novelty at different biological levels, such as measuring the isolation of novel phenotypes on the fitness landscape, locating the spatial and temporal origins of adaptive variation contributing to novelty, detecting dysfunctional gene regulation due to adaptive divergence, and connecting behaviors with novel traits. Evolutionary novelties are rare, almost by definition, and we conclude that integrative case studies can provide insights into this rarity relative to the dynamics of adaptation to more common ecological niches and repeated parallel speciation, such as the relative isolation of novel phenotypes on fitness landscapes and the transient availability of ecological, genetic, and behavioral opportunities.

Convergent phenotypic evolution of the visual system via different molecular routes: How Neotropical cichlid fishes adapt to novel light environments

How predictable is evolution? This remains a fundamental but contested issue in evolutionary biology. When independent lineages colonize the same environment, we are presented with a natural experiment that allows us to ask if genetic and ecological differences promote species-specific evolutionary outcomes or whether species phenotypically evolve in a convergent manner in response to shared selection pressures. If so, are the molecular mechanisms underlying phenotypic convergence the same? In Nicaragua, seven species of cichlid fishes concurrently colonized two novel photic environments. Hence, their visual system represents a compelling model to address these questions, particularly since the adaptive value of phenotypic changes is well-understood. By analyzing retinal transcriptomes, we found that differential expression of genes responsible for color vision (cone opsins and cyp27c1) produced rapid and mostly convergent changes of predicted visual sensitivities. Notably, these changes occurred in the same direction in all species although there were differences in underlying gene expression patterns illustrating nonconvergence at the molecular level. Adaptive phenotypes evolved deterministically, even when species differ substantially in ecology and genetic variation. This provides strong evidence that phenotypic evolution of the visual system occurred in response to similar selective forces of the photic environment.