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

The genome sequence of the endemic Mexican common mustached Bat, Pteronotus mexicanus. Miller, 1902 [Mormoopidae; Pteronotus]

We describe here the first characterization of the genome of the bat Pteronotus mexicanus, an endemic species of Mexico, as part of the Mexican Bat Genome Project which focuses on the characterization and assembly of the genomes of endemic bats in Mexico. The genome was assembled from a liver tissue sample of an adult male from Jalisco, Mexico provided by the Texas Tech University Museum tissue collection. The assembled genome size was 1.9 Gb. The assembly of the genome was fitted in a framework of 110,533 scaffolds and 1,659,535 contigs. The ecological importance of bats such as P. mexicanus, and their diverse ecological roles, underscores the value of having complete genomes in addressing information gaps and facing challenges regarding their function in ecosystems and their conservation.

Tracing histoplasmosis genomic epidemiology and species occurrence across the USA

ABSTRACTHistoplasmosis is an endemic mycosis in North America frequently reported along the Ohio and Mississippi River Valleys, although autochthonous cases occur in non-endemic areas. In the United States, the disease is provoked by two genetically distinct clades of Histoplasma capsulatum sensu lato, Histoplasma mississippiense (Nam1) and H. ohiense (Nam2). To bridge the molecular epidemiological gap, we genotyped 93 Histoplasma isolates (62 novel genomes) including clinical, environmental, and veterinarian samples from a broader geographical range by whole-genome sequencing, followed by evolutionary and species niche modelling analyses. We show that histoplasmosis is caused by two major lineages, H. ohiense and H. mississippiense; with sporadic cases caused by H. suramericanum in California and Texas. While H. ohiense is prevalent in eastern states, H. mississipiense was found to be prevalent in the central and western portions of the United States, but also geographically overlapping in some areas suggesting that these species might co-occur. Species Niche Modelling revealed that H. ohiense thrives in places with warmer and drier conditions, while H. mississippiense is endemic to areas with cooler temperatures and more precipitation. In addition, we predicted multiple areas of secondary contact zones where the two species co-occur, potentially facilitating gene exchange and hybridization. This study provides the most comprehensive understanding of the genomic epidemiology of histoplasmosis in the USA and lays a blueprint for the study of invasive fungal diseases.

A molecular phylogeny for all 21 families within Chiroptera (bats)

Bats, members of the Chiroptera order, rank as the second most diverse group among mammals. Recent molecular systematic studies on bats have successfully classified 21 families within two suborders: Yinpterochiroptera and Yangochiroptera. Nevertheless, the phylogeny within these 21 families has remained a subject of controversy. In this study, we have employed a balanced approach to establish a robust family-level phylogenetic hypothesis for bats, utilizing a more comprehensive molecular dataset. This dataset includes representative species from all 21 bat families, resulting in a reduced level of missing genetic information. The resulting phylogenetic tree comprises 21 lineages that are strongly supported, each corresponding to one of the bat families. Our findings support to place the Emballonuroidea superfamily as the basal lineage of Yangochiroptera, and that Myzopodidae should be situated as a basal lineage of Emballonuroidea, forming a sister relationship with the clade consisting of Nycteridae and Emballonuridae. Finally, we have conducted dating analyses on this newly resolved phylogenetic tree, providing divergence times for each bat family. Collectively, our study has employed a relatively comprehensive molecular dataset to establish a more robust phylogeny encompassing all 21 bat families. This improved phylogenetic framework will significantly contribute to our understanding of evolutionary processes, ecological roles, disease dynamics, and biodiversity conservation in the realm of bats.

Plethora of New Marsupial Genomes Informs Our Knowledge of Marsupial MHC Class II

The major histocompatibility complex (MHC) plays a vital role in the vertebrate immune system due to its role in infection, disease and autoimmunity, or recognition of "self". The marsupial MHC class II genes show divergence from eutherian MHC class II genes and are a unique taxon of therian mammals that give birth to altricial and immunologically naive young providing an opportune study system for investigating evolution of the immune system. Additionally, the MHC in marsupials has been implicated in disease associations, including susceptibility to Chlamydia pecorum infection in koalas. Due to the complexity of the gene family, automated annotation is not possible so here we manually annotate 384 class II MHC genes in 29 marsupial species. We find losses of key components of the marsupial MHC repertoire in the Dasyuromorphia order and the Pseudochiridae family. We perform PGLS analysis to show the gene losses we find are true gene losses and not artifacts of unresolved genome assembly. We investigate the associations between the number of loci and life history traits, including lifespan and reproductive output in lineages of marsupials and hypothesize that gene loss may be linked to the energetic cost and tradeoffs associated with pregnancy and reproduction. We found support for litter size being a significant predictor of the number of DBA and DBB loci, indicating a tradeoff between the energetic requirements of immunity and reproduction. Additionally, we highlight the increased susceptibility of Dasyuridae species to neoplasia and a potential link to MHC gene loss. Finally, these annotations provide a valuable resource to the immunogenetics research community to move forward and further investigate diversity in MHC genes in marsupials.

A haplotype-resolved reference genome of a long-distance migratory bat, Pipistrellus nathusii (Keyserling & Blasius, 1839)

We present a complete, chromosome-scale reference genome for the long-distance migratory bat Pipistrellus nathusii. The genome encompasses both haplotypic sets of autosomes and the separation of both sex chromosomes by utilizing highly accurate long-reads and preserving long-range phasing information through the use of three-dimensional chromatin conformation capture sequencing (Hi-C). This genome, accompanied by a comprehensive protein-coding sequence annotation, provides a valuable genomic resource for future investigations into the genomic bases of long-distance migratory flight in bats as well as uncovering the genetic architecture, population structure and evolutionary history of Pipistrellus nathusii. The reference-quality genome presented here gives a fundamental resource to further our understanding of bat genetics and evolution, adding to the growing number of high-quality genetic resources in this field. Here, we demonstrate its use in the phylogenetic reconstruction of the order Chiroptera, and in particular, we present the resources to allow detailed investigations into the genetic drivers and adaptations related to long-distance migration.

Unraveling genomic features and phylogenomics through the analysis of three Mexican endemic Myotis genomes

Background

Genomic resource development for non-model organisms is rapidly progressing, seeking to uncover molecular mechanisms and evolutionary adaptations enabling thriving in diverse environments. Limited genomic data for bat species hinder insights into their evolutionary processes, particularly within the diverse Myotis genus of the Vespertilionidae family. In Mexico, 15 Myotis species exist, with three-M. vivesi, M. findleyi, and M. planiceps-being endemic and of conservation concern.

Methods

We obtained samples of Myotis vivesi, M. findleyi, and M. planiceps for genomic analysis. Each of three genomic DNA was extracted, sequenced, and assembled. The scaffolding was carried out utilizing the M. yumanensis genome via a genome-referenced approach within the ntJoin program. GapCloser was employed to fill gaps. Repeat elements were characterized, and gene prediction was done via ab initio and homology methods with MAKER pipeline. Functional annotation involved InterproScan, BLASTp, and KEGG. Non-coding RNAs were annotated with INFERNAL, and tRNAscan-SE. Orthologous genes were clustered using Orthofinder, and a phylogenomic tree was reconstructed using IQ-TREE.

Results

We present genome assemblies of these endemic species using Illumina NovaSeq 6000, each exceeding 2.0 Gb, with over 90% representing single-copy genes according to BUSCO analyses. Transposable elements, including LINEs and SINEs, constitute over 30% of each genome. Helitrons, consistent with Vespertilionids, were identified. Values around 20,000 genes from each of the three assemblies were derived from gene annotation and their correlation with specific functions. Comparative analysis of orthologs among eight Myotis species revealed 20,820 groups, with 4,789 being single copy orthogroups. Non-coding RNA elements were annotated. Phylogenomic tree analysis supported evolutionary chiropterans' relationships. These resources contribute significantly to understanding gene evolution, diversification patterns, and aiding conservation efforts for these endangered bat species.

Complexin-1 enhances ultrasound neurotransmission in the mammalian auditory pathway

Unlike megabats, which rely on well-developed vision, microbats use ultrasonic echolocation to navigate and locate prey. To study ultrasound perception, here we compared the auditory cortices of microbats and megabats by constructing reference genomes and single-nucleus atlases for four species. We found that parvalbumin (PV)+ neurons exhibited evident cross-species differences and could respond to ultrasound signals, whereas their silencing severely affected ultrasound perception in the mouse auditory cortex. Moreover, megabat PV+ neurons expressed low levels of complexins (CPLX1-CPLX4), which can facilitate neurotransmitter release, while microbat PV+ neurons highly expressed CPLX1, which improves neurotransmission efficiency. Further perturbation of Cplx1 in PV+ neurons impaired ultrasound perception in the mouse auditory cortex. In addition, CPLX1 functioned in other parts of the auditory pathway in microbats but not megabats and exhibited convergent evolution between echolocating microbats and whales. Altogether, we conclude that CPLX1 expression throughout the entire auditory pathway can enhance mammalian ultrasound neurotransmission.

Different species of Chiroptera: Immune cells and molecules

The distinct composition and immune response characteristics of bats' innate and adaptive immune systems, which enable them to serve as host of numerous serious zoonotic viruses without falling ill, differ substantially from those of other mammals, it have garnered significant attention. In this article, we offer a systematic review of the names, attributes, and functions of innate and adaptive immune cells & molecules across different bat species. This includes descriptions of the differences shown by research between 71 bat species in 10 families, as well as comparisons between bats and other mammals. Studies of the immune cells & molecules of different bat species are necessary to understand the unique antiviral immunity of bats. By providing comprehensive information on these unique immune responses, it is hoped that new insights will be provided for the study of co-evolutionary dynamics between viruses and the bat immune system, as well as human antiviral immunity.

Bat Bites and Rabies PEP in the Croatian Reference Centre for Rabies 1995-2020

Seroprevalence of lyssaviruses in certain bat species has been proven in the Republic of Croatia, but there have been no confirmed positive bat brain isolates or human fatalities associated with bat injuries/bites. The study included a retrospective analysis of bat injuries/bites, post-exposure prophylaxis (PEP) and geographic distribution of bat injuries in persons examined at the Zagreb Antirabies Clinic, the Croatian Reference Centre for Rabies. In the period 1995-2020, we examined a total of 21,910 patients due to animal injuries, of which 71 cases were bat-related (0.32%). Of the above number of patients, 4574 received rabies PEP (20.87%). However, for bat injuries, the proportion of patients receiving PEP was significantly higher: 66 out of 71 patients (92.95%). Of these, 33 received only the rabies vaccine, while the other 33 patients received the vaccine with human rabies immunoglobulin (HRIG). In five cases, PEP was not administered, as there was no indication for treatment. Thirty-five of the injured patients were biologists or biology students (49.29%). The bat species was confirmed in only one of the exposure cases. This was a serotine bat (Eptesicus serotinus), a known carrier of Lyssavirus hamburg. The results showed that the bat bites were rather sporadic compared to other human injuries caused by animal bites. All bat injuries should be treated as if they were caused by a rabid animal, and according to WHO recommendations. People who come into contact with bats should be strongly advised to be vaccinated against rabies. Entering bat habitats should be done with caution and in accordance with current recommendations, and nationwide surveillance should be carried out by competent institutions and in close collaboration between bat experts, epidemiologists and rabies experts.

Adaptive Evolution of the Greater Horseshoe Bat AANAT: Insights into the Link between AANAT and Hibernation Rhythms

Arylalkylamine N-acetyltransferase (AANAT) is a crucial rate-limiting enzyme in the synthesis of melatonin. AANAT has been confirmed to be independently duplicated and inactivated in different animal taxa in order to adapt to the environment. However, the evolutionary forces associated with having a single copy of AANAT remain unclear. The greater horseshoe bat has a single copy of AANAT but exhibits different hibernation rhythms in various populations. We analyzed the adaptive evolution at the gene and protein levels of AANAT from three distinct genetic lineages in China: northeast (NE), central east (CE), and southwest (SW). The results revealed greater genetic diversity in the AANAT loci of the NE and CE lineage populations that have longer hibernation times, and there were two positive selection loci. The catalytic capacity of AANAT in the Liaoning population that underwent positive selection was significantly higher than that of the Yunnan population (p < 0.05). This difference may be related to the lower proportion of α helix and the variation in two interface residues. The adaptive evolution of AANAT was significantly correlated with climate and environment (p < 0.05). After controlling for geographical factors (latitude and altitude), the evolution of AANAT by the negative temperature factor was represented by the monthly mean temperature (r = -0.6, p < 0.05). The results identified the gene level variation, functional adaptation, and evolutionary driving factors of AANAT, provide an important foundation for further understanding the adaptive evolution of the single copy of AANAT in pteropods, and may offer evidence for adaptive hibernation rhythms in bats.

Translating genomic advances into biodiversity conservation

A key action of the new Global Biodiversity Framework is the maintenance of genetic diversity in all species to safeguard their adaptive potential. To achieve this goal, a translational mindset, which aims to convert results of basic research into direct practical benefits, needs to be applied to biodiversity conservation. Despite much discussion on the value of genomics to conservation, a disconnect between those generating genomic resources and those applying it to biodiversity management remains. As global efforts to generate reference genomes for non-model species increase, investment into practical biodiversity applications is critically important. Applications such as understanding population and multispecies diversity and longitudinal monitoring need support alongside education for policymakers on integrating the data into evidence-based decisions. Without such investment, the opportunity to revolutionize global biodiversity conservation using genomics will not be fully realized.

Landscape of IGH germline genes of Chiroptera and the pattern of Rhinolophus affinis bat IGH CDR3 repertoire

The emergence and re-emergence of abundant viruses from bats that impact human and animal health have resulted in a resurgence of interest in bat immunology. Characterizing the immune receptor repertoire is critical to understanding how bats coexist with viruses in the absence of disease and developing new therapeutics to target viruses in humans and susceptible livestock. In this study, IGH germline genes of Chiroptera including Rhinolophus ferrumequinum, Phyllostomus discolor, and Pipistrellus pipistrellus were annotated, and we profiled the characteristics of Rhinolophus affinis (RA) IGH CDR3 repertoire. The germline genes of Chiroptera are quite different from those of human, mouse, cow, and dog in evolution, but the three bat species have high homology. The CDR3 repertoire of RA is unique in many aspects including CDR3 subclass, V/J genes access and pairing, CDR3 clones, and somatic high-frequency mutation compared with that of human and mouse, which is an important point in understanding the asymptomatic nature of viral infection in bats. This study unveiled a detailed map of bat IGH germline genes on chromosome level and provided the first immune receptor repertoire of bat, which will stimulate new avenues of research that are directly relevant to human health and disease.IMPORTANCEThe intricate relationship between bats and viruses has been a subject of study since the mid-20th century, with more than 100 viruses identified, including those affecting humans. While preliminary investigations have outlined the innate immune responses of bats, the role of adaptive immunity remains unclear. This study presents a pioneering contribution to bat immunology by unveiling, for the first time, a detailed map of bat IGH germline genes at the chromosome level. This breakthrough not only provides a foundation for B cell receptor research in bats but also contributes to primer design and sequencing of the CDR3 repertoire. Additionally, we offer the first comprehensive immune receptor repertoire of bats, serving as a crucial library for future comparative analyses. In summary, this research significantly advances the understanding of bats' immune responses, providing essential resources for further investigations into viral tolerance and potential zoonotic threats.

Remarkably High Repeat Content in the Genomes of Sparrows: The Importance of Genome Assembly Completeness for Transposable Element Discovery

Transposable elements (TE) play critical roles in shaping genome evolution. Highly repetitive TE sequences are also a major source of assembly gaps making it difficult to fully understand the impact of these elements on host genomes. The increased capacity of long-read sequencing technologies to span highly repetitive regions promises to provide new insights into patterns of TE activity across diverse taxa. Here we report the generation of highly contiguous reference genomes using PacBio long-read and Omni-C technologies for three species of Passerellidae sparrow. We compared these assemblies to three chromosome-level sparrow assemblies and nine other sparrow assemblies generated using a variety of short- and long-read technologies. All long-read based assemblies were longer (range: 1.12 to 1.41 Gb) than short-read assemblies (0.91 to 1.08 Gb) and assembly length was strongly correlated with the amount of repeat content. Repeat content for Bell's sparrow (31.2% of genome) was the highest level ever reported within the order Passeriformes, which comprises over half of avian diversity. The highest levels of repeat content (79.2% to 93.7%) were found on the W chromosome relative to other regions of the genome. Finally, we show that proliferation of different TE classes varied even among species with similar levels of repeat content. These patterns support a dynamic model of TE expansion and contraction even in a clade where TEs were once thought to be fairly depauperate and static. Our work highlights how the resolution of difficult-to-assemble regions of the genome with new sequencing technologies promises to transform our understanding of avian genome evolution.

The impact of induced pluripotent stem cells in animal conservation

It is widely acknowledged that we are currently facing a critical tipping point with regards to global extinction, with human activities driving us perilously close to the brink of a devastating sixth mass extinction. As a promising option for safeguarding endangered species, induced pluripotent stem cells (iPSCs) hold great potential to aid in the preservation of threatened animal populations. For endangered species, such as the northern white rhinoceros (Ceratotherium simum cottoni), supply of embryos is often limited. After the death of the last male in 2019, only two females remained in the world. IPSC technology offers novel approaches and techniques for obtaining pluripotent stem cells (PSCs) from rare and endangered animal species. Successful generation of iPSCs circumvents several bottlenecks that impede the development of PSCs, including the challenges associated with establishing embryonic stem cells, limited embryo sources and immune rejection following embryo transfer. To provide more opportunities and room for growth in our work on animal welfare, in this paper we will focus on the progress made with iPSC lines derived from endangered and extinct species, exploring their potential applications and limitations in animal welfare research.

The translational potential of studying bat immunity

Molecular studies in bats have led to the discovery of antiviral adaptations that may explain how some bat species have evolved enhanced immune tolerance towards viruses. Accumulating data suggest that some bat species have also evolved remarkable features of longevity and low rates of cancer. Furthermore, recent research strongly suggests that discovering immune adaptations in bat models can be translated to develop immune modulators and recognize alternate therapeutic strategies for diseases affecting humans. We posit that research in bat immunology will lead to discoveries that can potentially be translated to improve health outcomes in humans.

Immunological features of bats: resistance and tolerance to emerging viruses

Bats are among the most diverse mammalian species, representing over 20% of mammalian diversity. The past two decades have witnessed a disproportionate spillover of viruses from bats to humans compared with other mammalian hosts, attributed to the viral richness within bats, their phylogenetic likeness to humans, and increased human contact with wildlife. Unique evolutionary adaptations in bat genomes, particularly in antiviral protection and immune tolerance genes, enable bats to serve as reservoirs for pandemic-inducing viruses. Here, we discuss current limitations and advances made in understanding the role of bats as drivers of pandemic zoonoses. We also discuss novel technologies that have revealed spatial, dynamic, and physiological factors driving virus and host coevolution.

Experimental evidence for cancer resistance in a bat species

Mammals exhibit different rates of cancer, with long-lived species generally showing greater resistance. Although bats have been suggested to be resistant to cancer due to their longevity, this has yet to be systematically examined. Here, we investigate cancer resistance across seven bat species by activating oncogenic genes in their primary cells. Both in vitro and in vivo experiments suggest that Myotis pilosus (MPI) is particularly resistant to cancer. The transcriptomic and functional analyses reveal that the downregulation of three genes (HIF1A, COPS5, and RPS3) largely contributes to cancer resistance in MPI. Further, we identify the loss of a potential enhancer containing the HIF1A binding site upstream of COPS5 in MPI, resulting in the downregulation of COPS5. These findings not only provide direct experimental evidence for cancer resistance in a bat species but also offer insights into the natural mechanisms of cancer resistance in mammals.

Evolution of guanylate binding protein genes shows a remarkable variability within bats (Chiroptera)

Background

GBPs (guanylate binding proteins), an evolutionary ancient protein family, play a key role in the host's innate immune response against bacterial, parasitic and viral infections. In Humans, seven GBP genes have been described (GBP1-7). Despite the interest these proteins have received over the last years, evolutionary studies have only been performed in primates, Tupaia and rodents. These have shown a pattern of gene gain and loss in each family, indicative of the birth-and-death evolution process.

Results

In this study, we analysed the evolution of this gene cluster in several bat species, belonging to the Yangochiroptera and Yinpterochiroptera sub-orders. Detailed analysis shows a conserved synteny and a gene expansion and loss history. Phylogenetic analysis showed that bats have GBPs 1,2 and 4-6. GBP2 has been lost in several bat families, being present only in Hipposideidae and Pteropodidae. GBPs1, 4 and 5 are present mostly as single-copy genes in all families but have suffered duplication events, particularly in Myotis myotis and Eptesicus fuscus. Most interestingly, we demonstrate that GBP6 duplicated in a Chiroptera ancestor species originating two genes, which we named GBP6a and GBP6b, with different subsequent evolutionary histories. GBP6a underwent several duplication events in all families while GBP6b is present as a single copy gene and has been lost in Pteropodidae, Miniopteridae and Desmodus rotundus, a Phyllostomidae. With 14 and 15 GBP genes, Myotis myotis and Eptesicus fuscus stand out as having far more copies than all other studied bat species. Antagonistically, Pteropodidae have the lowest number of GBP genes in bats.

Conclusion

Bats are important reservoirs of viruses, many of which have become zoonotic diseases in the last decades. Further functional studies on bats GBPs will help elucidate their function, evolutionary history, and the role of bats as virus reservoirs.

Sensommatic: an efficient pipeline to mine and predict sensory receptor genes in the era of reference-quality genomes

SUMMARY

Sensory receptor gene families have undergone extensive expansion and loss across vertebrate evolution, leading to significant variation in receptor counts between species. However, due to their species-specific nature, conventional reference-based annotation tools often underestimate the true number of sensory receptors in a given species. While there has been an exponential increase in the taxonomic diversity of publicly available genome assemblies in recent years, only ∼30% of vertebrate species on the NCBI database are currently annotated. To overcome these limitations, we developed 'Sensommatic', an automated and accessible sensory receptor annotation pipeline. Sensommatic implements BLAST and AUGUSTUS to mine and predict sensory receptor genes from whole genome assemblies, adopting a one-to-many gene mapping approach. While designed for vertebrates, Sensommatic can be extended to run on non-vertebrate species by generating customized reference files, making it a scalable and generalizable tool.

AVAILABILITY AND IMPLEMENTATION

Source code and associated files are available at: https://github.com/GMHughes/Sensommatic.

RERconverge Expansion: Using Relative Evolutionary Rates to Study Complex Categorical Trait Evolution

Comparative genomics approaches seek to associate evolutionary genetic changes with the evolution of phenotypes across a phylogeny. Many of these methods, including our evolutionary rates based method, RERconverge, lack the capability of analyzing non-ordinal, multicategorical traits. To address this limitation, we introduce an expansion to RERconverge that associates shifts in evolutionary rates with the convergent evolution of multi-categorical traits. The categorical RERconverge expansion includes methods for performing categorical ancestral state reconstruction, statistical tests for associating relative evolutionary rates with categorical variables, and a new method for performing phylogenetic permulations on multi-categorical traits. In addition to demonstrating our new method on a three-category diet phenotype, we compare its performance to naive pairwise binary RERconverge analyses and two existing methods for comparative genomic analyses of categorical traits: phylogenetic simulations and a phylogenetic signal based method. We also present a diagnostic analysis of the new permulations approach demonstrating how the method scales with the number of species and the number of categories included in the analysis. Our results show that our new categorical method outperforms phylogenetic simulations at identifying genes and enriched pathways significantly associated with the diet phenotype and that the new ancestral reconstruction drives an improvement in our ability to capture diet-related enriched pathways. Our categorical permulations were able to account for non-uniform null distributions and correct for non-independence in gene rank during pathway enrichment analysis. The categorical expansion to RERconverge will provide a strong foundation for applying the comparative method to categorical traits on larger data sets with more species and more complex trait evolution.

Bat-associated microbes: Opportunities and perils, an overview

The potential biotechnological uses of bat-associated bacteria are discussed briefly, indicating avenues for biotechnological applications of bat-associated microbes. The uniqueness of bats in terms of their lifestyle, genomes and molecular immunology may predispose bats to act as disease reservoirs. Molecular phylogenetic analysis has shown several instances of bats harbouring the ancestral lineages of bacterial (Bartonella), protozoal (Plasmodium, Trypanosoma cruzi) and viral (SARS-CoV2) pathogens infecting humans. Along with the transmission of viruses from bats, we also discuss the potential roles of bat-associated bacteria, fungi, and protozoan parasites in emerging diseases. Current evidence suggests that environmental changes and interactions between wildlife, livestock, and humans contribute to the spill-over of infectious agents from bats to other hosts. Domestic animals including livestock may act as intermediate amplifying hosts for bat-origin pathogens to transmit to humans. An increasing number of studies investigating bat pathogen diversity and infection dynamics have been published. However, whether or how these infectious agents are transmitted both within bat populations and to other hosts, including humans, often remains unknown. Metagenomic approaches are uncovering the dynamics and distribution of potential pathogens in bat microbiomes, which might improve the understanding of disease emergence and transmission. Here, we summarize the current knowledge on bat zoonoses of public health concern and flag the gaps in the knowledge to enable further research and allocation of resources for tackling future outbreaks.

An updated synthesis of and outstanding questions in the olfactory and vomeronasal systems in bats: Genetics asks questions only anatomy can answer

The extensive diversity observed in bat nasal chemosensory systems has been well-documented at the histological level. Understanding how this diversity evolved and developing hypotheses as to why particular patterns exist require a phylogenetic perspective, which was first outlined in the work of anatomist Kunwar Bhatnagar. With the onset of genetics and genomics, it might be assumed that the puzzling patterns observed in the morphological data have been clarified. However, there is still a widespread mismatch of genetic and morphological correlations among bat chemosensory systems. Novel genomic evidence has set up new avenues to explore that demand more evidence from anatomical structures. Here, we outline the progress that has been made in both morphological and molecular studies on the olfactory and vomeronasal systems in bats since the work of Bhatnagar. Genomic data of olfactory and vomeronasal receptors demonstrate the strong need for further morphological sampling, with a particular focus on receiving brain regions, glands, and ducts.

Functional comparisons of the virus sensor RIG-I from humans, the microbat Myotis daubentonii, and the megabat Rousettus aegyptiacus, and their response to SARS-CoV-2 infection

IMPORTANCE

A common hypothesis holds that bats (order Chiroptera) are outstanding reservoirs for zoonotic viruses because of a special antiviral interferon (IFN) system. However, functional studies about key components of the bat IFN system are rare. RIG-I is a cellular sensor for viral RNA signatures that activates the antiviral signaling chain to induce IFN. We cloned and functionally characterized RIG-I genes from two species of the suborders Yangochiroptera and Yinpterochiroptera. The bat RIG-Is were conserved in their sequence and domain organization, and similar to human RIG-I in (i) mediating virus- and IFN-activated gene expression, (ii) antiviral signaling, (iii) temperature dependence, and (iv) recognition of RNA ligands. Moreover, RIG-I of Rousettus aegyptiacus (suborder Yinpterochiroptera) and of humans were found to recognize SARS-CoV-2 infection. Thus, members of both bat suborders encode RIG-Is that are comparable to their human counterpart. The ability of bats to harbor zoonotic viruses therefore seems due to other features.

Long-Read Sequencing Reveals Rapid Evolution of Immunity- and Cancer-Related Genes in Bats

Bats are exceptional among mammals for their powered flight, extended lifespans, and robust immune systems and therefore have been of particular interest in comparative genomics. Using the Oxford Nanopore Technologies long-read platform, we sequenced the genomes of two bat species with key phylogenetic positions, the Jamaican fruit bat (Artibeus jamaicensis) and the Mesoamerican mustached bat (Pteronotus mesoamericanus), and carried out a comprehensive comparative genomic analysis with a diverse collection of bats and other mammals. The high-quality, long-read genome assemblies revealed a contraction of interferon (IFN)-α at the immunity-related type I IFN locus in bats, resulting in a shift in relative IFN-ω and IFN-α copy numbers. Contradicting previous hypotheses of constitutive expression of IFN-α being a feature of the bat immune system, three bat species lost all IFN-α genes. This shift to IFN-ω could contribute to the increased viral tolerance that has made bats a common reservoir for viruses that can be transmitted to humans. Antiviral genes stimulated by type I IFNs also showed evidence of rapid evolution, including a lineage-specific duplication of IFN-induced transmembrane genes and positive selection in IFIT2. In addition, 33 tumor suppressors and 6 DNA-repair genes showed signs of positive selection, perhaps contributing to increased longevity and reduced cancer rates in bats. The robust immune systems of bats rely on both bat-wide and lineage-specific evolution in the immune gene repertoire, suggesting diverse immune strategies. Our study provides new genomic resources for bats and sheds new light on the extraordinary molecular evolution in this critically important group of mammals.

Phylogeographic Aspects of Bat Lyssaviruses in Europe: A Review

During the last few decades, bat lyssaviruses have become the topic of intensive molecular and epidemiological investigations. Since ancient times, rhabdoviruses have caused fatal encephalitis in humans which has led to research into effective strategies for their eradication. Modelling of potential future cross-species virus transmissions forms a substantial component of the recent infection biology of rabies. In this article, we summarise the available data on the phylogeography of both bats and lyssaviruses in Europe and the adjacent reg ions, especially in the contact zone between the Palearctic and Ethiopian realms. Within these zones, three bat families are present with high potential for cross-species transmission and the spread of lyssaviruses in Phylogroup II to Europe (part of the western Palearctic). The lack of effective therapies for rabies viruses in Phylogroup II and the most divergent lyssaviruses generates impetus for additional phylogenetic and virological research within this geographical region.

The genome sequence of the northern bat, Eptesicus nilssonii (Keyserling & Blasius, 1839)

We present a genome assembly from an individual Eptesicus nilssonii (the northern bat; Chordata; Mammalia; Chiroptera; Vespertilionidae), derived from the placental tissue of a pregnancy that resulted a male pup. The genome sequence is 2,064.1 megabases in span. Most of the assembly is scaffolded into 26 chromosomal pseudomolecules, including the X and Y sex chromosomes. The mitochondrial genome has also been assembled and is 17.04 kilobases in length.

The genome sequence of the Soprano Pipistrelle, Pipistrellus pygmaeus (Leach, 1825)

We present a genome assembly from an individual male Pipistrellus pygmaeus (the Soprano Pipistrelle; Chordata; Mammalia; Chiroptera; Vespertilionidae). The genome sequence is 1,895.1 megabases in span. Most of the assembly is scaffolded into 23 chromosomal pseudomolecules, including the X and Y sex chromosomes. The mitochondrial genome has also been assembled and is 17.18 kilobases in length.

Comparison of antiviral responses in two bat species reveals conserved and divergent innate immune pathways

Bats host a range of disease-causing viruses without displaying clinical symptoms. The mechanisms behind this are a continuous source of interest. Here, we studied the antiviral response in the Egyptian fruit bat and Kuhl's pipistrelle, representing two subordinal clades. We profiled the antiviral response in fibroblasts using RNA sequencing and compared bat with primate and rodent responses. Both bats upregulate similar genes; however, a subset of these genes is transcriptionally divergent between them. These divergent genes also evolve rapidly in sequence, have specific promoter architectures, and are associated with programs underlying tolerance and resistance. Finally, we characterized antiviral genes that expanded in bats, with duplicates diverging in sequence and expression. Our study reveals a largely conserved antiviral program across bats and points to a set of genes that rapidly evolve through multiple mechanisms. These can contribute to bat adaptation to viral infection and provide directions to understanding the mechanisms behind it.

Population genomic analysis reveals distinct demographics and recent adaptation in the black flying fox (Pteropus alecto)

As the only mammalian group capable of powered flight, bats have many unique biological traits. Previous comparative genomic studies in bats have focused on long-term evolution. However, the micro-evolutionary processes driving recent evolution are largely under-explored. Using resequencing data from 50 black flying foxes (Pteropus alecto), one of the model species for bats, we find that black flying fox has much higher genetic diversity and lower levels of linkage disequilibrium than most of the mammalian species. Demographic inference reveals strong population fluctuations (>100 fold) coinciding with multiple historical events including the last glacial change and Toba super eruption, suggesting that the black flying fox is a very resilient species with strong recovery abilities. While long-term adaptation in the black flying fox is enriched in metabolic genes, recent adaptation in the black flying fox has a unique landscape where recently selected genes are not strongly enriched in any functional category. The demographic history and mode of adaptation suggest that black flying fox might be a well-adapted species with strong evolutionary resilience. Taken together, this study unravels a vibrant landscape of recent evolution for the black flying fox and sheds light on several unique evolutionary processes for bats comparing to other mammalian groups.

Diet and the evolution of ADH7 across seven orders of mammals

Dietary variation within and across species drives the eco-evolutionary responsiveness of genes necessary to metabolize nutrients and other components. Recent evidence from humans and other mammals suggests that sugar-rich diets of floral nectar and ripe fruit have favoured mutations in, and functional preservation of, the ADH7 gene, which encodes the ADH class 4 enzyme responsible for metabolizing ethanol. Here we interrogate a large, comparative dataset of ADH7 gene sequence variation, including that underlying the amino acid residue located at the key site (294) that regulates the affinity of ADH7 for ethanol. Our analyses span 171 mammal species, including 59 newly sequenced. We report extensive variation, especially among frugivorous and nectarivorous bats, with potential for functional impact. We also report widespread variation in the retention and probable pseudogenization of ADH7. However, we find little statistical evidence of an overarching impact of dietary behaviour on putative ADH7 function or presence of derived alleles at site 294 across mammals, which suggests that the evolution of ADH7 is shaped by complex factors. Our study reports extensive new diversity in a gene of longstanding ecological interest, offers new sources of variation to be explored in functional assays in future study, and advances our understanding of the processes of molecular evolution.

Museomics and phylogenomics of lovebirds (Psittaciformes, Psittaculidae, Agapornis) using low-coverage whole-genome sequencing

Natural history collections contain specimens that provide important insights into studies of ecology and evolution. With the advancement of high-throughput sequencing, historical DNA (hDNA) from museum specimens has become a valuable source of genomic data to study the evolutionary history of organisms. Low-coverage whole genome sequencing (WGS) has been increasingly applied to museum specimens for analyzing organelle genomes, but is still uncommon for genotyping the nuclear DNA fraction. In this study, we applied low-coverage WGS to phylogenomic analyses of parrots in the genus Agapornis by including both modern samples and historical specimens of ∼100-year-old. Agapornis are small-sized African and Malagasy parrots with diverse characters. Earlier phylogenetic studies failed to resolve the positions of some key lineages, prohibiting a robust interpretation of the biogeography and evolution of these African parrots. Here, we demonstrated the use of low-coverage WGS for generating both mitochondrial and nuclear genomic data, and evaluated data quality differences between modern and historical samples. Our resolved Agapornis phylogeny indicates the ancestor of Agapornis likely colonized Madagascar from Australasia by trans-oceanic dispersal events before dispersing to the African continent. Genome-wide SNPs also allowed us to identify the parental origins of hybrid Agapornis individuals. This study demonstrates the potential of applying low-coverage WGS to phylogenomics and population genomics analyses and illustrates how including historical museum specimens can address outstanding questions regarding the evolutionary history of contemporary lineages.

Next Generation Sequencing Revolutionizes Organismal Biology Research in Bats

The advent of next generation sequencing technologies (NGS) has greatly accelerated our understanding of critical aspects of organismal biology from non-model organisms. Bats form a particularly interesting group in this regard, as genomic data have helped unearth a vast spectrum of idiosyncrasies in bat genomes associated with bat biology, physiology, and evolution. Bats are important bioindicators and are keystone species to many eco-systems. They often live in proximity to humans and are frequently associated with emerging infectious diseases, including the COVID-19 pandemic. Nearly four dozen bat genomes have been published to date, ranging from drafts to chromosomal level assemblies. Genomic investigations in bats have also become critical towards our understanding of disease biology and host-pathogen coevolution. In addition to whole genome sequencing, low coverage genomic data like reduced representation libraries, resequencing data, etc. have contributed significantly towards our understanding of the evolution of natural populations, and their responses to climatic and anthropogenic perturbations. In this review, we discuss how genomic data have enhanced our understanding of physiological adaptations in bats (particularly related to ageing, immunity, diet, etc.), pathogen discovery, and host pathogen co-evolution. In comparison, the application of NGS towards population genomics, conservation, biodiversity assessment, and functional genomics has been appreciably slower. We reviewed the current areas of focus, identifying emerging topical research directions and providing a roadmap for future genomic studies in bats.

What do mammals have to say about the neurobiology of acoustic communication?

Auditory communication is crucial across taxa, including humans, because it enables individuals to convey information about threats, food sources, mating opportunities, and other social cues necessary for survival. Comparative approaches to auditory communication will help bridge gaps across taxa and facilitate our understanding of the neural mechanisms underlying this complex task. In this work, we briefly review the field of auditory communication processing and the classical champion animal, the songbird. In addition, we discuss other mammalian species that are advancing the field. In particular, we emphasize mice and bats, highlighting the characteristics that may inform how we think about communication processing.

Comparing genomic and epigenomic features across species using the WashU Comparative Epigenome Browser

Genome browsers have become an intuitive and critical tool to visualize and analyze genomic features and data. Conventional genome browsers display data/annotations on a single reference genome/assembly; there are also genomic alignment viewer/browsers that help users visualize alignment, mismatch, and rearrangement between syntenic regions. However, there is a growing need for a comparative epigenome browser that can display genomic and epigenomic data sets across different species and enable users to compare them between syntenic regions. Here, we present the WashU Comparative Epigenome Browser. It allows users to load functional genomic data sets/annotations mapped to different genomes and display them over syntenic regions simultaneously. The browser also displays genetic differences between the genomes from single-nucleotide variants (SNVs) to structural variants (SVs) to visualize the association between epigenomic differences and genetic differences. Instead of anchoring all data sets to the reference genome coordinates, it creates independent coordinates of different genome assemblies to faithfully present features and data mapped to different genomes. It uses a simple, intuitive genome-align track to illustrate the syntenic relationship between different species. It extends the widely used WashU Epigenome Browser infrastructure and can be expanded to support multiple species. This new browser function will greatly facilitate comparative genomic/epigenomic research, as well as support the recent growing needs to directly compare and benchmark the T2T CHM13 assembly and other human genome assemblies.

Integrating gene annotation with orthology inference at scale

Annotating coding genes and inferring orthologs are two classical challenges in genomics and evolutionary biology that have traditionally been approached separately, limiting scalability. We present TOGA (Tool to infer Orthologs from Genome Alignments), a method that integrates structural gene annotation and orthology inference. TOGA implements a different paradigm to infer orthologous loci, improves ortholog detection and annotation of conserved genes compared with state-of-the-art methods, and handles even highly fragmented assemblies. TOGA scales to hundreds of genomes, which we demonstrate by applying it to 488 placental mammal and 501 bird assemblies, creating the largest comparative gene resources so far. Additionally, TOGA detects gene losses, enables selection screens, and automatically provides a superior measure of mammalian genome quality. TOGA is a powerful and scalable method to annotate and compare genes in the genomic era.

Insights into mammalian TE diversity through the curation of 248 genome assemblies

We examined transposable element (TE) content of 248 placental mammal genome assemblies, the largest de novo TE curation effort in eukaryotes to date. We found that although mammals resemble one another in total TE content and diversity, they show substantial differences with regard to recent TE accumulation. This includes multiple recent expansion and quiescence events across the mammalian tree. Young TEs, particularly long interspersed elements, drive increases in genome size, whereas DNA transposons are associated with smaller genomes. Mammals tend to accumulate only a few types of TEs at any given time, with one TE type dominating. We also found association between dietary habit and the presence of DNA transposon invasions. These detailed annotations will serve as a benchmark for future comparative TE analyses among placental mammals.

Positive Selection and Duplication of Bat TRIM Family Proteins

Bats have received increasing attention because of some unique biological features they possess. TRIM is a large family of proteins that participate in diverse cellular functions, such as antiviral immunity, DNA damage repair, tumor suppression, and aging. These functional areas appear to be highly consistent with the special characteristics of bats, such as tolerance to viruses and DNA damage generated in flight, low cancer incidence, and longevity. However, there is still a lack of systematic study of the TRIM family in bats. Here, we explored the TRIM family of bats using the genomes of 16 representative species. The results showed that the bat TRIM family contains 70 members, with 24 under positive selection and 7 duplicated. Additional transcriptomic analysis revealed the tissue-specific expressions of TRIM9, 46, 54, 55, 63, and 72. Additionally, following interferon or viral stimulation, TRIM orthologs associated with antiviral immunity reported in humans were also upregulated in bat cells. The present study systematically analyzed the composition, evolution, and expression of bat TRIM genes. It may provide a theoretical basis for studies of bat TRIM in the fields of antiviral immunity, longevity, and tolerance to DNA damage.

Bat pluripotent stem cells reveal unusual entanglement between host and viruses

Bats are distinctive among mammals due to their ability to fly, use laryngeal echolocation, and tolerate viruses. However, there are currently no reliable cellular models for studying bat biology or their response to viral infections. Here, we created induced pluripotent stem cells (iPSCs) from two species of bats: the wild greater horseshoe bat (Rhinolophus ferrumequinum) and the greater mouse-eared bat (Myotis myotis). The iPSCs from both bat species showed similar characteristics and had a gene expression profile resembling that of cells attacked by viruses. They also had a high number of endogenous viral sequences, particularly retroviruses. These results suggest that bats have evolved mechanisms to tolerate a large load of viral sequences and may have a more intertwined relationship with viruses than previously thought. Further study of bat iPSCs and their differentiated progeny will provide insights into bat biology, virus host relationships, and the molecular basis of bats' special traits.

Oral Sampling of Little Brown Bat (Myotis lucifugus) Maternity Colonies for SARS-CoV-2 in the Northeast and Mid-Atlantic, USA

The potential introduction of SARS-CoV-2, the virus responsible for the COVID-19 pandemic, into North American bat populations is of interest to wildlife managers due to recent disease-mediated declines of several species. Populations of little brown bats (Myotis lucifugus) have collapsed due to white-nose syndrome (WNS), a disease caused by the introduction and spread of the fungal pathogen Pseudogymnoascus destructans (Pd). Throughout much of the United States and southern Canada, large colonies of the species routinely established diurnal roosts in anthropogenic structures, thereby creating the potential for direct human contact and cross-species disease transmission. Given recent declines and the potential for further disease impacts, we collected oral swabs from eight little brown bat maternity colonies to assess the presence and prevalence of SARS-CoV-2 by RT-qPCR analysis. Little brown bat colonies in Maryland (n = 1), New Hampshire (n = 1), New Jersey (n = 2), New York (n = 1), Rhode Island (n = 2), and Virginia (n = 1) were taken during May-August, 2022. From 235 assayed individuals, no bat tested positive for SARS-CoV-2. Our results indicate that little brown bats may not contract SARS-CoV-2 or that the virus persists at undetectable levels in populations of the Mid-Atlantic and Northeast during summer months. Nonetheless, continued monitoring and future work addressing other seasons may still be warranted to conclusively determine infection status.

Bat responses to climate change: a systematic review

Understanding how species respond to climate change is key to informing vulnerability assessments and designing effective conservation strategies, yet research efforts on wildlife responses to climate change fail to deliver a representative overview due to inherent biases. Bats are a species-rich, globally distributed group of organisms that are thought to be particularly sensitive to the effects of climate change because of their high surface-to-volume ratios and low reproductive rates. We systematically reviewed the literature on bat responses to climate change to provide an overview of the current state of knowledge, identify research gaps and biases and highlight future research needs. We found that studies are geographically biased towards Europe, North America and Australia, and temperate and Mediterranean biomes, thus missing a substantial proportion of bat diversity and thermal responses. Less than half of the published studies provide concrete evidence for bat responses to climate change. For over a third of studied bat species, response evidence is only based on predictive species distribution models. Consequently, the most frequently reported responses involve range shifts (57% of species) and changes in patterns of species diversity (26%). Bats showed a variety of responses, including both positive (e.g. range expansion and population increase) and negative responses (range contraction and population decrease), although responses to extreme events were always negative or neutral. Spatial responses varied in their outcome and across families, with almost all taxonomic groups featuring both range expansions and contractions, while demographic responses were strongly biased towards negative outcomes, particularly among Pteropodidae and Molossidae. The commonly used correlative modelling approaches can be applied to many species, but do not provide mechanistic insight into behavioural, physiological, phenological or genetic responses. There was a paucity of experimental studies (26%), and only a small proportion of the 396 bat species covered in the examined studies were studied using long-term and/or experimental approaches (11%), even though they are more informative about the effects of climate change. We emphasise the need for more empirical studies to unravel the multifaceted nature of bats' responses to climate change and the need for standardised study designs that will enable synthesis and meta-analysis of the literature. Finally, we stress the importance of overcoming geographic and taxonomic disparities through strengthening research capacity in the Global South to provide a more comprehensive view of terrestrial biodiversity responses to climate change.

Divergence in interspecific and intersubspecific gene expression between two closely related horseshoe bats (Rhinolophus)

Closely related species have been used as representative systems to investigate the genetic mechanisms involved in the early stages of species differentiation. Previous studies have indicated that variation in gene expression might be a sensitive indicator of initial species divergence, although the role of expression divergence, and especially that associated with phenotypic variation remained relatively undefined. For three organs (cochlea, brain, and liver) from two closely related bat species (Rhinolophus siamensis and R. episcopus), the interspecific and intersubspecific gene expression profiles were compared using transcriptomics in this study. Striking organ specificity of expression was observed, and expression profiles exhibited similarities between cochlea and brain tissues. Numerous differentially expressed genes (DEGs) were identified for each organ in the interspecific comparison (cochlea/brain/liver: 1,069/647/692) and intersubspecific comparison (608/528/368). Functional enrichment analysis indicated vital variation in expression related to the immune system, ion activities, neuronal function, and multisensory system regulation in both comparisons. DEGs relevant to the variation in echolocation calls (RF) were found, and some of them were involved in the pivotal patterns of expression variation. The regulation of immune, ion channel, neural activity, and sophisticated sensory functions at the expression level might be key mechanisms in the early species divergence of bats, and the expression variation related to acoustical signal could have played a crucial part. This study expands our knowledge of gene expression and patterns of variation for three key organs to echolocation at both the interspecific and intersubspecific levels. Further, the framework described here provides insight into the genetic basis of phenotypic variation during the incipient stage of species differentiation.

Taking flight: An ecological, evolutionary and genomic perspective on bat telomeres

Over 20% of all living mammals are bats (order Chiroptera). Bats possess extraordinary adaptations including powered flight, laryngeal echolocation and a unique immune system that enables them to tolerate a diversity of viral infections without presenting clinical disease symptoms. They occupy multiple trophic niches and environments globally. Significant physiological and ecological diversity occurs across the order. Bats also exhibit extreme longevity given their body size with many species showing few signs of ageing. The molecular basis of this extended longevity has recently attracted attention. Telomere maintenance potentially underpins bats' extended healthspan, although functional studies are still required to validate the causative mechanisms. In this review, we detail the current knowledge on bat telomeres, telomerase expression, and how these relate to ecology, longevity and life-history strategies. Patterns of telomere shortening and telomerase expression vary across species, and comparative genomic analyses suggest that alternative telomere maintenance mechanisms evolved in the longest-lived bats. We discuss the unique challenges faced when working with populations of wild bats and highlight ways to advance the field including expanding long-term monitoring across species that display contrasting life-histories and occupy different environmental niches. We further review how new high quality, chromosome-level genome assemblies can enable us to uncover the molecular mechanisms governing telomere dynamics and how phylogenomic analyses can reveal the adaptive significance of telomere maintenance and variation in bats.

Gene expression changes during the evolution of the tetrapod limb

Major changes in the vertebrate anatomy have preceded the conquest of land by the members of this taxon, and continuous changes in limb shape and use have occurred during the later radiation of tetrapods. While the main, conserved mechanisms of limb development have been discerned over the past century using a combination of classical embryological and molecular methods, only recent advances made it possible to identify and study the regulatory changes that have contributed to the evolution of the tetrapod appendage. These advances include the expansion of the model repertoire from traditional genetic model species to non-conventional ones, a proliferation of predictive mathematical models that describe gene interactions, an explosion in genomic data and the development of high-throughput methodologies. These revolutionary innovations make it possible to identify specific mutations that are behind specific transitions in limb evolution. Also, as we continue to apply them to more and more extant species, we can expect to gain a fine-grained view of this evolutionary transition that has been so consequential for our species as well.

Single-cell transcriptome analysis of the in vivo response to viral infection in the cave nectar bat Eonycteris spelaea

Bats are reservoir hosts of many zoonotic viruses with pandemic potential. We utilized single-cell transcriptome sequencing (scRNA-seq) to analyze the immune response in bat lungs upon in vivo infection with a double-stranded RNA virus, Pteropine orthoreovirus PRV3M. Bat neutrophils were distinguished by high basal IDO1 expression. NK cells and T cells were the most abundant immune cells in lung tissue. Three distinct CD8⁺ effector T cell populations could be delineated by differential expression of KLRB1, GFRA2, and DPP4. Select NK and T clusters increased expression of genes involved in T cell activation and effector function early after viral infection. Alveolar macrophages and classical monocytes drove antiviral interferon signaling. Infection expanded a CSF1R⁺ population expressing collagen-like genes, which became the predominant myeloid cell type post-infection. This work uncovers features relevant to viral disease tolerance in bats, lays a foundation for future experimental work, and serves as a resource for comparative immunology studies.

The pale spear-nosed bat: A neuromolecular and transgenic model for vocal learning

Vocal learning, the ability to produce modified vocalizations via learning from acoustic signals, is a key trait in the evolution of speech. While extensively studied in songbirds, mammalian models for vocal learning are rare. Bats present a promising study system given their gregarious natures, small size, and the ability of some species to be maintained in captive colonies. We utilize the pale spear-nosed bat (Phyllostomus discolor) and report advances in establishing this species as a tractable model for understanding vocal learning. We have taken an interdisciplinary approach, aiming to provide an integrated understanding across genomics (Part I), neurobiology (Part II), and transgenics (Part III). In Part I, we generated new, high-quality genome annotations of coding genes and noncoding microRNAs to facilitate functional and evolutionary studies. In Part II, we traced connections between auditory-related brain regions and reported neuroimaging to explore the structure of the brain and gene expression patterns to highlight brain regions. In Part III, we created the first successful transgenic bats by manipulating the expression of FoxP2, a speech-related gene. These interdisciplinary approaches are facilitating a mechanistic and evolutionary understanding of mammalian vocal learning and can also contribute to other areas of investigation that utilize P. discolor or bats as study species.

Best genome sequencing strategies for annotation of complex immune gene families in wildlife

BACKGROUND

The biodiversity crisis and increasing impact of wildlife disease on animal and human health provides impetus for studying immune genes in wildlife. Despite the recent boom in genomes for wildlife species, immune genes are poorly annotated in nonmodel species owing to their high level of polymorphism and complex genomic organisation. Our research over the past decade and a half on Tasmanian devils and koalas highlights the importance of genomics and accurate immune annotations to investigate disease in wildlife. Given this, we have increasingly been asked the minimum levels of genome quality required to effectively annotate immune genes in order to study immunogenetic diversity. Here we set out to answer this question by manually annotating immune genes in 5 marsupial genomes and 1 monotreme genome to determine the impact of sequencing data type, assembly quality, and automated annotation on accurate immune annotation.

RESULTS

Genome quality is directly linked to our ability to annotate complex immune gene families, with long reads and scaffolding technologies required to reassemble immune gene clusters and elucidate evolution, organisation, and true gene content of the immune repertoire. Draft-quality genomes generated from short reads with HiC or 10× Chromium linked reads were unable to achieve this. Despite mammalian BUSCOv5 scores of up to 94.1% amongst the 6 genomes, automated annotation pipelines incorrectly annotated up to 59% of manually annotated immune genes regardless of assembly quality or method of automated annotation.

CONCLUSIONS

Our results demonstrate that long reads and scaffolding technologies, alongside manual annotation, are required to accurately study the immune gene repertoire of wildlife species.

Genome assembly of the Pendlebury's roundleaf bat, Hipposideros pendleburyi, revealed the expansion of Tc1/Mariner DNA transposons in Rhinolophoidea

Bats (Chiroptera) constitute the second largest order of mammals and have several distinctive features, such as true self-powered flight and strong immunity. The Pendlebury's roundleaf bat, Hipposideros pendleburyi, is endemic to Thailand and listed as a vulnerable species. We employed the 10× Genomics linked-read technology to obtain a genome assembly of H. pendleburyi. The assembly size was 2.17 Gb with a scaffold N50 length of 15,398,518 bases. Our phylogenetic analysis placed H. pendleburyi within the rhinolophoid clade of the suborder Yinpterochiroptera. A synteny analysis showed that H. pendleburyi shared conserved chromosome segments (up to 105 Mb) with Rhinolophus ferrumequinum and Phyllostomus discolor albeit having different chromosome numbers and belonging different families. We found positive selection signals in genes involved in inflammation, spermatogenesis and Wnt signalling. The analyses of transposable elements suggested the contraction of short interspersed nuclear elements (SINEs) and the accumulation of young mariner DNA transposons in the analysed hipposiderids. Distinct mariners were likely horizontally transferred to hipposiderid genomes over the evolution of this family. The lineage-specific profiles of SINEs and mariners might involve in the evolution of hipposiderids and be associated with the phylogenetic separations of these bats from other bat families.

The new highest number of B chromosomes (Bs) in Leisler's bat Nyctalusleisleri (Kuhl, 1817)

B chromosomes (Bs) are supernumerary to the standard chromosome set, from which they prevalently derive. Variation in numbers both among individuals or populations and among cells within individuals is their constant feature. Leisler's bat Nyctalusleisleri (Kuhl, 1817) is one of only four species of Chiroptera with detected Bs. Four males of N.leisleri were collected from two localities on the territory of Serbia and cytogenetically analysed. All animals had Bs with interindividual variability ranging from two to five heterochromatic micro Bs. The highest number of Bs was detected in this species. Among mammals, Rodentia and Chiroptera are orders with the largest number of species, but Bs frequently appear in rodents and rarely in chiropterans. Possible explanations for this difference are offered.

Unique Evolution of Antiviral Tetherin in Bats

Bats are recognized as important reservoirs of viruses deadly to other mammals, including humans. These infections are typically nonpathogenic in bats, raising questions about host response differences that might exist between bats and other mammals. Tetherin is a restriction factor which inhibits the release of a diverse range of viruses from host cells, including retroviruses, coronaviruses, filoviruses, and paramyxoviruses, some of which are deadly to humans and transmitted by bats. Here, we characterize the tetherin genes from 27 bat species, revealing that they have evolved under strong selective pressure, and that fruit bats and vesper bats express unique structural variants of the tetherin protein. Tetherin was widely and variably expressed across fruit bat tissue types and upregulated in spleen tissue when stimulated with Toll-like receptor agonists. The expression of two computationally predicted splice isoforms of fruit bat tetherin was verified. We identified an additional third unique splice isoform which includes a C-terminal region that is not homologous to known mammalian tetherin variants but was functionally capable of restricting the release of filoviral virus-like particles. We also report that vesper bats possess and express at least five tetherin genes, including structural variants, more than any other mammal reported to date. These findings support the hypothesis of differential antiviral gene evolution in bats relative to other mammals.

IMPORTANCE Bats are an important host of various viruses which are deadly to humans and other mammals but do not cause outward signs of illness in bats. Furthering our understanding of the unique features of the immune system of bats will shed light on how they tolerate viral infections, potentially informing novel antiviral strategies in humans and other animals. This study examines the antiviral protein tetherin, which prevents viral particles from escaping their host cell. Analysis of tetherin from 27 bat species reveals that it is under strong evolutionary pressure, and we show that multiple bat species have evolved to possess more tetherin genes than other mammals, some of which encode structurally unique tetherins capable of activity against different viral particles. These data suggest that bat tetherin plays a potentially broad and important role in the management of viral infections in bats.