The genome of a songbird

The constitutive differential transcriptome of a brain circuit for vocal learning

Background

The ability to imitate the vocalizations of other organisms, a trait known as vocal learning, is shared by only a few organisms, including humans, where it subserves the acquisition of speech and language, and 3 groups of birds. In songbirds, vocal learning requires the coordinated activity of a set of specialized brain nuclei referred to as the song control system. Recent efforts have revealed some of the genes that are expressed in these vocal nuclei, however a thorough characterization of the transcriptional specializations of this system is still missing. We conducted a rigorous and comprehensive analysis of microarrays, and conducted a separate analysis of 380 genes by in situ hybridizations in order to identify molecular specializations of the major nuclei of the song system of zebra finches (Taeniopygia guttata), a songbird species.

Results

Our efforts identified more than 3300 genes that are differentially regulated in one or more vocal nuclei of adult male birds compared to the adjacent brain regions. Bioinformatics analyses provided insights into the possible involvement of these genes in molecular pathways such as cellular morphogenesis, intrinsic cellular excitability, neurotransmission and neuromodulation, axonal guidance and cela-to-cell interactions, and cell survival, which are known to strongly influence the functional properties of the song system. Moreover, an in-depth analysis of specific gene families with known involvement in regulating the development and physiological properties of neuronal circuits provides further insights into possible modulators of the song system.

Conclusion

Our study represents one of the most comprehensive molecular characterizations of a brain circuit that evolved to facilitate a learned behavior in a vertebrate. The data provide novel insights into possible molecular determinants of the functional properties of the song control circuitry. It also provides lists of compelling targets for pharmacological and genetic manipulations to elucidate the molecular regulation of song behavior and vocal learning.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4578-0) contains supplementary material, which is available to authorized users.

Karyotype Evolution in Birds: From Conventional Staining to Chromosome Painting

In the last few decades, there have been great efforts to reconstruct the phylogeny of Neoaves based mainly on DNA sequencing. Despite the importance of karyotype data in phylogenetic studies, especially with the advent of fluorescence in situ hybridization (FISH) techniques using different types of probes, the use of chromosomal data to clarify phylogenetic proposals is still minimal. Additionally, comparative chromosome painting in birds is restricted to a few orders, while in mammals, for example, virtually all orders have already been analyzed using this method. Most reports are based on comparisons using Gallus gallus probes, and only a small number of species have been analyzed with more informative sets of probes, such as those from Leucopternis albicollis and Gyps fulvus, which show ancestral macrochromosomes rearranged in alternative patterns. Despite this, it is appropriate to review the available cytogenetic information and possible phylogenetic conclusions. In this report, the authors gather both classical and molecular cytogenetic data and describe some interesting and unique characteristics of karyotype evolution in birds.

Comparative analysis examining patterns of genomic differentiation across multiple episodes of population divergence in birds

Heterogeneous patterns of genomic differentiation are commonly documented between closely related populations and there is considerable interest in identifying factors that contribute to their formation. These factors could include genomic features (e.g., areas of low recombination) that promote processes like linked selection (positive or purifying selection that affects linked neutral sites) at specific genomic regions. Examinations of repeatable patterns of differentiation across population pairs can provide insight into the role of these factors. Birds are well suited for this work, as genome structure is conserved across this group. Accordingly, we reestimated relative (F_ST ) and absolute (d_XY ) differentiation between eight sister pairs of birds that span a broad taxonomic range using a common pipeline. Across pairs, there were modest but significant correlations in window-based estimates of differentiation (up to 3% of variation explained for F_ST and 26% for d_XY ), supporting a role for processes at conserved genomic features in generating heterogeneous patterns of differentiation; processes specific to each episode of population divergence likely explain the remaining variation. The role genomic features play was reinforced by linear models identifying several genomic variables (e.g., gene densities) as significant predictors of F_ST and d_XY repeatability. F_ST repeatability was higher among pairs that were further along the speciation continuum (i.e., more reproductively isolated) providing further insight into how genomic differentiation changes with population divergence; early stages of speciation may be dominated by positive selection that is different between pairs but becomes integrated with processes acting according to shared genomic features as speciation proceeds.

CNVs are associated with genomic architecture in a songbird

Background

Understanding variation in genome structure is essential to understand phenotypic differences within populations and the evolutionary history of species. A promising form of this structural variation is copy number variation (CNV). CNVs can be generated by different recombination mechanisms, such as non-allelic homologous recombination, that rely on specific characteristics of the genome architecture. These structural variants can therefore be more abundant at particular genes ultimately leading to variation in phenotypes under selection. Detailed characterization of CNVs therefore can reveal evolutionary footprints of selection and provide insight in their contribution to phenotypic variation in wild populations.

Results

Here we use genotypic data from a long-term population of great tits (Parus major), a widely studied passerine bird in ecology and evolution, to detect CNVs and identify genomic features prevailing within these regions. We used allele intensities and frequencies from high-density SNP array data from 2,175 birds. We detected 41,029 CNVs concatenated into 8,008 distinct CNV regions (CNVRs). We successfully validated 93.75% of the CNVs tested by qPCR, which were sampled at different frequencies and sizes. A mother-daughter family structure allowed for the evaluation of the inheritance of a number of these CNVs. Thereby, only CNVs with 40 probes or more display segregation in accordance with Mendelian inheritance, suggesting a high rate of false negative calls for smaller CNVs. As CNVRs are a coarse-grained map of CNV loci, we also inferred the frequency of coincident CNV start and end breakpoints. We observed frequency-dependent enrichment of these breakpoints at homologous regions, CpG sites and AT-rich intervals. A gene ontology enrichment analyses showed that CNVs are enriched in genes underpinning neural, cardiac and ion transport pathways.

Conclusion

Great tit CNVs are present in almost half of the genes and prominent at repetitive-homologous and regulatory regions. Although overlapping genes under selection, the high number of false negatives make neutrality or association tests on CNVs detected here difficult. Therefore, CNVs should be further addressed in the light of their false negative rate and architecture to improve the comprehension of their association with phenotypes and evolutionary history.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4577-1) contains supplementary material, which is available to authorized users.

Multidirectional chromosome painting in Synallaxis frontalis (Passeriformes, Furnariidae) reveals high chromosomal reorganization, involving fissions and inversions

In this work we performed comparative chromosome painting using probes from Gallus gallus (GGA) Linnaeus, 1758 and Leucopternis albicollis (LAL) Latham, 1790 in Synallaxis frontalis Pelzeln, 1859 (Passeriformes, Furnariidae), an exclusively Neotropical species, in order to analyze whether the complex pattern of intrachromosomal rearrangements (paracentric and pericentric inversions) proposed for Oscines and Suboscines is shared with more basal species. S. frontalis has 82 chromosomes, similar to most Avian species, with a large number of microchromosomes and a few pairs of macrochromosomes. We found polymorphisms in pairs 1 and 3, where homologues were submetacentric and acrocentric. Hybridization of GGA probes showed syntenies in the majority of ancestral macrochromosomes, except for GGA1 and GGA2, which hybridized to more than one pair of chromosomes each. LAL probes confirmed the occurrence of intrachromosomal rearrangements in the chromosomes corresponding to GGA1q, as previously proposed for species from the order Passeriformes. In addition, LAL probes suggest that pericentric inversions or centromere repositioning were responsible for variations in the morphology of the heteromorphic pairs 1 and 3. Altogether, the analysis of our data on chromosome painting and the data published in other Passeriformes highlights chromosomal changes that have occurred during the evolution of Passeriformes.

Draft genome assembly of the Bengalese finch, Lonchura striata domestica, a model for motor skill variability and learning

Background

Vocal learning in songbirds has emerged as a powerful model for sensorimotor learning. Neurobehavioral studies of Bengalese finch (Lonchura striata domestica) song, naturally more variable and plastic than songs of other finch species, have demonstrated the importance of behavioral variability for initial learning, maintenance, and plasticity of vocalizations. However, the molecular and genetic underpinnings of this variability and the learning it supports are poorly understood.

Findings

To establish a platform for the molecular analysis of behavioral variability and plasticity, we generated an initial draft assembly of the Bengalese finch genome from a single male animal to 151× coverage and an N50 of 3.0 MB. Furthermore, we developed an initial set of gene models using RNA-seq data from 8 samples that comprise liver, muscle, cerebellum, brainstem/midbrain, and forebrain tissue from juvenile and adult Bengalese finches of both sexes.

Conclusions

We provide a draft Bengalese finch genome and gene annotation to facilitate the study of the molecular-genetic influences on behavioral variability and the process of vocal learning. These data will directly support many avenues for the identification of genes involved in learning, including differential expression analysis, comparative genomic analysis (through comparison to existing avian genome assemblies), and derivation of genetic maps for linkage analysis. Bengalese finch gene models and sequences will be essential for subsequent manipulation (molecular or genetic) of genes and gene products, enabling novel mechanistic investigations into the role of variability in learned behavior.

Rapid regulatory evolution of a nonrecombining autosome linked to divergent behavioral phenotypes

In the white-throated sparrow (Zonotrichia albicollis), the second chromosome bears a striking resemblance to sex chromosomes. First, within each breeding pair of birds, one bird is homozygous for the standard arrangement of the chromosome (ZAL2/ZAL2) and its mate is heterozygous for a different version (ZAL2/ZAL2^m). Second, recombination is profoundly suppressed between the two versions, leading to genetic differentiation between them. Third, the ZAL2^m version is linked with phenotypic traits, such as bright plumage, high aggression, and low parental behavior, which are usually associated with males. These similarities to sex chromosomes suggest that the evolutionary mechanisms that shape sex chromosomes, in particular genetic degeneration of the heterogametic version due to the suppression of recombination, are likely important in this system as well. Here, we investigated patterns of protein sequence evolution and gene expression evolution between the ZAL2 and ZAL2^m chromosomes by whole-genome sequencing and transcriptome analyses. Patterns of protein evolution exhibited only weak signals of genetic degeneration, and few genes harbored signatures of positive selection. We found substantial evidence of transcriptome evolution, such as significant expression divergence between ZAL2 and ZAL2^m alleles and signatures of dosage compensation for highly expressed genes. These results suggest that, early in the evolution of heteromorphic chromosomes, gene expression divergence and dosage compensation can prevail before large-scale genetic degeneration. Our results show further that suppression of recombination between heteromorphic chromosomes can lead to the evolution of alternative (sex-like) behavioral phenotypes before substantial genetic degeneration.

Repetitive DNAs and shrink genomes: A chromosomal analysis in nine Columbidae species (Aves, Columbiformes)

An extensive karyotype variation is found among species belonging to the Columbidae family of birds (Columbiformes), both in diploid number and chromosomal morphology. Although clusters of repetitive DNA sequences play an important role in chromosomal instability, and therefore in chromosomal rearrangements, little is known about their distribution and amount in avian genomes. The aim of this study was to analyze the distribution of 11 distinct microsatellite sequences, as well as clusters of 18S rDNA, in nine different Columbidae species, correlating their distribution with the occurrence of chromosomal rearrangements. We found 2n values ranging from 76 to 86 and nine out of 11 microsatellite sequences showed distinct hybridization signals among the analyzed species. The accumulation of microsatellite repeats was found preferentially in the centromeric region of macro and microchromosomes, and in the W chromosome. Additionally, pair 2 showed the accumulation of several microsatellites in different combinations and locations in the distinct species, suggesting the occurrence of intrachromosomal rearrangements, as well as a possible fission of this pair in Geotrygon species. Therefore, although birds have a smaller amount of repetitive sequences when compared to other Tetrapoda, these seem to play an important role in the karyotype evolution of these species.

Genomic associations with bill length and disease reveal drift and selection across island bird populations

Island species provide excellent models for investigating how selection and drift operate in wild populations, and for determining how these processes act to influence local adaptation and speciation. Here, we examine the role of selection and drift in shaping genomic and phenotypic variation across recently separated populations of Berthelot's pipit (Anthus berthelotii), a passerine bird endemic to three archipelagos in the Atlantic. We first characterized genetic diversity and population structuring that supported previous inferences of a history of recent colonizations and bottlenecks. We then tested for regions of the genome associated with the ecologically important traits of bill length and malaria infection, both of which vary substantially across populations in this species. We identified a SNP associated with variation in bill length among individuals, islands, and archipelagos; patterns of variation at this SNP suggest that both phenotypic and genotypic variation in bill length is largely shaped by founder effects. Malaria was associated with SNPs near/within genes involved in the immune response, but this relationship was not consistent among archipelagos, supporting the view that disease resistance is complex and rapidly evolving. Although we found little evidence for divergent selection at candidate loci for bill length and malaria resistance, genome scan analyses pointed to several genes related to immunity and metabolism as having important roles in divergence and adaptation. Our findings highlight the utility and challenges involved with combining association mapping and population genetic analysis in nonequilibrium populations, to disentangle the effects of drift and selection on shaping genotypes and phenotypes.

Primordial germ cell-mediated transgenesis and genome editing in birds

Transgenesis and genome editing in birds are based on a unique germline transmission system using primordial germ cells (PGCs), which is quite different from the mammalian transgenic and genome editing system. PGCs are progenitor cells of gametes that can deliver genetic information to the next generation. Since avian PGCs were first discovered in nineteenth century, there have been numerous efforts to reveal their origin, specification, and unique migration pattern, and to improve germline transmission efficiency. Recent advances in the isolation and in vitro culture of avian PGCs with genetic manipulation and genome editing tools enable the development of valuable avian models that were unavailable before. However, many challenges remain in the production of transgenic and genome-edited birds, including the precise control of germline transmission, introduction of exogenous genes, and genome editing in PGCs. Therefore, establishing reliable germline-competent PGCs and applying precise genome editing systems are critical current issues in the production of avian models. Here, we introduce a historical overview of avian PGCs and their application, including improved techniques and methodologies in the production of transgenic and genome-edited birds, and we discuss the future potential applications of transgenic and genome-edited birds to provide opportunities and benefits for humans.

A genetic signature of the evolution of loss of flight in the Galapagos cormorant

We have a limited understanding of the genetic and molecular basis of evolutionary changes in the size and proportion of limbs. We studied wing and pectoral skeleton reduction leading to flightlessness in the Galapagos cormorant (Phalacrocorax harrisi). We sequenced and de novo assembled the genomes of four cormorant species and applied a predictive and comparative genomics approach to find candidate variants that may have contributed to the evolution of flightlessness. These analyses and cross-species experiments in Caenorhabditis elegans and in chondrogenic cell lines implicated variants in genes necessary for transcriptional regulation and function of the primary cilium. Cilia are essential for Hedgehog signaling, and humans affected by skeletal ciliopathies suffer from premature bone growth arrest, mirroring skeletal features associated with loss of flight.

A dual transcript-discovery approach to improve the delimitation of gene features from RNA-seq data in the chicken model

The sequence of the chicken genome, like several other draft genome sequences, is presently not fully covered. Gaps, contigs assigned with low confidence and uncharacterized chromosomes result in gene fragmentation and imprecise gene annotation. Transcript abundance estimation from RNA sequencing (RNA-seq) data relies on read quality, library complexity and expression normalization. In addition, the quality of the genome sequence used to map sequencing reads, and the gene annotation that defines gene features, must also be taken into account. A partially covered genome sequence causes the loss of sequencing reads from the mapping step, while an inaccurate definition of gene features induces imprecise read counts from the assignment step. Both steps can significantly bias interpretation of RNA-seq data. Here, we describe a dual transcript-discovery approach combining a genome-guided gene prediction and a de novo transcriptome assembly. This dual approach enabled us to increase the assignment rate of RNA-seq data by nearly 20% as compared to when using only the chicken reference annotation, contributing therefore to a more accurate estimation of transcript abundance. More generally, this strategy could be applied to any organism with partial genome sequence and/or lacking a manually-curated reference annotation in order to improve the accuracy of gene expression studies.

Avian genomics lends insights into endocrine function in birds

The genomics era has brought along the completed sequencing of a large number of bird genomes that cover a broad range of the avian phylogenetic tree (>30 orders), leading to major novel insights into avian biology and evolution. Among recent findings, the discovery that birds lack a large number of protein coding genes that are organized in highly conserved syntenic clusters in other vertebrates is very intriguing, given the physiological importance of many of these genes. A considerable number of them play prominent endocrine roles, suggesting that birds evolved compensatory genetic or physiological mechanisms that allowed them to survive and thrive in spite of these losses. While further studies are needed to establish the exact extent of avian gene losses, these findings point to birds as potentially highly relevant model organisms for exploring the genetic basis and possible therapeutic approaches for a wide range of endocrine functions and disorders.

Synchrotron microtomography of a Nothosaurus marchicus skull informs on nothosaurian physiology and neurosensory adaptations in early Sauropterygia

Nothosaurs form a subclade of the secondarily marine Sauropterygia that was well represented in late Early to early Late Triassic marine ecosystems. Here we present and discuss the internal skull anatomy of the small piscivorous nothosaur Nothosaurus marchicus from coastal to shallow marine Lower Muschelkalk deposits (Anisian) of Winterswijk, The Netherlands, which represents the oldest sauropterygian endocast visualized to date. The cranial endocast is only partially encapsulated by ossified braincase elements. Cranial flattening and lateral constriction by hypertrophied temporal musculature grant the brain a straight, tubular geometry that lacks particularly well-developed cerebral lobes but does potentially involve distinguishable optic lobes, suggesting vision may have represented an important sense during life. Despite large orbit size, the circuitous muscular pathway linking the basisphenoidal and orbital regions indicates poor oculomotor performance. This suggests a rather fixed ocular orientation, although eye placement and neck manoeuvrability could have enabled binocular if not stereoscopic vision. The proportionally large dorsal projection of the braincase endocast towards the well-developed pineal foramen advocates substantial dependence on the corresponding pineal system in vivo. Structures corroborating keen olfactory or acoustic senses were not identified. The likely atrophied vomeronasal organ argues against the presence of a forked tongue in Nothosaurus, and the relative positioning of external and internal nares contrasts respiratory configurations proposed for pistosauroid sauropterygians. The antorbital domain furthermore accommodates a putative rostral sensory plexus and pronounced lateral nasal glands that were likely exapted as salt glands. Previously proposed nothosaurian 'foramina eustachii' arose from architectural constraints on braincase development rather than representing functional foramina. Several modifications to brain shape and accessory organs were achieved through heterochronic development of the cranium, particularly the braincase. In summary, the cranium of Nothosaurus marchicus reflects important physiological and neurosensory adaptations that enabled the group's explosive invasion of shallow marine habitats in the late Early Triassic.

Androgen and estrogen sensitivity of bird song: a comparative view on gene regulatory levels

Singing of songbirds is sensitive to testosterone and its androgenic and estrogenic metabolites in a species-specific way. The hormonal effects on song pattern are likely mediated by androgen receptors (AR) and estrogen receptor alpha (ERα), ligand activated transcription factors that are expressed in neurons of various areas of the songbirds' vocal control circuit. The distribution of AR in this circuit is rather similar between species while that of ERα is species variant and concerns a key vocal control area, the HVC (proper name). We discuss the regulation of the expression of the cognate AR and ERα and putative splice variants. In particular, we suggest that transcription factor binding sites in the promoter of these receptors differ between bird species. Further, we suggest that AR- and ERα-dependent gene regulation in vocal areas differs between species due to species-specific DNA binding sites of putative target genes that are required for the transcriptional activity of the receptors. We suggest that species differences in the distribution of AR and ERα in vocal areas and in the genomic sensitivity to these receptors contribute to species-specific hormonal regulation of the song.

Urban landscape genomics identifies fine-scale gene flow patterns in an avian invasive

Invasive species exert a serious impact on native fauna and flora and have been the target of many eradication and management efforts worldwide. However, a lack of data on population structure and history, exacerbated by the recency of many species introductions, limits the efficiency with which such species can be kept at bay. In this study we generated a novel genome of high assembly quality and genotyped 4735 genome-wide single nucleotide polymorphic (SNP) markers from 78 individuals of an invasive population of the Javan Myna Acridotheres javanicus across the island of Singapore. We inferred limited population subdivision at a micro-geographic level, a genetic patch size (~13-14 km) indicative of a pronounced dispersal ability, and barely an increase in effective population size since introduction despite an increase of four to five orders of magnitude in actual population size, suggesting that low population-genetic diversity following a bottleneck has not impeded establishment success. Landscape genomic analyses identified urban features, such as low-rise neighborhoods, that constitute pronounced barriers to gene flow. Based on our data, we consider an approach targeting the complete eradication of Javan Mynas across Singapore to be unfeasible. Instead, a mixed approach of localized mitigation measures taking into account urban geographic features and planning policy may be the most promising avenue to reducing the adverse impacts of this urban pest. Our study demonstrates how genomic methods can directly inform the management and control of invasive species, even in geographically limited datasets with high gene flow rates.

Reciprocal interactions between prostaglandin E2- and estradiol-dependent signaling pathways in the injured zebra finch brain

Background

Astrocytic aromatization and consequent increases in estradiol are neuroprotective in the injured brain. In zebra finches, cyclooxygenase-activity is necessary for injury-induced aromatase expression, and increased central estradiol lowers neuroinflammation. The mechanisms underlying these influences are unknown. Here, we document injury-induced, cyclooxygenase-dependent increases in glial aromatase expression and replicate previous work in our lab showing increases in central prostaglandin E2 and estradiol following brain damage. Further, we describe injury-dependent changes in E-prostanoid and estrogen receptor expression and reveal the necessity of E-prostanoid and estrogen receptors in the injury-dependent, reciprocal interactions of neuroinflammatory and neurosteroidogenic pathways.

Methods

Adult male and female birds were shams or received bilateral injections of the appropriate drug or vehicle into contralateral telencephalic lobes.

Results

Injuries sustained in the presence of indomethacin (a cyclooxygenase inhibitor) had fewer aromatase-expressing reactive astrocytes relative to injuries injected with vehicle suggesting that cyclooxygenase activity is necessary for the induction of glial aromatase around the site of damage. Injured hemispheres had higher prostaglandin E2 and estradiol content relative to shams. Importantly, injured hemispheres injected with E-prostanoid- or estrogen receptor-antagonists showed elevated prostaglandin E2 and estradiol, respectively, but lower prostaglandin E2 or estradiol-dependent downstream activity (protein kinase A or phosphoinositide-3-kinase mRNA) suggesting that receptor antagonism did not affect injury-induced prostaglandin E2 or estradiol, but inhibited the effects of these ligands. Antagonism of E-prostanoid receptors 3 or 4 prevented injury-induced increases in neural estradiol in males and females, respectively, albeit this apparent sex-difference needs to be tested more stringently. Further, estrogen receptor-α, but not estrogen receptor-β antagonism, exaggerated neural prostaglandin E2 levels relative to the contralateral lobe in both sexes.

Conclusion

These data suggest injury-induced, sex-specific prostaglandin E2-dependent estradiol synthesis, and estrogen receptor-α dependent decreases in neuroinflammation in the vertebrate brain.

De-novo emergence of SINE retroposons during the early evolution of passerine birds

BACKGROUND

Passeriformes ("perching birds" or passerines) make up more than half of all extant bird species. The genome of the zebra finch, a passerine model organism for vocal learning, was noted previously to contain thousands of short interspersed elements (SINEs), a group of retroposons that is abundant in mammalian genomes but considered largely inactive in avian genomes.

RESULTS

Here we resolve the deep phylogenetic relationships of passerines using presence/absence patterns of SINEs. The resultant retroposon-based phylogeny provides a powerful and independent corroboration of previous sequence-based analyses. Notably, SINE activity began in the common ancestor of Eupasseres (passerines excluding the New Zealand wrens Acanthisittidae) and ceased before the rapid diversification of oscine passerines (suborder Passeri - songbirds). Furthermore, we find evidence for very recent SINE activity within suboscine passerines (suborder Tyranni), following the emergence of a SINE via acquisition of a different tRNA head as we suggest through template switching.

CONCLUSIONS

We propose that the early evolution of passerines was unusual among birds in that it was accompanied by de-novo emergence and activity of SINEs. Their genomic and transcriptomic impact warrants further study in the light of the massive diversification of passerines.

Characterization of a de novo assembled transcriptome of the Common Blackbird (Turdus merula)

BACKGROUND

In recent years, next generation high throughput sequencing technologies have proven to be useful tools for investigations concerning the genomics or transcriptomics also of non-model species. Consequently, ornithologists have adopted these technologies and the respective bioinformatics tools to survey the genomes and transcriptomes of a few avian non-model species. The Common Blackbird is one of the most common bird species living in European cities, which has successfully colonized urban areas and for which no reference genome or transcriptome is publicly available. However, to target questions like genome wide gene expression analysis, a reference genome or transcriptome is needed.

METHODS

Therefore, in this study two Common Blackbirds were sacrificed, their mRNA was isolated and analyzed by RNA-Seq to de novo assemble a transcriptome and characterize it. Illumina reads (125 bp paired-end) and a Velvet/Oases pipeline led to 162,158 transcripts. For the annotation (using Blast+), an unfiltered protein database was used. SNPs were identified using SAMtools and BCFtools. Furthermore, mRNA from three single tissues (brain, heart and liver) of the same two Common Blackbirds were sequenced by Illumina (75 bp single-end reads). The draft transcriptome and the three single tissues were compared by their BLAST hits with the package VennDiagram in R.

RESULTS

Following the annotation against protein databases, we found evidence for 15,580 genes in the transcriptome (all well characterized hits after annotation). On 18% of the assembled transcripts, 144,742 SNPs were identified which are, consequently, 0.09% of all nucleotides in the assembled transcriptome. In the transcriptome and in the single tissues (brain, heart and liver), 10,182 shared genes were found.

DISCUSSION

Using a next-generation technology and bioinformatics tools, we made a first step towards the genomic investigation of the Common Blackbird. The de novo assembled transcriptome is usable for downstream analyses such as differential gene expression analysis and SNP identification. This study shows the importance of the approach to sequence single tissues to understand functions of tissues, proteins and the phenotype.

Mutations in different pigmentation genes are associated with parallel melanism in island flycatchers

The independent evolution of similar traits across multiple taxa provides some of the most compelling evidence of natural selection. Little is known, however, about the genetic basis of these convergent or parallel traits: are they mediated by identical or different mutations in the same genes, or unique mutations in different genes? Using a combination of candidate gene and reduced representation genomic sequencing approaches, we explore the genetic basis of and the evolutionary processes that mediate similar plumage colour shared by isolated populations of the Monarcha castaneiventris flycatcher of the Solomon Islands. A genome-wide association study (GWAS) that explicitly controlled for population structure revealed that mutations in known pigmentation genes are the best predictors of parallel plumage colour. That is, entirely black or melanic birds from one small island share an amino acid substitution in the melanocortin-1 receptor (MC1R), whereas similarly melanic birds from another small island over 100 km away share an amino acid substitution in a predicted binding site of agouti signalling protein (ASIP). A third larger island, which separates the two melanic populations, is inhabited by birds with chestnut bellies that lack the melanic MC1R and ASIP allelic variants. Formal FST outlier tests corroborated the results of the GWAS and suggested that strong, directional selection drives the near fixation of the MC1R and ASIP variants across islands. Our results, therefore, suggest that selection acting on different mutations with large phenotypic effects can drive the evolution of parallel melanism, despite the relatively small population size on islands.

Fundamental constraints in synchronous muscle limit superfast motor control in vertebrates

Superfast muscles (SFMs) are extremely fast synchronous muscles capable of contraction rates up to 250 Hz, enabling precise motor execution at the millisecond time scale. SFM phenotypes have been discovered in most major vertebrate lineages, but it remains unknown whether all SFMs share excitation-contraction coupling pathway adaptations for speed, and if SFMs arose once, or from independent evolutionary events. Here, we demonstrate that to achieve rapid actomyosin crossbridge kinetics bat and songbird SFM express myosin heavy chain genes that are evolutionarily and ontologically distinct. Furthermore, we show that all known SFMs share multiple functional adaptations that minimize excitation-contraction coupling transduction times. Our results suggest that SFM evolved independently in sound-producing organs in ray-finned fish, birds, and mammals, and that SFM phenotypes operate at a maximum operational speed set by fundamental constraints in synchronous muscle. Consequentially, these constraints set a fundamental limit to the maximum speed of fine motor control.

Across animals, different muscle types have evolved to perform vastly different tasks at different speeds. For example, tortoise leg muscles move slowly over several seconds, while the flight muscles of a hummingbird move quickly dozens of times per second. The speed record holders among vertebrates are the so-called superfast muscles, which can move up to 250 times per second. Superfast muscles power the alarming rattle of rattlesnakes, courtship calls in fish, rapid echolocation calls in bats and the elaborate vocal gymnastics of songbirds. Thus these extreme muscles are all around us and are always involved in sound production.

Did superfast muscles evolve from a common ancestor? And how do different superfast muscles achieve their extreme behavior? To answer these questions, Mead et al. studied the systems known to limit contraction speed in all currently known superfast muscles found in rattlesnakes, toadfish, bats and songbirds. This revealed that all the muscles share certain specific adaptations that allow superfast contractions. Furthermore, the three fastest examples – toadfish, songbird and bat – have nearly identical maximum speeds. Although this appears to support the idea that the adaptations all evolved from a shared ancestor, Mead et al. found evidence that suggests otherwise. Each of the three superfast muscles are powered by a different motor protein, which argues strongly in favor of the muscles evolving independently. The existence of such similar mechanisms and performance in independently evolved muscles raises the possibility that the fastest contraction rates measured by Mead et al. represent a maximum speed limit for all vertebrate muscles.

Any technical failure in a racecar most likely will slow it down, while the same failure in a robustly engineered family car may not be so noticeable. Similarly in superfast muscle many cellular and molecular systems need to perform maximally. Therefore by understanding how these extreme muscles work, we also gain a better understanding of how normal muscles contract.

Developmental song learning as a model to understand neural mechanisms that limit and promote the ability to learn

Songbirds famously learn their vocalizations. Some species can learn continuously, others seasonally, and still others just once. The zebra finch (Taeniopygia guttata) learns to sing during a single developmental "Critical Period," a restricted phase during which a specific experience has profound and permanent effects on brain function and behavioral patterns. The zebra finch can therefore provide fundamental insight into features that promote and limit the ability to acquire complex learned behaviors. For example, what properties permit the brain to come "on-line" for learning? How does experience become encoded to prevent future learning? What features define the brain in receptive compared to closed learning states? This piece will focus on epigenomic, genomic, and molecular levels of analysis that operate on the timescales of development and complex behavioral learning. Existing data will be discussed as they relate to Critical Period learning, and strategies for future studies to more directly address these questions will be considered. Birdsong learning is a powerful model for advancing knowledge of the biological intersections of maturation and experience. Lessons from its study not only have implications for understanding developmental song learning, but also broader questions of learning potential and the enduring effects of early life experience on neural systems and behavior.

Mosaic genome evolution in a recent and rapid avian radiation

Recent genomic analyses of evolutionary radiations suggest that ancestral or standing genetic variation may facilitate rapid diversification, particularly in cases involving convergence in ecological traits. Likewise, lateral transfer of alleles via hybridization may also facilitate adaptive convergence, but little is known about the role of ancestral variation in examples of explosive diversification that primarily involve the evolution of species recognition traits. Here, we show that genomic regions distinguishing sympatric species in an extraordinary radiation of small finches called munias (genus Lonchura) have phylogenetic histories that are discordant with each other, with the overall pattern of autosomal differentiation among species, and with sex-linked and mitochondrial components of the genome. Genome-wide data for 11 species sampled in Australia and Papua New Guinea indicate substantial autosomal introgression between sympatric species, but also identify a limited number of divergent autosomal regions, several of which overlap known colour genes (ASIP, EDN3, IGSF11, KITLG, MC1R and SOX10). Phylogenetic analysis of these outlier regions shows that different munia species have acquired unique combinations of alleles across a relatively small set of phenotypically relevant genes. Our results demonstrate that the recombination of ancestral genetic variation across multiple loci may be an important mechanism for generating phenotypic novelty and diversity.

Detailed temporal structure of communication networks in groups of songbirds

Animals in groups often exchange calls, in patterns whose temporal structure may be influenced by contextual factors such as physical location and the social network structure of the group. We introduce a model-based analysis for temporal patterns of animal call timing, originally developed for networks of firing neurons. This has advantages over cross-correlation analysis in that it can correctly handle common-cause confounds and provides a generative model of call patterns with explicit parameters for the influences between individuals. It also has advantages over standard Markovian analysis in that it incorporates detailed temporal interactions which affect timing as well as sequencing of calls. Further, a fitted model can be used to generate novel synthetic call sequences. We apply the method to calls recorded from groups of domesticated zebra finch (Taeniopygia guttata) individuals. We find that the communication network in these groups has stable structure that persists from one day to the next, and that ‘kernels’ reflecting the temporal range of influence have a characteristic structure for a calling individual's effect on itself, its partner and on others in the group. We further find characteristic patterns of influences by call type as well as by individual.

Northern Spotted Owl (Strix occidentalis caurina) Genome: Divergence with the Barred Owl (Strix varia) and Characterization of Light-Associated Genes

We report here the assembly of a northern spotted owl (Strix occidentalis caurina) genome. We generated Illumina paired-end sequence data at 90× coverage using nine libraries with insert lengths ranging from ∼250 to 9,600 nt and read lengths from 100 to 375 nt. The genome assembly is comprised of 8,108 scaffolds totaling 1.26 × 109 nt in length with an N50 length of 3.98 × 106 nt. We calculated the genome-wide fixation index (FST) of S. o. caurina with the closely related barred owl (Strix varia) as 0.819. We examined 19 genes that encode proteins with light-dependent functions in our genome assembly as well as in that of the barn owl (Tyto alba). We present genomic evidence for loss of three of these in S. o. caurina and four in T. alba. We suggest that most light-associated gene functions have been maintained in owls and their loss has not proceeded to the same extent as in other dim-light-adapted vertebrates.

On the independent gene trees assumption in phylogenomic studies

Multilocus coalescent methods for inferring species trees or historical demographic parameters typically require the assumption that gene trees for sampled SNPs or DNA sequence loci are conditionally independent given their species tree. In practice, researchers have used different criteria to delimit "independent loci." One criterion identifies sampled loci as being independent of each other if they undergo Mendelian independent assortment (IA criterion). O'Neill et al. (2013, Molecular Ecology, 22, 111-129) used this approach in their phylogeographic study of North American tiger salamander species complex. In two other studies, researchers developed a pair of related methods that employ an independent genealogies criterion (IG criterion), which considers the effects of population-level recombination on correlations between the gene trees of intrachromosomal loci. Here, I explain these three methods, illustrate their use with example data, and evaluate their efficacies. I show that the IA approach is more conservative, is simpler to use and requires fewer assumptions than the IG approaches. However, IG approaches can identify much larger numbers of independent loci than the IA method, which, in turn, allows researchers to obtain more precise and accurate estimates of species trees and historical demographic parameters. A disadvantage of the IG methods is that they require an estimate of the population recombination rate. Despite their drawbacks, IA and IG approaches provide molecular ecologists with promising a priori methods for selecting SNPs or DNA sequence loci that likely meet the independence assumption in coalescent-based phylogenomic studies.

Habitat Association Predicts Genetic Diversity and Population Divergence in Amazonian Birds

The ecological traits of organisms may predict their genetic diversity and population genetic structure and mediate the action of evolutionary processes important for speciation and adaptation. Making these ecological-evolutionary links is difficult because it requires comparable genetic estimates from many species with differing ecologies. In Amazonian birds, habitat association is an important component of ecological diversity. Here, we examine the link between habitat association and genetic parameters using 20 pairs of closely related Amazonian bird species in which one member of the pair occurs primarily in forest edge and floodplains and the other occurs in upland forest interior. We use standardized geographic sampling and data from 2,416 genomic markers to estimate genetic diversity, population genetic structure, and statistics reflecting demographic and evolutionary processes. We find that species of upland forest have greater genetic diversity and divergence across the landscape as well as signatures of older histories and less gene flow than floodplain species. Our results reveal that species ecology in the form of habitat association is an important predictor of genetic diversity and population divergence and suggest that differences in diversity between floodplain and upland avifaunas in the Amazon may be driven by differences in the demographic and evolutionary processes at work in the two habitats.

Annotated Draft Genome Assemblies for the Northern Bobwhite (Colinus virginianus) and the Scaled Quail (Callipepla squamata) Reveal Disparate Estimates of Modern Genome Diversity and Historic Effective Population Size

Northern bobwhite (Colinus virginianus; hereafter bobwhite) and scaled quail (Callipepla squamata) populations have suffered precipitous declines across most of their US ranges. Illumina-based first- (v1.0) and second- (v2.0) generation draft genome assemblies for the scaled quail and the bobwhite produced N50 scaffold sizes of 1.035 and 2.042 Mb, thereby producing a 45-fold improvement in contiguity over the existing bobwhite assembly, and ≥90% of the assembled genomes were captured within 1313 and 8990 scaffolds, respectively. The scaled quail assembly (v1.0 = 1.045 Gb) was ∼20% smaller than the bobwhite (v2.0 = 1.254 Gb), which was supported by kmer-based estimates of genome size. Nevertheless, estimates of GC content (41.72%; 42.66%), genome-wide repetitive content (10.40%; 10.43%), and MAKER-predicted protein coding genes (17,131; 17,165) were similar for the scaled quail (v1.0) and bobwhite (v2.0) assemblies, respectively. BUSCO analyses utilizing 3023 single-copy orthologs revealed a high level of assembly completeness for the scaled quail (v1.0; 84.8%) and the bobwhite (v2.0; 82.5%), as verified by comparison with well-established avian genomes. We also detected 273 putative segmental duplications in the scaled quail genome (v1.0), and 711 in the bobwhite genome (v2.0), including some that were shared among both species. Autosomal variant prediction revealed ∼2.48 and 4.17 heterozygous variants per kilobase within the scaled quail (v1.0) and bobwhite (v2.0) genomes, respectively, and estimates of historic effective population size were uniformly higher for the bobwhite across all time points in a coalescent model. However, large-scale declines were predicted for both species beginning ∼15-20 KYA.

An intra-population analysis of the indris' song dissimilarity in the light of genetic distance

The increasing interest in the evolution of human language has led several fields of research to focus on primate vocal communication. The 'singing primates', which produce elaborated and complex sequences of vocalizations, are of particular interest for this topic. Indris (Indri indri) are the only singing lemurs and emit songs whose most distinctive portions are "descending phrases" consisting of 2-5 units. We examined how the structure of the indris' phrases varied with genetic relatedness among individuals. We tested whether the acoustic structure could provide conspecifics with information about individual identity and group membership. When analyzing phrase dissimilarity and genetic distance of both sexes, we found significant results for males but not for females. We found that similarity of male song-phrases correlates with kin in both time and frequency parameters, while, for females, this information is encoded only in the frequency of a single type. Song phrases have consistent individual-specific features, but we did not find any potential for advertising group membership. We emphasize the fact that genetic and social factors may play a role in the acoustic plasticity of female indris. Altogether, these findings open a new perspective for future research on the possibility of vocal production learning in these primates.

Chromosomal inversion differences correlate with range overlap in passerine birds

Chromosomal inversions evolve frequently but the reasons for this remain unclear. We used cytological descriptions of 411 species of passerine birds to identify large pericentric inversion differences between species, based on the position of the centromere. Within 81 small clades comprising 284 of the species, we found 319 differences on the 9 largest autosomes combined, 56 on the Z chromosome, and 55 on the W chromosome. We also identified inversions present within 32 species. Using a new fossil-calibrated phylogeny, we examined the phylogenetic, demographic and genomic context in which these inversions have evolved. The number of inversion differences between closely related species is consistently predicted by whether the ranges of species overlap, even when time is controlled for as far as is possible. Fixation rates vary across the autosomes, but inversions are more likely to be fixed on the Z chromosome than the average autosome. Variable mutagenic input alone (estimated by chromosome size, map length, GC content or repeat density) cannot explain the differences between chromosomes in the number of inversions fixed. Together, these results support a model in which inversions increase because of their effects on recombination suppression in the face of hybridization. Other factors associated with hybridization may also contribute, including the possibility that inversions contain incompatibility alleles, making taxa less likely to collapse following secondary contact.

Lentiviral-Mediated Transgenesis in Songbirds

Transgenesis involves the insertion of an exogenous gene into an animal's genome, which allows the identification of the expressed phenotypes in brain function or behavior. Lentiviral-mediated transgenesis offers unique transduction potency making it possible to deliver and stably integrate transgenes into a wide variety of dividing and nondividing cells. The ability to establish long-term expression of such transgenes allows their use for transgenesis which is especially useful in organisms lacking quality pluripotent stem cell lines and which is otherwise difficult to produce via traditional pronuclear microinjection, such as songbirds. Here we describe a protocol to generate the transgenic songbird, the zebra finch, by producing and inserting lentiviral-mediated transgene into the blastoderm of freshly laid eggs. This protocol includes procedures for production of lentiviral vectors, injection of a virus into zebra finch embryos, and postinjection care. The implementation of the songbird transgenic approach provides a leap toward basic and translational neuroscience that uses an animal model for speech and language and their pathologies. Additionally, the highly quantifiable song behavior, combined with a well-characterized song circuitry, offers an exciting opportunity to develop therapeutic strategies for neurological disorders.

A sex-chromosome inversion causes strong overdominance for sperm traits that affect siring success

Male reproductive success depends on the competitive ability of sperm to fertilize the ova, which should lead to strong selection on sperm characteristics. This raises the question of how heritable variation in sperm traits is maintained. Here we show that in zebra finches (Taeniopygia guttata) nearly half of the variance in sperm morphology is explained by an inversion on the Z chromosome with a 40% allele frequency in the wild. The sperm of males that are heterozygous for the inversion had the longest midpieces and the highest velocity. Furthermore, such males achieved the highest fertility and the highest siring success, both within-pair and extra-pair. Males homozygous for the derived allele show detrimental sperm characteristics and the lowest siring success. Our results suggest heterozygote advantage as the mechanism that maintains the inversion polymorphism and hence variance in sperm design and in fitness.

Low diversity, activity, and density of transposable elements in five avian genomes

In this study, we conducted the activity, diversity, and density analysis of transposable elements (TEs) across five avian genomes (budgerigar, chicken, turkey, medium ground finch, and zebra finch) to explore the potential reason of small genome sizes of birds. We found that these avian genomes exhibited low density of TEs by about 10% of genome coverages and low diversity of TEs with the TE landscapes dominated by CR1 and ERV elements, and contrasting proliferation dynamics both between TE types and between species were observed across the five avian genomes. Phylogenetic analysis revealed that CR1 clade was more diverse in the family structure compared with R2 clade in birds; avian ERVs were classified into four clades (alpha, beta, gamma, and ERV-L) and belonged to three classes of ERV with an uneven distributed in these lineages. The activities of DNA and SINE TEs were very low in the evolution history of avian genomes; most LINEs and LTRs were ancient copies with a substantial decrease of activity in recent, with only LTRs and LINEs in chicken and zebra finch exhibiting weak activity in very recent, and very few TEs were intact; however, the recent activity may be underestimated due to the sequencing/assembly technologies in some species. Overall, this study demonstrates low diversity, activity, and density of TEs in the five avian species; highlights the differences of TEs in these lineages; and suggests that the current and recent activity of TEs in avian genomes is very limited, which may be one of the reasons of small genome sizes in birds.

Form of Dietary Methylmercury does not Affect Total Mercury Accumulation in the Tissues of Zebra Finch

Exposure to mercury in humans, other mammals, and birds is primarily dietary, with mercury in the methylated form and bound to cysteine in the tissues of prey items. Yet dosing studies are generally carried out using methylmercury chloride. Here we tested whether the accumulation of total mercury in zebra finch blood, egg, muscle, liver, kidney or brain differed depending on whether dietary mercury was complexed with chloride or cysteine. We found no effect of form of mercury on tissue accumulation. Some previous studies have found lower accumulation of mercury in tissues of animals fed complexed mercury. Much remains to be understood about what happens to ingested mercury once it enters the intestines, but our results suggest that dietary studies using methylmercury chloride in birds will produce similar tissue accumulation levels to those using methylmercury cysteine.

Transcriptional response to West Nile virus infection in the zebra finch (Taeniopygia guttata)

West Nile virus (WNV) is a widespread arbovirus that imposes a significant cost to both human and wildlife health. WNV exists in a bird-mosquito transmission cycle in which passerine birds act as the primary reservoir host. As a public health concern, the mammalian immune response to WNV has been studied in detail. Little, however, is known about the avian immune response to WNV. Avian taxa show variable susceptibility to WNV and what drives this variation is unknown. Thus, to study the immune response to WNV in birds, we experimentally infected captive zebra finches (Taeniopygia guttata). Zebra finches provide a useful model, as like many natural avian hosts they are moderately susceptible to WNV and thus provide sufficient viremia to infect mosquitoes. We performed RNAseq in spleen tissue during peak viremia to provide an overview of the transcriptional response. In general, we find strong parallels with the mammalian immune response to WNV, including upregulation of five genes in the Rig-I-like receptor signalling pathway, and offer insights into avian-specific responses. Together with complementary immunological assays, we provide a model of the avian immune response to WNV and set the stage for future comparative studies among variably susceptible populations and species.

Combination of short-read, long-read, and optical mapping assemblies reveals large-scale tandem repeat arrays with population genetic implications

Accurate and contiguous genome assembly is key to a comprehensive understanding of the processes shaping genomic diversity and evolution. Yet, it is frequently constrained by constitutive heterochromatin, usually characterized by highly repetitive DNA. As a key feature of genome architecture associated with centromeric and subtelomeric regions, it locally influences meiotic recombination. In this study, we assess the impact of large tandem repeat arrays on the recombination rate landscape in an avian speciation model, the Eurasian crow. We assembled two high-quality genome references using single-molecule real-time sequencing (long-read assembly [LR]) and single-molecule optical maps (optical map assembly [OM]). A three-way comparison including the published short-read assembly (SR) constructed for the same individual allowed assessing assembly properties and pinpointing misassemblies. By combining information from all three assemblies, we characterized 36 previously unidentified large repetitive regions in the proximity of sequence assembly breakpoints, the majority of which contained complex arrays of a 14-kb satellite repeat or its 1.2-kb subunit. Using whole-genome population resequencing data, we estimated the population-scaled recombination rate (ρ) and found it to be significantly reduced in these regions. These findings are consistent with an effect of low recombination in regions adjacent to centromeric or subtelomeric heterochromatin and add to our understanding of the processes generating widespread heterogeneity in genetic diversity and differentiation along the genome. By combining three different technologies, our results highlight the importance of adding a layer of information on genome structure that is inaccessible to each approach independently.

Blood transcriptomes and de novo identification of candidate loci for mating success in lekking great snipe (Gallinago media)

We assembled the great snipe blood transcriptome using data from fourteen lekking males, in order to de novo identify candidate genes related to sexual selection, and determined the expression profiles in relation to mating success. The three most highly transcribed genes were encoding different haemoglobin subunits. All tended to be overexpressed in males with high mating success. We also called single nucleotide polymorphisms (SNPs) from the transcriptome data and found considerable genetic variation for many genes expressed during lekking. Among these, we identified 14 polymorphic candidate SNPs that had a significant genotypic association with mating success (number of females mated with) and/or mating status (mated or not). Four of the candidate SNPs were found in HBAA (encoding the haemoglobin α-chain). Heterozygotes for one of these and one SNP in the gene PABPC1 appeared to enjoy higher mating success compared to males homozygous for either of the alleles. In a larger data set of individuals, we genotyped 38 of the identified SNPs but found low support for consistent selection as only one of the zygosities of previously identified candidate SNPs and none of their genotypes were associated with mating status. However, candidate SNPs generally showed lower levels of spatial genetic structure compared to noncandidate markers. We also scored the prevalence of avian malaria in a subsample of birds. Males infected with avian malaria parasites had lower mating success in the year of sampling than noninfected males. Parasite infection and its interaction with specific genes may thus affect performance on the lek.

Overexpression of human NR2B receptor subunit in LMAN causes stuttering and song sequence changes in adult zebra finches

Zebra finches (Taeniopygia guttata) learn to produce songs in a manner reminiscent of spoken language development in humans. One candidate gene implicated in influencing learning is the N-methyl-D-aspartate (NMDA) subtype 2B glutamate receptor (NR2B). Consistent with this idea, NR2B levels are high in the song learning nucleus LMAN (lateral magnocellular nucleus of the anterior nidopallium) during juvenile vocal learning, and decreases to low levels in adults after learning is complete and the song becomes more stereotyped. To test for the role of NR2B in generating song plasticity, we manipulated NR2B expression in LMAN of adult male zebra finches by increasing its protein levels to those found in juvenile birds, using a lentivirus containing the full-length coding sequence of the human NR2B subunit. We found that increased NR2B expression in adult LMAN induced increases in song sequence diversity and slower song tempo more similar to juvenile songs, but also increased syllable repetitions similar to stuttering. We did not observe these effects in control birds with overexpression of NR2B outside of LMAN or with the green fluorescent protein (GFP) in LMAN. Our results suggest that low NR2B subunit expression in adult LMAN is important in conserving features of stereotyped adult courtship song.

Mapping brain structure and function: cellular resolution, global perspective

A comprehensive understanding of the brain requires analysis, although from a global perspective, with cellular, and even subcellular, resolution. An important step towards this goal involves the establishment of three-dimensional high-resolution brain maps, incorporating brain-wide information about the cells and their connections, as well as the chemical architecture. The progress made in such anatomical brain mapping in recent years has been paralleled by the development of physiological techniques that enable investigators to generate global neural activity maps, also with cellular resolution, while simultaneously recording the organism's behavioral activity. Combination of the high-resolution anatomical and physiological maps, followed by theoretical systems analysis of the deduced network, will offer unprecedented opportunities for a better understanding of how the brain, as a whole, processes sensory information and generates behavior.