Comparison of Substitution Rates in ZFX and ZFY Introns of Sheep and Goat Related Species Supports the Hypothesis of Male-Biased Mutation Rates

There is a growing body of evidence that males serve as the major generators of mutations, due to the larger number of cell divisions involved in sperm compared to egg production. In mammals, this hypothesis (referred to as ``male-driven evolution') has been tested by comparison of nucleotide substitution rates on the X and Y sex chromosomes in a limited number of taxa, predominantly primates and rodents. This study asks whether male-driven evolution is a more general phenomenon among mammals, by comparison of paralogous ZFX and ZFY intron sequences in sheep and goat species (the tribe Caprini). The male-to-female mutation ratio, α_m, was estimated to be between 2.93 (95% CI, 1.51–8.61) and 3.94 (95% CI, 1.25–32.29) when calculated using pairwise distance and branch length, respectively, suggesting that the Caprini are subject to weak, male-driven evolution. Comparison to published values for primates, felids, and rodents implies that there may be some correlation with reproductive life span. However, this is difficult to test with current data because confidence intervals are large and overlapping. Nonindependent evolution of paralogous sequences and/or the presence of selective constraints could lead to inaccurate estimates of α_m. No evidence for gene conversion between the ZFX and the ZFY introns was found, and this suggests that they have evolved independently during the radiation of the Caprini. Finally, there was no apparent evidence that these introns are subject to selective constraints, although low levels of intraspecific polymorphism reduce the power of neutrality tests. Received: 13 February 2001 / Accepted: 23 May 2001     

The karyotype of Alligator mississippiensis,and chromosomal mapping of the ZFY/X homologue,Zfc

Comparative mapping studies of X-linked genes in mammals have provided insights into the evolution of the X chromosome. Many reptiles including the American alligator, Alligator mississippiensis, do not appear to possess heteromorphic sex chromosomes, and sex is determined by the incubation temperature of the egg during embryonic development. Mapping of homologues of mammalian X-linked genes in reptiles could lead to a greater understanding of the evolution of vertebrate sex chromosomes. One of the genes used in the mammalian mapping studies was ZFX, an X-linked copy of the human ZFY gene which was originally isolated as a candidate for the mammalian testis-determining factor (TDF). ZFX is X-linked in eutherians, but maps to two autosomal locations in marsupials and monotremes, close to other genes associated with the eutherian X. The alligator homologue of the ZFY/ZFX genes, Zfc, has been isolated and described previously. A detailed karyotype of A. mississippiensis is presented, together with chromosomal in situ hybridisation data localising the Zfc gene to chromosome 3. Further chromosomal mapping studies using eutherian X-linked genes may reveal conserved chromosomal regions in the alligator that have become part of the eutherian X chromosome during evolution.     

Cloning and mapping of bovine ZFX gene to the long arm of the X-chromosome (Xq34) and homologous mapping of ZFY gene to the distal region of the short arm of the bovine (Yp13), ovine (Yp12–p13), and caprine (Yp12–p13) Y Chromosome

A part of mouse Zfy-2 sequence was synthesized and used to screen a genomic library of the spinous country-rat (Tokudaia osimensis spp., 2n = 45). An isolated clone had the C-terminal region of Zfy, which consisted of 1190 bp, encoded 336 amino acid residues, and harbored 11 out of 13 zinc finger motifs. With this as a probe, a bovine testis cDNA library was screened. Two ZFX clones were isolated and their sequences combined. The short sequence, lacking part of the 5′ upstream region, was amplified by PCR or RT-PCR, cloned, and sequenced. A full-length ZFX was constructed by combining these three sequences. The bovine ZFX consisted of 5328 bp and encoded 800 amino acid residues, which contained 13 zinc finger motifs. ZFX was used as a probe for fluorescence in situ hybridization and was mapped to Xq34, different from its previously reported site at Xq21-q231. A SINE (short interspersed nuclear element) sequence consisting of 188 bp was found close to the end of the 3′-untranslated region of ZFX. The SINE sequence hybridized to all bovine chromosomes. ZFY is highly homologous with ZFX and, as a result, ZFY could be mapped simultaneously. ZFY was mapped to the distal region of the short arm of the Y Chromosome (Chr) (Yp13), contradicting the previously reported position Yq1. Ovine and caprine ZFY were also mapped with bovine ZFX. Both were mapped to the distal region of the short arm of the Y Chr (Yp12-p13). Ovine ZFX was mapped to a region close to the centromere of the X Chr (Xq13). Received: 23 July 1997 / Accepted: 30 September 1997     

Estimating the intensity of male-driven evolution in rodents by using X-linked and Y-linked Ube 1 genes and pseudogenes

Using sequence data from the last introns of ZFX and ZFY genes, we previously estimated the male-to-female ratio () of mutation rate to be close to 6 in higher primates and 1.8 in rodents. As the mutation rate may vary among different regions of the mammalian genome, it is interesting to see whether sequence data from other regions will give similar estimates. In this study, we have determined the partial genomic sequences of the ubiquitin-activating enzyme El genes (Ube 1x and Ube 1y for the X-linked and Y-linked homologues, respectively) of mice and rats and two mouse Ube 1y pseudogenes. From the intron sequences of the Ube 1 genes, we calculated the divergence of the Y-linked genes (Y = 0.161) and that of the X-linked genes (X = 0.107) between mouse and rat, and found the Y/X ratio to be 1.50. This ratio led to an estimate of = 2.0 with a 95% confidence interval of (1.0, 3.9). Similar estimates of were obtained if mouse Ube 1y pseudogenes were used instead of the mouse Ube 1y functional gene. These estimates are consistent with our previous estimate for rodents and suggest that the sex ratio of mutation rate in rodents is approximately only one-third of that in higher primates. Our estimate of the divergence time between Ube 1x and Ube 1y supports the view that the two genes separated before the eutherian radiation.Correspondence to: W.-H. Li     

Potential problems in estimating the male-to-female mutation rate ratio from DNA sequence data

It is commonly believed that the rate of mutation is much higher in males than in females because the number of germ-cell divisions per generation is much larger in males than in females. However, the precise magnitude of the male-to-female mutation rate ratio (α _m ) remains unknown. Recently there have been efforts to estimate α _m by using DNA sequence data from different species. We have studied the potential problems in such an approach. We found that the rate of synonymous substitution varies about fivefold among X-linked genes, as large as the variation among autosomal genes. This large variation makes the assumption of selective neutrality of synonymous changes dubious, so one should be cautious in using the synonymous rates in X-linked and autosomal genes to estimate α _m . A similar difficulty was also observed in using nonhomologous intron sequences to estimate α _m . Contrary to the expectation that X-linked sequences should evolve more slowly than autosomal sequences, theAlu repeat in the last intron of the X-linked zinc finger gene has evolved faster than the four autosomalAlu repeats used in this study. It appears that the best way to estimate α _m is to use homologous sequences. However, such sequences may be involved in gene conversion events. In fact, we found evidence that the Y-linked and X-linked zinc finger genes have been involved in multiple conversion events during primate evolution. Thus, the possibility of gene conversion should be considered when using homologous sequences to estimate α _m . Presented at the NATO Advanced Research Workshop onGenome Organization and Evolution, Spetsai, Greece, 16–22 September 1992     

AZFY-like sequence in fish,with comments on the evolution of theZFY family of genes in vertebrates

ZFY-like genes have been observed in a variety of vertebrate species. Although originally implicated as the primary testis-determining gene in humans and other placental mammals, more recent evidence indicates a role(s) outside that of testis determination. In this study, DNA from five species of fish,Carasius auratus, Rivulus marmoratus, Xiphophorus maculatus, X. milleri, andX. nigrensis was subjected to Southern blot analysis using a PCR-amplified fragment of mouseZFY-like sequence as a probe. Restriction fragment patterns were not polymorphic between sexes in any one species but showed a different pattern for each species. With one exception,Rivulus, a 3.1-kb band from theEcoRI digestion was common to all. Sequence and open reading frame analysis of this fragment showed a strong homology to other known vertebrateZFY-like genes. Of particular interest in this gene is a novel third finger domain similar to one human and one alligatorZFY-like gene. Our studies and others provide evidence for a family of vertebrateZFY genes, with those having this novel third finger being representative of the ancestral condition.     

Comparative mapping of ZFY in the hominoid apes

Summary Within our project of comparative mapping of candidate genes for sex-determination/testis differentiation, we used a cloned probe from the human ZFY locus for comparative hybridization studies in hominoids. As in the human, the ZFY probe detects X- and Y-specific restriction fragments in the chimpanzee, the gorilla, the orangutan, and the gibbon. Furthermore, the X-specific hybridization site in the great apes resides in Xp21.3, the same locus defining ZFX in the human. The Y-specific locus of ZFY maps closely to the early replicating pseudoautosomal segment in the telomeric or subtelomeric position of the Y chromosomes of the great apes, again as found in the human. Thus, despite cytogenetically visible structural alterations within the euchromatic parts of the Y chromosomes of the human species and the great apes, a segment of the Y chromosome defined by the pseudoautosomal region and ZFY seems to be more strongly conserved than the rest of the Y chromosome.     

The γ-crystallin gene families: Sequence and evolutionary patterns

Summary The -crystallin proteins consist of two topologically equivalent domains, each built up out of two similar motifs. They are encoded by a gene family, which already contained five members before the divergence of rodents and primates. A further gene duplication took place in each lineage. To analyze the pattern of evolution within this gene family, the coding sequences of six human genes, six rat genes, and four mouse genes were compared. Between species, a uniform rate of evolution of all regions of the protein is seen. The ratio of synonymous to nonsynonymous substitution in the human/rat or human/mouse comparison is much lower than the ratio when rat and mouse are compared indicating that the -crystallin proteins are better conserved in the rodent lineage. Within species, the regions encoding the two external motifs I and III of the protein show a greater extent of nonsynonymous substitution than the regions encoding the two internal protein motifs II and IV. The low extent of synonymous substitution between the second exons (encoding motifs I and II) of the rat -crystallin genes suggests the frequent occurrence of gene conversion. In contrast, a high extent of synonymous substitution is found in exon 3 (encoding motifs III and IV) of the rat genes. The same phenomenon is seen within the human gene family. The frequencies of occurrence of the various dinucleotides deviate less from those predicted from the frequencies of occurrence of each individual nucleotide in the second exons than in the third exons. The sequences of the third exons are significantly depleted in CpG, ApA, and GpT and enriched in CpT and GpA.     

Phylogeography and population history of the least weasel (Mustela nivalis) in the Palearctic based on multilocus analysis

The least weasel (Mustela nivalis) is one of the most widely distributed carnivorans. While previous studies have identified distinct western and eastern mitochondrial DNA (mtDNA) lineages of the species in the western Palearctic, their broader distributions across the Palearctic have remained unknown. To address the broad-scale phylogeographical structure, we expanded the sampling to populations in Eastern Europe, the Urals, the Russian Far East, and Japan, and analyzed the mtDNA control region and cytochrome b, the final intron of the zinc finger protein on Y chromosome (ZFY), and the autosomal agouti signaling protein gene (ASIP). The mtDNA data analysis exposed the previous western lineage (Clade I) but poorly supported assemblage extending across Palearctic, whereas the previous eastern lineage (Clade II) was reconfirmed and limited in the south western part of the Palearctic. The ZFY phylogeny showed a distinctive split that corresponding to the mtDNA lineage split, although less phylogeographical structure was seen in the ASIP variation. Our data concur with the previous inference of the Black Sea–Caspian Sea area having an ancestral character. The Urals region harbored high mitochondrial diversity, with an estimated coalescent time of around 100,000 years, suggesting this could have been a cryptic refugium. Based on the coalescent-based demographic reconstructions, the expansion of Clade I across the Palearctic was remarkably rapid, while Clade II was relatively stable for a longer time. It seems that Clade II has maintained a constant population size in the temperate region, and the expansive Clade I represents adaptation to the cold regions.     

Different rates of substitution may produce different phylogenies of the eutherian mammals

Summary In an attempt to resolve some points of branching order in the phylogeny of the eutherian mammals, a phylogenetic analysis of 26 nuclear and 6 mitochondrial genes was undertaken using a maximum likelihood method on a constant rate stochastic model of molecular evolution. Seventeen of the nuclear genes gave a primates/artiodactyls grouping highest support whereas three of the mitochondrial genes found a rodents/artiodactyls grouping to be best supported. The primates/rodents grouping was never the best supported. On the assumption that rodents are indeed an outgroup to primates and artiodactyls and that the latter taxa diverged 70 million years ago, an estimation was made, for each gene, of the time of divergence of the rodent lineage. In most cases such estimates were beyond the limits set by present interpretations of the paleontological record as were many estimates of the divergence time of mouse and rat. These results suggest that, although there is locus variation, the divergent position of the rodent lineage may be an artifact of an elevated rate of nucleotide substitution in this order.     

Analysis of mutation rates in the SMCY/SMCX genes shows that mammalian evolution is male driven

Mammalian evolution is believed to be male driven because the greater number of germ cell divisions per generation in males increases the opportunity for errors in DNA replication. Since the Y Chromosome (Chr) replicates exclusively in males, its genes should also evolve faster than X or autosomal genes. In addition, estimating the overall male-to-female mutation ratio (α_m) is of great importance as a large α_m implies that replication-independent mutagenic events play a relatively small role in evolution. A small α_m suggests that the impact of these factors may, in fact, be significant. In order to address this problem, we have analyzed the rates of evolution in the homologous X-Y common SMCX/SMCY genes from three different species—mouse, human, and horse. The SMC genes were chosen because the X and Y copies are highly homologous, well conserved in evolution, and in all probability functionally interchangeable. Sequence comparisons and analysis of synonymous substitutions in approximately 1kb of the 5′ coding region of the SMC genes reveal that the Y-linked copies are evolving approximately 1.8 times faster than their X homologs. The male-to-female mutation ratio α_m was estimated to be 3. These data support the hypothesis that mammalian evolution is male driven. However, the ratio value is far smaller than suggested in earlier works, implying significance of replication-independent mutagenic events in evolution. Received: 18 April 1996 / Accepted: 4 October 1996     

Interchromosomal gene conversion as a possible mechanism for explaining divergence patterns of ZFY-related genes

Summary The divergence pattern of mammalian ZFY-related genes from human (ZFY and ZFX) and mouse (Zfy-1 and Zfx) was reexamined on the basis of nucleotide substitutions at the synonymous codon-alternating positions. It is possible to explain the unusual divergence pattern of the mammalian Y-linked ZF genes by interchromosomal gene conversion by X-linked ZF genes. Furthermore, the rates of evolution of mammalian X- and Y-linked ZF genes were shown to agree well with those expected from our model. Offprint requests to: T. Miyata     

Molecular sexing assays in 114 mammalian species: In silico sequence reanalysis and a unified graphical visualization of diagnostic tests

Molecular‐based methods for identifying sex in mammals have a wide range of applications, from embryo manipulation to ecological studies. Various sex‐specific or homologous genes can be used for this purpose, PCR amplification being a common method. Over the years, the number of reported tests and the range of tested species have increased greatly. The aim of the present analysis was to retrieve PCR‐based sexing assays for a range of mammalian species, gathering the gene sequences from either the articles or online databases, and visualize the molecular design in a uniform manner. For nucleotide alignment and diagnostic test visualization, the following genomic databases and tools were used: NCBI, Ensembl Nucleotide BLAST, ClustalW2, and NEBcutter V2.0. In the 45 gathered articles, 59 different diagnostic tests based on eight different PCR‐based methods were developed for 114 mammalian species. Most commonly used genes for the analysis were ZFX, ZFY, AMELX, and AMELY. The tests were most commonly based on sex‐specific insertions and deletions (SSIndels) and sex‐specific sequence polymorphisms (SSSP). This review provides an overview of PCR‐based sexing methods developed for mammals. This information will facilitate more efficient development of novel molecular sexing assays and reuse of previously developed tests. Development of many novel and improvement of previously developed tests is also expected with the rapid increase in the quantity and quality of available genetic information.     

<Emphasis Type="Italic">Usp9y</Emphasis> (ubiquitin-specific protease 9 gene on the Y) is associated with a functional promoter and encodes an intact open reading frame homologous to <Emphasis Type="Italic">Usp9x</Emphasis> that is under selective constraint

Sequences complementary to the X-linked ubiquitin-specific protease gene Usp9x (Dffrx) have been shown to map to the Sxr^b interval of the mouse Y Chromosome (chr) and to be expressed in a testis-specific manner. In humans, ubiquitously expressed functional homologues (USP9Y and USP9X DFFRY/DFFRX) are present on both sex chromosomes, whereas in mouse it remains to be demonstrated that the Y-linked sequences encode a functional protein. In this paper, it is shown that the Usp9y gene encodes a potentially functional ubiquitin-specific protease possessing a core promoter region that shares several features characteristic of other testis-specific genes. Analysis of synonymous and nonsynonymous nucleotide changes suggests that there is constraint on the amino acid sequence of both the mouse Usp9x and Usp9y genes, a finding that mirrors similar analysis of the human orthologs. Thus, in both mouse and human, selection is acting to maintain the amino acid sequence of the X and Y-linked genes. This indicates that in both species the genes on each sex chromosome continue to encode an important function.     

Evidence for male-driven evolution in Drosophila

In several vertebrate taxa studied to date, mutation rates arehigher in males than females (male-driven evolution). The male-to-femalemutation rate () can be estimated by contrasting DNA divergencedata at X-linked, Y-linked, and autosomal loci. Previous studiesin Drosophila, comparing X-linked and autosomal divergence,have found no evidence for male-driven evolution in this genus.Here, I compare levels of nucleotide divergence between homologousX- and Y-linked loci in Drosophila miranda. Using divergenceat both synonymous sites and at short introns, I estimate tobe approximately 2. This study thus provides the first evidencefor male-biased mutation rates outside vertebrates, supportingthe view that DNA sequence evolution is male driven in a widevariety of taxa.     

Synonymous substitution rates in Drosophila: Mitochondrial versus nuclear genes

Synonymous substitution rates in mitochondrial and nuclear genes of Drosophila were compared. To make accurate comparisons, we considered the following: (1) relative synonymous rates, which do not require divergence time estimates, should be used; (2) methods estimating divergence should take into account base composition; (3) only very closely related species should be used to avoid effects of saturation; (4) the heterogeneity of rates should be examined. We modified the methods estimating synonymous substitution numbers to account for base composition bias. By using these methods, we found that mitochondrial genes have 1.7–3.4 times higher synonymous substitution rates than the fastest nuclear genes or 4.5–9.0 times higher rates than the average nuclear genes. The average rate of synonymous transversions was 2.7 (estimated from the melanogaster species subgroup) or 2.9 (estimated from the obscura group) times higher in mitochondrial genes than in nuclear genes. Synonymous transversions in mitochondrial genes occurred at an approximately equivalent rate to those in the fastest nuclear genes. This last result is not consistent with the hypothesis that the difference in turnover rates between mitochondrial and nuclear genomes is the major factor determining higher synonymous substitution rates in mtDNA. We conclude that the difference in synonymous substitution rates is due to a combination of two factors: a higher transitional mutation rate in mtDNA and constraints on nuclear genes due to selection for codon usage. Received: 27 November 1996 / Accepted: 8 May 1997     

Sexual differences in the foraging ecology of 19th century beluga whales (Delphinapterus leucas) from the Canadian High Arctic

Marine mammals often exhibit significant sexual segregation in their diet and habitat use but these differences have not been studied systematically in historic or ancient populations due to the difficulties associated with determining the sex of skeletal elements based on gross morphology. Using a combined ancient DNA and stable isotope approach, we document a sexual difference in the foraging ecology of late 19th century beluga whales (Delphinapterus leucas) from the Canadian High Arctic. Using two PCR assays that coamplify fragments of the Y-linked SRY and X-linked ZFX genes, we assigned reproducible sex identities to 35 beluga specimens. This provided a basis for investigating sex-specific differences in foraging ecology using stable carbon and nitrogen isotope analyses of bone collagen. These isotopic data demonstrate that although both males and females primarily consumed Arctic cod, males utilized a wider range of prey than females, feeding on high trophic level benthic prey (sculpins) to a greater extent. Because bone collagen integrates prey isotopic compositions over the course of several years these sex-based differences in beluga bone collagen isotopic compositions reflect long-term and sustained sexual differences in foraging.