A metaanalysis of bat phylogenetics and positive selection based on genomes and transcriptomes from 18 species

A genomic hotspot of diversifying selection and structural change in the hoary bat (Lasiurus cinereus)

Background

Previous work found that numerous genes positively selected within the hoary bat (Lasiurus cinereus) lineage are physically clustered in regions of conserved synteny. Here I further validate and expand on those finding utilizing an updated L. cinereus genome assembly and additional bat species as well as other tetrapod outgroups.

Methods

A chromosome-level assembly was generated by chromatin-contact mapping and made available by DNAZoo (www.dnazoo.org). The genomic organization of orthologous genes was extracted from annotation data for multiple additional bat species as well as other tetrapod clades for which chromosome-level assemblies were available from the National Center for Biotechnology Information (NCBI). Tests of branch-specific positive selection were performed for L. cinereus using PAML as well as with the HyPhy package for comparison.

Results

Twelve genes exhibiting significant diversifying selection in the L. cinereus lineage were clustered within a 12-Mb genomic window; one of these (Trpc4) also exhibited diversifying selection in bats generally. Ten of the 12 genes are landmarks of two distinct blocks of ancient synteny that are not linked in other tetrapod clades. Bats are further distinguished by frequent structural rearrangements within these synteny blocks, which are rarely observed in other Tetrapoda. Patterns of gene order and orientation among bat taxa are incompatible with phylogeny as presently understood, implying parallel evolution or subsequent reversals. Inferences of positive selection were found to be robust to alternative phylogenetic topologies as well as a strong shift in background nucleotide composition in some taxa.

Discussion

This study confirms and further localizes a genomic hotspot of protein-coding divergence in the hoary bat, one that also exhibits an increased tempo of structural change in bats compared with other mammals. Most genes in the two synteny blocks have elevated expression in brain tissue in humans and model organisms, and genetic studies implicate the selected genes in cranial and neurological development, among other functions.

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.

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.

Responses to predation pressure involve similar sets of genes in two divergent species of Daphnia

Species that are not closely related can express similar inducible traits, but molecular mechanisms underlying the observed responses are often unknown, nor is it known if these mechanisms are shared between such species. Here, we compared transcriptional profiles of two Daphnia species (D. mitsukuri and D. sinensis) from different subgenera, at both juvenile and adult developmental stages. Both species were exposed to the same predation threat (fish kairomones), and both showed similar induced morphological changes (reduced body length). At the early developmental stage, response to predation risk resulted in similar changes in expression levels of 23 orthologues in both species. These orthologues, involved in 107 GO categories, changed in the same direction in both species (over- or underexpressed), in comparison to non-exposed controls. Several of these orthologues were associated with DNA replication, structural constituents of cuticle or innate immune response. In both species, the differentially expressed (DE) genes on average had higher ω (dN /dS ) values than non-DE genes, suggesting that these genes had experienced greater positive selection or lower purifying selection than non-DE genes. Overall, our results suggest that similar suites of genes, responding in similar ways to predation pressure, have been retained in Daphnia for many millions of years.

Antiviral effects of interferon-stimulated genes in bats

The interferon pathway is the first line of defense in viral infection in all mammals, and its induction stimulates broad expression of interferon-stimulated genes (ISGs). In mice and also humans, the antiviral function of ISGs has been extensively studied. As an important viral reservoir in nature, bats can coexist with a variety of pathogenic viruses without overt signs of disease, yet only limited data are available for the role of ISGs in bats. There are multiple species of bats and work has begun deciphering the differences and similarities between ISG function of human/mouse and different bat species. This review summarizes the current knowledge of conserved and bat-specific-ISGs and their known antiviral effector functions.

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.

Bat ASC2 suppresses inflammasomes and ameliorates inflammatory diseases

Bats are special in their ability to live long and host many emerging viruses. Our previous studies showed that bats have altered inflammasomes, which are central players in aging and infection. However, the role of inflammasome signaling in combating inflammatory diseases remains poorly understood. Here, we report bat ASC2 as a potent negative regulator of inflammasomes. Bat ASC2 is highly expressed at both the mRNA and protein levels and is highly potent in inhibiting human and mouse inflammasomes. Transgenic expression of bat ASC2 in mice reduced the severity of peritonitis induced by gout crystals and ASC particles. Bat ASC2 also dampened inflammation induced by multiple viruses and reduced mortality of influenza A virus infection. Importantly, it also suppressed SARS-CoV-2-immune-complex-induced inflammasome activation. Four key residues were identified for the gain of function of bat ASC2. Our results demonstrate that bat ASC2 is an important negative regulator of inflammasomes with therapeutic potential in inflammatory diseases.

Reconstructing the dietary habits and trophic positions of the Longipterygidae (Aves: Enantiornithes) using neontological and comparative morphological methods

The Longipterygidae are a unique clade among the enantiornithines in that they exhibit elongate rostra (≥60% total skull length) with dentition restricted to the distal tip of the rostrum, and pedal morphologies suited for an arboreal lifestyle (as in other enantiornithines). This suite of features has made interpretations of this group’s diet and ecology difficult to determine due to the lack of analogous taxa that exhibit similar morphologies together. Many extant bird groups exhibit rostral elongation, which is associated with several disparate ecologies and diets (e.g., aerial insectivory, piscivory, terrestrial carnivory). Thus, the presence of rostral elongation in the Longipterygidae only somewhat refines trophic predictions of this clade. Anatomical morphologies do not function singularly but as part of a whole and thus, any dietary or ecological hypothesis regarding this clade must also consider other features such as their unique dentition. The only extant group of dentulous volant tetrapods are the chiropterans, in which tooth morphology and enamel thickness vary depending upon food preference. Drawing inferences from both avian bill proportions and variations in the dental morphology of extinct and extant taxa, we provide quantitative data to support the hypothesis that the Longipterygidae were animalivorous, with greater support for insectivory.

Genetic Characterization and Insular Habitat Enveloping of Endangered Leaf-Nosed Bat, Hipposideros nicobarulae (Mammalia: Chiroptera) in India: Phylogenetic Inference and Conservation Implication

The Nicobar leaf-nosed Bat (Hipposideros nicobarulae) was described in the early 20th century; however, its systematic classification has been debated for over 100 years. This endangered and endemic species has achieved species status through morphological data in the last 10 years. However, the genetic information and phylogenetic relationships of H. nicobarulae remain neglected. The generated mitochondrial cytochrome b gene (mtCytb) sequences (438 bp) of H. nicobarulae contains 53.42-53.65% AT composition and 1.82% variable sites. The studied species, H. nicobarulae maintains an 8.1% to 22.6% genetic distance from other Hipposideros species. The genetic divergence estimated in this study is congruent with the concept of gene speciation in bats. The Bayesian and Maximum-Likelihood phylogenies clearly discriminated all Hipposideros species and showed a sister relationship between H. nicobarulae and H. cf. antricola. Current mtCytb-based investigations of H. nicobarulae have confirmed the species status at the molecular level. Further, the MaxEnt-based species distribution modelling illustrates the most suitable habitat of H. nicobarulae (294 km²), of which the majority (171 km²) is located on Great Nicobar Island. The present study suggests rigorous sampling across the range, taxonomic coverage, the generation of multiple molecular markers (mitochondrial and nuclear), as well as more ecological information, which will help in understanding population genetic structure, habitat suitability, and the implementation of appropriate conservation action plans for H. nicobarulae and other Hipposideros species.

Exploring Evolutionary Relationships within Neodermata Using Putative Orthologous Groups of Proteins, with Emphasis on Peptidases

The phylogenetic relationships within Neodermata were examined based on putative orthologous groups of proteins (OGPs) from 11 species of Monogenea, Trematoda, and Cestoda. The dataset included OGPs from BUSCO and OMA. Additionally, peptidases were identified and evaluated as phylogenetic markers. Phylogenies were inferred using the maximum likelihood method. A network analysis and a hierarchical grouping analysis of the principal components (HCPC) of orthologous groups of peptidases were performed. The phylogenetic analyses showed the monopisthocotylean monogeneans as the sister-group of cestodes, and the polyopisthocotylean monogeneans as the sister-group of trematodes. However, the sister-group relationship between Monopisthocotylea and Cestoda was not statistically well supported. The network analysis and HCPC also showed a cluster formed by polyopisthocotyleans and trematodes. The present study supports the non-monophyly of Monogenea. An analysis of mutation rates indicated that secreted peptidases and inhibitors, and those with multiple copies, are under positive selection pressure, which could explain the expansion of some families such as C01, C19, I02, and S01. Whilst not definitive, our study presents another point of view in the discussion of the evolution of Neodermata, and we hope that our data drive further discussion and debate on this intriguing topic.

Evolutionary pressures of aerial insectivory reflected in anurognathid pterosaurs

Across the evolution of powered flight, the ecological niche of aerial insectivore has been occupied by members of the three volant vertebrate clades-Aves and Chiroptera, and the first known volant vertebrates, pterosaurs. However, morphological and quantitative evidence to support pterosaurs exhibiting this ecology remains scant. Anurognathids are an unusual group of pterosaurs in which the skull exhibits the unique morphology of being mediolaterally expanded, so much so that their skulls may be wider than rostrocaudally long. Here, we conduct quantitative comparative cranial measurements and dental morphology in anurognathids against extant avian and chiropteran taxa, respectively, with ecologies and behaviors that are similar to predicted putative behaviors of anurognathids. Comparative analyses of both skull and dental morphology suggest anurognathid specimens in similar morphospaces as insectivorous crepuscular and nocturnal extant volant taxa. Our results support that this unique group of pterosaurs likely occupied a niche of mid-flight insectivorous capture in low-light conditions.

Adaptive duplication and genetic diversification of protein kinase R contribute to the specificity of bat-virus interactions

Several bat species act as asymptomatic reservoirs for many viruses that are highly pathogenic in other mammals. Here, we have characterized the functional diversification of the protein kinase R (PKR), a major antiviral innate defense system. Our data indicate that PKR has evolved under positive selection and has undergone repeated genomic duplications in bats in contrast to all studied mammals that have a single copy of the gene. Functional testing of the relationship between PKR and poxvirus antagonists revealed how an evolutionary conflict with ancient pathogenic poxviruses has shaped a specific bat host-virus interface. We determined that duplicated PKRs of the Myotis species have undergone genetic diversification, allowing them to collectively escape from and enhance the control of DNA and RNA viruses. These findings suggest that viral-driven adaptations in PKR contribute to modern virus-bat interactions and may account for bat-specific immunity.

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.

Coronavirus Infections in Animals: Risks of Direct and Reverse Zoonoses

The publications on animal coronavirus infections that have the greatest emerging potential, as well as official data from the World Organization for Animal Health (OIE) on cases of animal infection with COVID-19, are analyzed. Like most infectious diseases common to humans, coronavirus infections were first discovered in animals. Due to the increased rate of replication and recombination activity compared to other viruses, mutations occur more often in the genome of coronaviruses, which contribute to the acquisition of new qualities in order to consolidate in the host organism. Examples of cross-species transmission are not only SARS-CoV, MERS-CoV, and SARS-CoV-2, which are dangerous to humans, but also coronaviruses of agricultural and domestic animals, between which there is a genetic relationship. There are several known cases of zoo, wild, domestic, and farm animals displaying symptoms characteristic of COVID-19 and identification of the genome of the SARS-CoV-2 virus in them. The issue of cross-species transmission of coronavirus infections, in particular the reverse zoonosis of SARS-CoV-2 from animals to humans, is widely discussed. According to the conclusions of many researchers, including OIE experts, there is no direct evidence base for infection of humans with COVID-19 from animals. However, people with suspected COVID-19 and with a confirmed diagnosis are still advised to isolate not only from people but also from animals. A number of methods for specific prevention, diagnosis, and immunization against a wide range of coronavirus infections are being developed at the All-Russia Research Institute for Animal Protection.

Chirohepevirus from Bats: Insights into Hepatitis E Virus Diversity and Evolution

Homologs of the human hepatitis E virus (HEV) have been identified in more than a dozen animal species. Some of them have been evidenced to cross species barriers and infect humans. Zoonotic HEV infections cause chronic liver diseases as well as a broad range of extrahepatic manifestations, which increasingly become significant clinical problems. Bats comprise approximately one-fifth of all named mammal species and are unique in their distinct immune response to viral infection. Most importantly, they are natural reservoirs of several highly pathogenic viruses, which have induced severe human diseases. Since the first discovery of HEV-related viruses in bats in 2012, multiple genetically divergent HEV variants have been reported in a total of 12 bat species over the last decade, which markedly expanded the host range of the HEV family and shed light on the evolutionary origin of human HEV. Meanwhile, bat-borne HEV also raised critical public health concerns about its zoonotic potential. Bat HEV strains resemble genomic features but exhibit considerable heterogeneity. Due to the close evolutionary relationships, bat HEV altogether has been recently assigned to an independent genus, Chirohepevirus. This review focuses on the current state of bat HEV and provides novel insights into HEV genetic diversity and molecular evolution.

Coevolutionary Analysis Implicates Toll-Like Receptor 9 in Papillomavirus Restriction

Upon infection, DNA viruses can be sensed by pattern recognition receptors (PRRs), leading to the activation of type I and III interferons to block infection. Therefore, viruses must inhibit these signaling pathways, avoid being detected, or both. Papillomavirus virions are trafficked from early endosomes to the Golgi apparatus and wait for the onset of mitosis to complete nuclear entry. This unique subcellular trafficking strategy avoids detection by cytoplasmic PRRs, a property that may contribute to the establishment of infection. However, as the capsid uncoats within acidic endosomal compartments, the viral DNA may be exposed to detection by Toll-like receptor 9 (TLR9). In this study, we characterized two new papillomaviruses from bats and used molecular archeology to demonstrate that their genomes altered their nucleotide compositions to avoid detection by TLR9, providing evidence that TLR9 acts as a PRR during papillomavirus infection. Furthermore, we showed that TLR9, like other components of the innate immune system, is under evolutionary selection in bats, providing the first direct evidence for coevolution between papillomaviruses and their hosts. Finally, we demonstrated that the cancer-associated human papillomaviruses show a reduction in CpG dinucleotides within a TLR9 recognition complex.

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

Background

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

Methods

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

Results

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

Next-Generation Sequencing for Confronting Virus Pandemics

Virus pandemics have happened, are happening and will happen again. In recent decades, the rate of zoonotic viral spillover into humans has accelerated, mirroring the expansion of our global footprint and travel network, including the expansion of viral vectors and the destruction of natural spaces, bringing humans closer to wild animals. Once viral cross-species transmission to humans occurs, transmission cannot be stopped by cement walls but by developing barriers based on knowledge that can prevent or reduce the effects of any pandemic. Controlling a local transmission affecting few individuals is more efficient that confronting a community outbreak in which infections cannot be traced. Genetic detection, identification, and characterization of infectious agents using next-generation sequencing (NGS) has been proven to be a powerful tool allowing for the development of fast PCR-based molecular assays, the rapid development of vaccines based on mRNA and DNA, the identification of outbreaks, transmission dynamics and spill-over events, the detection of new variants and treatment of vaccine resistance mutations, the development of direct-acting antiviral drugs, the discovery of relevant minority variants to improve knowledge of the viral life cycle, strengths and weaknesses, the potential for becoming dominant to take appropriate preventive measures, and the discovery of new routes of viral transmission.

TR Locus Annotation and Characteristics of Rhinolophus ferrumequinum

T-cell antigen receptors (TRs) in vertebrates can be divided into αβ or γδ, encoded by TRA/D, TRG, or TRB loci. TRs play a central role in mammal cellular immunity, which occurs by rearrangement of V, D, J, and C genes in the loci. The bat is the only mammal with flying ability and is considered the main host of zoonotic viruses, an important public health concern. However, at present, little is known about the composition of bat TR genes. Based on the whole genome sequence of the greater horseshoe bat (Rhinolophus ferrumequinum) and referring to the TR/IG annotation rules formulated by the international ImMunoGeneTics information system (IMGT), we present a complete annotation of TRA/D, TRG, and TRB loci of R. ferrumequinum. A total of 128 V segments, three D segments, 85 J segments, and 6 C segments were annotated and compared with other known mammalian data. The characteristics of the TR locus and germline genes of R. ferrumequinum are analyzed.

Of bats and men: Immunomodulatory treatment options for COVID-19 guided by the immunopathology of SARS-CoV-2 infection

[Figure: see text].

Positive Selection and Gene Expression Analyses from Salivary Glands Reveal Discrete Adaptations within the Ecologically Diverse Bat Family Phyllostomidae

The leaf-nosed bats (Phyllostomidae) are outliers among chiropterans with respect to the unusually high diversity of dietary strategies within the family. Salivary glands, owing to their functions and high ultrastructural variability among lineages, are proposed to have played an important role during the phyllostomid radiation. To identify genes underlying salivary gland functional diversification, we sequenced submandibular gland transcriptomes from phyllostomid species representative of divergent dietary strategies. From the assembled transcriptomes, we performed an array of selection tests and gene expression analyses to identify signatures of adaptation. Overall, we identified an enrichment of immunity-related gene ontology terms among 53 genes evolving under positive selection. Lineage-specific selection tests revealed several endomembrane system genes under selection in the vampire bat. Many genes that respond to insulin were under selection and differentially expressed genes pointed to modifications of amino acid synthesis pathways in plant-visitors. Results indicate salivary glands have diversified in various ways across a functional diverse clade of mammals in response to niche specializations.

Coronavirus, the King Who Wanted More Than a Crown: From Common to the Highly Pathogenic SARS-CoV-2, Is the Key in the Accessory Genes?

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), that emerged in late 2019, is the etiologic agent of the current "coronavirus disease 2019" (COVID-19) pandemic, which has serious health implications and a significant global economic impact. Of the seven human coronaviruses, all of which have a zoonotic origin, the pandemic SARS-CoV-2, is the third emerging coronavirus, in the 21st century, highly pathogenic to the human population. Previous human coronavirus outbreaks (SARS-CoV-1 and MERS-CoV) have already provided several valuable information on some of the common molecular and cellular mechanisms of coronavirus infections as well as their origin. However, to meet the new challenge caused by the SARS-CoV-2, a detailed understanding of the biological specificities, as well as knowledge of the origin are crucial to provide information on viral pathogenicity, transmission and epidemiology, and to enable strategies for therapeutic interventions and drug discovery. Therefore, in this review, we summarize the current advances in SARS-CoV-2 knowledges, in light of pre-existing information of other recently emerging coronaviruses. We depict the specificity of the immune response of wild bats and discuss current knowledge of the genetic diversity of bat-hosted coronaviruses that promotes viral genome expansion (accessory gene acquisition). In addition, we describe the basic virology of coronaviruses with a special focus SARS-CoV-2. Finally, we highlight, in detail, the current knowledge of genes and accessory proteins which we postulate to be the major keys to promote virus adaptation to specific hosts (bat and human), to contribute to the suppression of immune responses, as well as to pathogenicity.

Malignant Peripheral Nerve Sheath Tumour in a Seba's Short-Tailed Bat (Carollia perspicillata)

There are few reports of tumours in bats, despite the fact that members of the Chiroptera order include 22.4% of all currently known mammal species. We now describe a case of malignant peripheral nerve sheath tumour (MPNST) in an adult female Seba's short-tailed bat (Carollia perspicillata) from a private zoo in Switzerland. The macroscopical and histological findings indicated that the MPNST originated from a dorsal spinal nerve root sheath in the lumbar region and metastasized to the spleen. Immunohistochemical labelling for glial fibrillary acid protein, S100 protein, neuron specific enolase and synaptophysin was negative. The cytoplasm of a few neoplastic cells labelled for smooth muscle actin and desmin. To our knowledge, this is the first report of a PNST in a member of the Chiroptera order.

Lessons from the host defences of bats, a unique viral reservoir

There have been several major outbreaks of emerging viral diseases, including Hendra, Nipah, Marburg and Ebola virus diseases, severe acute respiratory syndrome (SARS) and Middle East respiratory syndrome (MERS)-as well as the current pandemic of coronavirus disease 2019 (COVID-19). Notably, all of these outbreaks have been linked to suspected zoonotic transmission of bat-borne viruses. Bats-the only flying mammal-display several additional features that are unique among mammals, such as a long lifespan relative to body size, a low rate of tumorigenesis and an exceptional ability to host viruses without presenting clinical disease. Here we discuss the mechanisms that underpin the host defence system and immune tolerance of bats, and their ramifications for human health and disease. Recent studies suggest that 64 million years of adaptive evolution have shaped the host defence system of bats to balance defence and tolerance, which has resulted in a unique ability to act as an ideal reservoir host for viruses. Lessons from the effective host defence of bats would help us to better understand viral evolution and to better predict, prevent and control future viral spillovers. Studying the mechanisms of immune tolerance in bats could lead to new approaches to improving human health. We strongly believe that it is time to focus on bats in research for the benefit of both bats and humankind.

Coronavirus Infections of Animals: Future Risks to Humans

Coronaviruses have tremendous evolutionary potential, and three major outbreaks of new human coronavirus infections have occurred in the recent history of humankind. In this paper, the patterns of occurrence of new zoonotic coronavirus infections and the role of bioveterinary control in preventing their potential outbreaks in the future are determined. The possibility of SARS-CoV-2 infection in companion animals is considered. Diverse human activities may trigger various interactions between animal species and their viruses, sometimes causing the emergence of new viral pathogens. In addition, the possibility of using probiotics for the control of viral infections in animals is discussed.

Fundamental Characteristics of Bat Interferon Systems

Interferons are an essential component of the innate arm of the immune system and are arguably one of the most important lines of defence against viruses. The human IFN system and its functionality has already been largely characterized and studied in detail. However, the IFN systems of bats have only been marginally examined to date up until the recent developments of the Bat1k project which have now opened new opportunities in research by identifying six new bat genomes to possess novel genes that are likely associated with viral tolerance exhibited in bats. Interestingly, bats have been hypothesized to possess the ability to establish a host-virus relationship where despite being infected, they exhibit limited signs of disease and still retain the ability to transmit the disease into other susceptible hosts. Bats are one of the most abundant and widespread vertebrates on the planet and host many zoonotic viruses that are highly pathogenic to humans. Several genomics, immunological, and biological features are thought to underlie novel antiviral mechanisms of bats. This review aims to explore the bat IFN system and developments in its diverse IFN features, focusing mainly on the model species, the Australian black flying fox (Pteropus alecto), while also highlighting bat innate immunity as an exciting and fruitful area of research to understand their ability to control viral-mediated pathogenesis.

Microchiropterans have a diminutive cerebral cortex, not an enlarged cerebellum, compared to megachiropterans and other mammals

Small echolocating bats are set apart from most other mammals by their relatively large cerebellum, a feature that has been associated to echolocation, as it is presumed to indicate a relatively enlarged number of neurons in the cerebellum in comparison to other brain structures. Here we quantify the neuronal composition of the cerebral cortex, cerebellum and remaining brain structures of seven species of large Pteropodid bats (formerly classified as megachiropterans), one of which echolocates, and six species of small bats (formerly classified as microchiropterans), all of which echolocate. This chiropteran data is compared to 60 mammalian species in our dataset to determine whether the relatively large cerebellum of the small echolocating bats, and possibly that of the echolocating Pteropodid, contains a relatively enlarged number of neurons. We find no evidence that the distribution of neurons differs between microchiropterans and megachiropterans, but our data indicate that microchiropterans, like the smallest shrew in our dataset, have diminutive cerebral cortices, which makes the cerebellum appear relatively large. We propose that, in agreement with the diminutive brain size of the earliest fossil mammals, this is a plesiomorphic, not a derived, feature of microchiropteran brains. The results of this study also reveal important neural characteristics related to the phylogenetic affinities and relationships of the chiropterans.

RTP4 Is a Potent IFN-Inducible Anti-flavivirus Effector Engaged in a Host-Virus Arms Race in Bats and Other Mammals

Among mammals, bats are particularly rich in zoonotic viruses, including flaviviruses. Certain bat species can be productively yet asymptomatically infected with viruses that cause overt disease in other species. However, little is known about the antiviral effector repertoire in bats relative to other mammals. Here, we report the black flying fox receptor transporter protein 4 (RTP4) as a potent interferon (IFN)-inducible inhibitor of human pathogens in the Flaviviridae family, including Zika, West Nile, and hepatitis C viruses. Mechanistically, RTP4 associates with the flavivirus replicase, binds viral RNA, and suppresses viral genome amplification. Comparative approaches revealed that RTP4 undergoes positive selection, that a flavivirus can mutate to escape RTP4-imposed restriction, and that diverse mammalian RTP4 orthologs exhibit striking patterns of specificity against distinct Flaviviridae members. Our findings reveal an antiviral mechanism that has likely adapted over 100 million years of mammalian evolution to accommodate unique host-virus genetic conflicts.

Complementary regulation of caspase-1 and IL-1β reveals additional mechanisms of dampened inflammation in bats

Bats have emerged as unique mammalian vectors harboring a diverse range of highly lethal zoonotic viruses with minimal clinical disease. Despite having sustained complete genomic loss of AIM2, regulation of the downstream inflammasome response in bats is unknown. AIM2 sensing of cytoplasmic DNA triggers ASC aggregation and recruits caspase-1, the central inflammasome effector enzyme, triggering cleavage of cytokines such as IL-1β and inducing GSDMD-mediated pyroptotic cell death. Restoration of AIM2 in bat cells led to intact ASC speck formation, but intriguingly resulted in a lack of caspase-1 or consequent IL-1β activation. We further identified two residues undergoing positive selection pressures in Pteropus alecto caspase-1 that abrogate its enzymatic function and are crucial in human caspase-1 activity. Functional analysis of another bat lineage revealed a targeted mechanism for loss of Myotis davidii IL-1β cleavage and elucidated an inverse complementary relationship between caspase-1 and IL-1β, resulting in overall diminished signaling across bats of both suborders. Thus we report strategies that additionally undermine downstream inflammasome signaling in bats, limiting an overactive immune response against pathogens while potentially producing an antiinflammatory state resistant to diseases such as atherosclerosis, aging, and neurodegeneration.

Evolutionary medical insights into the SARS-CoV-2 pandemic

The author apply concepts and tools from evolutionary medicine to understanding the SARS-CoV-2 pandemic. The pandemic represents a mismatched conflict, with dynamics and pathology apparently driven by three main factors: (i) bat immune systems that rely on low inflammation but high efficacy of interferon-based defenses; (ii) viral tactics that differentially target the human interferon system, leading to substantial asymptomatic and pre-symptomatic transmission; and (ii) high mortality caused by hyper-inflammatory and hyper-coagulatory phenotypes, that represent dysregulated tradeoffs whereby collateral immune-induced damage becomes systemic and severe. This framework can explain the association of mortality with age (which involves immune life-history shifts towards higher inflammation and coagulation and reduced adaptive immunity), and sex (since males senesce faster than females). Genetic-risk factors for COVID-19 mortality can be shown, from a phenome-wide association analysis of the relevant SNPs, to be associated with inflammation and coagulation; the phenome-wide association study also provides evidence, consistent with several previous studies, that the calcium channel blocking drug amlodipine mediates risk of mortality. Lay Summary: SARS-CoV-2 is a bat virus that jumped into humans. The virus is adapted to bat immune systems, where it evolved to suppress the immune defenses (interferons) that mammals use to tell that they are infected. In humans, the virus can apparently spread effectively in the body with a delay in the production of symptoms and the initiation of immune responses. This delay may then promote overactive immune responses, when the virus is detected, that damage the body as a side effect. Older people are more vulnerable to the virus because they are less adapted to novel infectious agents, and invest less in immune defense, compared to younger people. Genes that increase risk of mortality from SARS-CoV-2 are functionally associated with a drug called amlodipine, which may represent a useful treatment.

DGINN, an automated and highly-flexible pipeline for the detection of genetic innovations on protein-coding genes

Adaptive evolution has shaped major biological processes. Finding the protein-coding genes and the sites that have been subjected to adaptation during evolutionary time is a major endeavor. However, very few methods fully automate the identification of positively selected genes, and widespread sources of genetic innovations such as gene duplication and recombination are absent from most pipelines. Here, we developed DGINN, a highly-flexible and public pipeline to Detect Genetic INNovations and adaptive evolution in protein-coding genes. DGINN automates, from a gene's sequence, all steps of the evolutionary analyses necessary to detect the aforementioned innovations, including the search for homologs in databases, assignation of orthology groups, identification of duplication and recombination events, as well as detection of positive selection using five methods to increase precision and ranking of genes when a large panel is analyzed. DGINN was validated on nineteen genes with previously-characterized evolutionary histories in primates, including some engaged in host-pathogen arms-races. Our results confirm and also expand results from the literature, including novel findings on the Guanylate-binding protein family, GBPs. This establishes DGINN as an efficient tool to automatically detect genetic innovations and adaptive evolution in diverse datasets, from the user's gene of interest to a large gene list in any species range.

Echolocation call frequency variation in horseshoe bats: molecular basis revealed by comparative transcriptomics

Recently diverged taxa with contrasting phenotypes offer opportunities for unravelling the genetic basis of phenotypic variation in nature. Horseshoe bats are a speciose group that exhibit a derived form of high-duty cycle echolocation in which the inner ear is finely tuned to echoes of the narrowband call frequency. Here, by focusing on three recently diverged subspecies of the intermediate horseshoe bat (Rhinolophus affinis) that display divergent echolocation call frequencies, we aim to identify candidate loci putatively involved in hearing frequency variation. We used de novo transcriptome sequencing of two mainland taxa (himalayanus and macrurus) and one island taxon (hainanus) to compare expression profiles of thousands of genes. By comparing taxa with divergent call frequencies (around 15 kHz difference), we identified 252 differentially expressed genes, of which six have been shown to be involved in hearing or deafness in human/mouse. To obtain further validation of these results, we applied quantitative reverse transcription-PCR to the candidate gene FBXL15 and found a broad association between the level of expression and call frequency across taxa. The genes identified here represent strong candidate loci associated with hearing frequency variation in bats.

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

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

Analysis of Natural Selection of Immune Genes in Spinibarbus caldwelli by Transcriptome Sequencing

Spinibarbus caldwelli is an omnivorous cyprinid fish that is distributed widely in China. To investigate the adaptive evolution of S. caldwelli, the muscle transcriptome was sequenced by Illumina HiSeq 4000 platform. A total of 80,447,367 reads were generated by next-generation sequencing. Also, 211,386 unigenes were obtained by de novo assembly. Additionally, we calculated that the divergence time between S. caldwelli and Sinocyclocheilus grahami is 23.14 million years ago (Mya). And both of them diverged from Ctenopharyngodon idellus 46.95 Mya. Furthermore, 38 positive genes were identified by calculating Ka/Ks ratios from 9225 orthologs. Among them, several immune-related genes were identified as positively selected, such as POLR3B, PIK3C3, TOPORS, FASTKD3, CYPLP1A1, and UACA. Our results throw light on the nature of the natural selection of S. caldwelli and contribute to future immunological and transcriptome studies.

Transcriptome Ortholog Alignment Sequence Tools (TOAST) for phylogenomic dataset assembly

Background

Advances in next-generation sequencing technologies have reduced the cost of whole transcriptome analyses, allowing characterization of non-model species at unprecedented levels. The rapid pace of transcriptomic sequencing has driven the public accumulation of a wealth of data for phylogenomic analyses, however lack of tools aimed towards phylogeneticists to efficiently identify orthologous sequences currently hinders effective harnessing of this resource.

Results

We introduce TOAST, an open source R software package that can utilize the ortholog searches based on the software Benchmarking Universal Single-Copy Orthologs (BUSCO) to assemble multiple sequence alignments of orthologous loci from transcriptomes for any group of organisms. By streamlining search, query, and alignment, TOAST automates the generation of locus and concatenated alignments, and also presents a series of outputs from which users can not only explore missing data patterns across their alignments, but also reassemble alignments based on user-defined acceptable missing data levels for a given research question.

Conclusions

TOAST provides a comprehensive set of tools for assembly of sequence alignments of orthologs for comparative transcriptomic and phylogenomic studies. This software empowers easy assembly of public and novel sequences for any target database of candidate orthologs, and fills a critically needed niche for tools that enable quantification and testing of the impact of missing data. As open-source software, TOAST is fully customizable for integration into existing or novel custom informatic pipelines for phylogenomic inference.

Software, a detailed manual, and example data files are available through github carolinafishes.github.io

Novel Insights Into Immune Systems of Bats

In recent years, viruses similar to those that cause serious disease in humans and other mammals have been detected in apparently healthy bats. These include filoviruses, paramyxoviruses, and coronaviruses that cause severe diseases such as Ebola virus disease, Marburg haemorrhagic fever and severe acute respiratory syndrome (SARS) in humans. The evolution of flight in bats seem to have selected for a unique set of antiviral immune responses that control virus propagation, while limiting self-damaging inflammatory responses. Here, we summarize our current understanding of antiviral immune responses in bats and discuss their ability to co-exist with emerging viruses that cause serious disease in other mammals. We highlight how this knowledge may help us to predict viral spillovers into new hosts and discuss future directions for the field.