Ontogenetic stage-specific quantitative trait loci contribute to divergence in developmental trajectories of sexually dimorphic fins between medaka populations

The genetic regulatory architecture and epigenomic basis for age-related changes in rattlesnake venom

Developmental phenotypic changes can evolve under selection imposed by age- and size-related ecological differences. Many of these changes occur through programmed alterations to gene expression patterns, but the molecular mechanisms and gene-regulatory networks underlying these adaptive changes remain poorly understood. Many venomous snakes, including the eastern diamondback rattlesnake (Crotalus adamanteus), undergo correlated changes in diet and venom expression as snakes grow larger with age, providing models for identifying mechanisms of timed expression changes that underlie adaptive life history traits. By combining a highly contiguous, chromosome-level genome assembly with measures of expression, chromatin accessibility, and histone modifications, we identified cis-regulatory elements and trans-regulatory factors controlling venom ontogeny in the venom glands of C. adamanteus. Ontogenetic expression changes were significantly correlated with epigenomic changes within genes, immediately adjacent to genes (e.g., promoters), and more distant from genes (e.g., enhancers). We identified 37 candidate transcription factors (TFs), with the vast majority being up-regulated in adults. The ontogenetic change is largely driven by an increase in the expression of TFs associated with growth signaling, transcriptional activation, and circadian rhythm/biological timing systems in adults with corresponding epigenomic changes near the differentially expressed venom genes. However, both expression activation and repression contributed to the composition of both adult and juvenile venoms, demonstrating the complexity and potential evolvability of gene regulation for this trait. Overall, given that age-based trait variation is common across the tree of life, we provide a framework for understanding gene-regulatory-network-driven life-history evolution more broadly.

A Cryptic Sex-Linked Locus Revealed by the Elimination of a Master Sex-Determining Locus in Medaka Fish

AbstractSex chromosomes rapidly turn over in several taxonomic groups. Sex chromosome turnover is generally thought to start with the appearance of a new sex-determining gene on an autosome while an old sex-determining gene still exists, followed by the fixation of the new one. However, we do not know how prevalent the transient state is, where multiple sex-determining loci coexist within natural populations. Here, we removed a Y chromosome with a master male-determining gene DMY from medaka fish using high temperature-induced sex-reversed males. After four generations, the genomic characteristics of a sex chromosome were found on one chromosome, which was an autosome in the original population. Thus, the elimination of a master sex-determining locus can reveal a cryptic locus with a possible sex-determining effect, which can be the seed for sex chromosome turnover. Our results suggest that populations that seem to have a single-locus XY system may have other chromosomal regions with sex-determining effects. In conclusion, the coexistence of multiple sex-determining genes in a natural population may be more prevalent than previously thought. Experimental elimination of a master sex-determining locus may serve as a promising method for finding a locus that can be a protosex chromosome.

The transcriptional state and chromatin landscape of cichlid jaw shape variation across species and environments

Adaptive phenotypes are shaped by a combination of genetic and environmental forces, but how they interact remains poorly understood. Here, we utilize the cichlid oral jaw apparatus to better understand these gene-by-environment effects. First, we employed RNA-seq in bony and ligamentous tissues important for jaw opening to identify differentially expressed genes between species and across foraging environments. We used two Lake Malawi species adapted to different foraging habitats along the pelagic-benthic ecomorphological axis. Our foraging treatments were designed to force animals to employ either suction or biting/scraping, which broadly mimic pelagic or benthic modes of feeding. We found a large number of differentially expressed genes between species, and while we identified relatively few differences between environments, species differences were far more pronounced when they were challenged with a pelagic versus benthic foraging mode. Expression data carried the signature of genetic assimilation, and implicated cell cycle regulation in shaping the jaw across species and environments. Next, we repeated the foraging experiment and performed ATAC-seq procedures on nuclei harvested from the same tissues. Cross-referencing results from both analyses revealed subsets of genes that were both differentially expressed and differentially accessible. This reduced dataset implicated notable candidate genes including the Hedgehog effector, KIAA0586 and the ETS transcription factor, etv4, which connects environmental stress and craniofacial morphogenesis. Taken together, these data provide novel insights into the epigenetic, genetic and cellular bases of species- and environment-specific bone shapes.

Age-dependent genetic architecture across ontogeny of body size in sticklebacks

Heritable variation in traits under natural selection is a prerequisite for evolutionary response. While it is recognized that trait heritability may vary spatially and temporally depending on which environmental conditions traits are expressed under, less is known about the possibility that genetic variance contributing to the expected selection response in a given trait may vary at different stages of ontogeny. Specifically, whether different loci underlie the expression of a trait throughout development and thus providing an additional source of variation for selection to act on in the wild, is unclear. Here we show that body size, an important life-history trait, is heritable throughout ontogeny in the nine-spined stickleback (Pungitius pungitius). Nevertheless, both analyses of quantitative trait loci and genetic correlations across ages show that different chromosomes/loci contribute to this heritability in different ontogenic time-points. This suggests that body size can respond to selection at different stages of ontogeny but that this response is determined by different loci at different points of development. Hence, our study provides important results regarding our understanding of the genetics of ontogeny and opens an interesting avenue of research for studying age-specific genetic architecture as a source of non-parallel evolution.

The genetic basis of a novel reproductive strategy in Sulawesi ricefishes: How modularity and a low number of loci shape pelvic brooding

The evolution of complex phenotypes like reproductive strategies is challenging to understand, as they often depend on multiple adaptations that only jointly result in a specific functionality. Sulawesi ricefishes (Adrianichthyidae) evolved a reproductive strategy termed as pelvic brooding. In contrast to the more common transfer brooding, female pelvic brooders carry an egg bundle connected to their body for weeks until the fry hatches. To examine the genetic architecture of pelvic brooding, we crossed the pelvic brooding Oryzias eversi and the transfer brooding Oryzias nigrimas (species divergence time: ∼3.6 my). We hypothesize, that a low number of loci and modularity have facilitated the rapid evolution of pelvic brooding. Traits associated to pelvic brooding, like rib length, pelvic fin length, and morphology of the genital papilla, were correlated in the parental species but correlations were reduced or lost in their F1 and F2 hybrids. Using the Castle-Wright estimator, we found that generally few loci underlie the studied traits. Further, both parental species showed modularity in their body plans. In conclusion, morphological traits related to pelvic brooding were based on a few loci and the mid-body region likely could evolve independently from the remaining body parts. Both factors presumably facilitated the evolution of pelvic brooding.

Genome editing reveals fitness effects of a gene for sexual dichromatism in Sulawesian fishes

Sexual selection drives rapid phenotypic diversification of mating traits. However, we know little about the causative genes underlying divergence in sexually selected traits. Here, we investigate the genetic basis of male mating trait diversification in the medaka fishes (genus Oryzias) from Sulawesi, Indonesia. Using linkage mapping, transcriptome analysis, and genome editing, we identify csf1 as a causative gene for red pectoral fins that are unique to male Oryzias woworae. A cis-regulatory mutation enables androgen-induced expression of csf1 in male fins. csf1-knockout males have reduced red coloration and require longer for mating, suggesting that coloration can contribute to male reproductive success. Contrary to expectations, non-red males are more attractive to a predatory fish than are red males. Our results demonstrate that integrating genomics with genome editing enables us to identify causative genes underlying sexually selected traits and provides a new avenue for testing theories of sexual selection.

Latitudinal variation in sexual dimorphism in a freshwater fish group

Tropical animals are characterized by showy ornaments and conspicuous body colours as compared with their temperate relatives. Some recent studies have hypothesized that sexual selection pressures are stronger in the tropics than in the temperate zone. Although negative correlations between latitude and the degree of sexual dimorphism would support this hypothesis, phylogeny should be taken into account in such comparative studies. Comparisons of the degree of sexual dimorphism in body size and fin lengths among species of the Adrianichthyidae, a freshwater fish family having a wide geographical range throughout Southeast and East Asia, revealed that lower latitude species are sexually more dimorphic in all characters than higher latitude species. Phylogenetic generalized least squares analyses using a mitochondrial DNA phylogeny demonstrated that the negative correlations between latitude and the degree of sexual dimorphism become non-significant when phylogeny is considered, but that the variance in the degree of sexual dimorphism is explained not only by phylogeny but also almost equally by latitude. Ancestral state reconstruction indicated that sexual dimorphisms have evolved independently even within major clades. These findings are consistent with the view that tropical species are exposed to stronger sexual selection pressures than temperate species. We discuss possible causes of the latitudinal variation in sexual selection pressure.

QTL Fine Mapping for Sex Determination Region in Bighead Carp (Hypophthalmichthys nobilis) and Comparison with Silver Carp (Hypophthalmichthys molitrix)

Bighead carp (Hypophthalmichthys nobilis) and silver carp (Hypophthalmichthys molitrix) are genetically close aquaculture fish in the Cyprinidae, which have been confirmed to hold XX/XY sex determination. However, genomic locations of potential sex-related loci in these two fishes are still unknown. In this study, a high-resolution genetic linkage map was constructed by using 2976 SNP and 924 microsatellite markers in a F1 full-sib family of bighead carp, the length of which spanned 2022.34 cM with an average inter-marker distance of 0.52 cM. Comparative genomics revealed a high level of genomic synteny between bighead carp and zebrafish as well as grass carp. QTL fine mapping for sex trait was performed based on this linkage map of bighead carp and an unpublished linkage map of silver carp. A map distance of 3.863 cM (69.787-73.650 cM) on LG19 of bighead carp and 4.705 cM (79.096-83.801 cM) on LG21 of silver carp was significantly associated with sex phenotypes, and these two LGs are homologous between two fish species. Fourteen markers harboring in these regions were in strong linkage disequilibrium with the sex phenotype variance explained (PVE) varying from 89 to 100%. Two common markers were mapped on the QTL regions of bighead carp and silver carp, suggesting that these two carp species may have similar genetic bases for sex determination. Eleven potentially sex-related genes were identified within or near the sex QTL markers in two species. This study provided insights into elucidating mechanisms and evolution of sex determination in cyprinid fishes.

Contribution of gene flow to the evolution of recombination suppression in sex chromosomes

Polymorphism of alleles that benefit one sex but harm the other (sexually antagonistic alleles) generates selective pressures for reduced recombination between themselves and sex-determination loci. Such polymorphism can be maintained within a population when selection coefficients are sufficiently balanced between males and females. However, if regulatory mutations restrict gene expression only to one sex, these alleles become neutral in the other sex and easily fixed within a population, removing the selective pressures for recombination suppression in sex chromosomes. When there is spatial variation in selection regimes, however, alleles that are deleterious in one sex and neutral in the other can be maintained in other neighboring populations and gene flow may continuously supply deleterious alleles. We hypothesized that this maintenance of genetic variation may promote the establishment of recombination suppression in sex chromosomes even in cases where selection is limited to one sex. Using individual-based simulations, we show that spatial variation in male-limited selection and gene flow can promote the establishment of Y-autosome fusions, a special case of recombination suppression in sex chromosomes. This can be explained by the fact that fused Y-chromosomes that capture alleles that are beneficial for local males have a higher mean fitness compared to unfused Y chromosomes in the presence of deleterious gene flow. We also simulated the case of sex-concordant selection and found that gene flow of alleles that are deleterious in both sexes did not substantially increase the establishment rates of Y-autosome fusions across the parameter space examined. This can be accounted for by the fact that foreign alleles that are deleterious in both sexes can be efficiently removed from the population compared to alleles that are neutral in females. These results indicate that how gene flow affects the establishment rates of Y-autosome fusions depends largely on selection regimes. Spatial variation in sex-specific selection and gene flow should be appreciated as a factor affecting sex chromosome evolution.

Pleiotropy, constraint, and modularity in the evolution of life histories: insights from genomic analyses

Multicellular organisms display an enormous range of life history (LH) strategies and present an evolutionary conundrum; despite strong natural selection, LH traits are characterized by high levels of genetic variation. To understand the evolution of life histories and maintenance of this variation, the specific phenotypic effects of segregating alleles and the genetic networks in which they act need to be elucidated. In particular, the extent to which LH evolution is constrained by the pleiotropy of alleles contributing to LH variation is generally unknown. Here, we review recent empirical results that shed light on this question, with an emphasis on studies employing genomic analyses. While genome-scale analyses are increasingly practical and affordable, they face limitations of genetic resolution and statistical power. We describe new research approaches that we believe can produce new insights and evaluate their promise and applicability to different kinds of organisms. Two approaches seem particularly promising: experiments that manipulate selection in multiple dimensions and measure phenotypic and genomic response and analytical approaches that take into account genome-wide associations between markers and phenotypes, rather than applying a traditional marker-by-marker approach.

Variation in the ontogeny of sex steroid levels between latitudinal populations of the medaka

INTRODUCTION

Sex steroids mediate the expression of sexual dimorphism during ontogeny, and populations that differ in the magnitudes of sexual dimorphism may accordingly differ in the ontogenetic patterns of their sex steroid levels. The medaka, Oryzias latipes species complex, shows geographic variation in the magnitude of sexual dimorphism with respect to the lengths of their anal and dorsal fins; dimorphism is greater in low-latitude populations than in high-latitude populations. However, sexual differences in the ontogenetic dynamics of sex steroids, and its interpopulation variation, have not been examined.

RESULTS

We measured testosterone (T), estradiol-17β (E2), and 11-ketotestosterone (11-KT) concentrations throughout ontogeny of laboratory-reared fish from two latitudinal populations: Aomori (northern) and Okinawa (southern). In both populations, the levels of all three steroids were high during early ontogenetic stages and decreased with growth. After reaching about 15 mm in standard length, when sexual dimorphisms in fin lengths became apparent, steroid levels increased and tended to plateau. Sexual differences in the steroid levels were observed only in the later ontogenetic stages; T and 11-KT levels were higher in males, while E2 levels were higher in females. Accordingly, interpopulation differences also became clearer; the southern fish tended to show higher T levels and lower E2 levels than the northern fish.

CONCLUSIONS

The ontogenetic patterns of sex steroid levels paralleled the ontogeny of anal and dorsal fins in the two latitudinal populations, suggesting that interpopulation variation in the degree of sexual dimorphisms in fin lengths is mediated by sex steroid-dependent regulation of fin elongation.

Genetic Architecture of the Variation in Male-Specific Ossified Processes on the Anal Fins of Japanese Medaka

Traits involved in reproduction evolve rapidly and show great diversity among closely related species. However, the genetic mechanisms that underlie the diversification of courtship traits are mostly unknown. Japanese medaka fishes (Oryzias latipes) use anal fins to attract females and to grasp females during courtship; the males have longer anal fins with male-specific ossified papillary processes on the fin rays. However, anal fin morphology varies between populations: the southern populations tend to have longer anal fins and more processes than the northern populations. In the present study, we conducted quantitative trait locus (QTL) mapping to investigate the genetic architecture underlying the variation in the number of papillary processes of Japanese medaka fish and compared the QTL with previously identified QTL controlling anal fin length. First, we found that only a few QTL were shared between anal fin length and papillary process number. Second, we found that the numbers of papillary processes on different fin rays often were controlled by different QTL. Finally, we produced another independent cross and found that some QTL were repeatable between the two crosses, whereas others were specific to only one cross. These results suggest that variation in the number of papillary processes is polygenic and controlled by QTL that are distinct from those controlling anal fin length. Thus, different courtship traits in Japanese medaka share a small number of QTL and have the potential for independent evolution.

Phenotypic integration in the feeding system of the eastern diamondback rattlesnake (Crotalus adamanteus)

Selection can vary geographically across environments and temporally over the lifetime of an individual. Unlike geographic contexts, where different selective regimes can act on different alleles, age-specific selection is constrained to act on the same genome by altering age-specific expression. Snake venoms are exceptional traits for studying ontogeny because toxin expression variation directly changes the phenotype; relative amounts of venom components determine, in part, venom efficacy. Phenotypic integration is the dependent relationship between different traits that collectively produce a complex phenotype and, in venomous snakes, may include traits as diverse as venom, head shape and fang length. We examined the feeding system of the eastern diamondback rattlesnake (Crotalus adamanteus) across environments and over the lifetime of individuals and used a genotype-phenotype map approach, protein expression data and morphological data to demonstrate that: (i) ontogenetic effects explained more of the variation in toxin expression variation than geographic effects, (ii) both juveniles and adults varied geographically, (iii) toxin expression variation was a result of directional selection and (iv) different venom phenotypes covaried with morphological traits also associated with feeding in temporal (ontogenetic) and geographic (functional) contexts. These data are the first to demonstrate, to our knowledge, phenotypic integration between multiple morphological characters and a biochemical phenotype across populations and age classes. We identified copy number variation as the mechanism driving the difference in the venom phenotype associated with these morphological differences, and the parallel mitochondrial, venom and morphological divergence between northern and southern clades suggests that each clade may warrant classification as a separate evolutionarily significant unit.