Songbird genomics: analysis of 45 kb upstream of a polymorphic Mhc class II gene in red-winged blackbirds (Agelaius phoeniceus)

Here we present the sequence of a 45 kb cosmid containing a previously characterized poly-morphic Mhc class II B gene (Agph-DAB1) from the red-winged blackbird (Agelaius phoeniceus). We compared it with a previously sequenced cosmid from this species, revealing two regions of 7.5 kb and 13.0 kb that averaged greater than 97% similarity to each another, indicating a very recent shared duplication. We found 12 retroelements, including two chicken repeat 1 (CR1) elements, constituting 6.4% of the sequence and indicating a lower frequency of retroelements than that found in mammalian genomic DNA. Agph-DAB3, a new class II B gene discovered in the cosmid, showed a low rate of polymorphism and may be functional. In addition, we found a Mhc class II B gene fragment and three genes likely to be functional (encoding activin receptor type II, a zinc finger, and a putative gamma-filamin). Phylogenetic analysis of exon 2 alleles of all three known blackbird Mhc genes indicated strong clustering of alleles by locus, implying that large amounts of interlocus gene conversion have not occurred since these genes have been diverging. Despite this, interspecific comparisons indicate that all three blackbird Mhc genes diverged from one another less than 35 million years ago and are subject to concerted evolution in the long term. Comparison of blackbird and chicken Mhc promoter regions revealed songbird promoter elements for the first time. The high gene density of this cosmid confirms similar findings for the chicken Mhc, but the segment duplications and diversity of retroelements resembles mammalian sequences.     

Songbird Genomics: Analysis of 45 kb Upstream of a Polymorphic Mhc Class II Gene in Red-Winged Blackbirds (Agelaius phoeniceus)

Here we present the sequence of a 45 kb cosmid containing a previously characterized poly-morphic Mhc class II B gene (Agph-DAB1) from the red-winged blackbird (Agelaius phoeniceus). We compared it with a previously sequenced cosmid from this species, revealing two regions of 7.5 kb and 13.0 kb that averaged greater than 97% similarity to each another, indicating a very recent shared duplication. We found 12 retroelements, including two chicken repeat 1 (CR1) elements, constituting 6.4% of the sequence and indicating a lower frequency of retroelements than that found in mammalian genomic DNA. Agph-DAB3, a new class II B gene discovered in the cosmid, showed a low rate of polymorphism and may be functional. In addition, we found a Mhc class II B gene fragment and three genes likely to be functional (encoding activin receptor type II, a zinc finger, and a putative γ-filamin). Phylogenetic analysis of exon 2 alleles of all three known blackbird Mhc genes indicated strong clustering of alleles by locus, implying that large amounts of interlocus gene conversion have not occurred since these genes have been diverging. Despite this, interspecific comparisons indicate that all three blackbird Mhc genes diverged from one another less than 35 million years ago and are subject to concerted evolution in the long term. Comparison of blackbird and chicken Mhc promoter regions revealed songbird promoter elements for the first time. The high gene density of this cosmid confirms similar findings for the chicken Mhc, but the segment duplications and diversity of retroelements resembles mammalian sequences.     

Comparative genomic analysis of two avian (quail and chicken) MHC regions

We mapped two different quail Mhc haplotypes and sequenced one of them (haplotype A) for comparative genomic analysis with a previously sequenced haplotype of the chicken Mhc. The quail haplotype A spans 180 kb of genomic sequence, encoding a total of 41 genes compared with only 19 genes within the 92-kb chicken Mhc. Except for two gene families (B30 and tRNA), both species have the same basic set of gene family members that were previously described in the chicken "minimal essential" Mhc. The two Mhc regions have a similar overall organization but differ markedly in that the quail has an expanded number of duplicated genes with 7 class I, 10 class IIB, 4 NK, 6 lectin, and 8 B-G genes. Comparisons between the quail and chicken Mhc class I and class II gene sequences by phylogenetic analysis showed that they were more closely related within species than between species, suggesting that the quail Mhc genes were duplicated after the separation of these two species from their common ancestor. The proteins encoded by the NK and class I genes are known to interact as ligands and receptors, but unlike in the quail and the chicken, the genes encoding these proteins in mammals are found on different chromosomes. The finding of NK-like genes in the quail Mhc strongly suggests an evolutionary connection between the NK C-type lectin-like superfamily and the Mhc, providing support for future studies on the NK, lectin, class I, and class II interaction in birds.     

Mhc diversity in two passerine birds: no evidence for a minimal essential Mhc

Humans express an array of Mhc genes, while the chicken has an Mhc that is relatively small and compact with fewer expressed genes. Here we ask whether the "minimal essential Mhc" of the chicken is representative for birds. We investigated the RFLP genotypes in 55 great reed warblers Acrocephalus arundinaceus and 10 willow warblers Phylloscopus trochilus to obtain an overview of the number of class II B genes. There were 13-17 bands per individual in the great reed warblers and 25-30 in the willow warblers, and every individual had a unique RFLP genotype. The high number of RFLP bands indicates that both species have a large number of class II B genes although some may be pseudogenes. Seven different class II B sequences were detected in a great reed warbler cDNA library. There was considerable sequence divergence between the cDNA sequences in exon 2 (peptide-binding region, PBR), whereas they were very similar in exon 3. The cDNA sequences were easily alignable to a classical chicken class II B sequence, and balancing selection was acting in the PBR. One of the cDNA sequences had two deletions and is likely nonfunctional. Finally, the polymorphic class I and class II B RFLP fragments seemed to be linked in the five studied great reed warbler families. These and previous results suggest that birds of the order Passeriformes in general have more Mhc class I and II B genes than birds of the order Galliformes. This difference could be caused by their phylogenetic past, and/or by variance in the selection pressure for maintaining a high number of Mhc genes.     

Ancestral Polymorphism of Mhc Class II Genes in Mice: Implications for Balancing Selection and the Mammalian Molecular Clock

To investigate the evolutionary dynamics at Mhc class II DR genes of mice (genus Mus), we sequenced the peptide binding regions (PBRs) of 41 DRB (=Eβ) genes and eight DRA (=Eα) genes from 15 strains representing eight species. As expected trees of these PBR sequences imply extensive maintenance of ancestral DRB alleles across species. We use a coalescent simulation model to show that the number of interspecific coalescent events (c) observed on these trees was higher than the number expected for neutral genealogies and similar sample sizes and is more consistent with balancing selection than with neutrality. Patterns of ancestral polymorphism in mouse DRB alleles were also used to examine the tempo of synonymous substitution in the PBR of mouse class II genes. Both absolute and relative rate tests on DRA and DRB genes imply increased substitution rates at two- and fourfold degnerate sites of mice and rats relative to primates, and decreased rates for the DRB genes of primates relative to ungulate and carnivore relatives. Thus rates of synonymous substitution at Mhc DR genes in mammals appear to be subject to generation time effects in ways similar to those found at other mammalian genes.     

The chicken genome contains no HMG1 retropseudogenes but a functional HMG1 gene with long introns

We have cloned the genomic sequence coding for the high mobility group 1 (HMG1) protein in chickens. Multiple sequence alignment shows that the chicken HMG1 gene is highly homologous to the human and the mouse HMG1 genes. The gene structure of chicken HMG1 is similar to that of the mouse and the human HMG1 genes, with the same exon-intron boundaries. However, in contrast to other avian genes that have shorter introns, the chicken HMG1 gene has introns that are twice as long as their mammalian homologues. In addition to the functional, intron-containing HMG1 gene, all mammalian genomes contain more than 50 copies of HMG1 retropseudogenes each, while in the chicken genome there are no HMG1 retropseudogenes. This finding suggests that the HMG1 retropseudogenes arose in mammals after their divergence away from the birds.     

Isolation of chicken major histocompatibility complex class II (B-L) beta chain sequences: comparison with mammalian beta chains and expression in lymphoid organs

By cross-hybridization in low stringency conditions, using a probe derived from an HLA-DQ beta cDNA clone, we have isolated several chicken genomic DNA clones. These clones were mapped to the major histocompatibility complex (MHC) of the chick (B complex) by virtue of their ability to detect restriction enzyme length polymorphisms between congenic lines of chicken. Evidence was obtained for the presence of at least three B-L beta genes in the chicken genome. The B-L beta genes are transcribed specifically in tissues containing cells of the B lymphocyte and myeloid lineages and expressing the B-L antigens. Exons encoding the beta 1, beta 2 and transmembrane domains of a B-L beta chain have been identified with 63, 66 and 62% similarity with the HLA-DQ beta sequence. This first isolation of an MHC class II gene outside of the mammalian class provides insight into the evolution of MHC genes based on the comparison of avian and mammalian class II beta chain amino acid and nucleotide sequences.     

Sequence and gene organization of the chicken mitochondrial genome. A novel gene order in higher vertebrates

The 16,775 base-pair mitochondrial genome of the white Leghorn chicken has been cloned and sequenced. The avian genome encodes the same set of genes (13 proteins, 2 rRNAs and 22 tRNAs) as do other vertebrate mitochondrial DNAs and is organized in a very similar economical fashion. There are very few intergenic nucleotides and several instances of overlaps between protein or tRNA genes. The protein genes are highly similar to their mammalian and amphibian counterparts and are translated according to the same variant genetic code. Despite these highly conserved features, the chicken mitochondrial genome displays two distinctive characteristics. First, it exhibits a novel gene order, the contiguous tRNA(Glu) and ND6 genes are located immediately adjacent to the displacement loop region of the molecule, just ahead of the contiguous tRNA(Pro), tRNA(Thr) and cytochrome b genes, which border the displacement loop region in other vertebrate mitochondrial genomes. This unusual gene order is conserved among the galliform birds. Second, a light-strand replication origin, equivalent to the conserved sequence found between the tRNA(Cys) and tRNA(Asn) genes in all vertebrate mitochondrial genomes sequenced thus far, is absent in the chicken genome. These observations indicate that galliform mitochondrial genomes departed from their mammalian and amphibian counterparts during the course of evolution of vertebrate species. These unexpected characteristics represent useful markers for investigating phylogenetic relationships at a higher taxonomic level.     

The chicken gene nomenclature committee report

Comparative genomics is an essential component of the post-genomic era. The chicken genome is the first avian genome to be sequenced and it will serve as a model for other avian species. Moreover, due to its unique evolutionary niche, the chicken genome can be used to understand evolution of functional elements and gene regulation in mammalian species. However comparative biology both within avian species and within amniotes is hampered due to the difficulty of recognising functional orthologs. This problem is compounded as different databases and sequence repositories proliferate and the names they assign to functional elements proliferate along with them. Currently, genes can be published under more than one name and one name sometimes refers to unrelated genes. Standardized gene nomenclature is necessary to facilitate communication between scientists and genomic resources. Moreover, it is important that this nomenclature be based on existing nomenclature efforts where possible to truly facilitate studies between different species. We report here the formation of the Chicken Gene Nomenclature Committee (CGNC), an international and centralized effort to provide standardized nomenclature for chicken genes. The CGNC works in conjunction with public resources such as NCBI and Ensembl and in consultation with existing nomenclature committees for human and mouse. The CGNC will develop standardized nomenclature in consultation with the research community and relies on the support of the research community to ensure that the nomenclature facilitates comparative and genomic studies.     

Natural selection of the major histocompatibility complex (Mhc) in Hawaiian honeycreepers (Drepanidinae)

The native Hawaiian honeycreepers represent a classic example of adaptive radiation and speciation, but currently face one the highest extinction rates in the world. Although multiple factors have likely influenced the fate of Hawaiian birds, the relatively recent introduction of avian malaria is thought to be a major factor limiting honeycreeper distribution and abundance. We have initiated genetic analyses of class II beta chain Mhc genes in four species of honeycreepers using methods that eliminate the possibility of sequencing mosaic variants formed by cloning heteroduplexed polymerase chain reaction products. Phylogenetic analyses group the honeycreeper Mhc sequences into two distinct clusters. Variation within one cluster is high, with dN > dS and levels of diversity similar to other studies of Mhc (B system) genes in birds. The second cluster is nearly invariant and includes sequences from honeycreepers (Fringillidae), a sparrow (Emberizidae) and a blackbird (Emberizidae). This highly conserved cluster appears reminiscent of the independently segregating Rfp-Y system of genes defined in chickens. The notion that balancing selection operates at the Mhc in the honeycreepers is supported by transpecies polymorphism and strikingly high dN/dS ratios at codons putatively involved in peptide interaction. Mitochondrial DNA control region sequences were invariant in the i'iwi, but were highly variable in the 'amakihi. By contrast, levels of variability of class II beta chain Mhc sequence codons that are hypothesized to be directly involved in peptide interactions appear comparable between i'iwi and 'amakihi. In the i'iwi, natural selection may have maintained variation within the Mhc, even in the face of what appears to a genetic bottleneck.     

Concerted evolution within a trypsin gene cluster in Drosophila.

A cluster of four trypsin genes has previously been localized to cytological position 47D-F of the Drosophila melanogaster genome. One of these genes had been sequenced, and the presence of the other three genes was identified by cross-hybridization. Here, we present the DNA sequence of the entire genomic region encoding these four trypsin genes. In addition to the four previously inferred genes, we have identified a fifth trypsin-coding sequence located within this gene cluster. This new gene shows a high degree of sequence divergence (more than 30%) from the other four genes, although it retains all of the functional motifs that are characteristic of trypsin-coding sequences. In order to trace the molecular evolution of this gene cluster, we isolated and sequenced the homologous 7-kb region from the closely related species Drosophila erecta. A comparison of the DNA sequences between the two species provides strong evidence for the concerted evolution of some members of this gene family. Two genes within the cluster are evolving in concert, while a third gene appears to be evolving independently. The remaining two genes show an intermediate pattern of evolution. We propose a simple model, involving chromosome looping and gene conversion, to explain the relatively complex patterns of molecular evolution within this gene cluster.     

Characterization of the first growth hormone gene sequence for a passerine bird--the pied flycatcher (Ficedula hypoleuca).

While the growth hormone (GH) gene has been characterized in a broad range of vertebrates, surprisingly little is known about this gene in birds. In order to extend knowledge of the GH gene in avian species and non-domestic species, the pied flycatcher (Ficedula hypoleuca) GH gene has been sequenced in this study. The overall average pairwise sequence divergence level was 0.08 among all available avian sequences and 0.27 among other taxa. However, the overall genetic organization of the gene is quite conserved. The similarity of the GH gene sequence of pied flycatchers with those of chicken and duck suggests that the rapid bursts of molecular evolution observed in mammalian and fish GH have not occurred during the divergence of passerine and non-passerine birds.     

Multi-Platform Next-Generation Sequencing of the Domestic Turkey (Meleagris gallopavo): Genome Assembly and Analysis

A synergistic combination of two next-generation sequencing platforms with a detailed comparative BAC physical contig map provided a cost-effective assembly of the genome sequence of the domestic turkey (Meleagris gallopavo). Heterozygosity of the sequenced source genome allowed discovery of more than 600,000 high quality single nucleotide variants. Despite this heterozygosity, the current genome assembly (∼1.1 Gb) includes 917 Mb of sequence assigned to specific turkey chromosomes. Annotation identified nearly 16,000 genes, with 15,093 recognized as protein coding and 611 as non-coding RNA genes. Comparative analysis of the turkey, chicken, and zebra finch genomes, and comparing avian to mammalian species, supports the characteristic stability of avian genomes and identifies genes unique to the avian lineage. Clear differences are seen in number and variety of genes of the avian immune system where expansions and novel genes are less frequent than examples of gene loss. The turkey genome sequence provides resources to further understand the evolution of vertebrate genomes and genetic variation underlying economically important quantitative traits in poultry. This integrated approach may be a model for providing both gene and chromosome level assemblies of other species with agricultural, ecological, and evolutionary interest.     

A Revised Evolutionary History of the CYP1A Subfamily: Gene Duplication, Gene Conversion, and Positive Selection

Members of cytochrome P450 subfamily 1A (CYP1As) are involved in detoxification and bioactivation of common environmental pollutants. Understanding the functional evolution of these genes is essential to predicting and interpreting species differences in sensitivity to toxicity caused by such chemicals. The CYP1A gene subfamily comprises a single ancestral representative in most fish species and two paralogs in higher vertebrates, including birds and mammals. Phylogenetic analysis of complete coding sequences suggests that mammalian and bird paralog pairs (CYP1A1/2 and CYP1A4/5, respectively) are the result of independent gene duplication events. However, comparison of vertebrate genome sequences revealed that CYP1A genes lie within an extended region of conserved fine-scale synteny, suggesting that avian and mammalian CYP1A paralogs share a common genomic history. Algorithms designed to detect recombination between nucleotide sequences indicate that gene conversion has homogenized most of the length of the chicken CYP1A genes, as well as the 5′ end of mammalian CYP1As. Together, these data indicate that avian and mammalian CYP1A paralog pairs resulted from a single gene duplication event and that extensive gene conversion is responsible for the exceptionally high degree of sequence similarity between CYP1A4 and CYP1A5. Elevated nonsynonymous/synonymous substitution ratios within a putatively unconverted stretch of ∼250 bp suggests that positive selection may have reduced the effective rate of gene conversion in this region, which contains two substrate recognition sites. This work significantly alters our understanding of functional evolution in the CYP1A subfamily, suggesting that gene conversion and positive selection have been the dominant processes of sequence evolution. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Yves Van de Peer]     

Evolutionary constraint in flanking regions of avian genes

An important comprehension from comparative genomic analysis is that sequence conservation beyond neutral expectations is frequently found outside protein-coding regions, indicating important functional roles of noncoding DNA. Understanding the causes of constraint on noncoding sequence evolution forms an important area of research, not least in light of the importance for understanding the evolution of gene expression. We aligned all orthologous genes of chicken and zebra finch together with 5 kb of their upstream and downstream noncoding sequences, to study the evolution of gene flanking sequences in the avian genome. Using ancestral repeats as a neutral reference, we detected significant evolutionary constraint in the 3' flanking region, highest directly after termination (60%) and then gradually decreasing to about 20% 5 kb downstream. Constraint was higher in annotated 3' untranslated regions (UTRs) than in non-UTRs at the same distance from the stop codon and higher in sequences annotated as microRNA (miRNA)-binding sites than in non-miRNA-binding sites within 3' UTRs. Constraint was also higher when estimated for a smaller data set of genes from more closely related songbird species, indicating turnover of functional elements during avian evolution. On the 5' flanking side constraint was readily seen within the first 125 bp immediately upstream of the start codon (34%) and was about 10% for remaining sequence within 5 kb upstream. Analysis of chicken polymorphism data gave further support for the highest constraint directly before and after the translated region. Finally, we found that genes evolving under the highest constraint measured by d(N)/d(S) also had the highest level of constraint in the 3' flanking region. This study broadens the insights into gene flanking sequence evolution by adding new findings from a vertebrate lineage other than mammals.     

Genetic mapping in a natural population of collared flycatchers (Ficedula albicollis): conserved synteny but gene order rearrangements on the avian Z chromosome

Data from completely sequenced genomes are likely to open the way for novel studies of the genetics of nonmodel organisms, in particular when it comes to the identification and analysis of genes responsible for traits that are under selection in natural populations. Here we use the draft sequence of the chicken genome as a starting point for linkage mapping in a wild bird species, the collared flycatcher - one of the most well-studied avian species in ecological and evolutionary research. A pedigree of 365 flycatchers was established and genotyped for single nucleotide polymorphisms in 23 genes selected from (and spread over most of) the chicken Z chromosome. All genes were also found to be located on the Z chromosome in the collared flycatcher, confirming conserved synteny at the level of gene content across distantly related avian lineages. This high degree of conservation mimics the situation seen for the mammalian X chromosome and may thus be a general feature in sex chromosome evolution, irrespective of whether there is male or female heterogamety. Alternatively, such unprecedented chromosomal conservation may be characteristic of most chromosomes in avian genome evolution. However, several internal rearrangements were observed, meaning that the transfer of map information from chicken to nonmodel bird species cannot always assume conserved gene orders. Interestingly, the rate of recombination on the Z chromosome of collared flycatchers was only approximately 50% that of chicken, challenging the widely held view that birds generally have high recombination rates.     

Polymorphism across an exon-intron boundary in an avian Mhc class II B gene.

Twenty-three sequence haplotypes spanning the boundary of the second exon and intron of a red-winged blackbird Mhc class II B gene, Agph-DAB1, are presented. The polymorphism of the exon segment is distributed in two divergent allelic lineages which appear to be maintained by balancing selection. The silent nucleotide diversity of the exon (pi = 0.101) is more than five times that of the intron (pi = 0.018) and decays rapidly across the exon-intron boundary. Additionally, genealogical reconstruction indicates that divergence from a common ancestor in the exon sample is over four times that of the intron. The intron sequences reveal a pattern of polymorphism which is characteristic of directional selection, rather than a pattern expected from linkage to a balanced polymorphism. These results suggest that the evolutionary histories of these two adjacent regions have been disassociated by recombination or gene conversion. The estimated population recombination parameter between the exon and the intron is sufficiently high (4NeC = 8.545) to explain the homogenization of intron sequences. Compatibility analyses estimate that these events primarily occur from the exon-intron boundary to about 20-30 bases into the intron. Additionally, the observation that divergent exon alleles share identical intron sequence supports the conclusion of disassociation of exon and intron evolutionary histories by recombination.     

Analysis of the genomic structure of the porcine CD1 gene cluster

CD1 is an MHC class I-like protein that presents lipid antigens to T cell receptors. We determined 470,187 bp of the genomic sequence encompassing the region encoding porcine CD1 genes. We identified 16 genes in this region and newly identified CD1A2, CD1B, CD1C, CD1D, and CD1E. Porcine CD1 genes were located in clusters between KIRREL and olfactory receptor (OR) genes, as observed in humans, although they were divided into two regions by a region encoding OR genes. Comparison of the genomic sequences of CD1 gene loci in pigs with other mammals showed that separation of the CD1 gene cluster by ORs was observed only in pigs. CD1A duplication in the porcine genome was estimated to have occurred after the divergence of the human and porcine. This analysis of the genomic sequence of the porcine CD1 family will contribute to our understanding of the evolution of mammalian CD1 genes.     

The mtDNA sequence of the ostrich and the divergence between paleognathous and neognathous birds

The complete mitochondrial DNA (mtDNA) molecule of the ostrich, Struthio camelus, was sequenced. The size of the molecule is 16,591 nucleotides. Since the ostrich represents the paleognathous birds, comparison with the mtDNA of the neognathous chicken, the only avian species reported so far in databases, made it possible to identify common and, probably, general avian mtDNA characteristics. Relative to other vertebrates, the avian NADH6 and tRNA-Glu genes are positioned upstream of the control region rather than the cytochrome b gene. The NADH3 gene of the ostrich is terminated by a stop codon at position 207. Thus, the gene is about 140 nucleotides shorter than in other vertebrates. The sequence for L-strand origin of replication is missing in both birds, and four transfer RNA genes of the two avian mtDNAs deviate from common characteristics of tRNAs of vertebrate mtDNAs by having an adenine (and not a thymidine) at position 8. Due to the absence of suitable fossils, most paleontological datings of avian divergences are conjectural. Molecular dating of the divergence between the ostrich and the chicken indicates that these two avian lineages separated 80-90 MYA. Phylogenetic analysis of complete cytochrome b genes of six avian orders showed that Passeriformes represent the earliest divergence among recent birds, contradicting the commonly accepted notion of a basal position of the Palaeognathae among recent birds.