Dense sampling of bird diversity increases power of comparative genomics

Remarkably High Repeat Content in the Genomes of Sparrows: The Importance of Genome Assembly Completeness for Transposable Element Discovery

Transposable elements (TE) play critical roles in shaping genome evolution. Highly repetitive TE sequences are also a major source of assembly gaps making it difficult to fully understand the impact of these elements on host genomes. The increased capacity of long-read sequencing technologies to span highly repetitive regions promises to provide new insights into patterns of TE activity across diverse taxa. Here we report the generation of highly contiguous reference genomes using PacBio long-read and Omni-C technologies for three species of Passerellidae sparrow. We compared these assemblies to three chromosome-level sparrow assemblies and nine other sparrow assemblies generated using a variety of short- and long-read technologies. All long-read based assemblies were longer (range: 1.12 to 1.41 Gb) than short-read assemblies (0.91 to 1.08 Gb) and assembly length was strongly correlated with the amount of repeat content. Repeat content for Bell's sparrow (31.2% of genome) was the highest level ever reported within the order Passeriformes, which comprises over half of avian diversity. The highest levels of repeat content (79.2% to 93.7%) were found on the W chromosome relative to other regions of the genome. Finally, we show that proliferation of different TE classes varied even among species with similar levels of repeat content. These patterns support a dynamic model of TE expansion and contraction even in a clade where TEs were once thought to be fairly depauperate and static. Our work highlights how the resolution of difficult-to-assemble regions of the genome with new sequencing technologies promises to transform our understanding of avian genome evolution.

Improving the performance of supervised deep learning for regulatory genomics using phylogenetic augmentation

MOTIVATION

Supervised deep learning is used to model the complex relationship between genomic sequence and regulatory function. Understanding how these models make predictions can provide biological insight into regulatory functions. Given the complexity of the sequence to regulatory function mapping (the cis-regulatory code), it has been suggested that the genome contains insufficient sequence variation to train models with suitable complexity. Data augmentation is a widely used approach to increase the data variation available for model training, however current data augmentation methods for genomic sequence data are limited.

RESULTS

Inspired by the success of comparative genomics, we show that augmenting genomic sequences with evolutionarily related sequences from other species, which we term phylogenetic augmentation, improves the performance of deep learning models trained on regulatory genomic sequences to predict high-throughput functional assay measurements. Additionally, we show that phylogenetic augmentation can rescue model performance when the training set is down-sampled and permits deep learning on a real-world small dataset, demonstrating that this approach improves data efficiency. Overall, this data augmentation method represents a solution for improving model performance that is applicable to many supervised deep-learning problems in genomics.

AVAILABILITY AND IMPLEMENTATION

The open-source GitHub repository agduncan94/phylogenetic_augmentation_paper includes the code for rerunning the analyses here and recreating the figures.

Detection of evolutionary conserved and accelerated genomic regions related to adaptation to thermal niches in Anolis lizards

Understanding the genetic basis for adapting to thermal environments is important due to serious effects of global warming on ectothermic species. Various genes associated with thermal adaptation in lizards have been identified mainly focusing on changes in gene expression or the detection of positively selected genes using coding regions. Only a few comprehensive genome-wide analyses have included noncoding regions. This study aimed to identify evolutionarily conserved and accelerated genomic regions using whole genomes of eight Anolis lizard species that have repeatedly adapted to similar thermal environments in multiple lineages. Evolutionarily conserved genomic regions were extracted as regions with overall sequence conservation (regions with fewer base substitutions) across all lineages compared with the neutral model. Genomic regions that underwent accelerated evolution in the lineage of interest were identified as those with more base substitutions in the target branch than in the entire background branch. Conserved elements across all branches were relatively abundant in "intergenic" genomic regions among noncoding regions. Accelerated regions (ARs) of each lineage contained a significantly greater proportion of noncoding RNA genes than the entire multiple alignment. Common genes containing ARs within 5 kb of their vicinity in lineages with similar thermal habitats were identified. Many genes associated with circadian rhythms and behavior were found in hot-open and cool-shaded habitat lineages. These genes might play a role in contributing to thermal adaptation and assist future studies examining the function of genes involved in thermal adaptation via genome editing.

Landscape Heterogeneity Explains the Genetic Differentiation of a Forest Bird across the Sino-Himalayan Mountains

Mountains are the world's most important centers of biodiversity. The Sino-Himalayan Mountains are global biodiversity hotspot due to their extremely high species richness and endemicity. Ample research investigated the impact of the Qinghai-Tibet Plateau uplift and Quaternary glaciations in driving species diversification in plants and animals across the Sino-Himalayan Mountains. However, little is known about the role of landscape heterogeneity and other environmental features in driving diversification in this region. We utilized whole genomes and phenotypic data in combination with landscape genetic approaches to investigate population structure, demography, and genetic diversity in a forest songbird species native to the Sino-Himalayan Mountains, the red-billed leiothrix (Leiothrix lutea). We identified 5 phylogeographic clades, including 1 in the East of China, 1 in Yunnan, and 3 in Tibet, roughly consistent with differences in song and plumage coloration but incongruent with traditional subspecies boundaries. Isolation-by-resistance model best explained population differentiation within L. lutea, with extensive secondary contact after allopatric isolation leading to admixture among clades. Ecological niche modeling indicated relative stability in the extent of suitable distribution areas of the species across Quaternary glacial cycles. Our results underscore the importance of mountains in the diversification of this species, given that most of the distinct genetic clades are concentrated in a relatively small area in the Sino-Himalayan Mountain region, while a single shallow clade populates vast lower-lying areas to the east. This study highlights the crucial role of landscape heterogeneity in promoting differentiation and provides a deep genomic perspective on the mechanisms through which diversity hotspots form.

Brief Communication: The Predictable Network Topology of Evolutionary Genomic Constraint

Large-scale comparative genomics studies offer valuable resources for understanding both functional and evolutionary rate constraints. It is suggested that constraint aligns with the topology of genomic networks, increasing toward the center, with intermediate nodes combining relaxed constraint with higher contributions to the phenotype due to pleiotropy. However, this pattern has yet to be demonstrated in vertebrates. This study shows that constraint intensifies toward the network's center in placental mammals. Genes with rate changes associated with emergence of hibernation cluster mostly toward intermediate positions, with higher constraint in faster-evolving genes, which is indicative of a "sweet spot" for adaptation. If this trend holds universally, network node metrics could predict high-constraint regions even in clades lacking empirical constraint data.

Genetics and Evolution of Bird Migration

Bird migration has long been a subject of fascination for humankind and is a behavior that is both intricate and multifaceted. In recent years, advances in technology, particularly in the fields of genomics and animal tracking, have enabled significant progress in our understanding of this phenomenon. In this review, we provide an overview of the latest advancements in the genetics of bird migration, with a particular focus on genomics, and examine various factors that contribute to the evolution of this behavior, including climate change. Integration of research from the fields of genomics, ecology, and evolution can enhance our comprehension of the complex mechanisms involved in bird migration and inform conservation efforts in a rapidly changing world.

Adaptive evolution and loss of a putative magnetoreceptor in passerines

Migratory birds possess remarkable accuracy in orientation and navigation, which involves various compass systems including the magnetic compass. Identifying the primary magnetosensor remains a fundamental open question. Cryptochromes (Cry) have been shown to be magnetically sensitive, and Cry4a from a migratory songbird seems to show enhanced magnetic sensitivity in vitro compared to Cry4a from resident species. We investigate Cry and their potential involvement in magnetoreception in a phylogenetic framework, integrating molecular evolutionary analyses with protein dynamics modelling. Our analysis is based on 363 bird genomes and identifies different selection regimes in passerines. We show that Cry4a is characterized by strong positive selection and high variability, typical characteristics of sensor proteins. We identify key sites that are likely to have facilitated the evolution of an optimized sensory protein for night-time orientation in songbirds. Additionally, we show that Cry4 was lost in hummingbirds, parrots and Tyranni (Suboscines), and thus identified a gene deletion, which might facilitate testing the function of Cry4a in birds. In contrast, the other avian Cry (Cry1 and Cry2) were highly conserved across all species, indicating basal, non-sensory functions. Our results support a specialization or functional differentiation of Cry4 in songbirds which could be magnetosensation.

Ancestral genome reconstruction enhances transposable element annotation by identifying degenerate integrants

Growing evidence indicates that transposable elements (TEs) play important roles in evolution by providing genomes with coding and non-coding sequences. Identification of TE-derived functional elements, however, has relied on TE annotations in individual species, which limits its scope to relatively intact TE sequences. Here, we report a novel approach to uncover previously unannotated degenerate TEs (degTEs) by probing multiple ancestral genomes reconstructed from hundreds of species. We applied this method to the human genome and achieved a 10.8% increase in coverage over the most recent annotation. Further, we discovered that degTEs contribute to various cis-regulatory elements and transcription factor binding sites, including those of a known TE-controlling family, the KRAB zinc-finger proteins. We also report unannotated chimeric transcripts between degTEs and human genes expressed in embryos. This study provides a novel methodology and a freely available resource that will facilitate the investigation of TE co-option events on a full scale.

Past volcanic activity predisposes an endemic threatened seabird to negative anthropogenic impacts

Humans are regularly cited as the main driver of current biodiversity extinction, but the impact of historic volcanic activity is often overlooked. Pre-human evidence of wildlife abundance and diversity are essential for disentangling anthropogenic impacts from natural events. Réunion Island, with its intense and well-documented volcanic activity, endemic biodiversity, long history of isolation and recent human colonization, provides an opportunity to disentangle these processes. We track past demographic changes of a critically endangered seabird, the Mascarene petrel Pseudobulweria aterrima, using genome-wide SNPs. Coalescent modeling suggested that a large ancestral population underwent a substantial population decline in two distinct phases, ca. 125,000 and 37,000 years ago, coinciding with periods of major eruptions of Piton des Neiges. Subsequently, the ancestral population was fragmented into the two known colonies, ca. 1500 years ago, following eruptions of Piton de la Fournaise. In the last century, both colonies declined significantly due to anthropogenic activities, and although the species was initially considered extinct, it was rediscovered in the 1970s. Our findings suggest that the current conservation status of wildlife on volcanic islands should be firstly assessed as a legacy of historic volcanic activity, and thereafter by the increasing anthropogenic impacts, which may ultimately drive species towards extinction.

A chromosome-level genome assembly of a free-living white-crowned sparrow (Zonotrichia leucophrys gambelii)

The white-crowned sparrow, Zonotrichia leucophrys, is a passerine bird with a wide distribution and it is extensively adapted to environmental changes. It has historically acted as a model species in studies on avian ecology, physiology and behaviour. Here, we present a high-quality chromosome-level genome of Zonotrichia leucophrys using PacBio and OmniC sequencing data. Gene models were constructed by combining RNA-seq and Iso-seq data from liver, hypothalamus, and ovary. In total a 1,123,996,003 bp genome was generated, including 31 chromosomes assembled in complete scaffolds along with other, unplaced scaffolds. This high-quality genome assembly offers an important genomic resource for the research community using the white-crowned sparrow as a model for understanding avian genome biology and development, and provides a genomic basis for future studies, both fundamental and applied.

Unexpected divergence in magnetoreceptor MagR from robin and pigeon linked to two sequence variations

Birds exhibit extraordinary mobility and remarkable navigational skills, obtaining guidance cues from the Earth's magnetic field for orientation and long-distance movement. Bird species also show tremendous diversity in navigation strategies, with considerable differences even within the same taxa and among individuals from the same population. The highly conserved iron and iron-sulfur cluster binding magnetoreceptor (MagR) protein is suggested to enable animals, including birds, to detect the geomagnetic field and navigate accordingly. Notably, MagR is also implicated in other functions, such as electron transfer and biogenesis of iron-sulfur clusters, raising the question of whether variability exists in its biochemical and biophysical features among species, particularly birds. In the current study, we conducted a comparative analysis of MagR from two different bird species, including the migratory European robin and the homing pigeon. Sequence alignment revealed an extremely high degree of similarity between the MagRs of these species, with only three sequence variations. Nevertheless, two of these variations underpinned significant differences in metal binding capacity, oligomeric state, and magnetic properties. These findings offer compelling evidence for the marked differences in MagR between the two avian species, potentially explaining how a highly conserved protein can mediate such diverse functions.

Development-Associated Genes of the Epidermal Differentiation Complex (EDC)

The epidermal differentiation complex (EDC) is a cluster of genes that encode protein components of the outermost layers of the epidermis in mammals, reptiles and birds. The development of the stratified epidermis from a single-layered ectoderm involves an embryo-specific superficial cell layer, the periderm. An additional layer, the subperiderm, develops in crocodilians and over scutate scales of birds. Here, we review the expression of EDC genes during embryonic development. Several EDC genes are expressed predominantly or exclusively in embryo-specific cell layers, whereas others are confined to the epidermal layers that are maintained in postnatal skin. The S100 fused-type proteins scaffoldin and trichohyalin are expressed in the avian and mammalian periderm, respectively. Scaffoldin forms the so-called periderm granules, which are histological markers of the periderm in birds. Epidermal differentiation cysteine-rich protein (EDCRP) and epidermal differentiation protein containing DPCC motifs (EDDM) are expressed in the avian subperiderm where they are supposed to undergo cross-linking via disulfide bonds. Furthermore, a histidine-rich epidermal differentiation protein and feather-type corneous beta-proteins, also known as beta-keratins, are expressed in the subperiderm. The accumulating evidence for roles of EDC genes in the development of the epidermis has implications on the evolutionary diversification of the skin in amniotes.

Comparison of the DNBSEQ platform and Illumina HiSeq 2000 for bacterial genome assembly

The Illumina HiSeq platform has been a commonly used option for bacterial genome sequencing. Now the BGI DNA nanoball (DNB) nanoarrays platform may provide an alternative platform for sequencing of bacterial genomes. To explore the impact of sequencing platforms on bacterial genome assembly, quality assessment, sequence alignment, functional annotation, mutation detection, and metagenome mapping, we compared genome assemblies based on sequencing of cultured bacterial species using the HiSeq 2000 and BGISEQ-500 platforms. In addition, simulated reads were used to evaluate the impact of insert size on genome assembly. Genome assemblies based on BGISEQ-500 sequencing exhibited higher completeness and fewer N bases in high GC genomes, whereas HiSeq 2000 assemblies exhibited higher N50. The majority of assembly assessment parameters, sequences of 16S rRNA genes and genomes, numbers of single nucleotide variants (SNV), and mapping to metagenome data did not differ significantly between platforms. More insertions were detected in HiSeq 2000 genome assemblies, whereas more deletions were detected in BGISEQ-500 genome assemblies. Insert size had no significant impact on genome assembly. Taken together, our results suggest that DNBSEQ platforms would be a valid substitute for HiSeq 2000 for bacterial genome sequencing.

Complete Mitochondrial Genome and Phylogenetic Analysis of Tarsiger indicus (Aves: Passeriformes: Muscicapidae)

Tarsiger indicus (Vieillot, 1817), the White-browed Bush Robin, is a small passerine bird widely distributed in Asian countries. Here, we successfully sequenced its mitogenome using the Illumina Novaseq 6000 platform (Illumina, San Diego, CA, USA) for PE 2 × 150 bp sequencing. Combined with other published mitogenomes, we conducted the first comprehensive comparative mitogenome analysis of Muscicapidae birds and reconstructed the phylogenetic relationships between Muscicapidae and related groups. The T. indicus mitogenome was 16,723 bp in size, and it possessed the typical avian mitogenome structure and organization. Most PCGs of T. indicus were initiated strictly with the typical start codon ATG, while COX1 and ND2 were started with GTG. RSCU statistics showed that CUA, CGA, and GCC were relatively high frequency in the T. indicus mitogenome. T. cyanurus and T. indicus shared very similar mitogenomic features. All 13 PCGs of Muscicapidae mitogenomes had experienced purifying selection. Specifically, ATP8 had the highest rate of evolution (0.13296), whereas COX1 had the lowest (0.01373). The monophylies of Muscicapidae, Turdidae, and Paradoxornithidae were strongly supported. The clade of ((Muscicapidae + Turdidae) + Sturnidae) in Passeriformes was supported by both Bayesian Inference and Maximum likelihood analyses. The latest taxonomic status of many passerine birds with complex taxonomic histories were also supported. For example, Monticola gularis, T. indicus, and T. cyanurus were allocated to Turdidae in other literature; our phylogenetic topologies clearly supported their membership in Muscicapidae; Paradoxornis heudei, Suthora webbiana, S. nipalensis, and S. fulvifrons were formerly classified into Muscicapidae; we supported their membership in Paradoxornithidae; Culicicapa ceylonensis was originally classified as a member of Muscicapidae; our results are consistent with a position in Stenostiridae. Our study enriches the genetic data of T. indicus and provides new insights into the molecular phylogeny and evolution of passerine birds.

Genomic Architecture Predicts Tree Topology, Population Structuring, and Demographic History in Amazonian Birds

Geographic barriers are frequently invoked to explain genetic structuring across the landscape. However, inferences on the spatial and temporal origins of population variation have been largely limited to evolutionary neutral models, ignoring the potential role of natural selection and intrinsic genomic processes known as genomic architecture in producing heterogeneity in differentiation across the genome. To test how variation in genomic characteristics (e.g. recombination rate) impacts our ability to reconstruct general patterns of differentiation between species that cooccur across geographic barriers, we sequenced the whole genomes of multiple bird populations that are distributed across rivers in southeastern Amazonia. We found that phylogenetic relationships within species and demographic parameters varied across the genome in predictable ways. Genetic diversity was positively associated with recombination rate and negatively associated with species tree support. Gene flow was less pervasive in genomic regions of low recombination, making these windows more likely to retain patterns of population structuring that matched the species tree. We further found that approximately a third of the genome showed evidence of selective sweeps and linked selection, skewing genome-wide estimates of effective population sizes and gene flow between populations toward lower values. In sum, we showed that the effects of intrinsic genomic characteristics and selection can be disentangled from neutral processes to elucidate spatial patterns of population differentiation.

A Fast, Reproducible, High-throughput Variant Calling Workflow for Population Genomics

The increasing availability of genomic resequencing data sets and high-quality reference genomes across the tree of life present exciting opportunities for comparative population genomic studies. However, substantial challenges prevent the simple reuse of data across different studies and species, arising from variability in variant calling pipelines, data quality, and the need for computationally intensive reanalysis. Here, we present snpArcher, a flexible and highly efficient workflow designed for the analysis of genomic resequencing data in nonmodel organisms. snpArcher provides a standardized variant calling pipeline and includes modules for variant quality control, data visualization, variant filtering, and other downstream analyses. Implemented in Snakemake, snpArcher is user-friendly, reproducible, and designed to be compatible with high-performance computing clusters and cloud environments. To demonstrate the flexibility of this pipeline, we applied snpArcher to 26 public resequencing data sets from nonmammalian vertebrates. These variant data sets are hosted publicly to enable future comparative population genomic analyses. With its extensibility and the availability of public data sets, snpArcher will contribute to a broader understanding of genetic variation across species by facilitating the rapid use and reuse of large genomic data sets.

Gene regulation and speciation in a migratory divide between songbirds

Behavioral variation abounds in nature. This variation is important for adaptation and speciation, but its molecular basis remains elusive. Here, we use a hybrid zone between two subspecies of songbirds that differ in migration - an ecologically important and taxonomically widespread behavior---to gain insight into this topic. We measure gene expression in five brain regions. Differential expression between migratory states was dominated by circadian genes in all brain regions. The remaining patterns were largely brain-region specific. For example, expression differences between the subspecies that interact with migratory state likely help maintain reproductive isolation in this system and were documented in only three brain regions. Contrary to existing work on regulatory mechanisms underlying species-specific traits, two lines of evidence suggest that trans- (vs. cis) regulatory changes underlie these patterns - no evidence for allele-specific expression in hybrids and minimal associations between genomic differentiation and expression differences. Additional work with hybrids shows expression levels were often distinct (transgressive) from parental forms. Behavioral contrasts and functional enrichment analyses allowed us to connect these patterns to mitonuclear incompatibilities and compensatory responses to stress that could exacerbate selection on hybrids and contribute to speciation.

The MetaInvert soil invertebrate genome resource provides insights into below-ground biodiversity and evolution

Soil invertebrates are among the least understood metazoans on Earth. Thus far, the lack of taxonomically broad and dense genomic resources has made it hard to thoroughly investigate their evolution and ecology. With MetaInvert we provide draft genome assemblies for 232 soil invertebrate species, representing 14 common groups and 94 families. We show that this data substantially extends the taxonomic scope of DNA- or RNA-based taxonomic identification. Moreover, we confirm that theories of genome evolution cannot be generalised across evolutionarily distinct invertebrate groups. The soil invertebrate genomes presented here will support the management of soil biodiversity through molecular monitoring of community composition and function, and the discovery of evolutionary adaptations to the challenges of soil conditions.

The draft genome of the Temminck's tragopan (Tragopan temminckii) with evolutionary implications

BACKGROUND

High-quality genome data of birds play a significant role in the systematic study of their origin and adaptive evolution. The Temminck's tragopan (Tragopan temminckii) (Galliformes, Phasianidae), a larger pheasant, is one of the most abundant and widely distributed species of the genus Tragopan, and was defined as class II of the list of national key protected wild animals in China. The absence of a sequenced genome has restricted previous evolutionary trait studies of this taxa.

RESULTS

The whole genome of the Temminck's tragopan was sequenced using Illumina and PacBio platform, and then de novo assembled and annotated. The genome size was 1.06 Gb, with a contig N50 of 4.17 Mb. A total of 117.22 Mb (11.00%) repeat sequences were identified. 16,414 genes were predicted using three methods, with 16,099 (98.08%) annotated as functional genes based on five databases. In addition, comparative genome analyses were conducted across 12 Galliformes species. The results indicated that T. temminckii was the first species to branch off from the clade containing Lophura nycthemera, Phasianus colchicus, Chrysolophus pictus, Syrmaticus mikado, Perdix hodgsoniae, and Meleagris gallopavo, with a corresponding divergence time of 31.43 million years ago (MYA). Expanded gene families associated with immune response and energy metabolism were identified. Genes and pathways associated with plumage color and feather development, immune response, and energy metabolism were found in the list of positively selected genes (PSGs).

CONCLUSIONS

A genome draft of the Temminck's tragopan was reported, genome feature and comparative genome analysis were described, and genes and pathways related to plumage color and feather development, immune response, and energy metabolism were identified. The genomic data of the Temminck's tragopan considerably contribute to the genome evolution and phylogeny of the genus Tragopan and the whole Galliformes species underlying ecological adaptation strategies.

Leveraging genomic redundancy to improve inference and alignment of orthologous proteins

Identifying protein sequences with common ancestry is a core task in bioinformatics and evolutionary biology. However, methods for inferring and aligning such sequences in annotated genomes have not kept pace with the increasing scale and complexity of the available data. Thus, in this work, we implemented several improvements to the traditional methodology that more fully leverage the redundancy of closely related genomes and the organization of their annotations. Two highlights include the application of the more flexible k-clique percolation algorithm for identifying clusters of orthologous proteins and the development of a novel technique for removing poorly supported regions of alignments with a phylogenetic hidden Markov model (phylo-HMM). In making the latter, we wrote a fully documented Python package Homomorph that implements standard HMM algorithms and created a set of tutorials to promote its use by a wide audience. We applied the resulting pipeline to a set of 33 annotated Drosophila genomes, generating 22,813 orthologous groups and 8,566 high-quality alignments.

Low Mutation Load in a Supergene Underpinning Alternative Male Mating Strategies in Ruff (Calidris pugnax)

A paradox in evolutionary biology is how supergenes can maintain high fitness despite reduced effective population size, the suppression of recombination, and the expected accumulation of mutational load. The ruff supergene involves 2 rare inversion haplotypes (satellite and faeder). These are recessive lethals but with dominant effects on male mating strategies, plumage, and body size. Sequence divergence to the wild-type (independent) haplotype indicates that the inversion could be as old as 4 million years. Here, we have constructed a highly contiguous genome assembly of the inversion region for both the independent and satellite haplotypes. Based on the new data, we estimate that the recombination event(s) creating the satellite haplotype occurred only about 70,000 yr ago. Contrary to expectations for supergenes, we find no substantial expansion of repeats and only a modest mutation load on the satellite and faeder haplotypes despite high sequence divergence to the non-inverted haplotype (1.46%). The essential centromere protein N (CENPN) gene is disrupted by the inversion and is as well conserved on the inversion haplotypes as on the noninversion haplotype. These results suggest that the inversion may be much younger than previously thought. The low mutation load, despite recessive lethality, may be explained by the introgression of the inversion from a now extinct lineage.

Genomic divergence and introgression between cryptic species of a widespread North American songbird

Analysis of genomic variation among related populations can sometimes reveal distinct species that were previously undescribed due to similar morphological appearances, and close examination of such cases can provide much insight regarding speciation. Genomic data can also reveal the role of reticulate evolution in differentiation and speciation. White-breasted nuthatches (Sitta carolinensis) are widely distributed North American songbirds that are currently classified as a single species but have been suspected to represent a case of cryptic speciation. Previous genetic analyses suggested four divergent groups, but it was unclear whether these represented multiple reproductively isolated species. Using extensive genomic sampling of over 350 white-breasted nuthatches from across North America and a new chromosome-level reference genome, we asked if white-breasted nuthatches are comprised of multiple species and whether introgression has occurred between divergent populations. Genomic variation of over 300,000 loci revealed four highly differentiated populations (Pacific, n = 45; Eastern, n = 23; Rocky Mountains North, n = 138; and Rocky Mountains South, n = 150) with geographic ranges that are adjacent. We observed a moderate degree of admixture between Rocky Mountain populations but only a small number of hybrids between the Rockies and the Eastern population. The rarity of hybrids together with high levels of differentiation between populations is supportive of populations having some level of reproductive isolation. Between populations, we show evidence for introgression from a divergent ghost lineage of white-breasted nuthatches into the Rocky Mountains South population, which is otherwise closely related to Rocky Mountains North. We conclude that white-breasted nuthatches are best considered at least three species and that ghost lineage introgression has contributed to differentiation between the two Rocky Mountain populations. White-breasted nuthatches provide a dramatic case of morphological similarity despite high genomic differentiation, and the varying levels of reproductive isolation among the four groups provide an example of the speciation continuum.

Genetic Load and Adaptive Potential of a Recovered Avian Species that Narrowly Avoided Extinction

High genetic diversity is a good predictor of long-term population viability, yet some species persevere despite having low genetic diversity. Here we study the genomic erosion of the Seychelles paradise flycatcher (Terpsiphone corvina), a species that narrowly avoided extinction after having declined to 28 individuals in the 1960s. The species recovered unassisted to over 250 individuals in the 1990s and was downlisted from Critically Endangered to Vulnerable in the International Union for the Conservation of Nature Red List in 2020. By comparing historical, prebottleneck (130+ years old) and modern genomes, we uncovered a 10-fold loss of genetic diversity. Highly deleterious mutations were partly purged during the bottleneck, but mildly deleterious mutations accumulated. The genome shows signs of historical inbreeding during the bottleneck in the 1960s, but low levels of recent inbreeding after demographic recovery. Computer simulations suggest that the species long-term small Ne reduced the masked genetic load and made the species more resilient to inbreeding and extinction. However, the reduction in genetic diversity due to the chronically small Ne and the severe bottleneck is likely to have reduced the species adaptive potential to face environmental change, which together with a higher load, compromises its long-term population viability. Thus, small ancestral Ne offers short-term bottleneck resilience but hampers long-term adaptability to environmental shifts. In light of rapid global rates of population decline, our work shows that species can continue to suffer the effect of their decline even after recovery, highlighting the importance of considering genomic erosion and computer modeling in conservation assessments.

De novo assembly of a chromosome-level reference genome for the California Scrub-Jay, Aphelocoma californica

We announce the assembly of the first de novo reference genome for the California Scrub-Jay (Aphelocoma californica). The genus Aphelocoma comprises four currently recognized species including many locally adapted populations across Mesoamerica and North America. Intensive study of Aphelocoma has revealed novel insights into the evolutionary mechanisms driving diversification in natural systems. Additional insights into the evolutionary history of this group will require continued development of high-quality, publicly available genomic resources. We extracted high molecular weight genomic DNA from a female California Scrub-Jay from northern California and generated PacBio HiFi long-read data and Omni-C chromatin conformation capture data. We used these data to generate a de novo partially phased diploid genome assembly, consisting of two pseudo-haplotypes, and scaffolded them using inferred physical proximity information from the Omni-C data. The more complete pseudo-haplotype assembly (arbitrarily designated "Haplotype 1") is 1.35 Gb in total length, highly contiguous (contig N50 = 11.53 Mb), and highly complete (BUSCO completeness score = 97%), with comparable scaffold sizes to chromosome-level avian reference genomes (scaffold N50 = 66.14 Mb). Our California Scrub-Jay assembly is highly syntenic with the New Caledonian Crow reference genome despite ~10 million years of divergence, highlighting the temporal stability of the avian genome. This high-quality reference genome represents a leap forward in publicly available genomic resources for Aphelocoma, and the family Corvidae more broadly. Future work using Aphelocoma as a model for understanding the evolutionary forces generating and maintaining biodiversity across phylogenetic scales can now benefit from a highly contiguous, in-group reference genome.

Candidate genes under selection in song sparrows co-vary with climate and body mass in support of Bergmann's Rule

Ecogeographic rules denote spatial patterns in phenotype and environment that may reflect local adaptation as well as a species' capacity to adapt to change. To identify genes underlying Bergmann's Rule, which posits that spatial correlations of body mass and temperature reflect natural selection and local adaptation in endotherms, we compare 79 genomes from nine song sparrow (Melospiza melodia) subspecies that vary ~300% in body mass (17 - 50 g). Comparing large- and smaller-bodied subspecies revealed 9 candidate genes in three genomic regions associated with body mass. Further comparisons to the five smallest subspecies endemic to California revealed eight SNPs within four of the candidate genes (GARNL3, RALGPS1, ANGPTL2, and COL15A1) associated with body mass and varying as predicted by Bergmann's Rule. Our results support the hypothesis that co-variation in environment, body mass and genotype reflect the influence of natural selection on local adaptation and a capacity for contemporary evolution in this diverse species.

A Bird's-Eye View of Endangered Species Conservation: Avian Genomics and Stem Cell Approaches for Green Peafowl (Pavo muticus)

Aves ranks among the top two classes for the highest number of endangered and extinct species in the kingdom Animalia. Notably, the IUCN Red List classified the green peafowl as endangered. This highlights promising strategies using genetics and reproductive technologies for avian wildlife conservation. These platforms provide the capacity to predict population trends and enable the practical breeding of such species. The conservation of endangered avian species is facilitated through the application of genomic data storage and analysis. Storing the sequence is a form of biobanking. An analysis of sequence can identify genetically distinct individuals for breeding. Here, we reviewed avian genomics and stem cell approaches which not only offer hope for saving endangered species, such as the green peafowl but also for other birds threatened with extinction.

3D atlas of tinamou (Neornithes: Tinamidae) pectoral morphology: Implications for reconstructing the ancestral neornithine flight apparatus

Palaeognathae, the extant avian clade comprising the flightless ratites and flight-capable tinamous (Tinamidae), is the sister group to all other living birds, and recent phylogenetic studies illustrate that tinamous are phylogenetically nested within a paraphyletic assemblage of ratites. As the only extant palaeognaths that have retained the ability to fly, tinamous may provide key information on the nature of the flight apparatus of ancestral crown palaeognaths-and, in turn, crown birds-as well as insight into convergent modifications to the wing apparatus among extant ratite lineages. To reveal new information about the musculoskeletal anatomy of tinamous and facilitate development of computational biomechanical models of tinamou wing function, we generated a three-dimensional musculoskeletal model of the flight apparatus of the extant Andean tinamou (Nothoprocta pentlandii) using diffusible iodine-based contrast-enhanced computed tomography (diceCT). Origins and insertions of the pectoral flight musculature of N. pentlandii are generally consistent with those of other extant volant birds specialized for burst flight, and the entire suite of presumed ancestral neornithine flight muscles are present in N. pentlandii with the exception of the biceps slip. The pectoralis and supracoracoideus muscles are robust, similar to the condition in other extant burst-flying birds such as many extant Galliformes. Contrary to the condition in most extant Neognathae (the sister clade to Palaeognathae), the insertion of the pronator superficialis has a greater distal extent than the pronator profundus, although most other anatomical observations are broadly consistent with the conditions observed in extant neognaths. This work will help form a basis for future comparative studies of the avian musculoskeletal system, with implications for reconstructing the flight apparatus of ancestral crown birds and clarifying musculoskeletal modifications underlying the convergent origins of ratite flightlessness.

Phylogenomic insights into the polyphyletic nature of Altai falcons within eastern sakers (Falco cherrug) and the origins of gyrfalcons (Falco rusticolus)

The Altai falcon from Central Asia always attracted the attention of humans. Long considered a totemic bird in its native area, modern falconers still much appreciated this large-bodied and mighty bird of prey due to its rarity and unique look. The peculiar body characteristics halfway between the saker falcon (Falco cherrug) and the gyrfalcon (F. rusticolus) triggered debates about its contentious taxonomy. The weak phylogenetic signal associated with traditional genetic methods could not resolve this uncertainty. Here, we address the controversial evolutionary origin of Altai falcons by means of a genome-wide approach, Restriction-site Associated DNA sequencing, using sympatric eastern sakers falcons, allopatric western saker falcons and gyrfalcons as outgroup. This approach provided an unprecedented insight into the phylogenetic relationships of the studied populations by delivering 17,095 unlinked SNPs shedding light on the polyphyletic nature of Altai falcons within eastern sakers. Thus we concluded that the former must correspond to a low taxonomic rank, probably an ecotype or form of the latter. Also, we found that eastern sakers are paraphyletic without gyrfalcons, thus, these latter birds are best regarded as the direct sister lineage of the eastern sakers. This evolutionary relationship, corroborated also by re-analyzing the dataset with the inclusion of outgroup samples (F. biarmicus and F. peregrinus), put eastern sakers into a new light as the potential ancestral genetic source of high latitude and altitude adaptation in descendent populations. Finally, conservation genomic values hint at the stable genetic background of the studied saker populations.

Novel tRNA Gene Rearrangements in the Mitochondrial Genomes of Poneroid Ants and Phylogenetic Implication of Paraponerinae (Hymenoptera: Formicidae)

Ants (Formicidae) are the most diverse eusocial insects in Hymenoptera, distributed across 17 extant subfamilies grouped into 3 major clades, the Formicoid, Leptanilloid, and Poneroid. While the mitogenomes of Formicoid ants have been well studied, there is a lack of published data on the mitogenomes of Poneroid ants, which requires further characterization. In this study, we first present three complete mitogenomes of Poneroid ants: Paraponera clavata, the only extant species from the subfamily Paraponerinae, and two species (Harpegnathos venator and Buniapone amblyops) from the Ponerinae subfamily. Notable novel gene rearrangements were observed in the new mitogenomes, located in the gene blocks CR-trnM-trnI-trnQ-ND2, COX1-trnK-trnD-ATP8, and ND3-trnA-trnR-trnN-trnS1-trnE-trnF-ND5. We reported the duplication of tRNA genes for the first time in Formicidae. An extra trnQ gene was identified in H. venator. These gene rearrangements could be explained by the tandem duplication/random loss (TDRL) model and the slipped-strand mispairing model. Additionally, one large duplicated region containing tandem repeats was identified in the control region of P. clavata. Phylogenetic analyses based on protein-coding genes and rRNA genes via maximum likelihood and Bayes methods supported the monophyly of the Poneroid clade and the sister group relationship between the subfamilies Paraponerinae and Amblyoponinae. However, caution is advised in interpreting the positions of Paraponerinae due to the potential artifact of long-branch attraction.

PhyloMatcher: a tool for resolving conflicts in taxonomic nomenclature

Summary

Large-scale comparative studies rely on the application of both phylogenetic trees and phenotypic data, both of which come from a variety of sources, but due to the changing nature of phylogenetic classification over time, many taxon names in comparative datasets do not match the nomenclature in phylogenetic trees. Manual curation of taxonomic synonyms in large comparative datasets can be daunting. To address this issue, we introduce PhyloMatcher, a tool which allows for programmatic querying of the National Center for Biotechnology Information Taxonomy and Global Biodiversity Information Facility databases to find associated synonyms with given target species names.

Availability and implementation

PhyloMatcher is easily installed as a Python package with pip, or as a standalone GUI application. PhyloMatcher source code and documentation are freely available at https://github.com/Lswhiteh/PhyloMatcher, the GUI application can be downloaded from the Releases page.

The avian taste system

Taste or gustation is the sense evolving from the chemo-sensory system present in the oral cavity of avian species, which evolved to evaluate the nutritional value of foods by detecting relevant compounds including amino acids and peptides, carbohydrates, lipids, calcium, salts, and toxic or anti-nutritional compounds. In birds compared to mammals, due to the relatively low retention time of food in the oral cavity, the lack of taste papillae in the tongue, and an extremely limited secretion of saliva, the relevance of the avian taste system has been historically undermined. However, in recent years, novel data has emerged, facilitated partially by the advent of the genomic era, evidencing that the taste system is as crucial to avian species as is to mammals. Despite many similarities, there are also fundamental differences between avian and mammalian taste systems in terms of anatomy, distribution of taste buds, and the nature and molecular structure of taste receptors. Generally, birds have smaller oral cavities and a lower number of taste buds compared to mammals, and their distribution in the oral cavity appears to follow the swallowing pattern of foods. In addition, differences between bird species in the size, structure and distribution of taste buds seem to be associated with diet type and other ecological adaptations. Birds also seem to have a smaller repertoire of bitter taste receptors (T2Rs) and lack some taste receptors such as the T1R2 involved in sweet taste perception. This has opened new areas of research focusing on taste perception mechanisms independent of GPCR taste receptors and the discovery of evolutionary shifts in the molecular function of taste receptors adapting to ecological niches in birds. For example, recent discoveries have shown that the amino acid taste receptor dimer T1R1-T1R3 have mutated to sense simple sugars in almost half of the living bird species, or SGLT1 has been proposed as a part of a T1R2-independent sweet taste sensing in chicken. The aim of this review is to present the scientific data known to date related to the avian taste system across species and its impact on dietary choices including domestic and wild species.

Breaking Through the Bottleneck: Krogh's Principle in Behavioral Neuroendocrinology and the Potential of Gene Editing

In 1929, August Krogh wrote that for every question in biology, there is a species or collection of species in which pursuing such questions is the most appropriate for achieving the deepest insights. Referred to as "Krogh's Principle," these words are a guiding force for many biologists. In practice, Krogh's principle might guide a biologist interested in studying bi-parental care to choose not to use lab mice, in which the female does most of the parenting, but instead study species in which bi-parental care is present and clearly observable, such as in certain poison dart frogs. This approach to pursuing biological questions has been fruitful, with more in-depth insights achievable with new technologies. However, up until recently, an important limitation of Krogh's principle for biologists interested in the functions of certain genes, was certain techniques were only available for a few traditional model organisms such as lab mice, fruit flies (Drosophila melanogaster), zebrafish (Danio rerio) and C. elegans (Caenorhabditis elegans), in which testing the functions of molecular systems on biological processes can be achieved using genetic knockout (KO) and transgenic technology. These methods are typically more precise than other approaches (e.g., pharmacology) commonly used in nontraditional model organisms to address similar questions. Therefore, some of the most in-depth insights into our understanding of the molecular control of these mechanisms have come from a small number of genetically tractable species. Recent advances in gene editing technology such as CRISPR (Clustered Regularly Interspersed Short Palindromic Repeats)/Cas9 gene editing as a laboratory tool has changed the insights achievable for biologists applying Krogh's principle. In this review, we will provide a brief summary on how some researchers of nontraditional model organisms have been able to achieve different levels of experimental precision with limited genetic tractability in their non-traditional model organism in the field of behavioral neuroendocrinology, a field in which understanding tissue and brain-region specific actions of molecules of interest has been a major goal. Then, we will highlight the exciting potential of Krogh's principle using discoveries made in a popular model species of social behavior, the African cichlid fish Astatotilapia burtoni. Specifically, we will focus on insights gained from studies of the control of social status by sex steroid hormones (androgens and estrogens) in A. burtoni that originated during field observations during the 1970s, and have recently culminated in novel insights from CRISPR/Cas9 gene editing in laboratory studies. Our review highlighting discoveries in A. burtoni may function as a roadmap for others using Krogh's principle aiming to incorporate gene editing into their research program. Gene editing is thus a powerful complimentary laboratory tool researchers can use to yield novel insights into understanding the molecular mechanisms of physiology and behavior in non-traditional model organisms.

A highly contiguous reference genome for the Steller's jay (Cyanocitta stelleri)

The Steller's jay is a familiar bird of western forests from Alaska south to Nicaragua. Here, we report a draft reference assembly for the species generated from PacBio HiFi long-read and Omni-C chromatin-proximity sequencing data as part of the California Conservation Genomics Project (CCGP). Sequenced reads were assembled into 352 scaffolds totaling 1.16 Gb in length. Assembly metrics indicate a highly contiguous and complete assembly with a contig N50 of 7.8 Mb, scaffold N50 of 25.8 Mb, and BUSCO completeness score of 97.2%. Repetitive elements span 16.6% of the genome including nearly 90% of the W chromosome. Compared with high-quality assemblies from other members of the family Corvidae, the Steller's jay genome contains a larger proportion of repetitive elements than 4 crow species (Corvus), but a lower proportion of repetitive elements than the California scrub-jay (Aphelocoma californica). This reference genome will serve as an essential resource for future studies on speciation, local adaptation, phylogeography, and conservation genetics in this species of significant biological interest.

Genome mining shows that retroviruses are pervasively invading vertebrate genomes

Endogenous retroviruses (ERVs) record past retroviral infections, providing molecular archives for interrogating the evolution of retroviruses and retrovirus-host interaction. However, the vast majority of ERVs are not active anymore due to various disruptive mutations, and ongoing retroviral invasion of vertebrate genomes has been rarely documented. Here we analyze genomics data from 2004 vertebrates for mining invading ERVs (ERVi). We find that at least 412 ERVi elements representing 217 viral operational taxonomic units are invading the genomes of 123 vertebrates, 18 of which have been assessed to be threatened species. Our results reveal an unexpected prevalence of ongoing retroviral invasion in vertebrates and expand the diversity of retroviruses recently circulating in the wild. We characterize the pattern and nature of ERVi in the historical and biogeographical context of their hosts, for instance, the generation of model organisms, sympatric speciation, and domestication. We suspect that these ERVi are relevant to conservation of threatened species, zoonoses in the wild, and emerging infectious diseases in humans.

Considering admixture when producing draft genomes: an example in North American ratsnakes (Pantherophis alleghaniensis/Pantherophis obsoletus)

The number of reference genomes of snakes lags behind several other vertebrate groups (e.g. birds and mammals). However, in the last two years, a concerted effort by researchers from around the world has produced new genomes of snakes representing members from several new families. Here, we present a high-quality, annotated genome of the central ratsnake (Pantherophis alleghaniensis), a member of the most diverse snake lineage, Colubroidea. Pantherophis alleghaniensis is found in the central part of the Nearctic, east of the Mississippi River. This genome was sequenced using 10X Chromium synthetic long reads and polished using Illumina short reads. The final genome assembly had an N50 of 21.82 Mb and an L50 of 22 scaffolds with a maximum scaffold length of 82.078 Mb. The genome is composed of 49.24% repeat elements dominated by long interspersed elements. We annotated this genome using transcriptome assemblies from 14 tissue types and recovered 28,368 predicted proteins. Finally, we estimated admixture proportions between two species of ratsnakes and discovered that this specimen is an admixed individual containing genomes from the western (Pantherophis obsoletus) and central ratsnakes (P. alleghaniensis). We discuss the importance of considering interspecific admixture in downstream approaches for inferring demography and phylogeny.

PhyloMatcher: a tool for resolving conflicts in taxonomic nomenclature

Summary

Large-scale comparative studies rely on the application of both phylogenetic trees and phenotypic data, both of which come from a variety of sources, but due to the changing nature of phylogenetic classification over time, many taxon names in comparative datasets do not match the nomenclature in phylogenetic trees. Manual curation of taxonomic synonyms in large comparative datasets can be daunting. To address this issue, we introduce PhyloMatcher, a tool which allows for programmatic querying of two commonly used taxonomic databases to find associated synonyms with given target species names.

Availability and implementation

PhyloMatcher is easily installed as a Python package with pip, or as a standalone GUI application. PhyloMatcher source code and documentation are freely available at https://github.com/Lswhiteh/PhyloMatcher, the GUI application can be downloaded from the Releases page.

Contact

Lswhiteh@unc.edu.

Supplemental Information

We provide documentation for PhyloMatcher, including walkthrough instructions for the GUI application on the Releases page of https://github.com/Lswhiteh/PhyloMatcher.

A Chromosome-Level Reference Genome for the Black-Legged Kittiwake (Rissa tridactyla), a Declining Circumpolar Seabird

Amidst the current biodiversity crisis, the availability of genomic resources for declining species can provide important insights into the factors driving population decline. In the early 1990s, the black-legged kittiwake (Rissa tridactyla), a pelagic gull widely distributed across the arctic, subarctic, and temperate zones, suffered a steep population decline following an abrupt warming of sea surface temperature across its distribution range and is currently listed as Vulnerable by the International Union for the Conservation of Nature. Kittiwakes have long been the focus for field studies of physiology, ecology, and ecotoxicology and are primary indicators of fluctuating ecological conditions in arctic and subarctic marine ecosystems. We present a high-quality chromosome-level reference genome and annotation for the black-legged kittiwake using a combination of Pacific Biosciences HiFi sequencing, Bionano optical maps, Hi-C reads, and RNA-Seq data. The final assembly spans 1.35 Gb across 32 chromosomes, with a scaffold N50 of 88.21 Mb and a BUSCO completeness of 97.4%. This genome assembly substantially improves the quality of a previous draft genome, showing an approximately 5× increase in contiguity and a more complete annotation. Using this new chromosome-level reference genome and three more chromosome-level assemblies of Charadriiformes, we uncover several lineage-specific chromosome fusions and fissions, but find no shared rearrangements, suggesting that interchromosomal rearrangements have been commonplace throughout the diversification of Charadriiformes. This new high-quality genome assembly will enable population genomic, transcriptomic, and phenotype-genotype association studies in a widely studied sentinel species, which may provide important insights into the impacts of global change on marine systems.

Mutational Study of the Tryptophan Tetrad Important for Electron Transfer in European Robin Cryptochrome 4a

The ability of migratory birds to sense magnetic fields has been known for decades, although the understanding of the underlying mechanism is still elusive. Currently, the strongest magnetoreceptor candidate in birds is a protein called cryptochrome 4a. The cryptochrome 4a protein has changed through evolution, apparently endowing some birds with a more pronounced magnetic sensitivity than others. Using phylogenetic tools, we show that a specific tryptophan tetrad and a tyrosine residue predicted to be essential for cryptochrome activation are highly conserved in the avian clade. Through state-of-the-art molecular dynamics simulations and associated analyses, we also studied the role of these specific residues and the associated mutants on the overall dynamics of the protein. The analyses of the single residue mutations were used to judge how far a local change in the protein structure can impact specific dynamics of European robin cryptochrome 4a. We conclude that the replacements of each of the tryptophans one by one with a phenylalanine do not compromise the overall stability of the protein.

A high-quality genome assembly highlights the evolutionary history of the great bustard (Otis tarda, Otidiformes)

Conservation genomics often relies on non-invasive methods to obtain DNA fragments which limit the power of multi-omic analyses for threatened species. Here, we report multi-omic analyses based on a well-preserved great bustard individual (Otis tarda, Otidiformes) that was found dead in the mountainous region in Gansu, China. We generate a near-complete genome assembly containing only 18 gaps scattering in 8 out of the 40 assembled chromosomes. We characterize the DNA methylation landscape which is correlated with GC content and gene expression. Our phylogenomic analysis suggests Otidiformes and Musophagiformes are sister groups that diverged from each other 46.3 million years ago. The genetic diversity of great bustard is found the lowest among the four available Otidiformes genomes, possibly due to population declines during past glacial periods. As one of the heaviest migratory birds, great bustard possesses several expanded gene families related to cardiac contraction, actin contraction, calcium ion signaling transduction, as well as positively selected genes enriched for metabolism. Finally, we identify an extremely young evolutionary stratum on the sex chromosome, a rare case among birds. Together, our study provides insights into the conservation genomics, adaption and chromosome evolution of the great bustard.

Hybridisation has shaped a recent radiation of grass-feeding aphids

BACKGROUND

Aphids are common crop pests. These insects reproduce by facultative parthenogenesis involving several rounds of clonal reproduction interspersed with an occasional sexual cycle. Furthermore, clonal aphids give birth to live young that are already pregnant. These qualities enable rapid population growth and have facilitated the colonisation of crops globally. In several cases, so-called "super clones" have come to dominate agricultural systems. However, the extent to which the sexual stage of the aphid life cycle has shaped global pest populations has remained unclear, as have the origins of successful lineages. Here, we used chromosome-scale genome assemblies to disentangle the evolution of two global pests of cereals-the English (Sitobion avenae) and Indian (Sitobion miscanthi) grain aphids.

RESULTS

Genome-wide divergence between S. avenae and S. miscanthi is low. Moreover, comparison of haplotype-resolved assemblies revealed that the S. miscanthi isolate used for genome sequencing is likely a hybrid, with one of its diploid genome copies closely related to S. avenae (~ 0.5% divergence) and the other substantially more divergent (> 1%). Population genomics analyses of UK and China grain aphids showed that S. avenae and S. miscanthi are part of a cryptic species complex with many highly differentiated lineages that predate the origins of agriculture. The complex consists of hybrid lineages that display a tangled history of hybridisation and genetic introgression.

CONCLUSIONS

Our analyses reveal that hybridisation has substantially contributed to grain aphid diversity, and hence, to the evolutionary potential of this important pest species. Furthermore, we propose that aphids are particularly well placed to exploit hybridisation events via the rapid propagation of live-born "frozen hybrids" via asexual reproduction, increasing the likelihood of hybrid lineage formation.

The State of Squamate Genomics: Past, Present, and Future of Genome Research in the Most Speciose Terrestrial Vertebrate Order

Squamates include more than 11,000 extant species of lizards, snakes, and amphisbaenians, and display a dazzling diversity of phenotypes across their over 200-million-year evolutionary history on Earth. Here, we introduce and define squamates (Order Squamata) and review the history and promise of genomic investigations into the patterns and processes governing squamate evolution, given recent technological advances in DNA sequencing, genome assembly, and evolutionary analysis. We survey the most recently available whole genome assemblies for squamates, including the taxonomic distribution of available squamate genomes, and assess their quality metrics and usefulness for research. We then focus on disagreements in squamate phylogenetic inference, how methods of high-throughput phylogenomics affect these inferences, and demonstrate the promise of whole genomes to settle or sustain persistent phylogenetic arguments for squamates. We review the role transposable elements play in vertebrate evolution, methods of transposable element annotation and analysis, and further demonstrate that through the understanding of the diversity, abundance, and activity of transposable elements in squamate genomes, squamates can be an ideal model for the evolution of genome size and structure in vertebrates. We discuss how squamate genomes can contribute to other areas of biological research such as venom systems, studies of phenotypic evolution, and sex determination. Because they represent more than 30% of the living species of amniote, squamates deserve a genome consortium on par with recent efforts for other amniotes (i.e., mammals and birds) that aim to sequence most of the extant families in a clade.

A highly contiguous genome assembly for the California quail (Callipepla californica)

The California quail (Callipepla californica) is an iconic native bird of scrub and oak woodlands in California and the Baja Peninsula of Mexico. Here, we report a draft reference assembly for the species generated from PacBio HiFi long read and Omni-C chromatin-proximity sequencing data as part of the California Conservation Genomics Project (CCGP). Sequenced reads were assembled into 321 scaffolds totaling 1.08 Gb in length. Assembly metrics indicate a highly contiguous and complete assembly with a contig N50 of 5.5 Mb, scaffold N50 of 19.4 Mb, and BUSCO completeness score of 96.5%. Transposable elements (TEs) occupy 16.5% of the genome, more than previous Odontophoridae quail assemblies but in line with estimates of TE content for recent long-read assemblies of chicken and Peking duck. Together these metrics indicate that the present assembly is more complete than prior reference assemblies generated for Odontophoridae quail. This reference will serve as an essential resource for studies on local adaptation, phylogeography, and conservation genetics in this species of significant biological and recreational interest.

CNEReg Interprets Ruminant-specific Conserved Non-coding Elements by Developmental Gene Regulatory Network

The genetic information coded in DNA leads to trait innovation via a gene regulatory network (GRN) in development. Here, we developed a conserved non-coding element interpretation method to integrate multi-omics data into gene regulatory network (CNEReg) to investigate the ruminant multi-chambered stomach innovation. We generated paired expression and chromatin accessibility data during rumen and esophagus development in sheep, and revealed 1601 active ruminant-specific conserved non-coding elements (active-RSCNEs). To interpret the function of these active-RSCNEs, we defined toolkit transcription factors (TTFs) and modeled their regulation on rumen-specific genes via batteries of active-RSCNEs during development. Our developmental GRN revealed 18 TTFs and 313 active-RSCNEs regulating 7 rumen functional modules. Notably, 6 TTFs (OTX1, SOX21, HOXC8, SOX2, TP63, and PPARG), as well as 16 active-RSCNEs, functionally distinguished the rumen from the esophagus. Our study provides a systematic approach to understanding how gene regulation evolves and shapes complex traits by putting evo-devo concepts into practice with developmental multi-omics data.

Comparing genomic and epigenomic features across species using the WashU Comparative Epigenome Browser

Genome browsers have become an intuitive and critical tool to visualize and analyze genomic features and data. Conventional genome browsers display data/annotations on a single reference genome/assembly; there are also genomic alignment viewer/browsers that help users visualize alignment, mismatch, and rearrangement between syntenic regions. However, there is a growing need for a comparative epigenome browser that can display genomic and epigenomic data sets across different species and enable users to compare them between syntenic regions. Here, we present the WashU Comparative Epigenome Browser. It allows users to load functional genomic data sets/annotations mapped to different genomes and display them over syntenic regions simultaneously. The browser also displays genetic differences between the genomes from single-nucleotide variants (SNVs) to structural variants (SVs) to visualize the association between epigenomic differences and genetic differences. Instead of anchoring all data sets to the reference genome coordinates, it creates independent coordinates of different genome assemblies to faithfully present features and data mapped to different genomes. It uses a simple, intuitive genome-align track to illustrate the syntenic relationship between different species. It extends the widely used WashU Epigenome Browser infrastructure and can be expanded to support multiple species. This new browser function will greatly facilitate comparative genomic/epigenomic research, as well as support the recent growing needs to directly compare and benchmark the T2T CHM13 assembly and other human genome assemblies.

Integrating gene annotation with orthology inference at scale

Annotating coding genes and inferring orthologs are two classical challenges in genomics and evolutionary biology that have traditionally been approached separately, limiting scalability. We present TOGA (Tool to infer Orthologs from Genome Alignments), a method that integrates structural gene annotation and orthology inference. TOGA implements a different paradigm to infer orthologous loci, improves ortholog detection and annotation of conserved genes compared with state-of-the-art methods, and handles even highly fragmented assemblies. TOGA scales to hundreds of genomes, which we demonstrate by applying it to 488 placental mammal and 501 bird assemblies, creating the largest comparative gene resources so far. Additionally, TOGA detects gene losses, enables selection screens, and automatically provides a superior measure of mammalian genome quality. TOGA is a powerful and scalable method to annotate and compare genes in the genomic era.

Chromosome-level genome assembly of the critically endangered Baer's pochard (Aythya baeri)

Baer's pochard (Aythya baeri) is a critically endangered species historically widespread throughout East Asia, whose population according to a recent estimate has decreased to between 150 and 700 individuals, and faces a long-term risk of extinction. However, the lack of a reference genome limits the study of conservation management and molecular biology of this species. We therefore report the first high-quality genome assembly of Baer's pochard. The genome has a total length of 1.14 Gb with a scaffold N50 of 85,749,954 bp and a contig N50 of 29,098,202 bp. We anchored 97.88% of the scaffold sequences onto 35 chromosomes based on the Hi-C data. BUSCO assessment indicated that 97.00% of the highly conserved Aves genes were completely present in the genome assembly. Furthermore, a total of 157.06 Mb of repetitive sequences were identified and 18,581 protein-coding genes were predicted in the genome, of which 99.00% were functionally annotated. This genome will be useful for understanding Baer's pochard genetic diversity and facilitate the conservation planning of this species.

Evolution of the Growth Hormone Gene Duplication in Passerine Birds

Birds of the order Passeriformes represent the most speciose order of land vertebrates. Despite strong scientific interest in this super-radiation, genetic traits unique to passerines are not well characterized. A duplicate copy of growth hormone (GH) is the only gene known to be present in all major lineages of passerines, but not in other birds. GH genes plausibly influence extreme life history traits that passerines exhibit, including the shortest embryo-to-fledging developmental period of any avian order. To unravel the implications of this GH duplication, we investigated the molecular evolution of the ancestral avian GH gene (GH or GH1) and the novel passerine GH paralog (GH2), using 497 gene sequences extracted from 342 genomes. Passerine GH1 and GH2 are reciprocally monophyletic, consistent with a single duplication event from a microchromosome onto a macrochromosome in a common ancestor of extant passerines. Additional chromosomal rearrangements have changed the syntenic and potential regulatory context of these genes. Both passerine GH1 and GH2 display substantially higher rates of nonsynonymous codon change than non-passerine avian GH, suggesting positive selection following duplication. A site involved in signal peptide cleavage is under selection in both paralogs. Other sites under positive selection differ between the two paralogs, but many are clustered in one region of a 3D model of the protein. Both paralogs retain key functional features and are actively but differentially expressed in two major passerine suborders. These phenomena suggest that GH genes may be evolving novel adaptive roles in passerine birds.

Pathways to polar adaptation in fishes revealed by long-read sequencing

Long-read sequencing is driving a new reality for genome science in which highly contiguous assemblies can be produced efficiently with modest resources. Genome assemblies from long-read sequences are particularly exciting for understanding the evolution of complex genomic regions that are often difficult to assemble. In this study, we utilized long-read sequencing data to generate a high-quality genome assembly for an Antarctic eelpout, Ophthalmolycus amberensis, the first for the globally distributed family Zoarcidae. We used this assembly to understand how O. amberensis has adapted to the harsh Southern Ocean and compared it to another group of Antarctic fishes: the notothenioids. We showed that selection has largely acted on different targets in eelpouts relative to notothenioids. However, we did find some overlap; in both groups, genes involved in membrane structure, thermal tolerance and vision have evidence of positive selection. We found evidence for historical shifts of transposable element activity in O. amberensis and other polar fishes, perhaps reflecting a response to environmental change. We were specifically interested in the evolution of two complex genomic loci known to underlie key adaptations to polar seas: haemoglobin and antifreeze proteins (AFPs). We observed unique evolution of the haemoglobin MN cluster in eelpouts and related fishes in the suborder Zoarcoidei relative to other Perciformes. For AFPs, we identified the first species in the suborder with no evidence of afpIII sequences (Cebidichthys violaceus) in the genomic region where they are found in all other Zoarcoidei, potentially reflecting a lineage-specific loss of this cluster. Beyond polar fishes, our results highlight the power of long-read sequencing to understand genome evolution.

A de novo assembled genome of the Tibetan Partridge (Perdix hodgsoniae) and its high-altitude adaptation

The Tibetan Partridge (Perdix hodgsoniae) is an endemic species distributed in high-altitude areas of 3600-5600 m on the Qinghai-Tibet Plateau. To explore how the species is adapted to the high elevation environment, we assembled a draft genome based on both the Illumina and PacBio sequencing platforms with its population genetics and genomics analysis. In total, 134.74 Gb short reads and 30.81 Gb long reads raw data were generated. The 1.05-Gb assembled genome had a contig N50 of 4.56 Mb, with 91.94% complete BUSCOs. The 17 457 genes were annotated, and 11.35% of the genome was composed of repeat sequences. The phylogenetic tree showed that P. hodgsoniae was located at the basal position of the clade, including Golden Pheasant (Chrysolophus pictus), Common Pheasant (Phasianus colchicus), and Mikado Pheasant (Syrmaticus mikado). We found that 1014, 2595, and 2732 of the 6641 one-to-one orthologous genes were under positive selection in P. hodgsoniae, detected using PAML, BUSTED, and aBSREL programs, respectively, of which 965 genes were common under positive selection with 3 different programs. Several positively selected genes and immunity pathways relevant to high-altitude adaptation were detected. Gene family evolution showed that 99 gene families experienced significant expansion events, while 6 gene families were under contraction. The total number of olfactory receptor genes was relatively low in P. hodgsoniae. Genomic data provide an important resource for a further study on the evolutionary history of P. hodgsoniae, which provides a new insight into its high-altitude adaptation mechanisms.

DeepGWAS: Enhance GWAS Signals for Neuropsychiatric Disorders via Deep Neural Network

Genetic dissection of neuropsychiatric disorders can potentially reveal novel therapeutic targets. While genome-wide association studies (GWAS) have tremendously advanced our understanding, we approach a sample size bottleneck (i.e., the number of cases needed to identify >90% of all loci is impractical). Therefore, computationally enhancing GWAS on existing samples may be particularly valuable. Here, we describe DeepGWAS, a deep neural network-based method to enhance GWAS by integrating GWAS results with linkage disequilibrium and brain-related functional annotations. DeepGWAS enhanced schizophrenia (SCZ) loci by ~3X when applied to the largest European GWAS, and 21.3% enhanced loci were validated by the latest multi-ancestry GWAS. Importantly, DeepGWAS models can be transferred to other neuropsychiatric disorders. Transferring SCZ-trained models to Alzheimer's disease and major depressive disorder, we observed 1.3-17.6X detected loci compared to standard GWAS, among which 27-40% were validated by other GWAS studies. We anticipate DeepGWAS to be a powerful tool in GWAS studies.