Evolutionary Genetics of the Cavefish Astyanax mexicanus

Blind and depigmented fish belonging to the species Astyanax mexicanus are outstanding models for evolutionary genetics. During their evolution in the darkness of caves, they have undergone a number of changes at the morphological, physiological, and behavioral levels, but they can still breed with their river-dwelling conspecifics. The fertile hybrids between these two morphotypes allow forward genetic approaches, from the search of quantitative trait loci to the identification of the mutations underlying the evolution of troglomorphism. We review here the past 30years of evolutionary genetics on Astyanax: from the first crosses and the discovery of convergent evolution of different Astyanax cavefish populations to the most recent evolutionary transcriptomics and genomics studies that have provided researchers with potential candidate genes to be tested using functional genetic approaches. Although significant progress has been made and some genes have been identified, cavefish have not yet fully revealed the secret of their adaptation to the absence of light. In particular, the genetic determinism of their loss of eyes seems complex and still puzzles researchers. We also discuss future research directions, including searches for the origin of cave alleles and searches for selection genome-wide, as well as the necessary but missing information on the timing of cave colonization by surface fish.

Heterotachy, an important process of protein evolution

Because of functional constraints, substitution rates vary among the positions of a protein but are usually assumed to be constant at a given site during evolution. The distribution of the rates across the sequence positions generally fits a Gamma distribution. Models of sequence evolution were accordingly designed and led to improved phylogenetic reconstruction. However, it has been convincingly demonstrated that the evolutionary rate of a given position is not always constant throughout time. We called such within-site rate variations heterotachy (for "different speed" in Greek). Yet, heterotachy was found among homologous sequences of distantly related organisms, often with different functions. In such cases, the functional constraints are likely different, which would explain the different distribution of variable sites. To evaluate the importance of heterotachy, we focused on amino acid sequences of mitochondrial cytochrome b, for which the function is likely the same in all vertebrates. Using 2,038 sequences, we demonstrate that 95% of the variable positions are heterotachous, i.e., underwent dramatic variations of substitution rate among vertebrate lineages. Heterotachy even occurs at small evolutionary scale, and in these cases it is very unlikely to be related to functional changes. Since a large number of sequences are required to efficiently detect heterotachy, the extent of this phenomenon could not be estimated for all proteins yet. It could be as large as for cytochrome b, since this protein is not a peculiar case. The observations made here open several new avenues of research, such as the understanding of the evolution of functional constraints or the improvement of phylogenetic reconstruction methods.

Molecular evidence for precambrian origin of amelogenin, the major protein of vertebrate enamel

Although molecular dating of cladogenetic events is possible, no molecular method has been described to date the acquisition of various tissues. Taking into account the specificity of the major protein in enamel in formation (amelogenin), we were able to develop such a method for enamel. Indeed, because the amelogenin protein is exclusively involved in enamel formation and mineralization and because it lacks pleiotropic effects, this protein is a good candidate to estimate the date of acquisition of this highly mineralized tissue. We searched DNA banks for similarities between the amelogenin sequence and other sequences. Similarities were found only to exon 2 of SPARC (osteonectin) in two protostomians and in eight deuterostomians, and to exon 2 of three SPARC-related deuterostomian genes (SC1, hevin, and QR1). The other amelogenin exons did not reveal significant similarities to other sequences. In these proteins, exon 2 mainly encodes the peptide signal that plays the essential role in enabling the protein to be ultimately localized in the extracellular matrix. We tested the significance of the exon 2 similarities. The observed values were always significantly higher than the expected randomly generated similarities. This demonstrates a common evolutionary origin of this exon. The phylogenetic analyses of exon 2 sequences indicated that exon 2 was duplicated to amelogenin from an ancestral SPARC sequence in the deuterostomian lineage before the duplication of deuterostomian SPARC and SC1/hevin/QR1. We were able to date the origin of the latter duplication at approximately 630 MYA. Therefore, amelogenin exon 2 was acquired before this date, in the Proterozoic, long before the so-called "Cambrian explosion," the sudden appearance of several bilateralian phyla in the fossil record at the Proterozoic-Phanerozoic transition. This sudden appearance has been often suggested to reflect intensive cladogenesis during this period. However, molecular dating of protostomian-deuterostomian divergence and of the cladogenesis among several major clades of Bilateralia lead to a different conclusion: many bilateralian clades were already present during the late Proterozoic. It has previously been proposed that these bilateralians were not mineralized and that they had low fossilization potential. Our results strongly suggest that late Proterozoic fossils possessing a mineralized tissue homologous to enamel might be found in the future.

evx1 transcription in bony fin rays segment boundaries leads to a reiterated pattern during zebrafish fin development and regeneration

The dermoskeleton of zebrafish fins is composed of actinotrichia and segmented bony rays, or lepidotrichia, which grow by successive addition of distal segments. The present study shows that evx1, a new zebrafish even-skipped related gene (Thaëron et al., 2000) displays during bony ray morphogenesis, a unique repetitive expression pattern along the proximodistal axis of the fin. Whole-mount in situ hybridization performed on larvae and adult regenerating fins show that evx1 signal appears as parallel dash lines crossing the width of each developing and regenerating rays, in a ladder-like fashion. Cytological studies show that a subpopulation of bone forming cells (scleroblasts) expresses evx1 at the level of the joint between two adjacent segments except in the apical part of the differentiating ray where evx1 expression precedes the formation of the joint. This distal transcription is turned on again only when the latest differentiating segment reached its final size and might label the putative next segment boundary. This suggests the existence of a molecular mechanism controlling the periodic expression of evx1 which could be involved in the establishment of segment boundaries during fin ray morphogenesis, and could play a key role during dermal skeleton patterning.

Zebrafish evx1 is dynamically expressed during embryogenesis in subsets of interneurones, posterior gut and urogenital system

The even-skipped-related homeobox genes (evx) are widely distributed through animal kingdom and are thought to play key role in posterior body patterning and neurogenesis. We have cloned and analyzed the expression of evx1 in zebrafish (see also Borday et al. (Dev. Dyn. 220 (2001) in press) which displays a dynamic and restricted expression pattern during neurogenesis. In spinal cord, rhombencephalon, and epiphysis, evx1 is expressed in several subsets of emerging interneurones prior to their axonal outgrowth, identified as primary interneurones and a subset of Pax2.1(+) commissural interneurones. In the hindbrain, evx1 is expressed in reticulospinal interneurones of rhombomeres 5 and 6 as well as in rhombomere 7 interneurones. The latest emerging evx1(+) interneurones in the hindbrain correspond to commissural interneurones. evx1 is also dynamically transcribed during the formation of the posterior gut and the uro-genital system in mesenchymal cells that border the pronephric ducts, the wall of the pronephric duct, and later in the posterior gut and the wall of the uro-genital opening. In larvae, the ano-rectal epithelium and the muscular layer that surrounds the analia-genitalia region remain stained up to 27 days. In contrast other vertebrates, evx1displays no early nor caudal expression in zebrafish.

Human-chimpanzee DNA sequence variation in the four major genes of the renin angiotensin system

The renin angiotensin system (RAS) is involved in blood pressure control and water/sodium metabolism. The genes encoding the proteins of this system are candidate genes for essential hypertension. The RAS involves four main molecules: angiotensinogen, renin, angiotensin I-converting enzyme, and the angiotensin II type 1 receptor (encoded by the genes AGT, REN, DCP1, and AGTR1, respectively). We performed a molecular screening over 17,037 bp of the coding and 5' and 3' untranslated regions of these genes, from three to six common chimpanzees. We identified 44 single-nucleotide polymorphisms (SNPs) in chimpanzee samples, including 18 coding-region SNPs, 5 of which led to an amino acid replacement. We observed common and different features at various sites (synonymous, nonsynonymous, and noncoding) within and between the four chimpanzee genes: (1) the nucleotide diversity at noncoding sites was similar; (2) the nucleotide diversity at nonsynonymous sites was low, probably reflecting purifying selection, except for the AGT gene; (3) the nucleotide diversity at synonymous sites, which was dependent on the G+C content at the third position of the codon, was high, except for the AGTR1 gene. Comparison of the chimpanzee SNPs with those previously reported for humans identified 119 sites with fixed differences (including 62 coding sites, 17 of which resulted in amino acid differences between the species). Analysis of polymorphism within species and divergence between species shed light on the evolutionary constraints on these genes. In particular, comparison of the pattern of mutation at polymorphic and fixed sites between humans and chimpanzees suggested that the high G+C content of the DCP1 gene was maintained by positive selection at its silent sites. Finally, we propose 68 ancestral alleles for the human RAS genes and discuss the implications for their use in future hypertension-susceptibility association studies.

Mutation pattern variation among regions of the primate genome

We sequenced three argininosuccinate-synthetase-processed pseudogenes (PsiAS-A1, PsiAS-A3, PsiAS-3) and their noncoding flanking sequences in human, orangutan, baboon, and colobus. Our data showed that these pseudogenes were incorporated into the genome of the Old World monkeys after the divergence of the Old World and New World monkey lineages. These pseudogene flanking regions show variable mutation rates and patterns. The variation in the G/C to A/T mutation rate (u) can account for the unequal GC contents at equilibrium: 34.9, 36.9, and 41.7% in the pseudogene PsiAS-A1, PsiAS-A3, and PsiAS-3 flanking regions, respectively. The A/T to G/C mutation rate (v) seems stable and the u/v ratios equal 1.9, 1.7, and 1.4 in the flanking regions of PsiAS-A1, PsiAS-A3, and PsiAS-3, respectively. These "regional" variations of the mutation rate affect the evolution of the pseudogenes, too. The ratio u/v being greater than 1.0 in each case, the overall mutation rate in the GC-rich pseudogenes is, as expected, higher than in their GC-poor flanking regions. Moreover, a "sequence effect" has been found. In the three cases examined u and v are higher (at least 20%) in the pseudogene than in its flanking region-i.e., the pseudogene appears as mutation "hot" spots embedded in "cold" regions. This observation could be partly linked to the fact that the pseudogene flanking regions are long-standing unconstrained DNA sequences, whereas the pseudogenes were relieved of selection on their coding functions only around 30-40 million years ago. We suspect that relatively more mutable sites maintained unchanged during the evolution of the argininosuccinate gene are able to change in the pseudogenes, such sites being eliminated or rare in the flanking regions which have been void of strong selective constraints over a much longer period. Our results shed light on (1) the multiplicity of factors that tune the spontaneous mutation rate and (2) the impact of the genomic position of a sequence on its evolution.

Nonneutral evolution of tandem repeats in the mitochondrial DNA control region of lagomorphs

The mitochondrial DNA of the European rabbit (Oryctolagus cuniculus) contains a tandem array of 153-bp repeats in the vicinity of the replication origin of the H-stand. Variation among molecules in the number of these repeats results in inter- and intraindividual length polymorphism (heteroplasmy). Generally, in an individual, one predominant molecular type is observed, the others representing a low percentage of the mtDNA content. At the tissue level, we observe a particular distribution of this polymorphism in the gonads compared with liver, kidneys, or brain, implying a relationship between the differentiation status of the cells and the types of new mtDNA molecules which appear and accumulate during lifetime. Similar tandem repeats were also found in the mtDNA noncoding region of European hares (Lepus europaeus), a cottontail (Sylvilagus floridanus), and a pika (Ochotona rufescens). The lengths and the sequences of these units evolve rapidly and in a concerted way, but the number of repeats is maintained in a narrow range, and an internal 20-bp segment is highly conserved. Constraints restrict the evolution of the primary sequence of these repeated units, the number of which is probably controlled by a stabilizing selection.

Rabbit mitochondrial DNA diversity from prehistoric to modern times

The mitochondrial genetic variability in European rabbit (Oryctolagus cuniculus) populations present in Europe and North Africa from 11,000 years ago to the present day has been analyzed using ancient DNA techniques. DNA was extracted from 90 rabbit bones found in 22 archaeological sites dated between the Mesolithic and recent times. Nucleotide sequences present in a variable 233-bp domain of the cytochrome b gene were compared to those present in modern-day rabbits. The results show that the structure of ancient populations of wild rabbit exhibited remarkable stability over time until the Middle Ages. At this time, a novel type of mtDNA molecule abruptly appears into most wild populations studied from France. This mtDNA type corresponds to that currently present in the domestic breeds of rabbit examined so far. The relative rapidity by which this mtDNA type established and its absence in all sites examined before 1,700 years ago lend support to the hypothesis that between 2,000 and 1,000 years ago, man may have favored the development, into all regions of France, of animals carrying this particular mtDNA molecule. The origin of such animals has still to be found: animals previously living outside of France or within France but in very restricted areas? This event was concomitant with the documented establishment of warrens after the tenth century A.D. in Europe.

Genetic analysis of systematic mitochondrial heteroplasmy in rabbits

One unusual property of rabbit mitochondrial DNA (mtDNA) is the existence of repeated 153-bp motifs in the vicinity of the replication origin of its H strand. Furthermore, every individual is heteroplasmic: it carries mtDNA molecules with a variable number of repeats. A systematic study of 8 females and their progeny has been devised to analyze mtDNA transmission through generations. The results suggest that three mechanisms are acting simultaneously. (1) Genetic drift in the germ line is revealed by the evolution of heteroplasmy when two major molecular forms are present in a female. (2) A high mutation rate (around 10(-2) per animal generation) generating molecular diversity, by deletion and addition of repeated units, is required to explain the observation of heteroplasmy in every individual. Moreover, the rates of mutation from the most frequent type to the other types are unequal. The deletion of one unit is more frequent than a deletion of two units, which is in turn more frequent than a deletion of three. (3) Selection for shorter molecules in somatic cells is probable. The frequency distribution of mtDNA types depends on the organ analyzed (kidney-spleen and liver vs. gonads).

Ancient DNA from Bronze Age bones of European rabbit (Oryctolagus cuniculus)

The European rabbit (Oryctolagus cuniculus) is now widely distributed throughout the world as a result of transportation by man. The original populations, however, were confined to southern France and Spain. In order to investigate the role of human intervention in determining the genetic diversity of rabbit populations, we are studying the origin of rabbits introduced onto a small Mediterranean island (Zembra) near Tunis over 1400 years ago, by examining ancient DNA extracted from rabbit bones found both on Zembra and on the European mainland. Ancient DNA was successfully extracted from rabbit bones found at two archaeological sites dated to at least the Early Bronze Age (more than 3500 years ago) in south-central France, and compared to that found in modern mainland and island populations using a small variable region of the cytochrome b gene. The results confirm that the Zembra Island population is descended from that present over 1400 years ago. The technical aspects of DNA extraction from bones and the implications of this type of research for determining the origin of introduced rabbit populations are discussed.