Abstract
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.