Bat Biology, Genomes, and the Bat1K Project: To Generate Chromosome-Level Genomes for All Living Bat Species

Comparative transcriptome analysis reveals molecular adaptations underlying distinct immunity and inverted resting posture in bats

Understanding how natural selection shapes unique traits in mammals is a central topic in evolutionary biology. The mammalian order Chiroptera (bats) is attractive for biologists as well as the general public due to their specific traits of extraordinary immunity and inverted resting posture. However, genomic resources for bats that occupy key phylogenetic positions are not sufficient, which hinders comprehensive investigation of the molecular mechanisms underpinning the origin of specific traits in bats. Here, we sequenced the transcriptomes of five bats that are phylogenetically divergent and occupy key positions in the phylogenetic tree of bats. In combination with the available genomes of 19 bats and 21 other mammals, we built a database consisting of 10,918 one-to-one ortholog genes and reconstructed phylogenetic relationships of these mammals. We found that genes related to immunity, bone remodeling and cardiovascular system are targets of natural selection along the ancestral branch of bats. Further analyses revealed that the T cell receptor signaling pathway involved in immune adaptation is specifically enriched in bats. Moreover, molecular adaptations of bone remodeling, cardiovascular system, and balance sensing may help to explain the reverted resting posture in bats. Our study provides valuable transcriptome resources, enabling us to tentatively identify genetic changes associated with bat-specific traits. This work is among the first to advance our understanding of molecular underpinnings of inverted resting posture in bats, which could provide insight into healthcare applications such as hypertension in humans. This article is protected by copyright. All rights reserved.

Molecular, ecological, and behavioral drivers of the bat-virus relationship

Bats perform important ecological roles in our ecosystem. However, recent studies have demonstrated that bats are reservoirs of emerging viruses that have spilled over into humans and agricultural animals to cause severe diseases. These viruses include Hendra and Nipah paramyxoviruses, Ebola and Marburg filoviruses, and coronaviruses that are closely related to severe acute respiratory syndrome coronavirus (SARS-CoV), Middle East respiratory syndrome coronavirus (MERS-CoV), and the recently emerged SARS-CoV-2. Intriguingly, bats that are naturally or experimentally infected with these viruses do not show clinical signs of disease. Here we have reviewed ecological, behavioral, and molecular factors that may influence the ability of bats to harbor viruses. We have summarized known zoonotic potential of bat-borne viruses and stress on the need for further studies to better understand the evolutionary relationship between bats and their viruses, along with discovering the intrinsic and external factors that facilitate the successful spillover of viruses from bats.

Optimising predictive models to prioritise viral discovery in zoonotic reservoirs

Despite the global investment in One Health disease surveillance, it remains difficult and costly to identify and monitor the wildlife reservoirs of novel zoonotic viruses. Statistical models can guide sampling target prioritisation, but the predictions from any given model might be highly uncertain; moreover, systematic model validation is rare, and the drivers of model performance are consequently under-documented. Here, we use the bat hosts of betacoronaviruses as a case study for the data-driven process of comparing and validating predictive models of probable reservoir hosts. In early 2020, we generated an ensemble of eight statistical models that predicted host-virus associations and developed priority sampling recommendations for potential bat reservoirs of betacoronaviruses and bridge hosts for SARS-CoV-2. During a time frame of more than a year, we tracked the discovery of 47 new bat hosts of betacoronaviruses, validated the initial predictions, and dynamically updated our analytical pipeline. We found that ecological trait-based models performed well at predicting these novel hosts, whereas network methods consistently performed approximately as well or worse than expected at random. These findings illustrate the importance of ensemble modelling as a buffer against mixed-model quality and highlight the value of including host ecology in predictive models. Our revised models showed an improved performance compared with the initial ensemble, and predicted more than 400 bat species globally that could be undetected betacoronavirus hosts. We show, through systematic validation, that machine learning models can help to optimise wildlife sampling for undiscovered viruses and illustrates how such approaches are best implemented through a dynamic process of prediction, data collection, validation, and updating.

Characterization of Pipistrellus pygmaeus Bat Virome from Sweden

Increasing amounts of data indicate that bats harbor a higher viral diversity relative to other mammalian orders, and they have been recognized as potential reservoirs for pathogenic viruses, such as the Hendra, Nipah, Marburg, and SARS-CoV viruses. Here, we present the first viral metagenomic analysis of Pipistrellus pygmaeus from Uppsala, Sweden. Total RNA was extracted from the saliva and feces of individual bats and analyzed using Illumina sequencing. The results identified sequences related to 51 different viral families, including vertebrate, invertebrate, and plant viruses. These viral families include Coronaviridae, Picornaviridae, Dicistroviridae, Astroviridae, Hepeviridae, Reoviridae, Botourmiaviridae, Lispviridae, Totiviridae, Botoumiaviridae, Parvoviridae, Retroviridae, Adenoviridae, and Partitiviridae, as well as different unclassified viruses. We further characterized three near full-length genome sequences of bat coronaviruses. A phylogenetic analysis showed that these belonged to alphacoronaviruses with the closest similarity (78–99% at the protein level) to Danish and Finnish bat coronaviruses detected in Pipistrellus and Myotis bats. In addition, the full-length and the near full-length genomes of picornavirus were characterized. These showed the closest similarity (88–94% at the protein level) to bat picornaviruses identified in Chinese bats. Altogether, the results of this study show that Swedish Pipistrellus bats harbor a great diversity of viruses, some of which are closely related to mammalian viruses. This study expands our knowledge on the bat population virome and improves our understanding of the evolution and transmission of viruses among bats and to other species.

Genomewide Analysis and Biological Characterization of Cathelicidins with Potent Antimicrobial Activity and Low Cytotoxicity from Three Bat Species

Cathelicidins are potent antimicrobial peptides with broad spectrum antimicrobial activity in many vertebrates and an important component of the innate immune system. However, our understanding of the genetic variations and biological characteristics of bat cathelicidins is limited. In this study, we performed genome-level analysis of the antimicrobial peptide cathelicidins from seven bat species in the six families, listed 19 cathelicidin-like sequences, and showed that the number of functional cathelicidin genes differed among bat species. Based on the identified biochemical characteristics of bat cathelicidins, three cathelicidins, HA-CATH (from Hipposideros armiger), ML-CATH (from Myotis lucifugus), and PD-CATH (from Phyllostomus discolor), with clear antimicrobial signatures were chemically synthesized and evaluated antimicrobial activity. HA-CATH showed narrow-spectrum antibacterial activity against a panel of 12 reference bacteria, comprising 6 Gram-negative and 6 Gram-positive strains. However, ML-CATH and PD-CATH showed potent antibacterial activity against a broad spectrum of Gram-negative and Gram-positive bacteria with minimum inhibitory concentration (MIC) of 1 and 3 μg/mL, respectively, against Staphylococcus aureus. ML-CATH and PD-CATH also showed antifungal activities against Candida albicans and Cryptococcus cuniculi with MIC of 5 to 40 μg/mL, respectively, and 80% inhibition of the metabolism of Mucor hiemalis hyphae at 80 μg/mL, while displaying minimal cytotoxicity to HaCaT cells. Taken together, although the spectrum and efficacy of bat cathelicidins were species-dependent, the antimicrobial activity of ML-CATH and PD-CATH was comparable to that of other highly active cathelicidins in vertebrates while having negligible cytotoxicity to mammalian cells. ML-CATH and PD-CATH can be exploited as promising candidates for the development of antimicrobial therapeutics.

Interaction between Old World fruit bats and humans: From large scale ecosystem services to zoonotic diseases

The Old World tropical and subtropical frugivorous bat genus Rousettus (Pteropodidae) contains species with broad distributions, as well as those occurring in restricted geographical areas, particularly islands. Herein we review the role of Rousettus as a keystone species from a global "One Health" approach and related to ecosystem functioning, zoonotic disease and public health. Rousettus are efficient at dispersing seeds and pollinating flowers; their role in forest regeneration is related to their ability to fly considerable distances during nightly foraging bouts and their relatively small body size, which allows them to access fruits in forested areas with closed vegetation. Rousettus are also reservoirs for various groups of pathogens (viruses, bacteria, fungi, protozoa), which, by definition, are infectious agents causing disease. The study of day roosts of different species of Rousettus and the successful establishment of captive breeding colonies have provided important details related to the infection dynamics of their associated pathogens. Large-scale conversion of forested areas into agricultural landscapes has increased contact between humans and Rousettus, therefore augmenting the chances of infectious agent spillover. Many crucial scientific details are still lacking related to members of this genus, which have direct bearing on the prevention of emerging disease outbreaks, as well as the conservation of these bats. The public should be better informed on the capacity of fruit bats as keystone species for large scale forest regeneration and in spreading pathogens. Precise details on the transmission of zoonotic diseases of public health importance associated with Rousettus should be given high priority.

Gene prediction in the immunoglobulin loci

The V(D)J recombination process rearranges the variable (V), diversity (D), and joining (J) genes in the immunoglobulin loci to generate antibody repertoires. Annotation of these loci across various species and predicting the V, D, and J genes (IG genes) is critical for studies of the adaptive immune system. However, since the standard gene finding algorithms are not suitable for predicting IG genes, they have been semi-manually annotated in very few species. We developed the IGDetective algorithm for predicting IG genes and applied it to species with the assembled IG loci. IGDetective generated the first large collection of IG genes across many species and enabled their evolutionary analysis, including the analysis of the "bat IG diversity" hypothesis. This analysis revealed extremely conserved V genes in evolutionary distant species indicating that these genes may be subjected to the same selective pressure, e.g., pressure driven by common pathogens. IGDetective also revealed extremely diverged V genes and a new family of evolutionary conserved V genes in bats with unusual noncanonical cysteines. Moreover, in difference from all other previously reported antibodies, these cysteines are located within complementarity-determining regions. Since cysteines form disulfide bonds, we hypothesize that these cysteine-rich V genes might generate antibodies with noncanonical conformations and could potentially form a unique part of the immune repertoire in bats. We also analyzed the diversity landscape of the recombination signal sequences and revealed their features that trigger the high/low usage of the IG genes.

Emerging viruses: Cross-species transmission of Coronaviruses, Filoviruses, Henipaviruses and Rotaviruses from bats

Emerging infectious diseases, especially if caused by bat-borne viruses, significantly affect public health and the global economy. There is an urgent need to understand the mechanism of interspecies transmission, particularly to humans. Viral genetics; host factors, including polymorphisms in the receptors; and ecological, environmental, and population dynamics are major parameters to consider. Here, we describe the taxonomy, geographic distribution, and unique traits of bats associated with their importance as virus reservoirs. Then, we summarize the origin, intermediate hosts, and the current understanding of interspecies transmission of Middle East respiratory syndrome coronavirus (MERS-CoV), severe acute respiratory syndrome coronavirus (SARS-CoV), SARS-CoV-2, Nipah, Hendra, Ebola, Marburg virus, and rotaviruses. Finally, the molecular interactions of viral surface proteins with host cell receptors are examined, and a comparison of these interactions in humans, intermediate hosts, and bats is conducted. This uncovers adaptive mutations in virus spike protein that facilitate cross-species transmission and risk factors associated with the emergence of novel viruses from bats.

Inferring mammalian tissue-specific regulatory conservation by predicting tissue-specific differences in open chromatin

BACKGROUND

Evolutionary conservation is an invaluable tool for inferring functional significance in the genome, including regions that are crucial across many species and those that have undergone convergent evolution. Computational methods to test for sequence conservation are dominated by algorithms that examine the ability of one or more nucleotides to align across large evolutionary distances. While these nucleotide alignment-based approaches have proven powerful for protein-coding genes and some non-coding elements, they fail to capture conservation of many enhancers, distal regulatory elements that control spatial and temporal patterns of gene expression. The function of enhancers is governed by a complex, often tissue- and cell type-specific code that links combinations of transcription factor binding sites and other regulation-related sequence patterns to regulatory activity. Thus, function of orthologous enhancer regions can be conserved across large evolutionary distances, even when nucleotide turnover is high.

RESULTS

We present a new machine learning-based approach for evaluating enhancer conservation that leverages the combinatorial sequence code of enhancer activity rather than relying on the alignment of individual nucleotides. We first train a convolutional neural network model that can predict tissue-specific open chromatin, a proxy for enhancer activity, across mammals. Next, we apply that model to distinguish instances where the genome sequence would predict conserved function versus a loss of regulatory activity in that tissue. We present criteria for systematically evaluating model performance for this task and use them to demonstrate that our models accurately predict tissue-specific conservation and divergence in open chromatin between primate and rodent species, vastly out-performing leading nucleotide alignment-based approaches. We then apply our models to predict open chromatin at orthologs of brain and liver open chromatin regions across hundreds of mammals and find that brain enhancers associated with neuron activity have a stronger tendency than the general population to have predicted lineage-specific open chromatin.

CONCLUSION

The framework presented here provides a mechanism to annotate tissue-specific regulatory function across hundreds of genomes and to study enhancer evolution using predicted regulatory differences rather than nucleotide-level conservation measurements.

From bats to pangolins: new insights into species differences in the structure and function of the immune system

Species differences in the structure and function of the immune system of laboratory animals are known to exist and have been reviewed extensively. However, the number and diversity of wild and exotic species, along with their associated viruses, that come into contact with humans has increased worldwide sometimes with lethal consequences. Far less is known about the immunobiology of these exotic and wild species. Data suggest that species differences of the mechanisms of inflammation, innate immunity and adaptive immunity are all involved in the establishment and maintenance of viral infections across reservoir hosts. The current review attempts to collect relevant data concerning the basics of innate and adaptive immune functions of exotic and wild species followed by identification of those differences that may play a role in the maintenance of viral infections in reservoir hosts.

Positively selected genes in the hoary bat (Lasiurus cinereus) lineage: prominence of thymus expression, immune and metabolic function, and regions of ancient synteny

Background

Bats of the genus Lasiurus occur throughout the Americas and have diversified into at least 20 species among three subgenera. The hoary bat (Lasiurus cinereus) is highly migratory and ranges farther across North America than any other wild mammal. Despite the ecological importance of this species as a major insect predator, and the particular susceptibility of lasiurine bats to wind turbine strikes, our understanding of hoary bat ecology, physiology, and behavior remains poor.

Methods

To better understand adaptive evolution in this lineage, we used whole-genome sequencing to identify protein-coding sequence and explore signatures of positive selection. Gene models were predicted with Maker and compared to seven well-annotated and phylogenetically representative species. Evolutionary rate analysis was performed with PAML.

Results

Of 9,447 single-copy orthologous groups that met evaluation criteria, 150 genes had a significant excess of nonsynonymous substitutions along the L. cinereus branch (P < 0.001 after manual review of alignments). Selected genes as a group had biased expression, most strongly in thymus tissue. We identified 23 selected genes with reported immune functions as well as a divergent paralog of Steep1 within suborder Yangochiroptera. Seventeen genes had roles in lipid and glucose metabolic pathways, partially overlapping with 15 mitochondrion-associated genes; these adaptations may reflect the metabolic challenges of hibernation, long-distance migration, and seasonal variation in prey abundance. The genomic distribution of positively selected genes differed significantly from background expectation by discrete Kolmogorov-Smirnov test (P < 0.001). Remarkably, the top three physical clusters all coincided with islands of conserved synteny predating Mammalia, the largest of which shares synteny with the human cat-eye critical region (CECR) on 22q11. This observation coupled with the expansion of a novel Tbx1-like gene family may indicate evolutionary innovation during pharyngeal arch development: both the CECR and Tbx1 cause dosage-dependent congenital abnormalities in thymus, heart, and head, and craniodysmorphy is associated with human orthologs of other positively selected genes as well.

Fast and compact matching statistics analytics

MOTIVATION

Fast, lightweight methods for comparing the sequence of ever larger assembled genomes from ever growing databases are increasingly needed in the era of accurate long reads and pan-genome initiatives. Matching statistics is a popular method for computing whole-genome phylogenies and for detecting structural rearrangements between two genomes, since it is amenable to fast implementations that require a minimal setup of data structures. However, current implementations use a single core, take too much memory to represent the result, and do not provide efficient ways to analyze the output in order to explore local similarities between the sequences.

RESULTS

We develop practical tools for computing matching statistics between large-scale strings, and for analyzing its values, faster and using less memory than the state-of-the-art. Specifically, we design a parallel algorithm for shared-memory machines that computes matching statistics 30 times faster with 48 cores in the cases that are most difficult to parallelize. We design a lossy compression scheme that shrinks the matching statistics array to a bitvector that takes from 0.8 to 0.2 bits per character, depending on the dataset and on the value of a threshold, and that achieves 0.04 bits per character in some variants. And we provide efficient implementations of range-maximum and range-sum queries that take a few tens of milliseconds while operating on our compact representations, and that allow computing key local statistics about the similarity between two strings. Our toolkit makes construction, storage and analysis of matching statistics arrays practical for multiple pairs of the largest genomes available today, possibly enabling new applications in comparative genomics.

AVAILABILITY AND IMPLEMENTATION

Our C/C++ code is available at https://github.com/odenas/indexed_ms under GPL-3.0. The data underlying this article are available in NCBI Genome at https://www.ncbi.nlm.nih.gov/genome and in the International Genome Sample Resource (IGSR) at https://www.internationalgenome.org.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Alphavirus Identification in Neotropical Bats

Alphaviruses (Togaviridae) are arthropod-borne viruses responsible for several emerging diseases, maintained in nature through transmission between hematophagous arthropod vectors and susceptible vertebrate hosts. Although bats harbor many species of viruses, their role as reservoir hosts in emergent zoonoses has been verified only in a few cases. With bats being the second most diverse order of mammals, their implication in arbovirus infections needs to be elucidated. Reports on arbovirus infections in bats are scarce, especially in South American indigenous species. In this work, we report the genomic detection and identification of two different alphaviruses in oral swabs from bats captured in Northern Uruguay. Phylogenetic analysis identified Río Negro virus (RNV) in two different species: Tadarida brasiliensis (n = 6) and Myotis spp. (n = 1) and eastern equine encephalitis virus (EEEV) in Myotis spp. (n = 2). Previous studies of our group identified RNV and EEEV in mosquitoes and horse serology, suggesting that they may be circulating in enzootic cycles in our country. Our findings reveal that bats can be infected by these arboviruses and that chiropterans could participate in the viral natural cycle as virus amplifiers or dead-end hosts. Further studies are warranted to elucidate the role of these mammals in the biological cycle of these alphaviruses in Uruguay.

The European Variation Archive: a FAIR resource of genomic variation for all species

The European Variation Archive (EVA; https://www.ebi.ac.uk/eva/) is a resource for sharing all types of genetic variation data (SNPs, indels, and structural variants) for all species. The EVA was created in 2014 to provide FAIR access to genetic variation data and has since grown to be a primary resource for genomic variants hosting >3 billion records. The EVA and dbSNP have established a compatible global system to assign unique identifiers to all submitted genetic variants. The EVA is active within the Global Alliance of Genomics and Health (GA4GH), maintaining, contributing and implementing standards such as VCF, Refget and Variant Representation Specification (VRS). In this article, we describe the submission and permanent accessioning services along with the different ways the data can be retrieved by the scientific community.

Bats and their vital ecosystem services: a global review

Bats play crucial ecosystem services as seed dispersers, pollinators, controllers of insects, and nutrient recyclers. However, there has not been a thorough global review evaluating these roles in bats across all biogeographical regions of the world. We reviewed the literature published during the last two decades and identified 283 relevant studies: 78 dealt with the control of potential insect pests by bats, 80 related to the suppression of other arthropods, 60 on the dispersal of native or endemic seeds, 11 dealt with the dispersal of seeds of introduced plants, 29 on the pollination of native or endemic plants, 1 study on pollination of introduced plants, and 24 on the use of guano as fertilizer. Our literature search showed that queries combining the terms "seed dispersal," "insectivorous bats," "nectarivorous bats," "use of guano," and "ecosystem services" returned 577 studies, but half were experimental in nature. We found that the evaluation of ecosystem services by bats has been mostly conducted in the Neotropical and Palearctic regions. To detect differences across relevant studies, and to explain trends in the study of ecosystem services provided by bats, we performed generalized linear mixed models (GLMM) fitted with a Poisson distribution to analyze potential differences among sampling methods. We identified 409 bat species that provide ecosystem services, 752 insect species consumed by bats and 549 plant species either dispersed or pollinated by bats. Our review summarizes the importance of conserving bat populations and the ecological services they provide, which is especially important during the current pandemic.

Phenotyping in the era of genomics: MaTrics—a digital character matrix to document mammalian phenotypic traits

A new and uniquely structured matrix of mammalian phenotypes, MaTrics (Mammalian Traits for Comparative Genomics) in a digital form is presented. By focussing on mammalian species for which genome assemblies are available, MaTrics provides an interface between mammalogy and comparative genomics.

MaTrics was developed within a project aimed to find genetic causes of phenotypic traits of mammals using Forward Genomics. This approach requires genomes and comprehensive and recorded information on homologous phenotypes that are coded as discrete categories in a matrix. MaTrics is an evolving online resource providing information on phenotypic traits in numeric code; traits are coded either as absent/present or with several states as multistate. The state record for each species is linked to at least one reference (e.g., literature, photographs, histological sections, CT scans, or museum specimens) and so MaTrics contributes to digitalization of museum collections. Currently, MaTrics covers 147 mammalian species and includes 231 characters related to structure, morphology, physiology, ecology, and ethology and available in a machine actionable NEXUS-format*. Filling MaTrics revealed substantial knowledge gaps, highlighting the need for phenotyping efforts. Studies based on selected data from MaTrics and using Forward Genomics identified associations between genes and certain phenotypes ranging from lifestyles (e.g., aquatic) to dietary specializations (e.g., herbivory, carnivory). These findings motivate the expansion of phenotyping in MaTrics by filling research gaps and by adding taxa and traits. Only databases like MaTrics will provide machine actionable information on phenotypic traits, an important limitation to genomics. MaTrics is available within the data repository Morph·D·Base (www.morphdbase.de).

Retroviruses of Bats: a Threat Waiting in the Wings?

Bats are infamous reservoirs of deadly human viruses. While retroviruses, such as the human immunodeficiency virus (HIV), are among the most significant of virus families that have jumped from animals into humans, whether bat retroviruses have the potential to infect and cause disease in humans remains unknown. Recent reports of retroviruses circulating in bat populations builds on two decades of research describing the fossil records of retroviral sequences in bat genomes and of viral metagenomes extracted from bat samples. The impact of the global COVID-19 pandemic demands that we pay closer attention to viruses hosted by bats and their potential as a zoonotic threat. Here we review current knowledge of bat retroviruses and explore the question of whether they represent a threat to humans.

Evolution of pathogen tolerance and emerging infections: A missing experimental paradigm

Researchers worldwide are repeatedly warning us against future zoonotic diseases resulting from humankind's insurgence into natural ecosystems. The same zoonotic pathogens that cause severe infections in a human host frequently fail to produce any disease outcome in their natural hosts. What precise features of the immune system enable natural reservoirs to carry these pathogens so efficiently? To understand these effects, we highlight the importance of tracing the evolutionary basis of pathogen tolerance in reservoir hosts, while drawing implications from their diverse physiological and life-history traits, and ecological contexts of host-pathogen interactions. Long-term co-evolution might allow reservoir hosts to modulate immunity and evolve tolerance to zoonotic pathogens, increasing their circulation and infectious period. Such processes can also create a genetically diverse pathogen pool by allowing more mutations and genetic exchanges between circulating strains, thereby harboring rare alive-on-arrival variants with extended infectivity to new hosts (i.e., spillover). Finally, we end by underscoring the indispensability of a large multidisciplinary empirical framework to explore the proposed link between evolved tolerance, pathogen prevalence, and spillover in the wild.

Genomic landscape of mitochondrial DNA insertions in 23 bat genomes: characteristics, loci, phylogeny, and polymorphism

The transfer of mitochondrial DNA to the nuclear genome gives rise to the nuclear DNA sequences of mitochondrial origin (NUMTs), considered as a driving force in genome evolution. In this study, NUMTs in 23 bat genomes were investigated and compared systematically. The results showed that NUMTs existed in 22 genomes except for Noctilio leporinus, suggesting that mitochondrial fragment insertion in the nuclear genome was a common event in bat genomes. However, remarkable variations in NUMTs number, cumulative length, and proportion in the nuclear genome were discovered across bat species. Further orthologous NUMT loci analysis of the Phyllostomidae family indicated that the NUMTs insertion events in bat genomes were homoplasy-free. The NUMTs were mainly inserted into the intergenic regions, particularly, co-localized with repetitive sequences (especially transposable elements). However, several NUMTs were inserted into genes, some of which were in the exon region of functional genes. One NUMT in the genome of Pteropus alecto surprisingly matched with cDNA of ATP8B3 that provided evidence of NUMTs with coding function. Phylogenic analysis on NUMTs originating from COXI and COXII loci highlighted 2 clusters of Yinpterochiroptera and Yangochiroptera for Chiroptera. Seven NUMTs from Rhinolophus ferrumequinum were amplified, and the sequencing results confirmed the reliability of the NUMT analysis. One of them was polymorphic for the presence or absence of the NUMT insertion, and each genotype of NUMT loci showed a distinct regional distribution pattern. The information obtained in this study provides novel insights into the NUMT organization and features in bat genomes and establishes a basis for further studying of the evolution of bat species.

Novel de Novo Genome of Cynopterus brachyotis Reveals Evolutionarily Abrupt Shifts in Gene Family Composition across Fruit Bats

Major novel physiological or phenotypic adaptations often require accompanying modifications at the genic level. Conversely, the detection of considerable contractions and/or expansions of gene families can be an indicator of fundamental but unrecognized physiological change. We sequenced a novel fruit bat genome (Cynopterus brachyotis) and adopted a comparative approach to reconstruct the evolution of fruit bats, mapping contractions and expansions of gene families along their evolutionary history. Despite a radical change in life history as compared with other bats (e.g., loss of echolocation, large size, and frugivory), fruit bats have undergone surprisingly limited change in their genic composition, perhaps apart from a potentially novel gene family expansion relating to telomere protection and longevity. In sharp contrast, within fruit bats, the new Cynopterus genome bears the signal of unusual gene loss and gene family contraction, despite its similar morphology and lifestyle to two other major fruit bat lineages. Most missing genes are regulatory, immune-related, and olfactory in nature, illustrating the diversity of genomic strategies employed by bats to contend with responses to viral infection and olfactory requirements. Our results underscore that significant fluctuations in gene family composition are not always associated with obvious examples of novel physiological and phenotypic adaptations but may often relate to less-obvious shifts in immune strategies.

Adaptation-Driven Evolution of Sirtuin 1 (SIRT1), a Key Regulator of Metabolism and Aging, in Marmot Species

The sirtuin protein family plays a role in the lifespan of various species and is involved in numerous key metabolic processes. To understand the evolutionary role of sirtuins in marmots, a long-living rodent species group with remarkable metabolic shutdown during hibernation, we conducted a phylogeny-based substitution rate analysis of coding genes based on genetic information of seven marmot species. We show that sirtuin 1 (SIRT1) has evolved under positive selection in the marmot lineage. We pinpoint three amino acid changes in four different marmot species that underlie the signal of positive selection and that may favor increased longevity in marmots. Based on a computational structural analysis we can show that all three substitutions affect the secondary structure of the same region in human SIRT1. We propose that the identified region is close to the catalytic domain and that the potential structural changes may impact the catalytic activity of the enzyme and therefore might be playing a functional role in marmot's extended lifespan and metabolic shutdown.

Metagenomic analysis of fecal and tissue samples from 18 endemic bat species in Switzerland revealed a diverse virus composition including potentially zoonotic viruses

Many recent disease outbreaks in humans had a zoonotic virus etiology. Bats in particular have been recognized as reservoirs to a large variety of viruses with the potential to cross-species transmission. In order to assess the risk of bats in Switzerland for such transmissions, we determined the virome of tissue and fecal samples of 14 native and 4 migrating bat species. In total, sequences belonging to 39 different virus families, 16 of which are known to infect vertebrates, were detected. Contigs of coronaviruses, adenoviruses, hepeviruses, rotaviruses A and H, and parvoviruses with potential zoonotic risk were characterized in more detail. Most interestingly, in a ground stool sample of a Vespertilio murinus colony an almost complete genome of a Middle East respiratory syndrome-related coronavirus (MERS-CoV) was detected by Next generation sequencing and confirmed by PCR. In conclusion, bats in Switzerland naturally harbour many different viruses. Metagenomic analyses of non-invasive samples like ground stool may support effective surveillance and early detection of viral zoonoses.

Decoding bat immunity: the need for a coordinated research approach

Understanding antiviral immune responses in bats, which are reservoirs for many emerging viruses, could aid the response to future epidemics. Here, we discuss five key areas in which greater consensus among the bat research community is necessary to drive breakthroughs in the field.

The COVID-19 pandemic has stressed the importance of understanding species such as bats that can serve as reservoirs for emerging viral threats. Here, Wang and colleagues call for greater consensus among the bat immunology community in five key areas.

VIDHOP, viral host prediction with deep learning

Motivation

Zoonosis, the natural transmission of infections from animals to humans, is a far-reaching global problem. The recent outbreaks of Zikavirus, Ebolavirus, and Coronavirus are examples of viral zoonosis, which occur more frequently due to globalization. In case of a virus outbreak, it is helpful to know which host organism was the original carrier of the virus to prevent further spreading of viral infection. Recent approaches aim to predict a viral host based on the viral genome, often in combination with the potential host genome and arbitrarily selected features. These methods are limited in the number of different hosts they can predict or the accuracy of the prediction.

Results

Here, we present a fast and accurate deep learning approach for viral host prediction, which is based on the viral genome sequence only. We tested our deep neural network (DNN) on three different virus species (influenza A virus, rabies lyssavirus, rotavirus A). We achieved for each virus species an AUC between 0.93 and 0.98, allowing highly accurate predictions while using only fractions (100-400 bp) of the viral genome sequences. We show that deep neural networks are suitable to predict the host of a virus, even with a limited amount of sequences and highly unbalanced available data. The trained DNNs are the core of our virus-host prediction tool VIDHOP (VIrus Deep learning HOst Prediction). VIDHOP also allows the user to train and use models for other viruses.

Availability

VIDHOP is freely available under https://github.com/flomock/vidhop

Supplementary information

Available at DOI 10.17605/OSF.IO/UXT7

Genetic variation between long-lived versus short-lived bats illuminates the molecular signatures of longevity

Bats are the longest-lived mammals given their body size with majority of species exhibiting exceptional longevity. However, there are some short-lived species that do not exhibit extended lifespans. Here we conducted a comparative genomic and transcriptomic study on long-lived Myotis myotis (maximum lifespan = 37.1 years) and short-lived Molossus molossus (maximum lifespan = 5.6 years) to ascertain the genetic difference underlying their divergent longevities. Genome-wide selection tests on 12,467 single-copy genes between M. myotis and M. molossus revealed only three genes (CCDC175, FATE1 and MLKL) that exhibited significant positive selection. Although 97.96% of 12,467 genes underwent purifying selection, we observed a significant heterogeneity in their expression patterns. Using a linear mixed model, we obtained expression of 2,086 genes that may truly represent the genetic difference between M. myotis and M. molossus. Expression analysis indicated that long-lived M. myotis exhibited a transcriptomic profile of enhanced DNA repair and autophagy pathways, compared to M. molossus. Further investigation of the longevity-associated genes suggested that long-lived M. myotis have naturally evolved a diminished anti-longevity transcriptomic profile. Together with observations from other long-lived species, our results suggest that heightened DNA repair and autophagy activity may represent a universal mechanism to achieve longevity in long-lived mammals.

Evolutionary Constraint on Visual and Nonvisual Mammalian Opsins

Animals have evolved light-sensitive G protein-coupled receptors, known as opsins, to detect coherent and ambient light for visual and nonvisual functions. These opsins have evolved to satisfy the particular lighting niches of the organisms that express them. While many unique patterns of evolution have been identified in mammals for rod and cone opsins, far less is known about the atypical mammalian opsins. Using genomic data from over 400 mammalian species from 22 orders, unique patterns of evolution for each mammalian opsins were identified, including photoisomerases, RGR-opsin (RGR) and peropsin (RRH), as well as atypical opsins, encephalopsin (OPN3), melanopsin (OPN4), and neuropsin (OPN5). The results demonstrate that OPN5 and rhodopsin show extreme conservation across all mammalian lineages. The cone opsins, SWS1 and LWS, and the nonvisual opsins, OPN3 and RRH, demonstrate a moderate degree of sequence conservation relative to other opsins, with some instances of lineage-specific gene loss. Finally, the photoisomerase, RGR, and the best-studied atypical opsin, OPN4, have high sequence diversity within mammals. These conservation patterns are maintained in human populations. Importantly, all mammalian opsins retain key amino acid residues important for conjugation to retinal-based chromophores, permitting light sensitivity. These patterns of evolution are discussed along with known functions of each atypical opsin, such as in circadian or metabolic physiology, to provide insight into the observed patterns of evolutionary constraint.

Towards complete and error-free genome assemblies of all vertebrate species

High-quality and complete reference genome assemblies are fundamental for the application of genomics to biology, disease, and biodiversity conservation. However, such assemblies are available for only a few non-microbial species1-4. To address this issue, the international Genome 10K (G10K) consortium5,6 has worked over a five-year period to evaluate and develop cost-effective methods for assembling highly accurate and nearly complete reference genomes. Here we present lessons learned from generating assemblies for 16 species that represent six major vertebrate lineages. We confirm that long-read sequencing technologies are essential for maximizing genome quality, and that unresolved complex repeats and haplotype heterozygosity are major sources of assembly error when not handled correctly. Our assemblies correct substantial errors, add missing sequence in some of the best historical reference genomes, and reveal biological discoveries. These include the identification of many false gene duplications, increases in gene sizes, chromosome rearrangements that are specific to lineages, a repeated independent chromosome breakpoint in bat genomes, and a canonical GC-rich pattern in protein-coding genes and their regulatory regions. Adopting these lessons, we have embarked on the Vertebrate Genomes Project (VGP), an international effort to generate high-quality, complete reference genomes for all of the roughly 70,000 extant vertebrate species and to help to enable a new era of discovery across the life sciences.

Bioinformatic Analysis of Genome-Predicted Bat Cathelicidins

Bats are unique in their potential to serve as reservoir hosts for intracellular pathogens. Recently, the impact of COVID-19 has relegated bats from biomedical darkness to the frontline of public health as bats are the natural reservoir of many viruses, including SARS-Cov-2. Many bat genomes have been sequenced recently, and sequences coding for antimicrobial peptides are available in the public databases. Here we provide a structural analysis of genome-predicted bat cathelicidins as components of their innate immunity. A total of 32 unique protein sequences were retrieved from the NCBI database. Interestingly, some bat species contained more than one cathelicidin. We examined the conserved cysteines within the cathelin-like domain and the peptide portion of each sequence and revealed phylogenetic relationships and structural dissimilarities. The antibacterial, antifungal, and antiviral activity of peptides was examined using bioinformatic tools. The peptides were modeled and subjected to docking analysis with the region binding domain (RBD) region of the SARS-CoV-2 Spike protein. The appearance of multiple forms of cathelicidins verifies the complex microbial challenges encountered by these species. Learning more about antiviral defenses of bats and how they drive virus evolution will help scientists to investigate the function of antimicrobial peptides in these species.

Evolution of mammalian longevity: age-related increase in autophagy in bats compared to other mammals

Autophagy maintains cellular homeostasis and its dysfunction has been implicated in aging. Bats are the longest-lived mammals for their size, but the molecular mechanisms underlying their extended healthspan are not well understood. Here, drawing on >8 years of mark-recapture field studies, we report the first longitudinal analysis of autophagy regulation in bats. Mining of published population level aging blood transcriptomes (M. myotis, mouse and human) highlighted a unique increase of autophagy related transcripts with age in bats, but not in other mammals. This bat-specific increase in autophagy transcripts was recapitulated by the western blot determination of the autophagy marker, LC3II/I ratio, in skin primary fibroblasts (Myotis myotis,Pipistrellus kuhlii, mouse), that also showed an increase with age in both bat species. Further phylogenomic selection pressure analyses across eutherian mammals (n=70 taxa; 274 genes) uncovered 10 autophagy-associated genes under selective pressure in bat lineages. These molecular adaptations potentially mediate the exceptional age-related increase of autophagy signalling in bats, which may contribute to their longer healthspans.

Ecological and Conservation Significance of Herpesvirus Infection in Neotropical Bats

Bats are the second most diverse order of mammals and key species for ecosystem functioning, providing a wide range of ecosystem services, from pest control to seed dispersal. Chiropterans are known for hosting a large diversity of viruses, in some cases with little or no effect to their health. Here, we report on the results of a screening for DNA (Herpesviridae) and RNA viruses (Rhabdovirus and Pneumovirus), finding a high prevalence and wide diversity of both Beta- and Gamma-Herpesvirus in insectivorous and hematophagous bats of the southern cone of South America. Our findings suggest that bats in the southern neotropics harbor a high diversity of herpesviruses and, at least in some cases, the viral community in the bat species is more strongly associated with ecological traits of the hosts, rather than their taxonomy. The presence of a separate clade into the Gammaherpesvirinae subfamily in the common vampire bat suggests the independent circulation of herpesviruses in hematophagous and insectivorous bats and highlights the properness of these viruses to track vampire bats' population structure for rabies studies. Hence, we suggest that as other pathogens viruses may be used to track the population dynamics of their hosts, including movement and demographics.

Making a bat: The developmental basis of bat evolution

Bats are incredibly diverse, both morphologically and taxonomically. Bats are the only mammalian group to have achieved powered flight, an adaptation that is hypothesized to have allowed them to colonize various and diverse ecological niches. However, the lack of fossils capturing the transition from terrestrial mammal to volant chiropteran has obscured much of our understanding of bat evolution. Over the last 20 years, the emergence of evo-devo in non-model species has started to fill this gap by uncovering some developmental mechanisms at the origin of bat diversification. In this review, we highlight key aspects of studies that have used bats as a model for morphological adaptations, diversification during adaptive radiations, and morphological novelty. To do so, we review current and ongoing studies on bat evolution. We first investigate morphological specialization by reviewing current knowledge about wing and face evolution. Then, we explore the mechanisms behind adaptive diversification in various ecological contexts using vision and dentition. Finally, we highlight the emerging work into morphological novelties using bat wing membranes.

Bats and birds as viral reservoirs: A physiological and ecological perspective

The birds (class Aves) and bats (order Chiroptera, class Mammalia) are well known natural reservoirs of a diverse range of viruses, including some zoonoses. The only extant volant vertebrates, bats and birds have undergone dramatic adaptive radiations that have allowed them to occupy diverse ecological niches and colonize most of the planet. However, few studies have compared the physiology and ecology of these ecologically, and medically, important taxa. Here, we review convergent traits in the physiology, immunology, flight-related ecology of birds and bats that might enable these taxa to act as viral reservoirs and asymptomatic carriers. Many species of birds and bats are well adapted to urban environments and may host more zoonotic pathogens than species that do not colonize anthropogenic habitats. These convergent traits in birds and bats and their ecological interactions with domestic animals and humans increase the potential risk of viral spillover transmission and facilitate the emergence of novel viruses that most likely sources of zoonoses with the potential to cause global pandemics.

Bats in ecosystems and their Wide spectrum of viral infectious potential threats: SARS-CoV-2 and other emerging viruses

Bats have populated earth for approximately 52 million years, serving as natural reservoirs for a variety of viruses through the course of evolution. Transmission of highly pathogenic viruses from bats has been suspected or linked to a spectrum of potential emerging infectious diseases in humans and animals worldwide. Examples of such viruses include Marburg, Ebolavirus, Nipah, Hendra, Influenza A, Dengue, Equine Encephalitis viruses, Lyssaviruses, Madariaga and Coronaviruses, involving the now pandemic Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2). Herein, we provide a narrative review focused in selected emerging viral infectious diseases that have been reported from bats.

Identifying the causes and consequences of assembly gaps using a multiplatform genome assembly of a bird-of-paradise

Genome assemblies are currently being produced at an impressive rate by consortia and individual laboratories. The low costs and increasing efficiency of sequencing technologies now enable assembling genomes at unprecedented quality and contiguity. However, the difficulty in assembling repeat-rich and GC-rich regions (genomic "dark matter") limits insights into the evolution of genome structure and regulatory networks. Here, we compare the efficiency of currently available sequencing technologies (short/linked/long reads and proximity ligation maps) and combinations thereof in assembling genomic dark matter. By adopting different de novo assembly strategies, we compare individual draft assemblies to a curated multiplatform reference assembly and identify the genomic features that cause gaps within each assembly. We show that a multiplatform assembly implementing long-read, linked-read and proximity sequencing technologies performs best at recovering transposable elements, multicopy MHC genes, GC-rich microchromosomes and the repeat-rich W chromosome. Telomere-to-telomere assemblies are not a reality yet for most organisms, but by leveraging technology choice it is now possible to minimize genome assembly gaps for downstream analysis. We provide a roadmap to tailor sequencing projects for optimized completeness of both the coding and noncoding parts of nonmodel genomes.

A comparative genomics multitool for scientific discovery and conservation

The Zoonomia Project is investigating the genomics of shared and specialized traits in eutherian mammals. Here we provide genome assemblies for 131 species, of which all but 9 are previously uncharacterized, and describe a whole-genome alignment of 240 species of considerable phylogenetic diversity, comprising representatives from more than 80% of mammalian families. We find that regions of reduced genetic diversity are more abundant in species at a high risk of extinction, discern signals of evolutionary selection at high resolution and provide insights from individual reference genomes. By prioritizing phylogenetic diversity and making data available quickly and without restriction, the Zoonomia Project aims to support biological discovery, medical research and the conservation of biodiversity.

A whole-genome alignment of 240 phylogenetically diverse species of eutherian mammal—including 131 previously uncharacterized species—from the Zoonomia Project provides data that support biological discovery, medical research and conservation.

Wing: A suitable nonlethal tissue type for repeatable and rapid telomere length estimates in bats

Telomeres are used increasingly in ecology and evolution as biomarkers for ageing and environmental stress, and are typically measured from DNA extracted from nonlethally sampled blood. However, obtaining blood is not always possible in field conditions and only limited amounts can be taken from small mammals, such as bats, which moreover lack nucleated red blood cells and hence yield relatively low amounts of DNA. As telomere length can vary within species according to age and tissue, it is important to determine which tissues serve best as a representation of the organism as a whole. Here, we investigated whether wing tissue biopsies, a rapid and relatively noninvasive tissue collection method, could serve as a proxy for other tissues when measuring relative telomere length (rTL) in the Egyptian fruit bat (Rousettus aegyptiacus). Telomeres were measured from blood, brain, heart, kidney, liver lung, muscle and wing, and multiple wing biopsies were taken from the same individuals to determine intra-individual repeatability of rTL measured by using qPCR. Wing rTL correlated with rTL estimates from most tissues apart from blood. Blood rTL was not significantly correlated with rTL from any other tissue. Blood and muscle rTLs were significantly longer compared with other tissues, while lung displayed the shortest rTLs. Individual repeatability of rTL measures from wing tissue was high (>70%). Here we show the relationships between tissue telomere dynamics for the first time in a bat, and our results provide support for the use of wing tissue for rTL measurements.

Transcriptome sequencing of cochleae from constant-frequency and frequency-modulated echolocating bats

Echolocating bats are fascinating for their ability to ‘see’ the world in the darkness. Ultrahigh frequency hearing is essential for echolocation. In this study we collected cochlear tissues from constant-frequency (CF) bats (two subspecies of Rhinolophus affinis, Rhinolophidae) and frequency-modulated (FM) bats (Myotis ricketti, Vespertilionidae) and applied PacBio single-molecule real-time isoform sequencing (Iso-seq) technology to generate the full-length (FL) transcriptomes for the three taxa. In total of 10103, 9676 and 10504 non-redundant FL transcripts for R. a. hainanus, R. a. himalayanus and Myotis ricketti were obtained respectively. These data present a comprehensive list of transcripts involved in ultrahigh frequency hearing of echolocating bats including 26342 FL transcripts, 24833 of which are annotated by public databases. No further comparative analyses were performed on the current data in this study. This data can be reused to quantify gene or transcript expression, assess the level of alternative splicing, identify novel transcripts and improve genome annotation of bat species.

Measurement(s)	cochlea • transcriptome • sequence feature annotation
Technology Type(s)	isoform sequencing • sequence annotation
Factor Type(s)	species
Sample Characteristic - Organism	Rhinolophus affinis • Myotis ricketti

Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.12999614

Grand Challenges in Comparative Tooth Biology

Teeth are a model system for integrating developmental genomics, functional morphology, and evolution. We are at the cusp of being able to address many open issues in comparative tooth biology and we outline several of these newly tractable and exciting research directions. Like never before, technological advances and methodological approaches are allowing us to investigate the developmental machinery of vertebrates and discover both conserved and excitingly novel mechanisms of diversification. Additionally, studies of the great diversity of soft tissues, replacement teeth, and non-trophic functions of teeth are providing new insights into dental diversity. Finally, we highlight several emerging model groups of organisms that are at the forefront of increasing our appreciation of the mechanisms underlying tooth diversification.

Sequencing the Genome of Indian Flying Fox, Natural Reservoir of Nipah Virus, Using Hybrid Assembly and Conservative Secondary Scaffolding

Indian fruit bats, flying fox Pteropus medius was identified as an asymptomatic natural host of recently emerged Nipah virus, which is known to induce a severe infectious disease in humans. The absence of P. medius genome sequence presents an important obstacle for further studies of virus–host interactions and better understanding of mechanisms of zoonotic viral emergence. Generation of the high-quality genome sequence is often linked to a considerable effort associated to elevated costs. Although secondary scaffolding methods have reduced sequencing expenses, they imply the development of new tools for the integration of different data sources to achieve more reliable sequencing results. We initially sequenced the P. medius genome using the combination of Illumina paired-end and Nanopore sequencing, with a depth of 57.4x and 6.1x, respectively. Then, we introduced the novel scaff2link software to integrate multiple sources of information for secondary scaffolding, allowing to remove the association with discordant information among two sources. Different quality metrics were next produced to validate the benefits from secondary scaffolding. The P. medius genome, assembled by this method, has a length of 1,985 Mb and consists of 33,613 contigs and 16,113 scaffolds with an NG50 of 19 Mb. At least 22.5% of the assembled sequences is covered by interspersed repeats already described in other species and 19,823 coding genes are annotated. Phylogenetic analysis demonstrated the clustering of P. medius genome with two other Pteropus bat species, P. alecto and P. vampyrus, for which genome sequences are currently available. SARS-CoV entry receptor ACE2 sequence of P. medius was 82.7% identical with ACE2 of Rhinolophus sinicus bats, thought to be the natural host of SARS-CoV. Altogether, our results confirm that a lower depth of sequencing is enough to obtain a valuable genome sequence, using secondary scaffolding approaches and demonstrate the benefits of the scaff2link application. The genome sequence is now available to the scientific community to (i) proceed with further genomic analysis of P. medius, (ii) to characterize the underlying mechanism allowing Nipah virus maintenance and perpetuation in its bat host, and (iii) to monitor their evolutionary pathways toward a better understanding of bats’ ability to control viral infections.

Molecular adaptation and convergent evolution of frugivory in Old World and neotropical fruit bats

Although cases of independent adaptation to the same dietary niche have been documented in mammalian ecology, the molecular correlates of such shifts are seldom known. Here, we used genomewide analyses of molecular evolution to examine two lineages of bats that, from an insectivorous ancestor, have both independently evolved obligate frugivory: the Old World family Pteropodidae and the neotropical subfamily Stenodermatinae. New genome assemblies from two neotropical fruit bats (Artibeus jamaicensis and Sturnira hondurensis) provide a framework for comparisons with Old World fruit bats. Comparative genomics of 10 bat species encompassing dietary diversity across the phylogeny revealed convergent molecular signatures of frugivory in both multigene family evolution and single-copy genes. Evidence for convergent molecular adaptations associated with frugivorous diets includes the composition of three subfamilies of olfactory receptor genes, losses of three bitter taste receptor genes, losses of two digestive enzyme genes and convergent amino acid substitutions in several metabolic genes. By identifying suites of adaptations associated with the convergent evolution of frugivory, our analyses both reveal the extent of molecular mechanisms under selection in dietary shifts and will facilitate future studies of molecular ecology in mammals.

Six reference-quality genomes reveal evolution of bat adaptations

Bats possess extraordinary adaptations, including flight, echolocation, extreme longevity and unique immunity. High-quality genomes are crucial for understanding the molecular basis and evolution of these traits. Here we incorporated long-read sequencing and state-of-the-art scaffolding protocols1 to generate, to our knowledge, the first reference-quality genomes of six bat species (Rhinolophus ferrumequinum, Rousettus aegyptiacus, Phyllostomus discolor, Myotis myotis, Pipistrellus kuhlii and Molossus molossus). We integrated gene projections from our 'Tool to infer Orthologs from Genome Alignments' (TOGA) software with de novo and homology gene predictions as well as short- and long-read transcriptomics to generate highly complete gene annotations. To resolve the phylogenetic position of bats within Laurasiatheria, we applied several phylogenetic methods to comprehensive sets of orthologous protein-coding and noncoding regions of the genome, and identified a basal origin for bats within Scrotifera. Our genome-wide screens revealed positive selection on hearing-related genes in the ancestral branch of bats, which is indicative of laryngeal echolocation being an ancestral trait in this clade. We found selection and loss of immunity-related genes (including pro-inflammatory NF-κB regulators) and expansions of anti-viral APOBEC3 genes, which highlights molecular mechanisms that may contribute to the exceptional immunity of bats. Genomic integrations of diverse viruses provide a genomic record of historical tolerance to viral infection in bats. Finally, we found and experimentally validated bat-specific variation in microRNAs, which may regulate bat-specific gene-expression programs. Our reference-quality bat genomes provide the resources required to uncover and validate the genomic basis of adaptations of bats, and stimulate new avenues of research that are directly relevant to human health and disease1.

Near-Random Distribution of Chromosome-Derived Circular DNA in the Condensed Genome of Pigeons and the Larger, More Repeat-Rich Human Genome

Extrachromosomal circular DNA (eccDNA) elements of chromosomal origin are known to be common in a number of eukaryotic species. However, it remains to be addressed whether genomic features such as genome size, the load of repetitive elements within a genome, and/or animal physiology affect the number of eccDNAs. Here, we investigate the distribution and numbers of eccDNAs in a condensed and less repeat-rich genome compared with the human genome, using Columba livia domestica (domestic rock pigeon) as a model organism. By sequencing eccDNA in blood and breast muscle from three pigeon breeds at various ages and with different flight behavior, we characterize 30,000 unique eccDNAs. We identify genomic regions that are likely hotspots for DNA circularization in breast muscle, including genes involved in muscle development. We find that although eccDNA counts do not correlate with the biological age in pigeons, the number of unique eccDNAs in a nonflying breed (king pigeons) is significantly higher (9-fold) than homing pigeons. Furthermore, a comparison between eccDNA from skeletal muscle in pigeons and humans reveals ∼9-10 times more unique eccDNAs per human nucleus. The fraction of eccDNA sequences, derived from repetitive elements, exist in proportions to genome content, that is, human 72.4% (expected 52.5%) and pigeon 8.7% (expected 5.5%).

Overall, our results support that eccDNAs are common in pigeons, that the amount of unique eccDNA types per nucleus can differ between species as well as subspecies, and suggest that eccDNAs from repeats are found in proportions relative to the content of repetitive elements in a genome.

Bat-borne virus diversity, spillover and emergence

Most viral pathogens in humans have animal origins and arose through cross-species transmission. Over the past 50 years, several viruses, including Ebola virus, Marburg virus, Nipah virus, Hendra virus, severe acute respiratory syndrome coronavirus (SARS-CoV), Middle East respiratory coronavirus (MERS-CoV) and SARS-CoV-2, have been linked back to various bat species. Despite decades of research into bats and the pathogens they carry, the fields of bat virus ecology and molecular biology are still nascent, with many questions largely unexplored, thus hindering our ability to anticipate and prepare for the next viral outbreak. In this Review, we discuss the latest advancements and understanding of bat-borne viruses, reflecting on current knowledge gaps and outlining the potential routes for future research as well as for outbreak response and prevention efforts.

Testing the advantages and disadvantages of short- and long- read eukaryotic metagenomics using simulated reads

BACKGROUND

The first step in understanding ecological community diversity and dynamics is quantifying community membership. An increasingly common method for doing so is through metagenomics. Because of the rapidly increasing popularity of this approach, a large number of computational tools and pipelines are available for analysing metagenomic data. However, the majority of these tools have been designed and benchmarked using highly accurate short read data (i.e. Illumina), with few studies benchmarking classification accuracy for long error-prone reads (PacBio or Oxford Nanopore). In addition, few tools have been benchmarked for non-microbial communities.

RESULTS

Here we compare simulated long reads from Oxford Nanopore and Pacific Biosciences (PacBio) with high accuracy Illumina read sets to systematically investigate the effects of sequence length and taxon type on classification accuracy for metagenomic data from both microbial and non-microbial communities. We show that very generally, classification accuracy is far lower for non-microbial communities, even at low taxonomic resolution (e.g. family rather than genus). We then show that for two popular taxonomic classifiers, long reads can significantly increase classification accuracy, and this is most pronounced for non-microbial communities.

CONCLUSIONS

This work provides insight on the expected accuracy for metagenomic analyses for different taxonomic groups, and establishes the point at which read length becomes more important than error rate for assigning the correct taxon.

The Potential Role of Endogenous Viral Elements in the Evolution of Bats as Reservoirs for Zoonotic Viruses

Despite a small genome size, bats have comparable diversity of retroviral and non-retroviral endogenous sequences to other mammals. These include Class I and Class II retroviral sequences, foamy viruses, and deltaretroviruses, as well as filovirus, bornavirus, and parvovirus endogenous viral elements. Some of these endogenous viruses are sufficiently preserved in bat genomes to be expressed, with potential effects for host biology. It is clear that the bat immune system differs when compared with other mammals, yet the role that virus-derived endogenous elements may have played in the evolution of bat immunity is poorly understood. In this review, we discuss some of the bat-specific immune mechanisms that may have resulted in a virus-tolerant phenotype and link these to the long-standing virus-host coevolution that may have allowed a large diversity of endogenous retroviruses and other endogenous viral elements to colonize bat genomes. We also consider the possible effects of endogenization in the evolution of the bat immune system.

Ensembl 2019

The Ensembl project (https://www.ensembl.org) makes key genomic data sets available to the entire scientific community without restrictions. Ensembl seeks to be a fundamental resource driving scientific progress by creating, maintaining and updating reference genome annotation and comparative genomics resources. This year we describe our new and expanded gene, variant and comparative annotation capabilities, which led to a 50% increase in the number of vertebrate genomes we support. We have also doubled the number of available human variants and added regulatory regions for many mouse cell types and developmental stages. Our data sets and tools are available via the Ensembl website as well as a through a RESTful webservice, Perl application programming interface and as data files for download.

Proteomics of Long-Lived Mammals

Mammalian species differ up to 100-fold in their aging rates and maximum lifespans. Long-lived mammals appear to possess traits that extend lifespan and healthspan. Genomic analyses have not revealed a single pro-longevity function that would account for all longevity effects. In contrast, it appears that pro-longevity mechanisms may be complex traits afforded by connections between metabolism and protein functions that are impossible to predict by genomic approaches alone. Thus, metabolomics and proteomics studies will be required to understand the mechanisms of longevity. Several examples are reviewed that demonstrate the naked mole rat (NMR) shows unique proteomic signatures that contribute to longevity by overcoming several hallmarks of aging. SIRT6 is also discussed as an example of a protein that evolves enhanced enzymatic function in long-lived species. Finally, it is shown that several longevity-related proteins such as Cip1/p21, FOXO3, TOP2A, AKT1, RICTOR, INSR, and SIRT6 harbor posttranslational modification (PTM) sites that preferentially appear in either short- or long-lived species and provide examples of crosstalk between PTM sites. Prospects of enhancing lifespan and healthspan of humans by altering metabolism and proteoforms with drugs that mimic changes observed in long-lived species are discussed.