Ancient DNA is preserved in fish fossils from tropical lake sediments.

Tropical freshwater lakes are well known for their high biodiversity, and particularly the East African Great Lakes are renowned for their adaptive radiation of cichlid fishes. While comparative phylogenetic analyses of extant species flocks have revealed patterns and processes of their diversification, little is known about evolutionary trajectories within lineages, the impacts of environmental drivers, or the scope and nature of now-extinct diversity. Time-structured palaeodata from geologically young fossil records, such as fossil counts and particularly ancient DNA (aDNA) data, would help fill this large knowledge gap. High ambient temperatures can be detrimental to the preservation of DNA, but refined methodology now allows data generation even from very poorly preserved samples. Here, we show for the first time that fish fossils from tropical lake sediments yield endogenous aDNA. Despite generally low endogenous content and high sample dropout, the application of high-throughput sequencing and, in some cases, sequence capture allowed taxonomic assignment and phylogenetic placement of 17% of analysed fish fossils to family or tribe level, including remains which are up to 2700 years old or weigh less than 1 mg. The relationship between aDNA degradation and the thermal age of samples is similar to that described for terrestrial samples from cold environments when adjusted for elevated temperature. Success rates and aDNA preservation differed between the investigated lakes Chala, Kivu and Victoria, possibly caused by differences in bottom water oxygenation. Our study demonstrates that the sediment records of tropical lakes can preserve genetic information on rapidly diversifying fish taxa over time scales of millennia.

First characterization of PIWI-interacting RNA clusters in a cichlid fish with a B chromosome.

BACKGROUND: B chromosomes are extra elements found in several eukaryote species. Usually, they do not express a phenotype in the host. However, advances in bioinformatics over the last decades have allowed us to describe several genes and molecular functions related to B chromosomes. These advances enable investigations of the relationship between the B chromosome and the host to understand how this element has been preserved in genomes. However, considering that transposable elements (TEs) are highly abundant in this supernumerary chromosome, there is a lack of knowledge concerning the dynamics of TE control in B-carrying cells. Thus, the present study characterized PIWI-interacting RNA (piRNA) clusters and pathways responsible for silencing the mobilization of TEs in gonads of the cichlid fish Astatotilapia latifasciata carrying the B chromosome. RESULTS: Through small RNA-seq and genome assembly, we predicted and annotated piRNA clusters in the A. latifasciata genome for the first time. We observed that these clusters had biased expression related to sex and the presence of the B chromosome. Furthermore, three piRNA clusters, named curupira, were identified in the B chromosome. Two of them were expressed exclusively in gonads of samples with the B chromosome. The composition of these curupira sequences was derived from LTR, LINE, and DNA elements, representing old and recent transposition events in the A. latifasciata genome and the B chromosome. The presence of the B chromosome also affected the expression of piRNA pathway genes. The mitochondrial cardiolipin hydrolase-like (pld6) gene is present in the B chromosome, as previously reported, and an increase in its expression was detected in gonads with the B chromosome. CONCLUSIONS: Due to the high abundance of TEs in the B chromosome, it was possible to investigate the origin of piRNA from these jumping genes. We hypothesize that the B chromosome has evolved its own genomic guardians to prevent uncontrolled TE mobilization. Furthermore, we also detected an expression bias in the presence of the B chromosome over A. latifasciata piRNA clusters and pathway genes.

Feeding and locomotor rhythm of the ornamental cichlid fish Amatitlania sp. in different social groups.

The aim of this study was to evaluate the feeding activity and the rhythm of daily locomotor activity of the convict cichlid (Amatitlania sp.) kept in different social groups under a self-feeding system. A total of 120 animals was distributed among six repetitions of four social groups, as follows: group 1 with one male and one female per tank; group 2 with three males and three females per tank; group 3 with six males per tank; and group 4 with six females per tank. Feeding activity (FA) and locomotor activity (LA) were evaluated using photoelectric presence-sensors connected to automatic feeders. The fish were fed a commercial extruded diet (46% crude protein and 3600 kcal kg-1 of digestible energy). Animal growth was evaluated for all groups. After 30 days of experimentation, the fish stabilized their demands by adjusting their consumption. Amatitlania sp. showed predominantly diurnal FA and LA. All groups showed a peak of activity when the light was turned on and when it was turned off. In summary, FA and LA of Amatitlania sp. are predominantly diurnal and independent of social group. Pairs and groups of males and females together consume less food in relation to groups of one sex or the other due to reproductive behaviour. On the other hand, groups of only males or females consume more food because they lack reproductive stimuli and thus prioritize growth. These results may support good feeding management practices for this ornamental cichlid. Studies relating feeding behaviour with different social groups are of great importance for determining effective feeding strategies for this species in captivity. Thus, such a study assists in a more efficient production of Amatitlania sp.

Neural activity patterns differ between learning contexts in a social fish.

Learning and decision-making are greatly influenced by context. When navigating a complex social world, individuals must quickly ascertain where to gain important resources and which group members are useful sources of such information. Such dynamic behavioural processes require neural mechanisms that are flexible across contexts. Here we examine how the social context influences the learning response during a cue discrimination task and the neural activity patterns that underlie acquisition of this novel information. Using the cichlid fish, Astatotilapia burtoni, we show that learning of the task is faster in social groups than in a non-social context. We quantify the neural activity patterns by examining the expression of Fos, an immediate-early gene, across brain regions known to play a role in social behaviour and learning (such as the putative teleost homologues of the mammalian hippocampus, basolateral amygdala and medial amygdala/BNST complex). We find that neural activity patterns differ between social and non-social contexts. Taken together, our results suggest that while the same brain regions may be involved in the learning of a cue association, the activity in each region reflects an individual's social context.

Aggressive interactions and consistency of dominance hierarchies of the native and nonnative cichlid fishes of the Balsas basin.

Dominance hierarchies are generally established based on the levels of aggressiveness that animals present. Frequently, animals fight to establish a dominance hierarchy and obtain a disputed resource. The Mexican mojarra Cichlasoma istlanum is a native species of the Balsas river basin and coexists there with four nonnative cichlids: tilapia Oreochromis sp., convict cichlid Amatitlania nigrofasciata, spotcheek cichlid Thorichthys maculipinnis, and green terror Andinoacara rivulatus. These five cichlid species compete for spaces for reproduction, feeding, and shelter and frequently engage in aggressive interactions to obtain these resources. We quantified dominance indices to evaluate the hierarchical structure of dominance among these five cichlids and the duration of aggressive behaviors of the Mexican mojarra during experimental contests between the native species and each of the four nonnative species. The Mexican mojarra was consistently dominant over the other four cichlid species, performing a larger number of aggressive behaviors and investing more time in attacking than the nonnative cichlids, which resulted in a higher hierarchical position. Our results show that the native fish, Mexican mojarra, established dominance over all four nonnative cichlid fish of the Balsas basin. Thus, the establishment of nonnative cichlid species in the Balsas basin is likely associated with factors other than behavioral dominance.

Mapping of quantitative trait loci underlying a magic trait in ongoing ecological speciation.

BACKGROUND: Telmatochromis temporalis is a cichlid fish endemic to Lake Tanganyika. The normal and dwarf morphs of this fish are a clear example of ongoing ecological speciation, and body size plays an important role in this speciation event as a magic trait. However, the genetic basis underlying this trait has not been studied. RESULTS: Based on double-digested restriction-site associated DNA (ddRAD) sequencing of a hybrid cross between the morphs that includes F0 male, F0 female, and 206 F2 individuals, we obtained a linkage map consisting of 708 ddRAD markers in 22 linkage groups, which corresponded to the previously reported Oreochromis niloticus chromosomes, and identified one significant and five suggestive quantitative trait loci (QTL) for body size. From the body-size distribution pattern, the significant and three of the five suggestive QTL are possibly associated with genes responsible for the difference in body size between the morphs. CONCLUSIONS: The QTL analysis presented here suggests that multiple genes, rather than a single gene, control morph-specific body size. The present results provide further insights about the genes underlying the morph specific body size and evolution of the magic trait during ecological speciation.

Ancestral Hybridization Facilitated Species Diversification in the Lake Malawi Cichlid Fish Adaptive Radiation.

The adaptive radiation of cichlid fishes in East African Lake Malawi encompasses over 500 species that are believed to have evolved within the last 800,000 years from a common founder population. It has been proposed that hybridization between ancestral lineages can provide the genetic raw material to fuel such exceptionally high diversification rates, and evidence for this has recently been presented for the Lake Victoria region cichlid superflock. Here, we report that Lake Malawi cichlid genomes also show evidence of hybridization between two lineages that split 3-4 Ma, today represented by Lake Victoria cichlids and the riverine Astatotilapia sp. "ruaha blue." The two ancestries in Malawi cichlid genomes are present in large blocks of several kilobases, but there is little variation in this pattern between Malawi cichlid species, suggesting that the large-scale mosaic structure of the genomes was largely established prior to the radiation. Nevertheless, tens of thousands of polymorphic variants apparently derived from the hybridization are interspersed in the genomes. These loci show a striking excess of differentiation across ecological subgroups in the Lake Malawi cichlid assemblage, and parental alleles sort differentially into benthic and pelagic Malawi cichlid lineages, consistent with strong differential selection on these loci during species divergence. Furthermore, these loci are enriched for genes involved in immune response and vision, including opsin genes previously identified as important for speciation. Our results reinforce the role of ancestral hybridization in explosive diversification by demonstrating its significance in one of the largest recent vertebrate adaptive radiations.

De novo genome assembly of the cichlid fish Astatotilapia latifasciata reveals a higher level of genomic polymorphism and genes related to B chromosomes.

Supernumerary B chromosomes (Bs) are accessory elements to the regular chromosome set (As) and have been observed in a huge diversity of eukaryotic species. Although extensively investigated, the biological significance of Bs remains enigmatic. Here, we present de novo genome assemblies for the cichlid fish Astatotilapia latifasciata, a well-known model to study Bs. High coverage data with Illumina sequencing was obtained for males and females with 0B (B-), 1B, and 2B (B+) chromosomes to provide information regarding the diversity among these genomes. The draft assemblies comprised 771 Mb for the B- genome and 781 Mb for the B+ genome. Comparative analysis of the B+ and B- assemblies reveals syntenic discontinuity, duplicated blocks and several insertions, deletions, and inversions indicative of rearrangements in the B+ genome. Hundreds of transposable elements and 1546 protein coding sequences were annotated in the duplicated B+ regions. Our work contributes a list of thousands of genes harbored on the B chromosome, with functions in several biological processes, including the cell cycle.

Temporally consistent species differences in parasite infection but no evidence for rapid parasite-mediated speciation in Lake Victoria cichlid fish.

Parasites may have strong eco-evolutionary interactions with their hosts. Consequently, they may contribute to host diversification. The radiation of cichlid fish in Lake Victoria provides a good model to study the role of parasites in the early stages of speciation. We investigated patterns of macroparasite infection in a community of 17 sympatric cichlids from a recent radiation and 2 older species from 2 nonradiating lineages, to explore the opportunity for parasite-mediated speciation. Host species had different parasite infection profiles, which were only partially explained by ecological factors (diet, water depth). This may indicate that differences in infection are not simply the result of differences in exposure, but that hosts evolved species-specific resistance, consistent with parasite-mediated divergent selection. Infection was similar between sampling years, indicating that the direction of parasite-mediated selection is stable through time. We morphologically identified 6 Cichlidogyrus species, a gill parasite that is considered a good candidate for driving parasite-mediated speciation, because it is host species-specific and has radiated elsewhere in Africa. Species composition of Cichlidogyrus infection was similar among the most closely related host species (members of the Lake Victoria radiation), but two more distantly related species (belonging to nonradiating sister lineages) showed distinct infection profiles. This is inconsistent with a role for Cichlidogyrus in the early stages of divergence. To conclude, we find significant interspecific variation in parasite infection profiles, which is temporally consistent. We found no evidence that Cichlidogyrus-mediated selection contributes to the early stages of speciation. Instead, our findings indicate that species differences in infection accumulate after speciation.

Genome-wide effects of social status on DNA methylation in the brain of a cichlid fish, Astatotilapia burtoni.

BACKGROUND: Successful social behavior requires real-time integration of information about the environment, internal physiology, and past experience. The molecular substrates of this integration are poorly understood, but likely modulate neural plasticity and gene regulation. In the cichlid fish species Astatotilapia burtoni, male social status can shift rapidly depending on the environment, causing fast behavioral modifications and a cascade of changes in gene transcription, the brain, and the reproductive system. These changes can be permanent but are also reversible, implying the involvement of a robust but flexible mechanism that regulates plasticity based on internal and external conditions. One candidate mechanism is DNA methylation, which has been linked to social behavior in many species, including A. burtoni. But, the extent of its effects after A. burtoni social change were previously unknown. RESULTS: We performed the first genome-wide search for DNA methylation patterns associated with social status in the brains of male A. burtoni, identifying hundreds of Differentially Methylated genomic Regions (DMRs) in dominant versus non-dominant fish. Most DMRs were inside genes supporting neural development, synapse function, and other processes relevant to neural plasticity, and DMRs could affect gene expression in multiple ways. DMR genes were more likely to be transcription factors, have a duplicate elsewhere in the genome, have an anti-sense lncRNA, and have more splice variants than other genes. Dozens of genes had multiple DMRs that were often seemingly positioned to regulate specific splice variants. CONCLUSIONS: Our results revealed genome-wide effects of A. burtoni social status on DNA methylation in the brain and strongly suggest a role for methylation in modulating plasticity across multiple biological levels. They also suggest many novel hypotheses to address in mechanistic follow-up studies, and will be a rich resource for identifying the relationships between behavioral, neural, and transcriptional plasticity in the context of social status.

Interspecific hybridization can generate functional novelty in cichlid fish.

The role of interspecific hybridization in evolution is still being debated. Interspecific hybridization has been suggested to facilitate the evolution of ecological novelty, and hence the invasion of new niches and adaptive radiation when ecological opportunity is present beyond the parental species niches. On the other hand, hybrids between two ecologically divergent species may perform less well than parental species in their respective niches because hybrids would be intermediate in performance in both niches. The evolutionary consequences of hybridization may hence be context-dependent, depending on whether ecological opportunities, beyond those of the parental species, do or do not exist. Surprisingly, these complementary predictions may never have been tested in the same experiment in animals. To do so, we investigate if hybrids between ecologically distinct cichlid species perform less well than the parental species when feeding on food either parent is adapted to, and if the same hybrids perform better than their parents when feeding on food none of the species are adapted to. We generated two first-generation hybrid crosses between species of African cichlids. In feeding efficiency experiments we measured the performance of hybrids and parental species on food types representing both parental species niches and additional 'novel' niches, not used by either of the parental species but by other species in the African cichlid radiations. We found that hybrids can have higher feeding efficiencies on the 'novel' food types but typically have lower efficiencies on parental food types when compared to parental species. This suggests that hybridization can generate functional variation that can be of ecological relevance allowing the access to resources outside of either parental species niche. Hence, we provide support for the hypothesis of ecological context-dependency of the evolutionary impact of interspecific hybridization.

Specialized movement and laterality of fin-biting behaviour in Genyochoromis mento in Lake Malawi.

Several vertebrates, including fish, exhibit behavioural laterality and associated morphological asymmetry. Laterality may increase individual fitness as well as foraging strength, accuracy and speed. However, little is known about which behaviours are affected by laterality or what fish species exhibit obvious laterality. Previous research on the predatory behaviour of the scale-eating Lake Tanganyika cichlid Perissodus microlepis indicates behavioural laterality that reflects asymmetric jaw morphology. The Lake Malawi cichlid Genyochromis mento feeds on the fins of other fish, a behaviour that G. mento developed independently from the Tanganyikan Perissodini scale eaters. We investigated stomach contents and behavioural laterality of predation in aquarium to clarify the functional roles and evolution of laterality in cichlids. We also compared the behavioural laterality and mouth asymmetry of G. mento and P. microlepis The diet of G. mento mostly includes fin fragments, but also scales of several fish species. Most individual G. mento specimens showed significant attack bias favouring the skew mouth direction. However, there was no difference in success rate between attacks from the preferred side and those from the non-preferred side, and no lateralized kinetic elements in predation behaviour. Genyochromismento showed weaker laterality than P. microlepis, partly because of their different feeding habits, the phylogenetic constraints from their shorter evolutionary history and their origin from ancestor Haplochromini omnivorous/herbivorous species. Taken together, this study provides new insights into the functional roles of behavioural laterality: predatory fish aiming for prey that show escape behaviours frequently exhibit lateralized behaviour in predation.

Does a cichlid fish process face holistically? Evidence of the face inversion effect.

Faces are the most important body part for differentiating among human individuals by humans. Humans read the face as a whole, rather than looking at its parts, which makes it more difficult to recognise inverted faces than upright. Some other mammals also identify each other based on the upright face and take longer to recognise inverted faces. This effect is called the face inversion effect and is considered as evidence for face-specific perception. This ability has rarely been observed in animals other than mammals, but it was recently reported that some fish species could distinguish among individuals based on the face. For example, the cichlid fish Neolamprologus pulcher rapidly recognises familiar conspecifics by faces rather than other body parts. Here, we examined the face inversion effect in N. pulcher, by showing photographs of conspecific fish faces and objects in both upright and inverted orientations. Subjects gazed at novel faces longer than familiar faces in upright presentation, whereas they did not show such a tendency for inverted faces. Although the object discrimination was difficult, we did not observe the difference between upright and inverted object photographs. Our results indicate that fish exhibits the inversion effect for faces. These findings suggest that N. pulcher may process their conspecifics' face holistically, like humans.

The onset of ecological diversification 50 years after colonization of a crater lake by haplochromine cichlid fishes.

Adaptive radiation research typically relies on the study of evolution in retrospective, leaving the predictive value of the concept hard to evaluate. Several radiations, including the cichlid fishes in the East African Great Lakes, have been studied extensively, yet no study has investigated the onset of the intraspecific processes of niche expansion and differentiation shortly after colonization of an adaptive zone by cichlids. Haplochromine cichlids of one of the two lineages that seeded the Lake Victoria radiation recently arrived in Lake Chala, a lake perfectly suited for within-lake cichlid speciation. Here, we infer the colonization and demographic history, quantify phenotypic, ecological and genomic diversity and diversification, and investigate the selection regime to ask if the population shows signs of diversification resembling the onset of adaptive radiation. We find that since their arrival in the lake, haplochromines have colonized a wide range of depth habitats associated with ecological and morphological expansion and the beginning of phenotypic differentiation and potentially nascent speciation, consistent with the very early onset of an adaptive radiation process. Moreover, we demonstrate evidence of rugged phenotypic fitness surfaces, indicating that current ecological selection may contribute to the phenotypic diversification.

Continuous 1.3-million-year record of East African hydroclimate, and implications for patterns of evolution and biodiversity.

The transport of moisture in the tropics is a critical process for the global energy budget and on geologic timescales, has markedly influenced continental landscapes, migratory pathways, and biological evolution. Here we present a continuous, first-of-its-kind 1.3-My record of continental hydroclimate and lake-level variability derived from drill core data from Lake Malawi, East Africa (9-15° S). Over the Quaternary, we observe dramatic shifts in effective moisture, resulting in large-scale changes in one of the world's largest lakes and most diverse freshwater ecosystems. Results show evidence for 24 lake level drops of more than 200 m during the Late Quaternary, including 15 lowstands when water levels were more than 400 m lower than modern. A dramatic shift is observed at the Mid-Pleistocene Transition (MPT), consistent with far-field climate forcing, which separates vastly different hydroclimate regimes before and after ∼800,000 years ago. Before 800 ka, lake levels were lower, indicating a climate drier than today, and water levels changed frequently. Following the MPT high-amplitude lake level variations dominate the record. From 800 to 100 ka, a deep, often overfilled lake occupied the basin, indicating a wetter climate, but these highstands were interrupted by prolonged intervals of extreme drought. Periods of high lake level are observed during times of high eccentricity. The extreme hydroclimate variability exerted a profound influence on the Lake Malawi endemic cichlid fish species flock; the geographically extensive habitat reconfiguration provided novel ecological opportunities, enabling new populations to differentiate rapidly to distinct species.

Demographic modelling with whole-genome data reveals parallel origin of similar Pundamilia cichlid species after hybridization.

Modes and mechanisms of speciation are best studied in young species pairs. In older taxa, it is increasingly difficult to distinguish what happened during speciation from what happened after speciation. Lake Victoria cichlids in the genus Pundamilia encompass a complex of young species and polymorphic populations. One Pundamilia species pair, P. pundamilia and P. nyererei, is particularly well suited to study speciation because sympatric population pairs occur with different levels of phenotypic differentiation and reproductive isolation at different rocky islands within the lake. Genetic distances between allopatric island populations of the same nominal species often exceed those between the sympatric species. It thus remained unresolved whether speciation into P. nyererei and P. pundamilia occurred once, followed by geographical range expansion and interspecific gene flow in local sympatry, or if the species pair arose repeatedly by parallel speciation. Here, we use genomic data and demographic modelling to test these alternative evolutionary scenarios. We demonstrate that gene flow plays a strong role in shaping the observed patterns of genetic similarity, including both gene flow between sympatric species and gene flow between allopatric populations, as well as recent and early gene flow. The best supported model for the origin of P. pundamilia and P. nyererei population pairs at two different islands is one where speciation happened twice, whereby the second speciation event follows shortly after introgression from an allopatric P. nyererei population that arose earlier. Our findings support the hypothesis that very similar species may arise repeatedly, potentially facilitated by introgressed genetic variation.

One gene, multiple alleles: insights into the microevolution of pigmentation polymorphisms.

Understanding how phenotypic variation is generated and maintained, and the evolutionary forces that shape these processes is the main goal of evolutionary biology. Great progress has been made in uncovering the genetic basis of morphological diversity, yet little is known about both the genetics and developmental basis of discrete polymorphisms segregating in wild populations. Exploring variation in developmental mechanisms at the population level can address the long-standing question of whether the mechanisms of change are the same at the micro- and macroevolutionary scale. This integration has been difficult mainly because the study of the evolution of developmental mechanisms and population genetics remain separate (Genetics, 195, 625 and 2013). In this issue of Molecular Ecology, Roberts et al. (Molecular Ecology and 2017) make a significant contribution towards bridging this gap by studying the genetic and developmental basis of an extremely variable pigmentation pattern. A polymorphic blotched coloration is common among females of four genera of Lake Malawi cichlids. The presence of this phenotype associates with a noncoding SNP upstream of the transcription factor pax7a (Science, 326, 998 and 2009). The authors describe in detail the morphs' pigmentation development, showing that phenotypic differences result from alterations in pigment cell development and survival. Next, using controlled crosses and population genetics studies, they identified three putative pax7a dominant blotch alleles that are associated with specific morphs. These different alleles lead to higher levels of pax7a transcript that correlate with different pigment cell composition. Finally, sequence comparison of the locus within populations and between species revealed a common origin of the allele controlling the blotched morph followed by a pattern of sequential appearance of derived alleles that gave rise to morph diversity. The coupling of the evolutionary history of this allelic series with the developmental analysis of the phenotype paves the way for a mechanistic understanding of morphological innovation and diversification.

An allelic series at pax7a is associated with colour polymorphism diversity in Lake Malawi cichlid fish.

Despite long-standing interest in the evolution and maintenance of discrete phenotypic polymorphisms, the molecular genetic basis of such polymorphism in the wild is largely unknown. Female sex-associated blotched colour polymorphisms found in cichlids of Lake Malawi, East Africa, represent a highly successful polymorphic phenotype, found and maintained in four genera across the geographic expanse of the lake. Previously, we identified an association with an allelic variant of the paired-box transcription factor gene pax7a and blotched colour morphs in Lake Malawi cichlid fishes. Although a diverse range of blotched phenotypes are present in Lake Malawi cichlid species, they all appeared to result from an allele of pax7a that produces increased levels of transcript. Here, we examine the developmental and genetic basis of variation among blotched morphs. First, we confirm that pax7a-associated blotch morphs result primarily from modulation of melanophore development and survival. From laboratory crosses and natural population studies, we identify at least three alleles of pax7a associated with discrete subtypes of blotched morphs, in addition to the ancestral pax7a allele. Genotypes at pax7a support initial evolution of a novel pax7a allele to produce the blotched class of morphs, followed by subsequent evolution of that pax7a blotched allele to produce additional alleles associated with discrete colour morphs. Variant alleles of pax7a produce different levels of pax7a transcript, correlating with pigmentation phenotype at the cellular level. This naturally selected allelic series should serve as a case study for understanding the molecular genetic control of pax7a expression and the evolution of sex-associated alleles.

Towards understanding the genetic basis of mouth asymmetry in the scale-eating cichlid Perissodus microlepis.

How polymorphisms consisting in left-right asymmetries are produced and maintained in natural populations is a tantalizing question, which remains largely unanswered. The scale-eating cichlid fish Perissodus microlepis is a remarkable example of extreme ecological specialization achieved by morphological and behavioural laterality. Its asymmetric mouth is accompanied by a pronounced lateralized foraging behaviour, where a left-bending morph preferentially feeds on the scales of the right side of its prey, while the opposite is true for the right morph. This striking asymmetry made this fish a textbook example of the astounding degree of ecological specialization and negative frequency-dependent selection. Yet, the genetic basis underlying this spectacular laterality remains unknown. We addressed this question through analyses of wild-caught fish using high-throughput DNA sequencing data. A novel array of SNP markers was developed by ddRAD sequencing (ddRADseq) and the use of pooled DNA samples (PoolSeq). We obtained more than 155 000 SNPs using ddRADseq and 3 900 000 SNPs with PoolSeq. Among these, we identified one (ddRAD) SNP, and 38 or 378 (PoolSeq) windows that are differentiated between the left and right morphs accounting for spurious associations due to geographic structuring. This allowed us to uncover candidate genomic regions that potentially contain genes for this trait. Then, this interesting trait has a genetic basis that is likely to be influenced by multiple loci. This result contributes to a greater understanding of the genetic bases of left-right asymmetry and, ultimately, the evolutionary processes governing the maintenance of this striking case of laterality.

Molecular characterization of MHC class IIB genes of sympatric Neotropical cichlids.

BACKGROUND: The Major Histocompatibility Complex (MHC) is a key component of the adaptive immune system of all vertebrates and consists of the most polymorphic genes known to date. Due to this complexity, however, MHC remains to be characterized in many species including any Neotropical cichlid fish. Neotropical crater lake cichlids are ideal models to study evolutionary processes as they display one of the most convincing examples of sympatric and repeated parallel radiation events within and among isolated crater lakes. RESULTS: Here, we characterized the genes of MHC class IIB chain of the Midas cichlid species complex (Amphilophus cf. citrinellus) including fish from five lakes in Nicaragua. We designed 19 new specific primers anchored in a stepwise fashion in order to detect all alleles present. We obtained 866 genomic DNA (gDNA) sequences from thirteen individuals and 756 additional sequences from complementary DNA (cDNA) of seven of those individuals. We identified 69 distinct alleles with up to 25 alleles per individual. We also found considerable intron length variation and mismatches of alleles detected in cDNA and gDNA suggesting that some loci have undergone pseudogenization. Lastly, we created a model of protein structure homology for each allele and identified their key structural components. CONCLUSIONS: Overall, the Midas cichlid has one of the most diverse repertoires of MHC class IIB genes known, which could serve as a powerful tool to elucidate the process of divergent radiations, colonization and speciation in sympatry.