Molecular evolution of genes in avian genomes

Background

Obtaining a draft genome sequence of the zebra finch (Taeniopygia guttata), the second bird genome to be sequenced, provides the necessary resource for whole-genome comparative analysis of gene sequence evolution in a non-mammalian vertebrate lineage. To analyze basic molecular evolutionary processes during avian evolution, and to contrast these with the situation in mammals, we aligned the protein-coding sequences of 8,384 1:1 orthologs of chicken, zebra finch, a lizard and three mammalian species.

Results

We found clear differences in the substitution rate at fourfold degenerate sites, being lowest in the ancestral bird lineage, intermediate in the chicken lineage and highest in the zebra finch lineage, possibly reflecting differences in generation time. We identified positively selected and/or rapidly evolving genes in avian lineages and found an over-representation of several functional classes, including anion transporter activity, calcium ion binding, cell adhesion and microtubule cytoskeleton.

Conclusions

Focusing specifically on genes of neurological interest and genes differentially expressed in the unique vocal control nuclei of the songbird brain, we find a number of positively selected genes, including synaptic receptors. We found no evidence that selection for beneficial alleles is more efficient in regions of high recombination; in fact, there was a weak yet significant negative correlation between ω and recombination rate, which is in the direction predicted by the Hill-Robertson effect if slightly deleterious mutations contribute to protein evolution. These findings set the stage for studies of functional genetics of avian genes.     

Genetic mapping of the major histocompatibility complex in the zebra finch (<Emphasis Type="Italic">Taeniopygia guttata</Emphasis>)

Genes of the major histocompatibility complex (MHC) have received much attention in immunology, genetics, and ecology because they are highly polymorphic and play important roles in parasite resistance and mate choice. Until recently, the MHC of passerine birds was not well-described. However, the genome sequencing of the zebra finch (Taeniopygia guttata) has partially redressed this gap in our knowledge of avian MHC genes. Here, we contribute further to the understanding of the zebra finch MHC organization by mapping SNPs within or close to known MHC genes in the zebra finch genome. MHC class I and IIB genes were both mapped to zebra finch chromosome 16, and there was no evidence that MHC class I genes are located on chromosome 22 (as suggested by the genome assembly). We confirm the location in the MHC region on chromosome 16 for several other genes (BRD2, FLOT1, TRIM7.2, GNB2L1, and CSNK2B). Two of these (CSNK2B and FLOT1) have not previously been mapped in any other bird species. In line with previous results, we also find that orthologs to the immune-related genes B-NK and CLEC2D, which are part of the MHC region in chicken, are situated on zebra finch chromosome Z and not among other MHC genes in the zebra finch.     

Intrachromosomal rearrangements in avian genome evolution: evidence for regions prone to breakpoints

It is generally believed that the organization of avian genomes remains highly conserved in evolution as chromosome number is constant and comparative chromosome painting demonstrated there to be very few interchromosomal rearrangements. The recent sequencing of the zebra finch (Taeniopygia guttata) genome allowed an assessment of the number of intrachromosomal rearrangements between it and the chicken (Gallus gallus) genome, revealing a surprisingly high number of intrachromosomal rearrangements. With the publication of the turkey (Meleagris gallopavo) genome it has become possible to describe intrachromosomal rearrangements between these three important avian species, gain insight into the direction of evolutionary change and assess whether breakpoint regions are reused in birds. To this end, we aligned entire chromosomes between chicken, turkey and zebra finch, identifying syntenic blocks of at least 250 kb. Potential optimal pathways of rearrangements between each of the three genomes were determined, as was a potential Galliform ancestral organization. From this, our data suggest that around one-third of chromosomal breakpoint regions may recur during avian evolution, with 10% of breakpoints apparently recurring in different lineages. This agrees with our previous hypothesis that mechanisms of genome evolution are driven by hotspots of non-allelic homologous recombination.     

Reconstruction of gross avian genome structure,organization and evolution suggests that the chicken lineage most closely resembles the dinosaur avian ancestor

Background

The availability of multiple avian genome sequence assemblies greatly improves our ability to define overall genome organization and reconstruct evolutionary changes. In birds, this has previously been impeded by a near intractable karyotype and relied almost exclusively on comparative molecular cytogenetics of only the largest chromosomes. Here, novel whole genome sequence information from 21 avian genome sequences (most newly assembled) made available on an interactive browser (Evolution Highway) was analyzed.

Results

Focusing on the six best-assembled genomes allowed us to assemble a putative karyotype of the dinosaur ancestor for each chromosome. Reconstructing evolutionary events that led to each species’ genome organization, we determined that the fastest rate of change occurred in the zebra finch and budgerigar, consistent with rapid speciation events in the Passeriformes and Psittaciformes. Intra- and interchromosomal changes were explained most parsimoniously by a series of inversions and translocations respectively, with breakpoint reuse being commonplace. Analyzing chicken and zebra finch, we found little evidence to support the hypothesis of an association of evolutionary breakpoint regions with recombination hotspots but some evidence to support the hypothesis that microchromosomes largely represent conserved blocks of synteny in the majority of the 21 species analyzed. All but one species showed the expected number of microchromosomal rearrangements predicted by the haploid chromosome count. Ostrich, however, appeared to retain an overall karyotype structure of 2n = 80 despite undergoing a large number (26) of hitherto un-described interchromosomal changes.

Conclusions

Results suggest that mechanisms exist to preserve a static overall avian karyotype/genomic structure, including the microchromosomes, with widespread interchromosomal change occurring rarely (e.g., in ostrich and budgerigar lineages). Of the species analyzed, the chicken lineage appeared to have undergone the fewest changes compared to the dinosaur ancestor.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1060) contains supplementary material, which is available to authorized users.     

Cross-species overgo hybridization and comparative physical mapping within avian genomes

The chicken genome sequence facilitates comparative genomics within other avian species. We performed cross-species hybridizations using overgo probes designed from chicken genomic and zebra finch expressed sequence tags (ESTs) to turkey and zebra finch BAC libraries. As a result, 3772 turkey BACs were assigned to 336 markers or genes, and 1662 zebra finch BACs were assigned to 164 genes. As expected, cross-hybridization was more successful with overgos within coding sequences than within untranslated region, intron or flanking sequences and between chicken and turkey, when compared with chicken-zebra finch or zebra finch-turkey cross-hybridization. These data contribute to the comparative alignment of avian genome maps using a 'one sequence, multiple genomes' strategy.     

Disruption and pseudoautosomal localization of the major histocompatibility complex in monotremes

Background

The monotremes, represented by the duck-billed platypus and the echidnas, are the most divergent species within mammals, featuring a flamboyant mix of reptilian, mammalian and specialized characteristics. To understand the evolution of the mammalian major histocompatibility complex (MHC), the analysis of the monotreme genome is vital.

Results

We characterized several MHC containing bacterial artificial chromosome clones from platypus (Ornithorhynchus anatinus) and the short-beaked echidna (Tachyglossus aculeatus) and mapped them onto chromosomes. We discovered that the MHC of monotremes is not contiguous and locates within pseudoautosomal regions of two pairs of their sex chromosomes. The analysis revealed an MHC core region with class I and class II genes on platypus and echidna X3/Y3. Echidna X4/Y4 and platypus Y4/X5 showed synteny to the human distal class III region and beyond. We discovered an intron-containing class I pseudogene on platypus Y4/X5 at a genomic location equivalent to the human HLA-B,C region, suggesting ancestral synteny of the monotreme MHC. Analysis of male meioses from platypus and echidna showed that MHC chromosomes occupy different positions in the meiotic chains of either species.

Conclusion

Molecular and cytogenetic analyses reveal new insights into the evolution of the mammalian MHC and the multiple sex chromosome system of monotremes. In addition, our data establish the first homology link between chicken microchromosomes and the smallest chromosomes in the monotreme karyotype. Our results further suggest that segments of the monotreme MHC that now reside on separate chromosomes must once have been syntenic and that the complex sex chromosome system of monotremes is dynamic and still evolving.     

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.     

Reconstructing an ancestral mammalian immune supercomplex from a marsupial major histocompatibility complex

The first sequenced marsupial genome promises to reveal unparalleled insights into mammalian evolution. We have used theMonodelphis domestica (gray short-tailed opossum) sequence to construct the first map of a marsupial major histocompatibility complex (MHC). The MHC is the most gene-dense region of the mammalian genome and is critical to immunity and reproductive success. The marsupial MHC bridges the phylogenetic gap between the complex MHC of eutherian mammals and the minimal essential MHC of birds. Here we show that the opossum MHC is gene dense and complex, as in humans, but shares more organizational features with non-mammals. The Class I genes have amplified within the Class II region, resulting in a unique Class I/II region. We present a model of the organization of the MHC in ancestral mammals and its elaboration during mammalian evolution. The opossum genome, together with other extant genomes, reveals the existence of an ancestral “immune supercomplex” that contained genes of both types of natural killer receptors together with antigen processing genes and MHC genes.     

Genomic organization and molecular phylogenies of the beta (β) keratin multigene family in the chicken (<Emphasis Type="Italic">Gallus gallus</Emphasis>) and zebra finch (<Emphasis Type="Italic">Taeniopygia guttata</Emphasis>): implications for feather evolution

Background  

The epidermal appendages of reptiles and birds are constructed of beta (β) keratins. The molecular phylogeny of these keratins is important to understanding the evolutionary origin of these appendages, especially feathers. Knowing that the crocodilian β-keratin genes are closely related to those of birds, the published genomes of the chicken and zebra finch provide an opportunity not only to compare the genomic organization of their β-keratins, but to study their molecular evolution in archosaurians.     

Utilization of a zebra finch BAC library to determine the structure of an avian androgen receptor genomic region

The zebra finch (Taeniopygia guttata) is an important model organism for studying behavior, neuroscience, avian biology, and evolution. To support the study of its genome, we constructed a BAC library (TG__Ba) using DNA from livers of females. The BAC library consists of 147,456 clones with 98% containing inserts of an average size of 134 kb and represents 15.5 haploid genome equivalents. By sequencing a whole BAC, a full-length androgen receptor open reading frame was identified, the first in an avian species. Comparison of BAC end sequences and the whole BAC sequence with the chicken genome draft sequence showed a high degree of conserved synteny between the zebra finch and the chicken genome.     

Expression analysis of turkey (Meleagris gallopavo) toll-like receptors and molecular characterization of avian specific TLR15

Toll-like receptors (TLRs) constitute a multi-gene family, which plays a pivotal role in sensing invading pathogens by virtue of conserved microbial patterns. TLR repertoire of chicken and zebra finch has been well studied. However TLR family of other avian species is yet to be characterized. In the present study, we identified TLR repertoire of turkey, characterized avian specific receptor TLR15 in turkey and profiled the TLRs expressions in a range of tissues of turkey poults. All ten TLR genes orthologous to chicken TLR repertoire were found in turkey. Turkey TLR genes showed 81-93 % similarity at amino acid level to their chicken counter parts. Phylogenetic analysis confirmed the orthologous relationship of turkey TLRs with chicken and zebra finch TLRs. Open reading frame of turkey TLR15 was 2,607 bp long encoding 868 amino acids similar to that of broiler chicken and showed 92.4, 91.1 and 69.5 % identity at amino acid levels with chicken, Japanese quail and zebra finch TLR15 sequences respectively. Overall TLR expression was highest for TLR4 and lowest for TLR21. TLR1A, 2A, 2B and 21 were significantly higher in liver than other tissues investigated (P < 0.01). TLR3 expression was significantly higher in bone marrow (BM) and spleen in comparison to other tissues studied (P < 0.01). Furthermore, no significant differences in the expression levels of TLR1B, 4, 5, 7 and 15 genes were detected among the tissues studied. Our findings contribute to the characterization of innate immune system of birds and show the innate preparedness of young turkey poults to a range of pathogens.     

The zebra finch neuropeptidome: prediction,detection and expression

Background  

Among songbirds, the zebra finch (Taeniopygia guttata) is an excellent model system for investigating the neural mechanisms underlying complex behaviours such as vocal communication, learning and social interactions. Neuropeptides and peptide hormones are cell-to-cell signalling molecules known to mediate similar behaviours in other animals. However, in the zebra finch, this information is limited. With the newly-released zebra finch genome as a foundation, we combined bioinformatics, mass-spectrometry (MS)-enabled peptidomics and molecular techniques to identify the complete suite of neuropeptide prohormones and final peptide products and their distributions.     

Comparative genomics reveals molecular features unique to the songbird lineage

Background

Songbirds (oscine Passeriformes) are among the most diverse and successful vertebrate groups, comprising almost half of all known bird species. Identifying the genomic innovations that might be associated with this success, as well as with characteristic songbird traits such as vocal learning and the brain circuits that underlie this behavior, has proven difficult, in part due to the small number of avian genomes available until recently. Here we performed a comparative analysis of 48 avian genomes to identify genomic features that are unique to songbirds, as well as an initial assessment of function by investigating their tissue distribution and predicted protein domain structure.

Results

Using BLAT alignments and gene synteny analysis, we curated a large set of Ensembl gene models that were annotated as novel or duplicated in the most commonly studied songbird, the Zebra finch (Taeniopygia guttata), and then extended this analysis to 47 additional avian and 4 non-avian genomes. We identified 10 novel genes uniquely present in songbird genomes. A refined map of chromosomal synteny disruptions in the Zebra finch genome revealed that the majority of these novel genes localized to regions of genomic instability associated with apparent chromosomal breakpoints. Analyses of in situ hybridization and RNA-seq data revealed that a subset of songbird-unique genes is expressed in the brain and/or other tissues, and that 2 of these (YTHDC2L1 and TMRA) are highly differentially expressed in vocal learning-associated nuclei relative to the rest of the brain.

Conclusions

Our study reveals novel genes unique to songbirds, including some that may subserve their unique vocal control system, substantially improves the quality of Zebra finch genome annotations, and contributes to a better understanding of how genomic features may have evolved in conjunction with the emergence of the songbird lineage.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1082) contains supplementary material, which is available to authorized users.     

Towards precise classification of cancers based on robust gene functional expression profiles

Background

Despite the continuous production of genome sequence for a number of organisms, reliable, comprehensive, and cost effective gene prediction remains problematic. This is particularly true for genomes for which there is not a large collection of known gene sequences, such as the recently published chicken genome. We used the chicken sequence to test comparative and homology-based gene-finding methods followed by experimental validation as an effective genome annotation method.

Results

We performed experimental evaluation by RT-PCR of three different computational gene finders, Ensembl, SGP2 and TWINSCAN, applied to the chicken genome. A Venn diagram was computed and each component of it was evaluated. The results showed that de novo comparative methods can identify up to about 700 chicken genes with no previous evidence of expression, and can correctly extend about 40% of homology-based predictions at the 5' end.

Conclusions

De novo comparative gene prediction followed by experimental verification is effective at enhancing the annotation of the newly sequenced genomes provided by standard homology-based methods.     

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.     

DNA repeat arrays in chicken and human genomes and the adaptive evolution of avian genome size

Background

Birds have smaller average genome sizes than other tetrapod classes, and it has been proposed that a relatively low frequency of repeating DNA is one factor in reduction of avian genome sizes.

Results

DNA repeat arrays in the sequenced portion of the chicken (Gallus gallus) autosomes were quantified and compared with those in human autosomes. In the chicken 10.3% of the genome was occupied by DNA repeats, in contrast to 44.9% in human. In the chicken, the percentage of a chromosome occupied by repeats was positively correlated with chromosome length, but even the largest chicken chromosomes had repeat densities much lower than those in human, indicating that avoidance of repeats in the chicken is not confined to minichromosomes. When 294 simple sequence repeat types shared between chicken and human genomes were compared, mean repeat array length and maximum repeat array length were significantly lower in the chicken than in human.

Conclusions

The fact that the chicken simple sequence repeat arrays were consistently smaller than arrays of the same type in human is evidence that the reduction in repeat array length in the chicken has involved numerous independent evolutionary events. This implies that reduction of DNA repeats in birds is the result of adaptive evolution. Reduction of DNA repeats on minichromosomes may be an adaptation to permit chiasma formation and alignment of small chromosomes. However, the fact that repeat array lengths are consistently reduced on the largest chicken chromosomes supports the hypothesis that other selective factors are at work, presumably related to the reduction of cell size and consequent advantages for the energetic demands of flight.     

Genomic organization of zebra finch alpha and beta globin genes and their expression in primitive and definitive blood in comparison with globins in chicken

How alpha and beta globin genes are organized and expressed in amniotes is of interest to researchers in a wide variety of fields. Data regarding this from avian species have been scarce. Using genomic and proteomic approaches, we present here our analysis of alpha and beta globins of zebra finch, a passerine bird. We show that finch alpha globin gene cluster has three genes (alphas 1–3), each orthologous to its chicken counterpart. Finch beta globin gene cluster has three genes (betas 1–3), with an additional pseudogene at the 3′ end. Finch beta3 is orthologous to chicken betaA, but the orthology of beta1 and beta2 to chicken counterparts is less clear. All six finch globins are confirmed to encode functional proteins. Gene expression in both globin gene clusters is regulated developmentally. Adult finch blood has a globin profile similar to that of adult chicken, with high levels of beta3 and alpha3 and moderate levels of alpha2. Finch embryonic primitive blood exhibits a globin profile very different from that of equivalent stage chick embryos, with all six globins expressed at high levels. Overall, our data provide a valuable resource for future studies in avian globin gene evolution and globin switching during erythropoietic development.     

家禽MHC结构研究进展

禽主要组织相容性复合体(Major histocompatibility complex,MHC)的结构与禽病防控、禽免疫学、禽类遗传学研究密切相关。文章对鸡、火鸡、鹌鹑、鸭和鹅的MHC结构方面的研究进展进行了综述,表明其有以下共同特点:都有保守的MHC区域,包括MHC I基因和MHC II基因及一些功能未知基因;基因排列简单而紧凑;MHC I基因内含子的长度都比哺乳动物小;鸡、火鸡、鸭和鹅的MHC I基因组序列都有8个外显子和7个内含子,MHC IIβ基因组序列都有6个外显子和5个内含子;鸡、火鸡和鹌鹑的BG基因结构模式相同;都存在微卫星重复单元。但也存在种属差异:鸡的MHC I基因和MHC II基因是双拷贝,而鸭、鹅和鹌鹑有若干个拷贝;BG基因的拷贝数及其外显子数目不同。对主要家禽MHC结构进行分析比较,将有利于对禽病学及禽免疫遗传学的进究。