Dense sampling of bird diversity increases power of comparative genomics

Whole-genome resequencing landscape of adaptive evolution in Relict gull (Larus relictus)

BACKGROUND

The relict gull (Larus relictus, Charadriiformes, Laridae) classified as vulnerable in the IUCN Red List is defined as a first-class national protected bird in China. However, our knowledge of the evolutionary history of L. relictus is limited. Here, we performed whole-genome resequencing of L. relictus (n = 14) and L. brunnicephalus (n = 3) to explore the genetic relationships and population structures and understand their adaptive evolution.

RESULTS

The whole genome resequencing generated 667.55 Gb clean reads with an average sequencing depth of ~ 29×. The genomic variant analysis identified 13,717,267 heterozygous SNPs in the samples obtained from 17 individuals. Population genetic diversity analysis revealed that low nucleotide diversity (0.00029) and no obvious population structure in L. relictus. Demographic history revealed that from 180 to 5 kya (thousand years ago), the effective population size (Ne) of L. relictus exhibited declines (24,000 to 5,000), with a very low range population size (2,200 to 5,000). In contrast, from 100 to 80 kya, L. brunnicephalus peaked in ancestral Ne, followed by distinct declines at ~ 70 kya (100,000 to 16,000). The findings identified several genes associated with the correlated changed life-history traits of L. relictus, including BMP4 involved in beak adaptation; HAND2, NEUROG1, COL11A2, and EDNRB involved in the evolution of the palate, soft palate, and tongue; PIGR and PLCB2 involved in an enhanced response to bitter taste by sensing chemical secretions released by staple food substrate insects to activate protective mechanisms. Furthermore, thirty-four genes related to sperm development and activity, including KLHL10 and TEKT3, were identified in the expanded gene family. In addition, MED1, CNOT9, NR5A1, and PATZ1 were involved in enhanced male hormone secretion and a high density of candidate genes associated with embryonic development were identified. The findings indicated that the L. relictus population was in a male-biased diffusion mode; the function of the TEKT3 gene showed that males played a dominant role in brooding, which enhanced their attraction to females. Our study revealed that significant enrichment of olfactory signaling pathway genes, including OR14C36, OR14J1, OR14I1, and OR14A16; inner ear development-related, including PTN, PTPN11, GATA2, ATP8B1, and MYO15A; and those related to hypoxic adaptation to high-altitude breeding and iris colour.

CONCLUSIONS

Based on the results and the knowledge of this organism biology and habitat use, we infer that less adaptive evolutionary pressure on vision in L. relictus were related with their feeding behaviour and adaptation. In summary, this comprehensive analysis provides insights into the evolutionary features of L. relictus and a new perspective for scientific research on L. relictus to effectively determine its future survival viability.

A phylogenetic approach to comparative genomics

Comparative genomics, whereby the genomes of different species are compared, has the potential to address broad and fundamental questions at the intersection of genetics and evolution. However, species, genomes and genes cannot be considered as independent data points within statistical tests. Closely related species tend to be similar because they share genes by common descent, which must be accounted for in analyses. This problem of non-independence may be exacerbated when examining genomes or genes but can be addressed by applying phylogeny-based methods to comparative genomic analyses. Here, we review how controlling for phylogeny can change the conclusions of comparative genomics studies. We address common questions on how to apply these methods and illustrate how they can be used to test causal hypotheses. The combination of rapidly expanding genomic datasets and phylogenetic comparative methods is set to revolutionize the biological insights possible from comparative genomic studies.

The complete mitogenomes of all four Cryptospiza species (Aves: Estrildidae)

Crimsonwings are estrildid finches found in the understory of montane rainforests of sub-Saharan Africa. The genus includes four species: Cryptospiza jacksoni Sharpe 1902, C. shelleyi Sharpe 1902, C. reichenovii (Hartlaub 1874), and C. salvadorii Reichenow 1892. The first two are endemic to the Albertine Rift, while the latter two are more widespread. Despite being well-represented in museum collections, genetic resources are scarce. Here we provide complete mitogenomes for all four species, each containing the standard 37 avian genes. Analyses showed C. shelleyi as sister to the other three species, with C. reichenovii and C. salvadorii being highly similar (99.2%). Further research is needed to explore their evolutionary history.

Thrushes in Love: Extensive Gene Flow, With Differential Resistance and Selection, Obscures and Reveals the Evolutionary History of a Songbird Clade

The application of high-throughput sequencing to phylogenetic analyses is allowing authors to reconstruct the true evolutionary history of species. This work can illuminate specific mechanisms underlying divergence when combined with analyses of gene flow, recombination and selection. We conducted a phylogenomic analysis of Catharus, a songbird genus with considerable potential for gene flow, variation in migratory behaviour and genomic resources. We documented discordance among trees constructed for mitochondrial, autosomal and sex (Z) chromosome partitions. Two trees were recovered on the Z. Both trees differed from the autosomes, one matched the mitochondria, and the other was unique to the Z. Gene flow with one species likely generated much of this discordance; substantial admixture between ustulatus and the remaining species was documented and linked to at least two historic events. The tree unique to the Z likely reflects the true history of Catharus; local genomic analyses recovered the same tree in autosomal regions with reduced admixture and recombination. Genes previously connected to migration were enriched in these regions suggesting transitions between migratory and non-migratory states helped generate divergence. Migratory (vs. nonmigratory) Catharus formed a monophyletic clade in a subset of genomic regions. Gene flow was elevated in some of these regions suggesting adaptive introgression may have occurred, but the dominant pattern was of balancing selection maintaining ancestral polymorphisms important for olfaction and perhaps, by extension, adaptation to temperate climates. This work illuminates the evolutionary history of an important model in speciation and demonstrates how differential resistance to gene flow can affect local genomic patterns.

Pangenome graphs and their applications in biodiversity genomics

Complete datasets of genetic variants are key to biodiversity genomic studies. Long-read sequencing technologies allow the routine assembly of highly contiguous, haplotype-resolved reference genomes. However, even when complete, reference genomes from a single individual may bias downstream analyses and fail to adequately represent genetic diversity within a population or species. Pangenome graphs assembled from aligned collections of high-quality genomes can overcome representation bias by integrating sequence information from multiple genomes from the same population, species or genus into a single reference. Here, we review the available tools and data structures to build, visualize and manipulate pangenome graphs while providing practical examples and discussing their applications in biodiversity and conservation genomics across the tree of life.

Quality assessment of gene repertoire annotations with OMArk

In the era of biodiversity genomics, it is crucial to ensure that annotations of protein-coding gene repertoires are accurate. State-of-the-art tools to assess genome annotations measure the completeness of a gene repertoire but are blind to other errors, such as gene overprediction or contamination. We introduce OMArk, a software package that relies on fast, alignment-free sequence comparisons between a query proteome and precomputed gene families across the tree of life. OMArk assesses not only the completeness but also the consistency of the gene repertoire as a whole relative to closely related species and reports likely contamination events. Analysis of 1,805 UniProt Eukaryotic Reference Proteomes with OMArk demonstrated strong evidence of contamination in 73 proteomes and identified error propagation in avian gene annotation resulting from the use of a fragmented zebra finch proteome as a reference. This study illustrates the importance of comparing and prioritizing proteomes based on their quality measures.

High-Throughput Shotgun Metagenomics of Microbial Footprints Uncovers a Cocktail of Noxious Antibiotic Resistance Genes in the Winam Gulf of Lake Victoria, Kenya

Background: A diverse range of pollutants, including heavy metals, agrochemicals, pharmaceutical residues, illicit drugs, personal care products, and other anthropogenic contaminants, pose a significant threat to aquatic ecosystems. The Winam Gulf of Lake Victoria, heavily impacted by surrounding human activities, faces potential contamination from these pollutants. However, studies exploring the presence of antibiotic resistance genes (ARGs) in the lake remain limited. In the current study, a shotgun metagenomics approach was employed to identify ARGs and related pathways. Genomic DNA was extracted from water and sediment samples and sequenced using the high-throughput Illumina NovaSeq platform. Additionally, phenotypic antibiotic resistance was assessed using the disk diffusion method with commonly used antibiotics. Results: The analysis of metagenomes sequences from the Gulf ecosystem and Comprehensive Antibiotic Resistance Database (CARD) revealed worrying levels of ARGs in the lake. The study reported nine ARGs from the 37 high-risk resistant gene families previously documented by the World Health Organization (WHO). Proteobacteria had the highest relative abundance of antibiotic resistance (53%), Bacteriodes (4%), Verrucomicrobia (2%), Planctomycetes Chloroflexi, Firmicutes (2%), and other unclassified bacteria (39%). Genes that target protection, replacement, change, and antibiotic-resistant efflux were listed in order of dominance. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis revealed antibiotic resistance to beta-lactamase and vancomycin. Phenotypic resistance to vancomycin, tetracycline, sulfamethoxazole, erythromycin, trimethoprim, tetracycline, and penicillin was reported through the zone of inhibition. Conclusions: This study highlights that the Winam Gulf of Lake Victoria in Kenya harbors a diverse array of antibiotic-resistant genes, including those conferring multidrug resistance. These findings suggest that the Gulf could be serving as a reservoir for more antibiotic-resistant genes, posing potential risks to both human health and aquatic biodiversity. The insights gained from this research can guide policy development for managing antibiotic resistance in Kenya.

Response strategies to acute and chronic environmental stress in the arctic breeding Lapland longspur (Calcarius lapponicus)

The potentially devastating effects of climate change have raised awareness of the need to understand how the biology of wild animals is influenced by extreme-weather events. We investigate how a wild arctic-breeding bird, the Lapland longspur (Calcarius lapponicus), responds to different environmental perturbations and its coping strategies. We explore the transcriptomic response to environmental adversity during the transition from arrival at the breeding grounds to incubation on the Arctic tundra. The effects of an extremely cold spring on arrival and a severe storm during incubation are examined through RNA-seq analysis of pertinent tissues sampled across the breeding cycle. The stress response, circadian rhythms, reproduction, and metabolism are all affected. A key gene of the Hypothalamic-Pituitary-Adrenal axis, FKBP5, was significantly up-regulated in hypothalamus. The genome assembly and gene expression profiles provide comprehensive resources for future studies. Our findings on different coping strategies to chronic and acute stressors will contribute to understanding the interplay between changing environments and genomic regulation.

Museomics Sheds Light on Evolutionary Diversity in a Critically Endangered Cockatoo Species From Wallacea

Accurate identification of evolutionarily significant units of rare and threatened organisms provides a foundation for effective management and conservation. Up to seven subspecies of the critically endangered Yellow-crested Cockatoo (Cacatua sulphurea) have been described, four of which were commonly recognised pre-2014. In the absence of genotypic data, C. sulphurea subspecies delimitation has been based on morphology, behaviour and biogeography. To clarify genetic relationships and shed light on the diversification of this parrot radiation, whole genomes were sequenced for 16 museum specimens, covering the geographic range of the proposed seven subspecies as well as one C. galerita galerita. Combined with four museum-derived wild Cacatua sequences from NCBI, the results indicate there are three distinct C. sulphurea subspecies clusters centred in different biogeographic subregions of Wallacea (Timor; Sumba; as well as the Sulawesi Region and the main Lesser Sunda chain), separated by shallow genetic distances (da < 0.148%). The results raise questions about the recent species-level elevation of the phenotypically most distinct subspecies, C. s. citrinocristata, and about the origins of C. s. abbotti, the only subspecies west of Wallace's Line. Our analyses suggest C. s. abbotti is unlikely to be embedded within C. sulphurea, suggesting its origin on the remote Masalembu islands may be due to human translocation via historical trade routes. These genomic results inform the prioritisation and streamlining of conservation measures for the critically endangered C. sulphurea by identifying and delimiting likely conservation units.

Sexual selection promotes reproductive isolation in barn swallows

Despite the well-known effects of sexual selection on phenotypes, links between this evolutionary process and reproductive isolation, genomic divergence, and speciation have been difficult to establish. We unravel the genetic basis of sexually selected plumage traits to investigate their effects on reproductive isolation in barn swallows. The genetic architecture of sexual traits is characterized by 12 loci on two autosomes and the Z chromosome. Sexual trait loci exhibit signatures of divergent selection in geographic isolation and barriers to gene flow in secondary contact. Linkage disequilibrium between these genes has been maintained by selection in hybrid zones beyond what would be expected under admixture alone. Our findings reveal that selection on coupled sexual trait loci promotes reproductive isolation, providing key empirical evidence for the role of sexual selection in speciation.

Untangling the colonization history of the Australo-Pacific reed warblers, one of the world's great island radiations

Island radiations, such as those of the Australo-Pacific, offer unique insight into diversification, extinction, and early speciation processes. Yet, their speciation and colonization histories are often obscured by conflicting genomic signals from incomplete lineage sorting (ILS) or hybridization. Here, we integrated mitogenomes and genome-wide SNPs to unravel the evolutionary history of one of the world's most geographically widespread island radiations. The Australo-Pacific reed warblers (Acrocephalus luscinius complex) are a speciose lineage including five species that have become extinct since the 19th century and ten additional species of conservation concern. The radiation spans over 10,000 km across Australo-Papua, Micronesia and Polynesia, including the Mariana, Hawaii and Pitcairn Island archipelagos. Earlier mtDNA studies suggested a stepping-stone colonization process, resulting in archipelago-level secondary sympatry of divergent mtDNA lineages in the Mariana Islands and Marquesas. These studies hypothesized that morphologically similar species on neighboring islands arose from ecological convergence. Using DNA from historical museum specimens and modern genetic samples, we show that ILS and/or gene flow have shaped the radiation of Australo-Pacific reed warblers rather than secondary sympatry. The nuclear genome reconstructs a simpler biogeographic history than mtDNA, showing close relationships between species in the Mariana Islands and Marquesas despite their paraphyletic mtDNA lineages. Gene flow likely involved early and late colonizing waves of the radiation before the loss of ancestral dispersive ability. Our results highlight how collection genomics can elucidate evolutionary history and inform conservation efforts for threatened species.

Avian Island Radiations Shed Light on the Dynamics of Adaptive and Nonadaptive Radiation

Understanding the mechanisms underlying species formation and differentiation is a central goal of evolutionary biology and a formidable challenge. This understanding can provide valuable insights into the origins of the astonishing diversity of organisms living on our planet. Avian evolutionary radiations on islands have long fascinated biologists as they provide the ideal variation to study the ecological and evolutionary forces operating on the continuum between incipient lineages to complete speciation. In this review, we summarize the key insights gained from decades of research on adaptive and nonadaptive radiations of both extant and extinct insular bird species. We present a new comprehensive global list of potential avian radiations on oceanic islands, based on published island species checklists, taxonomic studies, and phylogenetic analyses. We demonstrate that our understanding of evolutionary processes is being greatly enhanced through the use of genomic tools. However, to advance the field, it is critical to complement this information with a solid understanding of the ecological and behavioral traits of both extinct and extant avian island species.

Genomic patterns in the dwarf kingfishers of northern Melanesia reveal a mechanistic framework explaining the paradox of the great speciators

The paradox of the great speciators describes a contradictory biogeographic pattern exhibited by numerous avian lineages in Oceania. Specifically, these lineages display broad geographic distributions across the region, implying strong over-water dispersal capabilities; yet, they also display repeated genetic and phenotypic divergence-even between geographically proximate islands-implying poor inter-island dispersal capabilities. One group originally cited as evidence for this paradox is the dwarf kingfishers of the genus Ceyx. Here, using genomic sequencing and comprehensive geographic sampling of the monophyletic Ceyx radiation from northern Melanesia, we find repeated, deep genetic divergence and no evidence for gene flow between lineages found on geographically proximate islands, providing an exceptionally clear example of the paradox of the great speciators. A dated phylogenetic reconstruction suggests a significant burst of diversification occurred rapidly after reaching northern Melanesia, between 3.9 and 2.9 MYA. This pattern supports a shift in net diversification rate, concordant with the expectations of the "colonization cycle" hypothesis, which implies a historical shift in dispersiveness among great speciator lineages during the evolutionary past. Here, we present a formalized framework that explains how repeated founder effects and shifting selection pressures on highly dispersive genotypes are the only ultimate causes needed to generate the paradox of the great speciators. Within this framework, we emphasize that lineage-specific traits and island-specific abiotic factors will result in varying levels of selection pressure against dispersiveness, caused by varying proximate eco-evolutionary mechanisms. Overall, we highlight how understanding patterns of diversification in the Ceyx dwarf kingfishers helped us generate a cohesive framework that provides a rigorous mechanistic explanation for patterns concordant with the paradox of the great speciators and the repeated emergence of geographic radiations in island archipelagoes across the globe.

A high-quality genome assembly of the Spectacled Fulvetta (Fulvetta ruficapilla) endemic to China

The Spectacled Fulvetta (Fulvetta ruficapilla) is the type species of Fulvetta, an evolutionarily distinct group whose species show a high degree of sympatry in distribution and phenotypic convergence. To pave the way for insights into their adaptive evolution and speciation, we have assembled the first high quality reference genome for F. ruficapilla using high-fidelity (HiFi) long-read and Hi-C sequencing technologies. The resulting assembly spans a total of ~1.21 Gb with a contig N50 of 18.8 Mb and scaffold N50 of 75.9 Mb, and has a BUSCO completeness of 97.0%. The quality assessment suggests a high standard in base accuracy, continuity, and completeness of the assembly, comparable or close to that of Vertebrate Genomes Project. On this basis, we have annotated 23,774 protein-coding genes, of which 18,832 are functionally identified. The availability of this high-quality genome provides a solid foundation for the future studies of evolution and local adaptation in birds.

Comparison of retinol binding protein 1 with cone specific G-protein as putative effector molecules in cryptochrome signalling

Vision and magnetoreception in navigating songbirds are strongly connected as recent findings link a light dependent radical-pair mechanism in cryptochrome proteins to signalling pathways in cone photoreceptor cells. A previous yeast-two-hybrid screening approach identified six putative candidate proteins showing binding to cryptochrome type 4a. So far, only the interaction of the cone specific G-protein transducin α-subunit was investigated in more detail. In the present study, we compare the binding features of the G-protein α-subunit with those of another candidate from the yeast-two-hybrid screen, cellular retinol binding protein. Purified recombinant European robin retinol binding protein bound retinol with high affinity, displaying an EC50 of less than 5 nM, thereby demonstrating its functional state. We applied surface plasmon resonance and a Förster resonance transfer analysis to test for interactions between retinol binding protein and cryptochrome 4a. In the absence of retinol, we observed no robust binding events, which contrasts the strong interaction we observed between cryptochrome 4a and the G-protein α-subunit. We conclude that retinol binding protein is unlikely to be involved in the primary magnetosensory signalling cascade.

Repeated evolution on oceanic islands: comparative genomics reveals species-specific processes in birds

Understanding the interplay between genetic drift, natural selection, gene flow, and demographic history in driving phenotypic and genomic differentiation of insular populations can help us gain insight into the speciation process. Comparing patterns across different insular taxa subjected to similar selective pressures upon colonizing oceanic islands provides the opportunity to study repeated evolution and identify shared patterns in their genomic landscapes of differentiation. We selected four species of passerine birds (Common Chaffinch Fringilla coelebs/canariensis, Red-billed Chough Pyrrhocorax pyrrhocorax, House Finch  Haemorhous mexicanus and Dark-eyed/island Junco Junco hyemalis/insularis) that have both mainland and insular populations. Changes in body size between island and mainland populations were consistent with the island rule. For each species, we sequenced whole genomes from mainland and insular individuals to infer their demographic history, characterize their genomic differentiation, and identify the factors shaping them. We estimated the relative (Fst) and absolute (dxy) differentiation, nucleotide diversity (π), Tajima's D, gene density and recombination rate. We also searched for selective sweeps and chromosomal inversions along the genome. All species shared a marked reduction in effective population size (Ne) upon island colonization. We found diverse patterns of differentiated genomic regions relative to the genome average in all four species, suggesting the role of selection in island-mainland differentiation, yet the lack of congruence in the location of these regions indicates that each species evolved differently in insular environments. Our results suggest that the genomic mechanisms involved in the divergence upon island colonization-such as chromosomal inversions, and historical factors like recurrent selection-differ in each species, despite the highly conserved structure of avian genomes and the similar selective factors involved. These differences are likely influenced by factors such as genetic drift, the polygenic nature of fitness traits and the action of case-specific selective pressures.

Cellular, Molecular, and Genetic Mechanisms of Avian Beak Development and Evolution

Diverse research programs employing complementary strategies have been uncovering cellular, molecular, and genetic mechanisms essential to avian beak development and evolution. In reviewing these discoveries, I offer an interdisciplinary perspective on bird beaks that spans their derivation from jaws of dinosaurian reptiles, their anatomical and ecological diversification across major taxonomic groups, their common embryonic origins, their intrinsic patterning processes, and their structural integration. I describe how descriptive and experimental approaches, including gene expression and cell lineage analyses, tissue recombinations, surgical transplants, gain- and loss-of-function methods, geometric morphometrics, comparative genomics, and genome-wide association studies, have identified key constituent parts and putative genes regulating beak morphogenesis and evolution. I focus throughout on neural crest mesenchyme, which generates the beak skeleton and other components, and describe how these embryonic progenitor cells mediate species-specific pattern and link form and function as revealed by 20 years of research using chimeras between quail and duck embryos.

Beyond genes-for-behaviour: The potential for genomics to resolve long-standing questions in avian brood parasitism

Behavioural ecology by definition of its founding 'Tinbergian framework' is an integrative field, however, it lags behind in incorporating genomic methods. 'Finding the gene/s for a behaviour' is still rarely feasible or cost-effective in the wild but as we show here, genomic data can be used to address broader questions. Here we use avian brood parasitism, a model system in behavioural ecology as a case study to highlight how behavioural ecologists could use the full potential of state-of-the-art genomic tools. Brood parasite-host interactions are one of the most easily observable and amenable natural laboratories of antagonistic coevolution, and as such have intrigued evolutionary biologists for decades. Using worked examples, we demonstrate how genomic data can be used to study the causes and mechanisms of (co)evolutionary adaptation and answer three key questions for the field: (i) Where and when should brood parasitism evolve?, (ii) When and how should hosts defend?, and (iii) Will coevolution persist with ecological change? In doing so, we discuss how behavioural and molecular ecologists can collaborate to integrate Tinbergen's questions and achieve the coherent science that he promoted to solve the mysteries of nature.

Genomic and Acoustic Biogeography of the Iconic Sulphur-crested Cockatoo Clarifies Species Limits and Patterns of Intraspecific Diversity

Many highly recognizable species lack genetic data important for conservation due to neglect over their hyperabundance. This likely applies to the Sulfur-crested Cockatoo (Cacatua galerita), one of the world's most iconic parrots. The species is native to Australia, New Guinea, and some surrounding Melanesian islands of the latter. Four subspecies are currently recognised based on morphology. Australian subspecies and populations are abundant, but several factors threaten those in New Guinea and Melanesia. Genetic data from natural populations are scarce-information that is vital to identifying evolutionarily significant units (ESUs) important for modern conservation planning. We used whole-genome resequencing to investigate patterns of differentiation, evolutionary affinities, and demographic history across C. galerita's distribution range to assess whether currently recognised subspecies represent ESUs. We complement this with an assessment of bioacoustic variation across the species' distribution landscape. Our results point to C. galerita sensu lato (s.l.) comprising two species. We restrict C. galerita sensu stricto (s.s.) to populations in Australia and the Trans-Fly ecodomain of southern New Guinea. The second species, recognised here as Cacatua triton, likely occurs over much of the rest of New Guinea. Restricting further discussion of intraspecific diversity in C. triton, we show that within C. galerita s.s. two ESUs exist, which align to Cacatua galerita galerita in eastern Australia and southern New Guinea and Cacatua galerita fitzroyi in northern and north-western Australia. We suggest that the evolution of these species and ESUs are linked to Middle and Late Pleistocene glacial cycles and their effects on sea level and preferential habitats. We argue that conservation assessments need updating, protection of preferential forest and woodland habitats are important and reintroductions require careful management to avoid possible negative hybridization effects of non-complementary lineages.

Insights into avian molecular cytogenetics-with reptilian comparisons

In last 100 years or so, much information has been accumulated on avian karyology, genetics, physiology, biochemistry and evolution. The chicken genome project generated genomic resources used in comparative studies, elucidating fundamental evolutionary processes, much of it funded by the economic importance of domestic fowl (which are also excellent model species in many areas). Studying karyotypes and whole genome sequences revealed population processes, evolutionary biology, and genome function, uncovering the role of repetitive sequences, transposable elements and gene family expansion. Knowledge of the function of many genes and non-expressed or identified regulatory components is however still lacking. Birds (Aves) are diverse, have striking adaptations for flight, migration and survival and inhabit all continents most islands. They also have a unique karyotype with ~ 10 macrochromosomes and ~ 30 microchromosomes that are smaller than other reptiles. Classified into Palaeognathae and Neognathae they are evolutionarily close, and a subset of reptiles. Here we overview avian molecular cytogenetics with reptilian comparisons, shedding light on their karyotypes and genome structure features. We consider avian evolution, then avian (followed by reptilian) karyotypes and genomic features. We consider synteny disruptions, centromere repositioning, and repetitive elements before turning to comparative avian and reptilian genomics. In this context, we review comparative cytogenetics and genome mapping in birds as well as Z- and W-chromosomes and sex determination. Finally, we give examples of pivotal research areas in avian and reptilian cytogenomics, particularly physical mapping and map integration of sex chromosomal genes, comparative genomics of chicken, turkey and zebra finch, California condor cytogenomics as well as some peculiar cytogenetic and evolutionary examples. We conclude that comparative molecular studies and improving resources continually contribute to new approaches in population biology, developmental biology, physiology, disease ecology, systematics, evolution and phylogenetic systematics orientation. This also produces genetic mapping information for chromosomes active in rearrangements during the course of evolution. Further insights into mutation, selection and adaptation of vertebrate genomes will benefit from these studies including physical and online resources for the further elaboration of comparative genomics approaches for many fundamental biological questions.

The systematic codon usage bias has an important effect on genetic adaption in native species

Random mutations increase genetic variety and natural selection enhances adaption over generations. Codon usage biases (CUB) provide clues about the genome adaptation mechanisms of native species and extremophile species. Significant numbers of gene (CDS) of nine classes of endangered, native species, including extremophiles and mesophiles were utilised to compute CUB. Codon usage patterns differ among the lineages of endangered and extremophiles with native species. Polymorphic usage of nucleotides with codon burial suggests parallelism of native species within relatively confined taxonomic groups. Utilizing the deviation pattern of CUB of endangered and native species, I present a calculation parameter to estimate the extinction risk of endangered species. Species diversity and extinction risk are both positively associated with the propensity of random mutation in CDS (Coding DNA sequence). Codon bias tenet profoundly selected and it governs to adaptive evolution of native species.

Visualizing genomic evolution in Caenorhabditis through WormSynteny

Understanding the syntenic relationships among genomes is crucial to elucidate the genomic mechanisms that drive the evolution of species. The nematode Caenorhabditis is a good model for studying genomic evolution due to the well-established biology of Caenorhabditis elegans and the availability of > 50 genomes in the genus. However, effective alignment of more than ten species in Caenorhabditis has not been conducted before, and there is currently no tool to visualize the synteny of more than two species. In this study, we used Progressive Cactus, a recently developed multigenome aligner, to align the genomes of eleven Caenorhabditis species. Through the progressive alignment, we reconstructed nine ancestral genomes, analyzed the mutational types that cause genomic rearrangement during speciation, and found that insertion and duplication are the major driving forces for genome expansion. Dioecious species appear to expand their genomes more than androdioecious species. We then built an online interactive app called WormSynteny to visualize the syntenic relationship among the eleven species. Users can search the alignment dataset using C. elegans query sequences, construct synteny plots at different genomic scales, and use a set of options to control alignment output and plot presentation. We showcased the use of WormSynteny to visualize the syntenic conservation of one-to-one orthologues among species, tandem and dispersed gene duplication in C. elegans, and the evolution of exon and intron structures. Importantly, the integration of orthogroup information with synteny linkage in WormSynteny allows the easy visualization of conserved genomic blocks and disruptive rearrangement. In conclusion, WormSynteny provides immediate access to the syntenic relationships among the most widely used Caenorhabditis species and can facilitate numerous comparative genomics studies. This pilot study with eleven species also serves as a proof-of-concept to a more comprehensive larger-scale analysis using hundreds of nematode genomes, which is expected to reveal mechanisms that drive genomic evolution in the Nematoda phylum. Finally, the WormSynteny software provides a generalizable solution for visualizing the output of Progressive Cactus with interactive graphics, which would be useful for a broad community of genome researchers.

Mitogenomics clarifies the position of the Nearctic magpies (Pica hudsonia and Pica nuttalli) within the Holarctic magpie radiation

Partial separation of a peripheral population may lead to its divergence and, potentially, speciation due to genetic drift followed by selection and geographic isolation. This process may cause taxonomic uncertainty because reproductive isolation in allopatry cannot be verified directly. The two Nearctic allopatric species of magpies (Aves, Corvidae: Pica) serve as a good example of these problems. The Black-billed magpie Pica hudsonia is widely distributed in North America, whereas the Yellow-billed Magpie Pica nuttalli is endemic to a restricted range in California. Their relationships with Palearctic species have been little studied. We obtained complete mitochondrial genomes of both Nearctic magpie species, along with the Eurasian Magpie (Pica pica) and the Oriental Magpie (Pica serica), 20 mitogenomes in total. Phylogenetic analysis reveals a basal position of P. serica, and P. pica as a sister clade to the two Nearctic species. P. hudsonia and P. nuttalli form reciprocal monophyletic subclades, showing recent divergence between and within them. Our data show that the Nearctic magpie lineage diverged from the common ancestor with P. pica, with a single migration wave via the Beringia. Within the Nearctic, we hypothesize a peripatric mode of speciation among Pica taxa due to the divergence and separation of the small marginal population in California below the Sierra-Nevada mountains. Diversifying amino acid substitutions in ND4-ND5-ND6 genes along the branch leading to the New World clade may indicate selection for heat-tolerance. Considering the clear phenotypic differences between P. hudsonia and P. nuttalli, our data, showing their reciprocal monophylies and genetic distinctness, is consistent with the two-species taxonomy.

Long-read sequencing and genome assembly of natural history collection samples and challenging specimens

Museum collections harbor millions of samples, largely unutilized for long-read sequencing. Here, we use ethanol-preserved samples containing kilobase-sized DNA to show that amplification-free protocols can yield contiguous genome assemblies. Additionally, using a modified amplification-based protocol, employing an alternative polymerase to overcome PCR bias, we assembled the 3.1 Gb maned sloth genome, surpassing the previous 500 Mb protocol size limit. Our protocol also improves assemblies of other difficult-to-sequence molluscs and arthropods, including millimeter-sized organisms. By highlighting collections as valuable sample resources and facilitating genome assembly of tiny and challenging organisms, our study advances efforts to obtain reference genomes of all eukaryotes.

A new high-quality genome assembly and annotation for the threatened Florida Scrub-Jay (Aphelocoma coerulescens)

The Florida Scrub-Jay (Aphelocoma coerulescens), a Federally Threatened, cooperatively-breeding bird, is an emerging model system in evolutionary biology and ecology. Extensive individual-based monitoring and genetic sampling for decades has yielded a wealth of data, allowing for the detailed study of social behavior, demography, and population genetics of this natural population. Here, we report a linkage map and a chromosome-level genome assembly and annotation for a female Florida Scrub-Jay made with long-read sequencing technology, chromatin conformation data, and the linkage map. We constructed a linkage map comprising 4,468 SNPs that had 34 linkage groups and a total sex-averaged autosomal genetic map length of 2446.78 cM. The new genome assembly is 1.33 Gb in length, consisting of 33 complete or near-complete autosomes and the sex chromosomes (ZW). This highly contiguous assembly has an NG50 of 68 Mb and a Benchmarking Universal Single-Copy Orthologs (BUSCO) completeness score of 97.1% with respect to the Aves database. The annotated gene set has a BUSCO transcriptome completeness score of 95.5% and 17,964 identified protein-coding genes, 92.5% of which have associated functional annotations. This new, high-quality genome assembly and linkage map of the Florida Scrub-Jay provides valuable tools for future research into the evolutionary dynamics of small, natural populations of conservation concern.

Sturnidae sensu lato Mitogenomics: Novel Insights into Codon Aversion, Selection, and Phylogeny

The Sturnidae family comprises 123 recognized species in 35 genera. The taxa Mimidae and Buphagidae were formerly treated as subfamilies within Sturnidae. The phylogenetic relationships among the Sturnidae and related taxa (Sturnidae sensu lato) remain unresolved due to high rates of morphological change and concomitant morphological homoplasy. This study presents five new mitogenomes of Sturnidae sensu lato and comprehensive mitogenomic analyses. The investigated mitogenomes exhibit an identical gene composition of 37 genes-including 13 protein-coding genes (PCGs), 2 rRNA genes, and 22 tRNA genes-and one control region (CR). The most important finding of this study is drawn from CAM analyses. The surprisingly unique motifs for each species provide a new direction for the molecular species identification of avian. Furthermore, the pervasiveness of the natural selection of PCGs is found in all examined species when analyzing their nucleotide composition and codon usage. We also determine the structures of mt-tRNA, mt-rRNA, and CR structures of Sturnidae sensu lato. Lastly, our phylogenetic analyses not only well support the monophyly of Sturnidae, Mimidae, and Buphagidae, but also define nine stable subclades. Taken together, our findings will enable the further elucidation of the evolutionary relationships within Sturnidae sensu lato.

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.

The complete mitochondrial genome of aglaeactis castelnaudii (bourcier & mulsant, 1848) (apodiformes: trochilidae: aglaeactis) and phylogenetic analysis

In this study, we employed high-throughput sequencing data to assemble the mitochondrial genome (mitogenome) of the White-tufted Sunbeam (Aglaeactis castelnaudii). The total length of the mitogenome was found to be 16,872 base pairs (bp), containing 13 protein-coding genes (PCGs), 22 transfer RNA genes, 2 ribosomal RNA genes, and 1 control region. The nucleotide composition was as follows: A 30.6%, T 24.0%, C 31.2%, and G 14.2%, resulting in a GC content of 45.4%. Phylogenetic analysis, utilizing the concatenation of the 13 mitochondrial PCGs, indicated a closer evolutionary relationship between the genus Aglaeactis and the genus Coeligena compared to other genera within the family Trochilidae investigated in this study. The mitogenome of A. castelnaudii not only contributes to species identification but also provides valuable insights for phylogenetic and conservation genetic analyses of A. castelnaudii.

Chromosome mapping of retrotransposon AviRTE in a neotropical bird species: Trogon surrucura (Trogoniformes; Trogonidae)

Avian genomes are characterized as being more compact than other amniotes, with less diversity and density of transposable elements (TEs). In addition, birds usually show bimodal karyotypes, exhibiting a great variation in diploid numbers. Some species present unusually large sex chromosomes, possibly due to the accumulation of repetitive sequences. Avian retrotransposon-like element (AviRTE) is a long interspersed nuclear element (LINE) recently discovered in the genomes of birds and nematodes, and it is still poorly characterized in terms of chromosomal mapping and phylogenetic relationships. In this study, we mapped AviRTE isolated from the Trogon surrucura genome into the T. surrucura (TSU) karyotype. Furthermore, we analyzed the phylogenetic relationships of this LINE in birds and other vertebrates. Our results showed that the distribution pattern of AviRTE is not restricted to heterochromatic regions, with accumulation on the W chromosome of TSU, yet another species with an atypical sex chromosome and TE hybridization. The phylogenetic analysis of AviRTE sequences in birds agreed with the proposed phylogeny of species in most clades, and allowed the detection of this sequence in other species, expanding the distribution of the element.

Spatiotemporal omics for biology and medicine

The completion of the Human Genome Project has provided a foundational blueprint for understanding human life. Nonetheless, understanding the intricate mechanisms through which our genetic blueprint is involved in disease or orchestrates development across temporal and spatial dimensions remains a profound scientific challenge. Recent breakthroughs in cellular omics technologies have paved new pathways for understanding the regulation of genomic elements and the relationship between gene expression, cellular functions, and cell fate determination. The advent of spatial omics technologies, encompassing both imaging and sequencing-based methodologies, has enabled a comprehensive understanding of biological processes from a cellular ecosystem perspective. This review offers an updated overview of how spatial omics has advanced our understanding of the translation of genetic information into cellular heterogeneity and tissue structural organization and their dynamic changes over time. It emphasizes the discovery of various biological phenomena, related to organ functionality, embryogenesis, species evolution, and the pathogenesis of diseases.

Genome and life-history evolution link bird diversification to the end-Cretaceous mass extinction

Complex patterns of genome evolution associated with the end-Cretaceous [Cretaceous-Paleogene (K-Pg)] mass extinction limit our understanding of the early evolutionary history of modern birds. Here, we analyzed patterns of avian molecular evolution and identified distinct macroevolutionary regimes across exons, introns, untranslated regions, and mitochondrial genomes. Bird clades originating near the K-Pg boundary exhibited numerous shifts in the mode of molecular evolution, suggesting a burst of genomic heterogeneity at this point in Earth's history. These inferred shifts in substitution patterns were closely related to evolutionary shifts in developmental mode, adult body mass, and patterns of metabolic scaling. Our results suggest that the end-Cretaceous mass extinction triggered integrated patterns of evolution across avian genomes, physiology, and life history near the dawn of the modern bird radiation.

Museum genomics approach to study the taxonomy and evolution of Woolly-necked storks using historic specimens

The accessibility of genomic tools in evolutionary biology has allowed for a thorough exploration of various evolutionary processes associated with adaptation and speciation. However, genomic studies in natural systems present numerous challenges, reflecting the inherent complexities of studying organisms in their native habitats. The utilization of museum specimens for genomics research has received increased attention in recent times, facilitated by advancements in ancient DNA techniques. In this study, we have utilized a museum genomics approach to analyze historic specimens of Woolly-necked storks (Ciconia spp.) and examine their genetic composition and taxonomic status and explore the evolutionary and adaptive trajectories of populations over the years. The Woolly-necked storks are distributed in Asia and Africa with a taxonomic classification that has been a matter of ambiguity. Asian and African Woollynecks were recently recognized as different species based on their morphological differences; however, their genomic validation was lacking. In this study, we have used ∼70-year-old museum samples for whole-genome population-scale sequencing. Our study has revealed that Asian and African Woollynecks are genetically distinct, consistent with the current taxonomic classification based on morphological features. However, we also found a high genetic divergence between the Asian subspecies Ciconia episcopus neglecta and Ciconia episcopus episcopus, suggesting this classification requires a detailed examination to explore processes of ongoing speciation. Because taxonomic classification directly impacts conservation efforts, and there is evidence of declining populations of Asian Woollynecks in Southeast Asia, our results highlight that population-scale studies are urgent to determine the genetic, ecological, and phylogenetic diversity of these birds.

Unanticipated Loss of Inflammasomes in Birds

Inflammasomes are multiprotein complexes that form in response to ligands originating from pathogens as well as alterations of normal cell physiology caused by infection or tissue damage. These structures engage a robust inflammatory immune response that eradicates environmental microbes before they cause disease, and slow the growth of bona fide pathogens. Despite their undeniable utility in immunity, inflammasomes are radically reduced in birds. Perhaps most surprising is that, within all birds, NLRP3 is retained, while its signaling adapter ASC is lost, suggesting that NLRP3 signals via a novel unknown adapter. Crocodilian reptiles and turtles, which share a more recent common ancestor with birds, retain many of the lost inflammasome components, indicating that the deletion of inflammasomes occurred after birds diverged from crocodiles. Some bird lineages have even more extensive inflammasome loss, with songbirds continuing to pare down their inflammasomes until only NLRP3 and CARD8 remain. Remarkably, songbirds have lost caspase-1 but retain the downstream targets of caspase-1: IL-1β, IL-18, and the YVAD-linker encoding gasdermin A. This suggests that inflammasomes can signal through alternative proteases to activate cytokine maturation and pyroptosis in songbirds. These observations may reveal new contexts of activation that may be relevant to mammalian inflammasomes and may suggest new avenues of research to uncover the enigmatic nature of the poorly understood NLRP3 inflammasome.

Please Mind the Gap: Indel-Aware Parsimony for Fast and Accurate Ancestral Sequence Reconstruction and Multiple Sequence Alignment Including Long Indels

Despite having important biological implications, insertion, and deletion (indel) events are often disregarded or mishandled during phylogenetic inference. In multiple sequence alignment, indels are represented as gaps and are estimated without considering the distinct evolutionary history of insertions and deletions. Consequently, indels are usually excluded from subsequent inference steps, such as ancestral sequence reconstruction and phylogenetic tree search. Here, we introduce indel-aware parsimony (indelMaP), a novel way to treat gaps under the parsimony criterion by considering insertions and deletions as separate evolutionary events and accounting for long indels. By identifying the precise location of an evolutionary event on the tree, we can separate overlapping indel events and use affine gap penalties for long indel modeling. Our indel-aware approach harnesses the phylogenetic signal from indels, including them into all inference stages. Validation and comparison to state-of-the-art inference tools on simulated data show that indelMaP is most suitable for densely sampled datasets with closely to moderately related sequences, where it can reach alignment quality comparable to probabilistic methods and accurately infer ancestral sequences, including indel patterns. Due to its remarkable speed, our method is well suited for epidemiological datasets, eliminating the need for downsampling and enabling the exploitation of the additional information provided by dense taxonomic sampling. Moreover, indelMaP offers new insights into the indel patterns of biologically significant sequences and advances our understanding of genetic variability by considering gaps as crucial evolutionary signals rather than mere artefacts.

Integrated genomics provides insights into the evolution of the polyphosphate accumulation trait of Ca. Accumulibacter

Candidatus Accumulibacter, a prominent polyphosphate-accumulating organism (PAO) in wastewater treatment, plays a crucial role in enhanced biological phosphorus removal (EBPR). The genetic underpinnings of its polyphosphate accumulation capabilities, however, remain largely unknown. Here, we conducted a comprehensive genomic analysis of Ca. Accumulibacter-PAOs and their relatives within the Rhodocyclaceae family, identifying 124 core genes acquired via horizontal gene transfer (HGT) at its least common ancestor. Metatranscriptomic analysis of an enrichment culture of Ca. Accumulibacter revealed active transcription of 44 of these genes during an EBPR cycle, notably including the polyphosphate kinase 2 (PPK2) gene instead of the commonly recognized polyphosphate kinase 1 (PPK1) gene. Intriguingly, the phosphate regulon (Pho) genes showed minimal transcriptions, pointing to a distinctive fact of Pho dysregulation, where PhoU, the phosphate signaling complex protein, was not regulating the high-affinity phosphate transport (Pst) system, resulting in continuous phosphate uptake. To prevent phosphate toxicity, Ca. Accumulibacter utilized the laterally acquired PPK2 to condense phosphate into polyphosphate, resulting in the polyphosphate-accumulating feature. This study provides novel insights into the evolutionary emergence of the polyphosphate-accumulating trait in Ca. Accumulibacter, offering potential advancements in understanding the PAO phenotype in the EBPR process.

Complete Mitochondrial Genome and Phylogenetic Analysis of the Blue Whistling Thrush (Myophonus caeruleus)

The blue whistling thrush (Myophonus caeruleus) is a bird belonging to the order Passeriformes and family Muscicapidae. M. caeruleus is widely distributed in China, Pakistan, India, and Myanmar and is a resident bird in the southern part of the Yangtze River in China and summer migratory bird in the northern part of the Yangtze River. At present, there are some controversies about the classification of M. caeruleus. We use complete mitochondrial genomes to provide insights into the phylogenetic position of M. caeruleus and its relationships among Muscicapidae. The mitochondrial genome (GenBank: MN564936) is 16,815 bp long and contains 13 protein-coding genes (PCGs), 2 rRNA genes, 22 tRNA genes, and a non-coding control region (D-loop). The thirteen PCGs started with GTG and ATG and ended with five types of stop codons. The nucleotide composition of T was 23.71%, that of C was 31.45%, that of A was 30.06%, and that of G was 14.78%. The secondary structures of 22 tRNAs were predicted, all of which could form typical cloverleaf structures. There were 24 mismatches, mainly G-U mismatches. Through phylogenetic tree reconstruction, it was found that Saxicola, Monticola, Oenanthe, and Phoenicurus were clustered into one clade, together with the sister group of Myophonus.

Paleoneurology of stem palaeognaths clarifies the plesiomorphic condition of the crown bird central nervous system

Lithornithidae, an assemblage of volant Palaeogene fossil birds, provide our clearest insights into the early evolutionary history of Palaeognathae, the clade that today includes the flightless ratites and volant tinamous. The neotype specimen of Lithornis vulturinus, from the early Eocene (approximately 53 million years ago) of Europe, includes a partial neurocranium that has never been thoroughly investigated. Here, we describe these cranial remains including the nearly complete digital endocasts of the brain and bony labyrinth. The telencephalon of Lithornis is expanded and its optic lobes are ventrally shifted, as is typical for crown birds. The foramen magnum is positioned caudally, rather than flexed ventrally as in some crown birds, with the optic lobes, cerebellum, and foramen magnum shifted further ventrally. The overall brain shape is similar to that of tinamous, the only extant clade of flying palaeognaths, suggesting that several aspects of tinamou neuroanatomy may have been evolutionarily conserved since at least the early Cenozoic. The estimated ratio of the optic lobe's surface area relative to the total brain suggests a diurnal ecology. Lithornis may provide the clearest insights to date into the neuroanatomy of the ancestral crown bird, combining an ancestrally unflexed brain with a caudally oriented connection with the spinal cord, a moderately enlarged telencephalon, and ventrally shifted, enlarged optic lobes.

ETV2 induces endothelial, but not hematopoietic, lineage specification in birds

Cardiovascular system develops from the lateral plate mesoderm. Its three primary cell lineages (hematopoietic, endothelial, and muscular) are specified by the sequential actions of conserved transcriptional factors. ETV2, a master regulator of mammalian hemangioblast development, however, is absent in the chicken genome and acts downstream of NPAS4L in zebrafish. Here, we investigated the epistatic relationship between NPAS4L and ETV2 in avian hemangioblast development. We showed that ETV2 is deleted in all 363 avian genomes analyzed. Mouse ETV2 induced LMO2, but not NPAS4L or SCL, expression in chicken mesoderm. Squamate (lizards, geckos, and snakes) genomes contain both NPAS4L and ETV2 In Madagascar ground gecko, both genes were expressed in developing hemangioblasts. Gecko ETV2 induced only LMO2 in chicken mesoderm. We propose that both NPAS4L and ETV2 were present in ancestral amniote, with ETV2 acting downstream of NPAS4L in endothelial lineage specification. ETV2 may have acted as a pioneer factor by promoting chromatin accessibility of endothelial-specific genes and, in parallel with NPAS4L loss in ancestral mammals, has gained similar function in regulating blood-specific genes.

Accurate, scalable, and fully automated inference of species trees from raw genome assemblies using ROADIES

Inference of species trees plays a crucial role in advancing our understanding of evolutionary relationships and has immense significance for diverse biological and medical applications. Extensive genome sequencing efforts are currently in progress across a broad spectrum of life forms, holding the potential to unravel the intricate branching patterns within the tree of life. However, estimating species trees starting from raw genome sequences is quite challenging, and the current cutting-edge methodologies require a series of error-prone steps that are neither entirely automated nor standardized. In this paper, we present ROADIES, a novel pipeline for species tree inference from raw genome assemblies that is fully automated, easy to use, scalable, free from reference bias, and provides flexibility to adjust the tradeoff between accuracy and runtime. The ROADIES pipeline eliminates the need to align whole genomes, choose a single reference species, or pre-select loci such as functional genes found using cumbersome annotation steps. Moreover, it leverages recent advances in phylogenetic inference to allow multi-copy genes, eliminating the need to detect orthology. Using the genomic datasets released from large-scale sequencing consortia across three diverse life forms (placental mammals, pomace flies, and birds), we show that ROADIES infers species trees that are comparable in quality with the state-of-the-art approaches but in a fraction of the time. By incorporating optimal approaches and automating all steps from assembled genomes to species and gene trees, ROADIES is poised to improve the accuracy, scalability, and reproducibility of phylogenomic analyses.

Role-reversed polyandry is associated with faster fast-Z in shorebirds

In birds, males are homogametic and carry two copies of the Z chromosome ('ZZ'), while females are heterogametic and exhibit a 'ZW' genotype. The Z chromosome evolves at a faster rate than similarly sized autosomes, a phenomenon termed 'fast-Z evolution'. This is thought to be caused by two independent processes-greater Z chromosome genetic drift owing to a reduced effective population size, and stronger Z chromosome positive selection owing to the exposure of partially recessive alleles to selection. Here, we investigate the relative contributions of these processes by considering the effect of role-reversed polyandry on fast-Z in shorebirds, a paraphyletic group of wading birds that exhibit unusually diverse mating systems. We find stronger fast-Z effects under role-reversed polyandry, which is consistent with particularly strong selection on polyandrous females driving the fixation of recessive beneficial alleles. This result contrasts with previous research in birds, which has tended to implicate a primary role of genetic drift in driving fast-Z variation. We suggest that this discrepancy can be interpreted in two ways-stronger sexual selection acting on polyandrous females overwhelms an otherwise central role of genetic drift, and/or sexual antagonism is also contributing significantly to fast-Z and is exacerbated in sexually dimorphic species.

A genomic hotspot of diversifying selection and structural change in the hoary bat (Lasiurus cinereus)

Background

Previous work found that numerous genes positively selected within the hoary bat (Lasiurus cinereus) lineage are physically clustered in regions of conserved synteny. Here I further validate and expand on those finding utilizing an updated L. cinereus genome assembly and additional bat species as well as other tetrapod outgroups.

Methods

A chromosome-level assembly was generated by chromatin-contact mapping and made available by DNAZoo (www.dnazoo.org). The genomic organization of orthologous genes was extracted from annotation data for multiple additional bat species as well as other tetrapod clades for which chromosome-level assemblies were available from the National Center for Biotechnology Information (NCBI). Tests of branch-specific positive selection were performed for L. cinereus using PAML as well as with the HyPhy package for comparison.

Results

Twelve genes exhibiting significant diversifying selection in the L. cinereus lineage were clustered within a 12-Mb genomic window; one of these (Trpc4) also exhibited diversifying selection in bats generally. Ten of the 12 genes are landmarks of two distinct blocks of ancient synteny that are not linked in other tetrapod clades. Bats are further distinguished by frequent structural rearrangements within these synteny blocks, which are rarely observed in other Tetrapoda. Patterns of gene order and orientation among bat taxa are incompatible with phylogeny as presently understood, implying parallel evolution or subsequent reversals. Inferences of positive selection were found to be robust to alternative phylogenetic topologies as well as a strong shift in background nucleotide composition in some taxa.

Discussion

This study confirms and further localizes a genomic hotspot of protein-coding divergence in the hoary bat, one that also exhibits an increased tempo of structural change in bats compared with other mammals. Most genes in the two synteny blocks have elevated expression in brain tissue in humans and model organisms, and genetic studies implicate the selected genes in cranial and neurological development, among other functions.

Reference Genome Choice and Filtering Thresholds Jointly Influence Phylogenomic Analyses

Molecular phylogenies are a cornerstone of modern comparative biology and are commonly employed to investigate a range of biological phenomena, such as diversification rates, patterns in trait evolution, biogeography, and community assembly. Recent work has demonstrated that significant biases may be introduced into downstream phylogenetic analyses from processing genomic data; however, it remains unclear whether there are interactions among bioinformatic parameters or biases introduced through the choice of reference genome for sequence alignment and variant calling. We address these knowledge gaps by employing a combination of simulated and empirical data sets to investigate the extent to which the choice of reference genome in upstream bioinformatic processing of genomic data influences phylogenetic inference, as well as the way that reference genome choice interacts with bioinformatic filtering choices and phylogenetic inference method. We demonstrate that more stringent minor allele filters bias inferred trees away from the true species tree topology, and that these biased trees tend to be more imbalanced and have a higher center of gravity than the true trees. We find the greatest topological accuracy when filtering sites for minor allele count (MAC) >3-4 in our 51-taxa data sets, while tree center of gravity was closest to the true value when filtering for sites with MAC >1-2. In contrast, filtering for missing data increased accuracy in the inferred topologies; however, this effect was small in comparison to the effect of minor allele filters and may be undesirable due to a subsequent mutation spectrum distortion. The bias introduced by these filters differs based on the reference genome used in short read alignment, providing further support that choosing a reference genome for alignment is an important bioinformatic decision with implications for downstream analyses. These results demonstrate that attributes of the study system and dataset (and their interaction) add important nuance for how best to assemble and filter short-read genomic data for phylogenetic inference.

Extreme elevational migration spurred cryptic speciation in giant hummingbirds

The ecoevolutionary drivers of species niche expansion or contraction are critical for biodiversity but challenging to infer. Niche expansion may be promoted by local adaptation or constrained by physiological performance trade-offs. For birds, evolutionary shifts in migratory behavior permit the broadening of the climatic niche by expansion into varied, seasonal environments. Broader niches can be short-lived if diversifying selection and geography promote speciation and niche subdivision across climatic gradients. To illuminate niche breadth dynamics, we can ask how "outlier" species defy constraints. Of the 363 hummingbird species, the giant hummingbird (Patagona gigas) has the broadest climatic niche by a large margin. To test the roles of migratory behavior, performance trade-offs, and genetic structure in maintaining its exceptional niche breadth, we studied its movements, respiratory traits, and population genomics. Satellite and light-level geolocator tracks revealed an >8,300-km loop migration over the Central Andean Plateau. This migration included a 3-wk, ~4,100-m ascent punctuated by upward bursts and pauses, resembling the acclimatization routines of human mountain climbers, and accompanied by surging blood-hemoglobin concentrations. Extreme migration was accompanied by deep genomic divergence from high-elevation resident populations, with decisive postzygotic barriers to gene flow. The two forms occur side-by-side but differ almost imperceptibly in size, plumage, and respiratory traits. The high-elevation resident taxon is the world's largest hummingbird, a previously undiscovered species that we describe and name here. The giant hummingbirds demonstrate evolutionary limits on niche breadth: when the ancestral niche expanded due to evolution (or loss) of an extreme migratory behavior, speciation followed.

A highly contiguous genome assembly for the Yellow Warbler (Setophaga petechia)

The Yellow Warbler (Setophaga petechia) is a small songbird in the wood-warbler family (Parulidae) that exhibits phenotypic and ecological differences across a widespread distribution and is important to California's riparian habitat conservation. Here, we present a high-quality de novo genome assembly of a vouchered female Yellow Warbler from southern California. Using HiFi long-read and Omni-C proximity sequencing technologies, we generated a 1.22 Gb assembly including 687 scaffolds with a contig N50 of 6.80 Mb, scaffold N50 of 21.18 Mb, and a BUSCO completeness score of 96.0%. This highly contiguous genome assembly provides an essential resource for understanding the history of gene flow, divergence, and local adaptation in Yellow Warblers and can inform conservation management of this charismatic bird 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.

Wildlife ecological risk assessment in the 21st century: Promising technologies to assess toxicological effects

Despite advances in toxicity testing and the development of new approach methodologies (NAMs) for hazard assessment, the ecological risk assessment (ERA) framework for terrestrial wildlife (i.e., air-breathing amphibians, reptiles, birds, and mammals) has remained unchanged for decades. While survival, growth, and reproductive endpoints derived from whole-animal toxicity tests are central to hazard assessment, nonstandard measures of biological effects at multiple levels of biological organization (e.g., molecular, cellular, tissue, organ, organism, population, community, ecosystem) have the potential to enhance the relevance of prospective and retrospective wildlife ERAs. Other factors (e.g., indirect effects of contaminants on food supplies and infectious disease processes) are influenced by toxicants at individual, population, and community levels, and need to be factored into chemically based risk assessments to enhance the "eco" component of ERAs. Regulatory and logistical challenges often relegate such nonstandard endpoints and indirect effects to postregistration evaluations of pesticides and industrial chemicals and contaminated site evaluations. While NAMs are being developed, to date, their applications in ERAs focused on wildlife have been limited. No single magic tool or model will address all uncertainties in hazard assessment. Modernizing wildlife ERAs will likely entail combinations of laboratory- and field-derived data at multiple levels of biological organization, knowledge collection solutions (e.g., systematic review, adverse outcome pathway frameworks), and inferential methods that facilitate integrations and risk estimations focused on species, populations, interspecific extrapolations, and ecosystem services modeling, with less dependence on whole-animal data and simple hazard ratios. Integr Environ Assess Manag 2024;20:725-748. © 2023 His Majesty the King in Right of Canada and The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental Toxicology & Chemistry (SETAC). Reproduced with the permission of the Minister of Environment and Climate Change Canada. This article has been contributed to by US Government employees and their work is in the public domain in the USA.

Complexity of avian evolution revealed by family-level genomes

Despite tremendous efforts in the past decades, relationships among main avian lineages remain heavily debated without a clear resolution. Discrepancies have been attributed to diversity of species sampled, phylogenetic method and the choice of genomic regions1-3. Here we address these issues by analysing the genomes of 363 bird species4 (218 taxonomic families, 92% of total). Using intergenic regions and coalescent methods, we present a well-supported tree but also a marked degree of discordance. The tree confirms that Neoaves experienced rapid radiation at or near the Cretaceous-Palaeogene boundary. Sufficient loci rather than extensive taxon sampling were more effective in resolving difficult nodes. Remaining recalcitrant nodes involve species that are a challenge to model due to either extreme DNA composition, variable substitution rates, incomplete lineage sorting or complex evolutionary events such as ancient hybridization. Assessment of the effects of different genomic partitions showed high heterogeneity across the genome. We discovered sharp increases in effective population size, substitution rates and relative brain size following the Cretaceous-Palaeogene extinction event, supporting the hypothesis that emerging ecological opportunities catalysed the diversification of modern birds. The resulting phylogenetic estimate offers fresh insights into the rapid radiation of modern birds and provides a taxon-rich backbone tree for future comparative studies.

Genome composition and GC content influence loci distribution in reduced representation genomic studies

BACKGROUND

Genomic architecture is a key evolutionary trait for living organisms. Due to multiple complex adaptive and neutral forces which impose evolutionary pressures on genomes, there is a huge variability of genomic features. However, their variability and the extent to which genomic content determines the distribution of recovered loci in reduced representation sequencing studies is largely unexplored.

RESULTS

Here, by using 80 genome assemblies, we observed that whereas plants primarily increase their genome size by expanding their intergenic regions, animals expand both intergenic and intronic regions, although the expansion patterns differ between deuterostomes and protostomes. Loci mapping in introns, exons, and intergenic categories obtained by in silico digestion using 2b-enzymes are positively correlated with the percentage of these regions in the corresponding genomes, suggesting that loci distribution mostly mirrors genomic architecture of the selected taxon. However, exonic regions showed a significant enrichment of loci in all groups regardless of the used enzyme. Moreover, when using selective adaptors to obtain a secondarily reduced loci dataset, the percentage and distribution of retained loci also varied. Adaptors with G/C terminals recovered a lower percentage of selected loci, with a further enrichment of exonic regions, while adaptors with A/T terminals retained a higher percentage of loci and slightly selected more intronic regions than expected.

CONCLUSIONS

Our results highlight how genome composition, genome GC content, RAD enzyme choice and use of base-selective adaptors influence reduced genome representation techniques. This is important to acknowledge in population and conservation genomic studies, as it determines the abundance and distribution of loci.

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.

The complete mitochondrial genome of Hypsipetes amaurotis (Passeriformes: Pycnonotidae)

The brown-eared bulbul (Hypsipetes amaurotis) is a medium-sized songbird native to East Asia and characterized by its prominent reddish-brown ear-coverts. Previous studies on it have primarily been from the taxonomic and morphological aspects, with limited research in the realm of molecular biology. In this study, we sequenced and annotated the complete mitochondrial genome of H. amaurotis, which was the first reported complete mitogenome of the genus Hypsipetes. The mitogenome of H. amaurotis is 17,871 bp in length and was predicted to encode 37 typical mitochondrial genes, including 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs). Specifically, this mitogenome contains two D-loop control regions that are of similar length and sequencing pattern. A total of 8 Pycnonotidae and six outgroup taxa were used to determine the phylogenetic placement with two methods: Maximum Likelihood Approximation (IQ-TREE) and Bayesian inference (MrBayes). Our findings reveal that H. amaurotis is phylogenetically closely related to Ixos mcclellandii. The outcomes are generally consistent with the phylogenetic trees constructed in previous studies. The data gathered from this research provides valuable insights for future genomic investigations into the evolution, ecology, and conservation of this species.