Search

The zebra finch is one of the most commonly studied songbirds in biology, particularly in genomics, neuroscience and vocal communication. However, this species lacks a robust cell line for molecular biology research and reagent optimization. We generated a cell line, designated CFS414, from zebra finch embryonic fibroblasts using the SV40 large and small T antigens. This cell line demonstrates an improvement over previous songbird cell lines through continuous and density-independent growth, allowing for indefinite culture and monoclonal line derivation. Cytogenetic, genomic, and transcriptomic profiling established the provenance of this cell line and identified the expression of genes relevant to ongoing songbird research. Using this cell line, we disrupted endogenous gene sequences using S.aureus Cas9 and confirmed a stress-dependent localization response of a song system specialized gene, SAP30L. The utility of CFS414 cells enhances the comprehensive molecular potential of the zebra finch and validates cell immortalization strategies in a songbird species., (© 2022. The Author(s).)

Zebra finch is a representative animal model for studying the molecular basis of human disorders of vocal development and communication. Accordingly, various functional studies of zebra finch have knocked down or introduced foreign genes in vivo; however, their germline transmission efficiency is remarkably low. The primordial germ cell (PGC)-mediated method is preferred for avian transgenic studies; however, use of this method is restricted in zebra finch due to the lack of an efficient gene transfer method for the germline. To target primary germ cells that are difficult to transfect and manipulate, an adenovirus-mediated gene transfer system with high efficiency in a wide range of cell types may be useful. Here, we isolated and characterized two types of primary germline-competent stem cells, PGCs and spermatogonial stem cells (SSCs), from embryonic and adult reproductive tissues of zebra finch and demonstrated that genes were most efficiently transferred into these cells using an adenovirus-mediated system. This system was successfully used to generate gene-edited PGCs in vitro. These results are expected to improve transgenic zebra finch production., (© 2021. The Author(s).)

Major histocompatibility complex (MHC) genes code for key proteins of the adaptive immune system, which present antigens from intra-cellular (MHC class I) and extra-cellular (MHC class II) pathogens. Because of their unprecedented diversity, MHC genes have long been an object of scientific interest, but due to methodological difficulties in genotyping of duplicated loci, our knowledge on the evolution of the MHC across different vertebrate lineages is still limited. Here, we compared the evolution of MHC class I and class II genes in three sister clades of common passerine birds, finches (Fringillinae and Carduelinae) and buntings (Emberizidae) using a uniform methodological (genotyping and data processing) approach and uniform sample sizes. Our analyses revealed contrasting evolutionary trajectories of the two MHC classes. We found a stronger signature of pervasive positive selection and higher allele diversity (allele numbers) at the MHC class I than class II. In contrast, MHC class II genes showed greater allele divergence (in terms of nucleotide diversity) and a much stronger recombination (gene conversion) signal. Gene copy numbers at both MHC class I and class II evolved via fluctuating selection and drift (Brownian Motion evolution), but the evolutionary rate was higher at class I. Our study constitutes one of few existing examples, where evolution of MHC class I and class II genes was directly compared using a multi-species approach. We recommend that re-focusing MHC research from single-species and single-class approaches towards multi-species analyses of both MHC classes can substantially increase our understanding MHC evolution in a broad phylogenetic context.

The fovea is a pit formed in the center of the retina that enables high-acuity vision in certain vertebrate species. While formation of the fovea fascinates many researchers, the molecular mechanisms underlying foveal development are poorly understood. In the current study, we histologically investigated foveal development in zebra finch (Taeniopygia guttata) and found that foveal pit formation begins just before post-hatch day 14 (P14). We next performed RNA-seq analysis to compare gene expression profiles between the central (foveal and parafoveal) and peripheral retina in zebra finch at P14. We found that the Arhgef33 expression is enriched in the middle layer of the inner nuclear layer at the parafovea, suggesting that Arhgef33 is dominantly expressed in Müller glial cells in the developing parafovea. We then performed a pull-down assay using Rhotekin-RBD and observed GEF activity of Arhgef33 against RhoA. We found that overexpression of Arhgef33 in HEK293 cells induces cell contraction and that Arhgef33 expression inhibits neurite extension in Neuro 2A cells, which is partially recovered by a Rho-kinase (ROCK) inhibitor. Taken together, we used zebra finch as a model animal to investigate foveal development and identified Arhgef33 as a candidate protein possibly involved in foveal development through modulating RhoA activity.

Whole-genome sequencing projects are increasingly populating the tree of life and characterizing biodiversity 1-4 . Sparse taxon sampling has previously been proposed to confound phylogenetic inference 5 , and captures only a fraction of the genomic diversity. Here we report a substantial step towards the dense representation of avian phylogenetic and molecular diversity, by analysing 363 genomes from 92.4% of bird families-including 267 newly sequenced genomes produced for phase II of the Bird 10,000 Genomes (B10K) Project. We use this comparative genome dataset in combination with a pipeline that leverages a reference-free whole-genome alignment to identify orthologous regions in greater numbers than has previously been possible and to recognize genomic novelties in particular bird lineages. The densely sampled alignment provides a single-base-pair map of selection, has more than doubled the fraction of bases that are confidently predicted to be under conservation and reveals extensive patterns of weak selection in predominantly non-coding DNA. Our results demonstrate that increasing the diversity of genomes used in comparative studies can reveal more shared and lineage-specific variation, and improve the investigation of genomic characteristics. We anticipate that this genomic resource will offer new perspectives on evolutionary processes in cross-species comparative analyses and assist in efforts to conserve species.

The transcription factor FOXP2 is crucial for the formation and function of cortico-striatal circuits. FOXP2 mutations are associated with specific speech and language impairments. In songbirds, experimentally altered FoxP2 expression levels in the striatal song nucleus Area X impair vocal learning and song production. Overall FoxP2 protein levels in Area X are low in adult zebra finches and decrease further with singing. However, some Area X medium spiny neurons (MSNs) express FoxP2 at high levels (FoxP2 high MSNs) and singing does not change this. Because Area X receives many new neurons throughout adulthood, we hypothesized that the FoxP2 high MSNs are newly recruited neurons, not yet integrated into the local Area X circuitry and thus not active during singing. Contrary to our expectation, FoxP2 protein levels did not predict whether new MSNs were active during singing, assayed via immediate early gene expression. However, new FoxP2 high MSNs had more complex dendrites, higher spine density and more mushroom spines than new FoxP2 low MSNs. In addition, FoxP2 expression levels correlated positively with nucleus size of new MSNs. Together, our data suggest that dynamic FoxP2 levels in new MSNs shape their morphology during maturation and their incorporation into a neural circuit that enables the maintenance and social modulation of adult birdsong.

All songbirds studied to date have an additional Germline Restricted Chromosome (GRC), which is not present in somatic cells. GRCs show a wide variation in genetic content and little homology between species. To check how this divergence affected the meiotic behavior of the GRC, we examined synapsis, recombination and copy number variation for GRCs in the closely related sand and pale martins (Riparia riparia and R. diluta) in comparison with distantly related estrildid finches. Using immunolocalization of meiotic proteins and FISH with GRC-specific DNA probes, we found a striking similarity in the meiotic behavior of GRCs between martins and estrildid finches despite the millions of years of independent evolution. GRCs are usually present in two copies in female and in one copy in male pachytene cells. However, we detected polymorphism in female and mosaicism in male martins for the number of GRCs. In martin and zebra finch females, two GRCs synapse along their whole length, but recombine predominately at their ends. We suggest that the shared features of the meiotic behavior of GRCs have been supported by natural selection in favor of a preferential segregation of GRCs to the eggs.

Developmental conditions can impact the adult phenotype via epigenetic changes that modulate gene expression. In mammals, methylation of the glucocorticoid receptor gene Nr3c1 has been implicated as mediator of long-term effects of developmental conditions, but this evidence is limited to humans and rodents, and few studies have simultaneously tested for associations between DNA methylation, gene expression and phenotype. Adverse environmental conditions during early life (large natal brood size) or adulthood (high foraging costs) exert multiple long-term phenotypic effects in zebra finches, and we here test for effects of these manipulations on DNA methylation and expression of the Nr3c1 gene in blood. Having been reared in a large brood induced higher DNA methylation of the Nr3c1 regulatory region in adulthood, and this effect persisted over years. Nr3c1 expression was negatively correlated with methylation at 2 out of 8 CpG sites, and was lower in hard foraging conditions, despite foraging conditions having no effect on Nr3c1 methylation at our target region. Nr3c1 expression also correlated with glucocorticoid traits: higher expression level was associated with lower plasma baseline corticosterone concentrations and enhanced corticosterone reactivity. Our results suggest that methylation of the Nr3c1 regulatory region can contribute to the mechanisms underlying the emergence of long-term effects of developmental conditions in birds, but in our system current adversity dominated over early life experiences with respect to receptor expression.

Estradiol provision via neural aromatization decreases neuro-inflammation and -degeneration, but almost nothing is known about the interactions between the peripheral immune system and brain aromatase. Given the vulnerability of the CNS we reasoned that brain aromatization may protect circuits from the threats of peripheral infection; perhaps shielding cells that are less resilient from the degeneration associated with peripheral infection or trauma. Lipopolysaccharide (LPS) or vehicle was administered peripherally to adult zebra finches and sickness behavior was recorded 2 or 24 hours later. The central transcription of cytokines and aromatase was measured, as were telencephalic aromatase activity and immunoreactive aromatase (24 hour time point only). Two hours post LPS, sickness-like behaviors increased, the transcription of IL-1β was higher in both sexes, and TNFα was elevated in females. 24 hours post-LPS, the behavior of LPS birds was similar to controls, and cytokines had returned to baseline, but aromatase mRNA and activity were elevated in both sexes. Immunocytochemistry revealed greater numbers of aromatase-expressing neurons in LPS birds. These data suggest that the activation of the immune system via peripheral endotoxin increases neuronal aromatase; a mechanism that may rapidly generate a potent anti-neuroinflammatory steroid in response to peripheral activation of the immune system.

The proportion of an individual's genome that is identical by descent (GWIBD) can be estimated from pedigrees (inbreeding coefficient 'Pedigree F') or molecular markers ('Marker F'), but both estimators come with error. Assuming unrelated pedigree founders, Pedigree F is the expected proportion of GWIBD given a specific inbreeding constellation. Meiotic recombination introduces variation around that expectation (Mendelian noise) and related pedigree founders systematically bias Pedigree F downward. Marker F is an estimate of the actual proportion of GWIBD but it suffers from the sampling error of markers plus the error that occurs when a marker is homozygous without reflecting common ancestry (identical by state). We here show via simulation of a zebra finch and a human linkage map that three aspects of meiotic recombination (independent assortment of chromosomes, number of crossovers and their distribution along chromosomes) contribute to variation in GWIBD and thus the precision of Pedigree and Marker F. In zebra finches, where the genome contains large blocks that are rarely broken up by recombination, the Mendelian noise was large (nearly twofold larger s.d. values compared with humans) and Pedigree F thus less precise than in humans, where crossovers are distributed more uniformly along chromosomes. Effects of meiotic recombination on Marker F were reversed, such that the same number of molecular markers yielded more precise estimates of GWIBD in zebra finches than in humans. As a consequence, in species inheriting large blocks that rarely recombine, even small numbers of microsatellite markers will often be more informative about inbreeding and fitness than large pedigrees.

Small heat shock proteins (sHSPs) range in size from 12 to 42 kDa and contain an α-crystalline domain. They have been proposed to play roles in the first line of defence against various stresses in an ATP-independent manner. In birds, a newly oviposited blastoderm can survive several weeks in a dormant state in low-temperature storage suggesting that blastoderm cells are basically tolerant of environmental stress. However, sHSPs in the stress-tolerant blastoderm have yet to be investigated. Thus, we characterised the expression and function of sHSPs in the chicken blastoderm. We found that chicken HSP25 was expressed especially in the blastoderm and was highly upregulated during low-temperature storage. Multiple alignments, phylogenetic trees, and expression in the blastoderms of Japanese quail and zebra finch showed homologues of HSP25 were conserved in other avian species. After knockdown of chicken HSP25, the expression of pluripotency marker genes decreased significantly. Furthermore, loss of function studies demonstrated that chicken HSP25 is associated with anti-apoptotic, anti-oxidant, and pro-autophagic effects in chicken blastoderm cells. Collectively, these results suggest avian HSP25 could play an important role in association with the first line of cellular defences against environmental stress and the protection of future embryonic cells in the avian blastoderm.

Colour polymorphism is known to facilitate speciation but the genetic basis of animal pigmentation and how colour polymorphisms contribute to speciation is poorly understood. Restricted recombination may promote linkage disequilibrium between the colour locus and incompatibility genes. Genomic rearrangement and the position of relevant loci within a chromosome are important factors that influence the frequency of recombination. Therefore, it is important to know the position of the colour locus, gene order and recombination landscape of the chromosome to understand the mechanism that generates incompatibilities between morphs. Recent studies showed remarkable pre- and postzygotic incompatibilities between sympatric colour morphs of the Gouldian finch (Erythrura gouldiae), in which head feather colour is genetically determined by a single sex-linked locus, Red. We constructed a genetic map for the Z chromosome of the Gouldian finch (male-specific map distance=131 cM), using 618 captive-bred birds and 34 microsatellite markers, to investigate the extent of inter- and intraspecific genomic rearrangements and variation in recombination rate within the Z chromosome. We refined the location of the Red locus to a ~7.2-cM interval in a region with a moderate recombination rate but outside the least-recombining, putative centromeric region. There was no evidence of chromosome-wide genomic rearrangements between the chromosomes carrying the red or black alleles with the current marker resolution. This work will contribute to identifying the causal gene, which will in turn enable alternative explanations for the association between incompatibility and colouration, such as fine-scale linkage disequilibrium, genomic rearrangements and pleiotropy, to be tested.

Learning and memory formation are known to require dynamic CpG (de)methylation and gene expression changes. Here, we aimed at establishing a genome-wide DNA methylation map of the zebra finch genome, a model organism in neuroscience, as well as identifying putatively epigenetically regulated genes. RNA- and MethylCap-seq experiments were performed on two zebra finch cell lines in presence or absence of 5-aza-2'-deoxycytidine induced demethylation. First, the MethylCap-seq methodology was validated in zebra finch by comparison with RRBS-generated data. To assess the influence of (variable) methylation on gene expression, RNA-seq experiments were performed as well. Comparison of RNA-seq and MethylCap-seq results showed that at least 357 of the 3,457 AZA-upregulated genes are putatively regulated by methylation in the promoter region, for which a pathway analysis showed remarkable enrichment for neurological networks. A subset of genes was validated using Exon Arrays, quantitative RT-PCR and CpG pyrosequencing on bisulfite-treated samples. To our knowledge, this study provides the first genome-wide DNA methylation map of the zebra finch genome as well as a comprehensive set of genes of which transcription is under putative methylation control.

Most molecular measures of inbreeding do not measure inbreeding at the scale that is most relevant for understanding inbreeding depression-namely the proportion of the genome that is identical-by-descent (IBD). The inbreeding coefficient FPed obtained from pedigrees is a valuable estimator of IBD, but pedigrees are not always available, and cannot capture inbreeding loops that reach back in time further than the pedigree. We here propose a molecular approach to quantify the realized proportion of the genome that is IBD (propIBD), and we apply this method to a wild and a captive population of zebra finches (Taeniopygia guttata). In each of 948 wild and 1057 captive individuals we analyzed available single-nucleotide polymorphism (SNP) data (260 SNPs) spread over four different genomic regions in each population. This allowed us to determine whether any of these four regions was completely homozygous within an individual, which indicates IBD with high confidence. In the highly nomadic wild population, we did not find a single case of IBD, implying that inbreeding must be extremely rare (propIBD=0-0.00094, 95% CI). In the captive population, a five-generation pedigree strongly underestimated the average amount of realized inbreeding (FPed=0.013

Darwin's finches, inhabiting the Galápagos archipelago and Cocos Island, constitute an iconic model for studies of speciation and adaptive evolution. Here we report the results of whole-genome re-sequencing of 120 individuals representing all of the Darwin's finch species and two close relatives. Phylogenetic analysis reveals important discrepancies with the phenotype-based taxonomy. We find extensive evidence for interspecific gene flow throughout the radiation. Hybridization has given rise to species of mixed ancestry. A 240 kilobase haplotype encompassing the ALX1 gene that encodes a transcription factor affecting craniofacial development is strongly associated with beak shape diversity across Darwin's finch species as well as within the medium ground finch (Geospiza fortis), a species that has undergone rapid evolution of beak shape in response to environmental changes. The ALX1 haplotype has contributed to diversification of beak shapes among the Darwin's finches and, thereby, to an expanded utilization of food resources.

Linking variation in quantitative traits to variation in the genome is an important, but challenging task in the study of life-history evolution. Linkage maps provide a valuable tool for the unravelling of such trait-gene associations. Moreover, they give insight into recombination landscapes and between-species karyotype evolution. Here we used genotype data, generated from a 10k single-nucleotide polymorphism (SNP) chip, of over 2000 individuals to produce high-density linkage maps of the great tit (Parus major), a passerine bird that serves as a model species for ecological and evolutionary questions. We created independent maps from two distinct populations: a captive F2-cross from The Netherlands (NL) and a wild population from the United Kingdom (UK). The two maps contained 6554 SNPs in 32 linkage groups, spanning 2010 cM and 1917 cM for the NL and UK populations, respectively, and were similar in size and marker order. Subtle levels of heterochiasmy within and between chromosomes were remarkably consistent between the populations, suggesting that the local departures from sex-equal recombination rates have evolved. This key and surprising result would have been impossible to detect if only one population was mapped. A comparison with zebra finch Taeniopygia guttata, chicken Gallus gallus and the green anole lizard Anolis carolinensis genomes provided further insight into the evolution of avian karyotypes.

Quantitative genetic analyses of basal metabolic rate (BMR) can inform us about the evolvability of the trait by providing estimates of heritability, and also of genetic correlations with other traits that may constrain the ability of BMR to respond to selection. Here, we studied a captive population of zebra finches (Taeniopygia guttata) in which selection lines for male courtship rate have been established. We measure BMR in these lines to see whether selection on male sexual activity would change BMR as a potentially correlated trait. We find that the genetic correlation between courtship rate and BMR is practically zero, indicating that the two traits can evolve independently of each other. Interestingly, we find that the heritability of BMR in our population (h(2)=0.45) is markedly higher than was previously reported for a captive zebra finch population from Norway. A comparison of the two studies shows that additive genetic variance in BMR has been largely depleted in the Norwegian population, especially the genetic variance in BMR that is independent of body mass. In our population, the slope of BMR increase with body mass differs not only between the sexes but also between the six selection lines, which we tentatively attribute to genetic drift and/or founder effects being strong in small populations. Our study therefore highlights two things. First, the evolvability of BMR may be less constrained by genetic correlations and lack of independent genetic variation than previously described. Second, genetic drift in small populations can rapidly lead to different evolvabilities across populations.

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

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.

We provide a first-generation linkage map of the blue tit (Cyanistes caeruleus), a passerine within the previously genetically uncharacterized family Paridae, which includes 91 orthologous loci with a single anchored position in the chicken (Gallus gallus) sequence assembly. The map consists of 18 linkage groups and covers 935 cM. There was highly conserved synteny between blue tit and chicken with the exception of a split on chromosome 1, potential splits on chromosome 4 and the translocation of two markers from chromosome 2 and 3, respectively, to chromosome 5. Gene order was very well conserved for the majority of chromosomes, an exception being chromosome 1 where multiple rearrangements were detected. Similar results were obtained in a comparison to the zebra finch (Taeniopygia guttata) genome assembly. The recombination rate in females was slightly higher than in males, implying a moderate degree of heterochiasmy in the blue tit. The map distance of the blue tit was approximately 78% of that of the Wageningen chicken broiler population, and very similar to the Uppsala chicken mapping population, over homologous genome regions. Apart from providing insights into avian recombination and genome evolution, our blue tit linkage map forms a valuable genetic resource for ecological and evolutionary research in Paridae.

A classic textbook example of adaptive radiation under natural selection is the evolution of 14 closely related species of Darwin's finches (Fringillidae, Passeriformes), whose primary diversity lies in the size and shape of their beaks. Thus, ground finches have deep and wide beaks, cactus finches have long and pointed beaks (low depth and narrower width), and warbler finches have slender and pointed beaks, reflecting differences in their respective diets. Previous work has shown that even small differences in any of the three major dimensions (depth, width and length) of the beak have major consequences for the overall fitness of the birds. Recently we used a candidate gene approach to explain one pathway involved in Darwin's finch beak morphogenesis. However, this type of analysis is limited to molecules with a known association with craniofacial and/or skeletogenic development. Here we use a less constrained, complementary DNA microarray analysis of the transcripts expressed in the beak primordia to find previously unknown genes and pathways whose expression correlates with specific beak morphologies. We show that calmodulin (CaM), a molecule involved in mediating Ca2+ signalling, is expressed at higher levels in the long and pointed beaks of cactus finches than in more robust beak types of other species. We validated this observation with in situ hybridizations. When this upregulation of the CaM-dependent pathway is artificially replicated in the chick frontonasal prominence, it causes an elongation of the upper beak, recapitulating the beak morphology of the cactus finches. Our results indicate that local upregulation of the CaM-dependent pathway is likely to have been a component of the evolution of Darwin's finch species with elongated beak morphology and provide a mechanistic explanation for the independence of beak evolution along different axes. More generally, our results implicate the CaM-dependent pathway in the developmental regulation of craniofacial skeletal structures.

Sperm design and function are important determinants of male reproductive success and are expected to be under strong selection. The way that spermatozoa phenotypes evolve is poorly understood, because there have been few studies of the quantitative genetics of sperm. Here we show, in the zebra finch Taeniopygia guttata, an extraordinary degree of inter-male variation in sperm design that is independent of sperm swimming velocity. A quantitative genetics study using data from over 900 zebra finches in a complex breeding experiment showed that sperm head, mid-piece and flagellum length are heritable, that negative genetic correlations exist between sperm traits, and that significant indirect (maternal) genetic effects exist. Selection on the zebra finch sperm phenotype may be low because sperm competition is infrequent in this species, and this, in combination with negative genetic correlations and maternal genetic effects, may account for the variation in sperm phenotype between males. These results have important implications for the evolution of sperm in other taxa.