An evolutionary test of the isoform switching hypothesis of functional progesterone withdrawal for parturition: humans have a weaker repressive effect of PR-A than mice

BACKGROUND

A decrease in maternal serum progesterone (P4) concentrations precedes the onset of labor in most placental mammals. Humans differ by maintaining high levels of P4 throughout birth. Parturition in humans probably includes mechanisms that undercut the pregnancy sustaining function of P4. One attractive hypothesis is the isoform switching hypothesis (ISH). ISH is supported by in vitro evidence that progesterone receptor isoform A (PR-A) inhibits PR-B and that the PR-A/PR-B ratio increases towards term.

MATERIALS AND METHODS

Here, we test the hypothesis that isoform switching is an adaptation to high levels of P4 at term, predicting that, in humans, PR-A mediated repression of PR-B is stronger than in mouse. We use reporter assays with human and mouse PRs to detect species differences in the repressive effects of PR-A.

RESULTS

We found that human PR-B is less sensitive to repression by human PR-A than mouse PR-B, contrary to our prediction. The difference between human and mouse PR-B sensitivity is most pronounced at PR-A/PR-B ratios typical for the preterm myometrium.

CONCLUSIONS

Our results are inconsistent with the ISH. We speculate that, instead, the lower sensitivity of human PR-B to PR-A may be relevant for the maintenance of pregnancy at high progesterone levels and increasing PR-A concentrations towards term.

The evolutionary dynamics of evolvability in a gene network model

Evolvability, the ability of populations to adapt, has recently emerged as a major unifying concept in biology. Although the study of evolvability offers new insights into many important biological questions, the conceptual bases of evolvability, and the mechanisms of its evolution, remain controversial. We used simulated evolution of a model of gene network dynamics to test the contentious hypothesis that natural selection can favour high evolvability, in particular in sexual populations. Our results conclusively demonstrate that fluctuating natural selection can increase the capacity of model gene networks to adapt to new environments. Detailed studies of the evolutionary dynamics of these networks establish a broad range of validity for this result and quantify the evolutionary forces responsible for changes in evolvability. Analysis of the genotype-phenotype map of these networks also reveals mechanisms connecting evolvability, genetic architecture and robustness. Our results suggest that the evolution of evolvability can have a pervasive influence on many aspects of organisms.

The topology of the possible: formal spaces underlying patterns of evolutionary change

The current implementation of the Neo-Darwinian model of evolution typically assumes that the set of possible phenotypes is organized into a highly symmetric and regular space equipped with a notion of distance, for example, a Euclidean vector space. Recent computational work on a biophysical genotype-phenotype model based on the folding of RNA sequences into secondary structures suggests a rather different picture. If phenotypes are organized according to genetic accessibility, the resulting space lacks a metric and is formalized by an unfamiliar structure, known as a pre-topology. Patterns of phenotypic evolution-such as punctuation, irreversibility, modularity--result naturally from the properties of this space. The classical framework, however, addresses these patterns by exclusively invoking natural selection on suitably imposed fitness landscapes. We propose to extend the explanatory level for phenotypic evolution from fitness considerations alone to include the topological structure of phenotype space as induced by the genotype-phenotype map. We introduce the mathematical concepts and tools necessary to formalize the notion of accessibility pre-topology relative to which we can speak of continuity in the genotype-phenotype map and in evolutionary trajectories. We connect the factorization of a pre-topology into a product space with the notion of phenotypic character and derive a condition for factorization. Based on anecdotal evidence from the RNA model, we conjecture that this condition is not globally fulfilled, but rather confined to regions where the genotype-phenotype map is continuous. Equivalently, local regions of genotype space on which the map is discontinuous are associated with the loss of character autonomy. This is consistent with the importance of these regions for phenotypic innovation. The intention of the present paper is to offer a perspective, a framework to implement this perspective, and a few results illustrating how this framework can be put to work. The RNA case is used as an example throughout the text.

The tetrapod limb: a hypothesis on its origin

A wrist joint and structures typical of the hand, such as digits, however, are absent in [Eustenopteron] (Andrews and Westoll, '68, p 240). Great changes must have been undergone during evolution of the ankle joint; the small number of large bones in the fin must somehow have developed into a large number of small bones, and it is very difficult to draw homologies in this region, or even be certain of what is being compared (Andrews and Westoll, '68, p 268). The tetrapod limb is one of the major morphological adaptations that facilitated the transition from an aquatic to a terrestrial lifestyle in vertebrate evolution. We review the paleontological evidence for the fin-limb transition and conclude that the innovation associated with evolution of the tetrapod limb is the zeugopodial-mesopodial transition, i.e., the evolution of the developmental mechanism that differentiates the distal parts of the limb (the autopodium, i.e., hand or foot) from the proximal parts. Based on a review of tetrapod limb and fish fin development, we propose a genetic hypothesis for the origin of the autopodium. In tetrapods the genes Hoxa-11 and Hoxa-13 have locally exclusive expression domains along the proximal-distal axis of the limb bud. The junction between the distal limit of Hoxa-11 expression and of the proximal limit of Hoxa-13 expression is involved in establishing the border between the zeugopodial and autopodial anlagen. In zebrafish, the expression domains of these genes are overlapping and there is no evidence for an autopodial equivalent in the fin skeleton. We propose that the evolution of the derived expression patterns of Hoxa-11 and Hoxa-13 may be causally involved in the origin of the tetrapod limb.

Epistasis and the mutation load: a measurement-theoretical approach

An approximate solution for the mean fitness in mutation-selection balance with arbitrary order of epistatic interaction is derived. The solution is based on the assumptions of coupling equilibrium and that the interaction effects are multilinear. We find that the effect of m-order epistatic interactions (i.e., interactions among groups of m loci) on the load is dependent on the total genomic mutation rate, U, to the mth power. Thus, higher-order gene interactions are potentially important if U is large and the interaction density among loci is not too low. The solution suggests that synergistic epistasis will decrease the mutation load and that variation in epistatic effects will elevate the load. Both of these results, however, are strictly true only if they refer to epistatic interaction strengths measured in the optimal genotype. If gene interactions are measured at mutation-selection equilibrium, only synergistic interactions among even numbers of genes will reduce the load. Odd-ordered synergistic interactions will then elevate the load. There is no systematic relationship between variation in epistasis and load at equilibrium. We argue that empirical estimates of gene interaction must pay attention to the genetic background in which the effects are measured and that it may be advantageous to refer to average interaction intensities as measured in mutation-selection equilibrium. We derive a simple criterion for the strength of epistasis that is necessary to overcome the twofold disadvantage of sex.

Modeling genetic architecture: a multilinear theory of gene interaction

The map from genotype to phenotype is an exceedingly complex function of central importance in biology. In this work we derive and analyze a mathematically tractable model of the genotype-phenotype map that allows for any order of gene interaction. By assuming that the alterations of the effect of a gene substitution due to changes in the genetic background can be described as a linear transformation, we show that the genotype-phenotype map is a sum of linear and multilinear terms of operationally defined "reference" effects at each locus. The "multilinear" model is used to study the effect of epistasis on quantitative genetic variation, on the response to selection, and on genetic canalization. It is shown how the model can be used to estimate the strength of "functional" epistasis from a variety of genetic experiments.

Evolution of Hoxa-11 in lineages phylogenetically positioned along the fin-limb transition

HOXA11 is a transcription factor implicated in paired appendage development. To identify signatures of evolutionary change in the structural, and putative functional, domains of HOXA11, we studied its evolution in tetrapod and nontetrapod lineages that represent approximately 1.5 billion years of evolutionary time. Here, Hoxa-11 gene proper sequences were determined for frog (Xenopus tropicalis), coelacanth (Latimeria chalumnae), common zebrafish (Danio rerio; Hoxa-11a and Hoxa-11b paralogs), and giant zebrafish (D. aequipinnatus; Hoxa-11b) and aligned against previously published Hoxa-11 sequences of human, mouse, chick, and newt. Based on aligned Hoxa-11 amino acid sequences, the protein was demarcated into three segments: Domains I (N-terminal) and III (homeobox + C-terminal), which varied slightly in rates and patterns of evolution, and a variable, overall hydrophilic region (HyD), which partially overlaps with Domain I. As judged by character reconstructions of HOXA11 Domains I and III, no significant changes in rates of coding sequence evolution occurred in tetrapods (frog and chick), relative to coelacanth (a lobe-finned fish), i.e., across the fin-limb transition. Accelerated rates of Hoxa-11 coding sequence evolution were observed for the mammalian and newt lineages. This was shown to be a gene-specific phenomenon. The duplicated Hoxa-11a and Hoxa-11b genes of zebrafish exhibited accelerated rates of evolution and accumulated substitutions at sites that are conserved among coelacanth and all tetrapods examined. Amino acid sequence comparisons of the HyD of HOXA11 suggested that a putative repressor subdomain, containing stretches of consecutive alanine residues, emerged within the tetrapods. A high degree of nucleotide conservation in the 5' half of the Hoxa-11 intron was observed for tetrapod and nontetrapod lineages. Using electrophoretic mobility shift assays, a 35-bp intron sequence, which is 100% conserved in all Hoxa-11 loci except for the zebrafish Hoxa-11a paralog, was found to bind protein(s) in HeLa and chick whole-cell extracts.

Is the genotype-phenotype map modular? A statistical approach using mouse quantitative trait loci data

Various theories about the evolution of complex characters make predictions about the statistical distribution of genetic effects on phenotypic characters, also called the genotype-phenotype map. With the advent of QTL technology, data about these distributions are becoming available. In this article, we propose simple tests for the prediction that functionally integrated characters have a modular genotype-phenotype map. The test is applied to QTL data on the mouse mandible. The results provide statistical support for the notion that the ascending ramus region of the mandible is modularized. A data set comprising the effects of QTL on a more extensive portion of the phenotype is required to determine if the alveolar region of the mandible is also modularized.

Population dependent Fourier decomposition of fitness landscapes over recombination spaces: evolvability of complex characters

The effect of recombination on genotypes can be represented in the form of P-structures, i.e., a map from the set of pairs of genotypes to the power set of genotypes. The interpretation is that the P-structure maps the pair of parental genotypes to the set of recombinant genotypes which result from the recombination of the parental genotypes. A recombination fitness landscape is then a function from the genotypes in a P-structure to the real numbers. In previous papers we have shown that the eigenfunctions of (a matrix associated with) the P-structure provide a basis for the Fourier decomposition of arbitrary recombination landscapes. Here we generalize this framework to include the effect of genotype frequencies, assuming linkage equilibrium. We find that the autocorrelation of the eigenfunctions of the population-weighted P-structure is independent of the population composition. As a consequence we can directly compare the performance of mutation and recombination operators by comparing the autocorrelations on the finite set of elementary landscapes. This comparison suggests that point mutation is a superior search strategy on landscapes with a low order and a moderate order of interaction p < n/3 (n is the number of loci). For more complex landscapes 1-point recombination is superior to both mutation and uniform recombination, but only if the distance among the interacting loci (defining length) is minimal. Furthermore we find that the autocorrelation on any landscape is increasing as the distribution of genotypes becomes more extreme, i.e., if the population occupies a location close to the boundary of the frequency simplex. Landscapes are smoother the more biased the distribution of genotype frequencies is. We suggest that this result explains the paradox that there is little epistatic interaction for quantitative traits detected in natural populations if one uses variance decomposition methods while there is evidence for strong interactions in molecular mapping studies for quantitative trait loci.

Modeling the evolution of genetic architecture: A continuum of alleles model with pairwise AxA epistasis

A continuum of alleles model with pair-wise AxA epistasis is proposed and its transmission genetic, and variational properties are analysed. The basic idea is that genes control the values of underlying variables, which affect the genotypic value of phenotypic characters proportional to a "scaling factor". Epistasis is the influence of one gene on the average effect of another gene. In this model, epistasis is introduced as a mutational effect of one gene on the scaling factors of another gene. In accordance with empirical results, the model assumes that the average direct effect of mutations is zero, as is the average epistatic effect. The model predicts that, on average, a mutation at one locus increases the expected mutational variance of mutations at another interacting locus. The increase in mutational variance is predicted to be equal to the variance of the pair-wise epistatic effects. This result is consistent with the observation that mutant phenotypes tend to be more variable than the wildtype phenotype. Another generic result of this model is that the frequency of canalizing mutations can at most be equal to the frequency of de-canalizing mutations. Furthermore, it is predicted that the mutational variance of a character increases at least linearly with the size of the character; hence this model is scale variant. In the case of two characters it is shown that the dimensionality of the locus-specific mutational effect distribution is invariant, i.e. the rank of the mutational covariance matrix M is invariant. While in additive models the mutational covariance matrix is always and entirely invariant, the invariance in the case of epistatic models is unexpected. Epistatic interactions can change the magnitude of the mutational (co)variances at a locus and can thus influence the structure of the mutational covariance matrix. However, in the present model the dimensionality of the mutational effect distribution remains the same. A consequence of this result is that, in this model, the genetic architecture of a set of characters is always evolvable i.e. no hard constraints can evolve.

A recombinogenic targeting method to modify large-inserts for cis-regulatory analysis in transgenic mice: construction and expression of a 100-kb, zebrafish Hoxa-11b-lacZ reporter gene

The identification of cis-sequences responsible for spatiotemporal patterns of gene expression often requires the functional analysis of large genomic regions. In this study a 100-kb zebrafish Hoxa-11b-lacZ reporter gene was constructed and expressed in transgenic mice. PAC clone 10-O19, containing a portion of the zebrafish HoxA-b cluster, was captured into the yeast-bacterial shuttle vector, pPAC-ResQ, by recombinogenic targeting. A lacZ reporter gene was then inserted in-frame into exon 1 of the zfHoxa-11b locus by a second round of recombinogenic targeting. Expression of the zfHoxa-11b-lacZ reporter gene in 10.5 d.p.f. transgenic mouse embryos was observed only in the posterior portion of the A-P axis, in the paraxial mesoderm, neural tube, and somites. These findings demonstrate the utility of recombinogenic targeting for the modification and expression of large inserts captured from P1/PAC clones.

Asymmetry of configuration space induced by unequal crossover: implications for a mathematical theory of evolutionary innovation

Evolution can be regarded as the exploration of genetic or morphological state space by populations. In traditional models of population and quantitative genetics, the state space can be formally represented as a configuration space with clearly defined concepts of neighborhood and distance, defined by the action of variational operators such as mutation and/or recombination. In this paper, we describe a process where no genetic configuration space closure (and hence, no non-arbitrary notion of distance and neighborhood) exists. The process is gene duplication by means of unequal crossover, which we regard as an example of an "innovation" process that changes the state space of the system rather than exploring a closed state space. We assert that such processes are qualitatively distinct from representations of the adaptation process, which occur on regular configuration spaces.

Duplication of modules facilitates the evolution of functional specialization

The evolution of simulated robots with three different architectures is studied in this article. We compare a nonmodular feed-forward network, a hardwired modular, and a duplication-based modular motor control network. We conclude that both modular architectures outperform the non-modular architecture, both in terms of rate of adaptation as well as the level of adaptation achieved. The main difference between the hardwired and duplication-based modular architectures is that in the latter the modules reached a much higher degree of functional specialization of their motor control units with regard to high-level behavioral functions. The hardwired architectures reach the same level of performance, but have a more distributed assignment of functional tasks to the motor control units. We conclude that the mechanism through which functional specialization is achieved is similar to the mechanism proposed for the evolution of duplicated genes. It is found that the duplication of multifunctional modules first leads to a change in the regulation of the module, leading to a differentiation of the functional context in which the module is used. Then the module adapts to the new functional context. After this second step the system is locked into a functionally specialized state. We suggest that functional specialization may be an evolutionary absorption state.

Is Hsp90 a regulator of evolvability?

In a recent paper, Rutherford and Lindquist (1998. Nature 396:336-342) identified mutations in the Hsp90 protein that act to unmask hidden genetic variation with a variety of phenotypic effects. The Hsp90 protein has a number of properties that suggest a role in regulating the expression of genetic variation and therefore in adjusting the evolvability of the organism. In this paper we reflect upon the evolutionary feasibility of such mechanisms and suggest some possible ways of testing the adaptation-for-evolvability hypothesis in more detail. We conclude that Hsp90 holds promise as a molecular model system for the evolution of evolvability. J. Exp. Zool. ( Mol. Dev. Evol. ) 285:116-118, 1999.

A research programme for testing the biological homology concept

The classical homology concept has served as a heuristic principle for organizing the enormous wealth of information on comparative anatomical patterns across a wide range of organisms. However, the classical homology concept reaches its limit as knowledge of the evolutionary, genetic and developmental processes that underlie these anatomical patterns increases. The biological homology concept places the known anatomical patterns into a mechanistic context and asserts that character identity is based on common variational properties in this chapter a research programme for testing the biological homology concept that involves the following steps is outlined: (1) identifying of two or more putative homologues is a clade; (2) determining the phylogenetic distribution of the putative homologues; (3) describing the intra- and interspecific variation patterns of each putative homologue; (4) describing the development of each putative homologue, and determining if modes of development and distribution of homologues are phylogenetically congruent; and (5) providing and testing a model of how differences in modes of development between putative homologues effect differences in variational tendencies. The goal is to demonstrate a link between developmental and variational differences of two homologues.

Evolution of chordate hox gene clusters

In this article, we consider the role of the Hox genes in chordate and vertebrate evolution from the viewpoints of molecular and developmental evolution. Models of Hox cluster duplication are considered with emphasis on a threefold duplication model. We also show that cluster duplication is consistent with a semiconservative model of duplication, where following duplication, one daughter cluster remains unmodified, while the other diverges and assumes a new architecture and presumably new functions. Evidence is reviewed, suggesting that Hox gene enhancers have played an important role in body plan evolution. Finally, we contrast the invertebrates and vertebrates in terms of genome and Hox cluster duplication which are present in the latter, but not the former. We question whether gene duplication has been important in vertebrates for the introduction of novel features such as limbs, a urogenital system, and specialized neuromuscular interactions.

1,2,3 = 2,3,4: a solution to the problem of the homology of the digits in the avian hand

Persistent contradictions in well supported empirical findings usually point to important scientific problems and may even lead to exciting new insights. One of the most enduring problems in evolutionary biology is the apparent conflict between paleontological and embryological evidence regarding the homology of the digits in the avian hand (1, 2). We propose that this problem highlights an important feature of morphological change: namely, the possible dissociation between the developmental origin of a particular repeated element and its subsequent individualization into a fully functional character. We argue that, although comparative embryological evidence correctly identifies the homology of the primordial condensations in avians as CII, CIII, and CIV, subsequent anatomical differentiation reflects a frame shift in the developmental identities of the avian digit anlagen in later ontogeny such that CII becomes DI, CIII becomes DII, and CIV becomes DIII.

Genetic measurement of theory of epistatic effects

Epistasis is defined as the influence of the genotype at one locus on the effect of a mutation at another locus. As such it plays a crucial role in a variety of evolutionary phenomena such as speciation, population bottle necks, and the evolution of genetic architecture (i.e., the evolution of dominance, canalization, and genetic correlations). In mathematical population genetics, however, epistasis is often represented as a mere noise term in an additive model of gene effects. In this paper it is argued that epistasis needs to be scaled in a way that is more directly related to the mechanisms of evolutionary change. A review of general measurement theory shows that the scaling of a quantitative concept has to reflect the empirical relationships among the objects. To apply these ideas to epistatic mutation effects, it is proposed to scale A x A epistatic effects as the change in the magnitude of the additive effect of a mutation at one locus due to a mutation at a second locus. It is shown that the absolute change in the additive effect at locus A due to a substitution at locus B is always identical to the absolute change in B due to the substitution at A. The absolute A x A epistatic effects of A on B and of B on A are identical, even if the relative effects can be different. The proposed scaling of A x A epistasis leads to particularly simple equations for the decomposition of genotypic variance. The Kacser Burns model of metabolic flux is analyzed for the presence of epistatic effects on flux. It is shown that the non-linearity of the Kacser Burns model is not sufficient to cause A x A epistasis among the genes coding for enzymes. It is concluded that non-linearity of the genotype-phenotype map is not sufficient to cause epistasis. Finally, it is shown that there exist correlations among the additive and epistatic effects among pairs of loci, caused by the inherent symmetries of Mendelian genetic systems. For instance, it is shown that a mutation that has a larger than average additive effect will tend to decrease the additive effect of a second mutation, i.e., it will tend to have a negative (canalizing) interaction with a subsequent gene substitution. This is confirmed in a preliminary analysis of QTL-data for adult body weight in mice.

Heterochronic differences of Hoxa-11 expression in Xenopus fore- and hind limb development: evidence for lower limb identity of the anuran ankle bones

The wrist (carpus) and ankle (tarsus) of most tetrapods, as well as the wrist of anurans, contains relatively small nodular skeletal elements. The anuran tarsus, however, comprises a pair of long bones, the proximal tarsals tibiale and fibulare, which resemble the lower leg bones, tibia and fibula (zeugopodium). In this paper we investigate whether the proximal tarsals of Xenopus are of zeugopodial character identity, i.e. whether they develop under the influence of the same genes that pattern the lower limb. We compare Hoxa-11 expression in the forelimb bud with that in the hind limb bud by whole-mount in situ hybridization. Hoxa-11 has been implicated in the development of the lower limb. In Xenopus we note three differences between Hoxa-11 expression in fore- and hind limb buds: (1) Hoxa-11 expression is maintained until the hind limb bud reaches a larger size (2 mm) than that of the forelimb bud (1.5 mm); (2) Hoxa-11 expression is maintained over larger spatial domains than in the forelimb; and (3) Hoxa-11 expression has a pronounced posterior polarity in the hind limb, but not in the forelimb. Hind limb expression of Hoxa-11 can be understood as a heterochronic prolonging of the expression dynamic in the forelimb. Finally we found that the proximal tarsals start to develop within the expression domain of Hoxa-11, while in the forelimb the lower arm elements reach the distal expression limit of Hoxa-11. The gene expression data presented here support the notion of a zeugopodial identity of the proximal tarsal elements in Xenopus.

Phylogenetic reconstruction of vertebrate Hox cluster duplications

In vertebrates and the cephalochordate, amphioxus, the closest vertebrate relative, Hox genes are linked in a single cluster. Accompanying the emergence of higher vertebrates, the Hox gene cluster duplicated in either a single step or multiple steps, resulting in the four-cluster state present in teleosts and tetrapods. Mammalian Hox clusters (designated A, B, C, and D) extend over 100 kb and are located on four different chromosomes. Reconstructing the history of the duplications and its relation to vertebrate evolution has been problematic due to the lack of alignable sequence information. In this study, the problem was approached by conducting a statistical analysis of sequences from the fibrillar-type collagens (I, II, III, and IV), genes closely linked to each Hox cluster which likely share the same duplication history as the Hox genes. We find statistical support for the hypothesis that the cluster duplication occurred as multiple distinct events and that the four-cluster situation arose by a three-step sequential process.

Algebraic theory of recombination spaces

A new mathematical representation is proposed for the configuration space structure induced by recombination, which we call "P-structure." It consists of a mapping of pairs of objects to the power set of all objects in the search space. The mapping assigns to each pair of parental "genotypes" the set of all recombinant genotypes obtainable from the parental ones. It is shown that this construction allows a Fourier decomposition of fitness landscapes into a superposition of "elementary landscapes." This decomposition is analogous to the Fourier decomposition of fitness landscapes on mutation spaces. The elementary landscapes are obtained as eigenfunctions of a Laplacian operator defined for P-structures. For binary string recombination, the elementary landscapes are exactly the p-spin functions (Walsh functions), that is, the same as the elementary landscapes of the string point mutation spaces (i.e., the hypercube). This supports the notion of a strong homomorphism between string mutation and recombination spaces. However, the effective nearest neighbor correlations on these elementary landscapes differ between mutation and recombination and among different recombination operators. On average, the nearest neighbor correlation is higher for one-point recombination than for uniform recombination. For one-point recombination, the correlations are higher for elementary landscapes with fewer interacting sites as well as for sites that have closer linkage, confirming the qualitative predictions of the Schema Theorem. We conclude that the algebraic approach to fitness landscape analysis can be extended to recombination spaces and provides an effective way to analyze the relative hardness of a landscape for a given recombination operator.

Evidence for four Hox clusters in the killifish Fundulus heteroclitus (teleostei)

We analyzed the HOM/Hox cluster composition of the teleost Fundulus heteroclitus (Killifish) using a PCR survey. We found a total number of 30 unique homeobox sequences which could be assigned to specific cognate groups of the Hox clusters by sequence comparisons. One sequence was of the msh group. Ten homeobox fragments could be identified as orthologs with specific mammalian Hox genes of cluster A-D. The number of representatives in cognate groups 1, 4, 5, 6, 7, 8 and 9 differed from those of human and mouse. For groups 1 and 9 we found four representatives, which provides the first evidence for possible four Hox clusters in a diploid teleost, a lower vertebrate. Furthermore, it demonstrates significant differences of the Fundulus clusters compared to those of mouse and humans. The implications of our data for the interpretation of the evolution of Hox clusters are discussed.

PCR-survey of Hox-genes of the zebrafish: new sequence information and evolutionary implications

We analyzed the Hox gene complement of the zebrafish Danio rerio using a PCR survey. We found 18 new zebrafish HOM/Hox type sequences and one sequence of the msh group. For groups 1-3 and 8-10 we could unambiguously assign the zebrafish fragments to cognate groups. The assignment for cognate groups 4-7 had to remain tentative due to insufficient sequence variation. The number of zebrafish Hox fragments classified as members of cognate groups 1-4, 8, and 9 is identical to the number of genes in corresponding cognate groups of the mouse and human genomes. We found only two differences between the zebrafish and mouse Hox gene complement: four putative genes in group 10 (three in mammals) and only seven in the medial groups 5 to 7 (eight in mammals). Together with the previously published Hox gene sequences of the killifish, the larger number of zebrafish genes in group 10 is positive evidence for variation in the Hox gene complements among bony fish. In contrast, the Hox gene complement appears to be highly conserved among all tetrapods.

Epistasis can facilitate the evolution of reproductive isolation by peak shifts: a two-locus two-allele model

The influence of epistasis on the evolution of reproductive isolation by peak shifts is studied in a two-locus two-allele model of a quantitative genetic character under stabilizing selection. Epistasis is introduced by a simple multiplicative term in the function that maps gene effects onto genotypic values. In the model with only additive effects on the trait, the probability of a peak shift and the amount of reproductive isolation are always inversely related, i.e., the higher the peak shift rate, the lower the amount of reproductive isolation caused by the peak shift. With epistatic characters there is no consistent relationship between these two values. Interestingly, there are cases where both transition rates as well as the amount of reproductive isolation are increased relative to the additive model. This effect has two main causes: a shift in the location of the transition point, and the hybrids between the two alternative optimal genotypes have lower average fitness in the epistatic case. A review of the empirical literature shows that the fitness relations resulting in higher peak shift rates and more reproductive isolation are qualitatively the same as those observed for genes causing hybrid inferiority.

Evolution and multi-functionality of the chitin system

Chitin, that is, the beta-1, 4 linked polysaccharide of N-acetylglucosamine, is best known as a cell wall component of fungi and as skeletal material of invertebrates. In recent years this simple picture has changed dramatically. Three developments have taken place: the discovery of chitinous tissues in vertebrates, the molecular analysis of the chitin-synthase genes, and the discovery that chitin derivatives play a crucial role in the interaction between higher plants and symbiotic bacteria. In this paper the methods for chitin detection and the current data on the evolution of chitin synthesis are reviewed. In addition, data is summarized which suggest that chitin synthesis may serve roles other than the production of skeletal material. In particular, anecdotal evidence suggests that chitin derivatives may play a role as signals in plant as well as animal development. Two major unresolved questions are identified: 1) Is there historical continuity of all the chitin synthesizing systems in protists, animals and, in particular, the deuterostome animals. 2) Are chitin derivatives actually involved in the development of plants and animals?

Regeneration in Salaria pavo (Blenniidae, Teleostei). Histogenesis of the regenerating pectoral fin suggests different mechanisms for morphogenesis and structural maintenance

The pectoral fin of blennies is differentiated into a dorsal field and a ventral hook field. A histogenetic analysis of the regenerating pectoral fin was related to two questions. First, are histological specializations of the hook field responsible for the impairment of the regenerative capacity of pectoral fins of blennies? Second, can analysis of the temporal sequence of histogenetic events be used to make testable predictions about the tissue interactions required to re-establish the adult pattern? Regeneration of pectoral fins was examined in Salaria pavo (Blenniidae, Teleostei). Approximately 80% of the length of the fin rays was amputated. Fin ray stumps were evaluated 7, 14, 24, 48 and 72 h after amputation, regenerates 4, 5, and 6 days after amputation and at length of about 30%, 50% and 60% regeneration of the original fin length. The regeneration process is subdivided into four stages: wound healing, blastema formation, fin ray formation and distal outgrowth and differentiation of hook characters. Analysis of the early events of regeneration, wound healing, blastema formation and distal outgrowth, yielded no profound differences from those of conventional fins in general. Impairment of regenerative capacity becomes manifested before histological differentiation of hook characters, and it is thus unlikely that their presence is the proximate cause of heteromorphic regeneration. The sequence in which the anatomical specializations characteristic of fin hooks (lepidotrichal cord, cuticle, fin web regression) appear was variable. Detailed analysis of older regenerates revealed a more regular pattern. In the first phase the characters appear to be largely independently organized, while they become locally correlated later. It is concluded that the anatomical differentiation passes through two stages, initiation of anatomical differentiation, and then mutual adjustment of character expression leading to spatially correlated expression of the lepidotrichal cord, the cuticle and the fin web regression.

Evolutionary modification of regenerative capability in vertebrates: a comparative study on teleost pectoral fin regeneration

The regenerative ability of the pectoral fins of 14 species from 6 euteleostean families was tested. Blastema formation and distal outgrowth was observed in all species, indicating the initiation of regeneration in all species tested. Interspecific variation exists with respect to the frequency of malformations and the patterns produced by heteromorphic regeneration. Taking into account published reports on pectoral fin regeneration, the systematic distribution of homo- and heteromorphic regeneration leads to the following conclusions: 1) regenerative ability of pectoral fins is a property inherited from the common ancestor of euteleosteans. Whether it is also the ancestral condition for the whole teleostean group cannot be determined, because reports on more primitive teleosteans like the herring and the osteoglossimorphs are missing. 2) A propensity to produce high frequencies of heteromorphic regenerates originated independently at least three times in Cypriniformes, Scorpaeniformes, and Perciformes. 3) Impaired regeneration is most commonly found in bottom fishes, although not all ground fish groups show heteromorphic regeneration. This suggests that impaired regeneration is not directly related to bottom dwelling, but most probably originated as a side effect of other adaptive changes. Hence, neither the presence nor the loss of faithful regeneration can be associated with particular adaptive scenarios in this group, since regeneration seems to be ancestral to all major euteleost groups and its loss has no clear adaptive significance. Whether there are adaptive reasons to maintain regenerative capability or whether there are cases of reestablishment of regeneration after it was lost cannot be decided on the basis of recent evidence. More observations on phylogenetically closely related species with variable regenerative capability are necessary to assess adaptive explanations of regeneration.

Multivariate mutation-selection balance with constrained pleiotropic effects

A multivariate quantitative genetic model is analyzed that is based on the assumption that the genetic variation at a locus j primarily influences an underlying physiological variable yj, while influence on the genotypic values is determined by a kind of "developmental function" which is not changed by mutations at this locus. Assuming additivity among loci the developmental function becomes a linear transformation of the underlying variables y onto the genotypic values x, x = By. In this way the pleiotropic effects become constrained by the structure of the B-matrix. The equilibrium variance under mutation-stabilizing selection balance in infinite and finite populations is derived by using the house of cards approximation. The results are compared to the predictions given by M. Turelli in 1985 for pleiotropic two-character models. It is shown that the B-matrix model gives the same results as Turelli's five-allele model, suggesting that the crucial factor determining the equilibrium variance in multivariate models with pleiotropy is the assumption about constraints on the pleiotropic effects, and not the number of alleles as proposed by Turelli. Finally it is shown that under Gaussian stabilizing selection the structure of the B-matrix has effectively no influence on the mean equilibrium fitness of an infinite population. Hence the B-matrix and consequently to some extent also the structure of the genetic correlation matrix is an almost neutral character. The consequences for the evolution of genetic covariance matrices are discussed.

Tangential organization of the infragranular fiber plexus in rat cerebral cortex

Cylindrical lesions (diameter 300-500 microns) were formed by poking needles into various parts of the cerebral cortex of adult albino rats. Degenerating axons were visualized in horizontal sections through the 'flattened' cortex using the silver impregnation method of Gallyas et al. [Stain Technol. 55: 291-297 (1980)] which stains degenerating axoplasm. The density and distribution of tangentially oriented axons were evaluated in the infragranular layers by TV image analysis. The sampling fields were concentrically arranged around the lesion at distances of 200, 400, 700 and 1,100 microns. The results indicate that the distribution patterns of degenerating (associational) axons covary with the cytoarchitectonic regions into which the lesions were placed. In the motor cortex, the majority of axons run in the antero-posterior direction. The density is generally lower around lesions in frontal regions than in parietal regions. The most extended degeneration was found around lesions near the border of or within the retrosplenial cortex, indicating an exceptionally strong internal connectivity in this area. Since only few degenerating axons were seen around lesions in the center of area 17, the high density of myelinated axons in the primary visual cortex seems to be due to fibers that originate in peristrate areas. It is concluded that the number and extension of fibers that degenerate tends to covary with some aspects of cortical architecture, but it is not area-specific.

Small lesions in the primary visual cortex of rats cause a specific reorganization of associational connections

Neuroplastic changes in associational connections were investigated 3 weeks after the intrinsic organization of the visual cortex of rats had been partially damaged by small cylindrical lesions (type I). These lesions caused the degeneration of short intracortical connections and associational connections that form patches in the primary and secondary visual areas. The resulting terminal degeneration disappeared within 20 days p.o. after which only some fiber degeneration was evident in the infragranular layers. Patches of terminal degeneration reappeared in the vicinity of the stab wounds, when the associational connections between the retrosplenial and the primary visual cortex had been secondarily interrupted by elongated lesions (type II), which penetrated the paramedian cortex and subcortical white matter. When type-II lesions were made in the intact cortex, patches of degeneration were absent, although in both cases some terminal degeneration was diffusely distributed in the primary visual cortex. Horseradish peroxidase (HRP) was applied to sites similar to those where type-I lesions were applied. In the intact cortex, HRP caused a granular labeling of numerous neurons in various positions including the retrosplenial cortex and patches of the postero-median visual cortex. HRP was also applied to type-I lesions that had been made 3 weeks earlier. In these cases, apparently HRP labeled the same subpopulations of neurons as it did in the intact cortex. However, a fraction of the labeled neurons showed a Golgi-like staining (e.g., 27% of the labeled neurons in the retrosplenial cortex) only when HRP was applied to stab wounds.(ABSTRACT TRUNCATED AT 250 WORDS)