The first whole genome and transcriptome of the cinereous vulture reveals adaptation in the gastric and immune defense systems and possible convergent evolution between the Old and New World vultures

Necrophages and necrophiles: a review of their antibacterial defenses and biotechnological potential

With antibiotic resistance on the rise, there is an urgent need for new antibacterial drugs and products to treat or prevent infection. Many such products in current use, for example human and veterinary antibiotics and antimicrobial food preservatives, were discovered and developed from nature. Natural selection acts on all living organisms and the presence of bacterial competitors or pathogens in an environment can favor the evolution of antibacterial adaptations. In this review, we ask if vultures, blow flies and other carrion users might be a good starting point for antibacterial discovery based on the selection pressure they are under from bacterial disease. Dietary details are catalogued for over 600 of these species, bacterial pathogens associated with the diets are described, and an overview of the antibacterial defenses contributing to disease protection is given. Biotechnological applications for these defenses are then discussed, together with challenges facing developers and possible solutions. Examples include use of (a) the antimicrobial peptide (AMP) gene sarcotoxin IA to improve crop resistance to bacterial disease, (b) peptide antibiotics such as serrawettin W2 as antibacterial drug leads, (c) lectins for targeted drug delivery, (d) bioconversion-generated chitin as an antibacterial biomaterial, (e) bacteriocins as antibacterial food preservatives and (f) mutualistic microbiota bacteria as alternatives to antibiotics in animal feed. We show that carrion users encounter a diverse range of bacterial pathogens through their diets and interactions, have evolved many antibacterial defenses, and are a promising source of genes, molecules, and microbes for medical, agricultural, and food industry product development.

The Role of Vulture (Accipitriformes) Cutaneous Microbiota in Infectious Disease Protection

Vultures (Accipitriformes), as obligate scavengers, are regularly exposed to a diverse array of pathogens present in decomposing carcasses. Nevertheless, they exhibit a remarkable ability to resist infections, suggesting a crucial role of skin microbiota in host defense. The microbial communities residing on necrophagic birds' skin create a protective barrier through competitive interactions, antimicrobial compound production, and immunity priming. Additionally, vultures contribute to ecosystem balance by reducing the spread of infectious agents. However, they may also serve as vectors for antimicrobial resistance (AMR) due to their exposure to contaminated food sources. Understanding the dynamics of their microbiota can provide valuable insights into host-microbe interactions, wildlife conservation, and public health. This review examines the composition and functional significance of vulture cutaneous microbiota. Specifically, it explores the role of necrophagic birds' skin microbiota in pathogen exclusion, immune system modulation, and environmental adaptation, with the aim of suggesting further research routes, besides clarifying the ecological implications of such birds.

Metagenomic comparison of gut communities between wild and captive Himalayan griffons

Introduction

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

Methods

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

Results

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

Discussion

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

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

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

First metagenomic analysis of the Andean condor (Vultur gryphus) gut microbiome reveals microbial diversity and wide resistome

Background

The Andean condor (Vultur gryphus) is the largest scavenger in South America. This predatory bird plays a crucial role in their ecological niche by removing carcasses. We report the first metagenomic analysis of the Andean condor gut microbiome.

Methods

This work analyzed shotgun metagenomics data from a mixture of fifteen captive Chilean Andean condors. To filter eukaryote contamination, we employed BWA-MEM v0.7. Taxonomy assignment was performed using Kraken2 and MetaPhlAn v2.0 and all filtered reads were assembled using IDBA-UD v1.1.3. The two most abundant species were used to perform a genome reference-guided assembly using MetaCompass. Finally, we performed a gene prediction using Prodigal and each gene predicted was functionally annotated. InterproScan v5.31-70.0 was additionally used to detect homology based on protein domains and KEGG mapper software for reconstructing metabolic pathways.

Results

Our results demonstrate concordance with the other gut microbiome data from New World vultures. In the Andean condor, Firmicutes was the most abundant phylum present, with Clostridium perfringens, a potentially pathogenic bacterium for other animals, as dominating species in the gut microbiome. We assembled all reads corresponding to the top two species found in the condor gut microbiome, finding between 94% to 98% of completeness for Clostridium perfringens and Plesiomonas shigelloides, respectively. Our work highlights the ability of the Andean condor to act as an environmental reservoir and potential vector for critical priority pathogens which contain relevant genetic elements. Among these genetic elements, we found 71 antimicrobial resistance genes and 1,786 virulence factors that we associated with several adaptation processes.

Evaluation of Oral Transmucosal Administration of Pentobarbital for Euthanasia of Conscious Wild Birds

This prospective study evaluated oral transmucosal pentobarbital sodium at three doses in 110 wild-caught wild birds requiring euthanasia. Birds received transmucosal pentobarbital at five (430 mg/kg), six (516 mg/kg), and seven times (602 mg/kg) the intravenous dose for mammals. Time to first effects and loss of consciousness, presence of pupillary light and corneal reflexes, apnea, and asystole were recorded each minute. When asystole was not achieved at 5 minutes, IV pentobarbital was administered. Combining data for all doses, loss of consciousness occurred at a median (range) of 2 minutes (0-4.75 min), apnea at 3 minutes (0-6 min), and asystole at 4 minutes (0.5-5 min). Loss of consciousness and apnea occurred significantly faster in the 602 mg/kg dose group than in the 430 mg/kg group (p = 0.009, difference of 0.6 ± 0.2 min; p = 0.024, difference of 0.7 ± 0.3 min), respectively. Apnea and asystole were achieved in 80/110 birds within 5 minutes. Oral transmucosal pentobarbital results in rapid loss of consciousness and respiratory arrest and provides a reliable alternative euthanasia method compared to intravenous administration.

Fungal signatures of oral disease reflect environmental degradation in a facultative avian scavenger

Degradation of natural ecosystems increases the risk of infections in wildlife due to microbiota dysbiosis. However, little is known about its influence on the development of fungal communities in predators and facultative avian scavengers. We evaluated the incidence of oral disease in wild nestling black kites (Milvus migrans) under contrasting environmental degradation conditions, and explored their oral fungal patterns using molecular methods and multivariate analysis. Oral lesions were found in 36.8% of the 38 nestlings examined in an anthropogenically altered habitat (southeastern Madrid, Spain), but in none of the 105 nestlings examined in a well-conserved natural area (Doñana National Park, Spain). In a subsample of 48 black kites, the composition of the oral fungal community differed among symptomatic nestlings from Madrid (SM) and asymptomatic nestlings from Madrid (AM) and Doñana (AD). Opportunistic fungal pathogens (e.g., Fusarium incarnatum-equiseti species complex, Mucor spp., Rhizopus oryzae) were more prevalent in SM and AM than in AD. Hierarchical clustering and principal component analyses revealed that fungal patterns were distinct between both study areas, and that anthropogenic and natural environmental factors had a greater impact on them than oral disease. Fungal signatures associated with anthropogenic and natural stresses harbored some taxa that could be used to flag oral infection (F. incarnatum-equiseti species complex and Alternaria), indicate environmental degradation (Alternaria) or provide protective benefits in degraded environments (Trichoderma, Epicoccum nigrum and Sordaria). Co-occurrence associations between potentially beneficial and pathogenic fungi were typical of AM and AD, hinting at a possible role in host health. This study shows that early-life exposure to highly degraded environments induces a shift towards a higher prevalence of pathogenic species in the oral cavity of black kites, favoring oral disease. Furthermore, our findings suggest potential ecological applications of the monitoring of oral mycobiome as a bioindication of oral disease and environmental degradation.

Long-read genome sequencing provides molecular insights into scavenging and societal complexity in spotted hyena Crocuta crocuta

The spotted hyena (Crocuta crocuta) is a large and unique terrestrial carnivore. It is a particularly fascinating species due to its distinct phenotypic traits, especially its complex social structure and scavenging lifestyle, with associated high dietary exposure to microbial pathogens. However, the underlying molecular mechanisms related to these phenotypes remain elusive. Here, we sequenced and assembled a high-quality long-read genome of the spotted hyena, with a contig N50 length of ∼13.75 Mb. Based on comparative genomics, immunoglobulin family members (e.g., IGKV4-1) showed significant adaptive duplications in the spotted hyena and striped hyena. Furthermore, immune-related genes (e.g., CD8A, LAG3, and TLR3) experienced species-specific positive selection in the spotted hyena lineage. These results suggest that immune tolerance between the spotted hyena and closely related striped hyena has undergone adaptive divergence to cope with prolonged dietary exposure to microbial pathogens from scavenging. Furthermore, we provided the potential genetic insights underlying social complexity, hinting at social behavior and cognition. Specifically, the RECNE-associated genes (e.g., UGP2 and ACTR2) in the spotted hyena genome are involved in regulation of social communication. Taken together, our genomic analyses provide molecular insights into the scavenging lifestyle and societal complexity of spotted hyenas.

Why Did the Bee Eat the Chicken? Symbiont Gain, Loss, and Retention in the Vulture Bee Microbiome

Diet and gut microbiomes are intricately linked on both short and long timescales. Changes in diet can alter the microbiome, while microbes in turn allow hosts to access novel diets. Bees are wasps that switched to a vegetarian lifestyle, and the vast majority of bees feed on pollen and nectar. Some stingless bee species, however, also collect carrion, and a few have fully reverted to a necrophagous lifestyle, relying on carrion for protein and forgoing flower visitation altogether. These “vulture” bees belong to the corbiculate apid clade, which is known for its ancient association with a small group of core microbiome phylotypes. Here, we investigate the vulture bee microbiome, along with closely related facultatively necrophagous and obligately pollinivorous species, to understand how these diets interact with microbiome structure. Via deep sequencing of the 16S rRNA gene and subsequent community analyses, we find that vulture bees have lost some core microbes, retained others, and entered into novel associations with acidophilic microbes found in the environment and on carrion. The abundance of acidophilic bacteria suggests that an acidic gut is important for vulture bee nutrition and health, as has been found in other carrion-feeding animals. Facultatively necrophagous bees have more variable microbiomes than strictly pollinivorous bees, suggesting that bee diet may interact with microbiomes on both short and long timescales. Further study of vulture bees promises to provide rich insights into the role of the microbiome in extreme diet switches.

Competition between obligate and facultative scavengers and infection: vulture-jackal-anthrax dynamics in Etosha National Park

Different species of scavengers may compete for the same food in an ecosystem. This case study considers the competition between jackals and vultures in Etosha National Park in Namibia. While jackals are facultative scavengers, able to hunt for food if necessary, vultures are obligate scavengers wholly dependent on carcasses of animals like zebras for persistence. This competition may be further affected by outbreaks of infections such as anthrax, which temporarily increase the number of carcasses but lower the zebra population, acting in some ways as a third competitor. We use a dynamical system to model the interplay between competition dynamics and infection dynamics, and how it is affected by the nature of the competition: indirect (exploitative) or direct (interference). A bifurcation analysis using reproduction numbers shows how vultures' survival may depend on their direct competitive edge in reaching carcasses faster than jackals, and how the infection and the scavengers complicate each other's persistence. Vultures' interference causes a backward bifurcation which enables them to persist. One possible outcome is a "strange bedfellows" bistability in which anthrax and vultures persist only together, not apart, despite being competitors.

Immune and hormonal regulation of the Boa constrictor (Serpentes; Boidae) in response to feeding

Feeding upregulates immune function and the systemic and local (gastrointestinal tract) concentrations of some immunoregulatory hormones, as corticosterone (CORT) and melatonin (MEL), in mammals and anurans. However, little is known about the immune and hormonal regulation in response to feeding in other ectothermic vertebrates, especially snakes, in which the postprandial metabolic changes are pronounced. Here, we investigated the effects feeding have on hormonal and innate immune responses in the snake, Boa constrictor. We divided juvenile males into two groups: fasting and fed with mice (30% of body mass). We measured the rates of oxygen consumption, plasma CORT levels, heterophil/lymphocyte ratio (HL ratio), plasma bacterial killing ability (BKA), and stomach and intestine MEL in fasting snakes and 48 h after meal intake. We observed increased rates of oxygen consumption, plasma CORT levels, and HL ratio, along with a tendency of decreased stomach and intestine MEL in fed snakes compared to fasting ones. BKA was not affected by feeding. Overall, we found that feeding modulates metabolic rates, CORT levels, and immune cell distribution in boas. Increased baseline CORT may be important to mobilize energy to support the metabolic increment during the postprandial period. Increased HL ratio might be an immunoregulatory effect of increased CORT, which has been shown in different physiological situations such as in response to immune challenge. Our results suggest that feeding activates the hypothalamic-pituitary-adrenal axis and modulates immune cell redistribution, possibly contributing to fighting potential injuries and infections derived from predation and from pathogens present in ingested food.

Genomic bases underlying the adaptive radiation of core landbirds

BACKGROUND

Core landbirds undergo adaptive radiation with different ecological niches, but the genomic bases that underlie their ecological diversification remain unclear.

RESULTS

Here we used the genome-wide target enrichment sequencing of the genes related to vision, hearing, language, temperature sensation, beak shape, taste transduction, and carbohydrate, protein and fat digestion and absorption to examine the genomic bases underlying their ecological diversification. Our comparative molecular phyloecological analyses show that different core landbirds present adaptive enhancement in different aspects, and two general patterns emerge. First, all three raptorial birds (Accipitriformes, Strigiformes, and Falconiformes) show a convergent adaptive enhancement for fat digestion and absorption, while non-raptorial birds tend to exhibit a promoted capability for protein and carbohydrate digestion and absorption. Using this as a molecular marker, our results show relatively strong support for the raptorial lifestyle of the common ancestor of core landbirds, consequently suggesting a single origin of raptors, followed by two secondary losses of raptorial lifestyle within core landbirds. In addition to the dietary niche, we find at temporal niche that diurnal birds tend to exhibit an adaptive enhancement in bright-light vision, while nocturnal birds show an increased adaption in dim-light vision, in line with previous findings.

CONCLUSIONS

Our molecular phyloecological study reveals the genome-wide adaptive differentiations underlying the ecological diversification of core landbirds.

Vulture Genomes Reveal Molecular Adaptations Underlying Obligate Scavenging and Low Levels of Genetic Diversity

Obligate scavenging on the dead and decaying animal matter is a rare dietary specialization that in extant vertebrates is restricted to vultures. These birds perform essential ecological services, yet many vulture species have undergone recent steep population declines and are now endangered. To test for molecular adaptations underlying obligate scavenging in vultures, and to assess whether genomic features might have contributed to their population declines, we generated high-quality genomes of the Himalayan and bearded vultures, representing both independent origins of scavenging within the Accipitridae, alongside a sister taxon, the upland buzzard. By comparing our data to published sequences from other birds, we show that the evolution of obligate scavenging in vultures has been accompanied by widespread positive selection acting on genes underlying gastric acid production, and immunity. Moreover, we find evidence of parallel molecular evolution, with amino acid replacements shared among divergent lineages of these scavengers. Our genome-wide screens also reveal that both the Himalayan and bearded vultures exhibit low levels of genetic diversity, equating to around a half of the mean genetic diversity of other bird genomes examined. However, demographic reconstructions indicate that population declines began at around the Last Glacial Maximum, predating the well-documented dramatic declines of the past three decades. Taken together, our genomic analyses imply that vultures harbor unique adaptations for processing carrion, but that modern populations are genetically depauperate and thus especially vulnerable to further genetic erosion through anthropogenic activities.

Genome-wide diversity in the California condor tracks its prehistoric abundance and decline

Due to their small population sizes, threatened and endangered species frequently suffer from a lack of genetic diversity, potentially leading to inbreeding depression and reduced adaptability.¹ During the latter half of the twentieth century, North America's largest soaring bird,² the California condor (Gymnogyps californianus; Critically Endangered³), briefly went extinct in the wild. Though condors once ranged throughout North America, by 1982 only 22 individuals remained. Following decades of captive breeding and release efforts, there are now >300 free-flying wild condors and ∼200 in captivity. The condor's recent near-extinction from lead poisoning, poaching, and loss of habitat is well documented,⁴ but much about its history remains obscure. To fill this gap and aid future management of the species, we produced a high-quality chromosome-length genome assembly for the California condor and analyzed its genome-wide diversity. For comparison, we also examined the genomes of two close relatives: the Andean condor (Vultur gryphus; Vulnerable³) and the turkey vulture (Cathartes aura; Least Concern³). The genomes of all three species show evidence of historic population declines. Interestingly, the California condor genome retains a high degree of variation, which our analyses reveal is a legacy of its historically high abundance. Correlations between genome-wide diversity and recombination rate further suggest a history of purifying selection against linked deleterious alleles, boding well for future restoration. We show how both long-term evolutionary forces and recent inbreeding have shaped the genome of the California condor, and provide crucial genomic resources to enable future research and conservation.

Molecular dating for phylogenies containing a mix of populations and species by using Bayesian and RelTime approaches

Simultaneous molecular dating of population and species divergences is essential in many biological investigations, including phylogeography, phylodynamics and species delimitation studies. In these investigations, multiple sequence alignments consist of both intra- and interspecies samples (mixed samples). As a result, the phylogenetic trees contain interspecies, interpopulation and within-population divergences. Bayesian relaxed clock methods are often employed in these analyses, but they assume the same tree prior for both inter- and intraspecies branching processes and require specification of a clock model for branch rates (independent vs. autocorrelated rates models). We evaluated the impact of a single tree prior on Bayesian divergence time estimates by analysing computer-simulated data sets. We also examined the effect of the assumption of independence of evolutionary rate variation among branches when the branch rates are autocorrelated. Bayesian approach with coalescent tree priors generally produced excellent molecular dates and highest posterior densities with high coverage probabilities. We also evaluated the performance of a non-Bayesian method, RelTime, which does not require the specification of a tree prior or a clock model. RelTime's performance was similar to that of the Bayesian approach, suggesting that it is also suitable to analyse data sets containing both populations and species variation when its computational efficiency is needed.

Vultures from different trophic guilds show distinct oral pathogenic yeast signatures and co-occurrence networks

Vultures have evolved adaptive mechanisms to prevent infections associated with their scavenging lifestyle. However, food-borne exposure to antimicrobial pharmaceuticals can promote opportunistic infections with adverse outcomes. Here, we used multivariate and network analyses to increase understanding of the behavior of the yeast communities causing oral mycosis outbreaks recently reported in wild nestling cinereous (Aegypius monachus), griffon (Gyps fulvus) and Egyptian (Neophron percnopterus) vultures (CV, GV and EV, respectively) exposed to antibiotics from livestock farming. Common and unique yeast signatures (of Candida, Debaromyces, Diutina, Meyerozyma, Naganishia, Pichia, Rhodotorula, Trichosporon and Yarrowia species) associated with oral mycoses were identified in the three vulture species. Hierarchical clustering analysis (HCA) and principal component analysis (PCA) highlighted that oral lesions from CV and GV shared similar yeast signatures (of major causative pathogens of opportunistic mycoses, such as Candida albicans, Candida parapsilosis and Candida tropicalis), while EV had a distinct yeast signature (of uncommon pathogenic species, such as Candida dubliniensis, Candida zeylanoides, Pichia fermentans and Rhodotorula spp.). Synergistic interactions between yeast species from distinct fungal phyla were found in lesions from CV and GV, but not in EV. These formed co-occurrence subnetworks with partially or fully connected topology. This study reveals that the composition, assembly and co-occurrence patterns of the yeast communities causing oral mycoses differ between vulture species with distinct feeding habits and scavenging lifestyles. Yeast species widely pathogenic to humans and animals, and yeast co-occurrence relationships, are distinctive hallmarks of oral mycoses in CV and GV. These vulture species are more exposed to antibiotics from intensively medicated livestock carcasses provided in supplementary feeding stations and show higher incidence of thrush-like oral lesions than EV. These findings may be useful for development of new initiatives or changes in the conservation of these avian scavengers affected by anthropogenic activities.

New genome assembly of the barn owl (Tyto alba alba)

New genomic tools open doors to study ecology, evolution, and population genomics of wild animals. For the Barn owl species complex, a cosmopolitan nocturnal raptor, a very fragmented draft genome was assembled for the American species (Tyto furcata pratincola) (Jarvis et al. 2014). To improve the genome, we assembled de novo Illumina and Pacific Biosciences (PacBio) long reads sequences of its European counterpart (Tyto alba alba). This genome assembly of 1.219 Gbp comprises 21,509 scaffolds and results in a N50 of 4,615,526 bp. BUSCO (Universal Single-Copy Orthologs) analysis revealed an assembly completeness of 94.8% with only 1.8% of the genes missing out of 4,915 avian orthologs searched, a proportion similar to that found in the genomes of the zebra finch (Taeniopygia guttata) or the collared flycatcher (Ficedula albicollis). By mapping the reads of the female American barn owl to the male European barn owl reads, we detected several structural variants and identified 70 Mbp of the Z chromosome. The barn owl scaffolds were further mapped to the chromosomes of the zebra finch. In addition, the completeness of the European barn owl genome is demonstrated with 94 of 128 proteins missing in the chicken genome retrieved in the European barn owl transcripts. This improved genome will help future barn owl population genomic investigations.

What can be learned by scanning the genome for molecular convergence in wild populations?

Convergent evolution, where independent lineages evolve similar phenotypes in response to similar challenges, can provide valuable insight into how selection operates and the limitations it encounters. However, it has only recently become possible to explore how convergent evolution is reflected at the genomic level. The overlapping outlier approach (OOA), where genome scans of multiple independent lineages are used to find outliers that overlap and therefore identify convergently evolving loci, is becoming popular. Here, we present a quantitative analysis of 34 studies that used this approach across many sampling designs, taxa, and sampling intensities. We found that OOA studies with increased biological sampling power within replicates have increased likelihood of finding overlapping, "convergent" signals of adaptation between them. When identifying convergent loci as overlapping outliers, it is tempting to assume that any false-positive outliers derived from individual scans will fail to overlap across replicates, but this cannot be guaranteed. We highlight how population demographics and genomic context can contribute toward both true convergence and false positives in OOA studies. We finish with an exploration of emerging methods that couple genome scans with phenotype and environmental measures, leveraging added information from genome data to more directly test hypotheses of the likelihood of convergent evolution.

Blood transcriptome resources of chinstrap (Pygoscelis antarcticus) and gentoo (Pygoscelis papua) penguins from the South Shetland Islands, Antarctica

The chinstrap (Pygoscelis antarcticus) and gentoo (P. papua) penguins are distributed throughout Antarctica and the sub-Antarctic islands. In this study, high-quality de novo assemblies of blood transcriptomes from these penguins were generated using the Illumina MiSeq platform. A total of 22.2 and 21.8 raw reads were obtained from chinstrap and gentoo penguins, respectively. These reads were assembled using the Oases assembly platform and resulted in 26,036 and 21,854 contigs with N50 values of 929 and 933 base pairs, respectively. Functional gene annotations through pathway analyses of the Gene Ontology, EuKaryotic Orthologous Groups, and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases were performed for each blood transcriptome, resulting in a similar compositional order between the two transcriptomes. Ortholog comparisons with previously published transcriptomes from the Adélie (P. adeliae) and emperor (Aptenodytes forsteri) penguins revealed that a high proportion of the four penguins’ transcriptomes had significant sequence homology. Because blood and tissues of penguins have been used to monitor pollution in Antarctica, immune parameters in blood could be important indicators for understanding the health status of penguins and other Antarctic animals. In the blood transcriptomes, KEGG analyses detected many essential genes involved in the major innate immunity pathways, which are key metabolic pathways for maintaining homeostasis against exogenous infections or toxins. Blood transcriptome studies such as this may be useful for checking the immune and health status of penguins without sacrifice.

Protective role of the vulture facial skin and gut microbiomes aid adaptation to scavenging

Background

Vultures have adapted the remarkable ability to feed on carcasses that may contain microorganisms that would be pathogenic to most other animals. The holobiont concept suggests that the genetic basis of such adaptation may not only lie within their genomes, but additionally in their associated microbes. To explore this, we generated shotgun DNA sequencing datasets of the facial skin and large intestine microbiomes of the black vulture (Coragyps atratus) and the turkey vulture (Cathartes aura). We characterized the functional potential and taxonomic diversity of their microbiomes, the potential pathogenic challenges confronted by vultures, and the microbial taxa and genes that could play a protective role on the facial skin and in the gut.

Results

We found microbial taxa and genes involved in diseases, such as dermatitis and pneumonia (more abundant on the facial skin), and gas gangrene and food poisoning (more abundant in the gut). Interestingly, we found taxa and functions with potential for playing beneficial roles, such as antilisterial bacteria in the gut, and genes for the production of antiparasitics and insecticides on the facial skin. Based on the identified phages, we suggest that phages aid in the control and possibly elimination, as in phage therapy, of microbes reported as pathogenic to a variety of species. Interestingly, we identified Adineta vaga in the gut, an invertebrate that feeds on dead bacteria and protozoans, suggesting a defensive predatory mechanism. Finally, we suggest a colonization resistance role through biofilm formation played by Fusobacteria and Clostridia in the gut.

Conclusions

Our results highlight the importance of complementing genomic analyses with metagenomics in order to obtain a clearer understanding of the host-microbial alliance and show the importance of microbiome-mediated health protection for adaptation to extreme diets, such as scavenging.

Electronic supplementary material

The online version of this article (10.1186/s13028-018-0415-3) contains supplementary material, which is available to authorized users.

Effects of vulture exclusion on carrion consumption by facultative scavengers

Vultures provide an essential ecosystem service through removal of carrion, but globally, many populations are collapsing and several species are threatened with extinction. Widespread declines in vulture populations could increase the availability of carrion to other organisms, but the ways facultative scavengers might respond to this increase have not been thoroughly explored. We aimed to determine whether facultative scavengers increase carrion consumption in the absence of vulture competition and whether they are capable of functionally replacing vultures in the removal of carrion biomass from the landscape. We experimentally excluded 65 rabbit carcasses from vultures during daylight hours and placed an additional 65 carcasses that were accessible to vultures in forested habitat in South Carolina, USA during summer (June-August). We used motion-activated cameras to compare carrion use by facultative scavenging species between the experimental and control carcasses. Scavenging by facultative scavengers did not increase in the absence of competition with vultures. We found no difference in scavenger presence between control carcasses and those from which vultures were excluded. Eighty percent of carcasses from which vultures were excluded were not scavenged by vertebrates, compared to 5% of carcasses that were accessible to vultures. At the end of the 7-day trials, there was a 10.1-fold increase in the number of experimental carcasses that were not fully scavenged compared to controls. Facultative scavengers did not functionally replace vultures during summer in our study. This finding may have been influenced by the time of the year in which the study took place, the duration of the trials, and the spacing of carcass sites. Our results suggest that under the warm and humid conditions of our study, facultative scavengers would not compensate for loss of vultures. Carcasses would persist longer in the environment and consumption of carrion would likely shift from vertebrates to decomposers. Such changes could have substantial implications for disease transmission, nutrient cycling, and ecosystem functioning.

Not in wilderness: African vulture strongholds remain in areas with high human density

Vultures constitute an important functional group in many ecosystems, providing crucial ecosystem services both in natural and humanized environments. These scavengers are facing massive declines worldwide, but in several African countries virtually nothing is known on populations’ status and threats, hampering the development of adequate conservation strategies. In Guinea-Bissau, globally important populations of Hooded Necrosyrtes monachus and African white-backed vultures Gyps africanus were recently reported. Using the country as a study area, we aim to characterize human-vulture interactions in West Africa applying a multidisciplinary approach. We assessed the status and distribution of vulture populations using data from 1711 km of roadside transects, examined predictors of their distribution, and produced a nationwide population estimate for the Hooded Vulture, using an innovative method based on the relationship between the size of human population in settlements and vulture numbers. We conducted 47 stakeholder interviews to assess perceived roles played by vultures, and to investigate potential anthropogenic threats. Hooded vultures were strongly associated with high human population densities, whereas no relation was found between African white-backed and Rüppell’s vultures and any of the tested predictors, which included cattle density, precipitation and Normalized Difference Vegetation Index, among others. We estimate a national population of 43347 Hooded vultures, the largest population reported in the species range. Respondents were generally aware of the services provided by vultures, especially waste and carcass removal, including in urban areas. Hunting for witchcraft and traditional medicine was the most frequently recognised threat, while poisoning was ranked as having the highest impact. We hypothesise that poisoning-related mortality may be affecting African white-backed and Rüppell’s vultures’ distribution and explain their scarcity in apparently highly suitable habitats. Our results suggest a mutualistic rather than a commensalistic relationship between vultures and humans, with important implications for designing and implementing conservation strategies.

Parasitaemia data and molecular characterization of Haemoproteus catharti from New World vultures (Cathartidae) reveals a novel clade of Haemosporida

Background

New World vultures (Cathartiformes: Cathartidae) are obligate scavengers comprised of seven species in five genera throughout the Americas. Of these, turkey vultures (Cathartes aura) and black vultures (Coragyps atratus) are the most widespread and, although ecologically similar, have evolved differences in morphology, physiology, and behaviour. Three species of haemosporidians have been reported in New World vultures to date: Haemoproteus catharti, Leucocytozoon toddi and Plasmodium elongatum, although few studies have investigated haemosporidian parasites in this important group of species. In this study, morphological and molecular methods were used to investigate the epidemiology and molecular biology of haemosporidian parasites of New World vultures in North America.

Methods

Blood and/or tissue samples were obtained from 162 turkey vultures and 95 black vultures in six states of the USA. Parasites were identified based on their morphology in blood smears, and sequences of the mitochondrial cytochrome b and nuclear adenylosuccinate lyase genes were obtained for molecular characterization.

Results

No parasites were detected in black vultures, whereas 24% of turkey vultures across all sampling locations were positive for H. catharti by blood smear analysis and/or PCR testing. The phylogenetic analysis of cytochrome b gene sequences revealed that H. catharti is closely related to MYCAMH1, a yet unidentified haemosporidian from wood storks (Mycteria americana) in southeastern USA and northern Brazil. Haemoproteus catharti and MYCAMH1 represent a clade that is unmistakably separate from all other Haemoproteus spp., being most closely related to Haemocystidium spp. from reptiles and to Plasmodium spp. from birds and reptiles.

Conclusions

Haemoproteus catharti is a widely-distributed parasite of turkey vultures in North America that is evolutionarily distinct from other haemosporidian parasites. These results reveal that the genetic diversity and evolutionary relationships of avian haemosporidians are still being uncovered, and future studies combining a comprehensive evaluation of morphological and life cycle characteristics with the analysis of multiple nuclear and mitochondrial genes will be useful to redefine the genus boundaries of these parasites and to re-evaluate the relationships amongst haemosporidians of birds, reptiles and mammals.

Electronic supplementary material

The online version of this article (10.1186/s12936-017-2165-5) contains supplementary material, which is available to authorized users.

Characterization of a de novo assembled transcriptome of the Common Blackbird (Turdus merula)

BACKGROUND

In recent years, next generation high throughput sequencing technologies have proven to be useful tools for investigations concerning the genomics or transcriptomics also of non-model species. Consequently, ornithologists have adopted these technologies and the respective bioinformatics tools to survey the genomes and transcriptomes of a few avian non-model species. The Common Blackbird is one of the most common bird species living in European cities, which has successfully colonized urban areas and for which no reference genome or transcriptome is publicly available. However, to target questions like genome wide gene expression analysis, a reference genome or transcriptome is needed.

METHODS

Therefore, in this study two Common Blackbirds were sacrificed, their mRNA was isolated and analyzed by RNA-Seq to de novo assemble a transcriptome and characterize it. Illumina reads (125 bp paired-end) and a Velvet/Oases pipeline led to 162,158 transcripts. For the annotation (using Blast+), an unfiltered protein database was used. SNPs were identified using SAMtools and BCFtools. Furthermore, mRNA from three single tissues (brain, heart and liver) of the same two Common Blackbirds were sequenced by Illumina (75 bp single-end reads). The draft transcriptome and the three single tissues were compared by their BLAST hits with the package VennDiagram in R.

RESULTS

Following the annotation against protein databases, we found evidence for 15,580 genes in the transcriptome (all well characterized hits after annotation). On 18% of the assembled transcripts, 144,742 SNPs were identified which are, consequently, 0.09% of all nucleotides in the assembled transcriptome. In the transcriptome and in the single tissues (brain, heart and liver), 10,182 shared genes were found.

DISCUSSION

Using a next-generation technology and bioinformatics tools, we made a first step towards the genomic investigation of the Common Blackbird. The de novo assembled transcriptome is usable for downstream analyses such as differential gene expression analysis and SNP identification. This study shows the importance of the approach to sequence single tissues to understand functions of tissues, proteins and the phenotype.

Myotis rufoniger genome sequence and analyses: M. rufoniger's genomic feature and the decreasing effective population size of Myotis bats

Myotis rufoniger is a vesper bat in the genus Myotis. Here we report the whole genome sequence and analyses of the M. rufoniger. We generated 124 Gb of short-read DNA sequences with an estimated genome size of 1.88 Gb at a sequencing depth of 66× fold. The sequences were aligned to M. brandtii bat reference genome at a mapping rate of 96.50% covering 95.71% coding sequence region at 10× coverage. The divergence time of Myotis bat family is estimated to be 11.5 million years, and the divergence time between M. rufoniger and its closest species M. davidii is estimated to be 10.4 million years. We found 1,239 function-altering M. rufoniger specific amino acid sequences from 929 genes compared to other Myotis bat and mammalian genomes. The functional enrichment test of the 929 genes detected amino acid changes in melanin associated DCT, SLC45A2, TYRP1, and OCA2 genes possibly responsible for the M. rufoniger's red fur color and a general coloration in Myotis. N6AMT1 gene, associated with arsenic resistance, showed a high degree of function alteration in M. rufoniger. We further confirmed that the M. rufoniger also has bat-specific sequences within FSHB, GHR, IGF1R, TP53, MDM2, SLC45A2, RGS7BP, RHO, OPN1SW, and CNGB3 genes that have already been published to be related to bat's reproduction, lifespan, flight, low vision, and echolocation. Additionally, our demographic history analysis found that the effective population size of Myotis clade has been consistently decreasing since ~30k years ago. M. rufoniger's effective population size was the lowest in Myotis bats, confirming its relatively low genetic diversity.

De novo transcriptome analysis shows differential expression of genes in salivary glands of edible bird's nest producing swiftlets

BACKGROUND

Edible bird's nest (EBN), produced from solidified saliva secretions of specific swiftlet species during the breeding season, is one of the most valuable animal by-products in the world. The composition and medicinal benefits of EBN have been extensively studied, however, genomic and transcriptomic studies of the salivary glands of these birds have not been conducted.

RESULTS

The study described the transcriptomes of salivary glands from three swiftlet species (28 samples) generated by RNASeq. A total of 14,835 annotated genes and 428 unmapped genes were cataloged. The current study investigated the genes and pathways that are associated with the development of salivary gland and EBN composition. Differential expression and pathway enrichment analysis indicated that the expression of CREB3L2 and several signaling pathways involved in salivary gland development, namely, the EGFR, BMP, and MAPK signaling pathways, were up-regulated in swiftlets producing white EBN (Aerodramus fuciphagus) and black EBN (Aerodramus maximus) compared with non-EBN-producing swiftlets (Apus affinis). Furthermore, MGAT, an essential gene for the biosynthesis of N-acetylneuraminic acid (sialic acid), was highly expressed in both white- and black-nest swiftlets compared to non-EBN-producing swiftlets. Interspecies comparison between Aerodramus fuciphagus and Aerodramus maximus indicated that the genes involved in N-acetylneuraminic and fatty acid synthesis were up-regulated in Aerodramus fuciphagus, while alanine and aspartate synthesis pathways were up-regulated in Aerodramus maximus. Furthermore, gender-based analysis revealed that N-glycan trimming pathway was significantly up-regulated in male Aerodramus fuciphagus from its natural habitat (cave) compared to their female counterpart.

CONCLUSIONS

Transcriptomic analysis of salivary glands of different swiftlet species reveal differential expressions of candidate genes that are involved in salivary gland development and in the biosynthesis of various bioactive compounds found in EBN.

Comparison of carnivore, omnivore, and herbivore mammalian genomes with a new leopard assembly

BACKGROUND

There are three main dietary groups in mammals: carnivores, omnivores, and herbivores. Currently, there is limited comparative genomics insight into the evolution of dietary specializations in mammals. Due to recent advances in sequencing technologies, we were able to perform in-depth whole genome analyses of representatives of these three dietary groups.

RESULTS

We investigated the evolution of carnivory by comparing 18 representative genomes from across Mammalia with carnivorous, omnivorous, and herbivorous dietary specializations, focusing on Felidae (domestic cat, tiger, lion, cheetah, and leopard), Hominidae, and Bovidae genomes. We generated a new high-quality leopard genome assembly, as well as two wild Amur leopard whole genomes. In addition to a clear contraction in gene families for starch and sucrose metabolism, the carnivore genomes showed evidence of shared evolutionary adaptations in genes associated with diet, muscle strength, agility, and other traits responsible for successful hunting and meat consumption. Additionally, an analysis of highly conserved regions at the family level revealed molecular signatures of dietary adaptation in each of Felidae, Hominidae, and Bovidae. However, unlike carnivores, omnivores and herbivores showed fewer shared adaptive signatures, indicating that carnivores are under strong selective pressure related to diet. Finally, felids showed recent reductions in genetic diversity associated with decreased population sizes, which may be due to the inflexible nature of their strict diet, highlighting their vulnerability and critical conservation status.

CONCLUSIONS

Our study provides a large-scale family level comparative genomic analysis to address genomic changes associated with dietary specialization. Our genomic analyses also provide useful resources for diet-related genetic and health research.