Rooting a phylogeny with homologous genes on opposite sex chromosomes (gametologs): a case study using avian CHD

Comparison of the chicken and zebra finch Z chromosomes shows evolutionary rearrangements

Using fluorescent in-situ hybridization (FISH) of zebra finch (Taeniopygia guttata) bacterial artificial chromosome (BAC) clones, we determined the chromosomal localizations of 14 zebra finch genes that are Z-linked in chickens: ATP5A1, CHD1, NR2F1, DMRT1, PAM, GHR, HSD17B4, NIPBL, ACO1, HINT1, SMAD2, SPIN, NTRK2 and UBE2R2. All 14 genes also map to the zebra finch Z chromosome, indicating substantial conservation of gene content on the Z chromosome in the two avian lineages. However, the physical order of these genes on the zebra finch Z chromosome differed from that of the chicken, in a pattern that would have required several inversions since the two lineages diverged. Eight of 14 zebra finch BAC DNA showed cross-hybridization to the W chromosome, usually to the entire W chromosome, suggesting that repetitive sequences are shared by the W and Z chromosomes. These repetitive sequences likely evolved in the finch lineage after it diverged from the Galliform lineage.

Comparison of the Z and W sex chromosomal architectures in elegant crested tinamou (Eudromia elegans) and ostrich (Struthio camelus) and the process of sex chromosome differentiation in palaeognathous birds

To clarify the process of avian sex chromosome differentiation in palaeognathous birds, we performed molecular and cytogenetic characterization of W chromosome-specific repetitive DNA sequences for elegant crested tinamou (Eudromia elegans, Tinamiformes) and constructed comparative cytogenetic maps of the Z and W chromosomes with nine chicken Z-linked gene homologues for E. elegans and ostrich (Struthio camelus, Struthioniformes). A novel family of W-specific repetitive sequences isolated from E. elegans was found to be composed of guanine- and cytosine-rich 293-bp elements that were tandemly arrayed in the genome as satellite DNA. No nucleotide sequence homologies were found for the Struthioniformes and neognathous birds. The comparative cytogenetic maps of the Z and W chromosomes of E. elegans and S. camelus revealed that there are partial deletions in the proximal regions of the W chromosomes in the two species, and the W chromosome is more differentiated in E. elegans than in S. camelus. These results suggest that a deletion firstly occurred in the proximal region close to the centromere of the acrocentric proto-W chromosome and advanced toward the distal region. In E. elegans, the W-specific repeated sequence elements were amplified site-specifically after deletion of a large part of the W chromosome occurred.

Evidence for different origin of sex chromosomes in snakes, birds, and mammals and step-wise differentiation of snake sex chromosomes

All snake species exhibit genetic sex determination with the ZZ/ZW type of sex chromosomes. To investigate the origin and evolution of snake sex chromosomes, we constructed, by FISH, a cytogenetic map of the Japanese four-striped rat snake (Elaphe quadrivirgata) with 109 cDNA clones. Eleven of the 109 clones were localized to the Z chromosome. All human and chicken homologues of the snake Z-linked genes were located on autosomes, suggesting that the sex chromosomes of snakes, mammals, and birds were all derived from different autosomal pairs of the common ancestor. We mapped the 11 Z-linked genes of E. quadrivirgata to chromosomes of two other species, the Burmese python (Python molurus bivittatus) and the habu (Trimeresurus flavoviridis), to investigate the process of W chromosome differentiation. All and 3 of the 11 clones were localized to both the Z and W chromosomes in P. molurus and E. quadrivirgata, respectively, whereas no cDNA clones were mapped to the W chromosome in T. flavoviridis. Comparative mapping revealed that the sex chromosomes are only slightly differentiated in P. molurus, whereas they are fully differentiated in T. flavoviridis, and E. quadrivirgata is at a transitional stage of sex-chromosome differentiation. The differentiation of sex chromosomes was probably initiated from the distal region on the short arm of the protosex chromosome of the common ancestor, and then deletion and heterochromatization progressed on the sex-specific chromosome from the phylogenetically primitive boids to the more advanced viperids.

A novel family of repetitive DNA sequences amplified site-specifically on the W chromosomes in Neognathous birds

A novel family of repetitive DNA sequences was molecularly cloned from ApaI-digested genomic DNA of two Galliformes species, Japanese quail (Coturnix japonica) and guinea fowl (Numida meleagris), and characterized by chromosome in-situ hybridization and filter hybridization. Both the repeated sequence elements produced intensely painted signals on the W chromosomes, whereas they weakly hybridized to whole chromosomal regions as interspersed-type repetitive sequences. The repeated elements of the two species had high similarity of nucleotide sequences, and cross-hybridized to chromosomes of two other Galliformes species, chicken (Gallus gallus) and blue-breasted quail (Coturnix chinensis). The nucleotide sequences were conserved in three other orders of Neognathous birds, the Strigiformes, Gruiformes and Falconiformes, but not in Palaeognathous birds, the Struthioniformes and Tinamiformes, indicating that the repeated sequence elements were amplified on the W chromosomes in the lineage of Neognathous birds after the common ancestor diverged into the Palaeognathae and Neognathae. They are components of the W heterochromatin in Neognathous birds, and a good molecular cytogenetic marker for estimating the phylogenetic relationships and for clarifying the origin of the sex chromosome heterochromatin and the process of sex chromosome differentiation in birds.

Inefficient purifying selection: the mammalian Y chromosome in the rodent genus Mus

Two related genes with potentially similar functions, one on the Y chromosome and one on the X chromosome, were examined to determine if they evolved differently because of their chromosomal positions. Six hundred fifty-seven base pairs of coding sequence of Jarid1d (Smcy) on the Y chromosome and Jarid1c (Smcx) on the X chromosome were sequenced in 13 rodent taxa. An analysis of replacement and silent substitutions, using a counting method designed for samples with small evolutionary distances, showed a significant difference between the two genes. The different patterns of replacement and silent substitutions within Jarid1d and Jarid1c may be a result of evolutionary mechanisms that are particularly strong on the Y chromosome because of its unique properties. These findings are similar to results of previous studies of Y chromosomal genes in these and other mammalian taxa, suggesting that genes on the mammalian Y evolve in a chromosome-specific manner.

The evolution of the spindlin gene in birds: sequence analysis of an intron of the spindlin W and Z gene reveals four major divisions of the Psittaciformes

The Psittaciformes (parrots, parakeets) are among the most widely held captive birds. Yet, their evolution and their phylogenetic relationships have been relatively little studied. This paper describes the phylogenetic relationships between a number of Psittaciformes as derived from the sequences of the third intron of the Z-chromosomal and W-chromosomal spindlin genes. The Z-chromosomal sequences of the kakapo (Strigops habroptilus), the kea (Nestor notabilis), and the kaka (Nestor meridionalis) from New Zealand form a cluster which is the sister group to all other Psittaciformes. The results show further that the Z-chromosomal sequences of the other species can be divided into two groups based on the occurrence of a sequence element ACCCT. The group with the insert (A) is mainly from species with an Australasian geographical distribution and includes such species as the Lories (Lorius, etc.), the budgerigar (Melospittacus undulatus), and the rosellas (Platycercus). It also includes the African lovebirds (Agapornidae), which are the only representative of group A outside Australasia. Group B, without the insert, includes the neotropical parrots and parakeets such as the amazons (Amazona, etc.), the macaws (Ara, etc.), and the conures (Aratinga, etc.), the Australian Cacatuini and the African species such as the African grey parrot (Psittacus erithacus) as well as Coracopsis vasa from Madagascar and Psittrichas fulgidus from New Guinea. The W-chromosomal sequence data show that another division of the Psittacidae is found in the replacement of a pyrimidine-rich segment occurring in many non-psittacines as well as the kakapo (S. habroptilus), the kea (N. notabilis), the kaka (N. meridionalis), and the Cacatuini by a microsatellite consisting of a variable number of TATTA monomers in the other Psittaciformes. The results support a Gondwanan origin of the Psittaciformes and the suggestion that paleogeographic events were a major force in psittacine divergence.

Phylogenetics of modern birds in the era of genomics

In the 14 years since the first higher-level bird phylogenies based on DNA sequence data, avian phylogenetics has witnessed the advent and maturation of the genomics era, the completion of the chicken genome and a suite of technologies that promise to add considerably to the agenda of avian phylogenetics. In this review, we summarize current approaches and data characteristics of recent higher-level bird studies and suggest a number of as yet untested molecular and analytical approaches for the unfolding tree of life for birds. A variety of comparative genomics strategies, including adoption of objective quality scores for sequence data, analysis of contiguous DNA sequences provided by large-insert genomic libraries, and the systematic use of retroposon insertions and other rare genomic changes all promise an integrated phylogenetics that is solidly grounded in genome evolution. The avian genome is an excellent testing ground for such approaches because of the more balanced representation of single-copy and repetitive DNA regions than in mammals. Although comparative genomics has a number of obvious uses in avian phylogenetics, its application to large numbers of taxa poses a number of methodological and infrastructural challenges, and can be greatly facilitated by a 'community genomics' approach in which the modest sequencing throughputs of single PI laboratories are pooled to produce larger, complementary datasets. Although the polymerase chain reaction era of avian phylogenetics is far from complete, the comparative genomics era-with its ability to vastly increase the number and type of molecular characters and to provide a genomic context for these characters-will usher in a host of new perspectives and opportunities for integrating genome evolution and avian phylogenetics.

Male-driven evolution in closely related species of the mouse genus Mus

Recently, other researchers have found that closely related primate species had a lower male-to-female mutation rate ratio (alpha) than distantly related species. To determine if this is a general phenomenon affecting other mammalian orders, eleven species or subspecies of the rodent genus Mus and two outgroup species were compared. Intron sequences from a gene in the nonrecombining region of the Y chromosome Jarid1d (Smcy) and its X chromosomal gametolog, Jarid1c (Smcx), were analyzed in a phylogenetic context. The male-to-female mutation rate ratio for all thirteen taxa is approximately 2.5, which is similar to previous estimates in more distantly related rodents. However, when branches with lengths of more than 2.5% were removed from the analysis, the male-to-female mutation rate ratio dropped to 0.9. Thus, in closely related rodents, as in closely related primates, the male-to-female mutation rate ratio is lower than expected.

Evolutionary strata on the chicken Z chromosome: implications for sex chromosome evolution

The human X chromosome exhibits four "evolutionary strata," interpreted to represent distinct steps in the process whereby recombination became arrested between the proto X and proto Y. To test if this is a general feature of sex chromosome evolution, we studied the Z-W sex chromosomes of birds, which have female rather than male heterogamety and evolved from a different autosome pair than the mammalian X and Y. Here we analyze all five known gametologous Z-W gene pairs to investigate the "strata" hypothesis in birds. Comparisons of the rates of synonymous substitution and intronic divergence between Z and W gametologs reveal the presence of at least two evolutionary strata spread over the p and q arms of the chicken Z chromosome. A phylogenetic analysis of intronic sequence data from different avian lineages indicates that Z-W recombination ceased in the oldest stratum (on Zq; CHD1Z, HINTZ, and SPINZ) 102-170 million years ago (MYA), before the split of the Neoaves and Eoaves. However, recombination continued in the second stratum (on Zp; UBAP2Z and ATP5A1Z) until after the divergence of extant avian orders, with Z and W diverging 58-85 MYA. Our data suggest that progressive and stepwise cessation of recombination is a general feature behind sex chromosome evolution.

Evolutionary strata on the mouse X chromosome correspond to strata on the human X chromosome

Lahn and Page previously observed that genes on the human X chromosome were physically arranged along the chromosome in "strata," roughly ordered by degree of divergence from related genes on the Y chromosome. They hypothesized that this ordering results from a historical series of suppressions of recombination along the mammalian Y chromosome, thereby allowing formerly recombining X and Y chromosomal genes to diverge independently. Here predictions of this hypothesis are confirmed in a nonprimate mammalian order, Rodentia, through an analysis of eight gene pairs from the X and Y chromosomes of the house mouse, Mus musculus. The mouse X chromosome has been rearranged relative to the human X, so strata were not found in the same physical order on the mouse X. However, based on synonymous evolutionary distances, X-linked genes in M. musculus fall into the same strata as orthologous genes in humans, as predicted. The boundary between strata 2 and 3 is statistically significant, but the boundary between strata 1 and 2 is not significant in mice. An analysis of smaller fragments of Smcy, Smcx, Zfy, and Zfx from seven species of Mus confirmed that the strata in Mus musculus were representative of the genus Mus.

Sex Differences in Structure and Expression of the Sex Chromosome Genes CHD1Z and CHD1W in Zebra Finches

Genes on the sex chromosomes are unique because of their sex-specific inheritance. One question is whether homologous gene pairs on the sex chromosomes, which have diverged in their sequence, have acquired different functions. We have analyzed the first homologous pair of genes (CHD1Z and CHD1W) discovered on the avian Z and W sex chromosomes of the zebra finch (Taeniopygia guttata) to examine whether functional differences may have evolved. Sequence analysis revealed that the two genes maintained a high degree of similarity especially within the C, H, and D domains, but outside of these regions larger differences were observed. Expression studies showed that CHD1W was unique to females and has the potential to produce a protein that CHD1Z does not. CHD1Z mRNA was expressed at a higher level in the male brain than in the female brain at various post-hatch ages. Reporter constructs containing the 5' flanking regions of each gene showed they had the ability to drive reporter expression in primary cell cultures. The 5' flanking region sequence of CHD1Z and CHD1W exhibited little homology, and differences in putative promoter elements were apparent. These differences between CHD1Z and CHD1W suggest that the two proteins may have diverged in their function.

Life history and the male mutation bias

If DNA replication is a major cause of mutation, then those life-history characters, which are expected to affect the number of male germline cell divisions, should also affect the male to female mutation bias (alpha(m)). We tested this hypothesis by comparing several clades of bird species, which show variation both in suitable life-history characters (generation time as measured by age at first breeding and sexual selection as measured by frequency of extrapair paternity) and in alpha(m), which was estimated by comparing Z-linked and W-linked substitution rates in gametologous introns. Alpha(m) differences between clades were found to positively covary with both generation time and sexual selection, as expected if DNA replication causes mutation. The effects of extrapair paternity frequency on alpha(m) suggests that increased levels of sexual selection cause higher mutation rates, which offers an interesting solution to the paradox of the loss of genetic variance associated with strong directional sexual selection. We also used relative rate tests to examine whether the observed differences in alpha(m) between clades were due to differences in W-linked or Z-linked substitution rates. In one case, a significant difference in alpha(m) between two clades was shown to be due to W-linked rates and not Z-linked rates, a result that suggests that mutation rates are not determined by replication alone.

Inter-familial relationships of the shorebirds (Aves: Charadriiformes) based on nuclear DNA sequence data

BACKGROUND

Phylogenetic hypotheses of higher-level relationships in the order Charadriiformes based on morphological data, partly disagree with those based on DNA-DNA hybridisation data. So far, these relationships have not been tested by analysis of DNA sequence data. Herein we utilize 1692 bp of aligned, nuclear DNA sequences obtained from 23 charadriiform species, representing 15 families. We also test earlier suggestions that bustards and sandgrouses may be nested with the charadriiforms. The data is analysed with methods based on the parsimony and maximum-likelihood criteria.

RESULTS

Several novel phylogenetic relationships were recovered and strongly supported by the data, regardless of which method of analysis was employed. These include placing the gulls and allied groups as a sistergroup to the sandpiper-like birds, and not to the plover-like birds. The auks clearly belong to the clade with the gulls and allies, and are not basal to most other charadriiform birds as suggested in analyses of morphological data. Pluvialis, which has been supposed to belong to the plover family (Charadriidae), represents a basal branch that constitutes the sister taxon to a clade with plovers, oystercatchers and avocets. The thick-knees and sheathbills unexpectedly cluster together.

CONCLUSION

The DNA sequence data contains a strong phylogenetic signal that results in a well-resolved phylogenetic tree with many strongly supported internodes. Taxonomically it is the most inclusive study of shorebird families that relies on nucleotide sequences. The presented phylogenetic hypothesis provides a solid framework for analyses of macroevolution of ecological, morphological and behavioural adaptations observed within the order Charadriiformes.

Two new avian mitochondrial genomes (penguin and goose) and a summary of bird and reptile mitogenomic features

We report complete mitochondrial (mt) genomes for a penguin (little blue, Eudyptula minor) and a goose (greater white-fronted, Anser albifrons). A revised annotation of avian and reptile mt genomes has been carried out, which improves consistency of labeling gene start and stop positions. In conjunction with this, a summary of mt gene features is presented and a number of conserved patterns and interesting differences identified. The protein-coding genes from the two new genomes were analysed together with those from 17 other birds plus outgroup (reptile) taxa. The unrooted amino acid tree from 19 avian genomes was locally stable with many high bootstrap values using several maximum likelihood methods. In particular, Anseriformes (goose and duck) grouped strongly with Galliformes (chicken) to form Gallianseres, while the penguin paired firmly with the stork. The position where the outgroup joined the avian tree varied with the combination of outgroup taxa used. The three best supported positions of the root were passerine, but the traditional rooting position between paleognaths and neognaths could not be excluded.

Sex and death: CHD1Z associated with high mortality in moorhens

Sex ratios in clutches of moorhens (Gallinula chloropus) in Britain were measured on 83 chicks using the sex-linked CHD1 gene (Chromo-helicase/ATPase-DNA binding protein 1). Among birds, the female is the heterogametic sex (Z and W chromosomes), and the male is homogametic (two copies of the Z chromosome). We report variation among the PCR-amplified fragments of the CHD1Z, and the death of nearly all heterozygous male chicks (92%). In contrast, survivorship among females and homozygote males was 54-60%. Mortality in male heterozygotes was significantly higher than that of male homozygotes (P < 0.001). Chick and egg biometrics were not significantly different between these males. The CHD1Z was unlikely to be directly responsible but may have been hitchhiked by the causal gene(s). The observations appear to follow a classic underdominance (heterozygote inferiority) pattern, but raise the paradoxical question of why one form of the Z chromosome has not been fixed, as is expected from evolutionary theory. We discuss possible explanations and include a survey of British populations based on skin specimens.

The effects of nucleotide substitution model assumptions on estimates of nonparametric bootstrap support

The use of parameter-rich substitution models in molecular phylogenetics has been criticized on the basis that these models can cause a reduction both in accuracy and in the ability to discriminate among competing topologies. We have explored the relationship between nucleotide substitution model complexity and nonparametric bootstrap support under maximum likelihood (ML) for six data sets for which the true relationships are known with a high degree of certainty. We also performed equally weighted maximum parsimony analyses in order to assess the effects of ignoring branch length information during tree selection. We observed that maximum parsimony gave the lowest mean estimate of bootstrap support for the correct set of nodes relative to the ML models for every data set except one. For several data sets, we established that the exact distribution used to model among-site rate variation was critical for a successful phylogenetic analysis. Site-specific rate models were shown to perform very poorly relative to gamma and invariable sites models for several of the data sets most likely because of the gross underestimation of branch lengths. The invariable sites model also performed poorly for several data sets where this model had a poor fit to the data, suggesting that addition of the gamma distribution can be critical. Estimates of bootstrap support for the correct nodes often increased under gamma and invariable sites models relative to equal rates models. Our observations are contrary to the prediction that such models cause reduced confidence in phylogenetic hypotheses. Our results raise several issues regarding the process of model selection, and we briefly discuss model selection uncertainty and the role of sensitivity analyses in molecular phylogenetics.

Dosage compensation: do birds do it as well?

In birds males carry ZZ and females ZW sex chromosomes, and it has been proposed that there is no dosage compensation in the expression of sex-linked genes. However, recent data suggest the opposite, indicating that male and female birds might demonstrate similar levels of expression of Z-linked genes. If they do, the equalization between the sexes is probably not achieved by inactivation of one of the male Z chromosomes. Other possible mechanisms include the transcription of Z-linked genes being upregulated in females or downregulated in males, or equalization at the translation stage in either sex. A recently identified hypermethylated region on the Z chromosome, with similarities to the X inactivation centre on the mammalian X chromosome, might play a part in this process or have a role in avian sex determination.

Contrasting levels of nucleotide diversity on the avian Z and W sex chromosomes

Sex chromosomes may provide a context for studying the local effects of mutation rate on molecular evolution, since the two types of sex chromosomes are generally exposed to different mutational environments in male and female germ lines. Importantly, recent studies of some vertebrates have provided evidence for a higher mutation rate among males than among females. Thus, in birds, the Z chromosome, which spends two thirds of its time in the male germ line, is exposed to more mutations than the female-specific W chromosome. We show here that levels of nucleotide diversity are drastically higher on the avian Z chromosome than in paralogous sequences on the W chromosome. In fact, no intraspecific polymorphism whatsoever was seen in about 3.4 kb of CHD1W intron sequence from a total of >150 W chromosome copies of seven different bird species. In contrast, the amount of genetic variability in paralogous sequences on the Z chromosome was significant, with an average pairwise nucleotide diversity (d) of 0.0020 between CHD1Z introns and with 37 segregating sites in a total of 3.8 kb of Z sequence. The contrasting levels of genetic variability on the avian sex chromosomes are thus in a direction predicted from a male-biased mutation rate. However, although a low gene number, as well as some other factors, argues against background selection and/or selective sweeps shaping the genetic variability of the avian W chromosome, we cannot completely exclude selection as a contributor to the low levels of variation on the W chromosome.

Evolutionary genetics. Clonal inheritance of avian mitochondrial DNA

We have taken a new approach to test the commonly accepted, but recently questioned, principle of clonal inheritance of vertebrate mitochondrial DNA (mtDNA) by relating its inheritance to a female-specific marker of nuclear DNA. Whereas this is impossible in organisms with male heterogamy (such as mammals), we show here that genealogies of mtDNA and the female-specific W chromosome of a bird species are completely concordant. Our results indicate that inheritance of mtDNA is free of detectable recombination effects over an evolutionary timescale.

Multiple and independent cessation of recombination between avian sex chromosomes

Birds are characterized by female heterogamety; females carry the Z and W sex chromosomes, while males have two copies of the Z chromosome. We suggest here that full differentiation of the Z and W sex chromosomes of birds did not take place until after the split of major contemporary lineages, in the late Cretaceous. The ATP synthase alpha-subunit gene is now present in one copy each on the nonrecombining part of the W chromosome (ATP5A1W) and on the Z chromosome (ATP5A1Z). This gene seems to have evolved on several independent occasions, in different lineages, from a state of free recombination into two sex-specific and nonrecombining variants. ATP5A1W and ATP5A1Z are thus more similar within orders, relative to what W (or Z) are between orders. Moreover, this cessation of recombination apparently took place at different times in different lineages (estimated at 13, 40, and 65 million years ago in Ciconiiformes, Galliformes, and Anseriformes, respectively). We argue that these observations are the result of recent and traceable steps in the process where sex chromosomes gradually cease to recombine and become differentiated. Our data demonstrate that this process, once initiated, may occur independently in parallel in sister lineages.