Recent secondary contact, genome-wide admixture, and asymmetric introgression of neo-sex chromosomes between two Pacific island bird species.

Secondary contact between closely related taxa represents a "moment of truth" for speciation-an opportunity to test the efficacy of reproductive isolation that evolved in allopatry and to identify the genetic, behavioral, and/or ecological barriers that separate species in sympatry. Sex chromosomes are known to rapidly accumulate differences between species, an effect that may be exacerbated for neo-sex chromosomes that are transitioning from autosomal to sex-specific inheritance. Here we report that, in the Solomon Islands, two closely related bird species in the honeyeater family-Myzomela cardinalis and Myzomela tristrami-carry neo-sex chromosomes and have come into recent secondary contact after ~1.1 my of geographic isolation. Hybrids of the two species were first observed in sympatry ~100 years ago. To determine the genetic consequences of hybridization, we use population genomic analyses of individuals sampled in allopatry and in sympatry to characterize gene flow in the contact zone. Using genome-wide estimates of diversity, differentiation, and divergence, we find that the degree and direction of introgression varies dramatically across the genome. For sympatric birds, autosomal introgression is bidirectional, with phenotypic hybrids and phenotypic parentals of both species showing admixed ancestry. In other regions of the genome, however, the story is different. While introgression on the Z/neo-Z-linked sequence is limited, introgression of W/neo-W regions and mitochondrial sequence (mtDNA) is highly asymmetric, moving only from the invading M. cardinalis to the resident M. tristrami. The recent hybridization between these species has thus enabled gene flow in some genomic regions but the interaction of admixture, asymmetric mate choice, and/or natural selection has led to the variation in the amount and direction of gene flow at sex-linked regions of the genome.

Genetic diversity of a recovering European roller (Coracias garrulus) population from Serbia.

The European Roller (Coracias garrulus), a long-distance migratory bird, faced a considerable decline in breeding pairs throughout Europe at the end of the 20th century. Due to conservation efforts and the installation of nesting boxes, the population of the European Roller in Serbia has made a remarkable recovery. Here, we used the variability of nucleotide sequences of the mitochondrial DNA (mtDNA) control region and 10 microsatellite loci to assess the genetic diversity and structuring, phylogeographic patterns and demographic history of this species using 224 individuals from Serbia. Our results showed moderate level of genetic diversity (HO = 0.392) and a slightly elevated level of inbreeding and homozygosity (FIS = 0.393). Genetic structuring based on microsatellite data indicated three genetic clusters, but without a clear spatial pattern. High haplotype diversity (Hd = 0.987) of the mtDNA control region sequences was detected, and neutrality tests indicated a recent demographic expansion. The phylogeographic analysis, which also included previously published sequences of the mtDNA control region, supported the subdivision into two distinct European and Asian haplogroups (ΦST = 0.712). However, the results of our study showed that a larger number of haplotypes sampled in Serbia are clustered in the Asian haplogroup as compared to previous studies, indicating a historically continuous distribution of this species and possibly a wider distribution of the subspecies Coracias garrulus semenovwi. Our results suggest that the European Roller population in Serbia is genetically stable, with no evidence of recent bottlenecks, and emphasize the importance of artificial nest boxes for promoting and maintaining population dynamics of European Rollers.

CARE-ing for Indigenous nonhuman genomic data - rethinking our approach.

Earlier this year, an international group of scientists published a paper in ScienceAdvances on the draft genome of the little bush moa (Anomalopteryx didiformis), one of about nine species of extinct flightless birds endemic to Aotearoa (New Zealand). The genome was sequenced from the ancient DNA of a "poorly provenanced" fossil bone acquired by the Royal Ontario Museum. It held important clues about the moa's evolutionary history and aspects of its biology.

A comparison of antibiotic resistance genes and mobile genetic elements in wild and captive Himalayan vultures.

As the most widely distributed scavenger birds on the Qinghai-Tibetan Plateau, Himalayan vultures (Gyps himalayensis) feed on the carcasses of various wild and domestic animals, facing the dual selection pressure of pathogens and antibiotics and are suitable biological sentinel species for monitoring antibiotic resistance genes (ARGs). This study used metagenomic sequencing to comparatively investigate the ARGs and mobile genetic elements (MGEs) of wild and captive Himalayan vultures. Overall, the resistome of Himalayan vultures contained 414 ARG subtypes resistant to 20 ARG types, with abundances ranging from 0.01 to 1,493.60 ppm. The most abundant resistance type was beta-lactam (175 subtypes), followed by multidrug resistance genes with 68 subtypes. Decreases in the abundance of macrolide-lincosamide-streptogramin (MLS) resistance genes were observed in the wild group compared with the zoo group. A total of 75 genera (five phyla) of bacteria were predicted to be the hosts of ARGs in Himalayan vultures, and the clinical (102 ARGs) and high-risk ARGs (35 Rank I and 56 Rank II ARGs) were also analyzed. Among these ARGs, twenty-two clinical ARGs, nine Rank I ARG subtypes, sixteen Rank II ARG subtypes were found to differ significantly between the two groups. Five types of MGEs (128 subtypes) were found in Himalayan vultures. Plasmids (62 subtypes) and transposases (44 subtypes) were found to be the main MGE types. Efflux pump and antibiotic deactivation were the main resistance mechanisms of ARGs in Himalayan vultures. Decreases in the abundance of cellular protection were identified in wild Himalayan vultures compared with the captive Himalayan vultures. Procrustes analysis and the co-occurrence networks analysis revealed different patterns of correlations among gut microbes, ARGs, and MGEs in wild and captive Himalayan vultures. This study is the first step in describing the characterization of the ARGs in the gut of Himalayan vultures and highlights the need to pay more attention to scavenging birds.

Evolutionary variation in gene conversion at the avian MHC is explained by fluctuating selection, gene copy numbers and life history.

The major histocompatibility complex (MHC) multigene family encodes key pathogen-recognition molecules of the vertebrate adaptive immune system. Hyper-polymorphism of MHC genes is de novo generated by point mutations, but new haplotypes may also arise by re-shuffling of existing variation through intra- and inter-locus gene conversion. Although the occurrence of gene conversion at the MHC has been known for decades, we still have limited understanding of its functional importance. Here, I took advantage of extensive genetic resources (~9000 sequences) to investigate broad scale macroevolutionary patterns in gene conversion processes at the MHC across nearly 200 avian species. Gene conversion was found to constitute a universal mechanism in birds, as 83% of species showed footprints of gene conversion at either MHC class and 25% of all allelic variants were attributed to gene conversion. Gene conversion processes were stronger at MHC-II than MHC-I, but inter-specific variation at both MHC classes was explained by similar evolutionary scenarios, reflecting fluctuating selection towards different optima and drift. Gene conversion showed uneven phylogenetic distribution across birds and was driven by gene copy number variation, supporting significant role of inter-locus gene conversion processes in the evolution of the avian MHC. Finally, MHC gene conversion was stronger in species with fast life histories (high fecundity) and in long-distance migrants, likely reflecting variation in population sizes and host-pathogen coevolutionary dynamics. The results provide a robust comparative framework for understanding macroevolutionary variation in gene conversion at the avian MHC and reinforce important contribution of this mechanism to functional MHC diversity.

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.

Digest: How environmental light conditions shape the evolution of visual systems in birds.

How do varying environmental light conditions influence the evolution of avian visual systems? Fröhlich et al. (2024) demonstrate that nocturnal birds evolved broader corneas and slightly longer axial lengths than their diurnal counterparts, increasing light capture efficiency. Nocturnal species also tended to maintain or reduce the size of brain regions responsible for vision, i.e., the optic tectum and the visual wulst. These results highlight adaptive trends in nocturnal species, where evolutionary improvement in low-light performance of eyes may be accompanied by compromised brain function.

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.

Diversity of major histocompatibility complex of II B gene and mate choice in a monogamous and long-lived seabird, the Little Auk (Alle alle).

The major histocompatibility complex (MHC) plays a key role in the adaptive immune system of vertebrates, and is known to influence mate choice in many species. In birds, the MHC has been extensively examined but mainly in galliforms and passerines while other taxa that represent specific ecological and evolutionary life-histories, like seabirds, are underexamined. Here, we characterized diversity of MHC Class II B exon 2 in a colonial pelagic seabird, the Little Auk (or Dovekie Alle alle). We further examined whether MHC variation could be maintained through balancing selection and disassortative mating. We found high polymorphism at the genotyped MHC fragment, characterizing 99 distinct alleles across 140 individuals from three populations. The alleles frequencies exhibited a similar skewed distribution in both sexes, with the four most commonly occurring alleles representing approximately 35% of allelic variation. The results of a Bayesian site-by-site selection analysis suggest evidence of balancing selection and no direct evidence for MHC-dependent disassortative mating preferences in the Little Auk. The latter result might be attributed to the high overall polymorphism of the examined fragment, which itself may be maintained by the large population size of the species.

Mitonuclear interactions impact aerobic metabolism in hybrids and may explain mitonuclear discordance in young, naturally hybridizing bird lineages.

Understanding genetic incompatibilities and genetic introgression between incipient species are major goals in evolutionary biology. Mitochondrial genes evolve rapidly and exist in dense gene networks with coevolved nuclear genes, suggesting that mitochondrial respiration may be particularly susceptible to disruption in hybrid organisms. Mitonuclear interactions have been demonstrated to contribute to hybrid dysfunction between deeply divergent taxa crossed in the laboratory, but there are few empirical examples of mitonuclear interactions between younger lineages that naturally hybridize. Here, we use controlled hybrid crosses and high-resolution respirometry to provide the first experimental evidence in a bird that inter-lineage mitonuclear interactions impact mitochondrial aerobic metabolism. Specifically, respiration capacity of the two mitodiscordant backcrosses (with mismatched mitonuclear combinations) differs from one another, although they do not differ significantly from the parental groups or mitoconcordant backcrosses as we would expect of mitonuclear disruptions. In the wild hybrid zone between these subspecies, the mitochondrial cline centre is shifted west of the nuclear cline centre, which is consistent with the direction of our experimental results. Our results therefore demonstrate asymmetric mitonuclear interactions that impact the capacity of cellular mitochondrial respiration and may help to explain the geographic discordance between mitochondrial and nuclear genomes observed in the wild.

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.

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.

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.

Mammalian and avian species quantification in homogenized foods: real time PCR and digital PCR as tools for label compliance controls.

Currently food fraud and authenticity of products composition are topics of great concern; ingredients quantification could allow to identify small amounts of contaminats or voluntary addition of improper components. Many molecular methods are available for species identification in foodstuffs but, for a better application, they should not be affected by the interference of other ingredients. The main purpose of this work was to verify the Real Time PCR and the Digital PCR (dPCR) quantification performances on baby food samples, specifically selected for their high miscibility to limit variability; chicken was selected as target to verify the performance of quantification of methods after having spiked the same quantity in different baby foods. The other aims were: (1) to verify a constant genome copies ratio existence between mammalian and avian species (2) to verify the dPCR performance, set up on housekeeping, to quantify mammalian and avian species in commercial products. Digital PCR showed fewer differences respect to Real Time PCR, at the same 15% w/w chicken spiking level. Despite the constant difference between mammalian and avian genome copies, in samples with the same spiking weight, the confidence intervals increasing towards the extreme values, made impossible to use genome copies ratio as a sort of correction factor between species. Finally, the dPCR system using the myostatin housekeeping gene to determine the chicken content seemed reliable to verify the labelling compliance in meat-based commercial products.

Detection and genomic characterization of Klebsiella pneumoniae and Escherichia coli harboring tet(X4) in black kites (Milvus migrans) in Pakistan.

The emergence of plasmid-mediated tigecycline resistance gene tet(X4) among clinically relevant bacteria has promoted significant concerns, as tigecycline is considered a last-resort drug against serious infections caused by multidrug-resistant bacteria. We herein focused on the isolation and molecular characterization of tet(X4)-positive Klebsiella pneumoniae (K. pneumoniae) and Escherichia coli (E. coli) in wild bird populations with anthropogenic interaction in Faisalabad, Pakistan. A total of 150 birds including black kites (Milvus migrans) and house crows (Corvus splendens) were screened for the presence of tigecycline resistance K. pneumoniae and E. coli. We found two K. pneumoniae and one E. coli isolate carrying tet(X4) originating from black kites. A combination of short- and long-read sequencing strategies showed that tet(X4) was located on a broad host range IncFII plasmid family in K. pneumoniae isolates whereas on an IncFII-IncFIB hybrid plasmid in E. coli. We also found an integrative and conjugative element ICEKp2 in K. pneumoniae isolate KP8336. We demonstrate the first description of tet(X4) gene in the WHO critical-priority pathogen K. pneumoniae among wild birds. The convergence of tet(X4) and virulence associated ICEKp2 in a wild bird with known anthropogenic contact should be further investigated to evaluate the potential epidemiological implications. The potential risk of global transmission of tet(X4)-positive K. pneumoniae and E. coli warrant comprehensive evaluation and emphasizes the need for effective mitigation strategies to reduce anthropogenic-driven dissemination of AMR in the environment.

Evolution and genetic architecture of sex-limited polymorphism in cuckoos.

Sex-limited polymorphism has evolved in many species including our own. Yet, we lack a detailed understanding of the underlying genetic variation and evolutionary processes at work. The brood parasitic common cuckoo (Cuculus canorus) is a prime example of female-limited color polymorphism, where adult males are monochromatic gray and females exhibit either gray or rufous plumage. This polymorphism has been hypothesized to be governed by negative frequency-dependent selection whereby the rarer female morph is protected against harassment by males or from mobbing by parasitized host species. Here, we show that female plumage dichromatism maps to the female-restricted genome. We further demonstrate that, consistent with balancing selection, ancestry of the rufous phenotype is shared with the likewise female dichromatic sister species, the oriental cuckoo (Cuculus optatus). This study shows that sex-specific polymorphism in trait variation can be resolved by genetic variation residing on a sex-limited chromosome and be maintained across species boundaries.

Family living and cooperative breeding in birds are associated with the number of avian predators.

Cooperative breeding occurs when individuals contribute parental care to offspring that are not their own. Numerous intra- and interspecific studies have aimed to explain the evolution of this behavior. Recent comparative work suggests that family living (i.e., when offspring remain with their parents beyond independence) is a critical stepping stone in the evolution of cooperative breeding. Thus, it is key to understand the factors that facilitate the evolution of family living. Within-species studies suggest that protection from predators is a critical function of group living, through both passive benefits such as dilution effects and active benefits such as prosocial antipredator behaviors in family groups. However, the association between predation risk and the formation and prevalence of family groups and cooperative breeding remains untested globally. Here, we use phylogenetic comparative analyses including 2,984 bird species to show that family living and cooperative breeding are associated with increased occurrence of avian predators. These cross-species findings lend support to previous suggestions based on intraspecific studies that social benefits of family living, such as protection against predation, could favor the evolution of delayed dispersal and cooperative breeding.

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.

Global assessment of effective population sizes: Consistent taxonomic differences in meeting the 50/500 rule.

Effective population size (Ne) is a particularly useful metric for conservation as it affects genetic drift, inbreeding and adaptive potential within populations. Current guidelines recommend a minimum Ne of 50 and 500 to avoid short-term inbreeding and to preserve long-term adaptive potential respectively. However, the extent to which wild populations reach these thresholds globally has not been investigated, nor has the relationship between Ne and human activities. Through a quantitative review, we generated a dataset with 4610 georeferenced Ne estimates from 3829 populations, extracted from 723 articles. These data show that certain taxonomic groups are less likely to meet 50/500 thresholds and are disproportionately impacted by human activities; plant, mammal and amphibian populations had a <54% probability of reaching N ̂ e = 50 and a <9% probability of reaching N ̂ e = 500. Populations listed as being of conservation concern according to the IUCN Red List had a smaller median N ̂ e than unlisted populations, and this was consistent across all taxonomic groups. N ̂ e was reduced in areas with a greater Global Human Footprint, especially for amphibians, birds and mammals, however relationships varied between taxa. We also highlight several considerations for future works, including the role that gene flow and subpopulation structure plays in the estimation of N ̂ e in wild populations, and the need for finer-scale taxonomic analyses. Our findings provide guidance for more specific thresholds based on Ne and help prioritise assessment of populations from taxa most at risk of failing to meet conservation thresholds.