Comparative genomics reveals insights into avian genome evolution and adaptation

Genome-wide analysis reveals the contributors to fast molecular evolution of the Chinese hook snout carp (Opsariichthys bidens)

Variations in molecular evolutionary rate have been widely investigated among lineages and genes. However, it remains an open question whether fast rate of molecular evolution is driven by natural selection or random drift, and how the fast rate is linked to metabolic rate. Additionally, previous studies on fast molecular evolution have been largely restricted to concatenated matrix of genes or a few specifically selected genes, but less is known for individual genes at the genome-wide level. Here we addressed these questions using more than 5000 single-copy orthologous (SCO) genes through comparative genomic and phylogenetic analyses among fishes, with a special focus on a newly-sequenced clupeocephalan fish the Chinese hook snout carp Opsariichthys bidens. We showed O. bidens displays significantly higher mean substitution rate and more fast-evolving SCO genes (2172 genes) than most fishes studied here. The rapidly evolving genes are enriched in highly conserved and very basic functions such as translation and ribosome that are critical for biological fitness. We further revealed that ∼25 % of these fast-evolving genes exhibit a constant increase of substitution rate from the common ancestor down to the present, suggesting a neglected but important contribution from ancestral states. Model fitting showed that ∼85 % of fast-evolving genes exclusive to O. bidens and related species follow the adaptive evolutionary model rather than random-drift model, and 7.6 % of fast-evolving genes identified in O. bidens have experienced positive selection, both indicating the reflection of adaptive selection. Finally, metabolic rate was observed to be linked with substitution rate in a gene-specific manner. Overall, our findings reveal fast molecular evolution of SCO genes at genome-wide level in O. bidens, and uncover the evolutionary and ecological contributors to it.

Application of multi-omics techniques to androgenetic alopecia: Current status and perspectives

The rapid advancement of sequencing technologies has enabled the generation of vast datasets, allowing for the in-depth analysis of sequencing data. This analysis has facilitated the validation of novel pathogenesis hypotheses for understanding and treating diseases through ex vivo and in vivo experiments. Androgenetic alopecia (AGA), a common hair loss disorder, has been a key focus of investigators attempting to uncover its underlying mechanisms. Abnormal changes in mRNA, proteins, and metabolites have been identified in individuals with AGA, and future developments in sequencing technologies may reveal new biomarkers for AGA. By integrating multiple omics analysis datasets such as genomics, transcriptomics, proteomics, and metabolomics-along with clinical phenotype data-we can achieve a comprehensive understanding of the molecular underpinnings of AGA. This review summarizes the data-mining studies conducted on various omics analysis datasets as related to AGA that have been adopted to interpret the biological data obtained from different omics layers. We herein discuss the challenges of integrative omics analyses, and suggest that collaborative multi-omics studies can enhance the understanding of the complete pathomechanism(s) of AGA by focusing on the interaction networks comprising DNA, RNA, proteins, and metabolites.

A reference genome for the Harpy Eagle reveals steady demographic decline and chromosomal rearrangements in the origin of Accipitriformes

The Harpy Eagle (Harpia harpyja) is an iconic species that inhabits forested landscapes in Neotropical regions, with decreasing population trends mainly due to habitat loss, and currently classified as vulnerable. Here, we report on a chromosome-scale genome assembly for a female individual combining long reads, optical mapping, and chromatin conformation capture reads. The final assembly spans 1.35 Gb, with N50scaffold equal to 58.1 Mb and BUSCO completeness of 99.7%. We built the first extensive transposable element (TE) library for the Accipitridae to date and identified 7,228 intact TEs. We found a burst of an unknown TE ~ 13-22 million years ago (MYA), coincident with the split of the Harpy Eagle from other Harpiinae eagles. We also report a burst of solo-LTRs and CR1 retrotransposons ~ 31-33 MYA, overlapping with the split of the ancestor to all Harpiinae from other Accipitridae subfamilies. Comparative genomics with other Accipitridae, the closely related Cathartidae and Galloanserae revealed major chromosome-level rearrangements at the basal Accipitriformes genome, in contrast to a conserved ancient genome architecture for the latter two groups. A historical demography reconstruction showed a rapid decline in effective population size over the last 20,000 years. This reference genome serves as a crucial resource for future conservation efforts towards the Harpy Eagle.

Reticulate evolution: Detection and utility in the phylogenomics era

Phylogenomics has enriched our understanding that the Tree of Life can have network-like or reticulate structures among some taxa and genes. Two non-vertical modes of evolution - hybridization/introgression and horizontal gene transfer - deviate from a strictly bifurcating tree model, causing non-treelike patterns. However, these reticulate processes can produce similar patterns to incomplete lineage sorting or recombination, potentially leading to ambiguity. Here, we present a brief overview of a phylogenomic workflow for inferring organismal histories and compare methods for distinguishing modes of reticulate evolution. We discuss how the timing of coalescent events can help disentangle introgression from incomplete lineage sorting and how horizontal gene transfer events can help determine the relative timing of speciation events. In doing so, we identify pitfalls of certain methods and discuss how to extend their utility across the Tree of Life. Workflows, methods, and future directions discussed herein underscore the need to embrace reticulate evolutionary patterns for understanding the timing and rates of evolutionary events, providing a clearer view of life's history.

Reconstructive evolutionary morphology: Tracing the historical process of modifications of complex systems driven by natural selection through changing ecological conditions

There is general consensus among evolutionary biologists that natural selection drives phenotypic modifications within populations over generational time. How to reconstruct this historical process, however, has been discussed mostly in theoretical terms, and recommendations and explanations on how to translate such theoretical insights into practice are needed. The present study aims at providing a theory-supported practical guide on how to reconstruct historical evolutionary processes by applying a morphology-centered approach through a series of interdependent steps of descriptive morphology, functional analysis, ecological observation, integration of paleoecological data, and evolutionary synthesis. Special attention is given to the development of tests regarding the accuracy, closeness to reality, and plausibility of the hypotheses at every level of the reconstructive process. This morphology-centered approach had its beginnings in the wake of the evolutionary synthesis and is part of the scientifically necessary process of reciprocal testing of hypotheses generated by different methods and data for the reconstruction of evolutionary history.

A reference genome, mitochondrial genome and associated transcriptomes for the critically endangered swift parrot ( Lathamus discolor)

The swift parrot ( Lathamus discolor) is a Critically Endangered migratory parrot that breeds in Tasmania and winters on the Australian mainland. Here we provide a reference genome assembly for the swift parrot. We sequence PacBio HiFi reads to create a high-quality reference assembly and identify a complete mitochondrial sequence. We also generate a reference transcriptome from five organs to inform genome annotation. The genome was 1.24 Gb in length and consisted of 847 contigs with a contig N50 of 18.97 Gb and L50 of 20 contigs. This study provides an annotated reference assembly and transcriptomic resources for the swift parrot to assist in future conservation genomic research.

Benefit-Sharing by Design: A Call to Action for Human Genomics Research

The ethical standards for the responsible conduct of human research have come a long way; however, concerns surrounding equity remain in human genetics and genomics research. Addressing these concerns will help society realize the full potential of human genomics research. One outstanding concern is the fair and equitable sharing of benefits from research on human participants. Several international bodies have recognized that benefit-sharing can be an effective tool for ethical research conduct, but international laws, including the Convention on Biological Diversity and its Nagoya Protocol on Access and Benefit-Sharing, explicitly exclude human genetic and genomic resources. These agreements face significant challenges that must be considered and anticipated if similar principles are applied in human genomics research. We propose that benefit-sharing from human genomics research can be a bottom-up effort and embedded into the existing research process. We propose the development of a "benefit-sharing by design" framework to address concerns of fairness and equity in the use of human genomic resources and samples and to learn from the aspirations and decade of implementation of the Nagoya Protocol.

High dispersal ability versus migratory traditions: Fine-scale population structure and post-glacial colonisation in bar-tailed godwits

In migratory animals, high mobility may reduce population structure through increased dispersal and enable adaptive responses to environmental change, whereas rigid migratory routines predict low dispersal, increased structure, and limited flexibility to respond to change. We explore the global population structure and phylogeographic history of the bar-tailed godwit, Limosa lapponica, a migratory shorebird known for making the longest non-stop flights of any landbird. Using nextRAD sequencing of 14,318 single-nucleotide polymorphisms and scenario-testing in an Approximate Bayesian Computation framework, we infer that bar-tailed godwits existed in two main lineages at the last glacial maximum, when much of their present-day breeding range persisted in a vast, unglaciated Siberian-Beringian refugium, followed by admixture of these lineages in the eastern Palearctic. Subsequently, population structure developed at both longitudinal extremes: in the east, a genetic cline exists across latitude in the Alaska breeding range of subspecies L. l. baueri; in the west, one lineage diversified into three extant subspecies L. l. lapponica, taymyrensis, and yamalensis, the former two of which migrate through previously glaciated western Europe. In the global range of this long-distance migrant, we found evidence of both (1) fidelity to rigid behavioural routines promoting fine-scale geographic population structure (in the east) and (2) flexibility to colonise recently available migratory flyways and non-breeding areas (in the west). Our results suggest that cultural traditions in highly mobile vertebrates can override the expected effects of high dispersal ability on population structure, and provide insights for the evolution and flexibility of some of the world's longest migrations.

Phylogenetic Conflict Between Species Tree and Maternally Inherited Gene Trees in a Clade of Emberiza Buntings (Aves: Emberizidae)

Different genomic regions may reflect conflicting phylogenetic topologies primarily due to incomplete lineage sorting and/or gene flow. Genomic data are necessary to reconstruct the true species tree and explore potential causes of phylogenetic conflict. Here, we investigate the phylogenetic relationships of 4 Emberiza species (Aves: Emberizidae) and discuss the potential causes of the observed mitochondrial non-monophyly of Emberiza godlewskii (Godlewski's bunting) using phylogenomic analyses based on whole genome resequencing data from 41 birds. Analyses based on both the whole mitochondrial genome and ~39 kilobases from the non-recombining W chromosome reveal sister relationships between each the northern and southern populations of E. godlewskii with E. cioides and E. cia, respectively. In contrast, the monophyly of E. godlewskii is reflected by the phylogenetic signal of autosomal and Z chromosomal sequence data as well as demographic inference analyses, which-in combination-support the following tree topology: ([{E. godlewskii, E. cia}, E. cioides], E. jankowskii). Using D-statistics, we detected multiple gene flow events among different lineages, indicating pervasive introgressive hybridization within this clade. Introgression from an unsampled lineage that is sister to E. cioides or introgression from an unsampled mitochondrial + W chromosomal lineage of E. cioides into northern E. godlewskii may explain the phylogenetic conflict between the species tree estimated from genome-wide data versus mtDNA/W tree topologies. These results underscore the importance of using genomic data for phylogenetic reconstruction and species delimitation.

Hybridization in birds-of-paradise: Widespread ancestral gene flow despite strong sexual selection in a lek-mating system

Sexual selection can directly contribute to reproductive isolation and is an important mechanism that can lead to speciation. Lek-mating is one of the most extreme forms of sexual selection, but surprisingly does not seem to preclude occasional hybridization in nature. However, hybridization among lekking species may still be trivial if selection against offspring with intermediate phenotypes prohibits introgression. Here we investigate this further by sequencing the genomes of nearly all bird-of-paradise (Paradisaeidae) species and 10 museum specimens of putative hybrid origin. We find that intergeneric hybridization indeed still takes place despite extreme differentiation in form, plumage, and behavior. In parallel, the genomes of contemporary species contain widespread signatures of past introgression, demonstrating that hybridization has repeatedly resulted in shared genetic variation despite strong sexual isolation. Our study raises important questions about extrinsic factors that modulate hybridization probability and the evolutionary consequences of introgressive hybridization between lekking species.

Chromosomal-level genome assembly and single-nucleotide polymorphism sites of black-faced spoonbill Platalea minor

Platalea minor, or black-faced spoonbill (Threskiornithidae), is a wading bird confined to coastal areas in East Asia. Due to habitat destruction, it was classified as globally endangered by the International Union for Conservation of Nature. However, the lack of genomic resources for this species hinders the understanding of its biology and diversity, and the development of conservation measures. Here, we report the first chromosomal-level genome assembly of P. minor using a combination of PacBio SMRT and Omni-C scaffolding technologies. The assembled genome (1.24 Gb) contains 95.33% of the sequences anchored to 31 pseudomolecules. The genome assembly has high sequence continuity with scaffold length N50 = 53 Mb. We predicted 18,780 protein-coding genes and measured high BUSCO score completeness (97.3%). Finally, we revealed 6,155,417 bi-allelic single nucleotide polymorphisms, accounting for ∼5% of the genome. This resource offers new opportunities for studying the black-faced spoonbill and developing conservation measures for this species.

The great tit HapMap project: A continental-scale analysis of genomic variation in a songbird

A major aim of evolutionary biology is to understand why patterns of genomic diversity vary within taxa and space. Large-scale genomic studies of widespread species are useful for studying how environment and demography shape patterns of genomic divergence. Here, we describe one of the most geographically comprehensive surveys of genomic variation in a wild vertebrate to date; the great tit (Parus major) HapMap project. We screened ca 500,000 SNP markers across 647 individuals from 29 populations, spanning ~30 degrees of latitude and 40 degrees of longitude - almost the entire geographical range of the European subspecies. Genome-wide variation was consistent with a recent colonisation across Europe from a South-East European refugium, with bottlenecks and reduced genetic diversity in island populations. Differentiation across the genome was highly heterogeneous, with clear 'islands of differentiation', even among populations with very low levels of genome-wide differentiation. Low local recombination rates were a strong predictor of high local genomic differentiation (FST), especially in island and peripheral mainland populations, suggesting that the interplay between genetic drift and recombination causes highly heterogeneous differentiation landscapes. We also detected genomic outlier regions that were confined to one or more peripheral great tit populations, probably as a result of recent directional selection at the species' range edges. Haplotype-based measures of selection were related to recombination rate, albeit less strongly, and highlighted population-specific sweeps that likely resulted from positive selection. Our study highlights how comprehensive screens of genomic variation in wild organisms can provide unique insights into spatio-temporal evolutionary dynamics.

Estrogen promotes gonadotropin-releasing hormone expression by regulating tachykinin 3 and prodynorphin systems in chicken

The "KNDy neurons" located in the hypothalamic arcuate nucleus (ARC) of mammals are known to co-express kisspeptin, neurokinin B (NKB), and dynorphin (DYN), and have been identified as key mediators of the feedback regulation of steroid hormones on gonadotropin-releasing hormone (GnRH). However, in birds, the genes encoding kisspeptin and its receptor GPR54 are genomic lost, leaving unclear mechanisms for feedback regulation of GnRH by steroid hormones. Here, the genes tachykinin 3 (TAC3) and prodynorphin (PDYN) encoding chicken NKB and DYN neuropeptides were successfully cloned. Temporal expression profiling indicated that TAC3, PDYN and their receptor genes (TACR3, OPRK1) were mainly expressed in the hypothalamus, with significantly higher expression at 30W than at 15W. Furthermore, overexpression or interference of TAC3 and PDYN can regulate the GnRH mRNA expression. In addition, in vivo and in vitro assays showed that estrogen (E2) could promote the mRNA expression of TAC3, PDYN, and GnRH, as well as the secretion of GnRH/LH. Mechanistically, E2 could dimerize the nuclear estrogen receptor 1 (ESR1) to regulate the expression of TAC3 and PDYN, which promoted the mRNA and protein expression of GnRH gene as well as the secretion of GnRH. In conclusion, these results revealed that E2 could regulate the GnRH expression through TAC3 and PDYN systems, providing novel insights for reproductive regulation in chickens.

Understanding microchromosomal organization and evolution in four representative woodpeckers (Picidae, Piciformes) through BAC-FISH analysis

The genome organization of woodpeckers has several distinctive features e.g., an uncommon accumulation of repetitive sequences, enlarged Z chromosomes, and atypical diploid numbers. Despite the large diversity of species, there is a paucity of detailed cytogenomic studies for this group and we thus aimed to rectify this. Genome organization patterns and hence evolutionary change in the microchromosome formation of four species (Colaptes campestris, Veniliornis spilogaster, Melanerpes candidus, and Picumnus nebulosus) was established through fluorescence in situ hybridization using bacterial artificial chromosomes originally derived from Gallus gallus and Taeniopygia guttata. Findings suggest that P. nebulosus (2n = 110), which was described for the first time, had the most basal karyotype among species of Picidae studied here, and probably arose as a result of fissions of avian ancestral macrochromosomes. We defined a new chromosomal number for V. spilogaster (2n = 88) and demonstrated microchromosomal rearrangements involving C. campestris plus a single, unique hitherto undescribed rearrangement in V. spilogaster. This comprised an inversion after a fusion involving the ancestral microchromosome 12 (homologous to chicken microchromosome 12). We also determined that the low diploid number of M. candidus is related to microchromosome fusions. Woodpeckers thus exhibit significantly rearranged karyotypes compared to the putative ancestral karyotype.

Genome Assembly of the Dyeing Poison Frog Provides Insights into the Dynamics of Transposable Element and Genome-Size Evolution

Genome size varies greatly across the tree of life and transposable elements are an important contributor to this variation. Among vertebrates, amphibians display the greatest variation in genome size, making them ideal models to explore the causes and consequences of genome size variation. However, high-quality genome assemblies for amphibians have, until recently, been rare. Here, we generate a high-quality genome assembly for the dyeing poison frog, Dendrobates tinctorius. We compare this assembly to publicly available frog genomes and find evidence for both large-scale conserved synteny and widespread rearrangements between frog lineages. Comparing conserved orthologs annotated in these genomes revealed a strong correlation between genome size and gene size. To explore the cause of gene-size variation, we quantified the location of transposable elements relative to gene features and find that the accumulation of transposable elements in introns has played an important role in the evolution of gene size in D. tinctorius, while estimates of insertion times suggest that many insertion events are recent and species-specific. Finally, we carry out population-scale mobile-element sequencing and show that the diversity and abundance of transposable elements in poison frog genomes can complicate genotyping from repetitive element sequence anchors. Our results show that transposable elements have clearly played an important role in the evolution of large genome size in D. tinctorius. Future studies are needed to fully understand the dynamics of transposable element evolution and to optimize primer or bait design for cost-effective population-level genotyping in species with large, repetitive genomes.

How structural variants shape avian phenotypes: Lessons from model systems

Despite receiving significant recent attention, the relevance of structural variation (SV) in driving phenotypic diversity remains understudied, although recent advances in long-read sequencing, bioinformatics and pangenomic approaches have enhanced SV detection. We review the role of SVs in shaping phenotypes in avian model systems, and identify some general patterns in SV type, length and their associated traits. We found that most of the avian SVs so far identified are short indels in chickens, which are frequently associated with changes in body weight and plumage colouration. Overall, we found that relatively short SVs are more frequently detected, likely due to a combination of their prevalence compared to large SVs, and a detection bias, stemming primarily from the widespread use of short-read sequencing and associated analytical methods. SVs most commonly involve non-coding regions, especially introns, and when patterns of inheritance were reported, SVs associated primarily with dominant discrete traits. We summarise several examples of phenotypic convergence across different species, mediated by different SVs in the same or different genes and different types of changes in the same gene that can lead to various phenotypes. Complex rearrangements and supergenes, which can simultaneously affect and link several genes, tend to have pleiotropic phenotypic effects. Additionally, SVs commonly co-occur with single-nucleotide polymorphisms, highlighting the need to consider all types of genetic changes to understand the basis of phenotypic traits. We end by summarising expectations for when long-read technologies become commonly implemented in non-model birds, likely leading to an increase in SV discovery and characterisation. The growing interest in this subject suggests an increase in our understanding of the phenotypic effects of SVs in upcoming years.

slc26a12-A novel member of the slc26 family, is located in tandem with slc26a2 in coelacanths, amphibians, reptiles, and birds

Solute carrier family 26 (Slc26) is a family of anion exchangers with 11 members in mammals (named Slc26a1-a11). Here, we identified a novel member of the slc26 family, slc26a12, located in tandem with slc26a2 in the genomes of several vertebrate lineages. BLAST and synteny analyses of various jawed vertebrate genome databases revealed that slc26a12 is present in coelacanths, amphibians, reptiles, and birds but not in cartilaginous fishes, lungfish, mammals, or ray-finned fishes. In some avian and reptilian lineages such as owls, penguins, egrets, and ducks, and most turtles examined, slc26a12 was lost or pseudogenized. Phylogenetic analysis showed that Slc26a12 formed an independent branch with the other Slc26 members and Slc26a12, Slc26a1 and Slc26a2 formed a single branch, suggesting that these three members formed a subfamily in Slc26. In jawless fish, hagfish have two genes homologous to slc26a2 and slc26a12, whereas lamprey has a single gene homologous to slc26a2. African clawed frogs express slc26a12 in larval gills, skin, and fins. These results show that slc26a12 was present at least before the separation of lobe-finned fish and tetrapods; the name slc26a12 is appropriate because the gene duplication occurred in the distant past.

Fast Bayesian Inference of Phylogenies from Multiple Continuous Characters

Time-scaled phylogenetic trees are an ultimate goal of evolutionary biology and a necessary ingredient in comparative studies. The accumulation of genomic data has resolved the tree of life to a great extent, yet timing evolutionary events remain challenging if not impossible without external information such as fossil ages and morphological characters. Methods for incorporating morphology in tree estimation have lagged behind their molecular counterparts, especially in the case of continuous characters. Despite recent advances, such tools are still direly needed as we approach the limits of what molecules can teach us. Here, we implement a suite of state-of-the-art methods for leveraging continuous morphology in phylogenetics, and by conducting extensive simulation studies we thoroughly validate and explore our methods' properties. While retaining model generality and scalability, we make it possible to estimate absolute and relative divergence times from multiple continuous characters while accounting for uncertainty. We compile and analyze one of the most data-type diverse data sets to date, comprised of contemporaneous and ancient molecular sequences, and discrete and continuous morphological characters from living and extinct Carnivora taxa. We conclude by synthesizing lessons about our method's behavior, and suggest future research venues.

The Limits of the Constant-rate Birth-Death Prior for Phylogenetic Tree Topology Inference

Birth-death models are stochastic processes describing speciation and extinction through time and across taxa and are widely used in biology for inference of evolutionary timescales. Previous research has highlighted how the expected trees under the constant-rate birth-death (crBD) model tend to differ from empirical trees, for example, with respect to the amount of phylogenetic imbalance. However, our understanding of how trees differ between the crBD model and the signal in empirical data remains incomplete. In this Point of View, we aim to expose the degree to which the crBD model differs from empirically inferred phylogenies and test the limits of the model in practice. Using a wide range of topology indices to compare crBD expectations against a comprehensive dataset of 1189 empirically estimated trees, we confirm that crBD model trees frequently differ topologically compared with empirical trees. To place this in the context of standard practice in the field, we conducted a meta-analysis for a subset of the empirical studies. When comparing studies that used Bayesian methods and crBD priors with those that used other non-crBD priors and non-Bayesian methods (i.e., maximum likelihood methods), we do not find any significant differences in tree topology inferences. To scrutinize this finding for the case of highly imbalanced trees, we selected the 100 trees with the greatest imbalance from our dataset, simulated sequence data for these tree topologies under various evolutionary rates, and re-inferred the trees under maximum likelihood and using the crBD model in a Bayesian setting. We find that when the substitution rate is low, the crBD prior results in overly balanced trees, but the tendency is negligible when substitution rates are sufficiently high. Overall, our findings demonstrate the general robustness of crBD priors across a broad range of phylogenetic inference scenarios but also highlight that empirically observed phylogenetic imbalance is highly improbable under the crBD model, leading to systematic bias in data sets with limited information content.

Geogenomic Predictors of Genetree Heterogeneity Explain Phylogeographic and Introgression History: A Case Study in an Amazonian Bird (Thamnophilus aethiops)

Can knowledge about genome architecture inform biogeographic and phylogenetic inference? Selection, drift, recombination, and gene flow interact to produce a genomic landscape of divergence wherein patterns of differentiation and genealogy vary nonrandomly across the genomes of diverging populations. For instance, genealogical patterns that arise due to gene flow should be more likely to occur on smaller chromosomes, which experience high recombination, whereas those tracking histories of geographic isolation (reduced gene flow caused by a barrier) and divergence should be more likely to occur on larger and sex chromosomes. In Amazonia, populations of many bird species diverge and introgress across rivers, resulting in reticulated genomic signals. Herein, we used reduced representation genomic data to disentangle the evolutionary history of 4 populations of an Amazonian antbird, Thamnophilus aethiops, whose biogeographic history was associated with the dynamic evolution of the Madeira River Basin. Specifically, we evaluate whether a large river capture event ca. 200 Ka, gave rise to reticulated genealogies in the genome by making spatially explicit predictions about isolation and gene flow based on knowledge about genomic processes. We first estimated chromosome-level phylogenies and recovered 2 primary topologies across the genome. The first topology (T1) was most consistent with predictions about population divergence and was recovered for the Z-chromosome. The second (T2), was consistent with predictions about gene flow upon secondary contact. To evaluate support for these topologies, we trained a convolutional neural network to classify our data into alternative diversification models and estimate demographic parameters. The best-fit model was concordant with T1 and included gene flow between non-sister taxa. Finally, we modeled levels of divergence and introgression as functions of chromosome length and found that smaller chromosomes experienced higher gene flow. Given that (1) genetrees supporting T2 were more likely to occur on smaller chromosomes and (2) we found lower levels of introgression on larger chromosomes (and especially the Z-chromosome), we argue that T1 represents the history of population divergence across rivers and T2 the history of secondary contact due to barrier loss. Our results suggest that a significant portion of genomic heterogeneity arises due to extrinsic biogeographic processes such as river capture interacting with intrinsic processes associated with genome architecture. Future phylogeographic studies would benefit from accounting for genomic processes, as different parts of the genome reveal contrasting, albeit complementary histories, all of which are relevant for disentangling the intricate geogenomic mechanisms of biotic diversification. [Amazonia; biogeography; demographic modeling; gene flow; gene tree; genome architecture; geogenomics; introgression; linked selection; neural network; phylogenomic; phylogeography; reproductive isolation; speciation; species tree.].

A nuclear genome assembly of an extinct flightless bird, the little bush moa

We present a draft genome of the little bush moa (Anomalopteryx didiformis)-one of approximately nine species of extinct flightless birds from Aotearoa, New Zealand-using ancient DNA recovered from a fossil bone from the South Island. We recover a complete mitochondrial genome at 249.9× depth of coverage and almost 900 megabases of a male moa nuclear genome at ~4 to 5× coverage, with sequence contiguity sufficient to identify more than 85% of avian universal single-copy orthologs. We describe a diverse landscape of transposable elements and satellite repeats, estimate a long-term effective population size of ~240,000, identify a diverse suite of olfactory receptor genes and an opsin repertoire with sensitivity in the ultraviolet range, show that the wingless moa phenotype is likely not attributable to gene loss or pseudogenization, and identify potential function-altering coding sequence variants in moa that could be synthesized for future functional assays. This genomic resource should support further studies of avian evolution and morphological divergence.

Metagenomic comparison of gut communities between wild and captive Himalayan griffons

Introduction

Himalayan griffons (Gyps himalayensis), known as the scavenger of nature, are large scavenging raptors widely distributed on the Qinghai-Tibetan Plateau and play an important role in maintaining the balance of the plateau ecosystem. The gut microbiome is essential for host health, helping to maintain homeostasis, improving digestive efficiency, and promoting the development of the immune system. Changes in environment and diet can affect the composition and function of gut microbiota, ultimately impacting the host health and adaptation. Captive rearing is considered to be a way to protect Himalayan griffons and increase their population size. However, the effects of captivity on the structure and function of the gut microbial communities of Himalayan griffons are poorly understood. Still, availability of sequenced metagenomes and functional information for most griffons gut microbes remains limited.

Methods

In this study, metagenome sequencing was used to analyze the composition and functional structures of the gut microbiota of Himalayan griffons under wild and captive conditions.

Results

Our results showed no significant differences in the alpha diversity between the two groups, but significant differences in beta diversity. Taxonomic classification revealed that the most abundant phyla in the gut of Himalayan griffons were Fusobacteriota, Proteobacteria, Firmicutes_A, Bacteroidota, Firmicutes, Actinobacteriota, and Campylobacterota. At the functional level, a series of Kyoto Encyclopedia of Genes and Genome (KEGG) functional pathways, carbohydrate-active enzymes (CAZymes) categories, virulence factor genes (VFGs), and pathogen-host interactions (PHI) were annotated and compared between the two groups. In addition, we recovered nearly 130 metagenome-assembled genomes (MAGs).

Discussion

In summary, the present study provided a first inventory of the microbial genes and metagenome-assembled genomes related to the Himalayan griffons, marking a crucial first step toward a wider investigation of the scavengers microbiomes with the ultimate goal to contribute to the conservation and management strategies for this near threatened bird.

Chromosome-level genome assembly of the European green woodpecker Picus viridis

The European green woodpecker, Picus viridis, is a widely distributed species found in the Western Palearctic region. Here, we assembled a highly contiguous genome assembly for this species using a combination of short- and long-read sequencing and scaffolded with chromatin conformation capture (Hi-C). The final genome assembly was 1.28 Gb and features a scaffold N50 of 37 Mb and a scaffold L50 of 39.165 Mb. The assembly incorporates 89.4% of the genes identified in birds in OrthoDB. Gene and repetitive content annotation on the assembly detected 15,805 genes and a ∼30.1% occurrence of repetitive elements, respectively. Analysis of synteny demonstrates the fragmented nature of the P. viridis genome when compared to the chicken (Gallus gallus). The assembly and annotations produced in this study will certainly help for further research into the genomics of P. viridis and the comparative evolution of woodpeckers. Five historical and seven contemporary samples have been resequenced and may give insights on the population history of this species.

Genome Assembly of Pyrocephalus nanus: A Step Toward the Genetic Conservation of the Endangered Little Vermilion Flycatcher of the Galapagos Islands

Incredibly powerful whole genome studies of conservation genetics, evolution, and biogeography become possible for non-model organisms when reference genomes are available. Here, we report the sequence and assembly of the whole genome of the little vermilion flycatcher (Pyrocephalus nanus; family Tyrannidae), which is an endemic, endangered, and declining species of the Galapagos Islands. Using PacBio HiFi reads to assemble long contigs and Hi-C reads for scaffolding, we assembled a genome of 1.07 Gb comprising 267 contigs in 152 scaffolds, scaffold N50 74 M, contig N50 17.8 M, with 98.9% assigned to candidate chromosomal sequences and 99.72% of the BUSCO passeriformes 10,844 single-copy orthologs present. In addition, we used the novel HiFiMiTie pipeline to fully assemble and verify all portions of the mitochondrial genome from HiFi reads, obtaining a mitogenome of 17,151 bases, containing 13 protein-coding genes, 22 tRNAs, 2 rRNAs, two control regions, and a unique structure of control region duplication and repeats. These genomes will be a critical tool for much-needed studies of phylogenetics, population genetics, biogeography, and conservation genetics of Pyrocephalus and related genera. This genome and other studies that use it will be able to provide recommendations for conservation management, taxonomic improvement, and to understand the evolution and diversification of this genus within the Galapagos Islands.

Comprehensive genome assembly reveals genetic diversity and carcass consumption insights in critically endangered Asian king vultures

The Asian king vulture (AKV), a vital forest scavenger, is facing globally critical endangerment. This study aimed to construct a reference genome to unveil the mechanisms underlying its scavenger abilities and to assess the genetic relatedness of the captive population in Thailand. A reference genome of a female AKV was assembled from sequencing reads obtained from both PacBio long-read and MGI short-read sequencing platforms. Comparative genomics with New World vultures (NWVs) and other birds in the Family Accipitridae revealed unique gene families in AKV associated with retroviral genome integration and feather keratin, contrasting with NWVs' genes related to olfactory reception. Expanded gene families in AKV were linked to inflammatory response, iron regulation and spermatogenesis. Positively selected genes included those associated with anti-apoptosis, immune response and muscle cell development, shedding light on adaptations for carcass consumption and high-altitude soaring. Using restriction site-associated DNA sequencing (RADseq)-based genome-wide single nucleotide polymorphisms (SNPs), genetic relatedness and inbreeding status of five captive AKVs were determined, revealing high genomic inbreeding in two females. In conclusion, the AKV reference genome was established, providing insights into its unique characteristics. Additionally, the potential of RADseq-based genome-wide SNPs for selecting AKV breeders was demonstrated.

A genomic basis of vocal rhythm in birds

Vocal rhythm plays a fundamental role in sexual selection and species recognition in birds, but little is known of its genetic basis due to the confounding effect of vocal learning in model systems. Uncovering its genetic basis could facilitate identifying genes potentially important in speciation. Here we investigate the genomic underpinnings of rhythm in vocal non-learning Pogoniulus tinkerbirds using 135 individual whole genomes distributed across a southern African hybrid zone. We find rhythm speed is associated with two genes that are also known to affect human speech, Neurexin-1 and Coenzyme Q8A. Models leveraging ancestry reveal these candidate loci also impact rhythmic stability, a trait linked with motor performance which is an indicator of quality. Character displacement in rhythmic stability suggests possible reinforcement against hybridization, supported by evidence of asymmetric assortative mating in the species producing faster, more stable rhythms. Because rhythm is omnipresent in animal communication, candidate genes identified here may shape vocal rhythm across birds and other vertebrates.

Adaptive expansion of ERVK solo-LTRs is associated with Passeriformes speciation events

Endogenous retroviruses (ERVs) are ancient retroviral remnants integrated in host genomes, and commonly deleted through unequal homologous recombination, leaving solitary long terminal repeats (solo-LTRs). This study, analysing the genomes of 362 bird species and their reptilian and mammalian outgroups, reveals an unusually higher level of solo-LTRs formation in birds, indicating evolutionary forces might have purged ERVs during evolution. Strikingly in the order Passeriformes, and especially the parvorder Passerida, endogenous retrovirus K (ERVK) solo-LTRs showed bursts of formation and recurrent accumulations coinciding with speciation events over past 22 million years. Moreover, our results indicate that the ongoing expansion of ERVK solo-LTRs in these bird species, marked by high transcriptional activity of ERVK retroviral genes in reproductive organs, caused variation of solo-LTRs between individual zebra finches. We experimentally demonstrated that cis-regulatory activity of recently evolved ERVK solo-LTRs may significantly increase the expression level of ITGA2 in the brain of zebra finches compared to chickens. These findings suggest that ERVK solo-LTRs expansion may introduce novel genomic sequences acting as cis-regulatory elements and contribute to adaptive evolution. Overall, our results underscore that the residual sequences of ancient retroviruses could influence the adaptive diversification of species by regulating host gene expression.

Different Host-Endogenous Retrovirus Relationships between Mammals and Birds Reflected in Genome-Wide Evolutionary Interaction Patterns

Mammals and birds differ largely in their average endogenous retrovirus loads, namely the proportion of endogenous retrovirus in the genome. The host-endogenous retrovirus relationships, including conflict and co-option, have been hypothesized among the causes of this difference. However, there has not been studies about the genomic evolutionary signal of constant host-endogenous retrovirus interactions in a long-term scale and how such interactions could lead to the endogenous retrovirus load difference. Through a phylogeny-controlled correlation analysis on ∼5,000 genes between the dN/dS ratio of each gene and the load of endogenous retrovirus in 12 mammals and 21 birds, separately, we detected genes that may have evolved in association with endogenous retrovirus loads. Birds have a higher proportion of genes with strong correlation between dN/dS and the endogenous retrovirus load than mammals. Strong evidence of association is found between the dN/dS of the coding gene for leucine-rich repeat-containing protein 23 and endogenous retrovirus load in birds. Gene set enrichment analysis shows that gene silencing rather than immunity and DNA recombination may have a larger contribution to the association between dN/dS and the endogenous retrovirus load for both mammals and birds. The above results together showing different evolutionary patterns between bird and mammal genes can partially explain the apparently lower endogenous retrovirus loads of birds, while gene silencing may be a universal mechanism that plays a remarkable role in the evolutionary interaction between the host and endogenous retrovirus. In summary, our study presents signals that the host genes might have driven or responded to endogenous retrovirus load changes in long-term evolution.

Species-aware DNA language models capture regulatory elements and their evolution

BACKGROUND

The rise of large-scale multi-species genome sequencing projects promises to shed new light on how genomes encode gene regulatory instructions. To this end, new algorithms are needed that can leverage conservation to capture regulatory elements while accounting for their evolution.

RESULTS

Here, we introduce species-aware DNA language models, which we trained on more than 800 species spanning over 500 million years of evolution. Investigating their ability to predict masked nucleotides from context, we show that DNA language models distinguish transcription factor and RNA-binding protein motifs from background non-coding sequence. Owing to their flexibility, DNA language models capture conserved regulatory elements over much further evolutionary distances than sequence alignment would allow. Remarkably, DNA language models reconstruct motif instances bound in vivo better than unbound ones and account for the evolution of motif sequences and their positional constraints, showing that these models capture functional high-order sequence and evolutionary context. We further show that species-aware training yields improved sequence representations for endogenous and MPRA-based gene expression prediction, as well as motif discovery.

CONCLUSIONS

Collectively, these results demonstrate that species-aware DNA language models are a powerful, flexible, and scalable tool to integrate information from large compendia of highly diverged genomes.

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.

Comparative Genomics Reveal Phylogenetic Relationship and Chromosomal Evolutionary Events of Eight Cervidae Species

Cervidae represents a family that is not only rich in species diversity but also exhibits a wide range of karyotypes. The controversies regarding the phylogeny and classification of Cervidae still persist. The flourishing development of the genomic era has made it possible to address these issues at the genomic level. Here, the genomes of nine species were used to explore the phylogeny and chromosomal evolutionary events of Cervidae. By conducting whole-genome comparisons, we identified single-copy orthologous genes across the nine species and constructed a phylogenetic tree based on the single-copy orthologous genes sequences, providing new insights into the phylogeny of Cervidae, particularly the phylogenetic relationship among sika deer, red deer, wapiti and Tarim red deer. Gene family analysis revealed contractions in the olfactory receptor gene family and expansions in the histone gene family across eight Cervidae species. Furthermore, synteny analysis was used to explore the chromosomal evolutionary events of Cervidae species, revealing six chromosomal fissions during the evolutionary process from Bovidae to Cervidae. Notably, specific chromosomal fusion events were found in four species of Cervus, and a unique chromosomal fusion event was identified in Muntiacus reevesi. Our study further completed the phylogenetic relationship within the Cervidae and demonstrated the feasibility of inferring species phylogeny at the whole-genome level. Additionally, our findings on gene family evolution and the chromosomal evolutionary events in eight Cervidae species lay a foundation for comprehensive research of the evolution of Cervidae.

Ancient and Modern Genomes Reveal Microsatellites Maintain a Dynamic Equilibrium Through Deep Time

Microsatellites are widely used in population genetics, but their evolutionary dynamics remain poorly understood. It is unclear whether microsatellite loci drift in length over time. This is important because the mutation processes that underlie these important genetic markers are central to the evolutionary models that employ microsatellites. We identify more than 27 million microsatellites using a novel and unique dataset of modern and ancient Adélie penguin genomes along with data from 63 published chordate genomes. We investigate microsatellite evolutionary dynamics over 2 timescales: one based on Adélie penguin samples dating to ∼46.5 ka and the other dating to the diversification of chordates aged more than 500 Ma. We show that the process of microsatellite allele length evolution is at dynamic equilibrium; while there is length polymorphism among individuals, the length distribution for a given locus remains stable. Many microsatellites persist over very long timescales, particularly in exons and regulatory sequences. These often retain length variability, suggesting that they may play a role in maintaining phenotypic variation within populations.

Genetic Basis and Evolution of Structural Color Polymorphism in an Australian Songbird

Island organisms often evolve phenotypes divergent from their mainland counterparts, providing a useful system for studying adaptation under differential selection. In the white-winged fairywren (Malurus leucopterus), subspecies on two islands have a black nuptial plumage whereas the subspecies on the Australian mainland has a blue nuptial plumage. The black subspecies have a feather nanostructure that could in principle produce a blue structural color, suggesting a blue ancestor. An earlier study proposed independent evolution of melanism on the islands based on the history of subspecies divergence. However, the genetic basis of melanism and the origin of color differentiation in this group are still unknown. Here, we used whole-genome resequencing to investigate the genetic basis of melanism by comparing the blue and black M. leucopterus subspecies to identify highly divergent genomic regions. We identified a well-known pigmentation gene ASIP and four candidate genes that may contribute to feather nanostructure development. Contrary to the prediction of convergent evolution of island melanism, we detected signatures of a selective sweep in genomic regions containing ASIP and SCUBE2 not in the black subspecies but in the blue subspecies, which possesses many derived SNPs in these regions, suggesting that the mainland subspecies has re-evolved a blue plumage from a black ancestor. This proposed re-evolution was likely driven by a preexisting female preference. Our findings provide new insight into the evolution of plumage coloration in island versus continental populations, and, importantly, we identify candidate genes that likely play roles in the development and evolution of feather structural coloration.

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.

Evolution and expression patterns of the neo-sex chromosomes of the crested ibis

Bird sex chromosomes play a unique role in sex-determination, and affect the sexual morphology and behavior of bird species. Core waterbirds, a major clade of birds, share the common characteristics of being sexually monomorphic and having lower levels of inter-sexual conflict, yet their sex chromosome evolution remains poorly understood. Here, by we analyse of a chromosome-level assembly of a female crested ibis (Nipponia nippon), a typical core waterbird. We identify neo-sex chromosomes resulting from fusion of microchromosomes with ancient sex chromosomes. These fusion events likely occurred following the divergence of Threskiornithidae and Ardeidae. The neo-W chromosome of the crested ibis exhibits the characteristics of slow degradation, which is reflected in its retention of abundant gametologous genes. Neo-W chromosome genes display an apparent ovary-biased gene expression, which is largely driven by genes that are retained on the crested ibis W chromosome but lost in other bird species. These results provide new insights into the evolutionary history and expression patterns for the sex chromosomes of bird species.

Rapid adaptation of cellular metabolic rate to the MicroRNA complements of mammals and its relevance to the evolution of endothermy

The metabolic efficiency of mammalian cells depends on the attenuation of intrinsic translation noise by microRNAs. We devised a metric of cellular metabolic rate (cMR), rMR/Mexp optimally fit to the number of microRNA families (mirFam), that is robust to variation in mass and sensitive to body temperature (Tb), consistent with the heat dissipation limit theory of Speakman and Król (2010). Using mirFam as predictor, an Ornstein-Uhlenbeck process of stabilizing selection, with an adaptive shift at the divergence of Boreoeutheria, accounted for 95% of the variation in cMR across mammals. Branchwise rates of evolution of cMR, mirFam and Tb concurrently increased 6- to 7-fold at the divergence of Boreoeutheria, independent of mass. Cellular MR variation across placental mammals was also predicted by the sum of model conserved microRNA-target interactions, revealing an unexpected degree of integration of the microRNA-target apparatus into the energy economy of the mammalian cell.

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.

De novo Phased Genome Assembly, Annotation and Population Genotyping of Alectoris Chukar

The Alectoris Chukar (chukar) is the most geographically widespread partridge species in the world, demonstrating exceptional adaptability to diverse ecological environments. However, the scarcity of genetic resources for chukar has hindered research into its adaptive evolution and molecular breeding. In this study, we have sequenced and assembled a high-quality, phased chukar genome that consists of 31 pairs of relatively complete diploid chromosomes. Our BUSCO analysis reported a high completeness score of 96.8% and 96.5%, with respect to universal single-copy orthologs and a low duplication rate (0.3% and 0.5%) for two assemblies. Through resequencing and population genomic analyses of six subspecies, we have curated invaluable genotype data that underscores the adaptive evolution of chukar in response to both arid and high-altitude environments. These data will significantly contribute to research on how chukars adaptively evolve to cope with desertification and alpine climates.

Genomic evolution of island birds from the view of the Swinhoe's pheasant (Lophura swinhoii)

Island endemic birds account for the majority of extinct vertebrates in the past few centuries. To date, the evolutionary characteristics of island endemic bird's is poorly known. In this research, we de novo assembled a high-quality chromosome-level reference genome for the Swinhoe's pheasant, which is a typical endemic island bird. Results of collinearity tests suggest rapid ancient chromosome rearrangement that may have contributed to the initial species radiation within Phasianidae, and a role for the insertions of CR1 transposable elements in rearranging chromosomes in Phasianidae. During the evolution of the Swinhoe's pheasant, natural selection positively selected genes involved in fecundity and body size functions, at both the species and population levels, which reflect genetic variation associated with island adaptation. We further tested for variation in population genomic traits between the Swinhoe's pheasant and its phylogenetically closely related mainland relative the silver pheasant, and found higher levels of genetic drift and inbreeding in the Swinhoe's pheasant genome. Divergent demographic histories of insular and mainland bird species during the last glacial period may reflect the differing impact of insular and continental climates on the evolution of species. Our research interprets the natural history and population genetic characteristics of the insular endemic bird the Swinhoe's pheasant, at a genome-wide scale, provides a broader perspective on insular speciation, and adaptive evolution and contributes to the genetic conservation of island endemic birds.

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.

Conserved chromatin and repetitive patterns reveal slow genome evolution in frogs

Frogs are an ecologically diverse and phylogenetically ancient group of anuran amphibians that include important vertebrate cell and developmental model systems, notably the genus Xenopus. Here we report a high-quality reference genome sequence for the western clawed frog, Xenopus tropicalis, along with draft chromosome-scale sequences of three distantly related emerging model frog species, Eleutherodactylus coqui, Engystomops pustulosus, and Hymenochirus boettgeri. Frog chromosomes have remained remarkably stable since the Mesozoic Era, with limited Robertsonian (i.e., arm-preserving) translocations and end-to-end fusions found among the smaller chromosomes. Conservation of synteny includes conservation of centromere locations, marked by centromeric tandem repeats associated with Cenp-a binding surrounded by pericentromeric LINE/L1 elements. This work explores the structure of chromosomes across frogs, using a dense meiotic linkage map for X. tropicalis and chromatin conformation capture (Hi-C) data for all species. Abundant satellite repeats occupy the unusually long (~20 megabase) terminal regions of each chromosome that coincide with high rates of recombination. Both embryonic and differentiated cells show reproducible associations of centromeric chromatin and of telomeres, reflecting a Rabl-like configuration. Our comparative analyses reveal 13 conserved ancestral anuran chromosomes from which contemporary frog genomes were constructed.

Fine-Scale Map Reveals Highly Variable Recombination Rates Associated with Genomic Features in the Eurasian Blackcap

Recombination is responsible for breaking up haplotypes, influencing genetic variability, and the efficacy of selection. Bird genomes lack the protein PR domain-containing protein 9, a key determinant of recombination dynamics in most metazoans. Historical recombination maps in birds show an apparent stasis in positioning recombination events. This highly conserved recombination pattern over long timescales may constrain the evolution of recombination in birds. At the same time, extensive variation in recombination rate is observed across the genome and between different species of birds. Here, we characterize the fine-scale historical recombination map of an iconic migratory songbird, the Eurasian blackcap (Sylvia atricapilla), using a linkage disequilibrium-based approach that accounts for population demography. Our results reveal variable recombination rates among and within chromosomes, which associate positively with nucleotide diversity and GC content and negatively with chromosome size. Recombination rates increased significantly at regulatory regions but not necessarily at gene bodies. CpG islands are associated strongly with recombination rates, though their specific position and local DNA methylation patterns likely influence this relationship. The association with retrotransposons varied according to specific family and location. Our results also provide evidence of heterogeneous intrachromosomal conservation of recombination maps between the blackcap and its closest sister taxon, the garden warbler. These findings highlight the considerable variability of recombination rates at different scales and the role of specific genomic features in shaping this variation. This study opens the possibility of further investigating the impact of recombination on specific population-genomic features.

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.

Web of venom: exploration of big data resources in animal toxin research

Research on animal venoms and their components spans multiple disciplines, including biology, biochemistry, bioinformatics, pharmacology, medicine, and more. Manipulating and analyzing the diverse array of data required for venom research can be challenging, and relevant tools and resources are often dispersed across different online platforms, making them less accessible to nonexperts. In this article, we address the multifaceted needs of the scientific community involved in venom and toxin-related research by identifying and discussing web resources, databases, and tools commonly used in this field. We have compiled these resources into a comprehensive table available on the VenomZone website (https://venomzone.expasy.org/10897). Furthermore, we highlight the challenges currently faced by researchers in accessing and using these resources and emphasize the importance of community-driven interdisciplinary approaches. We conclude by underscoring the significance of enhancing standards, promoting interoperability, and encouraging data and method sharing within the venom research community.

The genome of a globally invasive passerine, the common myna, Acridotheres tristis

In an era of global climate change, biodiversity conservation is receiving increased attention. Conservation efforts are greatly aided by genetic tools and approaches, which seek to understand patterns of genetic diversity and how they impact species health and their ability to persist under future climate regimes. Invasive species offer vital model systems in which to investigate questions regarding adaptive potential, with a particular focus on how changes in genetic diversity and effective population size interact with novel selection regimes. The common myna (Acridotheres tristis) is a globally invasive passerine and is an excellent model species for research both into the persistence of low-diversity populations and the mechanisms of biological invasion. To underpin research on the invasion genetics of this species, we present the genome assembly of the common myna. We describe the genomic landscape of this species, including genome wide allelic diversity, methylation, repeats, and recombination rate, as well as an examination of gene family evolution. Finally, we use demographic analysis to identify that some native regions underwent a dramatic population increase between the two most recent periods of glaciation, and reveal artefactual impacts of genetic bottlenecks on demographic analysis.

Comparison of Endogenous Alpharetroviruses (ALV-like) across Galliform Species: New Distant Proviruses

The Genus Alpharetrovirus contains viruses pathogenic mainly for chickens, forming the Avian Sarcoma and Leukosis Virus group (ASLV). Cells of most Galliform species, besides chickens, contain genetic elements (endogenous retroviruses, ERVs) that could recombine with other alpharetroviruses or express proteins, complementing defective ASLV, which may successfully replicate and cause disease. However, they are quite unknown, and only ALV-F, from ring-necked pheasants, has been partially published. Upon scrutiny of 53 genomes of different avian species, we found Alpharetrovirus-like sequences only in 12 different Galliformes, including six full-length (7.4-7.6 Kbp) and 27 partial sequences. Phylogenetic studies of the regions studied (LTR, gag, pol, and env) consistently resulted in five almost identical clades containing the same ERVs: Clade I (presently known ASLVs); Clade II (Callipepla spp. ERVs); Clade IIIa (Phasianus colchicus ERVs); Clade IIIb (Alectoris spp. ERVs); and Clade IV (Centrocercus spp. ERVs). The low pol identity scores suggested that each of these Clades may be considered a different species. ORF analysis revealed that putatively encoded proteins would be very similar in length and domains to those of other alpharetroviruses and thus potentially functional. This will undoubtedly contribute to better understanding the biology of defective viruses, especially in wild Galliformes, their evolution, and the danger they may represent for other wild species and the poultry industry.

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.

RERconverge Expansion: Using Relative Evolutionary Rates to Study Complex Categorical Trait Evolution

Comparative genomics approaches seek to associate evolutionary genetic changes with the evolution of phenotypes across a phylogeny. Many of these methods, including our evolutionary rates based method, RERconverge, lack the capability of analyzing non-ordinal, multicategorical traits. To address this limitation, we introduce an expansion to RERconverge that associates shifts in evolutionary rates with the convergent evolution of multi-categorical traits. The categorical RERconverge expansion includes methods for performing categorical ancestral state reconstruction, statistical tests for associating relative evolutionary rates with categorical variables, and a new method for performing phylogenetic permulations on multi-categorical traits. In addition to demonstrating our new method on a three-category diet phenotype, we compare its performance to naive pairwise binary RERconverge analyses and two existing methods for comparative genomic analyses of categorical traits: phylogenetic simulations and a phylogenetic signal based method. We also present a diagnostic analysis of the new permulations approach demonstrating how the method scales with the number of species and the number of categories included in the analysis. Our results show that our new categorical method outperforms phylogenetic simulations at identifying genes and enriched pathways significantly associated with the diet phenotype and that the new ancestral reconstruction drives an improvement in our ability to capture diet-related enriched pathways. Our categorical permulations were able to account for non-uniform null distributions and correct for non-independence in gene rank during pathway enrichment analysis. The categorical expansion to RERconverge will provide a strong foundation for applying the comparative method to categorical traits on larger data sets with more species and more complex trait evolution.

A chromosome-level genome assembly for the Silkie chicken resolves complete sequences for key chicken metabolic, reproductive, and immunity genes

A set of high-quality pan-genomes would help identify important genes that are still hidden/incomplete in bird reference genomes. In an attempt to address these issues, we have assembled a de novo chromosome-level reference genome of the Silkie (Gallus gallus domesticus), which is an important avian model for unique traits, like fibromelanosis, with unclear genetic foundation. This Silkie genome includes the complete genomic sequences of well-known, but unresolved, evolutionarily, endocrinologically, and immunologically important genes, including leptin, ovocleidin-17, and tumor-necrosis factor-α. The gap-less and manually annotated MHC (major histocompatibility complex) region possesses 38 recently identified genes, with differentially regulated genes recovered in response to pathogen challenges. We also provide whole-genome methylation and genetic variation maps, and resolve a complex genetic region that may contribute to fibromelanosis in these animals. Finally, we experimentally show leptin binding to the identified leptin receptor in chicken, confirming an active leptin ligand-receptor system. The Silkie genome assembly not only provides a rich data resource for avian genome studies, but also lays a foundation for further functional validation of resolved genes.

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.