Towards complete and error-free genome assemblies of all vertebrate species.

High-quality and complete reference genome assemblies are fundamental for the application of genomics to biology, disease, and biodiversity conservation. However, such assemblies are available for only a few non-microbial species1-4. To address this issue, the international Genome 10K (G10K) consortium5,6 has worked over a five-year period to evaluate and develop cost-effective methods for assembling highly accurate and nearly complete reference genomes. Here we present lessons learned from generating assemblies for 16 species that represent six major vertebrate lineages. We confirm that long-read sequencing technologies are essential for maximizing genome quality, and that unresolved complex repeats and haplotype heterozygosity are major sources of assembly error when not handled correctly. Our assemblies correct substantial errors, add missing sequence in some of the best historical reference genomes, and reveal biological discoveries. These include the identification of many false gene duplications, increases in gene sizes, chromosome rearrangements that are specific to lineages, a repeated independent chromosome breakpoint in bat genomes, and a canonical GC-rich pattern in protein-coding genes and their regulatory regions. Adopting these lessons, we have embarked on the Vertebrate Genomes Project (VGP), an international effort to generate high-quality, complete reference genomes for all of the roughly 70,000 extant vertebrate species and to help to enable a new era of discovery across the life sciences.

Layered evolution of gene expression in "superfast" muscles for courtship.

Identifying the molecular process of complex trait evolution is a core goal of biology. However, pinpointing the specific context and timing of trait-associated changes within the molecular evolutionary history of an organism remains an elusive goal. We study this topic by exploring the molecular basis of elaborate courtship evolution, which represents an extraordinary example of trait innovation. Within the behaviorally diverse radiation of Central and South American manakin birds, species from two separate lineages beat their wings together using specialized "superfast" muscles to generate a "snap" that helps attract mates. Here, we develop an empirical approach to analyze phylogenetic lineage-specific shifts in gene expression in the key snap-performing muscle and then integrate these findings with comparative transcriptomic sequence analysis. We find that rapid wing displays are associated with changes to a wide range of molecular processes that underlie extreme muscle performance, including changes to calcium trafficking, myocyte homeostasis and metabolism, and hormone action. We furthermore show that these changes occur gradually in a layered manner across the species history, wherein which ancestral genetic changes to many of these molecular systems are built upon by later species-specific shifts that likely finalized the process of display performance adaptation. Our study demonstrates the potential for combining phylogenetic modeling of tissue-specific gene expression shifts with phylogenetic analysis of lineage-specific sequence changes to reveal holistic evolutionary histories of complex traits.

Neurogenomic landscape associated with status-dependent cooperative behaviour.

The neurogenomic mechanisms mediating male-male reproductive cooperative behaviours remain unknown. We leveraged extensive transcriptomic and behavioural data on a neotropical bird species (Pipra filicauda) that performs cooperative courtship displays to understand these mechanisms. In this species, the cooperative display is modulated by testosterone, which promotes cooperation in non-territorial birds, but suppresses cooperation in territory holders. We sought to understand the neurogenomic underpinnings of three related traits: social status, cooperative display behaviour and testosterone phenotype. To do this, we profiled gene expression in 10 brain nuclei spanning the social decision-making network (SDMN), and two key endocrine tissues that regulate social behaviour. We associated gene expression with each bird's behavioural and endocrine profile derived from 3 years of repeated measures taken from free-living birds in the Ecuadorian Amazon. We found distinct landscapes of constitutive gene expression were associated with social status, testosterone phenotype and cooperation, reflecting the modular organization and engagement of neuroendocrine tissues. Sex-steroid and neuropeptide signalling appeared to be important in mediating status-specific relationships between testosterone and cooperation, suggesting shared regulatory mechanisms with male aggressive and sexual behaviours. We also identified differentially regulated genes involved in cellular activity and synaptic potentiation, suggesting multiple mechanisms underpin these genomic states. Finally, we identified SDMN-wide gene expression differences between territorial and floater males that could form the basis of 'status-specific' neurophysiological phenotypes, potentially mediated by testosterone and growth hormone. Overall, our findings provide new, systems-level insights into the mechanisms of cooperative behaviour and suggest that differences in neurogenomic state are the basis for individual differences in social behaviour.

Systems biology as a framework to understand the physiological and endocrine bases of behavior and its evolution-From concepts to a case study in birds.

Organismal behavior, with its tremendous complexity and diversity, is generated by numerous physiological systems acting in coordination. Understanding how these systems evolve to support differences in behavior within and among species is a longstanding goal in biology that has captured the imagination of researchers who work on a multitude of taxa, including humans. Of particular importance are the physiological determinants of behavioral evolution, which are sometimes overlooked because we lack a robust conceptual framework to study mechanisms underlying adaptation and diversification of behavior. Here, we discuss a framework for such an analysis that applies a "systems view" to our understanding of behavioral control. This approach involves linking separate models that consider behavior and physiology as their own networks into a singular vertically integrated behavioral control system. In doing so, hormones commonly stand out as the links, or edges, among nodes within this system. To ground our discussion, we focus on studies of manakins (Pipridae), a family of Neotropical birds. These species have numerous physiological and endocrine specializations that support their elaborate reproductive displays. As a result, manakins provide a useful example to help imagine and visualize the way systems concepts can inform our appreciation of behavioral evolution. In particular, manakins help clarify how connectedness among physiological systems-which is maintained through endocrine signaling-potentiate and/or constrain the evolution of complex behavior to yield behavioral differences across taxa. Ultimately, we hope this review will continue to stimulate thought, discussion, and the emergence of research focused on integrated phenotypes in behavioral ecology and endocrinology.

Molecular mechanisms of local adaptation for salt-tolerance in a treefrog.

Salinization is a global phenomenon affecting ecosystems and forcing freshwater organisms to deal with increasing levels of ionic stress. However, our understanding of mechanisms that permit salt tolerance in amphibians is limited. This study investigates mechanisms of salt tolerance in locally adapted, coastal populations of a treefrog, Hyla cinerea. Using a common garden experiment, we (i) determine the extent that environment (i.e., embryonic and larval saltwater exposure) or genotype (i.e., coastal vs. inland) affects developmental benchmarks and transcriptome expression, and (ii) identify genes that may underpin differences in saltwater tolerance. Differences in gene expression, survival, and plasma osmolality were most strongly associated with genotype. Population genetic analyses on expressed genes also delineated coastal and inland groups based on genetic similarity. Coastal populations differentially expressed osmoregulatory genes including ion transporters (atp1b1, atp6V1g2, slc26a), cellular adhesion components (cdh26, cldn1, gjb3, ocln), and cytoskeletal components (odc1-a, tgm3). Several of these genes are the same genes expressed by euryhaline fish after exposure to freshwater, which is a novel finding for North American amphibians and suggests that these genes may be associated with local salinity adaptation. Coastal populations also highly expressed glycerol-3-phosphate dehydrogenase 1 (gpd1), which indicates they use glycerol as a compatible osmolyte to reduce water loss - another mechanism of saltwater tolerance previously unknown in frogs. These data signify that Hyla cinerea inhabiting coastal, brackish wetlands have evolved a salt-tolerant ecotype, and highlights novel candidate pathways that can lead to salt tolerance in freshwater organisms facing habitat salinization.

Light pollution increases West Nile virus competence of a ubiquitous passerine reservoir species.

Among the many anthropogenic changes that impact humans and wildlife, one of the most pervasive but least understood is light pollution. Although detrimental physiological and behavioural effects resulting from exposure to light at night are widely appreciated, the impacts of light pollution on infectious disease risk have not been studied. Here, we demonstrate that artificial light at night (ALAN) extends the infectious-to-vector period of the house sparrow (Passer domesticus), an urban-dwelling avian reservoir host of West Nile virus (WNV). Sparrows exposed to ALAN maintained transmissible viral titres for 2 days longer than controls but did not experience greater WNV-induced mortality during this window. Transcriptionally, ALAN altered the expression of gene regulatory networks including key hubs (OASL, PLBD1 and TRAP1) and effector genes known to affect WNV dissemination (SOCS). Despite mounting anti-viral immune responses earlier, transcriptomic signatures indicated that ALAN-exposed individuals probably experienced pathogen-induced damage and immunopathology, potentially due to evasion of immune effectors. A simple mathematical modelling exercise indicated that ALAN-induced increases of host infectious-to-vector period could increase WNV outbreak potential by approximately 41%. ALAN probably affects other host and vector traits relevant to transmission, and additional research is needed to advise the management of zoonotic diseases in light-polluted areas.

Parent and offspring genotypes influence gene expression in early life.

Parents can have profound effects on offspring fitness. Little, however, is known about the mechanisms through which parental genetic variation influences offspring physiology in natural systems. White-throated sparrows (Zonotrichia albicollis, WTSP) exist in two genetic morphs, tan and white, controlled by a large polymorphic supergene. Morphs mate disassortatively, resulting in two pair types: tan male × white female (T × W) pairs, which provide biparental care and white male × tan female (W × T) pairs, which provide female-biased care. To investigate how parental composition impacts offspring, we performed RNA-seq on whole blood of WTSP nestlings sampled from nests of both pair types. Parental pair type had a large effect on nestling gene expression, with 881 genes differentially expressed (DE) and seven correlated gene coexpression modules. The DE genes and modules expressed at higher levels in W × T nests with female-biased parental care function in metabolism and stress-related pathways resulting from the overrepresentation of proteolysis and stress-response genes (e.g., SOD2, NR3C1). These results show that parental genotypes and/or associated behaviours influence nestling physiology, and highlight avenues of further research investigating the ultimate implications for the maintenance of this polymorphism. Nestlings also exhibited morph-specific gene expression, with 92 differentially expressed genes, comprising immunity genes and genes encompassed by the supergene. Remarkably, we identified the same regulatory hub genes in these blood-derived expression networks as were previously identified in adult WTSP brains (EPM2A, BPNT1, TAF5L). These hub genes were located within the supergene, highlighting the importance of this gene complex in structuring regulatory networks across diverse tissues.

Mechanisms of Assortative Mating in Speciation with Gene Flow: Connecting Theory and Empirical Research.

The large body of theory on speciation with gene flow has brought to light fundamental differences in the effects of two types of mating rules on speciation: preference/trait rules, in which divergence in both (female) preferences and (male) mating traits is necessary for assortment, and matching rules, in which individuals mate with like individuals on the basis of the presence of traits or alleles that they have in common. These rules can emerge from a variety of behavioral or other mechanisms in ways that are not always obvious. We discuss the theoretical properties of both types of rules and explain why speciation is generally thought to be more likely under matching rather than preference/trait rules. We furthermore discuss whether specific assortative mating mechanisms fall under a preference/trait or matching rule, present empirical evidence for these mechanisms, and propose empirical tests that could distinguish between them. The synthesis of the theoretical literature on these assortative mating rules with empirical studies of the mechanisms by which they act can provide important insights into the occurrence of speciation with gene flow. Finally, by providing a clear framework we hope to inspire greater alignment in the ways that both theoreticians and empiricists study mating rules and how these rules affect speciation through maintaining or eroding barriers to gene flow among closely related species or populations.

A systems approach to animal communication.

Why animal communication displays are so complex and how they have evolved are active foci of research with a long and rich history. Progress towards an evolutionary analysis of signal complexity, however, has been constrained by a lack of hypotheses to explain similarities and/or differences in signalling systems across taxa. To address this, we advocate incorporating a systems approach into studies of animal communication--an approach that includes comprehensive experimental designs and data collection in combination with the implementation of systems concepts and tools. A systems approach evaluates overall display architecture, including how components interact to alter function, and how function varies in different states of the system. We provide a brief overview of the current state of the field, including a focus on select studies that highlight the dynamic nature of animal signalling. We then introduce core concepts from systems biology (redundancy, degeneracy, pluripotentiality, and modularity) and discuss their relationships with system properties (e.g. robustness, flexibility, evolvability). We translate systems concepts into an animal communication framework and accentuate their utility through a case study. Finally, we demonstrate how consideration of the system-level organization of animal communication poses new practical research questions that will aid our understanding of how and why animal displays are so complex.

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.

High major histocompatibility complex class I polymorphism despite bottlenecks in wild and domesticated populations of the zebra finch (Taeniopygia guttata).

BACKGROUND: Two subspecies of zebra finch, Taeniopygia guttata castanotis and T. g. guttata are native to Australia and the Lesser Sunda Islands, respectively. The Australian subspecies has been domesticated and is now an important model system for research. Both the Lesser Sundan subspecies and domesticated Australian zebra finches have undergone population bottlenecks in their history, and previous analyses using neutral markers have reported reduced neutral genetic diversity in these populations. Here we characterize patterns of variation in the third exon of the highly variable major histocompatibility complex (MHC) class I α chain. As a benchmark for neutral divergence, we also report the first mitochondrial NADH dehydrogenase 2 (ND2) sequences in this important model system. RESULTS: Despite natural and human-mediated population bottlenecks, we find that high MHC class I polymorphism persists across all populations. As expected, we find higher levels of nucleotide diversity in the MHC locus relative to neutral loci, and strong evidence of positive selection acting on important residues forming the peptide-binding region (PBR). Clear population differentiation of MHC allele frequencies is also evident, and this may be due to adaptation to new habitats and associated pathogens and/or genetic drift. Whereas the MHC Class I locus shows broad haplotype sharing across populations, ND2 is the first locus surveyed to date to show reciprocal monophyly of the two subspecies. CONCLUSIONS: Despite genetic bottlenecks and genetic drift, all surveyed zebra finch populations have maintained high MHC Class I diversity. The diversity at the MHC Class I locus in the Lesser Sundan subspecies contrasts sharply with the lack of diversity in previously examined neutral loci, and may thus be a result of selection acting to maintain polymorphism. Given uncertainty in historical population demography, however, it is difficult to rule out neutral processes in maintaining the observed diversity. The surveyed populations also differ in MHC Class I allele frequencies, and future studies are needed to assess whether these changes result in functional immune differences.

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.

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.

Gene expression in male and female stickleback from populations with convergent and divergent throat coloration.

Understanding of genetic mechanisms underlying variation in sexual dichromatism remains limited, especially for carotenoid-based colors. We addressed this knowledge gap in a gene expression study with threespine stickleback. We compared male and female throat tissues across five populations, including two in which female red coloration has evolved convergently. We found that the expression of individual genes, gene ontologies, and coexpression networks associated with red female color within a population differed between California and British Columbia populations, suggesting differences in underlying mechanisms. Comparing females from each of these populations to females from populations dominated by dull females, we again found extensive expression differences. For each population, genes and networks associated with female red color showed the same patterns for males only inconsistently. The functional roles of genes showing correlated expression with female color are unclear within populations, whereas genes highlighted through inter-population comparisons include some previously suggested to function in carotenoid pathways. Among these, the most consistent patterns involved TTC39B (Tetratricopeptide Repeat Domain 39B), which is within a known red coloration QTL in stickleback and implicated in red coloration in other taxa.

Comparative transcriptome analysis of Indian domestic duck reveals candidate genes associated with egg production.

Egg production is an important economic trait and a key indicator of reproductive performance in ducks. Egg production is regulated by several factors including genes. However the genes involved in egg production in duck remain unclear. In this study, we compared the ovarian transcriptome of high egg laying (HEL) and low egg laying (LEL) ducks using RNA-Seq to identify the genes involved in egg production. The HEL ducks laid on average 433 eggs while the LEL ducks laid 221 eggs over 93 weeks. A total of 489 genes were found to be significantly differentially expressed out of which 310 and 179 genes were up and downregulated, respectively, in the HEL group. Thirty-eight differentially expressed genes (DEGs), including LHX9, GRIA1, DBH, SYCP2L, HSD17B2, PAR6, CAPRIN2, STC2, and RAB27B were found to be potentially related to egg production and folliculogenesis. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis suggested that DEGs were enriched for functions related to glutamate receptor activity, serine-type endopeptidase activity, immune function, progesterone mediated oocyte maturation and MAPK signaling. Protein-protein interaction network analysis (PPI) showed strong interaction between 32 DEGs in two distinct clusters. Together, these findings suggest a mix of genetic and immunological factors affect egg production, and highlights candidate genes and pathways, that provides an understanding of the molecular mechanisms regulating egg production in ducks and in birds more broadly.

Dynamic molecular evolution of a supergene with suppressed recombination in white-throated sparrows.

In white-throated sparrows, two alternative morphs differing in plumage and behavior segregate with a large chromosomal rearrangement. As with sex chromosomes such as the mammalian Y, the rearranged version of chromosome two (ZAL2m) is in a near-constant state of heterozygosity, offering opportunities to investigate both degenerative and selective processes during the early evolutionary stages of 'supergenes.' Here, we generated, synthesized, and analyzed extensive genome-scale data to better understand the forces shaping the evolution of the ZAL2 and ZAL2m chromosomes in this species. We found that features of ZAL2m are consistent with substantially reduced recombination and low levels of degeneration. We also found evidence that selective sweeps took place both on ZAL2m and its standard counterpart, ZAL2, after the rearrangement event. Signatures of positive selection were associated with allelic bias in gene expression, suggesting that antagonistic selection has operated on gene regulation. Finally, we discovered a region exhibiting long-range haplotypes inside the rearrangement on ZAL2m. These haplotypes appear to have been maintained by balancing selection, retaining genetic diversity within the supergene. Together, our analyses illuminate mechanisms contributing to the evolution of a young chromosomal polymorphism, revealing complex selective processes acting concurrently with genetic degeneration to drive the evolution of supergenes.

Gene duplication and fragmentation in the zebra finch major histocompatibility complex.

BACKGROUND: Due 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. RESULTS: The 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. CONCLUSION: The 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.

Highly Contiguous Genomes Improve the Understanding of Avian Olfactory Receptor Repertoires.

Third-generation (long-read-based) sequencing technologies are reshaping our understanding of genome structure and function. One of the most persistent challenges in genome biology has been confidently reconstructing radiations of complex gene families. Olfactory receptors (ORs) represent just such a gene family with upward of thousands of receptors in some mammalian taxa. Whereas in birds olfaction was historically an overlooked sensory modality, new studies have revealed an important role for smell. Chromosome-level assemblies for birds allow a new opportunity to characterize patterns of OR diversity among major bird lineages. Previous studies of short-read-based (second-generation) genome assemblies have associated OR gene family size with avian ecology, but such conclusions could be premature especially when new assembly methods reshape our understanding of avian OR evolution. Here we provide a fundamental characterization of OR repertoires in five recent genome assemblies, including the most recent assembly of golden-collared manakin (Manacus vitellinus). We find that short read-based assemblies systematically undercount the avian-specific gamma-c OR subfamily, a subfamily that comprises over 65% of avian OR diversity. Therefore, in contrast to previous studies, we find a high diversity of gamma-c ORs across the avian tree of life. Building on these findings, ongoing sequencing efforts and improved genome assemblies will clarify the relationship between OR diversity and avian ecology.

Digital gene expression analysis of the zebra finch genome.

BACKGROUND: In 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). RESULTS: Almost 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. CONCLUSIONS: Our 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.

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 (pi = 0.010), a rapid decay of linkage disequilibrium (LD), and a high population recombination rate (rho approximately 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 (K(ST) = 0.193), reduced genetic diversity (pi = 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.