Distinct genetic origins of eumelanin intensity and barring patterns in cichlid fishes.

Pigment patterns are incredibly diverse across vertebrates and are shaped by multiple selective pressures from predator avoidance to mate choice. A common pattern across fishes, but for which we know little about the underlying mechanisms, is repeated melanic vertical bars. In order to understand genetic factors that modify the level or pattern of vertical barring, we generated a genetic cross of 322 F2 hybrids between two cichlid species with distinct barring patterns, Aulonocara koningsi and Metriaclima mbenjii. We identify 48 significant quantitative trait loci that underlie a series of seven phenotypes related to the relative pigmentation intensity, and four traits related to patterning of the vertical bars. We find that genomic regions that generate variation in the level of eumelanin produced are largely independent of those that control the spacing of vertical bars. Candidate genes within these intervals include novel genes and those newly-associated with vertical bars, which could affect melanophore survival, fate decisions, pigment biosynthesis, and pigment distribution. Together, this work provides insights into the regulation of pigment diversity, with direct implications for an animal's fitness and the speciation process.

Genetic basis of ecologically relevant body shape variation among four genera of cichlid fishes.

Divergence in body shape is one of the most widespread and repeated patterns of morphological variation in fishes and is associated with habitat specification and swimming mechanics. Such ecological diversification is the first stage of the explosive adaptive radiation of cichlid fishes in the East African Rift Lakes. We use two hybrid crosses of cichlids (Metriaclima sp. × Aulonocara sp. and Labidochromis sp. × Labeotropheus sp., >975 animals total) to determine the genetic basis of body shape diversification that is similar to benthic-pelagic divergence across fishes. Using a series of both linear and geometric shape measurements, we identified 34 quantitative trait loci (QTL) that underlie various aspects of body shape variation. These QTL are spread throughout the genome, each explaining 3.2-8.6% of phenotypic variation, and are largely modular. Further, QTL are distinct both between these two crosses of Lake Malawi cichlids and compared to previously identified QTL for body shape in fishes such as sticklebacks. We find that body shape is controlled by many genes of small effect. In all, we find that convergent body shape phenotypes commonly observed across fish clades are most likely due to distinct genetic and molecular mechanisms.

Morphometric and Genetic Description of Trophic Adaptations in Cichlid Fishes.

Since Darwin, biologists have sought to understand the evolution and origins of phenotypic adaptations. The skull is particularly diverse due to intense natural selection on feeding biomechanics. We investigated the genetic and molecular origins of trophic adaptation using Lake Malawi cichlids, which have undergone an exemplary evolutionary radiation. We analyzed morphological differences in the lateral and ventral head shape among an insectivore that eats by suction feeding, an obligate biting herbivore, and their F2 hybrids. We identified variation in a series of morphological traits-including mandible width, mandible length, and buccal length-that directly affect feeding kinematics and function. Using quantitative trait loci (QTL) mapping, we found that many genes of small effects influence these craniofacial adaptations. Intervals for some traits were enriched in genes related to potassium transport and sensory systems, the latter suggesting co-evolution of feeding structures and sensory adaptations for foraging. Despite these indications of co-evolution of structures, morphological traits did not show covariation. Furthermore, phenotypes largely mapped to distinct genetic intervals, suggesting that a common genetic basis does not generate coordinated changes in shape. Together, these suggest that craniofacial traits are mostly inherited as separate modules, which confers a high potential for the evolution of morphological diversity. Though these traits are not restricted by genetic pleiotropy, functional demands of feeding and sensory structures likely introduce constraints on variation. In all, we provide insights into the quantitative genetic basis of trophic adaptation, identify mechanisms that influence the direction of morphological evolution, and provide molecular inroads to craniofacial variation.

Polygenic sex determination produces modular sex polymorphism in an African cichlid fish.

For many vertebrates, a single genetic locus initiates a cascade of developmental sex differences in the gonad and throughout the organism, resulting in adults with two phenotypically distinct sexes. Species with polygenic sex determination (PSD) have multiple interacting sex determination alleles segregating within a single species, allowing for more than two genotypic sexes and scenarios where sex genotype at a given locus can be decoupled from gonadal sex. Here we investigate the effects of PSD on secondary sexual characteristics in the cichlid fish Metriaclima mbenjii, where one female (W) and one male (Y) sex determination allele interact to produce siblings with four possible sex classes: ZZXX females, ZWXX females, ZWXY females, and ZZXY males. We find that PSD in M. mbenjii produces an interplay of sex linkage and sex limitation resulting in modular variation in morphological and behavioral traits. Further, the evolution or introgression of a newly acquired sex determiner creates additional axes of phenotypic variation for varied traits, including genital morphology, craniofacial morphology, gastrointestinal morphology, and home tank behaviors. In contrast to single-locus sex determination, which broadly results in sexual dimorphism, polygenic sex determination can induce higher-order sexual polymorphism. The modularity of secondary sexual characteristics produced by PSD provides context for understanding the evolutionary causes and consequences of maintenance, gain, or loss of sex determination alleles in populations.

Origin of a Giant Sex Chromosome.

Chromosome size and morphology vary within and among species, but little is known about the proximate or ultimate causes of these differences. Cichlid fish species in the tribe Oreochromini share an unusual giant chromosome that is ∼3 times longer than the other chromosomes. This giant chromosome functions as a sex chromosome in some of these species. We test two hypotheses of how this giant sex chromosome may have evolved. The first hypothesis proposes that it evolved by accumulating repetitive elements as recombination was reduced around a dominant sex determination locus, as suggested by canonical models of sex chromosome evolution. An alternative hypothesis is that the giant sex chromosome originated via the fusion of an autosome with a highly repetitive B chromosome, one of which carried a sex determination locus. We test these hypotheses using comparative analysis of chromosome-scale cichlid and teleost genomes. We find that the giant sex chromosome consists of three distinct regions based on patterns of recombination, gene and transposable element content, and synteny to the ancestral autosome. The WZ sex determination locus encompasses the last ∼105 Mb of the 134-Mb giant chromosome. The last 47 Mb of the giant chromosome shares no obvious homology to any ancestral chromosome. Comparisons across 69 teleost genomes reveal that the giant sex chromosome contains unparalleled amounts of endogenous retroviral elements, immunoglobulin genes, and long noncoding RNAs. The results favor the B chromosome fusion hypothesis for the origin of the giant chromosome.

Chromosome-scale assemblies reveal the structural evolution of African cichlid genomes.

BACKGROUND: African cichlid fishes are well known for their rapid radiations and are a model system for studying evolutionary processes. Here we compare multiple, high-quality, chromosome-scale genome assemblies to elucidate the genetic mechanisms underlying cichlid diversification and study how genome structure evolves in rapidly radiating lineages. RESULTS: We re-anchored our recent assembly of the Nile tilapia (Oreochromis niloticus) genome using a new high-density genetic map. We also developed a new de novo genome assembly of the Lake Malawi cichlid, Metriaclima zebra, using high-coverage Pacific Biosciences sequencing, and anchored contigs to linkage groups (LGs) using 4 different genetic maps. These new anchored assemblies allow the first chromosome-scale comparisons of African cichlid genomes. Large intra-chromosomal structural differences (∼2-28 megabase pairs) among species are common, while inter-chromosomal differences are rare (<10 megabase pairs total). Placement of the centromeres within the chromosome-scale assemblies identifies large structural differences that explain many of the karyotype differences among species. Structural differences are also associated with unique patterns of recombination on sex chromosomes. Structural differences on LG9, LG11, and LG20 are associated with reduced recombination, indicative of inversions between the rock- and sand-dwelling clades of Lake Malawi cichlids. M. zebra has a larger number of recent transposable element insertions compared with O. niloticus, suggesting that several transposable element families have a higher rate of insertion in the haplochromine cichlid lineage. CONCLUSION: This study identifies novel structural variation among East African cichlid genomes and provides a new set of genomic resources to support research on the mechanisms driving cichlid adaptation and speciation.

The intestinal environment as an evolutionary adaptation to mouthbrooding in the Astatotilapia burtoni cichlid.

Many of the various parental care strategies displayed by animals are accompanied by a significant reduction in food intake that imposes a substantial energy trade-off. Mouthbrooding, as seen in several species of fish in which the parent holds the developing eggs and fry in the buccal cavity, represents an extreme example of reduced food intake during parental investment and is accompanied by a range of physiological adaptations. In this study we use 16S sequencing to characterize the gut microbiota of female Astatotilapia burtoni cichlid fish throughout the obligatory phase of self-induced starvation during the brooding cycle in comparison to stage-matched females that have been denied food for the same duration. In addition to a reduction of gut epithelial turnover, we find a dramatic reduction in species diversity in brooding stages that recovers upon release of fry and refeeding that is not seen in females that are simply starved. Based on overall species diversity as well as differential abundance of specific bacterial taxa, we suggest that rather than reflecting a simple deprivation of caloric intake, the gut microbiota is more strongly influenced by physiological changes specific to mouthbrooding including the reduced epithelial turnover and possible production of antimicrobial agents.

Genetic sex determination in Astatotilapia calliptera, a prototype species for the Lake Malawi cichlid radiation.

East African cichlids display extensive variation in sex determination systems. The species Astatotilapia calliptera is one of the few cichlids that reside both in Lake Malawi and in surrounding waterways. A. calliptera is of interest in evolutionary studies as a putative immediate outgroup species for the Lake Malawi species flock and possibly as a prototype ancestor-like species for the radiation. Here, we use linkage mapping to test association of sex in A. calliptera with loci that have been previously associated with genetic sex determination in East African cichlid species. We identify a male heterogametic XY system segregating at linkage group (LG) 7 in an A. calliptera line that originated from Lake Malawi, at a locus previously shown to act as an XY sex determination system in multiple species of Lake Malawi cichlids. Significant association of genetic markers and sex produce a broad genetic interval of approximately 26 megabases (Mb) using the Nile tilapia genome to orient markers; however, we note that the marker with the strongest association with sex is near a gene that acts as a master sex determiner in other fish species. We demonstrate that alleles of the marker are perfectly associated with sex in Metriaclima mbenjii, a species from the rock-dwelling clade of Lake Malawi. While we do not rule out the possibility of other sex determination loci in A. calliptera, this study provides a foundation for fine mapping of the cichlid sex determination gene on LG7 and evolutionary context regarding the origin and persistence of the LG7 XY across diverse, rapidly evolving lineages.

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.

Polygenic sex determination in the cichlid fish Astatotilapia burtoni.

BACKGROUND: The East African riverine cichlid species Astatotilapia burtoni serves as an important laboratory model for sexually dimorphic physiology and behavior, and also serves as an outgroup species for the explosive adaptive radiations of cichlid species in Lake Malawi and Lake Victoria. An astounding diversity of genetic sex determination systems have been revealed within the adaptive radiation of East African cichlids thus far, including polygenic sex determination systems involving the epistatic interaction of multiple, independently segregating sex determination alleles. However, sex determination has remained unmapped in A. burtoni. Here we present mapping results supporting the presence of multiple, novel sex determination alleles, and thus the presence of polygenic sex determination in A. burtoni. RESULTS: Using mapping in small families in conjunction with restriction-site associated DNA sequencing strategies, we identify associations with sex at loci on linkage group 13 and linkage group 5-14. Inheritance patterns support an XY sex determination system on linkage group 5-14 (a chromosome fusion relative to other cichlids studied), and an XYW system on linkage group 13, and these associations are replicated in multiple families. Additionally, combining our genetic data with comparative genomic analysis identifies another fusion that is unassociated with sex, with linkage group 8-24 and linkage group 16-21 fused in A. burtoni relative to other East African cichlid species. CONCLUSIONS: We identify genetic signals supporting the presence of three previously unidentified sex determination alleles at two loci in the species A. burtoni, strongly supporting the presence of polygenic sex determination system in the species. These results provide a foundation for future mapping of multiple sex determination genes and their interactions. A better understanding of sex determination in A. burtoni provides important context for their use in behavioral studies, as well as studies of the evolution of genetic sex determination and sexual conflicts in East African cichlids.

Structure and decay of a proto-Y region in Tilapia, Oreochromis niloticus.

BACKGROUND: Sex-determination genes drive the evolution of adjacent chromosomal regions. Sexually antagonistic selection favors the accumulation of inversions that reduce recombination in regions adjacent to the sex-determination gene. Once established, the clonal inheritance of sex-linked inversions leads to the accumulation of deleterious alleles, repetitive elements and a gradual decay of sex-linked genes. This in turn creates selective pressures for the evolution of mechanisms that compensate for the unequal dosage of gene expression. Here we use whole genome sequencing to characterize the structure of a young sex chromosome and quantify sex-specific gene expression in the developing gonad. RESULTS: We found an 8.8 Mb block of strong differentiation between males and females that corresponds to the location of a previously mapped sex-determiner on linkage group 1 of Oreochromis niloticus. Putatively disruptive mutations are found in many of the genes within this region. We also found a significant female-bias in the expression of genes within the block of differentiation compared to those outside the block of differentiation. Eight candidate sex-determination genes were identified within this region. CONCLUSIONS: This study demonstrates a block of differentiation on linkage group 1, suggestive of an 8.8 Mb inversion encompassing the sex-determining locus. The enrichment of female-biased gene expression inside the proposed inversion suggests incomplete dosage compensation. This study helps establish a model for studying the early-to-intermediate stages of sex chromosome evolution.

Genetic basis of continuous variation in the levels and modular inheritance of pigmentation in cichlid fishes.

Variation in pigmentation type and levels is a hallmark of myriad evolutionary radiations, and biologists have long been fascinated by the factors that promote and maintain variation in coloration across populations. Here, we provide insights into the genetic basis of complex and continuous patterns of colour variation in cichlid fishes, which offer a vast diversity of pigmentation patterns that have evolved in response to both natural and sexual selection. Specifically, we crossed two divergent cichlid species to generate an F2 mapping population that exhibited extensive variation in pigmentation levels and patterns. Our experimental design is robust in that it combines traditional quantitative trait locus (QTL) analysis with population genomics, which has allowed us to move efficiently from QTL interval to candidate gene. In total, we detected 41 QTL and 13 epistatic interactions that underlie melanocyte- and xanthophore-based coloration across the fins and flanks of these fishes. We also identified 2 QTL and 1 interaction for variation in the magnitude of integration among these colour traits. This finding in particular is notable as there are marked differences both within and between species with respect to the complexity of pigmentation patterns. While certain individuals are characterized by more uniform 'integrated' colour patterns, others exhibit many more degrees of freedom with respect to the distribution of colour 'modules' across the fins and flank. Our data reveal, for the first time, a genetic basis for this difference. Finally, we implicate pax3a as a mediator of continuous variation in the levels of xanthophore-based colour along the cichlid flank.

Origins of shared genetic variation in African cichlids.

Cichlid fishes have evolved tremendous morphological and behavioral diversity in the waters of East Africa. Within each of the Great Lakes Tanganyika, Malawi, and Victoria, the phenomena of hybridization and retention of ancestral polymorphism explain allele sharing across species. Here, we explore the sharing of single nucleotide polymorphisms (SNPs) between the major East African cichlid assemblages. A set of approximately 200 genic and nongenic SNPs was ascertained in five Lake Malawi species and genotyped in a diverse collection of ~160 species from across Africa. We observed segregating polymorphism outside of the Malawi lineage for more than 50% of these loci; this holds similarly for genic versus nongenic SNPs, as well as for SNPs at putative CpG versus non-CpG sites. Bayesian and principal component analyses of genetic structure in the data demonstrate that the Lake Malawi endemic flock is not monophyletic and that river species have likely contributed significantly to Malawi genomes. Coalescent simulations support the hypothesis that river cichlids have transported polymorphism between lake assemblages. We observed strong genetic differentiation between Malawi lineages for approximately 8% of loci, with contributions from both genic and nongenic SNPs. Notably, more than half of these outlier loci between Malawi groups are polymorphic outside of the lake. Cichlid fishes have evolved diversity in Lake Malawi as new mutations combined with standing genetic variation shared across East Africa.

A small number of genes underlie male pigmentation traits in Lake Malawi cichlid fishes.

Pigmentation patterns are one of the most recognizable forms of phenotypic diversity and an important component of organismal fitness. While much progress has been made in understanding the genes controlling pigmentation in model systems, many questions remain about the genetic basis of pigment traits observed in nature. Lake Malawi cichlid fishes are known for their diversity of male pigmentation patterns, which have been shaped by sexual selection. To begin the process of identifying the genes underlying this diversity, we quantified the number of pigment cells on the body and fins of two species of the genus Metriaclima and their hybrids. We then used the Castle-Wright equation to estimate that differences in individual pigmentation traits between these species are controlled by one to four genes each. Different pigmentation traits are highly correlated in the F(2) , suggesting shared developmental pathways and genetic pleiotropy. Melanophore and xanthophore traits fall on opposite ends of the first principal component axis of the F(2) phenotypes, suggesting a tradeoff during the development of these two pigment cell types.

Craniofacial divergence and ongoing adaptation via the hedgehog pathway.

Adaptive variation in craniofacial structure contributes to resource specialization and speciation, but the genetic loci that underlie craniofacial adaptation remain unknown. Here we show that alleles of the hedgehog pathway receptor Patched1 (Ptch1) gene are responsible for adaptive variation in the shape of the lower jaw both within and among genera of Lake Malawi cichlid fish. The evolutionarily derived allele of Ptch1 reduces the length of the retroarticular (RA) process of the lower jaw, a change predicted to increase speed of jaw rotation for improved suction-feeding. The alternate allele is associated with a longer RA and a more robustly mineralized jaw, typical of species that use a biting mode of feeding. Genera with the most divergent feeding morphologies are nearly fixed for different Ptch1 alleles, whereas species with intermediate morphologies still segregate variation at Ptch1. Thus, the same alleles that help to define macroevolutionary divergence among genera also contribute to microevolutionary fine-tuning of adaptive traits within some species. Variability of craniofacial morphology mediated by Ptch1 polymorphism has likely contributed to niche partitioning and ecological speciation of these fishes.

Multiple interacting loci control sex determination in lake Malawi cichlid fish.

Several models have been proposed to suggest how the evolution of sex-determining mechanisms might contribute to speciation. Here, we describe the inheritance of sex in 19 fish species from the rapidly evolving flock of cichlids in Lake Malawi, Africa. We found that many of these species have a male heterogametic (XY) system on linkage group 7. Some species also segregate for a female heterogametic (ZW) system on linkage group 5 that is coincident with a dominant orange-blotch (OB) color pattern in females. The ZW system is epistatically dominant to the XY system when both are segregating within a family. Several lines of evidence suggest that additional sex-determining loci are segregating in some species. These results are consistent with the idea that genetic conflicts play an important role in the evolution of these species flocks and suggest that evolution of sex-determining mechanisms has contributed to the radiation of cichlid fish in East Africa.

Sexual conflict resolved by invasion of a novel sex determiner in Lake Malawi cichlid fishes.

Sex determination mechanisms differ among animal species, but it is not clear how these differences evolve. New sex determiners may arise in response to sexual conflicts, which occur when traits benefit one sex but hinder the other. We identified the genetic basis for the orange-blotch (OB) color pattern, a trait under sexually antagonistic selection in the cichlid fish of Lake Malawi, East Africa. The OB phenotype is due to a cis-regulatory mutation in the Pax7 gene. OB provides benefits of camouflage to females but disrupts the species-specific male color patterns used for mate recognition. We suggest that the resulting sexual conflict over the OB allele has been resolved by selection for a novel female sex determination locus that has invaded populations with an ancestral male sex determination system.

Reduced EGFR causes abnormal valvular differentiation leading to calcific aortic stenosis and left ventricular hypertrophy in C57BL/6J but not 129S1/SvImJ mice.

Epidermal growth factor receptor (EGFR) signaling contributes to aortic valve development in mice. Because developmental phenotypes in Egfr-null mice are dependent on genetic background, the hypomorphic Egfr(wa2) allele was made congenic on C57BL/6J (B6) and 129S1/SvImJ (129) backgrounds and used to identify the underlying cellular cause of EGFR-related aortic valve abnormalities. Egfr(wa2/wa2) mice on both genetic backgrounds develop aortic valve hyperplasia. Many B6-Egfr(wa2/wa2) mice die before weaning, and those surviving to 3 mo of age or older develop severe left ventricular hypertrophy and heart failure. The cardiac phenotype was accompanied by significantly thicker aortic cusps and larger transvalvular gradients in B6-Egfr(wa2/wa2) mice compared with heterozygous controls and age-matched Egfr(wa2) homozygous mice on either 129 or B6129F1 backgrounds. Histological analysis revealed cellular changes in B6-Egfr(wa2/wa2) aortic valves underlying elevated pressure gradients and progression to heart failure, including increased cellular proliferation, ectopic cartilage formation, extensive calcification, and inflammatory infiltrate, mimicking changes seen in human calcific aortic stenosis. Despite having congenitally enlarged valves, 129 and B6129F1-Egfr(wa2/wa2) mice have normal lifespans, absence of left ventricular hypertrophy, and normal systolic function. These results show the requirement of EGFR activity for normal valvulogenesis and demonstrate that dominantly acting genetic modifiers curtail pathological changes in congenitally deformed valves. These studies provide a novel model of aortic sclerosis and stenosis and suggest that long-term inhibition of EGFR signaling for cancer therapy may have unexpected consequences on aortic valves in susceptible individuals.

Transcriptional recapitulation and subversion of embryonic colon development by mouse colon tumor models and human colon cancer.

BACKGROUND: The expression of carcino-embryonic antigen by colorectal cancer is an example of oncogenic activation of embryonic gene expression. Hypothesizing that oncogenesis-recapitulating-ontogenesis may represent a broad programmatic commitment, we compared gene expression patterns of human colorectal cancers (CRCs) and mouse colon tumor models to those of mouse colon development embryonic days 13.5-18.5. RESULTS: We report here that 39 colon tumors from four independent mouse models and 100 human CRCs encompassing all clinical stages shared a striking recapitulation of embryonic colon gene expression. Compared to normal adult colon, all mouse and human tumors over-expressed a large cluster of genes highly enriched for functional association to the control of cell cycle progression, proliferation, and migration, including those encoding MYC, AKT2, PLK1 and SPARC. Mouse tumors positive for nuclear beta-catenin shifted the shared embryonic pattern to that of early development. Human and mouse tumors differed from normal embryonic colon by their loss of expression modules enriched for tumor suppressors (EDNRB, HSPE, KIT and LSP1). Human CRC adenocarcinomas lost an additional suppressor module (IGFBP4, MAP4K1, PDGFRA, STAB1 and WNT4). Many human tumor samples also gained expression of a coordinately regulated module associated with advanced malignancy (ABCC1, FOXO3A, LIF, PIK3R1, PRNP, TNC, TIMP3 and VEGF). CONCLUSION: Cross-species, developmental, and multi-model gene expression patterning comparisons provide an integrated and versatile framework for definition of transcriptional programs associated with oncogenesis. This approach also provides a general method for identifying pattern-specific biomarkers and therapeutic targets. This delineation and categorization of developmental and non-developmental activator and suppressor gene modules can thus facilitate the formulation of sophisticated hypotheses to evaluate potential synergistic effects of targeting within- and between-modules for next-generation combinatorial therapeutics and improved mouse models.