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 regions 1-3. Here, we address these issues by analyzing genomes of 363 bird species 4 (218 taxonomic families, 92% of total). Using intergenic regions and coalescent methods, we present a well-supported tree but also a remarkable degree of discordance. The tree confirms that Neoaves experienced rapid radiation at or near the Cretaceous-Paleogene (K-Pg) boundary. Sufficient loci rather than extensive taxon sampling were more effective in resolving difficult nodes. Remaining recalcitrant nodes involve species that challenge modeling due to extreme GC content, variable substitution rates, incomplete lineage sorting, or complex evolutionary events such as ancient hybridization. Assessment of the impacts 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 K-Pg extinction event, supporting the hypothesis that emerging ecological opportunities catalyzed the diversification of modern birds. The resulting phylogenetic estimate offers novel insights into the rapid radiation of modern birds and provides a taxon-rich backbone tree for future comparative studies.

Late Pleistocene stickleback environmental genomes reveal the chronology of freshwater adaptation

Directly observing the chronology and tempo of adaptation in response to ecological change is rarely possible in natural ecosystems. Sedimentary ancient DNA (sedaDNA) has been shown to be a tractable source of genome-scale data of long-dead organisms1,2,3 and to thereby potentially provide an understanding of the evolutionary histories of past populations.4,5 To date, time series of ecosystem biodiversity have been reconstructed from sedaDNA, typically using DNA metabarcoding or shotgun sequence data generated from less than 1 g of sediment.6,7 Here, we maximize sequence coverage by extracting DNA from ∼50× more sediment per sample than the majority of previous studies1,2,3 to achieve genotype resolution. From a time series of Late Pleistocene sediments spanning from a marine to freshwater ecosystem, we compare adaptive genotypes reconstructed from the environmental genomes of three-spined stickleback at key time points of this transition. We find a staggered temporal dynamic in which freshwater alleles at known loci of large effect in marine-freshwater divergence of three-spined stickleback (e.g., EDA)8 were already established during the brackish phase of the formation of the isolation basin. However, marine alleles were still detected across the majority of marine-freshwater divergence-associated loci, even after the complete isolation of the lake from marine ingression. Our retrospective approach to studying adaptation from environmental genomes of three-spined sticklebacks at the end of the last glacial period complements contemporary experimental approaches9,10,11 and highlights the untapped potential for retrospective "evolve and resequence" natural experiments using sedaDNA.

Host-gut microbiota interactions shape parasite infections in farmed Atlantic salmon

Animals and their associated microbiota share long evolutionary histories. However, it is not always clear how host genotype and microbiota interact to affect phenotype. We applied a hologenomic approach to explore how host-microbiota interactions shape lifetime growth and parasite infection in farmed Atlantic salmon (Salmo salar). Multi-omics data sets were generated from the guts of 460 salmon, 82% of which were naturally infected with an intestinal cestode. A single Mycoplasma bacterial strain, MAG01, dominated the gut metagenome of large, non-parasitized fish, consistent with previous studies showing high levels of Mycoplasma in the gut microbiota of healthy salmon. While small and/or parasitized salmon also had high abundance of MAG01, we observed increased alpha diversity in these individuals, driven by increased frequency of low-abundance Vibrionaceae and other Mycoplasma species that carried known virulence genes. Colonization by one of these cestode-associated Mycoplasma strains was associated with host individual genomic variation in long non-coding RNAs. Integrating the multi-omic data sets revealed coordinated changes in the salmon gut mRNA transcriptome and metabolome that correlated with shifts in the microbiota of smaller, parasitized fish. Our results suggest that the gut microbiota of small and/or parasitized fish is in a state of dysbiosis that partly depends on the host genotype, highlighting the value of using a hologenomic approach to incorporate the microbiota into the study of host-parasite dynamics.IMPORTANCEStudying host-microbiota interactions through the perspective of the hologenome is gaining interest across all life sciences. Intestinal parasite infections are a huge burden on human and animal health; however, there are few studies investigating the role of the hologenome during parasite infections. We address this gap in the largest multi-omics fish microbiota study to date using natural cestode infection of farmed Atlantic salmon. We find a clear association between cestode infection, salmon lifetime growth, and perturbation of the salmon gut microbiota. Furthermore, we provide the first evidence that the genetic background of the host may partly determine how the gut microbiota changes during parasite-associated dysbiosis. Our study therefore highlights the value of a hologenomic approach for gaining a more in-depth understanding of parasitism.

Exploring the microbial diversity of novel misos with metagenomics

Interest in fermented foods, especially plant-based ones, has increased considerably in the last decade. Miso-a Japanese paste traditionally fermented with soybeans, salt, and kōji (Aspergillus oryzae grown on grains or beans)-has gained attention among chefs for its rich flavour and versatility. Some chefs have even been experimenting with making novel misos with untraditional substrates to create new flavours. Such novel fermented foods also offer new scientific opportunities. To explore the microbial diversity of these new traditional foods, we sampled six misos made by the team at a leading restaurant called Noma in Copenhagen (Denmark), using yellow peas (including a nixtamalised treatment), lupin seeds, Swedish Vreta peas, grey peas, and Gotland lentils as substrates. All misos were made with the same recipe and fermented for 3 months at 28 °C. Samples were collected at the end of fermentation for subsequent shotgun metagenomic sequencing and a genome-resolved metagenomic analysis. The taxonomic profile of the samples revealed the presence of kōji mould (A. oryzae) and Bacillus amyloliquefaciens in all misos. Various species of the genera Latilactobacillus, Lactiplantibacillus, Pediococcus and Staphylococcus were also detected. The Metagenome-Assembled Genomes (MAGs) revealed genomic sequences belonging to 12 different species and functional analyses of these MAGs were performed. Notably, we detected the presence of Exiguobacterium-the first reported instance of the genus in miso-and Average Nucleotide Identity (ANI) analyses suggest a potentially new species. We hope these results will improve the scientific literature on misos and contribute to developing novel fermented plant-based foods.

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

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

A pangenome graph reference of 30 chicken genomes allows genotyping of large and complex structural variants

BACKGROUND

The red junglefowl, the wild outgroup of domestic chickens, has historically served as a reference for genomic studies of domestic chickens. These studies have provided insight into the etiology of traits of commercial importance. However, the use of a single reference genome does not capture diversity present among modern breeds, many of which have accumulated molecular changes due to drift and selection. While reference-based resequencing is well-suited to cataloging simple variants such as single-nucleotide changes and short insertions and deletions, it is mostly inadequate to discover more complex structural variation in the genome.

METHODS

We present a pangenome for the domestic chicken consisting of thirty assemblies of chickens from different breeds and research lines.

RESULTS

We demonstrate how this pangenome can be used to catalog structural variants present in modern breeds and untangle complex nested variation. We show that alignment of short reads from 100 diverse wild and domestic chickens to this pangenome reduces reference bias by 38%, which affects downstream genotyping results. This approach also allows for the accurate genotyping of a large and complex pair of structural variants at the K feathering locus using short reads, which would not be possible using a linear reference.

CONCLUSIONS

We expect that this new paradigm of genomic reference will allow better pinpointing of exact mutations responsible for specific phenotypes, which will in turn be necessary for breeding chickens that meet new sustainability criteria and are resilient to quickly evolving pathogen threats.

Redefining the Evolutionary History of the Rock Dove, Columba livia, Using Whole Genome Sequences

The domestic pigeon's exceptional phenotypic diversity was key in developing Darwin's Theory of Evolution and establishing the concept of artificial selection. However, unlike its domestic counterpart, its wild progenitor, the rock dove Columba livia has received considerably less attention. Therefore, questions regarding its domestication, evolution, taxonomy, and conservation status remain unresolved. We generated whole-genome sequencing data from 65 historical rock doves that represent all currently recognized subspecies and span the species' original geographic distribution. Our dataset includes 3 specimens from Darwin's collection, and the type specimens of 5 different taxa. We characterized their population structure, genomic diversity, and gene-flow patterns. Our results show the West African subspecies C. l. gymnocyclus is basal to rock doves and domestic pigeons, and suggests gene-flow between the rock dove's sister species C. rupestris, and the ancestor of rock doves after its split from West African populations. These genomes allowed us to propose a model for the evolution of the rock dove in light of the refugia theory. We propose that rock dove genetic diversity and introgression patterns derive from a history of allopatric cycles and dispersion waves during the Quaternary glacial and interglacial periods. To explore the rock dove domestication history, we combined our new dataset with available genomes from domestic pigeons. Our results point to at least 1 domestication event in the Levant that gave rise to all domestic breeds analysed in this study. Finally, we propose a species-level taxonomic arrangement to reflect the evolutionary history of the West African rock dove populations.

Ancient DNA reveals genetic admixture in China during tiger evolution

The tiger (Panthera tigris) is a charismatic megafauna species that originated and diversified in Asia and probably experienced population contraction and expansion during the Pleistocene, resulting in low genetic diversity of modern tigers. However, little is known about patterns of genomic diversity in ancient populations. Here we generated whole-genome sequences from ancient or historical (100-10,000 yr old) specimens collected across mainland Asia, including a 10,600-yr-old Russian Far East specimen (RUSA21, 8× coverage) plus six ancient mitogenomes, 14 South China tigers (0.1-12×) and three Caspian tigers (4-8×). Admixture analysis showed that RUSA21 clustered within modern Northeast Asian phylogroups and partially derived from an extinct Late Pleistocene lineage. While some of the 8,000-10,000-yr-old Russian Far East mitogenomes are basal to all tigers, one 2,000-yr-old specimen resembles present Amur tigers. Phylogenomic analyses suggested that the Caspian tiger probably dispersed from an ancestral Northeast Asian population and experienced gene flow from southern Bengal tigers. Lastly, genome-wide monophyly supported the South China tiger as a distinct subspecies, albeit with mitochondrial paraphyly, hence resolving its longstanding taxonomic controversy. The distribution of mitochondrial haplogroups corroborated by biogeographical modelling suggested that Southwest China was a Late Pleistocene refugium for a relic basal lineage. As suitable habitat returned, admixture between divergent lineages of South China tigers took place in Eastern China, promoting the evolution of other northern subspecies. Altogether, our analysis of ancient genomes sheds light on the evolutionary history of tigers and supports the existence of nine modern subspecies.

Unravelling animal-microbiota evolution on a chip

Whether and how microorganisms have shaped the evolution of their animal hosts is a major question in biology. Although many animal evolutionary processes appear to correlate with changes in their associated microbial communities, the mechanistic processes leading to these patterns and their causal relationships are still far from being resolved. Gut-on-a-chip models provide an innovative approach that expands beyond the potential of conventional microbiome profiling to study how different animals sense and react to microbes by comparing responses of animal intestinal tissue models to different microbial stimuli. This complementary knowledge can contribute to our understanding of how host genetic features facilitate or prevent different microbiomes from being assembled, and in doing so elucidate the role of host-microbiota interactions in animal evolution.

Assessing the genetic composition of cotton-top tamarins (Saguinus oedipus) before sweeping anthropogenic impact

During the last century, the critically endangered cotton-top tamarin (Saguinus oedipus) has been threatened by multiple anthropogenic factors that drastically affected their habitat and population size. As the genetic impact of these pressures is largely unknown, this study aimed to establish a genetic baseline with the use of temporal sampling to determine the genetic makeup before detrimental anthropogenic impact. Genomes were resequenced from a combination of historical museum samples and modern wild samples at low-medium coverage, to unravel how the cotton-top tamarin population structure and genomic diversity may have changed during this period. Our data suggest two populations can be differentiated, probably separated historically by the mountain ranges of the Paramillo Massif in Colombia. Although this population structure persists in the current populations, modern samples exhibit genomic signals consistent with recent inbreeding, such as long runs of homozygosity and a reduction in genome-wide heterozygosity especially in the greater northeast population. This loss is likely the consequence of the population reduction following the mass exportation of cotton-top tamarins for biomedical research in the 1960s, coupled with the habitat loss this species continues to experience. However, current populations have not experienced an increase in genetic load. We propose that the historical genetic baseline established in this study can be used to provide insight into alteration in the modern population influenced by a drastic reduction in population size as well as providing background information to be used for future conservation decision-making for the species.

Field and laboratory guidelines for reliable bioinformatic and statistical analysis of bacterial shotgun metagenomic data

Shotgun metagenomics is an increasingly cost-effective approach for profiling environmental and host-associated microbial communities. However, due to the complexity of both microbiomes and the molecular techniques required to analyze them, the reliability and representativeness of the results are contingent upon the field, laboratory, and bioinformatic procedures employed. Here, we consider 15 field and laboratory issues that critically impact downstream bioinformatic and statistical data processing, as well as result interpretation, in bacterial shotgun metagenomic studies. The issues we consider encompass intrinsic properties of samples, study design, and laboratory-processing strategies. We identify the links of field and laboratory steps with downstream analytical procedures, explain the means for detecting potential pitfalls, and propose mitigation measures to overcome or minimize their impact in metagenomic studies. We anticipate that our guidelines will assist data scientists in appropriately processing and interpreting their data, while aiding field and laboratory researchers to implement strategies for improving the quality of the generated results.

The ancestry and geographical origins of St Helena's liberated Africans

The island of St Helena played a crucial role in the suppression of the transatlantic slave trade. Strategically located in the middle of the South Atlantic, it served as a staging post for the Royal Navy and reception point for enslaved Africans who had been "liberated" from slave ships intercepted by the British. In total, St Helena received approximately 27,000 liberated Africans between 1840 and 1867. Written sources suggest that the majority of these individuals came from West Central Africa, but their precise origins are unknown. Here, we report the results of ancient DNA analyses that we conducted as part of a wider effort to commemorate St Helena's liberated Africans and to restore knowledge of their lives and experiences. We generated partial genomes (0.1-0.5×) for 20 individuals whose remains had been recovered during archaeological excavations on the island. We compared their genomes with genotype data for over 3,000 present-day individuals from 90 populations across sub-Saharan Africa and conclude that the individuals most likely originated from different source populations within the general area between northern Angola and Gabon. We also find that the majority (17/20) of the individuals were male, supporting a well-documented sex bias in the latter phase of the transatlantic slave trade. The study expands our understanding of St Helena's liberated African community and illustrates how ancient DNA analyses can be used to investigate the origins and identities of individuals whose lives were bound up in the story of slavery and its abolition.

Historic Sampling of a Vanishing Beast: Population Structure and Diversity in the Black Rhinoceros

The black rhinoceros (Diceros bicornis L.) is a critically endangered species historically distributed across sub-Saharan Africa. Hunting and habitat disturbance have diminished both its numbers and distribution since the 19th century, but a poaching crisis in the late 20th century drove them to the brink of extinction. Genetic and genomic assessments can greatly increase our knowledge of the species and inform management strategies. However, when a species has been severely reduced, with the extirpation and artificial admixture of several populations, it is extremely challenging to obtain an accurate understanding of historic population structure and evolutionary history from extant samples. Therefore, we generated and analyzed whole genomes from 63 black rhinoceros museum specimens collected between 1775 and 1981. Results showed that the black rhinoceros could be genetically structured into six major historic populations (Central Africa, East Africa, Northwestern Africa, Northeastern Africa, Ruvuma, and Southern Africa) within which were nested four further subpopulations (Maasailand, southwestern, eastern rift, and northern rift), largely mirroring geography, with a punctuated north-south cline. However, we detected varying degrees of admixture among groups and found that several geographical barriers, most prominently the Zambezi River, drove population discontinuities. Genomic diversity was high in the middle of the range and decayed toward the periphery. This comprehensive historic portrait also allowed us to ascertain the ancestry of 20 resequenced genomes from extant populations. Lastly, using insights gained from this unique temporal data set, we suggest management strategies, some of which require urgent implementation, for the conservation of the remaining black rhinoceros diversity.

The extinct Sicilian wolf shows a complex history of isolation and admixture with ancient dogs

The Sicilian wolf remained isolated in Sicily from the end of the Pleistocene until its extermination in the 1930s-1960s. Given its long-term isolation on the island and distinctive morphology, the genetic origin of the Sicilian wolf remains debated. We sequenced four nuclear genomes and five mitogenomes from the seven existing museum specimens to investigate the Sicilian wolf ancestry, relationships with extant and extinct wolves and dogs, and diversity. Our results show that the Sicilian wolf is most closely related to the Italian wolf but carries ancestry from a lineage related to European Eneolithic and Bronze Age dogs. The average nucleotide diversity of the Sicilian wolf was half of the Italian wolf, with 37-50% of its genome contained in runs of homozygosity. Overall, we show that, by the time it went extinct, the Sicilian wolf had high inbreeding and low-genetic diversity, consistent with a population in an insular environment.

Conservation implications of elucidating the Korean wolf taxonomic ambiguity through whole-genome sequencing

The taxonomic status of the now likely extirpated Korean Peninsula wolf has been extensively debated, with some arguing it represents an independent wolf lineage, Canis coreanus. To investigate the Korean wolf's genetic affiliations and taxonomic status, we sequenced and analysed the genomes of a Korean wolf dated to the beginning of the 20th century, and a captive wolf originally from the Pyongyang Central Zoo. Our results indicated that the Korean wolf bears similar genetic ancestry to other regional East Asian populations, therefore suggesting it is not a distinct taxonomic lineage. We identified regional patterns of wolf population structure and admixture in East Asia with potential conservation consequences in the Korean Peninsula and on a regional scale. We find that the Korean wolf has similar genomic diversity and inbreeding to other East Asian wolves. Finally, we show that, in contrast to the historical sample, the captive wolf is genetically more similar to wolves from the Tibetan Plateau; hence, Korean wolf conservation programmes might not benefit from the inclusion of this specimen.

Genomics of adaptive evolution in the woolly mammoth

Ancient genomes provide a tool to investigate the genetic basis of adaptations in extinct organisms. However, the identification of species-specific fixed genetic variants requires the analysis of genomes from multiple individuals. Moreover, the long-term scale of adaptive evolution coupled with the short-term nature of traditional time series data has made it difficult to assess when different adaptations evolved. Here, we analyze 23 woolly mammoth genomes, including one of the oldest known specimens at 700,000 years old, to identify fixed derived non-synonymous mutations unique to the species and to obtain estimates of when these mutations evolved. We find that at the time of its origin, the woolly mammoth had already acquired a broad spectrum of positively selected genes, including ones associated with hair and skin development, fat storage and metabolism, and immune system function. Our results also suggest that these phenotypes continued to evolve during the last 700,000 years, but through positive selection on different sets of genes. Finally, we also identify additional genes that underwent comparatively recent positive selection, including multiple genes related to skeletal morphology and body size, as well as one gene that may have contributed to the small ear size in Late Quaternary woolly mammoths.

Species-specific traits mediate avian demographic responses under past climate change

Anticipating species' responses to environmental change is a pressing mission in biodiversity conservation. Despite decades of research investigating how climate change may affect population sizes, historical context is lacking, and the traits that mediate demographic sensitivity to changing climate remain elusive. We use whole-genome sequence data to reconstruct the demographic histories of 263 bird species over the past million years and identify networks of interacting morphological and life history traits associated with changes in effective population size (N_e) in response to climate warming and cooling. Our results identify direct and indirect effects of key traits representing dispersal, reproduction and survival on long-term demographic responses to climate change, thereby highlighting traits most likely to influence population responses to ongoing climate warming.

Ancient DNA from a lost Negev Highlands desert grape reveals a Late Antiquity wine lineage

Recent excavations of Late Antiquity settlements in the Negev Highlands of southern Israel uncovered a society that established commercial-scale viticulture in an arid environment [D. Fuks et al., Proc. Natl. Acad. Sci. U.S.A. 117, 19780-19791 (2020)]. We applied target-enriched genome-wide sequencing and radiocarbon dating to examine grapevine pips that were excavated at three of these sites. Our analyses revealed centuries long and continuous grape cultivation in the Southern Levant. The genetically diverse pips also provided clues to ancient cultivation strategies aimed at improving agricultural productivity and ensuring food security. Applying genomic prediction analysis, a pip dated to the eighth century CE was determined to likely be from a white grape, to date the oldest to be identified. In a kinship analysis, another pip was found to be descendant from a modern Greek cultivar and was thus linked with several popular historic wines that were once traded across the Byzantine Empire. These findings shed light on historical Byzantine trading networks and on the genetic contribution of Levantine varieties to the classic Aegean landscape.

Intestinal epigenotype of Atlantic salmon (Salmo salar) associates with tenacibaculosis and gut microbiota composition

It remains a challenge to obtain the desired phenotypic traits in aquacultural production of Atlantic salmon, and part of the challenge might come from the effect that host-associated microorganisms have on the fish phenotype. To manipulate the microbiota towards the desired host traits, it is critical to understand the factors that shape it. The bacterial gut microbiota composition can vary greatly among fish, even when reared in the same closed system. While such microbiota differences can be linked to diseases, the molecular effect of disease on host-microbiota interactions and the potential involvement of epigenetic factors remain largely unknown. The aim of this study was to investigate the DNA methylation differences associated with a tenacibaculosis outbreak and microbiota displacement in the gut of Atlantic salmon. Using Whole Genome Bisulfite Sequencing (WGBS) of distal gut tissue from 20 salmon, we compared the genome-wide DNA methylation levels between uninfected individuals and sick fish suffering from tenacibaculosis and microbiota displacement. We discovered >19,000 differentially methylated cytosine sites, often located in differentially methylated regions, and aggregated around genes. The 68 genes connected to the most significant regions had functions related to the ulcerous disease such as epor and slc48a1a but also included prkcda and LOC106590732 whose orthologs are linked to microbiota changes in other species. Although the expression level was not analysed, our epigenetic analysis suggests specific genes potentially involved in host-microbiota interactions and more broadly it highlights the value off considering epigenetic factors in efforts to manipulate the microbiota of farmed fish.

Evolution of the germline mutation rate across vertebrates

The germline mutation rate determines the pace of genome evolution and is an evolving parameter itself¹. However, little is known about what determines its evolution, as most studies of mutation rates have focused on single species with different methodologies². Here we quantify germline mutation rates across vertebrates by sequencing and comparing the high-coverage genomes of 151 parent-offspring trios from 68 species of mammals, fishes, birds and reptiles. We show that the per-generation mutation rate varies among species by a factor of 40, with mutation rates being higher for males than for females in mammals and birds, but not in reptiles and fishes. The generation time, age at maturity and species-level fecundity are the key life-history traits affecting this variation among species. Furthermore, species with higher long-term effective population sizes tend to have lower mutation rates per generation, providing support for the drift barrier hypothesis³. The exceptionally high yearly mutation rates of domesticated animals, which have been continually selected on fecundity traits including shorter generation times, further support the importance of generation time in the evolution of mutation rates. Overall, our comparative analysis of pedigree-based mutation rates provides ecological insights on the mutation rate evolution in vertebrates.

On the origin and evolution of RNA editing in metazoans

Extensive adenosine-to-inosine (A-to-I) editing of nuclear-transcribed mRNAs is the hallmark of metazoan transcriptional regulation. Here, by profiling the RNA editomes of 22 species that cover major groups of Holozoa, we provide substantial evidence supporting A-to-I mRNA editing as a regulatory innovation originating in the last common ancestor of extant metazoans. This ancient biochemistry process is preserved in most extant metazoan phyla and primarily targets endogenous double-stranded RNA (dsRNA) formed by evolutionarily young repeats. We also find intermolecular pairing of sense-antisense transcripts as an important mechanism for forming dsRNA substrates for A-to-I editing in some but not all lineages. Likewise, recoding editing is rarely shared across lineages but preferentially targets genes involved in neural and cytoskeleton systems in bilaterians. We conclude that metazoan A-to-I editing might first emerge as a safeguard mechanism against repeat-derived dsRNA and was later co-opted into diverse biological processes due to its mutagenic nature.

Co-diversification of an intestinal Mycoplasma and its salmonid host

Understanding the evolutionary relationships between a host and its intestinal resident bacteria can transform how we understand adaptive phenotypic traits. The interplay between hosts and their resident bacteria inevitably affects the intestinal environment and, thereby, the living conditions of both the host and the microbiota. Thereby this co-existence likely influences the fitness of both bacteria and host. Whether this co-existence leads to evolutionary co-diversification in animals is largely unexplored, mainly due to the complexity of the environment and microbial communities and the often low host selection. We present the gut metagenome from wild Atlantic salmon (Salmo salar), a new wild organism model with an intestinal microbiota of low complexity and a well-described population structure, making it well-suited for investigating co-evolution. Our data reveal a strong host selection of a core gut microbiota dominated by a single Mycoplasma species. We found a clear co-diversification between the population structure of Atlantic salmon and nucleotide variability of the intestinal Mycoplasma populations conforming to expectations from co-evolution between host and resident bacteria. Our results show that the stable microbiota of Atlantic salmon has evolved with its salmonid host populations while potentially providing adaptive traits to the salmon host populations, including defence mechanisms, biosynthesis of essential amino acids, and metabolism of B vitamins. We highlight Atlantic salmon as a novel model for studying co-evolution between vertebrate hosts and their resident bacteria.

Genomic diversity and differentiation between island and mainland populations of white-tailed eagles (Haliaeetus albicilla)

Divergence in the face of high dispersal capabilities is a documented but poorly understood phenomenon. The white-tailed eagle (Haliaeetus albicilla) has a large geographic dispersal capability and should theoretically be able to maintain genetic homogeneity across its dispersal range. However, following analysis of the genomic variation of white-tailed eagles, from both historical and contemporary samples, clear signatures of ancient biogeographic substructure across Europe and the North-East Atlantic is observed. The greatest genomic differentiation was observed between island (Greenland and Iceland) and mainland (Denmark, Norway and Estonia) populations. The two island populations share a common ancestry from a single mainland population, distinct from the other sampled mainland populations, and despite the potential for high connectivity between Iceland and Greenland they are well separated from each other and are characterized by inbreeding and little variation. Temporal differences also highlight a pattern of regional populations persisting despite the potential for admixture. All sampled populations generally showed a decline in effective population size over time, which may have been shaped by four historical events: (1) Isolation of refugia during the last glacial period 110-115,000 years ago, (2) population divergence following the colonization of the deglaciated areas ~10,000 years ago, (3) human population expansion, which led to the settlement in Iceland ~1100 years ago, and (4) human persecution and exposure to toxic pollutants during the last two centuries.

The genetic history of Scandinavia from the Roman Iron Age to the present

We investigate a 2,000-year genetic transect through Scandinavia spanning the Iron Age to the present, based on 48 new and 249 published ancient genomes and genotypes from 16,638 modern individuals. We find regional variation in the timing and magnitude of gene flow from three sources: the eastern Baltic, the British-Irish Isles, and southern Europe. British-Irish ancestry was widespread in Scandinavia from the Viking period, whereas eastern Baltic ancestry is more localized to Gotland and central Sweden. In some regions, a drop in current levels of external ancestry suggests that ancient immigrants contributed proportionately less to the modern Scandinavian gene pool than indicated by the ancestry of genomes from the Viking and Medieval periods. Finally, we show that a north-south genetic cline that characterizes modern Scandinavians is mainly due to the differential levels of Uralic ancestry and that this cline existed in the Viking Age and possibly earlier.

A simplified protocol for DNA extraction from FTA cards for faecal microbiome studies

As metagenomic studies continue to increase in size and complexity, they are often required to incorporate data from geographically isolated locations or longitudinal time samples. This represents a technical challenge, given that many of the commonly used methods used for sample collection, storage, and DNA extraction are sensitive to differences related to the time, storage and chemistry involved. FTA cards have been previously proposed as a simple, reliable and cost-efficient method for the preservation of animal faecal microbiomes. In this study, we report a simplified extraction methodology for recovering microbiome DNA from faeces stored on FTA cards and compare its performance to a common alternative means of characterising such microbiomes; namely, immediate freezing of the faeces followed by DNA extraction using the Qiagen PowerSoil DNA isolation kit. Our results show that overall the application of our simplified DNA extraction methodology yields microbial community results that have higher diversity and an expanded core microbiome than that found using the PowerSoil methodology. This suggests that the FTA card extraction method presented here is a viable alternative for metagenomic studies using faecal material when traditional freeze-based storage methods are not feasible.

Enriching captivity conditions with natural elements does not prevent the loss of wild‐like gut microbiota but shapes its compositional variation in two small mammals

As continued growth in gut microbiota studies in captive and model animals elucidates the importance of their role in host biology, further pursuit of how to retain a wild‐like microbial community is becoming increasingly important to obtain representative results from captive animals. In this study, we assessed how the gut microbiota of two wild‐caught small mammals, namely Crocidura russula (Eulipotyphla, insectivore) and Apodemus sylvaticus (Rodentia, omnivore), changed when bringing them into captivity. We analyzed fecal samples of 15 A. sylvaticus and 21 C. russula, immediately after bringing them into captivity and 5 weeks later, spread over two housing treatments: a “natural” setup enriched with elements freshly collected from nature and a “laboratory” setup with sterile artificial elements. Through sequencing of the V3–V4 region of the 16S recombinant RNA gene, we found that the initial microbial diversity dropped during captivity in both species, regardless of treatment. Community composition underwent a change of similar magnitude in both species and under both treatments. However, we did observe that the temporal development of the gut microbiome took different trajectories (i.e., changed in different directions) under different treatments, particularly in C. russula, suggesting that C. russula may be more susceptible to environmental change. The results of this experiment do not support the use of microbially enriched environments to retain wild‐like microbial diversities and compositions, yet show that specific housing conditions can significantly affect the drift of microbial communities under captivity.

The population genomic legacy of the second plague pandemic

Human populations have been shaped by catastrophes that may have left long-lasting signatures in their genomes. One notable example is the second plague pandemic that entered Europe in ca. 1,347 CE and repeatedly returned for over 300 years, with typical village and town mortality estimated at 10%-40%.¹ It is assumed that this high mortality affected the gene pools of these populations. First, local population crashes reduced genetic diversity. Second, a change in frequency is expected for sequence variants that may have affected survival or susceptibility to the etiologic agent (Yersinia pestis).² Third, mass mortality might alter the local gene pools through its impact on subsequent migration patterns. We explored these factors using the Norwegian city of Trondheim as a model, by sequencing 54 genomes spanning three time periods: (1) prior to the plague striking Trondheim in 1,349 CE, (2) the 17^th-19^th century, and (3) the present. We find that the pandemic period shaped the gene pool by reducing long distance immigration, in particular from the British Isles, and inducing a bottleneck that reduced genetic diversity. Although we also observe an excess of large F_ST values at multiple loci in the genome, these are shaped by reference biases introduced by mapping our relatively low genome coverage degraded DNA to the reference genome. This implies that attempts to detect selection using ancient DNA (aDNA) datasets that vary by read length and depth of sequencing coverage may be particularly challenging until methods have been developed to account for the impact of differential reference bias on test statistics.

Novel strategies to improve chicken performance and welfare by unveiling host-microbiota interactions through hologenomics

Fast optimisation of farming practices is essential to meet environmental sustainability challenges. Hologenomics, the joint study of the genomic features of animals and the microbial communities associated with them, opens new avenues to obtain in-depth knowledge on how host-microbiota interactions affect animal performance and welfare, and in doing so, improve the quality and sustainability of animal production. Here, we introduce the animal trials conducted with broiler chickens in the H2020 project HoloFood, and our strategy to implement hologenomic analyses in light of the initial results, which despite yielding negligible effects of tested feed additives, provide relevant information to understand how host genomic features, microbiota development dynamics and host-microbiota interactions shape animal welfare and performance. We report the most relevant results, propose hypotheses to explain the observed patterns, and outline how these questions will be addressed through the generation and analysis of animal-microbiota multi-omic data during the HoloFood project.

Uncovering the genomic basis of an extraordinary plant invasion

Invasive species are a key driver of the global biodiversity crisis, but the drivers of invasiveness, including the role of pathogens, remain debated. We investigated the genomic basis of invasiveness in Ambrosia artemisiifolia (common ragweed), introduced to Europe in the late 19th century, by resequencing 655 ragweed genomes, including 308 herbarium specimens collected up to 190 years ago. In invasive European populations, we found selection signatures in defense genes and lower prevalence of disease-inducing plant pathogens. Together with temporal changes in population structure associated with introgression from closely related Ambrosia species, escape from specific microbial enemies likely favored the plant's remarkable success as an invasive species.

Reference genomes for conservation

High-quality reference genomes for non-model species can benefit conservation.

Grey wolf genomic history reveals a dual ancestry of dogs

The grey wolf (Canis lupus) was the first species to give rise to a domestic population, and they remained widespread throughout the last Ice Age when many other large mammal species went extinct. Little is known, however, about the history and possible extinction of past wolf populations or when and where the wolf progenitors of the present-day dog lineage (Canis familiaris) lived^1–8. Here we analysed 72 ancient wolf genomes spanning the last 100,000 years from Europe, Siberia and North America. We found that wolf populations were highly connected throughout the Late Pleistocene, with levels of differentiation an order of magnitude lower than they are today. This population connectivity allowed us to detect natural selection across the time series, including rapid fixation of mutations in the gene IFT88 40,000–30,000 years ago. We show that dogs are overall more closely related to ancient wolves from eastern Eurasia than to those from western Eurasia, suggesting a domestication process in the east. However, we also found that dogs in the Near East and Africa derive up to half of their ancestry from a distinct population related to modern southwest Eurasian wolves, reflecting either an independent domestication process or admixture from local wolves. None of the analysed ancient wolf genomes is a direct match for either of these dog ancestries, meaning that the exact progenitor populations remain to be located.

DNA from ancient wolves spanning 100,000 years sheds light on wolves’ evolutionary history and the genomic origin of dogs.

Phylogenomics of the world's otters

Comparative whole-genome analyses hold great power to illuminate commonalities and differences in the evolution of related species that share similar ecologies. The mustelid subfamily Lutrinae includes 13 currently recognized extant species of otters,^1-5 a semiaquatic group whose evolutionary history is incompletely understood. We assembled a dataset comprising 24 genomes from all living otter species, 14 of which were newly sequenced. We used this dataset to infer phylogenetic relationships and divergence times, to characterize patterns of genome-wide genealogical discordance, and to investigate demographic history and current genomic diversity. We found that genera Lutra, Aonyx, Amblonyx, and Lutrogale form a coherent clade that should be synonymized under Lutra, simplifying the taxonomic structure of the subfamily. The poorly known tropical African Aonyx congicus and the more widespread Aonyx capensis were found to be reciprocally monophyletic (having diverged 440,000 years ago), supporting the validity of the former as a distinct species. We observed variable changes in effective population sizes over time among otters within and among continents, although several species showed similar trends of expansions and declines during the last 100,000 years. This has led to different levels of genomic diversity assessed by overall heterozygosity, genome-wide SNV density, and run of homozygosity burden. Interestingly, there were cases in which diversity metrics were consistent with the current threat status (mostly based on census size), highlighting the potential of genomic data for conservation assessment. Overall, our results shed light on otter evolutionary history and provide a framework for further in-depth comparative genomic studies targeting this group.

The genomic basis of the plant island syndrome in Darwin's giant daisies

The repeated, rapid and often pronounced patterns of evolutionary divergence observed in insular plants, or the ‘plant island syndrome’, include changes in leaf phenotypes, growth, as well as the acquisition of a perennial lifestyle. Here, we sequence and describe the genome of the critically endangered, Galápagos-endemic species Scalesia atractyloides Arnot., obtaining a chromosome-resolved, 3.2-Gbp assembly containing 43,093 candidate gene models. Using a combination of fossil transposable elements, k-mer spectra analyses and orthologue assignment, we identify the two ancestral genomes, and date their divergence and the polyploidization event, concluding that the ancestor of all extant Scalesia species was an allotetraploid. There are a comparable number of genes and transposable elements across the two subgenomes, and while their synteny has been mostly conserved, we find multiple inversions that may have facilitated adaptation. We identify clear signatures of selection across genes associated with vascular development, growth, adaptation to salinity and flowering time, thus finding compelling evidence for a genomic basis of the island syndrome in one of Darwin’s giant daisies.

Many island plant species share a syndrome of characteristic phenotype and life history. Cerca et al. find the genomic basis of the plant island syndrome in one of Darwin’s giant daisies, while separating ancestral genomes in a chromosome-resolved polyploid assembly.

Travel Tales of a Worldwide Weed: Genomic Signatures of Plantago major L. Reveal Distinct Genotypic Groups With Links to Colonial Trade Routes

Retracing pathways of historical species introductions is fundamental to understanding the factors involved in the successful colonization and spread, centuries after a species' establishment in an introduced range. Numerous plants have been introduced to regions outside their native ranges both intentionally and accidentally by European voyagers and early colonists making transoceanic journeys; however, records are scarce to document this. We use genotyping-by-sequencing and genotype-likelihood methods on the selfing, global weed, Plantago major, collected from 50 populations worldwide to investigate how patterns of genomic diversity are distributed among populations of this global weed. Although genomic differentiation among populations is found to be low, we identify six unique genotype groups showing very little sign of admixture and low degree of outcrossing among them. We show that genotype groups are latitudinally restricted, and that more than one successful genotype colonized and spread into the introduced ranges. With the exception of New Zealand, only one genotype group is present in the Southern Hemisphere. Three of the most prevalent genotypes present in the native Eurasian range gave rise to introduced populations in the Americas, Africa, Australia, and New Zealand, which could lend support to the hypothesis that P. major was unknowlingly dispersed by early European colonists. Dispersal of multiple successful genotypes is a likely reason for success. Genomic signatures and phylogeographic methods can provide new perspectives on the drivers behind the historic introductions and the successful colonization of introduced species, contributing to our understanding of the role of genomic variation for successful establishment of introduced taxa.

Extreme freeze-tolerance in cryophilic tardigrades relies on controlled ice formation but does not involve significant change in transcription

Subzero temperatures are among the most significant factors defining the distribution of organisms, yet, certain taxa have evolved to overcome this barrier. The microscopic tardigrades are among the most freeze-tolerant animals, with selected species reported to survive milli-Kelvin temperatures. Here, we estimate survival of fully hydrated eutardigrades of the species Ramazzottius varieornatus following exposures to -20 °C and  -80 °C as well as -196 °C with or without initial cooling to -80 °C. The tardigrades easily survive these temperatures, yet with a significant decrease in viability following rapid cooling by direct exposure to -196 °C. Hence, post-freeze recovery of R. varieornatus seems to rely on cooling rate and thus controlled ice formation. Cryophilic organisms are renowned for having cold-active enzymes that secure appropriate reaction rates at low temperatures. Hence, extreme freeze-tolerance in R. varieornatus could potentially involve syntheses of cryoprotectants and de novo transcription. We therefore generated a reference transcriptome for this cryophilic R. varieornatus population and explored for differential gene expression patterns following cooling to -80 °C as compared to active 5 °C controls. Specifically, we tested for fast transcription potentially occurring within 25 min of cooling from room temperature to a supercooling point of ca. -20 °C, at which the tardigrades presumably freeze and enter into the ametabolic state of cryobiosis. Our analyses revealed no evidence for differential gene expression. We, therefore, conclude that extreme freeze-tolerance in R. varieornatus relies on controlled extracellular freezing with any freeze-tolerance related genes being constitutively expressed.

Disentangling host-microbiota complexity through hologenomics

Research on animal-microbiota interactions has become a central topic in biological sciences because of its relevance to basic eco-evolutionary processes and applied questions in agriculture and health. However, animal hosts and their associated microbial communities are still seldom studied in a systemic fashion. Hologenomics, the integrated study of the genetic features of a eukaryotic host alongside that of its associated microbes, is becoming a feasible - yet still underexploited - approach that overcomes this limitation. Acknowledging the biological and genetic properties of both hosts and microbes, along with the advantages and disadvantages of implemented techniques, is essential for designing optimal studies that enable some of the major questions in biology to be addressed.

Differential expression profiling of heat stressed tardigrades reveals major shift in the transcriptome

Tardigrades are renowned for their extreme stress tolerance, which includes the ability to endure complete desiccation, high levels of radiation and very low sub-zero temperatures. Nevertheless, tardigrades appear to be vulnerable to high temperatures and thus the potential effects of global warming. Here, we provide the first analysis of transcriptome data obtained from heat stressed specimens of the eutardigrade Ramazzottius varieornatus, with the aim of providing new insights into the molecular processes affected by high temperatures. Specifically, we compare RNA-seq datasets obtained from active, heat-exposed (35 °C) tardigrades to that of active controls kept at 5 °C. Our data reveal a surprising shift in transcription, involving 9634 differentially expressed transcripts, corresponding to >35% of the transcriptome. The latter data are in striking contrast to the hitherto observed constitutive expression underlying tardigrade extreme stress tolerance and entrance into the latent state of life, known as cryptobiosis. Thus, when examining the molecular response, heat-stress appears to be more stressful for R. varieornatus than extreme conditions, such as desiccation or freezing. A gene ontology analysis reveals that the heat stress response involves a change in transcription and presumably translation, including an adjustment of metabolism, and, putatively, preparation for encystment and subsequent diapause. Among the differentially expressed transcripts we find heat-shock proteins as well as the eutardigrade specific proteins (CAHS, SAHS, MAHS, RvLEAM, and Dsup). The latter proteins thus seem to contribute to a general stress response, and may not be directly related to cryptobiosis.

The Australian dingo is an early offshoot of modern breed dogs

Dogs are uniquely associated with human dispersal and bring transformational insight into the domestication process. Dingoes represent an intriguing case within canine evolution being geographically isolated for thousands of years. Here, we present a high-quality de novo assembly of a pure dingo (CanFam_DDS). We identified large chromosomal differences relative to the current dog reference (CanFam3.1) and confirmed no expanded pancreatic amylase gene as found in breed dogs. Phylogenetic analyses using variant pairwise matrices show that the dingo is distinct from five breed dogs with 100% bootstrap support when using Greenland wolf as the outgroup. Functionally, we observe differences in methylation patterns between the dingo and German shepherd dog genomes and differences in serum biochemistry and microbiome makeup. Our results suggest that distinct demographic and environmental conditions have shaped the dingo genome. In contrast, artificial human selection has likely shaped the genomes of domestic breed dogs after divergence from the dingo.

A multi-omics approach unravels metagenomic and metabolic alterations of a probiotic and synbiotic additive in rainbow trout (Oncorhynchus mykiss)

BACKGROUND

Animal protein production is increasingly looking towards microbiome-associated services such as the design of new and better probiotic solutions to further improve gut health and production sustainability. Here, we investigate the functional effects of bacteria-based pro- and synbiotic feed additives on microbiome-associated functions in relation to growth performance in the commercially important rainbow trout (Oncorhynchus mykiss). We combine complementary insights from multiple omics datasets from gut content samples, including 16S bacterial profiling, whole metagenomes, and untargeted metabolomics, to investigate bacterial metagenome-assembled genomes (MAGs) and their molecular interactions with host metabolism.

RESULTS

Our findings reveal that (I) feed additives changed the microbiome and that rainbow trout reared with feed additives had a significantly reduced relative abundance of the salmonid related Candidatus Mycoplasma salmoninae in both the mid and distal gut content, (II) genome resolved metagenomics revealed that alterations of microbial arginine biosynthesis and terpenoid backbone synthesis pathways were directly associated with the presence of Candidatus Mycoplasma salmoninae, and (III) differences in the composition of intestinal microbiota among feed types were directly associated with significant changes of the metabolomic landscape, including lipids and lipid-like metabolites, amino acids, bile acids, and steroid-related metabolites.

CONCLUSION

Our results demonstrate how the use of multi-omics to investigate complex host-microbiome interactions enable us to better evaluate the functional potential of probiotics compared to studies that only measure overall growth performance or that only characterise the microbial composition in intestinal environments. Video Abstract.

To bee, or not to bee? One leg is the question

Temporal genomic studies that utilise museum insects are invaluable for understanding changes in ecological processes in which insects are essential, such as wild and agricultural pollination, seed dispersal, nutrient cycling, and food web architecture, to name a few. However, given such analyses come at the cost of physical damage to museum specimens required for such work, there is a natural interest in the development and/or application of methods to minimise the damage incurred. We explored the efficacy of a recently published single stranded library construction protocol, on DNA extracted from single legs taken from eight dry-preserved historic bee specimens collected 150 years ago. Specifically, the DNA was extracted using a 'minimally destructive' method that leaves the samples' exterior intact. Our sequencing data revealed not only that the endogenous DNA recovered from some of the samples was at a relatively high level (up to 58%), but that the complexity of the libraries was sufficient in the best samples to theoretically allow deeper sequencing to a predicted level of 69x genome coverage. As such, these combined protocols offer the possibility to generate sequencing data at levels that are suitable for many common evolutionary genomic analyses, while simultaneously minimising the damage conferred to the valuable dried museum bee samples. Furthermore, we anticipate that these methods may have much wider application on many other invertebrate taxa stored in a similar way. We hope that the results from this research may be able to contribute to the increased willingness of museums to loan much needed dry-preserved insects for future genomic studies.

Gut Microbiota Linked with Reduced Fear of Humans in Red Junglefowl Has Implications for Early Domestication

Domestication of animals can lead to profound phenotypic modifications within short evolutionary time periods, and for many species behavioral selection is likely at the forefront of this process. Animal studies have strongly implicated that the gut microbiome plays a major role in host behavior and cognition through the microbiome-gut-brain axis. Consequently, herein, it is hypothesized that host gut microbiota may be one of the earliest phenotypes to change as wild animals were domesticated. Here, the gut microbiome community in two selected lines of red junglefowl that are selected for either high or low fear of humans up to eight generations is examined. Microbiota profiles reveal taxonomic differences in gut bacteria known to produce neuroactive compounds between the two selection lines. Gut-brain module analysis by means of genome-resolved metagenomics identifies enrichment in the microbial synthesis and degradation potential of metabolites associated with fear extinction and reduces anxiety-like behaviors in low fear fowls. In contrast, high fear fowls are enriched in gut-brain modules from the butyrate and glutamate pathways, metabolites associated with fear conditioning. Overall, the results identify differences in the composition and functional potential of the gut microbiota across selection lines that may provide insights into the mechanistic explanations of the domestication process.

Kouprey (Bos sauveli) genomes unveil polytomic origin of wild Asian Bos

The evolution of the genera Bos and Bison, and the nature of gene flow between wild and domestic species, is poorly understood, with genomic data of wild species being limited. We generated two genomes from the likely extinct kouprey (Bos sauveli) and analyzed them alongside other Bos and Bison genomes. We found that B. sauveli possessed genomic signatures characteristic of an independent species closely related to Bos javanicus and Bos gaurus. We found evidence for extensive incomplete lineage sorting across the three species, consistent with a polytomic diversification of the major ancestry in the group, potentially followed by secondary gene flow. Finally, we detected significant gene flow from an unsampled Asian Bos-like source into East Asian zebu cattle, demonstrating both that the full genomic diversity and evolutionary history of the Bos complex has yet to be elucidated and that museum specimens and ancient DNA are valuable resources to do so.

•

We generated two genomes from the likely extinct kouprey (Bos sauveli)

•

Extensive mt and nuclear-genome-wide incomplete lineage sorting across wild Asian Bos

•

Initial polytomic diversification of the wild Asian Bos—kouprey, banteng, and gaur

The role of the gut microbiota in the dietary niche expansion of fishing bats

Background

Due to its central role in animal nutrition, the gut microbiota is likely a relevant factor shaping dietary niche shifts. We analysed both the impact and contribution of the gut microbiota to the dietary niche expansion of the only four bat species that have incorporated fish into their primarily arthropodophage diet.

Results

We first compared the taxonomic and functional features of the gut microbiota of the four piscivorous bats to that of 11 strictly arthropodophagous species using 16S rRNA targeted amplicon sequencing. Second, we increased the resolution of our analyses for one of the piscivorous bat species, namely Myotis capaccinii, and analysed multiple populations combining targeted approaches with shotgun sequencing. To better understand the origin of gut microorganisms, we also analysed the gut microbiota of their fish prey (Gambusia holbrooki). Our analyses showed that piscivorous bats carry a characteristic gut microbiota that differs from that of their strict arthropodophagous counterparts, in which the most relevant bacteria have been directly acquired from their fish prey. This characteristic microbiota exhibits enrichment of genes involved in vitamin biosynthesis, as well as complex carbohydrate and lipid metabolism, likely providing their hosts with an enhanced capacity to metabolise the glycosphingolipids and long-chain fatty acids that are particularly abundant in fish.

Conclusions

Our results depict the gut microbiota as a relevant element in facilitating the dietary transition from arthropodophagy to piscivory.

Supplementary Information

The online version contains supplementary material available at 10.1186/s42523-021-00137-w.