Deciphering bat influenza H18N11 infection dynamics in male Jamaican fruit bats on a single-cell level.

Jamaican fruit bats (Artibeus jamaicensis) naturally harbor a wide range of viruses of human relevance. These infections are typically mild in bats, suggesting unique features of their immune system. To better understand the immune response to viral infections in bats, we infected male Jamaican fruit bats with the bat-derived influenza A virus (IAV) H18N11. Using comparative single-cell RNA sequencing, we generated single-cell atlases of the Jamaican fruit bat intestine and mesentery. Gene expression profiling showed that H18N11 infection resulted in a moderate induction of interferon-stimulated genes and transcriptional activation of immune cells. H18N11 infection was predominant in various leukocytes, including macrophages, B cells, and NK/T cells. Confirming these findings, human leukocytes, particularly macrophages, were also susceptible to H18N11, highlighting the zoonotic potential of this bat-derived IAV. Our study provides insight into a natural virus-host relationship and thus serves as a fundamental resource for future in-depth characterization of bat immunology.

A chromosome-level genome assembly of an avivorous bat species (Nyctalus aviator).

Currently, three carnivorous bat species, namely Ia io, Nyctalus lasiopterus, and Nyctalus aviator, are known to actively prey on seasonal migratory birds (hereinafter referred to as "avivorous bats"). However, the absence of reference genomes impedes a thorough comprehension of the molecular adaptations of avivorous bat species. Herein, we present the high-quality chromosome-scale reference genome of N. aviator based on PacBio subreads, DNBSEQ short-reads and Hi-C sequencing data. The genome assembly size of N. aviator is 1.77 Gb, with a scaffold N50 of 102 Mb, of which 99.8% assembly was anchored into 21 pseudo-chromosomes. After masking 635.1 Mb repetitive sequences, a total of 19,412 protein-coding genes were identified, of which 99.3% were functionally annotated. The genome assembly and gene prediction reached 96.1% and 96.1% completeness of Benchmarking Universal Single-Copy Orthologs (BUSCO), respectively. This chromosome-level reference genome of N. aviator fills a gap in the existing information on the genomes of carnivorous bats, especially avivorous ones, and will be valuable for mechanism of adaptations to dietary niche expansion in bat species.

A comparative study on eDNA-based detection of Siamese bat catfish (Oreoglanis siamensis) in wet and dry conditions.

The use of environmental DNA (eDNA) analysis has demonstrated notable efficacy in detecting the existence of freshwater species, including those that are endangered or uncommon. This application holds significant potential for enhancing environmental monitoring and management efforts. However, the efficacy of eDNA-based detection relies on several factors. In this study, we assessed the impact of rainfall on the detection of eDNA for the Siamese bat catfish (Oreoglanis siamensis). Quantitative polymerase chain reaction (qPCR) analysis indicated that samples from days with average rainfall exceeding 35 mm (classified as heavy and very heavy rain) yielded negative results. While eDNA detection remains feasible on light or moderate rainy days, a noteworthy reduction in eDNA concentration and qPCR-positive likelihood was observed. Analysis across 12 sampling sites established a statistically significant negative relationship (p < 0.001) between eDNA detection and rainfall. Specifically, for each 1 mm increase in rainfall, there was an observed drop in eDNA concentration of 0.19 copies/mL (±0.14). The findings of this study provide definitive evidence that precipitation has a significant impact on the detection of eDNA in Siamese bat catfish. However, in the case of adverse weather conditions occurring on the day of sampling, our research indicates that it is acceptable to continue with the task, as long as the rainfall is not heavy or very heavy. To enhance the effectiveness of an eDNA survey, it is crucial to consider many factors related to climatic conditions. The aforementioned factor holds significant importance not only for the specific species under scrutiny but also for the broader dynamics of the climate.

Comparative mitochondrial genomics of endemic Mexican vesper yellow bats genus Rhogeessa (Chiroptera: Vespertilionidae) and insights into internal relationships in the family Vespertilionidae.

In the species-rich family Vespertilionidae, vesper yellow bats in the genus Rhogeessa include eleven species, three of them endemic to Mexico. These insectivorous bats provide important ecosystem services, including pest control. Even though some aspects of their biology are well- known, only a few genomic resources are available for these species, which limits our understanding of their biology. In this study, we assembled and annotated the mitochondrial genome of four species: R. aenea, R. genowaysi, R. mira, and R. parvula. We generated a phylomitogenomic hypothesis based on translated protein-coding genes for a total of 52 species in the family Vespertilionidae and examined the phylogenetic position of the genus Rhogeessa and species within the family. The AT-rich mitogenomes of R. aenea, R. genowaysi, R. mira, and R. parvula are 16,763, 16,781, 16,807, and 16,794 pb in length, respectively. Each studied mitogenome encodes 13 Protein Coding Genes (PCGs), 22 transfer RNA genes, and 2 rRNA genes, and contains a putative control region (CR). All tRNAs exhibit a 'cloverleaf' secondary structure, except tRNA-Serine-1 that lacked the DHU arm in all studied mitogenomes. Selective pressure analyses indicated that all protein-coding genes are exposed to purifying selection. The phylomitogenomic analysis supported the monophyletic status of the family Vespertilionidae, confirmed the placement of Rhogeessa within the tribe Antrozoini, and clarified phylogenetic relationships within and among subfamilies and tribes in this family. Our results indicate that phylomitogenomics are useful to explore the evolutionary history of vesper bats. The assembly and comprehensive analysis of mitochondrial genomes offer the potential to generate molecular references and resources beneficial for genetic analyses aimed at understanding the ecology and evolution of these remarkable bats.

Molecular adaptations underlying high-frequency hearing in the brain of CF bats species.

BACKGROUND: The majority of bat species have developed remarkable echolocation ability, especially for the laryngeally echolocating bats along with high-frequency hearing. Adaptive evolution has been widely detected for the cochleae in the laryngeally echolocating bats, however, limited understanding for the brain which is the central to echolocation signal processing in the auditory perception system, the laryngeally echolocating bats brain may also undergo adaptive changes. RESULT: In order to uncover the molecular adaptations related with high-frequency hearing in the brain of laryngeally echolocating bats, the genes expressed in the brain of Rhinolophus ferrumequinum (CF bat) and Myotis pilosus (FM bat) were both detected and also compared. A total of 346,891 genes were detected and the signal transduction mechanisms were annotated by the most abundant genes, followed by the transcription. In hence, there were 3,088 DEGs were found between the two bat brains, with 1,426 highly expressed in the brain of R. ferrumequinum, which were significantly enriched in the neuron and neurodevelopmental processes. Moreover, we found a key candidate hearing gene, ADCY1, playing an important role in the R. ferrumequinum brain and undergoing adaptive evolution in CF bats. CONCLUSIONS: Our study provides a new insight to the molecular bases of high-frequency hearing in two laryngeally echolocating bats brain and revealed different nervous system activities during auditory perception in the brain of CF bats.

The Role of (Co)variation in Shaping the Response to Selection in New World Leaf-Nosed Bats.

AbstractUnderstanding and predicting the evolutionary responses of complex morphological traits to selection remains a major challenge in evolutionary biology. Because traits are genetically correlated, selection on a particular trait produces both direct effects on the distribution of that trait and indirect effects on other traits in the population. The correlations between traits can strongly impact evolutionary responses to selection and may thus impose constraints on adaptation. Here, we used museum specimens and comparative quantitative genetic approaches to investigate whether the covariation among cranial traits facilitated or constrained the response to selection during the major dietary transitions in one of the world's most ecologically diverse mammalian families-the phyllostomid bats. We reconstructed the set of net selection gradients that would have acted on each cranial trait during the major transitions to feeding specializations and decomposed the selection responses into their direct and indirect components. We found that for all transitions, most traits capturing craniofacial length evolved toward adaptive directions owing to direct selection. Additionally, we showed instances of dietary transitions in which the complex interaction between the patterns of covariation among traits and the strength and direction of selection either constrained or facilitated evolution. Our work highlights the importance of considering the within-species covariation estimates to quantify evolvability and to disentangle the relative contribution of variational constraints versus selective causes for observed patterns.

Interferon-induced protein with tetratricopeptide repeats 5 of black fruit bat (Pteropus alecto) displays a broad inhibition of RNA viruses.

Bats are natural host reservoirs and have adapted a unique innate immune system that permits them to host many viruses without exhibiting symptoms. Notably, bat interferon stimulated genes (ISGs) have been shown to play antiviral roles. Interferon induced protein with tetratricopeptide repeats 5 (IFIT5) is a well-characterised ISG in humans with antiviral activities against negative-sense RNA viruses via inhibiting viral transcription. Here, we aim to investigate if Pteropus alecto (pa) IFIT5 (paIFIT5) possess the ability to inhibit negative-sense RNA viruses. Initially, gene syntenic and comparative structural analyses of multiple animals highlighted a high level of similarity between Pteropus alecto and human IFIT5 proteins. Our results showed that paIFIT5 was significantly inducible by viral and dsRNA stimulation. Transient overexpression of paIFIT5 inhibited the replication of vesicular stomatitis virus (VSV). Using minireplicon and transcription reporter assays, we demonstrated the ability of paIFIT5 specifically to inhibit H17N10 polymerase activity. Mechanistically, we noticed that the antiviral potential of paIFIT5 against negative sense RNA viruses was retributed to its interaction with 5'ppp containing RNA. Taken together, these findings highlight the genetic and functional conservation of IFIT5 among mammals.

The hierarchical radiation of phyllostomid bats as revealed by adaptive molar morphology.

Adaptive radiations are bursts in biodiversity that generate new evolutionary lineages and phenotypes. However, because they typically occur over millions of years, it is unclear how their macroevolutionary dynamics vary through time and among groups of organisms. Phyllostomid bats radiated extensively for diverse diets-from insects to vertebrates, fruit, nectar, and blood-and we use their molars as a model system to examine the dynamics of adaptive radiations. Three-dimensional shape analyses of lower molars of Noctilionoidea (Phyllostomidae and close relatives) indicate that different diet groups exhibit distinct morphotypes. Comparative analyses further reveal that phyllostomids are a striking example of a hierarchical radiation; phyllostomids' initial, higher-level diversification involved an "early burst" in molar morphological disparity as lineages invaded new diet-affiliated adaptive zones, followed by subsequent lower-level diversifications within adaptive zones involving less dramatic morphological changes. We posit that strong selective pressures related to initial shifts to derived diets may have freed molars from morpho-functional constraints associated with the ancestral molar morphotype. Then, lineages with derived diets (frugivores and nectarivores) diversified within broad adaptive zones, likely reflecting finer-scale niche partitioning. Importantly, the observed early burst pattern is only evident when examining molar traits that are strongly linked to diet, highlighting the value of ecomorphological traits in comparative studies. Our results support the hypothesis that adaptive radiations are commonly hierarchical and involve different tempos and modes at different phylogenetic levels, with early bursts being more common at higher levels.

Parallel evolution in an island archipelago revealed by genomic sequencing of Hipposideros leaf-nosed bats.

Body size is a key morphological attribute, often used to delimit species boundaries among closely related taxa. But body size can evolve in parallel, reaching similar final states despite independent evolutionary and geographic origins, leading to faulty assumptions of evolutionary history. Here, we document parallel evolution in body size in the widely distributed leaf-nosed bat genus Hipposideros, which has misled both taxonomic and evolutionary inference. We sequenced reduced representation genomic loci and measured external morphological characters from three closely related species from the Solomon Islands archipelago, delimited by body size. Species tree reconstruction confirms the paraphyly of two morphologically designated species. The nonsister relationship between large-bodied H. dinops lineages found on different islands indicates that large-bodied ecomorphs have evolved independently at least twice in the history of this radiation. A lack of evidence for gene flow between sympatric, closely related taxa suggests the rapid evolution of strong reproductive isolating barriers between morphologically distinct populations. Our results position Solomon Islands Hipposideros as a novel vertebrate system for studying the repeatability of parallel evolution under natural conditions. We conclude by offering testable hypotheses for how geography and ecology could be mediating the repeated evolution of large-bodied Hipposideros lineages in the Solomon Islands.

Adapting to stress: The effects of hibernation and hibernacula temperature on the hepatic transcriptome of Rhinolophus pusillus.

Hibernation, a survival strategy in mammals for extreme climates, induces physiological phenomena such as ischemia-reperfusion and metabolic shifts that hold great potential for advancements in modern medicine. Despite this, the molecular mechanisms underpinning hibernation remain largely unclear. This study used RNA-seq and Iso-seq techniques to investigate the changes in liver transcriptome expression of Rhinolophus pusillus during hibernation and active periods, as well as under different microhabitat temperatures. We identified 11 457 differentially expressed genes during hibernation and active periods, of which 395 showed significant differential expression. Genes associated with fatty acid catabolism were significantly upregulated during hibernation, whereas genes related to carbohydrate metabolism and glycogen synthesis were downregulated. Conversely, immune-related genes displayed differential expression patterns: genes tied to innate immunity were significantly upregulated, while those linked to adaptive immunity and inflammatory response were downregulated. The analysis of transcriptomic data obtained from different microhabitat temperatures revealed that R. pusillus exhibited an upregulation of genes associated with lipid metabolism in lower microhabitat temperature. This upregulation facilitated an enhanced utilization rate of triglyceride, ultimately resulting in increased energy provision for the organism. Additionally, R. pusillus upregulated gluconeogenesis-related genes regardless of the microhabitat temperature, demonstrating the importance of maintaining blood glucose levels during hibernation. Our transcriptomic data reveal that these changes in liver gene expression optimize energy allocation during hibernation, suggesting that liver tissue adaptively responds to the inherent stress of its function during hibernation. This study sheds light on the role of differential gene expression in promoting more efficient energy allocation during hibernation. It contributes to our understanding of how liver tissue adapts to the stressors associated with this state.

Phylogenomic analyses unraveled the evolution of viral tolerance in bats.

Bats host a range of viruses, exhibiting a coevolution process with many virus genera and a special capacity for viral tolerance. Foley et al.1 performed phylogenomic analyses for 60 bat species, finding that swarming behavior might facilitate cross-species introgression and the spread of anti-virus immunity gene loci across species.

Molecular and morphological data suggest a new species of big-eared bat (Vespertilionidae: Corynorhinus) endemic to northeastern Mexico.

Corynorhinus mexicanus is an insectivorous bat endemic to Mexico that inhabits the high and humid regions of the Sierra Madre Oriental (SMO), the Trans-Mexican Volcanic Belt (TMVB), and the Sierra Madre Occidental (SMOC). A previous study suggested that C. mexicanus could be a cryptic species complex due to the genetic divergence observed between specimens from the TMVB and SMOC. The present study implemented phylogenetic, population genetics, and morphological analyses to evaluate the hypothesis that C. mexicanus is a species complex. The phylogenetic analysis indicated that C. mexicanus is a polyphyletic species composed of three indirectly related lineages. The estimated divergence times for the lineages suggest that they first originated during the Pliocene, while the second and third shared a common ancestor with C. townsendii 1.55 million years ago, and diverged 600,000 years ago during the Middle Pleistocene. The population genetics analysis reveals the SMO lineage of C. mexicanus is an isolated genetic group and highly diverged from the rest of lineages (SMOC and TMVB). The morphological analyses showed variation in the skull and mandible associated with the lineages and sex of the specimens, highlighting a difference in mandible shape between the specimens of the SMO and the rest of C. mexicanus. The results of this study suggest the presence of an undescribed species of the genus Corynorhinus.

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.

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.

Resting posture drives the evolution of agonistic displays in bats.

Agonistic displays are one of the most diverse social behaviors that have important functions in animal's life history. However, their origin and driving factors have largely been unexplored. Here, we evaluated agonistic displays of 71 bat species across 10 families and classified these displays into two categories: (a) boxing displays where a bat attacks its opponent with its wrist and thumb and (b) pushing displays where a bat uses its head or body to hit a rival. We estimated the strength of the phylogenetic signal of the agonistic displays, revealed their origin, and tested the potential evolutionary relationships between agonistic behaviors and body size or resting posture (free hanging vs. contact hanging where the bat is in contact with some surface). We found that agonistic displays were phylogenetically conserved and that boxing displays are the ancestral state. Moreover, we found that bats with a free-hanging resting posture were more likely to exhibit boxing displays than pushing displays. In addition, bats with longer forearms do not have a higher propensity for boxing displays. This study expands our limited knowledge of the evolution of agonistic displays and highlights the importance of resting posture as a driving force in the diversity of agonistic displays.

Vocal learning-associated convergent evolution in mammalian proteins and regulatory elements.

Vocal production learning ("vocal learning") is a convergently evolved trait in vertebrates. To identify brain genomic elements associated with mammalian vocal learning, we integrated genomic, anatomical, and neurophysiological data from the Egyptian fruit bat (Rousettus aegyptiacus) with analyses of the genomes of 215 placental mammals. First, we identified a set of proteins evolving more slowly in vocal learners. Then, we discovered a vocal motor cortical region in the Egyptian fruit bat, an emergent vocal learner, and leveraged that knowledge to identify active cis-regulatory elements in the motor cortex of vocal learners. Machine learning methods applied to motor cortex open chromatin revealed 50 enhancers robustly associated with vocal learning whose activity tended to be lower in vocal learners. Our research implicates convergent losses of motor cortex regulatory elements in mammalian vocal learning evolution.

Karyotypic stasis and swarming influenced the evolution of viral tolerance in a species-rich bat radiation.

The emergence of COVID-19 and severe acute respiratory syndrome (SARS) has prioritized understanding bats' viral tolerance. Myotis bats are exceptionally species rich and have evolved viral tolerance. They also exhibit swarming, a cryptic behavior where large, multi-species assemblages gather for mating, which has been hypothesized to promote interspecific hybridization. To resolve the coevolution of genome architecture and their unusual antiviral tolerance, we undertook a phylogenomic analysis of 60 Old World Myotis genomes. We demonstrate an extensive history of introgressive hybridization that has replaced the species phylogeny across 17%-93% of the genome except for pericentromeric regions of macrochromosomes. Introgression tracts were enriched on microchromosome regions containing key antiviral pathway genes overexpressed during viral challenge experiments. Together, these results suggest that the unusual Myotis karyotype may have evolved to selectively position immune-related genes in high recombining genomic regions prone to introgression of divergent alleles, including a diversity of interleukin loci responsible for the release of pro-inflammatory cytokines.

Species-specific SNP arrays for non-invasive genetic monitoring of a vulnerable bat.

Genetic tagging from scats is one of the minimally invasive sampling (MIS) monitoring approaches commonly used to guide management decisions and evaluate conservation efforts. Microsatellite markers have traditionally been used but are prone to genotyping errors. Here, we present a novel method for individual identification in the Threatened ghost bat Macroderma gigas using custom-designed Single Nucleotide Polymorphism (SNP) arrays on the MassARRAY system. We identified 611 informative SNPs from DArTseq data from which three SNP panels (44-50 SNPs per panel) were designed. We applied SNP genotyping and molecular sexing to 209 M. gigas scats collected from seven caves in the Pilbara, Western Australia, employing a two-step genotyping protocol and identifying unique genotypes using a custom-made R package, ScatMatch. Following data cleaning, the average amplification rate was 0.90 ± 0.01 and SNP genotyping errors were low (allelic dropout 0.003 ± 0.000) allowing clustering of scats based on one or fewer allelic mismatches. We identified 19 unique bats (9 confirmed/likely males and 10 confirmed/likely females) from a maternity and multiple transitory roosts, with two male bats detected using roosts, 9 km and 47 m apart. The accuracy of our SNP panels enabled a high level of confidence in the identification of individual bats. Targeted SNP genotyping is a valuable tool for monitoring and tracking of non-model species through a minimally invasive sampling approach.

Musculoskeletal morphogenesis supports the convergent evolution of bat laryngeal echolocation.

The order Chiroptera (bats) is the second largest group of mammals. One of the essential adaptations that have allowed bats to dominate the night skies is laryngeal echolocation, where bats emit ultrasonic pulses and listen to the returned echo to produce high-resolution 'images' of their surroundings. There are two possible scenarios for the evolutionary origin of laryngeal echolocation in bats: (1) a single origin in a common ancestor followed by the secondary loss in Pteropodidae, or (2) two convergent origins in Rhinolophoidea and Yangochiroptera. Although data from palaeontological, anatomical, developmental and genomic studies of auditory apparatuses exist, they remain inconclusive concerning the evolutionary origin of bat laryngeal echolocation. Here we compared musculoskeletal morphogenesis of the larynx in several chiropteran lineages and found distinct laryngeal modifications in two echolocating lineages, rhinolophoids and yangochiropterans. Our findings support the second scenario that rhinolophoids and yangochiropterans convergently evolved advanced laryngeal echolocation through anatomical modifications of the larynx for ultrasonic sound generation and refinement of the auditory apparatuses for more detailed sound perception.

A minimally invasive, field-applicable CRISPR/Cas biosensor to aid in the detection of Pseudogymnoascus destructans, the causative fungal agent of white-nose syndrome in bats.

The accessibility to CRISPR/Cas (Clustered Regularly Interspaced Short Palindromic Repeats/CRISPR-associated protein) genetic tools has given rise to applications beyond site-directed genome editing for the detection of DNA and RNA. These tools include precise diagnostic detection of human disease pathogens, such as SARS-CoV-2 and Zika virus. Despite the technology being rapid and cost-effective, the use of CRISPR/Cas tools in the surveillance of the causative agents of wildlife diseases has not been prominent. This study presents the development of a minimally invasive, field-applicable and user-friendly CRISPR/Cas-based biosensor for the detection of Pseudogymnoascus destructans (Pd), the causative fungal agent of white-nose syndrome (WNS), an infectious disease that has killed more than five million bats in North America since its discovery in 2006. The biosensor assay combines a recombinase polymerase amplification (RPA) step followed by CRISPR/Cas12a nuclease cleavage to detect Pd DNA from bat dermal swab and guano samples. The biosensor had similar detection results when compared to quantitative PCR in distinguishing Pd-positive versus negative field samples. Although bat dermal swabs could be analysed with the biosensor without nucleic acid extraction, DNA extraction was needed when screening guano samples to overcome inhibitors. This assay can be applied to help with more rapid delineation of Pd-positive sites in the field to inform management decisions. With further optimization, this technology has broad translation potential to wildlife disease-associated pathogen detection and monitoring applications.