The genome of a songbird

A novel neo-sex chromosome in Sylvietta brachyura (Macrosphenidae) adds to the extraordinary avian sex chromosome diversity among Sylvioidea songbirds

We report the discovery of a novel neo-sex chromosome in an African warbler, Sylvietta brachyura (northern crombec; Macrosphenidae). This species is part of the Sylvioidea superfamily, where four separate autosome-sex chromosome translocation events have previously been discovered via comparative genomics of 11 of the 22 families in this clade. Our discovery here resulted from analyses of genomic data of single species-representatives from three additional Sylvioidea families (Macrosphenidae, Pycnonotidae and Leiothrichidae). In all three species, we confirmed the translocation of a part of chromosome 4A to the sex chromosomes, which originated basally in Sylvioidea. In S. brachyura, we found that a part of chromosome 8 has been translocated to the sex chromosomes, forming a unique neo-sex chromosome in this lineage. Furthermore, the non-recombining part of 4A in S. brachyura is smaller than in other Sylvioidea species, which suggests that recombination continued along this region after the fusion event in the Sylvioidea ancestor. These findings reveal additional sex chromosome diversity among the Sylvioidea, where five separate translocation events are now confirmed.

The Rattlesnake W Chromosome: A GC-Rich Retroelement Refugium with Retained Gene Function Across Ancient Evolutionary Strata

Sex chromosomes diverge after the establishment of recombination suppression, resulting in differential sex-linkage of genes involved in genetic sex determination and dimorphic traits. This process produces systems of male or female heterogamety wherein the Y and W chromosomes are only present in one sex and are often highly degenerated. Sex-limited Y and W chromosomes contain valuable information about the evolutionary transition from autosomes to sex chromosomes, yet detailed characterizations of the structure, composition, and gene content of sex-limited chromosomes are lacking for many species. In this study, we characterize the female-specific W chromosome of the prairie rattlesnake (Crotalus viridis) and evaluate how recombination suppression and other processes have shaped sex chromosome evolution in ZW snakes. Our analyses indicate that the rattlesnake W chromosome is over 80% repetitive and that an abundance of GC-rich mdg4 elements has driven an overall high degree of GC-richness despite a lack of recombination. The W chromosome is also highly enriched for repeat sequences derived from endogenous retroviruses and likely acts as a "refugium" for these and other retroelements. We annotated 219 putatively functional W-linked genes across at least two evolutionary strata identified based on estimates of sequence divergence between Z and W gametologs. The youngest of these strata is relatively gene-rich, however gene expression across strata suggests retained gene function amidst a greater degree of degeneration following ancient recombination suppression. Functional annotation of W-linked genes indicates a specialization of the W chromosome for reproductive and developmental function since recombination suppression from the Z chromosome.

A conserved molecular template underlies color pattern diversity in estrildid finches

The color patterns that adorn animals' coats not only exhibit extensive diversity linked to various ecological functions but also display recurrences in geometry, orientation, or body location. How processes of pattern formation shape such phenotypic trends remains a mystery. Here, we surveyed plumage color patterns in passerine finches displaying extreme apparent variation and identified a conserved set of color domains. We linked these domains to putative embryonic skin regions instructed by early developmental tissues and outlined by the combinatory expression of few genetic markers. We found that this embryonic prepattern is largely conserved in birds displaying drastic color differences in the adult, interspecies variation resulting from the masking or display of each domain depending on their coloration. This work showed that a simple molecular landscape serves as common spatial template to extensive color pattern variation in finches, revealing that early conserved landmarks and molecular pathways are a major cause of phenotypic trends.

Widespread genomic signatures of reproductive isolation and sex-specific selection in the Eastern Yellow Robin, Eopsaltria australis

How new species evolve is one of the most fundamental questions in biology. Population divergence, which may lead to speciation, may be occurring in the Eastern Yellow Robin, a common passerine that lives along the eastern coast of Australia. This species is composed of 2 parapatric lineages that have highly divergent mitochondrial DNA; however, similar levels of divergence have not been observed in the nuclear genome. Here we re-examine the nuclear genomes of these mitolineages to test potential mechanisms underlying the discordance between nuclear and mitochondrial divergence. We find that nuclear admixture occurs in a narrow hybrid zone, although the majority of markers across the genome show evidence of reproductive isolation between populations of opposing mitolineages. There is an 8 MB section of a previously identified putative neo-sex chromosome that is highly diverged between allopatric but not parapatric populations, which may be the result of a chromosomal inversion. The neo-sex chromosomal nature of this region, as well as the geographic patterns in which it exhibits divergence, suggest it is unlikely to be contributing to reproductive isolation through mitonuclear incompatibilities as reported in earlier studies. In addition, there are sex differences in the number of markers that are differentiated between populations of opposite mitolineages, with greater differentiation occurring in females, which are heterozygous, than males. These results suggest that, despite the absence of previously observed assortative mating, mitolineages of Eastern Yellow Robin experience at least some postzygotic isolation from each other, in a pattern consistent with Haldane’s Rule.

Benchmarking ultra-high molecular weight DNA preservation methods for long-read and long-range sequencing

BACKGROUND

Studies in vertebrate genomics require sampling from a broad range of tissue types, taxa, and localities. Recent advancements in long-read and long-range genome sequencing have made it possible to produce high-quality chromosome-level genome assemblies for almost any organism. However, adequate tissue preservation for the requisite ultra-high molecular weight DNA (uHMW DNA) remains a major challenge. Here we present a comparative study of preservation methods for field and laboratory tissue sampling, across vertebrate classes and different tissue types.

RESULTS

We find that storage temperature was the strongest predictor of uHMW fragment lengths. While immediate flash-freezing remains the sample preservation gold standard, samples preserved in 95% EtOH or 20-25% DMSO-EDTA showed little fragment length degradation when stored at 4°C for 6 hours. Samples in 95% EtOH or 20-25% DMSO-EDTA kept at 4°C for 1 week after dissection still yielded adequate amounts of uHMW DNA for most applications. Tissue type was a significant predictor of total DNA yield but not fragment length. Preservation solution had a smaller but significant influence on both fragment length and DNA yield.

CONCLUSION

We provide sample preservation guidelines that ensure sufficient DNA integrity and amount required for use with long-read and long-range sequencing technologies across vertebrates. Our best practices generated the uHMW DNA needed for the high-quality reference genomes for phase 1 of the Vertebrate Genomes Project, whose ultimate mission is to generate chromosome-level reference genome assemblies of all ∼70,000 extant vertebrate species.

Improved zebra finch brain transcriptome identifies novel proteins with sex differences

The zebra finch (Taeniopygia guttata), a representative oscine songbird species, has been widely studied to investigate behavioral neuroscience, most notably the neurobiological basis of vocal learning, a rare trait shared in only a few animal groups including humans. In 2019, an updated zebra finch genome annotation (bTaeGut1_v1.p) was released from the Ensembl database and is substantially more comprehensive than the first version published in 2010. In this study, we utilized the publicly available RNA-seq data generated from Illumina-based short-reads and PacBio single-molecule real-time (SMRT) long-reads to assess the bird transcriptome. To analyze the high-throughput RNA-seq data, we adopted a hybrid bioinformatic approach combining short and long-read pipelines. From our analysis, we added 220 novel genes and 8,134 transcript variants to the Ensembl annotation, and predicted a new proteome based on the refined annotation. We further validated 18 different novel proteins by using mass-spectrometry data generated from zebra finch caudal telencephalon tissue. Our results provide additional resources for future studies of zebra finches utilizing this improved bird genome annotation and proteome.

De Novo Whole-Genome Sequencing and Assembly of the Yellow-Throated Bunting (Emberiza elegans) Provides Insights into Its Evolutionary Adaptation

Simple Summary

We report the genomic sequence of Emberiza elegans for understanding the evolutionary mechanisms of environmental adaptation and for studying a more effective genetic monitoring of this species. The E. elegans assembly was approximately 1.14 Gb, with a scaffold N50 of 28.94 Mb. About 15,868 protein-coding genes were predicted, and 16.62% of the genome was identified as having repetitive elements. Our genomic evolution analyses found considerable numbers of adaptive genes that may help the yellow-throated bunting cope with migratory behavior and environmental stressors of diseases. These results provide us with new insights into genomic evolution and adaptation, thus providing a valuable resource for further studies of population genetic diversity and genome evolution in this species.

Abstract

Yellow-throated bunting is a small migratory songbird unique to the Palearctic region. However, the genetic studies of this species remain limited, with no nuclear genomic sequence reported to date. In this study, the genomic DNA from the bird was sequenced in long reads using Nanopore sequencing technology. Combining short-read sequencing, the genome was well-assembled and annotated. The final length of the assembly is approximately 1.14 Gb, with a scaffold N50 of 28.94 Mb. About 15,868 protein-coding genes were predicted, and 16.62% of the genome was identified as having repetitive elements. Comparative genomic analysis showed numerous expanded gene families and positively selected genes significantly enriched in those KEGG pathways that are associated with migratory behavior adaptation and immune response. Here, this newly generated de novo genome of the yellow-throated bunting using long reads provide the research community with a valuable resource for further studies of population genetic diversity and genome evolution in this species.

The evolutionary history and mechanistic basis of female ornamentation in a tropical songbird

Ornamentation, such as the showy plumage of birds, is widespread among female vertebrates, yet the evolutionary pressures shaping female ornamentation remain uncertain. In part this is due to a poor understanding of the mechanistic route to ornamentation in females. To address this issue, we evaluated the evolutionary history of ornament expression in a tropical passerine bird, the White-shouldered Fairywren, whose females, but not males, strongly vary between populations in occurrence of ornamented black-and-white plumage. We first use phylogenomic analysis to demonstrate that female ornamentation is derived and that female ornamentation evolves independently of changes in male plumage. We then use exogenous testosterone in a field experiment to induce partial ornamentation in naturally unornamented females. By sequencing the transcriptome of experimentally induced ornamented and natural feathers, we identify genes expressed during ornament production and evaluate the degree to which female ornamentation in this system is associated with elevated testosterone, as is common in males. We reveal that some ornamentation in females is linked to testosterone and that sexes differ in ornament-linked gene expression. Lastly, using genomic outlier analysis we identify a candidate melanogenesis gene that lies in a region of high genomic divergence among populations that is also differentially expressed in feather follicles of different female plumages. Taken together, these findings are consistent with sex-specific selection favoring the evolution of female ornaments and demonstrate a key role for testosterone in generating population divergence in female ornamentation through gene regulation. More broadly, our work highlights similarities and differences in how ornamentation evolves in the sexes.

Fine-scale population genetic structure and barriers to gene flow in a widespread seabird (Ardenna pacifica)

Pelagic seabirds are highly mobile, reducing opportunities for population isolation that might promote differentiation and speciation. At the same time, many seabirds are philopatric, and their tendency to return to their natal islands to breed might reduce gene flow sufficiently to permit local adaptation and differentiation. To test the net impact of these competing processes, estimates of differentiation and gene flow based on comprehensive geographical sampling are required. We leveraged diverse source material to achieve comprehensive geographical sampling in a widespread seabird, the Wedge-tailed Shearwater (Ardenna pacifica). Using data from sequence capture and high-throughput sequencing of 2402 loci containing 20 780 single nucleotide polymorphisms, we tested for population differentiation and gene flow among breeding areas. We found little evidence of deep divergences within A. pacifica but were able to resolve fine-scale differentiation across island groups. This differentiation was sufficient to assign individuals sampled away from breeding areas to their likely source populations. Estimated effective migration surfaces revealed reduced migration between the Indian Ocean and Pacific Ocean, presumably owing to land barriers, and across the equatorial Pacific Ocean, perhaps associated with differences in breeding schedule. Our results reveal that, despite their mobility, pelagic seabirds can exhibit fine-scale population differentiation and reduced gene flow among ocean basins.

Extreme variation in recombination rate and genetic diversity along the Sylvioidea neo-sex chromosome

Recombination strongly impacts sequence evolution by affecting the extent of linkage and the efficiency of selection. Here, we study recombination over the Z chromosome in great reed warblers (Acrocephalus arundinaceus) using pedigree-based linkage mapping. This species has extended Z and W chromosomes ("neo-sex chromosomes") formed by a fusion between a part of chromosome 4A and the ancestral sex chromosomes, which provides a unique opportunity to assess recombination and sequence evolution in sex-linked regions of different ages. We assembled an 87.54 Mbp and 90.19 cM large Z with a small pseudoautosomal region (0.89 Mbp) at one end and the fused Chr4A-part at the other end of the chromosome. A prominent feature in our data was an extreme variation in male recombination rate along Z with high values at both chromosome ends, but an apparent lack of recombination over a substantial central section, covering 78% of the chromosome. The nonrecombining region showed a drastic loss of genetic diversity and accumulation of repeats compared to the recombining parts. Thus, our data emphasize a key role of recombination in affecting local levels of polymorphism. Nonetheless, the evolutionary rate of genes (dN/dS) did not differ between high and low recombining regions, suggesting that the efficiency of selection on protein-coding sequences can be maintained also at very low levels of recombination. Finally, the Chr4A-derived part showed a similar recombination rate as the part of the ancestral Z that did recombine, but its sequence characteristics reflected both its previous autosomal, and current Z-linked, recombination patterns.

A high-quality genome assembly and annotation of the dark-eyed junco Junco hyemalis, a recently diversified songbird

The dark-eyed junco (Junco hyemalis) is one of the most common passerines of North America, and has served as a model organism in studies related to ecophysiology, behavior, and evolutionary biology for over a century. It is composed of at least 6 distinct, geographically structured forms of recent evolutionary origin, presenting remarkable variation in phenotypic traits, migratory behavior, and habitat. Here, we report a high-quality genome assembly and annotation of the dark-eyed junco generated using a combination of shotgun libraries and proximity ligation Chicago and Dovetail Hi-C libraries. The final assembly is ∼1.03 Gb in size, with 98.3% of the sequence located in 30 full or nearly full chromosome scaffolds, and with a N50/L50 of 71.3 Mb/5 scaffolds. We identified 19,026 functional genes combining gene prediction and similarity approaches, of which 15,967 were associated to GO terms. The genome assembly and the set of annotated genes yielded 95.4% and 96.2% completeness scores, respectively when compared with the BUSCO avian dataset. This new assembly for J. hyemalis provides a valuable resource for genome evolution analysis, and for identifying functional genes involved in adaptive processes and speciation.

How thermal challenges change gene regulation in the songbird brain and gonad: implications for sexual selection in our changing world

In a rapidly warming world, exposure to high temperatures may impact fitness, but the gene regulatory mechanisms that link sublethal heat to sexually selected traits are not well understood, particularly in endothermic animals. Our experiment used zebra finches (Taeniopygia guttata), songbirds that experience extreme temperature fluctuations in their native Australia. We exposed captive males to an acute thermal challenge (43°C) compared with thermoneutral (35°C) and lower (27°C) temperatures. We found significantly more heat dissipation behaviors at 43°C, a temperature previously shown to reduce song production and fertility, and more heat retention behaviors at 27°C. Next, we characterized transcriptomic responses in tissues important for mating effort - the posterior telencephalon, for its role in song production, and the testis, for its role in fertility and hormone production. Differential expression of hundreds of genes in the testes, but few in the brain, suggest the brain is less responsive to extreme temperatures. Nevertheless, gene network analyses revealed that expression related to dopaminergic signaling in the brain co-varied with heat dissipation behaviors, providing a mechanism by which temporary thermal challenges may alter motivational circuits for song production. In both brain and testis, we observed correlations between thermally sensitive gene networks and individual differences in thermoregulatory behavior. Although we cannot directly relate these gene regulatory changes to mating success, our results suggest that individual variation in response to thermal challenges could impact sexually selected traits in a warming world.

Genetic architecture facilitates then constrains adaptation in a host-parasite coevolutionary arms race

In coevolutionary arms races, interacting species impose selection on each other, generating reciprocal adaptations and counter adaptations. This process is typically enhanced by genetic recombination and heterozygosity, but these sources of evolutionary novelty may be secondarily lost when uniparental inheritance evolves to ensure the integrity of sex-linked adaptations. We demonstrate that host-specific egg mimicry in the African cuckoo finch Anomalospiza imberbis is maternally inherited, confirming the validity of an almost century-old hypothesis. We further show that maternal inheritance not only underpins the mimicry of different host species but also additional mimetic diversification that approximates the range of polymorphic egg “signatures” that have evolved within host species as an escalated defense against parasitism. Thus, maternal inheritance has enabled the evolution and maintenance of nested levels of mimetic specialization in a single parasitic species. However, maternal inheritance and the lack of sexual recombination likely disadvantage cuckoo finches by stifling further adaptation in the ongoing arms races with their individual hosts, which we show have retained biparental inheritance of egg phenotypes. The inability to generate novel genetic combinations likely prevents cuckoo finches from mimicking certain host phenotypes that are currently favored by selection (e.g., the olive-green colored eggs laid by some tawny-flanked prinia, Prinia subflava, females). This illustrates an important cost of coding coevolved adaptations on the nonrecombining sex chromosome, which may impede further coevolutionary change by effectively reversing the advantages of sexual reproduction in antagonistic coevolution proposed by the Red Queen hypothesis.

Genomic associations with poxvirus across divergent island populations in Berthelot's pipit

Understanding the mechanisms and genes that enable animal populations to adapt to pathogens is important from an evolutionary, health and conservation perspective. Berthelot's pipit (Anthus berthelotii) experiences extensive and consistent spatial heterogeneity in avian pox infection pressure across its range of island populations, thus providing an excellent system with which to examine how pathogen-mediated selection drives spatial variation in immunogenetic diversity. Here we test for evidence of genetic variation associated with avian pox at both an individual and population-level. At the individual level, we find no evidence that variation in MHC class I and TLR4 (both known to be important in recognising viral infection) was associated with pox infection within two separate populations. However, using genotype-environment association (Bayenv) in conjunction with genome-wide (ddRAD-seq) data, we detected strong associations between population-level avian pox prevalence and allele frequencies of single nucleotide polymorphisms (SNPs) at a number of sites across the genome. These sites were located within genes involved in cellular stress signalling and immune responses, many of which have previously been associated with responses to viral infection in humans and other animals. Consequently, our analyses indicates that pathogen-mediated selection may play a role in shaping genomic variation among relatively recently colonised island bird populations and highlights the utility of genotype-environment associations for identifying candidate genes potentially involved in host-pathogen interactions.

Chromosome-length genome assembly and linkage map of a critically endangered Australian bird: the helmeted honeyeater

BACKGROUND

The helmeted honeyeater (Lichenostomus melanops cassidix) is a Critically Endangered bird endemic to Victoria, Australia. To aid its conservation, the population is the subject of genetic rescue. To understand, monitor, and modulate the effects of genetic rescue on the helmeted honeyeater genome, a chromosome-length genome and a high-density linkage map are required.

RESULTS

We used a combination of Illumina, Oxford Nanopore, and Hi-C sequencing technologies to assemble a chromosome-length genome of the helmeted honeyeater, comprising 906 scaffolds, with length of 1.1 Gb and scaffold N50 of 63.8 Mb. Annotation comprised 57,181 gene models. Using a pedigree of 257 birds and 53,111 single-nucleotide polymorphisms, we obtained high-density linkage and recombination maps for 25 autosomes and Z chromosome. The total sex-averaged linkage map was 1,347 cM long, with the male map being 6.7% longer than the female map. Recombination maps revealed sexually dimorphic recombination rates (overall higher in males), with average recombination rate of 1.8 cM/Mb. Comparative analyses revealed high synteny of the helmeted honeyeater genome with that of 3 passerine species (e.g., 32 Hi-C scaffolds mapped to 30 zebra finch autosomes and Z chromosome). The genome assembly and linkage map suggest that the helmeted honeyeater exhibits a fission of chromosome 1A into 2 chromosomes relative to zebra finch. PSMC analysis showed a ∼15-fold decline in effective population size to ∼60,000 from mid- to late Pleistocene.

CONCLUSIONS

The annotated chromosome-length genome and high-density linkage map provide rich resources for evolutionary studies and will be fundamental in guiding conservation efforts for the helmeted honeyeater.

The genomic architecture of the passerine MHC region: high repeat content and contrasting evolutionary histories of single copy and tandemly duplicated MHC genes

The Major Histocompatibility Complex (MHC) is of central importance to the immune system, and an optimal MHC diversity is believed to maximize pathogen elimination. Birds show substantial variation in MHC diversity, ranging from few genes in most bird orders to very many genes in passerines. Our understanding of the evolutionary trajectories of the MHC in passerines is hampered by lack of data on genomic organization. Therefore, we assemble and annotate the MHC genomic region of the great reed warbler (Acrocephalus arundinaceus), using long-read sequencing and optical mapping. The MHC region is large (>5.5Mb), characterized by structural changes compared to hitherto investigated bird orders and shows higher repeat content than the genome average. These features were supported by analyses in three additional passerines. MHC genes in passerines are found in two different chromosomal arrangements, either as single copy MHC genes located among non-MHC genes, or as tandemly duplicated tightly linked MHC genes. Some single copy MHC genes are old and putative orthologs among species. In contrast tandemly duplicated MHC genes are monophyletic within species and have evolved by simultaneous gene duplication of several MHC genes. Structural differences in the MHC genomic region among bird orders seem substantial compared to mammals and have possibly been fuelled by clade-specific immune system adaptations. Our study provides methodological guidance in characterizing complex genomic regions, constitutes a resource for MHC research in birds, and calls for a revision of the general belief that avian MHC has a conserved gene order and small size compared to mammals.

Induction of an immortalized songbird cell line allows for gene characterization and knockout by CRISPR-Cas9

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.

Varied diversification patterns and distinct demographic trajectories in Ethiopian montane forest bird (Aves: Passeriformes) populations separated by the Great Rift Valley

Taxon-specific characteristics and extrinsic climatic and geological forces may both shape population differentiation and speciation. In geographically and taxonomically focused investigations, differentiation may occur synchronously as species respond to the same external conditions. Conversely, when evolution is investigated in taxa with largely varying traits, population differentiation and speciation is complex and shaped by interactions of Earth's template and species-specific traits. As such, it is important to characterize evolutionary histories broadly across the tree of life, especially in geographic regions that are exceptionally diverse and under pressures from human activities such as in biodiversity hotspots. Here, using whole-genome sequencing data, we characterize genomic variation in populations of six Ethiopian Highlands forest bird species separated by a lowland biogeographic barrier, the Great Rift Valley (GRV). In all six species, populations on either side of the GRV exhibited significant but varying levels of genetic differentiation. Species' dispersal ability was negatively correlated with levels of population differentiation. Isolation with migration models indicated varied patterns of population differentiation and connectivity among populations of the focal species. We found that demographic histories-estimated for each individual-varied by both species and population but were consistent between individuals of the same species and sampling region. We found that genomic diversity varied by half an order of magnitude across species, and that this variation could largely be explained by the harmonic mean of effective population size over the past 200,000 years. Overall, we found that even in highly dispersive species like birds, the GRV acts as a substantial biogeographic barrier.

Reciprocal processes of sensory perception and social bonding: an integrated social-sensory framework of social behavior

Organisms filter the complexity of natural stimuli through their individual sensory and perceptual systems. Such perceptual filtering is particularly important for social stimuli. A shared "social umwelt" allows individuals to respond appropriately to the expected diversity of cues and signals during social interactions. In this way, the behavioral and neurobiological mechanisms of sociality and social bonding cannot be disentangled from perceptual mechanisms and sensory processing. While a degree of embeddedness between social and sensory processes is clear, our dominant theoretical frameworks favor treating the social and sensory processes as distinct. An integrated social-sensory framework has the potential to greatly expand our understanding of the mechanisms underlying individual variation in social bonding and sociality more broadly. Here we leverage what is known about sensory processing and pair bonding in two common study systems with significant species differences in their umwelt (rodent chemosensation and avian acoustic communication). We primarily highlight that (1) communication is essential for pair bond formation and maintenance, (2) the neural circuits underlying perception, communication and social bonding are integrated, and (3) candidate neuromodulatory mechanisms that regulate pair bonding also impact communication and perception. Finally, we propose approaches and frameworks that more fully integrate sensory processing, communication, and social bonding across levels of analysis: behavioral, neurobiological, and genomic. This perspective raises two key questions: (1) how is social bonding shaped by differences in sensory processing?, and (2) to what extent is sensory processing and the saliency of signals shaped by social interactions and emerging relationships?

Inhibition of miR-128 Enhances Vocal Sequence Organization in Juvenile Songbirds

The molecular mechanisms underlying learned vocal communication are not well characterized. This is a major barrier for developing treatments for conditions affecting social communication, such as autism spectrum disorder (ASD). Our group previously generated an activity-dependent gene expression network in the striatopallidal song control nucleus, Area X, in adult zebra finches to identify master regulators of learned vocal behavior. This dataset revealed that the two host genes for microRNA-128, ARPP21 and R3HDM1, are among the top genes whose expression correlates to how much birds sing. Here we examined whether miR-128 itself is behaviorally regulated in Area X and found that its levels decline with singing. We hypothesized that reducing miR-128 during the critical period for vocal plasticity would enhance vocal learning. To test this, we bilaterally injected an antisense miR-128 construct (AS miR-128) or a control scrambled sequence into Area X at post-hatch day 30 (30 d) using sibling-matched experimental and control pupils. The juveniles were then returned to their home cage and raised with their tutors. Strikingly, inhibition of miR-128 in young birds enhanced the organization of learned vocal sequences. Tutor and pupil stereotypy scores were positively correlated, though the correlation was stronger between tutors and control pupils compared to tutors and AS miR-128 pupils. This difference was driven by AS miR-128 pupils achieving higher stereotypy scores despite their tutors’ lower syntax scores. AS miR-128 birds with tutors on the higher end of the stereotypy spectrum were more likely to produce songs with faster tempos relative to sibling controls. Our results suggest that low levels of miR-128 facilitate vocal sequence stereotypy. By analogy, reducing miR-128 could enhance the capacity to learn to speak in patients with non-verbal ASD. To our knowledge, this study is the first to directly link miR-128 to learned vocal communication and provides support for miR-128 as a potential therapeutic target for ASD.

Evaluation of four methods to identify the homozygotic sex chromosome in small populations

BACKGROUND

Whole genomes are commonly assembled into a collection of scaffolds and often lack annotations of autosomes, sex chromosomes, and organelle genomes (i.e., mitochondrial and chloroplast). As these chromosome types differ in effective population size and can have highly disparate evolutionary histories, it is imperative to take this information into account when analysing genomic variation. Here we assessed the accuracy of four methods for identifying the homogametic sex chromosome in a small population using two whole genome sequences (WGS) and 133 RAD sequences of white-tailed eagles (Haliaeetus albicilla): i) difference in read depth per scaffold in a male and a female, ii) heterozygosity per scaffold in a male and a female, iii) mapping to the reference genome of a related species (chicken) with annotated sex chromosomes, and iv) analysis of SNP-loadings from a principal components analysis (PCA), based on the low-depth RADseq data.

RESULTS

The best performing approach was the reference mapping (method iii), which identified 98.12% of the expected homogametic sex chromosome (Z). Read depth per scaffold (method i) identified 86.41% of the homogametic sex chromosome with few false positives. SNP-loading scores (method iv) identified 78.6% of the Z-chromosome and had a false positive discovery rate of more than 10%. Heterozygosity per scaffold (method ii) did not provide clear results due to a lack of diversity in both the Z and autosomal chromosomes, and potential interference from the heterogametic sex chromosome (W). The evaluation of these methods also revealed 10 Mb of putative PAR and gametologous regions.

CONCLUSION

Identification of the homogametic sex chromosome in a small population is best accomplished by reference mapping or examining differences in read depth between sexes.

Recurrent chromosome reshuffling and the evolution of neo-sex chromosomes in parrots

The karyotype of most birds has remained considerably stable during more than 100 million years' evolution, except for some groups, such as parrots. The evolutionary processes and underlying genetic mechanism of chromosomal rearrangements in parrots, however, are poorly understood. Here, using chromosome-level assemblies of four parrot genomes, we uncover frequent chromosome fusions and fissions, with most of them occurring independently among lineages. The increased activities of chromosomal rearrangements in parrots are likely associated with parrot-specific loss of two genes, ALC1 and PARP3, that have known functions in the repair of double-strand breaks and maintenance of genome stability. We further find that the fusion of the ZW sex chromosomes and chromosome 11 has created a pair of neo-sex chromosomes in the ancestor of parrots, and the chromosome 25 has been further added to the sex chromosomes in monk parakeet. Together, the combination of our genomic and cytogenetic analyses characterizes the complex evolutionary history of chromosomal rearrangements and sex chromosomes in parrots.

Occasional paternal inheritance of the germline-restricted chromosome in songbirds

Songbirds have one special accessory chromosome, the so-called germline-restricted chromosome (GRC), which is only present in germline cells and absent from all somatic tissues. Earlier work on the zebra finch (Taeniopygia guttata castanotis) showed that the GRC is inherited only through the female line-like the mitochondria-and is eliminated from the sperm during spermatogenesis. Here, we show that the GRC has the potential to be paternally inherited. Confocal microscopy using GRC-specific fluorescent in situ hybridization probes indicated that a considerable fraction of sperm heads (1 to 19%) in zebra finch ejaculates still contained the GRC. In line with these cytogenetic data, sequencing of ejaculates revealed that individual males from two families differed strongly and consistently in the number of GRCs in their ejaculates. Examining a captive-bred male hybrid of the two zebra finch subspecies (T. g. guttata and T. g. castanotis) revealed that the mitochondria originated from a castanotis mother, whereas the GRC came from a guttata father. Moreover, analyzing GRC haplotypes across nine castanotis matrilines, estimated to have diverged for up to 250,000 y, showed surprisingly little variability among GRCs. This suggests that a single GRC haplotype has spread relatively recently across all examined matrilines. A few diagnostic GRC mutations that arose since this inferred spreading suggest that the GRC has continued to jump across matriline boundaries. Our findings raise the possibility that certain GRC haplotypes could selfishly spread through the population via occasional paternal transmission, thereby outcompeting other GRC haplotypes that were limited to strict maternal inheritance, even if this was partly detrimental to organismal fitness.

Neo-sex chromosome evolution and phenotypic differentiation across an elevational gradient in horned larks (Eremophila Alpestris)

Genetic structure and phenotypic variation among populations is affected by both geographic distance and environmental variation across species' distributions. Understanding the relative contributions of isolation by distance (IBD) and isolation by environment (IBE) is important for elucidating population dynamics across habitats and ecological gradients. In this study, we compared phenotypic and genetic variation among Horned Lark (Eremophila alpestris) populations from 10 sites encompassing an elevational gradient from low-elevation desert scrub in Death Valley (285 a.s.l.) to high-elevation meadows in the White Mountains of the Sierra Nevada of California (greater than 3000 m a.s.l.). Using a ddRAD dataset of 28,474 SNPs aligned to a high-quality reference genome, we compared genetic structure with elevational, environmental, and spatial distance to quantify how different aspects of the landscape drive genomic and phenotypic differentiation in Horned Larks. We found larger-bodied birds were associated with sites that had less seasonality and higher annual precipitation, and longer spurs occurred in soils with more clay and silt content, less sand, and finer fragments. Larks have large neo-sex chromosomes, and we found that associations with elevation and environmental variation were much stronger among neo-sex chromosomes compared to autosomes. Furthermore, we found that putative chromosomal translocations, fusions, and inversions were associated with elevation and may underlie local adaptation across an elevational gradient in Horned Larks. Our results suggest that genetic variation in Horned Larks is affected more by IBD than IBE, but specific phenotypes and genomic regions-particually on neo-sex chromosomes-bear stronger associations with the environment.

Zebra finch song is a very short-range signal in the wild: evidence from an integrated approach

Birdsong is typically seen as a long-range signal functioning in mate attraction and territory defense. Among birds, the zebra finch is the prime model organism in bioacoustics, yet almost exclusively studied in the lab. In the wild, however, zebra finch song differs strikingly from songbirds commonly studied in the wild as zebra finch males sing most after mating and in the absence of territoriality. Using data from the wild, we here provide an ecological context for a wealth of laboratory studies. By integrating calibrated sound recordings, sound transmission experiments and social ecology of zebra finches in the wild with insights from hearing physiology we show that wild zebra finch song is a very short-range signal with an audible range of about nine meters and that even the louder distance calls do not carry much farther (up to about fourteen meters). These integrated findings provide an ecological context for the interpretation of laboratory studies of this species and indicate that the vocal communication distance of the main laboratory species for avian acoustics contrasts strikingly with songbirds that use their song as a long-range advertisement signal.

Single-Cell RNA Sequencing Revealed the Heterogeneity of Gonadal Primordial Germ Cells in Zebra Finch (Taeniopygia guttata)

Primordial germ cells (PGCs) are undifferentiated gametes with heterogeneity, an evolutionarily conserved characteristic across various organisms. Although dynamic selection at the level of early germ cell populations is an important biological feature linked to fertility, the heterogeneity of PGCs in avian species has not been characterized. In this study, we sought to evaluate PGC heterogeneity in zebra finch using a single-cell RNA sequencing (scRNA-seq) approach. Using scRNA-seq of embryonic gonadal cells from male and female zebra finches at Hamburger and Hamilton (HH) stage 28, we annotated nine cell types from 20 cell clusters. We found that PGCs previously considered a single population can be separated into three subtypes showing differences in apoptosis, proliferation, and other biological processes. The three PGC subtypes were specifically enriched for genes showing expression patterns related to germness or pluripotency, suggesting functional differences in PGCs according to the three subtypes. Additionally, we discovered a novel biomarker, SMC1B, for gonadal PGCs in zebra finch. The results provide the first evidence of substantial heterogeneity in PGCs previously considered a single population in birds. This discovery expands our understanding of PGCs to avian species, and provides a basis for further research.

Avian Neo-Sex Chromosomes Reveal Dynamics of Recombination Suppression and W Degeneration

How the avian sex chromosomes first evolved from autosomes remains elusive as 100 million years (My) of divergence and degeneration obscure their evolutionary history. The Sylvioidea group of songbirds is interesting for understanding avian sex chromosome evolution because a chromosome fusion event ∼24 Ma formed "neo-sex chromosomes" consisting of an added (new) and an ancestral (old) part. Here, we report the complete female genome (ZW) of one Sylvioidea species, the great reed warbler (Acrocephalus arundinaceus). Our long-read assembly shows that the added region has been translocated to both Z and W, and whereas the added-Z has retained its gene order the added-W part has been heavily rearranged. Phylogenetic analyses show that recombination between the homologous added-Z and -W regions continued after the fusion event, and that recombination suppression across this region took several million years to be completed. Moreover, recombination suppression was initiated across multiple positions over the added-Z, which is not consistent with a simple linear progression starting from the fusion point. As expected following recombination suppression, the added-W show signs of degeneration including repeat accumulation and gene loss. Finally, we present evidence for nonrandom maintenance of slowly evolving and dosage-sensitive genes on both ancestral- and added-W, a process causing correlated evolution among orthologous genes across broad taxonomic groups, regardless of sex linkage.

Oxytocin variation and brain region-specific gene expression in a domesticated avian species

The Bengalese finch was domesticated more than 250 years ago from the wild white-rumped munia (WRM). Similar to other domesticated species, Bengalese finches show a reduced fear response and have lower corticosterone levels, compared to WRMs. Bengalese finches and munias also have different song types. Since oxytocin (OT) has been found to be involved in stress coping and auditory processing, we tested whether the OT sequence and brain expression pattern and content differ in wild munias and domesticated Bengalese finches. We sequenced the OT from 10 wild munias and 11 Bengalese finches and identified intra-strain variability in both the untranslated and protein-coding regions of the sequence, with all the latter giving rise to synonymous mutations. Several of these changes fall in specific transcription factor-binding sites, and show either a conserved or a relaxed evolutionary trend in the avian lineage, and in vertebrates in general. Although in situ hybridization in several hypothalamic nuclei did not reveal significant differences in the number of cells expressing OT between the two strains, real-time quantitative PCR showed a significantly higher OT mRNA expression in the cerebrum of the Bengalese finches relative to munias, but a significantly lower expression in their diencephalon. Our study thus points to a brain region-specific pattern of neurochemical expression in domesticated and wild avian strains, which could be linked to domestication and the behavioral changes associated with it.

Genome Stability Is in the Eye of the Beholder: CR1 Retrotransposon Activity Varies Significantly across Avian Diversity

Since the sequencing of the zebra finch genome it has become clear that avian genomes, while largely stable in terms of chromosome number and gene synteny, are more dynamic at an intrachromosomal level. A multitude of intrachromosomal rearrangements and significant variation in transposable element (TE) content have been noted across the avian tree. TEs are a source of genome plasticity, because their high similarity enables chromosomal rearrangements through nonallelic homologous recombination, and they have potential for exaptation as regulatory and coding sequences. Previous studies have investigated the activity of the dominant TE in birds, chicken repeat 1 (CR1) retrotransposons, either focusing on their expansion within single orders, or comparing passerines with nonpasserines. Here, we comprehensively investigate and compare the activity of CR1 expansion across orders of birds, finding levels of CR1 activity vary significantly both between and within orders. We describe high levels of TE expansion in genera which have speciated in the last 10 Myr including kiwis, geese, and Amazon parrots; low levels of TE expansion in songbirds across their diversification, and near inactivity of TEs in the cassowary and emu for millions of years. CR1s have remained active over long periods of time across most orders of neognaths, with activity at any one time dominated by one or two families of CR1s. Our findings of higher TE activity in species-rich clades and dominant families of TEs within lineages mirror past findings in mammals and indicate that genome evolution in amniotes relies on universal TE-driven processes.

Comparative analysis of transposable elements provides insights into genome evolution in the genus Camelus

BACKGROUND

Transposable elements (TEs) are common features in eukaryotic genomes that are known to affect genome evolution critically and to play roles in gene regulation. Vertebrate genomes are dominated by TEs, which can reach copy numbers in the hundreds of thousands. To date, details regarding the presence and characteristics of TEs in camelid genomes have not been made available.

RESULTS

We conducted a genome-wide comparative analysis of camelid TEs, focusing on the identification of TEs and elucidation of transposition histories in four species: Camelus dromedarius, C. bactrianus, C. ferus, and Vicugna pacos. Our TE library was created using both de novo structure-based and homology-based searching strategies ( https://github.com/kacst-bioinfo-lab/TE_ideintification_pipeline ). Annotation results indicated a similar proportion of each genomes comprising TEs (35-36%). Class I LTR retrotransposons comprised 16-20% of genomes, and mostly consisted of the endogenous retroviruses (ERVs) groups ERVL, ERVL-MaLR, ERV_classI, and ERV_classII. Non-LTR elements comprised about 12% of genomes and consisted of SINEs (MIRs) and the LINE superfamilies LINE1, LINE2, L3/CR1, and RTE clades. Least represented were the Class II DNA transposons (2%), consisting of hAT-Charlie, TcMar-Tigger, and Helitron elements and comprising about 1-2% of each genome.

CONCLUSIONS

The findings of the present study revealed that the distribution of transposable elements across camelid genomes is approximately similar. This investigation presents a characterization of TE content in four camelid to contribute to developing a better understanding of camelid genome architecture and evolution.

Cytogenetic Evidence Clarifies the Phylogeny of the Family Rhynchocyclidae (Aves: Passeriformes)

The phylogenetic position and taxonomic status of Rhynchocyclidae (Aves: Passeriformes) have been the subject of debate since their first description. In most models, Rhynchocyclidae represents a subfamily-level taxon placed within the Tyrant Flycatchers (Tyrannidae). Considering that this classification does not include cytotaxonomic characters, we tested the hypothesis that the chromosome organization of Rhynchocyclidae members differs from that of Tyrannidae. Hence, we selected two species, Tolmomyias sulphurescens, and Pitangus sulphuratus, representing Rhynchocyclidae and Tyrannidae, respectively. Results revealed a diploid number (2n) of 60 in T. sulphurescens and 2n = 80 in P. sulphuratus, indicating significant chromosomal differences. Chromosome mapping of Gallus gallus (GGA) and Taeniopygia guttata bacterial artificial chromosome (BAC) corresponding to chromosomes GGA1-28 (except 16) revealed that the genome evolution of T. sulphurescens involved extensive chromosome fusions of macrochromosomes and microchromosomes. On the other hand, P. sulphuratus retained the ancestral pattern of organization of macrochromosomes (except the centric fission involving GGA1) and microchromosomes. In conclusion, comparing our results with previous studies in Tyrant Flycatchers and allies indicates that P. sulphuratus has similar karyotypes to other Tyrannidae members. However, T. sulphurescens does not resemble the Tyrannidae family, reinforcing family status to the clade named Rhynchocyclidae.

First constitutive heterochromatin characterization and Karyotype of white stork Ciconia ciconia (Aves: Ciconiidae)

The karyotype and constitutive heterochromatin pattern of the white stork Ciconia ciconia samples obtained from Manzala lake, Dimiaat, Egypt was described. Somatic cells of Ciconia ciconia samples have diploid number 2n= 68 chromosomes. Out of 68 chromosomes, 11 pairs including sex chromosomes were macrochromosomes and the remaining pairs were microchromosomes. Of the 11 macrochromosome pairs, no.1, 2, 4 and 5 were submetacentric and pairs no. 6, 7 and 8 were described as metacentric. In addition, the autosome pair no.3 was subtelocentric, while autosome pair no.9 was acrocentric. Also, the sex chromosome Z represents the fourth one in size and it was classified as submetacentric while, W chromosome appeared as medium size and was acrocentric. Furthermore, C-banding pattern (constitutive heterochromatin) revealed variation in their sizes and occurrence between macrochromosomes. Pairs no. 7 and 8 of autosomes exhibited unusual distribution of heterochromatin, where they appeared as entirely heterochromatic. This may be related to the origin of sex chromosomes Z and W. However, there is no sufficient evidence illustrate the appearance of entirely heterochromatic autosomes. Therefore, there is no available cytogenetic literature that describes the C-banding and karyotype of Ciconia Ciconia, so the results herein are important and may assist in cytogenetic study and evolutionary pattern of Ciconiiformes.

The avian W chromosome is a refugium for endogenous retroviruses with likely effects on female-biased mutational load and genetic incompatibilities

It is a broadly observed pattern that the non-recombining regions of sex-limited chromosomes (Y and W) accumulate more repeats than the rest of the genome, even in species like birds with a low genome-wide repeat content. Here, we show that in birds with highly heteromorphic sex chromosomes, the W chromosome has a transposable element (TE) density of greater than 55% compared to the genome-wide density of less than 10%, and contains over half of all full-length (thus potentially active) endogenous retroviruses (ERVs) of the entire genome. Using RNA-seq and protein mass spectrometry data, we were able to detect signatures of female-specific ERV expression. We hypothesize that the avian W chromosome acts as a refugium for active ERVs, probably leading to female-biased mutational load that may influence female physiology similar to the 'toxic-Y' effect in Drosophila males. Furthermore, Haldane's rule predicts that the heterogametic sex has reduced fertility in hybrids. We propose that the excess of W-linked active ERVs over the rest of the genome may be an additional explanatory variable for Haldane's rule, with consequences for genetic incompatibilities between species through TE/repressor mismatches in hybrids. Together, our results suggest that the sequence content of female-specific W chromosomes can have effects far beyond sex determination and gene dosage. This article is part of the theme issue 'Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part II)'.

A sex chromosome inversion is associated with copy number variation of mitochondrial DNA in zebra finch sperm

The propulsion of sperm cells via movement of the flagellum is of vital importance for successful fertilization. While the exact mechanism of energy production for this movement varies between species, in avian species energy is thought to come predominantly from the mitochondria located in the sperm midpiece. Larger midpieces may contain more mitochondria, which should enhance the energetic capacity and possibly promote mobility. Due to an inversion polymorphism on their sex chromosome TguZ, zebra finches (Taeniopygia guttata castanotis) exhibit large within-species variation in sperm midpiece length, and those sperm with the longest midpieces swim the fastest. Here, we test through quantitative real-time PCR in zebra finch ejaculates whether the inversion genotype has an effect on the copy number of mitochondrial DNA (mtDNA). We find that zebra finches carrying the derived allele (correlated with longer sperm midpieces) have more copies of the mtDNA in their ejaculates than those homozygous for the ancestral allele (shorter midpieces). We suggest downstream effects of mtDNA copy number variation on the rate of adenosine triphosphate production, which in turn may influence sperm swimming speed and fertilization success. Central components of gamete energy metabolism may thus be the proximate cause for a fitness-relevant genetic polymorphism, stabilizing a megabase-scale inversion at an intermediate allele frequency in the wild.

PAX2 + Mesenchymal Origin of Gonadal Supporting Cells Is Conserved in Birds

During embryonic gonadal development, the supporting cell lineage is the first cell type to differentiate, giving rise to Sertoli cells in the testis and pre-granulosa cells in the ovary. These cells are thought to direct other gonadal cell lineages down the testis or ovarian pathways, including the germline. Recent research has shown that, in contrast to mouse, chicken gonadal supporting cells derive from a PAX2/OSR1/DMRT1/WNT4 positive mesenchymal cell population. These cells colonize the undifferentiated genital ridge during early gonadogenesis, around the time that germ cells migrate into the gonad. During the process of somatic gonadal sex differentiation, PAX2 expression is down-regulated in embryonic chicken gonads just prior to up-regulation of testis- and ovary-specific markers and prior to germ cell differentiation. Most research on avian gonadal development has focused on the chicken model, and related species from the Galloanserae clade. There is a lack of knowledge on gonadal sex differentiation in other avian lineages. Comparative analysis in birds is required to fully understand the mechanisms of avian sex determination and gonadal differentiation. Here we report the first comparative molecular characterization of gonadal supporting cell differentiation in birds from each of the three main clades, Galloanserae (chicken and quail), Neoaves (zebra finch) and Palaeognathe (emu). Our analysis reveals conservation of PAX2+ expression and a mesenchymal origin of supporting cells in each clade. Moreover, down-regulation of PAX2 expression precisely defines the onset of gonadal sex differentiation in each species. Altogether, these results indicate that gonadal morphogenesis is conserved among the major bird clades.

Geographic contact drives increased reproductive isolation in two cryptic Empidonax flycatchers

Geographic contact between sister lineages often occurs near the final stages of speciation, but its role in speciation's completion remains debated. Reproductive isolation may be essentially complete prior to secondary contact. Alternatively, costly interactions between partially reproductively isolated species - such as maladaptive hybridization or competition for resources - may select for divergence, increasing reproductive isolation and driving speciation toward completion. Here, we use coalescent demographic modelling and whole-genome data sets to show that a period of contact and elevated hybridization between sympatric eastern North American populations of two cryptic bird species preceded a major increase in reproductive isolation between these populations within the last 10,000 years. In contrast, substantial introgression continues to the present in a western contact zone where geographic overlap is much narrower and probably of more recent origin. In the sympatric eastern region where reproductive isolation has increased, it is not accompanied by character displacement in key morphometric traits, plumage coloration, or ecological traits. While the precise trait and underlying mechanism driving increased reproductive isolation remains unknown, we discuss several possibilities and outline avenues for future research. Overall, our results highlight how demographic models can reveal the geographic context in which reproductive isolation was completed, and demonstrate how contact can accelerate the final stages of speciation.

The Genomic Signature of Allopatric Speciation in a Songbird Is Shaped by Genome Architecture (Aves: Certhia americana)

The genomic signature of speciation with gene flow is often attributed to the strength of divergent selection and recombination rate in regions harboring targets for selection. In contrast, allopatric speciation provides a different geographic context and evolutionary scenario, whereby introgression is limited by isolation rather than selection against gene flow. Lacking shared divergent selection or selection against hybridization, we would predict the genomic signature of allopatric speciation would largely be shaped by genomic architecture-the nonrandom distribution of functional elements and chromosomal characteristics-through its role in affecting the processes of selection and drift. Here, we built and annotated a chromosome-scale genome assembly for a songbird (Passeriformes: Certhia americana). We show that the genomic signature of allopatric speciation between its two primary lineages is largely shaped by genomic architecture. Regionally, gene density and recombination rate variation explain a large proportion of variance in genomic diversity, differentiation, and divergence. We identified a heterogeneous landscape of selection and neutrality, with a large portion of the genome under the effects of indirect selection. We found higher proportions of small chromosomes under the effects of indirect selection, likely because they have relatively higher gene density. At the chromosome scale, differential genomic architecture of macro- and microchromosomes shapes the genomic signatures of speciation: chromosome size has: 1) a positive relationship with genetic differentiation, genetic divergence, rate of lineage sorting in the contact zone, and proportion neutral evolution and 2) a negative relationship with genetic diversity and recombination rate.

Autosomal, sex-linked and mitochondrial loci resolve evolutionary relationships among wrens in the genus Campylorhynchus

Although there is general consensus that sampling of multiple genetic loci is critical in accurate reconstruction of species trees, the exact numbers and the best types of molecular markers remain an open question. In particular, the phylogenetic utility of sex-linked loci is underexplored. Here, we sample all species and 70% of the named diversity of the New World wren genus Campylorhynchus using sequences from 23 loci, to evaluate the effects of linkage on efficiency in recovering a well-supported tree for the group. At a tree-wide level, we found that most loci supported fewer than half the possible clades and that sex-linked loci produced similar resolution to slower-coalescing autosomal markers, controlling for locus length. By contrast, we did find evidence that linkage affected the efficiency of recovery of individual relationships; as few as two sex-linked loci were necessary to resolve a selection of clades with long to medium subtending branches, whereas 4-6 autosomal loci were necessary to achieve comparable results. These results support an expanded role for sampling of the avian Z chromosome in phylogenetic studies, including target enrichment approaches. Our concatenated and species tree analyses represent significant improvements in our understanding of diversification in Campylorhynchus, and suggest a relatively complex scenario for its radiation across the Miocene/Pliocene boundary, with multiple invasions of South America.

Highly contiguous genomes improve the understanding of avian olfactory receptor repertoires

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

Neurogenomic insights into the behavioral and vocal development of the zebra finch

The zebra finch (Taeniopygia guttata) is a socially monogamous and colonial opportunistic breeder with pronounced sexual differences in singing and plumage coloration. Its natural history has led to it becoming a model species for research into sex differences in vocal communication, as well as behavioral, neural and genomic studies of imitative auditory learning. As scientists tap into the genetic and behavioral diversity of both wild and captive lineages, the zebra finch will continue to inform research into culture, learning, and social bonding, as well as adaptability to a changing climate.

Molecular changes associated with migratory departure from wintering areas in obligate songbird migrants

Day length regulates the development of spring migratory and subsequent reproductive phenotypes in avian migrants. This study used molecular approaches, and compared mRNA and proteome-wide expression in captive redheaded buntings that were photostimulated under long-day (LD) conditions for 4 days (early stimulated, LD-eS) or for ∼3 weeks until each bird had shown 4 successive nights of Zugunruhe (stimulated, LD-S); controls were maintained under short days. After ∼3 weeks of LD, photostimulated indices of the migratory preparedness (fattening, weight gain and Zugunruhe) were paralleled with upregulated expression of acc, dgat2 and apoa1 genes in the liver, and of cd36, fabp3 and cpt1 genes in the flight muscle, suggesting enhanced fatty acid (FA) synthesis and transport in the LD-S state. Concurrently, elevated expression of genes involved in the calcium ion signalling and transport (camk1 and atp2a2; camk2a in LD-eS), cellular stress (hspa8 and sod1, not nos2) and metabolic pathways (apoa1 and sirt1), but not of genes associated with migratory behaviour (adcyap1 and vps13a), were found in the mediobasal hypothalamus (MBH). Further, MBH-specific quantitative proteomics revealed that out of 503 annotated proteins, 28 were differentially expressed (LD-eS versus LD-S: 21 up-regulated and 7 down-regulated) and they enriched five physiological pathways that are associated with FA transport and metabolism. These first comprehensive results on gene and protein expression suggest that changes in molecular correlates of FA transport and metabolism may aid the decision for migratory departure from wintering areas in obligate songbird migrants.

The Gene Expression Profile of the Song Control Nucleus HVC Shows Sex Specificity, Hormone Responsiveness, and Species Specificity Among Songbirds

Singing occurs in songbirds of both sexes, but some species show typical degrees of sex-specific performance. We studied the transcriptional sex differences in the HVC, a brain nucleus critical for song pattern generation, of the forest weaver (Ploceus bicolor), the blue-capped cordon-bleu (Uraeginthus cyanocephalus), and the canary (Serinus canaria), which are species that show low, medium, and high levels of sex-specific singing, respectively. We observed persistent sex differences in gene expression levels regardless of the species-specific sexual singing phenotypes. We further studied the HVC transcriptomes of defined phenotypes of canary, known for its testosterone-sensitive seasonal singing. By studying both sexes of canaries during both breeding and non-breeding seasons, non-breeding canaries treated with testosterone, and spontaneously singing females, we found that the circulating androgen levels and sex were the predominant variables associated with the variations in the HVC transcriptomes. The comparison of natural singing with testosterone-induced singing in canaries of the same sex revealed considerable differences in the HVC transcriptomes. Strong transcriptional changes in the HVC were detected during the transition from non-singing to singing in canaries of both sexes. Although the sex-specific genes of singing females shared little resemblance with those of males, our analysis showed potential functional convergences. Thus, male and female songbirds achieve comparable singing behaviours with sex-specific transcriptomes.

Karyotype Evolution and Genomic Organization of Repetitive DNAs in the Saffron Finch, Sicalis flaveola (Passeriformes, Aves)

The Saffron finch (Sicalis flaveola), a semi-domestic species, is tolerant of human proximity and nesting in roof spaces. Considering the importance of cytogenomic approaches in revealing different aspects of genomic organization and evolution, we provide detailed cytogenetic data for S. flaveola, including the standard Giemsa karyotype, C- and G-banding, repetitive DNA mapping, and bacterial artificial chromosome (BAC) FISH. We also compared our results with the sister groups, Passeriformes and Psittaciformes, bringing new insights into the chromosome and genome evolution of birds. The results revealed contrasting rates of intrachromosomal changes, highlighting the role of SSR (simple short repetition probes) accumulation in the karyotype reorganization. The SSRs showed scattered hybridization, but brighter signals were observed in the microchromosomes and the short arms of Z chromosome in S. flaveola. BACs probes showed conservation of ancestral syntenies of macrochromosomes (except GGA1), as well as the tested microchromosomes. The comparison of our results with previous studies indicates that the great biological diversity observed in Passeriformes was not likely accompanied by interchromosomal changes. In addition, although repetitive sequences often act as hotspots of genome rearrangements, Passeriformes species showed a higher number of signals when compared with the sister group Psittaciformes, indicating that these sequences were not involved in the extensive karyotype reorganization seen in the latter.

Adaptive Radiation Genomics of Two Ecologically Divergent Hawai'ian Honeycreepers: The 'akiapōlā'au and the Hawai'i 'amakihi

The Hawai'ian honeycreepers (drepanids) are a classic example of adaptive radiation: they adapted to a variety of novel dietary niches, evolving a wide range of bill morphologies. Here we investigated genomic diversity, demographic history, and genes involved in bill morphology phenotypes in 2 honeycreepers: the 'akiapōlā'au (Hemignathus wilsoni) and the Hawai'i 'amakihi (Chlorodrepanis virens). The 'akiapōlā'au is an endangered island endemic, filling the "woodpecker" niche by using a unique bill morphology, while the Hawai'i 'amakihi is a dietary generalist common on the islands of Hawai'i and Maui. We de novo sequenced the 'akiapōlā'au genome and compared it to the previously sequenced 'amakihi genome. The 'akiapōlā'au is far less heterozygous and has a smaller effective population size than the 'amakihi, which matches expectations due to its smaller census population and restricted ecological niche. Our investigation revealed genomic islands of divergence, which may be involved in the honeycreeper radiation. Within these islands of divergence, we identified candidate genes (including DLK1, FOXB1, KIF6, MAML3, PHF20, RBP1, and TIMM17A) that may play a role in honeycreeper adaptations. The gene DLK1, previously shown to influence Darwin's finch bill size, may be related to honeycreeper bill morphology evolution, while the functions of the other candidates remain unknown.

As above, so below: Whole transcriptome profiling demonstrates strong molecular similarities between avian dorsal and ventral pallial subdivisions

Over the last two decades, beginning with the Avian Brain Nomenclature Forum in 2000, major revisions have been made to our understanding of the organization and nomenclature of the avian brain. However, there are still unresolved questions on avian pallial organization, particularly whether the cells above the vestigial ventricle represent distinct populations to those below it or similar populations. To test these two hypotheses, we profiled the transcriptomes of the major avian pallial subdivisions dorsal and ventral to the vestigial ventricle boundary using RNA sequencing and a new zebra finch genome assembly containing about 22,000 annotated, complete genes. We found that the transcriptomes of neural populations above and below the ventricle were remarkably similar. Each subdivision in dorsal pallium (Wulst) had a corresponding molecular counterpart in the ventral pallium (dorsal ventricular ridge). In turn, each corresponding subdivision exhibited shared gene co‐expression modules that contained gene sets enriched in functional specializations, such as anatomical structure development, synaptic transmission, signaling, and neurogenesis. These findings are more in line with the continuum hypothesis of avian brain subdivision organization above and below the vestigial ventricle space, with the pallium as a whole consisting of four major cell populations (intercalated pallium, mesopallium, hyper‐nidopallium, and arcopallium) instead of seven (hyperpallium apicale, interstitial hyperpallium apicale, intercalated hyperpallium, hyperpallium densocellare, mesopallium, nidopallium, and arcopallium). We suggest adopting a more streamlined hierarchical naming system that reflects the robust similarities in gene expression, neural connectivity motifs, and function. These findings have important implications for our understanding of overall vertebrate brain evolution.

An island-hopping bird reveals how founder events shape genome-wide divergence

When populations colonize new areas, both strong selection and strong drift can be experienced due to novel environments and small founding populations, respectively. Empirical studies have predominantly focused on the phenotype when assessing the role of selection, and limited neutral-loci when assessing founder-induced loss of diversity. Consequently, the extent to which processes interact to influence evolutionary trajectories is difficult to assess. Genomic-level approaches provide the opportunity to simultaneously consider these processes. Here, we examine the roles of selection and drift in shaping genomic diversity and divergence in historically documented sequential island colonizations by the silvereye (Zosterops lateralis). We provide the first empirical demonstration of the rapid appearance of highly diverged genomic regions following population founding, the position of which are highly idiosyncratic. As these regions rarely contained loci putatively under selection, it is most likely that these differences arise via the stochastic nature of the founding process. However, selection is required to explain rapid evolution of larger body size in insular silvereyes. Reconciling our genomic data with these phenotypic patterns suggests there may be many genomic routes to the island phenotype, which vary across populations. Finally, we show that accelerated divergence associated with multiple founding steps is the product of genome-wide rather than localized differences, and that diversity erodes due to loss of rare alleles. However, even multiple founder events do not result in divergence and diversity levels seen in evolutionary older subspecies, and therefore do not provide a shortcut to speciation as proposed by founder-effect speciation models.

Chromosome-Level Genome Assembly of the Common Chaffinch (Aves: Fringilla coelebs): A Valuable Resource for Evolutionary Biology

The common chaffinch, Fringilla coelebs, is one of the most common, widespread, and well-studied passerines in Europe, with a broad distribution encompassing Western Europe and parts of Asia, North Africa, and the Macaronesian archipelagos. We present a high-quality genome assembly of the common chaffinch generated using Illumina shotgun sequencing in combination with Chicago and Hi-C libraries. The final genome is a 994.87-Mb chromosome-level assembly, with 98% of the sequence data located in chromosome scaffolds and a N50 statistic of 69.73 Mb. Our genome assembly shows high completeness, with a complete BUSCO score of 93.9% using the avian data set. Around 7.8% of the genome contains interspersed repetitive elements. The structural annotation yielded 17,703 genes, 86.5% of which have a functional annotation, including 7,827 complete universal single-copy orthologs out of 8,338 genes represented in the BUSCO avian data set. This new annotated genome assembly will be a valuable resource as a reference for comparative and population genomic analyses of passerine, avian, and vertebrate evolution.

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

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

Comparative Analysis of Annotation Pipelines Using the First Japanese White-Eye (Zosterops japonicus) Genome

Introduced into Hawaii in the early 1900s, the Japanese white-eye or warbling white-eye (Zosterops japonicus) is now the most abundant land bird in the archipelago. Here, we present the first Z. japonicus genome, sequenced from an individual in its invasive range. This genome provides an important resource for future studies in invasion genomics. We annotated the genome using two workflows-standalone AUGUSTUS and BRAKER2. We found that AUGUSTUS was more conservative with gene predictions when compared with BRAKER2. The final number of annotated gene models was similar between the two workflows, but standalone AUGUSTUS had over 70% of gene predictions with Blast2GO annotations versus under 30% using BRAKER2. Additionally, we tested whether using RNA-seq data from 47 samples had a significant impact on annotation quality when compared with data from a single sample, as generating RNA-seq data for genome annotation can be expensive and requires well preserved tissue. We found that more data did not significantly change the number of annotated genes using AUGUSTUS but using BRAKER2 the number increased substantially. The results presented here will aid researchers in annotating draft genomes of nonmodel species as well as those studying invasion success.

Snake Recombination Landscapes Are Concentrated in Functional Regions despite PRDM9

Meiotic recombination in vertebrates is concentrated in hotspots throughout the genome. The location and stability of hotspots have been linked to the presence or absence of PRDM9, leading to two primary models for hotspot evolution derived from mammals and birds. Species with PRDM9-directed recombination have rapid turnover of hotspots concentrated in intergenic regions (i.e., mammals), whereas hotspots in species lacking PRDM9 are concentrated in functional regions and have greater stability over time (i.e., birds). Snakes possess PRDM9, yet virtually nothing is known about snake recombination. Here, we examine the recombination landscape and test hypotheses about the roles of PRDM9 in rattlesnakes. We find substantial variation in recombination rate within and among snake chromosomes, and positive correlations between recombination rate and gene density, GC content, and genetic diversity. Like mammals, snakes appear to have a functional and active PRDM9, but rather than being directed away from genes, snake hotspots are concentrated in promoters and functional regions-a pattern previously associated only with species that lack a functional PRDM9. Snakes therefore provide a unique example of recombination landscapes in which PRDM9 is functional, yet recombination hotspots are associated with functional genic regions-a combination of features that defy existing paradigms for recombination landscapes in vertebrates. Our findings also provide evidence that high recombination rates are a shared feature of vertebrate microchromosomes. Our results challenge previous assumptions about the adaptive role of PRDM9 and highlight the diversity of recombination landscape features among vertebrate lineages.