The genome of a songbird

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.

Prenatal testosterone triggers long-term behavioral changes in male zebra finches: unravelling the neurogenomic mechanisms

BACKGROUND

Maternal hormones, like testosterone, can strongly influence developing offspring, even generating long-term organizational effects on adult behavior; yet, the mechanisms facilitating these effects are still unclear. Here, we experimentally elevated prenatal testosterone in the eggs of zebra finches (Taeniopygia guttata) and measured male aggression in adulthood along with patterns of neural gene expression (RNA-seq) and DNA methylation (MethylC-Seq) in two socially relevant brain regions (hypothalamus and nucleus taenia of the amygdala). We used enrichment analyses and protein-protein interaction networks to find candidate processes and hub genes potentially affected by the treatment. We additionally identified differentially expressed genes that contained differentially methylated regions.

RESULTS

We found that males from testosterone-injected eggs displayed more aggressive behaviors compared to males from control eggs. Hundreds of genes were differentially expressed, particularly in the hypothalamus, including potential aggression-related hub genes (e.g., brain derived neurotrophic factor). There were also enriched processes with well-established links to aggressive phenotypes (e.g., somatostatin and glutamate signaling). Furthermore, several highly connected genes identified in protein-protein interaction networks also showed differential methylation, including adenylate cyclase 2 and proprotein convertase 2.

CONCLUSIONS

These results highlight genes and processes that may play an important role in mediating the effects of prenatal testosterone on long-term phenotypic outcomes, thereby providing insights into the molecular mechanisms that facilitate hormone-mediated maternal effects.

Analysis of multiple chromosomal rearrangements in the genome of Willisornis vidua using BAC-FISH and chromosome painting on a supposed conserved karyotype

BACKGROUND

Thamnophilidae birds are the result of a monophyletic radiation of insectivorous Passeriformes. They are a diverse group of 225 species and 45 genera and occur in lowlands and lower montane forests of Neotropics. Despite the large degree of diversity seen in this family, just four species of Thamnophilidae have been karyotyped with a diploid number ranging from 76 to 82 chromosomes. The karyotypic relationships within and between Thamnophilidae and another Passeriformes therefore remain poorly understood. Recent studies have identified the occurrence of intrachromosomal rearrangements in Passeriformes using in silico data and molecular cytogenetic tools. These results demonstrate that intrachromosomal rearrangements are more common in birds than previously thought and are likely to contribute to speciation events. With this in mind, we investigate the apparently conserved karyotype of Willisornis vidua, the Xingu Scale-backed Antbird, using a combination of molecular cytogenetic techniques including chromosome painting with probes derived from Gallus gallus (chicken) and Burhinus oedicnemus (stone curlew), combined with Bacterial Artificial Chromosome (BAC) probes derived from the same species. The goal was to investigate the occurrence of rearrangements in an apparently conserved karyotype in order to understand the evolutionary history and taxonomy of this species. In total, 78 BAC probes from the Gallus gallus and Taeniopygia guttata (the Zebra Finch) BAC libraries were tested, of which 40 were derived from Gallus gallus macrochromosomes 1-8, and 38 from microchromosomes 9-28.

RESULTS

The karyotype is similar to typical Passeriformes karyotypes, with a diploid number of 2n = 80. Our chromosome painting results show that most of the Gallus gallus chromosomes are conserved, except GGA-1, 2 and 4, with some rearrangements identified among macro- and microchromosomes. BAC mapping revealed many intrachromosomal rearrangements, mainly inversions, when comparing Willisornis vidua karyotype with Gallus gallus, and corroborates the fissions revealed by chromosome painting.

CONCLUSIONS

Willisornis vidua presents multiple chromosomal rearrangements despite having a supposed conservative karyotype, demonstrating that our approach using a combination of FISH tools provides a higher resolution than previously obtained by chromosome painting alone. We also show that populations of Willisornis vidua appear conserved from a cytogenetic perspective, despite significant phylogeographic structure.

Estrogen and sex-dependent loss of the vocal learning system in female zebra finches

Sex hormones alter the organization of the brain during early development and coordinate various behaviors throughout life. In zebra finches, song learning is limited to males, with the associated song learning brain pathways only maturing in males and atrophying in females. While this atrophy can be prevented by treating females with exogenous estrogen during early post-hatch development, the requirement of estrogen during normal male song system development is uncertain. For the first time in songbirds, we administered exemestane, a potent third generation estrogen synthesis inhibitor, from the day of hatching until adulthood in order to reassess the role of estrogen in song circuit development. We examined the behavior, brain anatomy, and transcriptomes of individual song nuclei in these pharmacologically manipulated animals. We found that males with long-term exemestane treatment had diminished male-specific plumage and impaired song learning, but minimal effect on song nuclei sizes and their specialized transcriptome. Consistent with prior findings, females with long-term estrogen treatment retained a functional song system with song nuclei that had specialized gene expression similar, but not identical to males. We also observed that different song nuclei responded to estrogen manipulation differently, with Area X in the striatum being the most altered by estrogen modulation. These findings support the hypothesis that song learning is an ancestral trait in both sexes that was subsequently suppressed in females of some species and that estrogen has come to play a critical role in modulating this suppression as well as refinement of song learning.

Genome Size Reduction and Transposon Activity Impact tRNA Gene Diversity While Ensuring Translational Stability in Birds

As a highly diverse vertebrate class, bird species have adapted to various ecological systems. How this phenotypic diversity can be explained genetically is intensively debated and is likely grounded in differences in the genome content. Larger and more complex genomes could allow for greater genetic regulation that results in more phenotypic variety. Surprisingly, avian genomes are much smaller compared to other vertebrates but contain as many protein-coding genes as other vertebrates. This supports the notion that the phenotypic diversity is largely determined by selection on non-coding gene sequences. Transfer RNAs (tRNAs) represent a group of non-coding genes. However, the characteristics of tRNA genes across bird genomes have remained largely unexplored. Here, we exhaustively investigated the evolution and functional consequences of these crucial translational regulators within bird species and across vertebrates. Our dense sampling of 55 avian genomes representing each bird order revealed an average of 169 tRNA genes with at least 31% being actively used. Unlike other vertebrates, avian tRNA genes are reduced in number and complexity but are still in line with vertebrate wobble pairing strategies and mutation-driven codon usage. Our detailed phylogenetic analyses further uncovered that new tRNA genes can emerge through multiplication by transposable elements. Together, this study provides the first comprehensive avian and cross-vertebrate tRNA gene analyses and demonstrates that tRNA gene evolution is flexible albeit constrained within functional boundaries of general mechanisms in protein translation.

Genomic variation, population history and within-archipelago adaptation between island bird populations

Oceanic island archipelagos provide excellent models to understand evolutionary processes. Colonization events and gene flow can interact with selection to shape genetic variation at different spatial scales. Landscape-scale variation in biotic and abiotic factors may drive fine-scale selection within islands, while long-term evolutionary processes may drive divergence between distantly related populations. Here, we examine patterns of population history and selection between recently diverged populations of the Berthelot's pipit (Anthus berthelotii), a passerine endemic to three North Atlantic archipelagos. First, we use demographic trees and f3 statistics to show that genome-wide divergence across the species range is largely shaped by colonization and bottlenecks, with evidence of very weak gene flow between populations. Then, using a genome scan approach, we identify signatures of divergent selection within archipelagos at single nucleotide polymorphisms (SNPs) in genes potentially associated with craniofacial development and DNA repair. We did not detect within-archipelago selection at the same SNPs as were detected previously at broader spatial scales between archipelagos, but did identify signatures of selection at loci associated with similar biological functions. These findings suggest that similar ecological factors may repeatedly drive selection between recently separated populations, as well as at broad spatial scales across varied landscapes.

The evolution of a tropical biodiversity hotspot

The tropics are the source of most biodiversity yet inadequate sampling obscures answers to fundamental questions about how this diversity evolves. We leveraged samples assembled over decades of fieldwork to study diversification of the largest tropical bird radiation, the suboscine passerines. Our phylogeny, estimated using data from 2389 genomic regions in 1940 individuals of 1283 species, reveals that peak suboscine species diversity in the Neotropics is not associated with high recent speciation rates but rather with the gradual accumulation of species over time. Paradoxically, the highest speciation rates are in lineages from regions with low species diversity, which are generally cold, dry, unstable environments. Our results reveal a model in which species are forming faster in environmental extremes but have accumulated in moderate environments to form tropical biodiversity hotspots.

K. VP, Vijay N, Sharma VK. Genomic Insights Into the Molecular Basis of Sexual Selection in Birds

Sexual selection is a well-known biological process, yet the genomic basis and patterns of sexual selection are not fully understood. The extravagant ornamental plumage of peacock (Pavo cristatus) was instrumental in shaping Charles Darwin's theory of sexual selection and is considered to be an honest signal of its immunocompetence. Here, we used the recently generated draft genome sequence of peafowl (Pavo cristatus) and carried out a comparative analysis across 11 bird genomes that encompass a range of sexual selection and also had high-quality genomic and phenotypic data publically available to study the genomic basis of sexual selection. We found that varying degree of purifying selection was the predominant mechanism of action for sexual selection at the genome-wide scale and observed that sexual selection mostly influences genes regulating gene expression and protein processing. Specifically, the genome-wide phylogenetically corrected regression analysis supported the continuous or ongoing model of sexual selection. Genes involved in nucleic acid binding and gene expression regulation, including a specific regulator of sex-determination known as TRA2A to be under positive selection in the species with high post-copulatory sexual selection manifested as high sperm competition. We also detected specific feather-related and immune-related gene-pairs evolving under similar selection pressures across the 11 species, including peacock (Pavo cristatus), which is consistent with the Hamilton-Zuk hypothesis. The comparative genomics analysis of 11 avian taxa has provided new insights on the molecular underpinnings of sexual selection and identifies specific genomic regions for future in-depth analysis.

Island songbirds as windows into evolution in small populations

Due to their limited ranges and inherent isolation, island species have long been recognized as crucial systems for tackling a range of evolutionary questions, including in the early study of speciation.^1,2 Such species have been less studied in the understanding of the evolutionary forces driving DNA sequence evolution. Island species usually have lower census population sizes (N) than continental species and, supposedly, lower effective population sizes (Ne). Given that both the rates of change caused by genetic drift and by selection are dependent upon Ne, island species are theoretically expected to exhibit (1) lower genetic diversity, (2) less effective natural selection against slightly deleterious mutations,^3,4 and (3) a lower rate of adaptive evolution.^5-8 Here, we have used a large set of newly sequenced and published whole-genome sequences of Passerida species (14 insular and 11 continental) to test these predictions. We confirm that island species exhibit lower census size and Ne, supporting the hypothesis that the smaller area available on islands constrains the upper bound of Ne. In the insular species, we find lower nucleotide diversity in coding regions, higher ratios of non-synonymous to synonymous polymorphisms, and lower adaptive substitution rates. Our results provide robust evidence that the lower Ne experienced by island species has affected both the ability of natural selection to efficiently remove weakly deleterious mutations and also the adaptive potential of island species, therefore providing considerable empirical support for the nearly neutral theory. We discuss the implications for both evolutionary and conservation biology.

Avian Population Genomics Taking Off: Latest Findings and Future Prospects

Birds are one of the most recognizable and diverse groups of organisms on earth. This group has played an important role in many fields, including the development of methods in behavioral ecology and evolutionary theory. The use of population genomics took off following the advent of high-throughput sequencing in various taxa. Several features of avian genomes make them particularly amenable for work in this field, including their nucleated red blood cells permitting easy DNA extraction and small, compact genomes. We review the latest findings in the population genomics of birds here, emphasizing questions related to behavior, ecology, evolution, and conservation. Additionally, we include insights in trait mapping and the ability to obtain accurate estimates of important summary statistics for conservation (e.g., genetic diversity and inbreeding). We highlight roadblocks that will need to be overcome in order to advance work on the population genomics of birds and prospects for future work. Roadblocks include the assembly of more contiguous reference genomes using long-reads and optical mapping. Prospects include the integration of population genomics with additional fields (e.g., landscape genetics, phylogeography, and genomic mapping) along with studies beyond genetic variants (e.g., epigenetics).

Identifying the causes and consequences of assembly gaps using a multiplatform genome assembly of a bird-of-paradise

Genome assemblies are currently being produced at an impressive rate by consortia and individual laboratories. The low costs and increasing efficiency of sequencing technologies now enable assembling genomes at unprecedented quality and contiguity. However, the difficulty in assembling repeat-rich and GC-rich regions (genomic "dark matter") limits insights into the evolution of genome structure and regulatory networks. Here, we compare the efficiency of currently available sequencing technologies (short/linked/long reads and proximity ligation maps) and combinations thereof in assembling genomic dark matter. By adopting different de novo assembly strategies, we compare individual draft assemblies to a curated multiplatform reference assembly and identify the genomic features that cause gaps within each assembly. We show that a multiplatform assembly implementing long-read, linked-read and proximity sequencing technologies performs best at recovering transposable elements, multicopy MHC genes, GC-rich microchromosomes and the repeat-rich W chromosome. Telomere-to-telomere assemblies are not a reality yet for most organisms, but by leveraging technology choice it is now possible to minimize genome assembly gaps for downstream analysis. We provide a roadmap to tailor sequencing projects for optimized completeness of both the coding and noncoding parts of nonmodel genomes.

The potential role of Arhgef33 RhoGEF in foveal development in the zebra finch retina

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.