Christopher N. Balakrishnan

The genome of a songbird.

The zebra finch is an important model organism in several fields with unique relevance to human neuroscience. Like other songbirds, the zebra finch communicates through learned vocalizations, an ability otherwise documented only in humans and a few other animals and lacking in the chicken—the only bird with a sequenced genome until now. Here we present a structural, functional and comparative analysis of the genome sequence of the zebra finch (Taeniopygia guttata), which is a songbird belonging to the large avian order Passeriformes. We find that the overall structures of the genomes are similar in zebra finch and chicken, but they differ in many intrachromosomal rearrangements, lineage-specific gene family expansions, the number of long-terminal-repeat-based retrotransposons, and mechanisms of sex chromosome dosage compensation. We show that song behaviour engages gene regulatory networks in the zebra finch brain, altering the expression of long non-coding RNAs, microRNAs, transcription factors and their targets. We also show evidence for rapid molecular evolution in the songbird lineage of genes that are regulated during song experience. These results indicate an active involvement of the genome in neural processes underlying vocal communication and identify potential genetic substrates for the evolution and regulation of this behaviour.

Gene duplication and fragmentation in the zebra finch major histocompatibility complex

BackgroundDue to its high polymorphism and importance for disease resistance, the major histocompatibility complex (MHC) has been an important focus of many vertebrate genome projects. Avian MHC organization is of particular interest because the chicken Gallus gallus, the avian species with the best characterized MHC, possesses a highly streamlined minimal essential MHC, which is linked to resistance against specific pathogens. It remains unclear the extent to which this organization describes the situation in other birds and whether it represents a derived or ancestral condition. The sequencing of the zebra finch Taeniopygia guttata genome, in combination with targeted bacterial artificial chromosome (BAC) sequencing, has allowed us to characterize an MHC from a highly divergent and diverse avian lineage, the passerines.ResultsThe zebra finch MHC exhibits a complex structure and history involving gene duplication and fragmentation. The zebra finch MHC includes multiple Class I and Class II genes, some of which appear to be pseudogenes, and spans a much more extensive genomic region than the chicken MHC, as evidenced by the presence of MHC genes on each of seven BACs spanning 739 kb. Cytogenetic (FISH) evidence and the genome assembly itself place core MHC genes on as many as four chromosomes with TAP and Class I genes mapping to different chromosomes. MHC Class II regions are further characterized by high endogenous retroviral content. Lastly, we find strong evidence of selection acting on sites within passerine MHC Class I and Class II genes.ConclusionThe zebra finch MHC differs markedly from that of the chicken, the only other bird species with a complete genome sequence. The apparent lack of synteny between TAP and the expressed MHC Class I locus is in fact reminiscent of a pattern seen in some mammalian lineages and may represent convergent evolution. Our analyses of the zebra finch MHC suggest a complex history involving chromosomal fission, gene duplication and translocation in the history of the MHC in birds, and highlight striking differences in MHC structure and organization among avian lineages.

Stable recombination hotspots in birds

Recombination: The birds and the yeast Apes and mice have a specific gene, PRDM9, that is associated with genomic regions with high rates of recombination, called hotspots. In species with PRDM9, hotspots move rapidly within the genome, varying among populations and closely related species (see the Perspective by Lichten). To investigate recombination hotspots in species lacking PRDM9, Singhal et al. examined bird genomes, which lack a PRDM9 gene. They looked closely at the genomes of finch species and found that recombination was localized to the promoter regions of genes and highly conserved over millions of years. Similarly, Lam and Keeney examined recombination localization within yeast, which also lacks PRDM9. They found a similar more-or-less fixed pattern of hotspots. Thus, recombination in species lacking a PRDM9 gene shows similar patterns of hotspot localization and evolution. Science, this issue p. 913, p. 928; see also p. 932 DNA recombination hotspots in birds are localized to gene regulatory sequences and show evolutionary conservation. [Also see Perspective by Lichten] The DNA-binding protein PRDM9 has a critical role in specifying meiotic recombination hotspots in mice and apes, but it appears to be absent from other vertebrate species, including birds. To study the evolution and determinants of recombination in species lacking the gene that encodes PRDM9, we inferred fine-scale genetic maps from population resequencing data for two bird species: the zebra finch, Taeniopygia guttata, and the long-tailed finch, Poephila acuticauda. We found that both species have recombination hotspots, which are enriched near functional genomic elements. Unlike in mice and apes, most hotspots are shared between the two species, and their conservation seems to extend over tens of millions of years. These observations suggest that in the absence of PRDM9, recombination targets functional features that both enable access to the genome and constrain its evolution.

Clade-Limited Colonization in Brood Parasitic Finches (Vidua spp.)

The African brood parasitic finches (Vidua spp.) are host specialists that mimic the songs and nestling mouth markings of their finch hosts (family Estrildidae). Although recent molecular analyses suggest rapid speciation associated with host switches in some members of this group, the association of different Vidua lineages with particular host genera suggests the possibility of cospeciation at higher levels in the host and parasite phylogenies. We compared a phylogeny of all Vidua species with a phylogeny of their estrildid finch hosts and compared divergence time estimates for the two groups. Basal divergences among extant members of the Vidulidae and among Vidua species are more recent than those among host genera and species, respectively, allowing a model of cospeciation to be rejected at most or all levels of the Vidua phylogeny. Nonetheless, some tests for cospeciation indicated significant congruence between host and parasite tree topologies. This result may be an artifact of clade-limited colonization. Host switches in parasitic finches have most often involved new hosts in the same or a closely related genus, an effect that increases the apparent congruence of host and parasites trees.

Comparative genomics based on massive parallel transcriptome sequencing reveals patterns of substitution and selection across 10 bird species

Next‐generation sequencing technology provides an attractive means to obtain large‐scale sequence data necessary for comparative genomic analysis. To analyse the patterns of mutation rate variation and selection intensity across the avian genome, we performed brain transcriptome sequencing using Roche 454 technology of 10 different non‐model avian species. Contigs from de novo assemblies were aligned to the two available avian reference genomes, chicken and zebra finch. In total, we identified 6499 different genes across all 10 species, with ∼1000 genes found in each full run per species. We found evidence for a higher mutation rate of the Z chromosome than of autosomes (male‐biased mutation) and a negative correlation between the neutral substitution rate (dS) and chromosome size. Analyses of the mean dN/dS ratio (ω) of genes across chromosomes supported the Hill–Robertson effect (the effect of selection at linked loci) and point at stochastic problems with ω as an independent measure of selection. Overall, this study demonstrates the usefulness of next‐generation sequencing for obtaining genomic resources for comparative genomic analysis of non‐model organisms.

Nucleotide Variation, Linkage Disequilibrium and Founder-facilitated Speciation in Wild Populations of the Zebra Finch (Taeniopygia guttata)

The zebra finch has long been an important model system for the study of vocal learning, vocal production, and behavior. With the imminent sequencing of its genome, the zebra finch is now poised to become a model system for population genetics. Using a panel of 30 noncoding loci, we characterized patterns of polymorphism and divergence among wild zebra finch populations. Continental Australian populations displayed little population structure, exceptionally high levels of nucleotide diversity (π = 0.010), a rapid decay of linkage disequilibrium (LD), and a high population recombination rate (ρ ≈ 0.05), all of which suggest an open and fluid genomic background that could facilitate adaptive variation. By contrast, substantial divergence between the Australian and Lesser Sunda Island populations (KST = 0.193), reduced genetic diversity (π = 0.002), and higher levels of LD in the island population suggest a strong but relatively recent founder event, which may have contributed to speciation between these populations as envisioned under founder-effect speciation models. Consistent with this hypothesis, we find that under a simple quantitative genetic model both drift and selection could have contributed to the observed divergence in six quantitative traits. In both Australian and Lesser Sundas populations, diversity in Z-linked loci was significantly lower than in autosomal loci. Our analysis provides a quantitative framework for studying the role of selection and drift in shaping patterns of molecular evolution in the zebra finch genome.

Molecular evolution of genes in avian genomes

BackgroundObtaining 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.ResultsWe 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.ConclusionsFocusing 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.

Phylogeography and conservation genetics of Eld's deer (Cervus eldi)

Elds deer (Cervus eldi) is a highly endangered cervid, distributed historically throughout much of South Asia and Indochina. We analysed variation in the mitochondrial DNA (mtDNA) control region for representatives of all three Elds deer subspecies to gain a better understanding of the genetic population structure and evolutionary history of this species. A phylogeny of mtDNA haplotypes indicates that the critically endangered and ecologically divergent C. eldi eldi is related more closely to C. e. thamin than to C. e. siamensis, a result that is consistent with biogeographic considerations. The results also suggest a strong degree of phylogeographic structure both between subspecies and among populations within subspecies, suggesting that dispersal of individuals between populations has been very limited historically. Haplotype diversity was relatively high for two of the three subspecies (thamin and siamensis), indicating that recent population declines have not yet substantially eroded genetic diversity. In contrast, we found no haplotype variation within C. eldi eldi or the Hainan Island population of C. eldi siamensis, two populations which are known to have suffered severe population bottlenecks. We also compared levels of haplotype and nucleotide diversity in an unmanaged captive population, a managed captive population and a relatively healthy wild population. Diversity indices were higher in the latter two, suggesting the efficacy of well‐designed breeding programmes for maintaining genetic diversity in captivity. Based on significant genetic differentiation among Elds deer subspecies, we recommend the continued management of this species in three distinct evolutionarily significant units (ESUs). Where possible, it may be advisable to translocate individuals between isolated populations within a subspecies to maintain levels of genetic variation in remaining Elds deer populations.

Integrating Genomes, Brain and Behavior in the Study of Songbirds

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

Digital gene expression analysis of the zebra finch genome

BackgroundIn order to understand patterns of adaptation and molecular evolution it is important to quantify both variation in gene expression and nucleotide sequence divergence. Gene expression profiling in non-model organisms has recently been facilitated by the advent of massively parallel sequencing technology. Here we investigate tissue specific gene expression patterns in the zebra finch (Taeniopygia guttata) with special emphasis on the genes of the major histocompatibility complex (MHC).ResultsAlmost 2 million 454-sequencing reads from cDNA of six different tissues were assembled and analysed. A total of 11,793 zebra finch transcripts were represented in this EST data, indicating a transcriptome coverage of about 65%. There was a positive correlation between the tissue specificity of gene expression and non-synonymous to synonymous nucleotide substitution ratio of genes, suggesting that genes with a specialised function are evolving at a higher rate (or with less constraint) than genes with a more general function. In line with this, there was also a negative correlation between overall expression levels and expression specificity of contigs. We found evidence for expression of 10 different genes related to the MHC. MHC genes showed relatively tissue specific expression levels and were in general primarily expressed in spleen. Several MHC genes, including MHC class I also showed expression in brain. Furthermore, for all genes with highest levels of expression in spleen there was an overrepresentation of several gene ontology terms related to immune function.ConclusionsOur study highlights the usefulness of next-generation sequence data for quantifying gene expression in the genome as a whole as well as in specific candidate genes. Overall, the data show predicted patterns of gene expression profiles and molecular evolution in the zebra finch genome. Expression of MHC genes in particular, corresponds well with expression patterns in other vertebrates.