Genetic analysis of vertebral regionalization and number in medaka (Oryzias latipes) inbred lines

Adaptive tail-length evolution in deer mice is associated with differential Hoxd13 expression in early development

Variation in the size and number of axial segments underlies much of the diversity in animal body plans. Here we investigate the evolutionary, genetic and developmental mechanisms driving tail-length differences between forest and prairie ecotypes of deer mice (Peromyscus maniculatus). We first show that long-tailed forest mice perform better in an arboreal locomotion assay, consistent with tails being important for balance during climbing. We then identify six genomic regions that contribute to differences in tail length, three of which associate with caudal vertebra length and the other three with vertebra number. For all six loci, the forest allele increases tail length, indicative of the cumulative effect of natural selection. Two of the genomic regions associated with variation in vertebra number contain Hox gene clusters. Of those, we find an allele-specific decrease in Hoxd13 expression in the embryonic tail bud of long-tailed forest mice, consistent with its role in axial elongation. Additionally, we find that forest embryos have more presomitic mesoderm than prairie embryos and that this correlates with an increase in the number of neuromesodermal progenitors, which are modulated by Hox13 paralogues. Together, these results suggest a role for Hoxd13 in the development of natural variation in adaptive morphology on a microevolutionary timescale.

High-speed system to generate congenic strains in medaka

The congenic strain, an inbred strain containing a small genomic region from another strain, is a powerful tool to assess the phenotypic effect of polymorphisms and/or mutations in the substituted genomic region. Recent substantial progress in the genetic studies of complex traits increases the necessity of congenic strains and, therefore, a quick breeding system for congenic strains has become increasingly important in model organisms such as mouse and medaka. Traditionally, more than ten generations are necessary to produce a congenic strain. In contrast, a quick method has been reported previously for the mouse, in which the use of genetic markers reduces the required number of backcross generations to about a half that of the traditional method, so that it would take around six generations to obtain a congenic strain. Here, we present an even quicker congenic production system, which takes only about four generations. The system can produce medaka congenic strains having part of the HNI-II (an inbred medaka strain derived from the northern Japanese population, Oryzias sakaizumii) genome in the HdrR-II1 (another inbred strain from the southern Japanese population, O. latipes) background. In this system, the availability of frozen sperm and genotype data of the BC1 male population makes it possible to start marker-assisted congenic production after obtaining the BC2 population. Our evaluation revealed that the system could work well to increase the percentage of recipient genome as expected, so that a congenic strain may be obtained in about one year.

A Baseline for Skeletal Investigations in Medaka (Oryzias latipes): The Effects of Rearing Density on the Postcranial Phenotype

Oryzias latipes is increasingly used as a model in biomedical skeletal research. The standard approach is to generate genetic variants with particular skeletal phenotypes which resemble skeletal diseases in humans. The proper diagnosis of skeletal variation is key for this type of research. However, even laboratory rearing conditions can alter skeletal phenotypes. The subject of this study is the link between skeletal phenotypes and rearing conditions. Thus, wildtype medaka were reared from hatching to an early juvenile stage at low (LD: 5 individuals/L), medium (MD: 15 individuals/L), and high (HD: 45 individuals/L) densities. The objectives of the study are: (I) provide a comprehensive overview of the postcranial skeletal elements in medaka; (II) evaluate the effects of rearing density on specific meristic counts and on the variability in type and incidence of skeletal anomalies; (III) define the best laboratory settings to obtain a skeletal reference for a sound evaluation of future experimental conditions; (IV) contribute to elucidating the structural and cellular changes related to the onset of skeletal anomalies. The results from this study reveal that rearing densities greater than 5 medaka/L reduce the animals' growth. This reduction is related to decreased mineralization of dermal (fin rays) and perichondral (fin supporting elements) bone. Furthermore, high density increases anomalies affecting the caudal fin endoskeleton and dermal rays, and the preural vertebral centra. A series of static observations on Alizarin red S whole mount-stained preural fusions provide insights into the etiology of centra fusion. The fusion of preural centra involves the ectopic formation of bony bridges over the intact intervertebral ligament. An apparent consequence is the degradation of the intervertebral ligaments and the remodeling and reshaping of the fused vertebral centra into a biconoid-shaped centrum. From this study it can be concluded that it is paramount to take into account the rearing conditions, natural variability, skeletal phenotypic plasticity, and the genetic background along with species-specific peculiarities when screening for skeletal phenotypes of mutant or wildtype medaka.

Eight in One: Hidden Diversity of the Bagrid Catfish Tachysurus albomarginatus s.l. (Rendhal, 1928) Widespread in Lowlands of South China

There is increasing evidence that species diversity is underestimated in the current taxonomy of widespread freshwater fishes. The bagrid species T. albomarginatus s.l. is mainly distributed in the lowlands of South China, as currently identified. A total of 40 localities (including the type locality), which covers most of its known range, were sampled. Molecular phylogenetic analyses based on concatenated mtDNA and nuclear genes recover nine highly supported lineages clustering into eight geographic populations. The integration of molecular evidence, morphological data, and geographic distribution demonstrates the delineation of T. albomarginatus s.l. as eight putative species. Four species, namely, T. albomarginatus, T. lani, T. analis, and T. zhangfei sp. nov. and the T. similis complex are taxonomically recognized herein. Moreover, T. zhangfei sp. nov. comprises two genetically distinct lineages with no morphological and geographical difference. This study also reveals aspects of estimation of divergence time, distribution, and ecological adaption within the T. albomarginatus group. The unraveling of the hidden species diversity of this lowland bagrid fish highlights the need for not only the molecular scrutiny of widely distributed species of South China but also the adjustment of current biodiversity conservation strategies to protect the largely overlooked diversity of fishes from low-elevation rapids.

Genetic analysis of body weight in wild populations of medaka fish from different latitudes

The genetic bases of growth and body weight are of economic and scientific interest, and teleost fish models have proven useful in such investigations. The Oryzias latipes species complex (medaka) is an abundant freshwater fish in Japan and suitable for genetic studies. We compared two wild medaka stocks originating from different latitudes. The Maizuru population from higher latitudes weighed more than the Ginoza population. We investigated the genetic basis of body weight, using quantitative trait locus (QTL) analysis of the F2 offspring of these populations. We detected one statistically significant QTL for body weight on medaka chromosome 4 and identified 12 candidate genes that might be associated with body weight or growth. Nine of these 12 genes had at least one single nucleotide polymorphism that caused amino acid substitutions in protein-coding regions, and we estimated the effects of these substitutions. The present findings might contribute to the marker-assisted selection of economically important aquaculture species.

Effects of temperature and water turbulence on vertebral number and body shape in Astyanax mexicanus (Teleostei: Characidae)

Environmental changes can modify the phenotypic characteristics of populations, which in turn can influence their evolutionary trajectories. In ectotherms like fishes, temperature is a particularly important environmental variable that is known to have significant impacts on the phenotype. Here, we raised specimens of the surface ecomorph of Astyanax mexicanus at temperatures of 20°C, 23°C, 25°C, and 28°C to examine how temperature influenced vertebral number and body shape. To increase biological realism, specimens were also subjected to two water turbulence regimes. Vertebral number was counted from x-rays and body shape variation was analysed using geometric morphometric methods. Temperature significantly impacted mean total vertebral number, which increased at the lowest and highest temperatures. Fish reared at lower temperatures had relatively more precaudal vertebrae while fish reared at higher temperatures had relatively more caudal vertebrae. Vertebral anomalies, especially vertebral fusions, were most frequent at the extreme temperature treatments. Temperature significantly impacted body shape as well, with fish reared at 20°C being particularly divergent. Water turbulence also impacted body shape in a generally predictable manner, with specimens reared in high turbulence environments being more streamlined and having extended dorsal and anal fin bases. Variation in environmental variables thus resulted in significant changes in morphological traits known to impact fish fitness, indicating that A. mexicanus has the capacity to exhibit a range of phenotypic plasticity when challenged by environmental change. Understanding the biochemical mechanisms underlying this plasticity and whether adaptive plasticity has influenced the evolutionary radiation of the Characidae, are major directions for future research.

Exaggerated heterochiasmy in a fish with sex-linked male coloration polymorphisms

It is often stated that polymorphisms for mutations affecting fitness of males and females in opposite directions [sexually antagonistic (SA) polymorphisms] are the main selective force for the evolution of recombination suppression between sex chromosomes. However, empirical evidence to discriminate between different hypotheses is difficult to obtain. We report genetic mapping results in laboratory-raised families of the guppy (Poecilia reticulata), a sexually dimorphic fish with SA polymorphisms for male coloration genes, mostly on the sex chromosomes. Comparison of the genetic and physical maps shows that crossovers are distributed very differently in the two sexes (heterochiasmy); in male meiosis, they are restricted to the termini of all four chromosomes studied, including chromosome 12, which carries the sex-determining locus. Genome resequencing of male and female guppies from a population also indicates sex linkage of variants across almost the entire chromosome 12. More than 90% of the chromosome carrying the male-determining locus is therefore transmitted largely through the male lineage. A lack of heterochiasmy in a related fish species suggests that it originated recently in the lineage leading to the guppy. Our findings do not support the hypothesis that suppressed recombination evolved in response to the presence of SA polymorphisms. Instead, a low frequency of recombination on a chromosome that carries a male-determining locus and has not undergone genetic degeneration has probably facilitated the establishment of male-beneficial coloration polymorphisms.

pnp4a Is the Causal Gene of the Medaka Iridophore Mutant guanineless

See-through medaka lines are suitable for observing internal organs throughout life. They were bred by crossing multiple color mutants. However, some of the causal genes for these mutants have not been identified. The medaka has four pigment cell types: black melanophores, yellow xanthophores, white leucophores, and silvery iridophores. The causal genes of melanophore, xanthophore, and leucophore mutants have been elucidated, but the causal gene for the iridophore mutant remains unknown. Here, we describe the iridophore mutant, guanineless (gu), which exhibits a strong reduction in visible iridophores throughout its larval to adult stages. The gu locus was previously mapped to chromosome 5, but was located near the telomeric region, making it difficult to integrate into the chromosome. We sought the causal gene of gu using synteny analysis with the zebrafish genome and found a strong candidate, purine nucleoside phosphorylase 4a (pnp4a). Gene targeting and complementation testing showed that pnp4a is the causal gene of gu. This result will allow the establishment of inbred medaka strains or other useful strains with see-through phenotypes without major disruption in the genetic background of each strain.

Genetic control of startle behavior in medaka fish

Genetic polymorphisms are thought to generate intraspecific behavioral diversities, both within and among populations. The mechanisms underlying genetic control of behavioral properties, however, remain unclear in wild-type vertebrates, including humans. To explore this issue, we used diverse inbred strains of medaka fish (Oryzias latipes) established from the same and different local populations. Medaka exhibit a startle response to a visual stimulus (extinction of illumination) by rapidly bending their bodies (C-start) 20-ms after the stimulus presentation. We measured the rates of the response to repeated stimuli (1-s interval, 40 times) among four inbred strains, HNI-I, HNI-II, HO5, and Hd-rR-II1, and quantified two properties of the startle response: sensitivity (response rate to the first stimulus) and attenuation of the response probability with repeated stimulus presentation. Among the four strains, the greatest differences in these properties were detected between HNI-II and Hd-rR-II1. HNI-II exhibited high sensitivity (approximately 80%) and no attenuation, while Hd-rR-II1 exhibited low sensitivity (approximately 50%) and almost complete attenuation after only five stimulus presentations. Our findings suggested behavioral diversity of the startle response within a local population as well as among different populations. Linkage analysis with F2 progeny between HNI-II and Hd-rR-II1 detected quantitative trait loci (QTL) highly related to attenuation, but not to sensitivity, with a maximum logarithm of odds score of 11.82 on linkage group 16. The three genotypes (homozygous for HNI-II and Hd-rR-II1 alleles, and heterozygous) at the marker nearest the QTL correlated with attenuation. Our findings are the first to suggest that a single genomic region might be sufficient to generate individual differences in startle behavior between wild-type strains. Further identification of genetic polymorphisms that define the behavioral trait will contribute to our understanding of the neural mechanisms underlying behavioral diversity, allowing us to investigate the adaptive significance of intraspecific behavioral polymorphisms of the startle response.

Analysis of a novel gene, Sdgc, reveals sex chromosome-dependent differences of medaka germ cells prior to gonad formation

In vertebrates that have been examined to date, the sexual identity of germ cells is determined by the sex of gonadal somatic cells. In the teleost fish medaka, a sex-determination gene on the Y chromosome, DMY/dmrt1bY, is expressed in gonadal somatic cells and regulates the sexual identity of germ cells. Here, we report a novel mechanism by which sex chromosomes cell-autonomously confer sexually different characters upon germ cells prior to gonad formation in a genetically sex-determined species. We have identified a novel gene, Sdgc (sex chromosome-dependent differential expression in germ cells), whose transcripts are highly enriched in early XY germ cells. Chimeric analysis revealed that sexually different expression of Sdgc is controlled in a germ cell-autonomous manner by the number of Y chromosomes. Unexpectedly, DMY/dmrt1bY was expressed in germ cells prior to gonad formation, but knockdown and overexpression of DMY/dmrt1bY did not affect Sdgc expression. We also found that XX and XY germ cells isolated before the onset of DMY/dmrt1bY expression in gonadal somatic cells behaved differently in vitro and were affected by Sdgc. Sdgc maps close to the sex-determination locus, and recombination around the two loci appears to be repressed. Our results provide important insights into the acquisition and plasticity of sexual differences at the cellular level even prior to the developmental stage of sex determination.

Regionalization of the axial skeleton in the 'ambush predator' guild--are there developmental rules underlying body shape evolution in ray-finned fishes?

Background

A long, slender body plan characterized by an elongate antorbital region and posterior displacement of the unpaired fins has evolved multiple times within ray-finned fishes, and is associated with ambush predation. The axial skeleton of ray-finned fishes is divided into abdominal and caudal regions, considered to be evolutionary modules. In this study, we test whether the convergent evolution of the ambush predator body plan is associated with predictable, regional changes in the axial skeleton, specifically whether the abdominal region is preferentially lengthened relative to the caudal region through the addition of vertebrae. We test this hypothesis in seven clades showing convergent evolution of this body plan, examining abdominal and caudal vertebral counts in over 300 living and fossil species. In four of these clades, we also examined the relationship between the fineness ratio and vertebral regionalization using phylogenetic independent contrasts.

Results

We report that in five of the clades surveyed, Lepisosteidae, Esocidae, Belonidae, Sphyraenidae and Fistulariidae, vertebrae are added preferentially to the abdominal region. In Lepisosteidae, Esocidae, and Belonidae, increasing abdominal vertebral count was also significantly related to increasing fineness ratio, a measure of elongation. Two clades did not preferentially add abdominal vertebrae: Saurichthyidae and Aulostomidae. Both of these groups show the development of a novel caudal region anterior to the insertion of the anal fin, morphologically differentiated from more posterior caudal vertebrae.

Conclusions

The preferential addition of abdominal vertebrae in fishes with an elongate body shape is consistent with the existence of a conservative positioning module formed by the boundary between the abdominal and caudal vertebral regions and the anterior insertion of the anal fin. Dissociation of this module is possible, although less probable than changes in the independently evolving abdominal region. Dissociation of the axial skeleton-median fin module leads to increased regionalization within the caudal vertebral column, something that has evolved several times in bony fishes, and may be homologous with the sacral region of tetrapods. These results suggest that modularity of the axial skeleton may result in somewhat predictable evolutionary outcomes in bony fishes.