Genetic basis of male colour dimorphism in a Lake Tanganyika cichlid fish

Stages of embryonic development in the live-bearing fish, Gambusia holbrooki

BACKGROUND

Divergent morphology and placentation of Poeciliids make them suitable model for investigating how evolutionary selection has altered and conserved the developmental mechanisms. However, there is limited description of their embryonic staging, despite representing a key evolutionary node that shares developmental strategy with placental vertebrates. Here, we describe the embryonic developmental stages of Gambusia holbrooki from zygote to parturition using freshly harvested embryos.

RESULTS

We defined 40 embryonic stages using a numbered (stages 0-39; zygote to parturition, respectively) and named (grouped into seven periods, ie, zygote, cleavage, blastula, gastrula, segmentation, pharyngula, and parturition) staging system. Two sets of quantitative (ie, egg diameter, embryonic total length, otic vesicle closure index, heart rates, the number of caudal fin rays and elements) and qualitative (ie, three-dimensional analysis of images and key morphological criteria) data were acquired and used in combination to describe each stage. All 40 stages are separated by well-defined morphological traits, revealing developmental novelties that are influenced by narrow perivitelline space, placentation, internal gestation, and sex differentiation.

CONCLUSIONS

The principal diagnostic features described are quick, reliable, and easy to apply. This system will benefit researchers investigating molecular ontogeny, particularly sexual differentiation mechanisms in G. holbrooki.

Novel Sex Chromosomes in 3 Cichlid Fishes from Lake Tanganyika

African cichlids are well known for their adaptive radiations, but it is now apparent that they also harbor an extraordinary diversity of sex chromosome systems. In this study, we sequenced pools of males and females from species in 3 different genera of cichlids from Lake Tanganyika. We then searched for regions that were differentiated following the patterns expected for sex chromosomes. We report 2 novel sex chromosomes systems, an XY system on LG19 in Tropheus sp. "black" and a ZW system on LG7 in Hemibates stenosoma. We also identify a ZW system on LG5 in Cyprichromis leptosoma that may be convergent with a system previously described in Lake Malawi cichlids. Our data also identify candidate single nucleotide polymorphisms for the blue/yellow tail color polymorphism observed among male C. leptosoma.

Population genomics of rapid evolution in natural populations: polygenic selection in response to power station thermal effluents

BACKGROUND

Examples of rapid evolution are common in nature but difficult to account for with the standard population genetic model of adaptation. Instead, selection from the standing genetic variation permits rapid adaptation via soft sweeps or polygenic adaptation. Empirical evidence of this process in nature is currently limited but accumulating.

RESULTS

We provide genome-wide analyses of rapid evolution in Fundulus heteroclitus populations subjected to recently elevated temperatures due to coastal power station thermal effluents using 5449 SNPs across two effluent-affected and four reference populations. Bayesian and multivariate analyses of population genomic structure reveal a substantial portion of genetic variation that is most parsimoniously explained by selection at the site of thermal effluents. An FST outlier approach in conjunction with additional conservative requirements identify significant allele frequency differentiation that exceeds neutral expectations among exposed and closely related reference populations. Genomic variation patterns near these candidate loci reveal that individuals living near thermal effluents have rapidly evolved from the standing genetic variation through small allele frequency changes at many loci in a pattern consistent with polygenic selection on the standing genetic variation.

CONCLUSIONS

While the ultimate trajectory of selection in these populations is unknown and we survey only a minority of genomic loci, our findings suggest that polygenic models of adaptation may play important roles in large, natural populations experiencing recent selection due to environmental changes that cause broad physiological impacts.

Morphometric and genetic differentiation among populations of flat-headed cusimanse (Crossarchus platycephalus) in Nigeria

Geographic barriers can partition genetic diversity among populations and drive evolutionary divergence between populations, promoting the speciation process and affecting conservation goals. We integrated morphological and genomic data to assess the distribution of variation in the flat-headed cusimanse (Crossarchus platycephalus), a species of least conservation concern, on either side of the River Niger in Nigeria. Ecological disturbances affect the conservation status of many other animals in this region. The two populations were differentiated in the snout and fore limbs, with greater morphological diversity in the western population. We used Restriction site Associated DNA sequencing (RAD-seq) and identified two genotypic clusters in a STRUCTURE analysis. Individuals from the eastern population are almost entirely assigned to one cluster, whereas genotypes from the western population are a mixture of the two clusters. The population from west of the River Niger also had higher heterozygosity. The morphological and population genetic data are therefore in agreement that the population from west of the River Niger is more diverse than the eastern population, and the eastern population contains a subset of the genetic variation found in the western population. Our results demonstrate that combining morphological and genotypic measures of diversity can provide a congruent picture of the distribution of intraspecific variation. The results also suggest that future work should explore the role of the River Niger as a natural barrier to migration in Nigeria.

The Colorful Sex Chromosomes of Teleost Fish

Teleost fish provide some of the most intriguing examples of sexually dimorphic coloration, which is often advantageous for only one of the sexes. Mapping studies demonstrated that the genetic loci underlying such color patterns are frequently in tight linkage to the sex-determining locus of a species, ensuring sex-specific expression of the corresponding trait. Several genes affecting color synthesis and pigment cell development have been previously described, but the color loci on the sex chromosomes have mostly remained elusive as yet. Here, we summarize the current knowledge about the genetics of such color loci in teleosts, mainly from studies on poeciliids and cichlids. Further studies on these color loci will certainly provide important insights into the evolution of sex chromosomes.

SNP Discovery from Single and Multiplex Genome Assemblies of Non-model Organisms

Population genetic studies of non-model organisms often rely on initial ascertainment of genetic markers from a single individual or a small pool of individuals. This initial screening has been a significant barrier to beginning population studies on non-model organisms (Aitken et al., Mol Ecol 13:1423-1431, 2004; Morin et al., Trends Ecol Evol 19:208-216, 2004). As genomic data become increasingly available for non-model species, SNP ascertainment from across the genome can be performed directly from published genome contigs and short-read archive data. Alternatively, low to medium genome coverage from shotgun NGS library sequencing of single or pooled samples, or from reduced-representation libraries (e.g., capture enrichment; see Ref. "Hancock-Hanser et al., Mol Ecol Resour 13:254-268, 2013") can produce sufficient new data for SNP discovery with limited investment. We describe protocols for assembly of short read data to reference or related species genome contig sequences, followed by SNP discovery and filtering to obtain an optimal set of SNPs for population genotyping using a variety of downstream high-throughput genotyping methods.

Using NextRAD sequencing to infer movement of herbivores among host plants

Herbivores often move among spatially interspersed host plants, tracking high-quality resources through space and time. This dispersal is of particular interest for vectors of plant pathogens. Existing molecular tools to track such movement have yielded important insights, but often provide insufficient genetic resolution to infer spread at finer spatiotemporal scales. Here, we explore the use of Nextera-tagmented reductively-amplified DNA (NextRAD) sequencing to infer movement of a highly-mobile winged insect, the potato psyllid (Bactericera cockerelli), among host plants. The psyllid vectors the pathogen that causes zebra chip disease in potato (Solanum tuberosum), but understanding and managing the spread of this pathogen is limited by uncertainty about the insect's host plant(s) outside of the growing season. We identified 1,978 polymorphic loci among psyllids separated spatiotemporally on potato or in patches of bittersweet nightshade (S. dulcumara), a weedy plant proposed to be the source of potato-colonizing psyllids. A subset of the psyllids on potato exhibited genetic similarity to insects on nightshade, consistent with regular movement between these two host plants. However, a second subset of potato-collected psyllids was genetically distinct from those collected on bittersweet nightshade; this suggests that a currently unrecognized source, i.e., other nightshade patches or a third host-plant species, could be contributing to psyllid populations in potato. Oftentimes, dispersal of vectors of pathogens must be tracked at a fine scale in order to understand, predict, and manage disease spread. We demonstrate that emerging sequencing technologies that detect genome-wide SNPs of a vector can be used to infer such localized movement.

Integration of Next Generation Sequencing and EPR Analysis to Uncover Molecular Mechanism Underlying Shell Color Variation in Scallops

The Yesso scallop Patinopecten yessoensis displays polymorphism in shell colors, which is of great interest for the scallop industry. To identify genes involved in the shell coloration, in the present study, we investigate the transcriptome differences by Illumina digital gene expression (DGE) analysis in two extreme color phenotypes, Red and White. Illumina sequencing yields a total of 62,715,364 clean sequence reads, and more than 85% reads are mapped into our previously sequenced transcriptome. There are 25 significantly differentially expressed genes between Red and White scallops. EPR (Electron paramagnetic resonance) analysis has identified EPR spectra of pheomelanin and eumelanin in the red shells, but not in the white shells. Compared to the Red scallops, the White scallops have relatively higher mRNA expression in tyrosinase genes, but lower expression in other melanogensis-associated genes. Meantime, the relatively lower tyrosinase protein and decreased tyrosinase activity in White scallops are suggested to be associated with the lack of melanin in the white shells. Our findings highlight the functional roles of melanogensis-associated genes in the melanization process of scallop shells, and shed new lights on the transcriptional and post-transcriptional mechanisms in the regulation of tyrosinase activity during the process of melanin synthesis. The present results will assist our molecular understanding of melanin synthesis underlying shell color polymorphism in scallops, as well as other bivalves, and also help the color-based breeding in shellfish aquaculture.

Extrapolating from local ecological processes to genus-wide patterns in colour polymorphism in South African Protea

Polymorphic traits are central to many fundamental discoveries in evolution, yet why they are found in some species and not others remains poorly understood. We use the African genus Protea-within which more than 40% of species have co-occurring pink and white floral colour morphs-to ask whether convergent evolution and ecological similarity could explain the genus-wide pattern of polymorphism. First, we identified environmental correlates of pink morph frequency across 28 populations of four species. Second, we determined whether the same correlates could predict species-level polymorphism and monomorphism across 31 species. We found that pink morph frequency increased with elevation in Protea repens and three section Exsertae species, increased eastward in P. repens, and increased with seed predation intensity in section Exsertae. For cross-species comparisons, populations of monomorphic pink species occurred at higher elevations than populations of monomorphic white species, and 18 polymorphic species spanned broader elevational gradients than 13 monomorphic species. These results suggest that divergent selection along elevational clines has repeatedly favoured polymorphism, and that more uniform selection in altitudinally restricted species may promote colour monomorphism. Our findings are, to our knowledge, the first to link selection acting within species to the presence and absence of colour polymorphism at broader phylogenetic scales.

Quantitative Genetic Analyses of Male Color Pattern and Female Mate Choice in a Pair of Cichlid Fishes of Lake Malawi, East Africa

The traits involved in sexual selection, such as male secondary sexual characteristics and female mate choice, often co-evolve which can promote population differentiation. However, the genetic architecture of these phenotypes can influence their evolvability and thereby affect the divergence of species. The extraordinary diversity of East African cichlid fishes is often attributed to strong sexual selection and thus this system provides an excellent model to test predictions regarding the genetic architecture of sexually selected traits that contribute to reproductive isolation. In particular, theory predicts that rapid speciation is facilitated when male sexual traits and female mating preferences are controlled by a limited number of linked genes. However, few studies have examined the genetic basis of male secondary sexual traits and female mating preferences in cichlids and none have investigated the genetic architecture of both jointly. In this study, we artificially hybridized a pair of behaviorally isolated cichlid fishes from Lake Malawi and quantified both melanistic color pattern and female mate choice. We investigated the genetic architecture of both phenotypes using quantitative genetic analyses. Our results suggest that 1) many non-additively acting genetic factors influence melanistic color patterns, 2) female mate choice may be controlled by a minimum of 1-2 non-additive genetic factors, and 3) F2 female mate choice is not influenced by male courting effort. Furthermore, a joint analysis of color pattern and female mate choice indicates that the genes underlying these two traits are unlikely to be physically linked. These results suggest that reproductive isolation may evolve rapidly owing to the few genetic factors underlying female mate choice. Hence, female mate choice likely played an important role in the unparalleled speciation of East African cichlid fish.

Carotenoid-based coloration in cichlid fishes

Animal colors play important roles in communication, ecological interactions and speciation. Carotenoid pigments are responsible for many yellow, orange and red hues in animals. Whereas extensive knowledge on the proximate mechanisms underlying carotenoid coloration in birds has led to testable hypotheses on avian color evolution and signaling, much less is known about the expression of carotenoid coloration in fishes. Here, we promote cichlid fishes (Perciformes: Cichlidae) as a system in which to study the physiological and evolutionary significance of carotenoids. Cichlids include some of the best examples of adaptive radiation and color pattern diversification in vertebrates. In this paper, we examine fitness correlates of carotenoid pigmentation in cichlids and review hypotheses regarding the signal content of carotenoid-based ornaments. Carotenoid-based coloration is influenced by diet and body condition and is positively related to mating success and social dominance. Gaps in our knowledge are discussed in the last part of this review, particularly in the understanding of carotenoid metabolism pathways and the genetics of carotenoid coloration. We suggest that carotenoid metabolism and transport are important proximate mechanisms responsible for individual and population-differences in cichlid coloration that may ultimately contribute to diversification and speciation.

Frequency-dependent and correlational selection pressures have conflicting consequences for assortative mating in a color-polymorphic lizard, Uta stansburiana

Genetically determined polymorphisms incorporating multiple traits can persist in nature under chronic, fluctuating, and sometimes conflicting selection pressures. Balancing selection among morphs preserves equilibrium frequencies, while correlational selection maintains favorable trait combinations within each morph. Under negative frequency-dependent selection, females should mate (often disassortatively) with rare male morphotypes to produce conditionally fit offspring. Conversely, under correlational selection, females should mate assortatively to preserve coadapted gene complexes and avoid ontogenetic conflict. Using controlled breeding designs, we evaluated consequences of assortative mating patterns in color-polymorphic side-blotched lizards (Uta stansburiana), to identify conflict between these sources of selection. Females who mated disassortatively, and to conditionally high-quality males in the context of frequency-dependent selection, experienced highest fertility rates. In contrast, assortatively mated females experienced higher fetal viability rates. The trade-off between fertility and egg viability resulted in no overall fitness benefit to either assortative or disassortative mating patterns. These results suggest that ongoing conflict between correlational and frequency dependent selection in polymorphic populations may generate a trade-off between rare-morph advantage and phenotypic integration and between assortative and disassortative mating decisions. More generally, interactions among multiple sources of diversity-promoting selection can alter adaptations and dynamics predicted to arise under any of these regimes alone.

Colour variation in cichlid fish: developmental mechanisms, selective pressures and evolutionary consequences

Cichlid fishes constitute one of the most species-rich families of vertebrates. In addition to complex social behaviour and morphological versatility, they are characterised by extensive diversity in colouration, both within and between species. Here, we review the cellular and molecular mechanisms underlying colour variation in this group and the selective pressures responsible for the observed variation. We specifically address the evidence for the hypothesis that divergence in colouration is associated with the evolution of reproductive isolation between lineages. While we conclude that cichlid colours are excellent models for understanding the role of animal communication in species divergence, we also identify taxonomic and methodological biases in the current research effort. We suggest that the integration of genomic approaches with ecological and behavioural studies, across the entire cichlid family and beyond it, will contribute to the utility of the cichlid model system for understanding the evolution of biological diversity.

Mapping of pigmentation QTL on an anchored genome assembly of the cichlid fish, Metriaclima zebra

Background

Pigmentation patterns are one of the most recognizable phenotypes across the animal kingdom. They play an important role in camouflage, communication, mate recognition and mate choice. Most progress on understanding the genetics of pigmentation has been achieved via mutational analysis, with relatively little work done to understand variation in natural populations. Pigment patterns vary dramatically among species of cichlid fish from Lake Malawi, and are thought to be important in speciation. In this study, we crossed two species, Metriaclima zebra and M. mbenjii, that differ in several aspects of their body and fin color. We genotyped 798 SNPs in 160 F₂ male individuals to construct a linkage map that was used to identify quantitative trait loci (QTL) associated with the pigmentation traits of interest. We also used the linkage map to anchor portions of the M. zebra genome assembly.

Results

We constructed a linkage map consisting of 834 markers in 22 linkage groups that spanned over 1,933 cM. QTL analysis detected one QTL each for dorsal fin xanthophores, caudal fin xanthophores, and pelvic fin melanophores. Dorsal fin and caudal fin xanthophores share a QTL on LG12, while pelvic fin melanophores have a QTL on LG11. We used the mapped markers to anchor 66.5% of the M. zebra genome assembly. Within each QTL interval we identified several candidate genes that might play a role in pigment cell development.

Conclusion

This is one of a few studies to identify QTL for natural variation in fish pigmentation. The QTL intervals we identified did not contain any pigmentation genes previously identified by mutagenesis studies in other species. We expect that further work on these intervals will identify new genes involved in pigment cell development in natural populations.

RAD-Seq derived markers flank the shell colour and banding loci of the Cepaea nemoralis supergene

Studies on the classic shell colour and banding polymorphism of the land snail Cepaea played a crucial role in establishing the importance of natural selection in maintaining morphological variation. Cepaea is also a pre-eminent model for ecological genetics because the outward colour and banding phenotype is entirely genetically determined, primarily by a 'supergene' of at least five loci. Unfortunately, progress in understanding the evolution and maintenance of the Cepaea polymorphism stalled, partly because of a lack of genetic markers. With a view to re-establish Cepaea as a prominent model of molecular ecology, we made six laboratory crosses of Cepaea nemoralis, five of which segregated for shell ground colour (C) and the presence or absence of bands (B). First, scoring of colour and banding in 323 individuals found no recombination between the C and B loci of the supergene. Second, using restriction site-associated DNA sequencing (RAD-Seq) of two parents and 22 offspring, we identified 44 anonymous markers putatively linked to the colour (C) and banding (B) loci. The genotype of eleven of the most promising RAD-Seq markers was independently validated in the same 22 offspring, then up to a further 146 offspring were genotyped. The closest RAD-Seq markers scored are within ~0.6 centimorgan (cM) of the C-B supergene linkage group, with the combined loci together forming a 35.8 cM linkage map of markers that flank both sides of the Cepaea C-B supergene.