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

A pair of gametologous genes provides further insights into avian comparative cytogenomics

Exploration of avian gametologous genes, i.e., homologous genes located on both the Z and W chromosomes, provides a crucial information about the underlying mechanism pertaining to the evolution of these chromosomes. The domestic chicken (Gallus gallus (Linnaeus 1758); GGA) traditionally serves as the primary reference subject of these comparative cytogenomic studies. Using bioinformatic, molecular (overgo BAC library scanning), and cytogenetic (BAC-based FISH) techniques, we have investigated in detail a pair of UBE2R2/UBE2R2L gametologs. By screening a gridded genomic jungle fowl BAC library, CHORI-261, with a short labeled UBE2R2L gene fragment called overgo probe, we detected seven specific clones. For three of them, CH261-019I23, CH261-105E16, and CH261-114G22, we identified their precise cytogenetic location on the Gallus gallus W chromosome (GGAW). They also co-localized with the UBAP2L2 gene on the, as was shown previously, along with the CH261-053P09 BAC clone also containing the GGAW-specific UBE2R2L DNA sequence. The fine mapping of the UBE2R2/UBE2R2L homologs in the chicken genome also shed the light on comparative cytogenetic aspects in birds. Our findings provided further evidence that bird genomes moderately changed only during evolution and are suitable for successful use of interspecies hybridization using both overgo-based BAC library screen and BAC-based FISH.

Dynamic Patterns of Sex Chromosome Evolution in Neognath Birds: Many Independent Barriers to Recombination at the ATP5F1A Locus

Avian sex chromosomes evolved after the divergence of birds and crocodilians from their common ancestor, so they are younger than the better-studied chromosomes of mammals. It has long been recognized that there may have been several stages to the evolution of avian sex chromosomes. For example, the CHD1 undergoes recombination in paleognaths but not neognaths. Genome assemblies have suggested that there may be variation in the timing of barriers to recombination among Neognathae, but there remains little understanding of the extent of this variability. Here, we look at partial sequences of ATP5F1A, which is on the avian Z and W chromosomes. It is known that recombination of this gene has independently ceased in Galliformes, Anseriformes, and at least five neoavian orders, but whether there are other independent cessations of recombination among Neoaves is not understood. We analyzed a combination of data extracted from published chromosomal-level genomes with data collected using PCR and cloning to identify Z and W copies in 22 orders. Our results suggest that there may be at least 19 independent cessations of recombination within Neognathae, and 3 clades that may still be undergoing recombination (or have only recently ceased recombination). Analyses of ATP5F1A protein sequences revealed an increased amino acid substitution rate for W chromosome gametologs, suggesting relaxed purifying selection on the W chromosome. Supporting this hypothesis, we found that the increased substitution rate was particularly pronounced for buried residues, which are expected to be more strongly constrained by purifying selection. This highlights the dynamic nature of avian sex chromosomes, and that this level of variation among clades means they should be a good system to understand sex chromosome evolution.

Avian Neo-Sex Chromosomes Reveal Dynamics of Recombination Suppression and W Degeneration

How the avian sex chromosomes first evolved from autosomes remains elusive as 100 million years (My) of divergence and degeneration obscure their evolutionary history. The Sylvioidea group of songbirds is interesting for understanding avian sex chromosome evolution because a chromosome fusion event ∼24 Ma formed "neo-sex chromosomes" consisting of an added (new) and an ancestral (old) part. Here, we report the complete female genome (ZW) of one Sylvioidea species, the great reed warbler (Acrocephalus arundinaceus). Our long-read assembly shows that the added region has been translocated to both Z and W, and whereas the added-Z has retained its gene order the added-W part has been heavily rearranged. Phylogenetic analyses show that recombination between the homologous added-Z and -W regions continued after the fusion event, and that recombination suppression across this region took several million years to be completed. Moreover, recombination suppression was initiated across multiple positions over the added-Z, which is not consistent with a simple linear progression starting from the fusion point. As expected following recombination suppression, the added-W show signs of degeneration including repeat accumulation and gene loss. Finally, we present evidence for nonrandom maintenance of slowly evolving and dosage-sensitive genes on both ancestral- and added-W, a process causing correlated evolution among orthologous genes across broad taxonomic groups, regardless of sex linkage.

Putative Epigenetic Biomarkers of Stress in Red Blood Cells of Chickens Reared Across Different Biomes

Production animals are constantly subjected to early adverse environmental conditions that influence the adult phenotype and produce epigenetic effects. CpG dinucleotide methylation in red blood cells (RBC) could be a useful epigenetic biomarker to identify animals subjected to chronic stress in the production environment. Here we compared a reduced fraction of the RBC methylome of chickens exposed to social isolation to non-exposed. These experiments were performed in two different locations: Brazil and Sweden. The aim was to identify stress-associated DNA methylation profiles in RBC across these populations, in spite of the variable conditions to which birds are exposed in each facility and their different lineages. Birds were increasingly exposed to a social isolation treatment, combined with food and water deprivation, at random periods of the day from weeks 1-4 after hatching. We then collected the RBC DNA from individuals and compared a reduced fraction of their methylome between the experimental groups using two bioinformatic approaches to identify differentially methylated regions (DMRs): one using fixed-size windows and another that preselected differential peaks with MACS2. Three levels of significance were used (P ≤ 0.05, P ≤ 0.005, and P ≤ 0.0005) to identify DMRs between experimental groups, which were then used for different analyses. With both of the approaches more DMRs reached the defined significance thresholds in BR individuals compared to SW. However, more DMRs had higher fold change values in SW compared to BR individuals. Interestingly, ChrZ was enriched above expectancy for the presence of DMRs. Additionally, when analyzing the locations of these DMRs in relation to the transcription starting site (TSS), we found three peaks with high DMR presence: 10 kb upstream, the TSS itself, and 20-40 kb downstream. Interestingly, these peaks had DMRs with a high presence (>50%) of specific transcription factor binding sites. Three overlapping DMRs were found between the BR and SW population using the most relaxed p-value (P ≤ 0.05). With the most stringent p-value (P ≤ 0.0005), we found 7 and 4 DMRs between treatments in the BR and SW populations, respectively. This study is the first approximation to identify epigenetic biomarkers of long-term exposure to stress in different lineages of production animals.

Transcriptome analysis of genes related to gonad differentiation and development in Muscovy ducks

BACKGROUND

Sex-related genes play a crucial role in gonadal differentiation into testes or ovaries. However, the genetic control of gonadal differentiation in Muscovy ducks remains unknown. Therefore, the objective of our study was to screen new candidate genes associated with ovarian and testicular development.

RESULTS

In this study, 24 males before gonadal differentiation (MB), 24 females before gonadal differentiation (FB), 24 males after gonadal differentiation (MA) and 24 females after gonadal differentiation (FA) were selected from Putian Muscovy ducks, forming 4 groups. RNA-Seq revealed 101.76 Gb of clean reads and 2800 differentially expressed genes (DEGs), including 46 in MB vs FB, 609 in MA vs FA, 1027 in FA vs FB, and 1118 in MA vs MB. A total of 146 signalling pathways were enriched by KEGG analysis, among which 20, 108, 108 and 116 signalling pathways were obtained in MB vs FB, MA vs MB, MA vs FA and FA vs FB, respectively. In further GO and KEGG analyses, a total of 21 candidate genes related to gonad differentiation and development in Muscovy ducks were screened. Among these, 9 genes were involved in the differentiation and development of the testes, and 12 genes were involved in the differentiation and development of the ovaries. In addition, RNA-Seq data revealed 2744 novel genes.

CONCLUSIONS

RNA-Seq data revealed 21 genes related to gonadal differentiation and development in Muscovy ducks. We further identified 12 genes, namely, WNT5B, HTRA3, RSPO3, BMP3, HNRNPK, NIPBL, CREB3L4, DKK3, UBE2R2, UBPL3KCMF1, ANXA2, and OSR1, involved in the differentiation and development of ovaries. Moreover, 9 genes, namely, TTN, ATP5A1, DMRT1, DMRT3, AMH, MAP3K1, PIK3R1, AGT and ADAMTSL1, were related to the differentiation and development of testes. Moreover, after gonadal differentiation, DMRT3, AMH, PIK3R1, ADAMTSL1, AGT and TTN were specifically highly expressed in males. WNT5B, ANXA2 and OSR1 were specifically highly expressed in females. These results provide valuable information for studies on the sex control of Muscovy ducks and reveal novel candidate genes for the differentiation and development of testes and ovaries.

Linkage mapping of a polymorphic plumage locus associated with intermorph incompatibility in the Gouldian finch (Erythrura gouldiae)

Colour polymorphism is known to facilitate speciation but the genetic basis of animal pigmentation and how colour polymorphisms contribute to speciation is poorly understood. Restricted recombination may promote linkage disequilibrium between the colour locus and incompatibility genes. Genomic rearrangement and the position of relevant loci within a chromosome are important factors that influence the frequency of recombination. Therefore, it is important to know the position of the colour locus, gene order and recombination landscape of the chromosome to understand the mechanism that generates incompatibilities between morphs. Recent studies showed remarkable pre- and postzygotic incompatibilities between sympatric colour morphs of the Gouldian finch (Erythrura gouldiae), in which head feather colour is genetically determined by a single sex-linked locus, Red. We constructed a genetic map for the Z chromosome of the Gouldian finch (male-specific map distance=131 cM), using 618 captive-bred birds and 34 microsatellite markers, to investigate the extent of inter- and intraspecific genomic rearrangements and variation in recombination rate within the Z chromosome. We refined the location of the Red locus to a ~7.2-cM interval in a region with a moderate recombination rate but outside the least-recombining, putative centromeric region. There was no evidence of chromosome-wide genomic rearrangements between the chromosomes carrying the red or black alleles with the current marker resolution. This work will contribute to identifying the causal gene, which will in turn enable alternative explanations for the association between incompatibility and colouration, such as fine-scale linkage disequilibrium, genomic rearrangements and pleiotropy, to be tested.

Sex differences and rapid estrogen signaling: A look at songbird audition

The actions of estrogens have been associated with brain differentiation and sexual dimorphism in a wide range of vertebrates. Here we consider the actions of brain-derived 'neuroestrogens' in the forebrain and the accompanying differences and similarities observed between males and females in a variety of species. We summarize recent evidence showing that baseline and fluctuating levels of neuroestrogens within the auditory forebrain of male and female zebra finches are largely similar, and that neuroestrogens enhance auditory representations in both sexes. With a comparative perspective we review evidence that non-genomic mechanisms of neuroestrogen actions are sexually differentiated, and we propose a working model for nonclassical estrogen signaling via the MAPK intracellular signaling cascade in the songbird auditory forebrain that is informed by the way sex differences may be compensated. This view may lead to a more comprehensive understanding of how sex influences estradiol-dependent modulation of sensorimotor representations.

Chaos of Rearrangements in the Mating-Type Chromosomes of the Anther-Smut Fungus Microbotryum lychnidis-dioicae

Sex chromosomes in plants and animals and fungal mating-type chromosomes often show exceptional genome features, with extensive suppression of homologous recombination and cytological differentiation between members of the diploid chromosome pair. Despite strong interest in the genetics of these chromosomes, their large regions of suppressed recombination often are enriched in transposable elements and therefore can be challenging to assemble. Here we show that the latest improvements of the PacBio sequencing yield assembly of the whole genome of the anther-smut fungus, Microbotryum lychnidis-dioicae (the pathogenic fungus causing anther-smut disease of Silene latifolia), into finished chromosomes or chromosome arms, even for the repeat-rich mating-type chromosomes and centromeres. Suppressed recombination of the mating-type chromosomes is revealed to span nearly 90% of their lengths, with extreme levels of rearrangements, transposable element accumulation, and differentiation between the two mating types. We observed no correlation between allelic divergence and physical position in the nonrecombining regions of the mating-type chromosomes. This may result from gene conversion or from rearrangements of ancient evolutionary strata, i.e., successive steps of suppressed recombination. Centromeres were found to be composed mainly of copia-like transposable elements and to possess specific minisatellite repeats identical between the different chromosomes. We also identified subtelomeric motifs. In addition, extensive signs of degeneration were detected in the nonrecombining regions in the form of transposable element accumulation and of hundreds of gene losses on each mating-type chromosome. Furthermore, our study highlights the potential of the latest breakthrough PacBio chemistry to resolve complex genome architectures.

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.

Advancing avian behavioral neuroendocrinology through genomics

Genome technologies are transforming all areas of biology, including the study of hormones, brain and behavior. Annotated reference genome assemblies are rapidly being produced for many avian species. Here we briefly review the basic concepts and tools used in genomics. We then consider how these are informing the study of avian behavioral neuroendocrinology, focusing in particular on lessons from the study of songbirds. We discuss the impact of having a complete "parts list" for an organism; the transformational potential of studying large sets of genes at once instead one gene at a time; the growing recognition that environmental and behavioral signals trigger massive shifts in gene expression in the brain; and the prospects for using comparative genomics to uncover the genetic roots of behavioral variation. Throughout, we identify promising new directions for bolstering the application of genomic information to further advance the study of avian brain and behavior.

Genome-brain-behavior interdependencies as a framework to understand hormone effects on learned behavior

Hormones have profound effects on the maturation and function of the zebra finch song system. Hormones often signal through receptors that directly or indirectly regulate transcription. In this way, hormones and the genome are functionally connected. Genome-brain-behavior interdependencies are often studied on evolutionary timescales but we can now apply and test these relationships on short timescales, relevant to an individual. Here, we begin to place patterns of hormone-related gene expression into the timeframe of an individual's lifespan to consider how hormones contribute to organization of neural systems necessary for learned behavior, and how they might signal during experience in ways that affect future behavior. This framework illustrates both how much investigations into genome and hormone function are intertwined, and how much we still need to learn.

TECHNICAL ADVANCES: A microarray for large-scale genomic and transcriptional analyses of the zebra finch (Taeniopygia guttata) and other passerines

The microarray technology has revolutionized biological research in the last decade. By monitoring the expression of many genes simultaneously, microarrays can elucidate gene function, as well as scan entire genomes for candidate genes encoding complex traits. However, because of high costs of sequencing and design, microarrays have largely been restricted to a few model species. Cross-species microarray (CSM) analyses, where microarrays are used for other species than the one they were designed for, have had varied success. We have conducted a CSM analysis by hybridizing genomic DNA from the common whitethroat (Sylvia communis) on a newly developed Affymetrix array designed for the zebra finch (Taeniopygia guttata), the Lund-zf array. The results indicate a very high potential for the zebra finch array to act as a CSM utility in other passerine birds. When hybridizing zebra finch genomic DNA, 98% of the gene representatives had higher signal intensities than the background cut-off, and for the common whitethroat, we found the equivalent proportion to be as high as 96%. This was surprising given the fact that finches and warblers diverged 25-50 million years ago, but may be explained by a relatively low sequence divergence between passerines (89-93%). Passerine birds are widely used in studies of ecology and evolution, and a zebra finch array that can be used for many species may have a large impact on future research directions.

RNA-seq transcriptome analysis of male and female zebra finch cell lines

The derivation of stably cultured cell lines has been critical to the advance of molecular biology. We profiled gene expression in the first two generally available cell lines derived from the zebra finch. Using Illumina RNA-seq, we generated ~93 million reads and mapped the majority to the recently assembled zebra finch genome. Expression of most Ensembl-annotated genes was detected, but over half of the mapped reads aligned outside annotated genes. The male-derived G266 line expressed Z-linked genes at a higher level than did the female-derived ZFTMA line, indicating persistence in culture of the distinctive lack of avian sex chromosome dosage compensation. Although these cell lines were not derived from neural tissue, many neurobiologically relevant genes were expressed, although typically at lower levels than in a reference sample from auditory forebrain. These cell lines recapitulate fundamental songbird biology and will be useful for future studies of songbird gene regulation and function.

A genetic polymorphism in the sex-linked ATP5A1 gene is associated with individual fitness in Ovenbirds (Seiurus aurocapilla)

While testing genetic sexing techniques in Ovenbirds (Seiurus aurocapilla), we found a genetic polymorphism in the ATP5A1 gene in 38% of individuals. The Z′ allele included changes in both intronic and exonic portions of the sequenced region, but there was no evidence that this changed the resulting ATP synthase product. Males that had one or more copies of this allele had higher relative body mass (mass corrected for size) than other genotypes. This allele was unrelated to stable isotope signatures, and so was not a useful predictor of latitude within the eastern portion of the Ovenbird breeding range. Future studies are needed to determine whether this polymorphism may be a useful geographic marker. This study is the first to link polymorphisms in the sex-linked ATP5A1 gene with fitness effects.

Patterns of molecular evolution of an avian neo-sex chromosome

Newer parts of sex chromosomes, neo-sex chromosomes, offer unique possibilities for studying gene degeneration and sequence evolution in response to loss of recombination and population size decrease. We have recently described a neo-sex chromosome system in Sylvioidea passerines that has resulted from a fusion between the first half (10 Mb) of chromosome 4a and the ancestral sex chromosomes. In this study, we report the results of molecular analyses of neo-Z and neo-W gametologs and intronic parts of neo-Z and autosomal genes on the second half of chromosome 4a in three species within different Sylvioidea lineages (Acrocephalidea, Timaliidae, and Alaudidae). In line with hypotheses of neo-sex chromosome evolution, we observe 1) lower genetic diversity of neo-Z genes compared with autosomal genes, 2) moderate synonymous and weak nonsynonymous sequence divergence between neo-Z and neo-W gametologs, and 3) lower GC content on neo-W than neo-Z gametologs. Phylogenetic reconstruction of eight neo-Z and neo-W gametologs suggests that recombination continued after the split of Alaudidae from the rest of the Sylvioidea lineages (i.e., after ~42.2 Ma) and with some exceptions also after the split of Acrocephalidea and Timaliidae (i.e., after ~39.4 Ma). The Sylvioidea neo-sex chromosome shares classical evolutionary features with the ancestral sex chromosomes but, as expected from its more recent origin, shows weaker divergence between gametologs.

The long and the short of avian W chromosomes: no evidence for gradual W shortening

The well-established view of the evolution of sex chromosome dimorphism is of a gradual genetic and morphological degeneration of the hemizygous chromosome. Yet, no large-scale comparative analysis exists to support this view. Here, we analysed karyotypes of 200 bird species to test whether the supposed directional changes occur in bird sex chromosomes. We found no support for the view that W chromosomes gradually become smaller over evolutionary time. On the contrary, the length of the W chromosome can fluctuate over short time scales, probably involving both shortening and elongation of non-coding regions. Recent discoveries of near-identical palindromes and neo-sex chromosomes in birds may also contribute to the observed variation. Further studies are now needed to investigate how chromosome morphology relates to its gene content, and whether the changes in size were driven by selection.

Microsatellite resources for Passeridae species: a predicted microsatellite map of the house sparrow Passer domesticus

We identified microsatellite sequences of potential utility in the house sparrow (Passer domesticus) and assigned their predicted genome locations. These sequences included newly isolated house sparrow loci, which we fully characterized. Many of the newly isolated loci were polymorphic in two other species of Passeridae: Berthelot's pipit Anthus berthelotii and zebra finch Taeniopygia guttata. In total, we identified 179 microsatellite markers that were either isolated directly from, or are of known utility in, the house sparrow. Sixty-seven of these markers were designed from unique sequences that we isolated from a house sparrow genomic library. These new markers were combined with 36 house sparrow markers isolated by other studies and 76 markers isolated from other passerine species but known to be polymorphic in the house sparrow. We utilized sequence homology to assign chromosomal locations for these loci in the assembled zebra finch genome. One hundred and thirty-four loci were assigned to 25 different autosomes and eight loci to the Z chromosome. Examination of the genotypes of known-sex house sparrows for 37 of the new loci revealed a W-linked locus and an additional Z-linked locus. Locus Pdoμ2, previously reported as autosomal, was found to be Z-linked. These loci enable the creation of powerful and cost-effective house sparrow multiplex primer sets for population and parentage studies. They can be used to create a house sparrow linkage map and will aid the identification of quantitative trait loci in passerine species.

Strong conservation of the bird Z chromosome in reptilian genomes is revealed by comparative painting despite 275 million years divergence

The divergence of lineages leading to extant squamate reptiles (lizards, snakes, and amphisbaenians) and birds occurred about 275 million years ago. Birds, unlike squamates, have karyotypes that are typified by the presence of a number of very small chromosomes. Hence, a number of chromosome rearrangements might be expected between bird and squamate genomes. We used chromosome-specific DNA from flow-sorted chicken (Gallus gallus) Z sex chromosomes as a probe in cross-species hybridization to metaphase spreads of 28 species from 17 families representing most main squamate lineages and single species of crocodiles and turtles. In all but one case, the Z chromosome was conserved intact despite very ancient divergence of sauropsid lineages. Furthermore, the probe painted an autosomal region in seven species from our sample with characterized sex chromosomes, and this provides evidence against an ancestral avian-like system of sex determination in Squamata. The avian Z chromosome synteny is, therefore, conserved albeit it is not a sex chromosome in these squamate species.

Comparative BAC-based mapping in the white-throated sparrow, a novel behavioral genomics model, using interspecies overgo hybridization

Background

The genomics era has produced an arsenal of resources from sequenced organisms allowing researchers to target species that do not have comparable mapping and sequence information. These new "non-model" organisms offer unique opportunities to examine environmental effects on genomic patterns and processes. Here we use comparative mapping as a first step in characterizing the genome organization of a novel animal model, the white-throated sparrow (Zonotrichia albicollis), which occurs as white or tan morphs that exhibit alternative behaviors and physiology. Morph is determined by the presence or absence of a complex chromosomal rearrangement. This species is an ideal model for behavioral genomics because the association between genotype and phenotype is absolute, making it possible to identify the genomic bases of phenotypic variation.

Findings

We initiated a genomic study in this species by characterizing the white-throated sparrow BAC library via filter hybridization with overgo probes designed for the chicken, turkey, and zebra finch. Cross-species hybridization resulted in 640 positive sparrow BACs assigned to 77 chicken loci across almost all macro-and microchromosomes, with a focus on the chromosomes associated with morph. Out of 216 overgos, 36% of the probes hybridized successfully, with an average number of 3.0 positive sparrow BACs per overgo.

Conclusions

These data will be utilized for determining chromosomal architecture and for fine-scale mapping of candidate genes associated with phenotypic differences. Our research confirms the utility of interspecies hybridization for developing comparative maps in other non-model organisms.

The sex-biased brain: sexual dimorphism in gene expression in two species of songbirds

Background

Despite virtually identical DNA sequences between the sexes, sexual dimorphism is a widespread phenomenon in nature. To a large extent the systematic differences between the sexes must therefore arise from processes involving gene regulation. In accordance, sexual dimorphism in gene expression is common and extensive. Genes with sexually dimorphic regulation are known to evolve rapidly, both in DNA sequence and in gene expression profile. Studies of gene expression in related species can shed light on the flexibility, or degree of conservation, of the gene expression profiles underlying sexual dimorphism.

Results

We have studied the extent of sexual dimorphism in gene expression in the brain of two species of songbirds, the zebra finch (Taeniopygia guttata) and the common whitethroat (Sylvia communis), using large-scale microarray technology. Sexual dimorphism in gene expression was extensive in both species, and predominantly sex-linked: most genes identified were male-biased and Z-linked. Interestingly, approximately 50% of the male-biased Z-linked genes were sex-biased only in one of the study species.

Conclusion

Our results corroborate the results of recent studies in chicken and zebra finch which have been interpreted as caused by a low degree of dosage compensation in female birds (i.e. the heterogametic sex). Moreover, they suggest that zebra finches and common whitethroats dosage compensate partly different sets of genes on the Z chromosome. It is possible that this pattern reflects differences in either the essentiality or the level of sexual antagonism of these genes in the respective species. Such differences might correspond to genes with different rates of evolution related to sexual dimorphism in the avian brain, and might therefore be correlated with differences between the species in sex-specific behaviours.

Karyotypic polymorphism of the zebra finch Z chromosome

We describe a karyotypic polymorphism on the zebra finch Z chromosome. This polymorphism was discovered because of a difference in the position of the centromere and because it occurs at varying frequencies in domesticated colonies in the USA and Germany and among two zebra finch subspecies. Using DNA fluorescent in situ hybridization to map specific Z genes and measurements of DNA replication, we show that this polymorphism is the result of a large pericentric inversion involving the majority of the chromosome. We sequenced a likely breakpoint for the inversion and found many repetitive sequences. Around the breakpoint, there are numerous repetitive sequences and several copies of PAK3 (p21-activated kinase 3)-related sequences (PAK3Z) which showed testes-specific expression by RT-PCR. Our findings further suggest that the sequenced genome of the zebra finch may be derived from a male heterozygote for the Z chromosome polymorphism. This finding, in combination with regional differences in the frequency of the polymorphism, has important consequences for future studies using zebra finches.

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.

No evidence for Z-chromosome rearrangements between the pied flycatcher and the collared flycatcher as judged by gene-based comparative genetic maps

Revealing the genetic basis of reproductive isolation is fundamental for understanding the speciation process. Chromosome speciation models propose a role for chromosomal rearrangements in promoting the build up of reproductive isolation between diverging populations and empirical data from several animal and plant taxa support these models. The pied flycatcher and the collared flycatcher are two closely related species that probably evolved reproductive isolation during geographical separation in Pleistocene glaciation refugia. Despite the short divergence time and current hybridization, these two species demonstrate a high degree of intrinsic post-zygotic isolation and previous studies have shown that traits involved in mate choice and hybrid viability map to the Z-chromosome. Could rearrangements of the Z-chromosome between the species explain their reproductive isolation? We developed high coverage Z-chromosome linkage maps for both species, using gene-based markers and large-scale SNP genotyping. Best order maps contained 57-62 gene markers with an estimated average density of one every 1-1.5 Mb. We estimated the recombination rates in flycatcher Z-chromosomes to 1.1-1.3 cM/Mb. A comparison of the maps of the two species revealed extensive co-linearity with no strong evidence for chromosomal rearrangements. This study does therefore not provide support the idea that sex chromosome rearrangements have caused the relatively strong post-zygotic reproductive isolation between these two Ficedula species.

Genomic and neural analysis of the estradiol-synthetic pathway in the zebra finch

Background

Steroids are small molecule hormones derived from cholesterol. Steroids affect many tissues, including the brain. In the zebra finch, estrogenic steroids are particularly interesting because they masculinize the neural circuit that controls singing and their synthesis in the brain is modulated by experience. Here, we analyzed the zebra finch genome assembly to assess the content, conservation, and organization of genes that code for components of the estrogen-synthetic pathway and steroid nuclear receptors. Based on these analyses, we also investigated neural expression of a cholesterol transport protein gene in the context of song neurobiology.

Results

We present sequence-based analysis of twenty steroid-related genes using the genome assembly and other resources. Generally, zebra finch genes showed high homology to genes in other species. The diversity of steroidogenic enzymes and receptors may be lower in songbirds than in mammals; we were unable to identify all known mammalian isoforms of the 3β-hydroxysteroid dehydrogenase and 17β-hydroxysteroid dehydrogenase families in the zebra finch genome assembly, and not all splice sites described in mammals were identified in the corresponding zebra finch genes. We did identify two factors, Nobox and NR1H2-RXR, that may be important for coordinated transcription of multiple steroid-related genes. We found very little qualitative overlap in predicted transcription factor binding sites in the genes for two cholesterol transport proteins, the 18 kDa cholesterol transport protein (TSPO) and steroidogenic acute regulatory protein (StAR). We therefore performed in situ hybridization for TSPO and found that its mRNA was not always detected in brain regions where StAR and steroidogenic enzymes were previously shown to be expressed. Also, transcription of TSPO, but not StAR, may be regulated by the experience of hearing song.

Conclusions

The genes required for estradiol synthesis and action are represented in the zebra finch genome assembly, though the complement of steroidogenic genes may be smaller in birds than in mammals. Coordinated transcription of multiple steroidogenic genes is possible, but results were inconsistent with the hypothesis that StAR and TSPO mRNAs are co-regulated. Integration of genomic and neuroanatomical analyses will continue to provide insights into the evolution and function of steroidogenesis in the songbird brain.

Sex bias and dosage compensation in the zebra finch versus chicken genomes: general and specialized patterns among birds

We compared global patterns of gene expression between two bird species, the chicken and zebra finch, with regard to sex bias of autosomal versus Z chromosome genes, dosage compensation, and evolution of sex bias. Both species appear to lack a Z chromosome-wide mechanism of dosage compensation, because both have a similar pattern of significantly higher expression of Z genes in males relative to females. Unlike the chicken Z chromosome, which has female-specific expression of the noncoding RNA MHM (male hypermethylated) and acetylation of histone 4 lysine 16 (H4K16) near MHM, the zebra finch Z chromosome appears to lack the MHM sequence and acetylation of H4K16. The zebra finch also does not show the reduced male-to-female (M:F) ratio of gene expression near MHM similar to that found in the chicken. Although the M:F ratios of Z chromosome gene expression are similar across tissues and ages within each species, they differ between the two species. Z genes showing the greatest species difference in M:F ratio were concentrated near the MHM region of the chicken Z chromosome. This study shows that the zebra finch differs from the chicken because it lacks a specialized region of greater dosage compensation along the Z chromosome, and shows other differences in sex bias. These patterns suggest that different avian taxa may have evolved specific compensatory mechanisms.

Neural expression and post-transcriptional dosage compensation of the steroid metabolic enzyme 17beta-HSD type 4

Background

Steroids affect many tissues, including the brain. In the zebra finch, the estrogenic steroid estradiol (E₂) is especially effective at promoting growth of the neural circuit specialized for song. In this species, only the males sing and they have a much larger and more interconnected song circuit than females. Thus, it was surprising that the gene for 17β-hydroxysteroid dehydrogenase type 4 (HSD17B4), an enzyme that converts E₂ to a less potent estrogen, had been mapped to the Z sex chromosome. As a consequence, it was likely that HSD17B4 was differentially expressed in males (ZZ) and females (ZW) because dosage compensation of Z chromosome genes is incomplete in birds. If a higher abundance of HSD17B4 mRNA in males than females was translated into functional enzyme in the brain, then contrary to expectation, males could produce less E₂ in their brains than females.

Results

Here, we used molecular and biochemical techniques to confirm the HSD17B4 Z chromosome location in the zebra finch and to determine that HSD17B4 mRNA and activity were detectable in the early developing and adult brain. As expected, HSD17B4 mRNA expression levels were higher in males compared to females. This provides further evidence of the incomplete Z chromosome inactivation mechanisms in birds. We detected HSD17B4 mRNA in regions that suggested a role for this enzyme in the early organization and adult function of song nuclei. We did not, however, detect significant sex differences in HSD17B4 activity levels in the adult brain.

Conclusions

Our results demonstrate that the HSD17B4 gene is expressed and active in the zebra finch brain as an E₂ metabolizing enzyme, but that dosage compensation of this Z-linked gene may occur via post-transcriptional mechanisms.

Female heterogamety and speciation: reduced introgression of the Z chromosome between two species of nightingales

Several lines of evidence suggest that the X chromosome plays a large role in intrinsic postzygotic isolation. The role of the Z chromosome in speciation is much less understood. To explore the role of the Z chromosome in reproductive isolation, we studied nucleotide variation in two closely related bird species, the Thrush Nightingale (Luscinia luscinia) and the Common Nightingale (L. megarhynchos). These species are isolated by incomplete prezygotic isolation and female hybrid sterility. We sequenced introns of four Z-linked and eight autosomal loci and analyzed patterns of polymorphism and divergence using a divergence-with-gene flow framework. Our results suggest that the nightingale species diverged approximately 1.8 Mya. We found strong evidence of gene flow after divergence in both directions, although more introgression occurred from L. megarhynchos into L. luscinia. Gene flow was significantly higher on the autosomes than on the Z chromosome. Our results support the idea that the Z chromosome plays an important role in intrinsic postzygotic isolation in birds, although it may also contribute to the evolution of prezygotic isolation through sexual selection. This highlights the similarities in the genetic basis of reproductive isolation between organisms with heterogametic males and organisms with heterogametic females during the early stages of speciation.

Molecular cytogenetics of the california condor: evolutionary and conservation implications

Evolutionary cytogenetic comparisons involved 5 species of birds (California condor, chicken, zebra finch, collared flycatcher and black stork) belonging to divergent taxonomic orders. Seventy-four clones from a condor BAC library containing 80 genes were mapped to condor chromosomes using FISH, and 15 clones containing 16 genes were mapped to the stork Z chromosome. Maps for chicken and finch were derived from genome sequence databases, and that for flycatcher from the published literature. Gene content and gene order were highly conserved when individual condor, chicken, and zebra finch autosomes were compared, confirming that these species largely retain karyotypes close to the ancestral condition for neognathous birds. However, several differences were noted: zebra finch chromosomes 1 and 1A are homologous to condor and chicken chromosomes 1, the CHUNK1 gene appears to have transposed on condor chromosome 1, condor chromosomes 4 and 9 and zebra finch chromosomes 4 and 4A are homologous to chicken chromosome arms 4q and 4p, and novel inversions on chromosomes 4, 12 and 13 were found. Condor and stork Z chromosome gene orders are collinear and differentiated by a series of inversions/transpositions when compared to chicken, zebra finch, or flycatcher; phylogenetic analyses suggest independent rearrangement along the chicken, finch, and flycatcher lineages.

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.

Genomic evidence for a large-Z effect

The 'large-X effect' suggests that sex chromosomes play a disproportionate role in adaptive evolution. Theoretical work indicates that this effect may be most pronounced in genetic systems with female heterogamety under both good-genes and Fisher's runaway models of sexual selection (males ZZ, females ZW). Here, I use a comparative genomic approach (alignments of several thousands of chicken-zebra finch-human-mouse-opossum orthologues) to show that avian Z-linked genes are highly overrepresented among those bird-mammalian orthologues that show evidence of accelerated rate of functional evolution in birds relative to mammals; the data suggest a twofold excess of such genes on the Z chromosome. A reciprocal analysis of genes accelerated in mammals found no evidence for an excess of X-linkage. This would be compatible with theoretical expectations for differential selection on sex-linked genes under male and female heterogamety, although the power in this case was not sufficient to statistically show that 'large-Z' was more pronounced than 'large-X'. Accelerated Z-linked genes include a variety of functional categories and are characterized by higher non-synonymous to synonymous substitution rate ratios than both accelerated autosomal and non-accelerated genes. This points at a genomic 'large-Z effect', which is widespread and of general significance for adaptive divergence in birds.

A linkage map of the zebra finch Taeniopygia guttata provides new insights into avian genome evolution

Passeriformes are the largest order of birds and one of the most widely studied groups in evolutionary biology and ecology. Until recently genomic tools in passerines relied on chicken genomic resources. Here we report the construction and analysis of a whole-genome linkage map for the zebra finch (Taeniopygia guttata) using a 354-bird pedigree. The map contains 876 SNPs dispersed across 45 linkage groups and we found only a few instances of interchromosomal rearrangement between the zebra finch and the chicken genomes. Interestingly, there was a greater than expected degree of intrachromosomal rearrangements compared to the chicken, suggesting that gene order is not conserved within avian chromosomes. At 1068 cM the map is approximately only one quarter the length of the chicken linkage map, providing further evidence that the chicken has an unusually high recombination rate. Male and female linkage-map lengths were similar, suggesting no heterochiasmy in the zebra finch. This whole-genome map is the first for any passerine and a valuable tool for the zebra finch genome sequence project and for studies of quantitative trait loci.

The chicken (Gallus gallus) Z chromosome contains at least three nonlinear evolutionary strata

Birds have female heterogamety with Z and W sex chromosomes. These evolved from different autosomal precursor chromosomes than the mammalian X and Y. However, previous work has suggested that the pattern and process of sex chromosome evolution show many similarities across distantly related organisms. Here we show that stepwise restriction of recombination between the protosex chromosomes of birds has resulted in regions of the chicken Z chromosome showing discrete levels of divergence from W homologs (gametologs). The 12 genes analyzed fall into three levels of estimated divergence values, with the most recent divergence (d(S) = 0.18-0.21) displayed by 6 genes in a region on the Z chromosome corresponding to the interval 1-11 Mb of the assembled genome sequence. Another 4 genes show intermediate divergence (d(S) = 0.27-0.38) and are located in the interval 16-53 Mb. Two genes (at positions 42 and 50 Mb) with higher d(S) values are located proximal to the most distal of the 4 genes with intermediate divergence, suggesting an inversion event. The distribution of genes and their divergence indicate at least three evolutionary strata, with estimated times for cessation of recombination between Z and W of 132-150 (stratum 1), 71-99 (stratum 2), and 47-57 (stratum 3) million years ago. An inversion event, or some other form of intrachromosomal rearrangement, subsequent to the formation of strata 1 and 2 has scrambled the gene order to give rise to the nonlinear arrangement of evolutionary strata currently seen on the chicken Z chromosome. These observations suggest that the progressive restriction of recombination is an integral feature of sex chromosome evolution and occurs also in systems of female heterogamety.

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

Pseudoautosomal regions (PARs) shared by avian Z and W sex chromosomes are typically small homologous regions within which recombination still occurs and are hypothesized to share the properties of autosomes. We capitalized on the unusual structure of the sex chromosomes of emus, Dromaius novaehollandiae, which consist almost entirely of PAR shared by both sex chromosomes, to test this hypothesis. We compared recombination, linkage disequilibrium (LD), GC content, and nucleotide diversity between pseudoautosomal and autosomal loci derived from 11 emu bacterial artificial chromosome (BAC) clones that were mapped to chromosomes by fluorescent in situ hybridization. Nucleotide diversity (pi = 4N(e)mu) was not significantly lower in pseudoautosomal loci (14 loci, 1.9 +/- 2.4 x 10(-3)) than autosomal loci (8 loci, 4.2 +/- 6.1 x 10(-3)). By contrast, recombination per site within BAC-end sequences (rho = 4Nc) (pseudoautosomal, 3.9 +/- 6.9 x 10(-2); autosomal, 2.3 +/- 3.7 x 10(-2)) was higher and average LD (D') (pseudoautosomal, 4.2 +/- 0.2 x 10(-1); autosomal, 4.7 +/- 0.5 x 10(-1)) slightly lower in pseudoautosomal sequences. We also report evidence of deviation from a simple neutral model in the PAR and in autosomal loci, possibly caused by departures from demographic equilibrium, such as population growth. This study provides a snapshot of the population genetics of avian sex chromosomes at an early stage of differentiation.

Relationship between physical and genetic distances along the zebra finch Z chromosome

Nine bacterial artificial chromosomes containing genes linked to the Z chromosome of the zebra finch (Taeniopygia guttata) were localized using FISH on synaptonemal complex spreads. Their positions were correlated with those previously reported on the mitotic Z chromosome, showing a linear relationship between positions along the mitotic chromosome and its synaptonemal complex. Distances in cM between the genes were calculated using a cytological map of the crossing-over based on the distribution of MLH1 foci along the ZZ synaptonemal complex (MLH1-cM map). It is shown that physical and genetic distances lack a linear relationship along most of the chromosome length, due to clustering of crossover events around the telomeres. This relationship departs strongly from that observed in the chicken Z chromosome and reflects the existence of different recombination rates and patterns among birds in spite of wide genomic conservation.

A gene-based genetic linkage map of the collared flycatcher (Ficedula albicollis) reveals extensive synteny and gene-order conservation during 100 million years of avian evolution

By taking advantage of a recently developed reference marker set for avian genome analysis we have constructed a gene-based genetic map of the collared flycatcher, an important "ecological model" for studies of life-history evolution, sexual selection, speciation, and quantitative genetics. A pedigree of 322 birds from a natural population was genotyped for 384 single nucleotide polymorphisms (SNPs) from 170 protein-coding genes and 71 microsatellites. Altogether, 147 gene markers and 64 microsatellites form 33 linkage groups with a total genetic distance of 1787 cM. Male recombination rates are, on average, 22% higher than female rates (total distance 1982 vs. 1627 cM). The ability to anchor the collared flycatcher map with the chicken genome via the gene-based SNPs revealed an extraordinary degree of both synteny and gene-order conservation during avian evolution. The great majority of chicken chromosomes correspond to a single linkage group in collared flycatchers, with only a few cases of inter- and intrachromosomal rearrangements. The rate of chromosomal diversification, fissions/fusions, and inversions combined is thus considerably lower in birds (0.05/MY) than in mammals (0.6-2.0/MY). A dearth of repeat elements, known to promote chromosomal breakage, in avian genomes may contribute to their stability. The degree of genome stability is likely to have important consequences for general evolutionary patterns and may explain, for example, the comparatively slow rate by which genetic incompatibility among lineages of birds evolves.

The chromosomal polymorphism linked to variation in social behavior in the white-throated sparrow (Zonotrichia albicollis) is a complex rearrangement and suppressor of recombination

Variation in social behavior and plumage in the white-throated sparrow (Zonotrichia albicollis) is linked to an inversion polymorphism on chromosome 2. Here we report the results of our comparative cytogenetic mapping efforts and population genetics studies focused on the genomic characterization of this balanced chromosomal polymorphism. Comparative chromosome painting and cytogenetic mapping of 15 zebra finch BAC clones to the standard (ZAL2) and alternative (ZAL2(m)) arrangements revealed that this chromosome is orthologous to chicken chromosome 3, and that at a minimum, ZAL2 and ZAL2(m) differ by a pair of included pericentric inversions that we estimate span at least 98 Mb. Population-based sequencing and genotyping of multiple loci demonstrated that ZAL2(m) suppresses recombination in the heterokaryotype and is evolving as a rare nonrecombining autosomal segment of the genome. In addition, we estimate that the first inversion within the ZAL2(m) arrangement originated 2.2+/-0.3 million years ago. Finally, while previously recognized as a genetic model for the evolution of social behavior, we found that the ZAL2/ZAL2(m) polymorphism also shares genetic and phenotypic features with the mouse t complex and we further suggest that the ZAL2/ZAL2(m) polymorphism is a heretofore unrecognized model for the early stages of sex chromosome evolution.

Genomics of natural bird populations: a gene-based set of reference markers evenly spread across the avian genome

Although there is growing interest to take genomics into the complex realms of natural populations, there is a general shortage of genomic resources and tools available for wild species. This applies not at least to birds, for which genomic approaches should be helpful to questions such as adaptation, speciation and population genetics. In this study, we describe a genome-wide reference set of conserved avian gene markers, broadly applicable across birds. By aligning protein-coding sequences from the recently assembled chicken genome with orthologous sequences in zebra finch, we identified particularly conserved exonic regions flanking introns of suitable size for subsequent amplification and sequencing. Primers were designed for 242 gene markers evenly distributed across the chicken genome, with a mean inter-marker interval of 4.2 Mb. Between 78% and 93% of the markers amplified a specific product in five species tested (chicken, peregrine falcon, collared flycatcher, great reed warbler and blue tit). Two hundred markers were sequenced in collared flycatcher, yielding a total of 122.41 kb of genomic DNA sequence (12096 bp coding sequence and 110 314 bp noncoding). Intron size of collared flycatcher and chicken was highly correlated, as was GC content. A polymorphism screening using these markers in a panel of 10 unrelated collared flycatchers identified 871 single nucleotide polymorphisms (pi = 0.0029) and 33 indels (mainly very short). Avian genome characteristics such as uniform genome size and low rate of syntenic rearrangements suggest that this marker set will find broad utility as a genome-wide reference resource for molecular ecological and population genomic analysis of birds. We envision that it will be particularly useful for obtaining large-scale orthologous targets in different species--important in, for instance, phylogenetics--and for large-scale identification of evenly distributed single nucleotide polymorphisms needed in linkage mapping or in studies of gene flow and hybridization.

Avian genomics in the 21st century

The chicken has long been an important model organism for developmental biology, as well as a major source of protein with billions of birds used in meat and egg production each year. Chicken genomics has been transformed in recent years, with the characterisation of large EST collections and most recently with the assembly of the chicken genome sequence. As the first livestock genome to be fully sequenced it leads the way for others to follow--with zebra finch later this year. The genome sequence and the availability of three million genetic polymorphisms are expected to aid the identification of genes that control traits of importance in poultry. As the first bird genome to be sequenced it is a model for the remaining 9,600 species thought to exist today. Many of the features of avian biology and organisation of the chicken genome make it an ideal model organism for phylogenetics and embryology, along with applications in agriculture and medicine. The availability of new tools such as whole-genome gene expression arrays and SNP panels, coupled with information resources on the genes and proteins are likely to enhance this position.