Recombination and nucleotide diversity in the sex chromosomal pseudoautosomal region of the emu, Dromaius novaehollandiae

The evolutionary maintenance of ancient recombining sex chromosomes in the ostrich

Sex chromosomes have evolved repeatedly across the tree of life and often exhibit extreme size dimorphism due to genetic degeneration of the sex-limited chromosome (e.g. the W chromosome of some birds and Y chromosome of mammals). However, in some lineages, ancient sex-limited chromosomes have escaped degeneration. Here, we study the evolutionary maintenance of sex chromosomes in the ostrich (Struthio camelus), where the W remains 65% the size of the Z chromosome, despite being more than 100 million years old. Using genome-wide resequencing data, we show that the population scaled recombination rate of the pseudoautosomal region (PAR) is higher than similar sized autosomes and is correlated with pedigree-based recombination rate in the heterogametic females, but not homogametic males. Genetic variation within the sex-linked region (SLR) (π = 0.001) was significantly lower than in the PAR, consistent with recombination cessation. Conversely, genetic variation across the PAR (π = 0.0016) was similar to that of autosomes and dependent on local recombination rates, GC content and to a lesser extent, gene density. In particular, the region close to the SLR was as genetically diverse as autosomes, likely due to high recombination rates around the PAR boundary restricting genetic linkage with the SLR to only ~50Kb. The potential for alleles with antagonistic fitness effects in males and females to drive chromosome degeneration is therefore limited. While some regions of the PAR had divergent male-female allele frequencies, suggestive of sexually antagonistic alleles, coalescent simulations showed this was broadly consistent with neutral genetic processes. Our results indicate that the degeneration of the large and ancient sex chromosomes of the ostrich may have been slowed by high recombination in the female PAR, reducing the scope for the accumulation of sexually antagonistic variation to generate selection for recombination cessation.

Whole-genome phylogeography of the blue-faced honeyeater (Entomyzon cyanotis) and discovery and characterization of a neo-Z chromosome

Whole-genome surveys of genetic diversity and geographic variation often yield unexpected discoveries of novel structural variation, which long-read DNA sequencing can help clarify. Here, we report on whole-genome phylogeography of a bird exhibiting classic vicariant geographies across Australia and New Guinea, the blue-faced honeyeater (Entomyzon cyanotis), and the discovery and characterization of a novel neo-Z chromosome by long-read sequencing. Using short-read genome-wide SNPs, we inferred population divergence events within E. cyanotis across the Carpentarian and other biogeographic barriers during the Pleistocene (~0.3-1.7 Ma). Evidence for introgression between nonsister populations supports a hypothesis of reticulate evolution around a triad of dynamic barriers around Pleistocene Lake Carpentaria between Australia and New Guinea. During this phylogeographic survey, we discovered a large (134 Mbp) neo-Z chromosome and we explored its diversity, divergence and introgression landscape. We show that, as in some sylvioid passerine birds, a fusion occurred between chromosome 5 and the Z chromosome to form a neo-Z chromosome; and in E. cyanotis, the ancestral pseudoautosomal region (PAR) appears nonrecombinant between Z and W, along with most of the fused chromosome 5. The added recombination-suppressed portion of the neo-Z (~37.2 Mbp) displays reduced diversity and faster population genetic differentiation compared with the ancestral-Z. Yet, the new PAR (~17.4 Mbp) shows elevated diversity and reduced differentiation compared to autosomes, potentially resulting from introgression. In our case, long-read sequencing helped clarify the genomic landscape of population divergence on autosomes and sex chromosomes in a species where prior knowledge of genome structure was still incomplete.

Evolutionary Dynamics of Sex Chromosomes of Paleognathous Birds

Standard models of sex chromosome evolution propose that recombination suppression leads to the degeneration of the heterogametic chromosome, as is seen for the Y chromosome in mammals and the W chromosome in most birds. Unlike other birds, paleognaths (ratites and tinamous) possess large nondegenerate regions on their sex chromosomes (PARs or pseudoautosomal regions). It remains unclear why these large PARs are retained over >100 Myr, and how this retention impacts the evolution of sex chromosomes within this system. To address this puzzle, we analyzed Z chromosome evolution and gene expression across 12 paleognaths, several of whose genomes have recently been sequenced. We confirm at the genomic level that most paleognaths retain large PARs. As in other birds, we find that all paleognaths have incomplete dosage compensation on the regions of the Z chromosome homologous to degenerated portions of the W (differentiated regions), but we find no evidence for enrichments of male-biased genes in PARs. We find limited evidence for increased evolutionary rates (faster-Z) either across the chromosome or in differentiated regions for most paleognaths with large PARs, but do recover signals of faster-Z evolution in tinamou species with mostly degenerated W chromosomes, similar to the pattern seen in neognaths. Unexpectedly, in some species, PAR-linked genes evolve faster on average than genes on autosomes, suggested by diverse genomic features to be due to reduced efficacy of selection in paleognath PARs. Our analysis shows that paleognath Z chromosomes are atypical at the genomic level, but the evolutionary forces maintaining largely homomorphic sex chromosomes in these species remain elusive.

Broad-scale recombination pattern in the primitive bird Rhea americana (Ratites, Palaeognathae)

Birds have genomic and chromosomal features that make them an attractive group to analyze the evolution of recombination rate and the distribution of crossing over. Yet, analyses are biased towards certain species, especially domestic poultry and passerines. Here we analyze for the first time the recombination rate and crossover distribution in the primitive ratite bird, Rhea americana (Rheiformes, Palaeognathae). Using a cytogenetic approach for in situ mapping of crossovers we found that the total genetic map is 3050 cM with a global recombination rate of 2.1 cM/Mb for female rheas. In the five largest macrobivalents there were 3 or more crossovers in most bivalents. Recombination rates for macrobivalents ranges between 1.8-2.1 cM/Mb and the physical length of their synaptonemal complexes is highly predictive of their genetic lengths. The crossover rate at the pseudoautosomal region is 2.1 cM/Mb, similar to those of autosomal pairs 5 and 6 and only slightly higher compared to other macroautosomes. It is suggested that the presence of multiple crossovers on the largest macrobivalents is a feature common to many avian groups, irrespective of their position throughout phylogeny. These data provide new insights to analyze the heterogeneous recombination landscape of birds.

The Pseudoautosomal Regions of the U/V Sex Chromosomes of the Brown Alga Ectocarpus Exhibit Unusual Features

The recombining regions of sex chromosomes (pseudoautosomal regions, PARs) are predicted to exhibit unusual features due to their being genetically linked to the nonrecombining, sex-determining region. This phenomenon is expected to occur in both diploid (XY, ZW) and haploid (UV) sexual systems, with slightly different consequences for UV sexual systems because of the absence of masking during the haploid phase (when sex is expressed) and because there is no homozygous sex in these systems. Despite a considerable amount of theoretical work on PAR genetics and evolution, these genomic regions have remained poorly characterized empirically. We show here that although the PARs of the U/V sex chromosomes of the brown alga Ectocarpus recombine at a similar rate to autosomal regions of the genome, they exhibit many genomic features typical of nonrecombining regions. The PARs were enriched in clusters of genes that are preferentially, and often exclusively, expressed during the sporophyte generation of the life cycle, and many of these genes appear to have evolved since the Ectocarpales diverged from other brown algal lineages. A modeling-based approach was used to investigate possible evolutionary mechanisms underlying this enrichment in sporophyte-biased genes. Our results are consistent with the evolution of the PAR in haploid systems being influenced by differential selection pressures in males and females acting on alleles that are advantageous during the sporophyte generation of the life cycle.

Genomic structures of the kW1 loci on the Z and W chromosomes in ratite birds: structural changes at an early stage of W chromosome differentiation

The W chromosome of ratite birds shows minimal morphological differentiation and retains homology of genetic linkage and gene order with a substantial stretch of the Z chromosome; however, the molecular structure in the differentiated region is still not well known. The kW1 sequence was isolated from the kiwi as a W-specific DNA marker for PCR-based molecular sexing of ratite birds. In ratite W chromosomes, this sequence commonly contains a ∼200-bp deletion. To characterize the very early event of avian sex chromosome differentiation, we performed molecular cytogenetic analyses of kW1 and its flanking sequences in paleognathous and neognathous birds and reptiles. Female-specific repeats were found in the kW1-flanking sequence of the cassowary (Casuarius casuarius), and the repeats have been amplified in the pericentromeric region of the W chromosomes of ratites, which may have resulted from the cessation of meiotic recombination between the Z and W chromosomes at an early stage of sex chromosome differentiation. The presence of the kW1 sequence in neognathous birds and a crocodilian species suggests that the kW1 sequence was present in the ancestral genome of Archosauria; however, it disappeared in other reptilian taxa and several lineages of neognathous birds.

The effect of variable frequency of sexual reproduction on the genetic structure of natural populations of a cyclical parthenogen

Cyclical parthenogens are a valuable system in which to empirically test theoretical predictions as to the genetic consequences of sexual reproduction in natural populations, particularly if the frequency of sexual relative to asexual reproduction can be quantified. In this study, we used a series of lake populations of the cyclical parthenogen, Daphnia pulicaria, that vary consistently in their investment in sexual reproduction, to address the questions of whether the ecological variation in investment in sex is detectable at the genetic level, and if so, whether the genetic patterns seen are consistent with theoretical predictions. We show that there is variation in the genetic structure of these populations in a manner consistent with their investment in sexual reproduction. Populations engaging in a high frequency of sex were in Hardy-Weinberg and gametic phase equilibrium, and showed little genotypic differentiation across sampled years. In contrast, populations with a lower frequency of sex deviated widely from equilibrium, had reduced multilocus clonal diversity, and showed significant temporal genotypic deviation.

Diverse stages of sex-chromosome differentiation in tinamid birds: evidence from crossover analysis in Eudromia elegans and Crypturellus tataupa

All extant birds share the same sex-chromosome system: ZZ males and ZW females with striking differences in the stages of sex-chromosome differentiation between the primitive palaeognathus ratites and the large majority of avian species grouped within neognaths. Evolutionarily close to ratites is the neotropical order Tinamiformes that has been scarcely explored regarding their ZW pair morphology and constitution. Tinamous, when compared to ratites, constitute a large group among Palaeognathae, therefore, exploring the extent of homology between the Z and W chromosomes in this group might reveal key features on the evolution of the avian sex chromosomes. We mapped MLH1 foci that are crossover markers on pachytene bivalents to determine the size and localization of the homologous region shared by the Z and W chromosomes in two tinamous: Eudromia elegans and Crypturellus tataupa. We found that the homologous (pseudoautosomal) region differ significantly in size between these two species. They both have a single recombination event on the long arm of the acrocentric Z and W chromosomes. However, in E. elegans the pseudoautosomal region occupies one-fourth of the W chromosome, while in C. tataupa it is restricted to the tip of the long arm of the W. The W chromosomes in these two species differ in their heterochromatin content: in E. elegans it shows a terminal euchromatic segment and in C. tataupa is completely heterochromatic. These results show that tinamous have ZW pairs with more diversified stages of differentiation compared to ratites. Finally, the idea that the avian proto-sex chromosomes started to diverge from the end of the long arm towards the centromere of an acrocentric pair is discussed.

Sex chromosome evolution in amniotes: applications for bacterial artificial chromosome libraries

Variability among sex chromosome pairs in amniotes denotes a dynamic history. Since amniotes diverged from a common ancestor, their sex chromosome pairs and, more broadly, sex-determining mechanisms have changed reversibly and frequently. These changes have been studied and characterized through the use of many tools and experimental approaches but perhaps most effectively through applications for bacterial artificial chromosome (BAC) libraries. Individual BAC clones carry 100–200 kb of sequence from one individual of a target species that can be isolated by screening, mapped onto karyotypes, and sequenced. With these techniques, researchers have identified differences and similarities in sex chromosome content and organization across amniotes and have addressed hypotheses regarding the frequency and direction of past changes. Here, we review studies of sex chromosome evolution in amniotes and the ways in which the field of research has been affected by the advent of BAC libraries.

Genome evolution in Reptilia, the sister group of mammals

The genomes of birds and nonavian reptiles (Reptilia) are critical for understanding genome evolution in mammals and amniotes generally. Despite decades of study at the chromosomal and single-gene levels, and the evidence for great diversity in genome size, karyotype, and sex chromosome diversity, reptile genomes are virtually unknown in the comparative genomics era. The recent sequencing of the chicken and zebra finch genomes, in conjunction with genome scans and the online publication of the Anolis lizard genome, has begun to clarify the events leading from an ancestral amniote genome--predicted to be large and to possess a diverse repeat landscape on par with mammals and a birdlike sex chromosome system--to the small and highly streamlined genomes of birds. Reptilia exhibit a wide range of evolutionary rates of different subgenomes and, from isochores to mitochondrial DNA, provide a critical contrast to the genomic paradigms established in mammals.

Sex-specific population structure, natural selection, and linkage disequilibrium in a wild bird population as revealed by genome-wide microsatellite analyses

Background

Sexual dimorphism in ecologically important traits is widespread, yet the differences in the genomic architecture between the two sexes are largely unexplored. We employed a genome-wide multilocus approach to examine the sexual differences in population subdivision, natural selection and linkage disequilibrium (LD) in a wild Siberian jay (Perisoreus infaustus) population, using genotypes at a total of 107 autosomal and Z-chromosomal microsatellites.

Results

Mean observed heterozygosity was significantly higher in females (H_O = 0.567) than in males (H_O = 0.532), and autosomal markers (H_O = 0.561) were more variable than Z-chromosomal markers (H_O = 0.512). Genetic differentiation (F_ST = 0.002, P < 0.05) between the two sexes was low but significant and males were on average significantly more genetically related to each other than females. Genomescan analyses revealed that 3 out of 101 (3%) autosomal loci were under directional selection, while 4 out of 6 (67%) Z-chromosomal markers were indicated to be under balancing selection. This suggests a significantly greater but contrasting selection force on the Z-chromosome in comparison to autosomes, which is consistent with an overall significantly (P < 0.05) lower F_STvalue for Z-chromosomal (-0.014, 95% CI: -0.025 - -0.011) than for the autosomal loci (0.003, 95% CI: 0.001 - 0.004). Analysis of syntenic marker pairs revealed high levels of LD in both sexes but significantly (P < 0.05) lower levels of LD in the females both on autosomes and Z-chromosome, probably due to the higher rate of dispersal and the higher recombination rates on autosomes, as well as the pseudoautosomal markers. In both sexes LD decayed rapidly with genetic distance in a similar fashion on autosomes, while a more rapid decay of LD in Z-chromosome was detected in females than in males.

Conclusion

We conclude that there are many clear differences in genomic architecture between the sexes studied here which can be at least partly understood in the light of higher dispersal rate of females as compared to males and the unusual structure of the Z-chromosome of the species.

Variability in sex-determining mechanisms influences genome complexity in reptilia

In this review, we describe the history of amniote sex determination as a classic example of Darwinian evolution. We suggest that evolutionary changes in sex determination provide a foundation for understanding important aspects of chromosome and genome organization that otherwise appear haphazard in their origins and contents. Species with genotypic sex determination often possess heteromorphic sex chromosomes, whereas species with environmental sex determination lack them. Through a series of mutations followed by selection at key genes, sex-determining mechanisms have turned over many times throughout the amniote lineage. As a consequence, amniote genomes have undergone gains or losses of sex chromosomes. We review the genomic and ecological contexts in which either temperature-dependent or genotypic sex determination has evolved. Once genotypic sex determination emerges in a lineage, viviparity and heteromorphic sex chromosomes become more likely to evolve. For example, in extinct marine reptiles, genotypic sex determination apparently led to viviparity, which in turn facilitated their pelagic radiation. Sex chromosomes comprise genome regions that differ from autosomes in recombination rate, mutation rate, levels of polymorphism, and the presence of sex-determining and sexually antagonistic genes. In short, many aspects of amniote genome complexity, life history, and adaptive radiation appear contingent on evolutionary changes in sex-determining mechanisms.

Patterns of neutral genetic variation on recombining sex chromosomes

Many animals and plants have sex chromosomes that recombine over much of their length. Here we develop coalescent models for neutral sites on these chromosomes. The emphasis is on expected coalescence times (proportional to the expected amount of neutral genetic polymorphism), but we also derive some results for linkage disequilibria between neutral sites. We analyze the standard neutral model, a model with polymorphic Y chromosomes under balancing selection, and the invasion of a neo-Y chromosome. The results may be useful for testing hypotheses regarding how new sex chromosomes originate and how selection acts upon them.

Integrating genomes, brain and behavior in the study of songbirds

Songbirds share some essential traits but are extraordinarily diverse, allowing comparative analyses aimed at identifying specific genotype-phenotype associations. This diversity encompasses traits like vocal communication and complex social behaviors that are of great interest to humans, but that are not well represented in other accessible research organisms. Many songbirds are readily observable in nature and thus afford unique insight into the links between environment and organism. The distinctive organization of the songbird brain will facilitate analysis of genomic links to brain and behavior. Access to the zebra finch genome sequence will, therefore, prompt new questions and provide the ability to answer those questions.

Genome evolution in Reptilia: in silico chicken mapping of 12,000 BAC-end sequences from two reptiles and a basal bird

BACKGROUND

With the publication of the draft chicken genome and the recent production of several BAC clone libraries from non-avian reptiles and birds, it is now possible to undertake more detailed comparative genomic studies in Reptilia. Of interest in particular are the genomic events that transformed the large, repeat-rich genomes of mammals and non-avian reptiles into the minimalist chicken genome. We have used paired BAC end sequences (BESs) from the American alligator (Alligator mississippiensis), painted turtle (Chrysemys picta) and emu (Dromaius novaehollandiae) to investigate patterns of sequence divergence, gene and retroelement content, and microsynteny between these species and chicken.

RESULTS

From a total of 11,967 curated BESs, we successfully mapped 725, 773 and 2597 sequences in alligator, turtle, and emu, respectively, to sites in the draft chicken genome using a stringent BLAST protocol. Most commonly, sequences mapped to a single site in the chicken genome. Of 1675, 1828 and 2936 paired BESs obtained for alligator, turtle, and emu, respectively, a total of 34 (alligator, 2%), 24 (turtle, 1.3%) and 479 (emu, 16.3%) pairs were found to map with high confidence and in the correct orientation and with BAC-sized intermarker distances to single chicken chromosomes, including 25 such paired hits in emu mapping to the chicken Z chromosome. By determining the insert sizes of a subset of BAC clones from these three species, we also found a significant correlation between the intermarker distance in alligator and turtle and in chicken, with slopes as expected on the basis of the ratio of the genome sizes.

CONCLUSION

Our results suggest that a large number of small-scale chromosomal rearrangements and deletions in the lineage leading to chicken have drastically reduced the number of detected syntenies observed between the chicken and alligator, turtle, and emu genomes and imply that small deletions occurring widely throughout the genomes of reptilian and avian ancestors led to the ~50% reduction in genome size observed in birds compared to reptiles. We have also mapped and identified likely gene regions in hundreds of new BAC clones from these species.