Differential Evolution of Antiretroviral Restriction Factors in Pteropid Bats as Revealed by APOBEC3 Gene Complexity

RNA-Seq analysis reveals the long noncoding RNAs associated with immunity in wild Myotis myotis bats

BACKGROUND

Bats possess a uniquely adapted immune system that enables them to live with viral infections without the expected maladies. The molecular basis and regulation of bats' immune response is still not fully understood. Long non-coding RNAs (lncRNAs) represent an emerging class of molecules with critical regulatory roles in multiple biological processes, including immunity. We hypothesise that lncRNA-based regulation in bats may enable them to limit disease and live with viral pathogens.

RESULTS

We developed a lncRNA prediction pipeline to annotate the long non-coding transcriptome across multiple bat tissues and at the population level. Characterisation of our lncRNA dataset based on 100 blood transcriptomes from wild Myotis myotis bats revealed lower and more tissue-specific expression compared with coding genes, reduced GC content and shorter length distributions, consistent with lncRNA profiles observed in other species. Using WGCNA network analyses and gene ontology, we identified two mRNA-lncRNA co-expression modules in Myotis myotis associated with distinct immune response: one linked to T-cell activation and vial processes, and the other to inflammation. From these immune-related lncRNAs, we selected four candidates with high translational potential for regulating viral infections and inflammation. These include a newly identified lncRNA, BatLnc1, with potential antiviral functions; the M. myotis ortholog of TUG1, implicated in viral-host interactions; and well-known lncRNAs MALAT1 and NEAT1, recognised for their roles in inflammatory regulation.

CONCLUSIONS

We conducted the first ab initio prediction of lncRNAs in a non-model bat species, the wild-caught M. myotis. Our network analysis revealed significant variation in immune status among a subset of individuals, potentially due to pathogenic conditions. From these variations, we identified lncRNAs most likely associated with immune response in bats. This initial exploration lays the groundwork for future experimental validations of lncRNA functions, offering promising insights into their role in bat immunity.

Coevolution of Lentiviral Vif with Host A3F and A3G: Insights from Computational Modelling and Ancestral Sequence Reconstruction

The evolutionary arms race between host restriction factors and viral antagonists provides crucial insights into immune system evolution and viral adaptation. This study investigates the structural and evolutionary dynamics of the double-domain restriction factors A3F and A3G and their viral inhibitor, Vif, across diverse primate species. By constructing 3D structural homology models and integrating ancestral sequence reconstruction (ASR), we identified patterns of sequence diversity, structural conservation, and functional adaptation. Inactive CD1 (Catalytic Domain 1) domains displayed greater sequence diversity and more positive surface charges than active CD2 domains, aiding nucleotide chain binding and intersegmental transfer. Despite variability, the CD2 DNA-binding grooves remained structurally consistent with conserved residues maintaining critical functions. A3F and A3G diverged in loop 7' interaction strategies, utilising distinct molecular interactions to facilitate their roles. Vif exhibited charge variation linked to host species, reflecting its coevolution with A3 proteins. These findings illuminate how structural adaptations and charge dynamics enable both restriction factors and their viral antagonists to adapt to selective pressures. Our results emphasize the importance of studying structural evolution in host-virus interactions, with implications for understanding immune defense mechanisms, zoonotic risks, and viral evolution. This work establishes a foundation for further exploration of restriction factor diversity and coevolution across species.

[The diversification of the protein kinase R contributes to the specificity of bat-virus interactions]

Several bat species asymptomatically harbor certain viruses that are highly pathogenic in other mammals. The underlying mechanisms involve an evolutionary balance between tolerance and immune resistance to viral infections. However, how bats innate immunity has evolved in response to viruses remains to be elucidated. Here, we review the evolution of the protein kinase R (PKR) in bats, a major antiviral protein of vertebrate innate antiviral defense. Our recent results indicate that PKR has evolved under selective pressure and has undergone genomic duplications in bats, in contrast to all mammals studied, which possess only a single copy of the gene. The genetic changes in bat PKR are probably partly the result of genetic conflicts with ancient pathogenic poxviruses, shaping a bat-specific host–virus interface. Furthermore, the duplicated PKRs in Myotis species enable them to collectively escape viruses and enhance their viral control. These results suggest that viral adaptations of PKR contribute to the specificity of modern interactions between viruses and bats, and may explain unique antiviral mechanisms in bats.

APOBEC3-Related Editing and Non-Editing Determinants of HIV-1 and HTLV-1 Restriction

The apolipoprotein B mRNA editing enzyme catalytic polypeptide-like 3 (APOBEC3/A3) family of cytosine deaminases serves as a key innate immune barrier against invading retroviruses and endogenous retroelements. The A3 family's restriction activity against these parasites primarily arises from their ability to catalyze cytosine-to-uracil conversions, resulting in genome editing and the accumulation of lethal mutations in viral genomes. Additionally, non-editing mechanisms, including deaminase-independent pathways, such as blocking viral reverse transcription, have been proposed as antiviral strategies employed by A3 family proteins. Although viral factors can influence infection progression, the determinants that govern A3-mediated restriction are critical in shaping retroviral infection outcomes. This review examines the interactions between retroviruses, specifically human immunodeficiency virus type 1 and human T-cell leukemia virus type 1, and A3 proteins to better understand how editing and non-editing activities contribute to the trajectory of these retroviral infections.

Human and bats genome robustness under COSMIC mutational signatures

Carcinogenesis is an evolutionary process, and mutations can fix the selected phenotypes in selective microenvironments. Both normal and neoplastic cells are robust to the mutational stressors in the microenvironment to the extent that secure their fitness. To test the robustness of genes under a range of mutagens, we developed a sequential mutation simulator, Sinabro, to simulate single base substitution under a given mutational process. Then, we developed a pipeline to measure the robustness of genes and cells under those mutagenesis processes. We discovered significant human genome robustness to the APOBEC mutational signature SBS2, which is associated with viral defense mechanisms and is implicated in cancer. Robustness evaluations across over 70,000 sequences against 41 signatures showed higher resilience under signatures predominantly causing C-to-T (G-to-A) mutations. Principal component analysis indicates the GC content at the codon's wobble position significantly influences robustness, with increased resilience noted under transition mutations compared to transversions. Then, we tested our results in bats at extremes of the lifespan-to-mass relationship and found the long-lived bat is more robust to APOBEC than the short-lived one. By revealing robustness to APOBEC ranked highest in human (and bats with much more than number of APOBEC) genome, this work bolsters the key potential role of APOBECs in aging and cancer, as well as evolved countermeasures to this innate mutagenic process. It also provides the baseline of the human and bat genome robustness under mutational processes associated with aging and cancer.

A tale of endurance: bats, viruses and immune dynamics

The emergence of highly zoonotic viral infections has propelled bat research forward. The viral outbreaks including Hendra virus, Nipah virus, Marburg virus, Ebola virus, Rabies virus, Middle East respiratory syndrome coronavirus, SARS-CoV and the latest SARS-CoV-2 have been epidemiologically linked to various bat species. Bats possess unique immunological characteristics that allow them to serve as a potential viral reservoir. Bats are also known to protect themselves against viruses and maintain their immunity. Therefore, there is a need for in-depth understanding into bat-virus biology to unravel the major factors contributing to the coexistence and spread of viruses.

Evolvability of cancer-associated genes under APOBEC3A/B selection

Evolvability is an emergent hallmark of cancer that depends on intra-tumor heterogeneity and genetic variation. Mutations generated by APOBEC3 contribute to genetic variation and tumor evolvability. However, the influence of APOBEC3 on the evolvability of the genome and its differential impact on cancer genes versus non-cancer genes remains unclear. Analyzing over 40,000 human protein-coding transcripts, we identified distinct distribution patterns of APOBEC3A/B TC motifs between cancer and non-cancer genes, suggesting unique associations with cancer. Studying a bat species with numerous APOBEC3 genes, we found distinct motif patterns in orthologs of cancer genes compared to non-cancer genes, as in humans, suggesting APOBEC3 evolution to reduce impacts on the genome rather than the converse. Simulations confirmed that APOBEC3-induced heterogeneity enhances cancer evolution through bimodal patterns of mutations in certain classes of genes. Our results suggest the bimodal distribution of APOBEC-induced mutations can significantly increase cancer heterogeneity.

Screening for Viruses in Indigenous Greek Black Pigs

The successful advancement of xenotransplantation has led to the development of highly sensitive detection systems for the screening of potentially zoonotic viruses in donor pigs and preventing their transmission to the recipient. To validate these methods, genetically modified pigs generated for xenotransplantation, numerous minipigs and other pig breeds have been tested, thereby increasing our knowledge concerning the pig virome and the distribution of pig viruses. Of particular importance are the porcine cytomegalovirus, a porcine roseolovirus (PCMV/PRV) and the hepatitis E virus genotype 3 (HEV3). PCMV/PRV has been shown to reduce the survival time of pig transplants in non-human primates and was also transmitted in the first pig heart transplantation to a human patient. The main aim of this study was to determine the sensitivities of our methods to detect PCMV/PRV, HEV3, porcine lymphotropic herpesvirus-1 (PLHV-1), PLHV-2, PLHV-3, porcine circovirus 2 (PCV2), PCV3, PCV4 and porcine parvovirus 1 (PPV1) and to apply the methods to screen indigenous Greek black pigs. The high number of viruses found in these animals allowed for the evaluation of numerous detection methods. Since porcine endogenous retroviruses (PERVs) type A and B are integrated in the genome of all pigs, but PERV-C is not, the animals were screened for PERV-C and PERV-A/C. Our detection methods were sensitive and detected PCMV/PRV, PLHV-1, PLHV-1, PLHV-3, PVC3 and PERV-C in most animals. PPV1, HEV3, PCV4 and PERV-A/C were not detected. These data are of great interest since the animals are healthy and resistant to diseases.

Engineered deaminases as a key component of DNA and RNA editing tools

Over recent years, zinc-dependent deaminases have attracted increasing interest as key components of nucleic acid editing tools that can generate point mutations at specific sites in either DNA or RNA by combining a targeting module (such as a catalytically impaired CRISPR-Cas component) and an effector module (most often a deaminase). Deaminase-based molecular tools are already being utilized in a wide spectrum of therapeutic and research applications; however, their medical and biotechnological potential seems to be much greater. Recent reports indicate that the further development of nucleic acid editing systems depends largely on our ability to engineer the substrate specificity and catalytic activity of the editors themselves. In this review, we summarize the current trends and achievements in deaminase engineering. The presented data indicate that the potential of these enzymes has not yet been fully revealed or understood. Several examples show that even relatively minor changes in the structure of deaminases can give them completely new and unique properties.

Evolvability of cancer-associated genes under APOBEC3A/B selection

Evolvability is an emergent hallmark of cancer that depends on intra-tumor heterogeneity and, ultimately, genetic variation. Mutations generated by APOBEC3 cytidine deaminases can contribute to genetic variation and the consequences of APOBEC activation differ depending on the stage of cancer, with the most significant impact observed during the early stages. However, how APOBEC activity shapes evolutionary patterns of genes in the host genome and differential impacts on cancer-associated and non-cancer genes remain unclear. Analyzing over 40,000 human protein-coding transcripts, we identified distinct distribution patterns of APOBEC3A/B TC motifs between cancer-related genes and controls, suggesting unique associations with cancer. Studying a bat species with many more APOBEC3 genes, we found diverse motif patterns in orthologs of cancer genes compared to controls, similar to humans and suggesting APOBEC evolution to reduce impacts on the genome rather than the converse. Simulations confirmed that APOBEC-induced heterogeneity enhances cancer evolution, shaping clonal dynamics through bimodal introduction of mutations in certain classes of genes. Our results suggest that a major consequence of the bimodal distribution of APOBEC affects greater cancer heterogeneity.

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.

Reservoir host immunology and life history shape virulence evolution in zoonotic viruses

The management of future pandemic risk requires a better understanding of the mechanisms that determine the virulence of emerging zoonotic viruses. Meta-analyses suggest that the virulence of emerging zoonoses is correlated with but not completely predictable from reservoir host phylogeny, indicating that specific characteristics of reservoir host immunology and life history may drive the evolution of viral traits responsible for cross-species virulence. In particular, bats host viruses that cause higher case fatality rates upon spillover to humans than those derived from any other mammal, a phenomenon that cannot be explained by phylogenetic distance alone. In order to disentangle the fundamental drivers of these patterns, we develop a nested modeling framework that highlights mechanisms that underpin the evolution of viral traits in reservoir hosts that cause virulence following cross-species emergence. We apply this framework to generate virulence predictions for viral zoonoses derived from diverse mammalian reservoirs, recapturing trends in virus-induced human mortality rates reported in the literature. Notably, our work offers a mechanistic hypothesis to explain the extreme virulence of bat-borne zoonoses and, more generally, demonstrates how key differences in reservoir host longevity, viral tolerance, and constitutive immunity impact the evolution of viral traits that cause virulence following spillover to humans. Our theoretical framework offers a series of testable questions and predictions designed to stimulate future work comparing cross-species virulence evolution in zoonotic viruses derived from diverse mammalian hosts.

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.

Promoter evolution of mammalian gene duplicates

BACKGROUND

Gene duplication is thought to be a central process in evolution to gain new functions. The factors that dictate gene retention following duplication as well paralog gene divergence in sequence, expression and function have been extensively studied. However, relatively little is known about the evolution of promoter regions of gene duplicates and how they influence gene duplicate divergence. Here, we focus on promoters of paralog genes, comparing their similarity in sequence, in the sets of transcription factors (TFs) that bind them, and in their overall promoter architecture.

RESULTS

We observe that promoters of recent duplications display higher sequence similarity between them and that sequence similarity rapidly declines between promoters of more ancient paralogs. In contrast, similarity in cis-regulation, as measured by the set of TFs that bind promoters of both paralogs, does not simply decrease with time from duplication and is instead related to promoter architecture-paralogs with CpG Islands (CGIs) in their promoters share a greater fraction of TFs, while CGI-less paralogs are more divergent in their TF binding set. Focusing on recent duplication events and partitioning them by their duplication mechanism enables us to uncover promoter properties associated with gene retention, as well as to characterize the evolution of promoters of newly born genes: In recent retrotransposition-mediated duplications, we observe asymmetry in cis-regulation of paralog pairs: Retrocopy genes are lowly expressed and their promoters are bound by fewer TFs and are depleted of CGIs, in comparison with the original gene copy. Furthermore, looking at recent segmental duplication regions in primates enable us to compare successful retentions versus loss of duplicates, showing that duplicate retention is associated with fewer TFs and with CGI-less promoter architecture.

CONCLUSIONS

In this work, we profiled promoters of gene duplicates and their inter-paralog divergence. We also studied how their characteristics are associated with duplication time and duplication mechanism, as well as with the fate of these duplicates. These results underline the importance of cis-regulatory mechanisms in shaping the evolution of new genes and their fate following duplication.

Characterization of Megabat-Favored, CA-Dependent Susceptibility to Retrovirus Infection

The isolation of the Koala retrovirus-like virus from Australian megabats and the identification of endogenous retroviruses in the bat genome have raised questions on bat susceptibility to retroviruses in general. To answer this, we studied the susceptibility of 12 cell lines from 11 bat species to four well-studied retroviruses (human and simian immunodeficiency viruses [HIV and SIV] and murine leukemia viruses [B- and N-MLV]). Systematic comparison of retroviral susceptibility among bats revealed that megabat cell lines were overall less susceptible to the four retroviruses than microbat cell lines, particularly to HIV-1 infection, whereas lineage-specific differences were observed for MLV susceptibility. Quantitative PCR of reverse transcription (RT) products, infection in heterokaryon cells, and point mutation analysis of the capsid (CA) revealed that (i) HIV-1 and MLV replication were blocked at the nuclear transport of the pre-integration complexes and before and/or during RT, respectively, and (ii) the observed lineage-specific restriction can be attributed to a dominant cellular factor constrained by specific positions in CA. Investigation of bat homologs of the three previously reported post-entry restriction factors constrained by the same residues in CA, tripartite motif-protein 5α (TRIM5α), myxovirus resistance 2/B (Mx2/MxB), and carboxy terminus-truncated cleavage and polyadenylation factor 6 (CPSF6-358), demonstrated poor anti-HIV-1 activity in megabat cells, whereas megabat TRIM5α restricted MLV infection, suggesting that the major known CA-dependent restriction factors were not dominant in the observed lineage-specific susceptibility to HIV-1 in bat cells. Therefore, HIV-1 susceptibility of megabat cells may be determined in a manner distinct from that of primate cells. IMPORTANCE Recent studies have demonstrated the circulation of gammaretroviruses among megabats in Australia and the bats' resistance to HIV-1 infection; however, the origins of these viruses in megabats and the contribution of bats to retrovirus spread to other mammalian species remains unclear. To determine the intrinsic susceptibility of bat cells to HIV-1 infection, we investigated 12 cell lines isolated from 11 bat species. We report that lineage-specific retrovirus restriction in the bat cell lines can be attributed to CA-dependent factors. However, in the megabat cell lines examined, factors known to bind capsid and block infection in primate cell culture, including homologs of TRIM5α, Mx2/MxB, and CPSF6, failed to exhibit significant anti-HIV-1 activities. These results suggested that the HIV-1 susceptibility of megabat cells occurs in a manner distinct from that of primate cells, where cellular factors, other than major known CA-dependent restriction factors, with lineage-specific functions could recognize retroviral proteins in megabats.

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.

Adaptive Evolution of the OAS Gene Family Provides New Insights into the Antiviral Ability of Laurasiatherian Mammals

Many mammals risk damage from virus invasion due to frequent environmental changes. The oligoadenylate synthesis (OAS) gene family, which is an important component of the immune system, provides an essential response to the antiviral activities of interferons by regulating immune signal pathways. However, little is known about the evolutionary characteristics of OASs in Laurasiatherian mammals. Here, we examined the evolution of the OAS genes in 64 mammals to explore the accompanying molecular mechanisms of the antiviral ability of Laurasiatherian mammals living in different environments. We found that OAS2 and OAS3 were found to be pseudogenes in Odontoceti species. This may be related to the fact that they live in water. Some Antilopinae, Caprinae, and Cervidae species lacked the OASL gene, which may be related to their habitats being at higher altitudes. The OASs had a high number of positive selection sites in Cetartiodactyla, which drove the expression of strong antiviral ability. The OAS gene family evolved in Laurasiatherian mammals at different rates and was highly correlated with the species' antiviral ability. The gene evolution rate in Cetartiodactyla was significantly higher than that in the other orders. Compared to other species of the Carnivora family, the higher selection pressure on the OAS gene and the absence of positive selection sites in Canidae may be responsible for its weak resistance to rabies virus. The OAS gene family was relatively conserved during evolution. Conserved genes are able to provide better maintenance of gene function. The rate of gene evolution and the number of positively selected sites combine to influence the resistance of a species to viruses. The positive selection sites demonstrate the adaptive evolution of the OAS gene family to the environment. Adaptive evolution combined with conserved gene function improves resistance to viruses. Our findings offer insights into the molecular and functional evolution of the antiviral ability of Laurasian mammals.

Evidence linking APOBEC3B genesis and evolution of innate immune antagonism by gamma-herpesvirus ribonucleotide reductases

Viruses have evolved diverse mechanisms to antagonize host immunity such as direct inhibition and relocalization of cellular APOBEC3B (A3B) by the ribonucleotide reductase (RNR) of Epstein-Barr virus. Here, we investigate the mechanistic conservation and evolutionary origin of this innate immune counteraction strategy. First, we find that human gamma-herpesvirus RNRs engage A3B via largely distinct surfaces. Second, we show that RNR-mediated enzymatic inhibition and relocalization of A3B depend upon binding to different regions of the catalytic domain. Third, we show that the capability of viral RNRs to antagonize A3B is conserved among gamma-herpesviruses that infect humans and Old World monkeys that encode this enzyme but absent in homologous viruses that infect New World monkeys that naturally lack the A3B gene. Finally, we reconstruct the ancestral primate A3B protein and demonstrate that it is active and similarly engaged by the RNRs from viruses that infect humans and Old World monkeys but not by the RNRs from viruses that infect New World monkeys. These results combine to indicate that the birth of A3B at a critical branchpoint in primate evolution may have been a driving force in selecting for an ancestral gamma-herpesvirus with an expanded RNR functionality through counteraction of this antiviral enzyme.

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.

Functions and consequences of AID/APOBEC-mediated DNA and RNA deamination

The AID/APOBEC polynucleotide cytidine deaminases have historically been classified as either DNA mutators or RNA editors based on their first identified nucleic acid substrate preference. DNA mutators can generate functional diversity at antibody genes but also cause genomic instability in cancer. RNA editors can generate informational diversity in the transcriptome of innate immune cells, and of cancer cells. Members of both classes can act as antiviral restriction factors. Recent structural work has illuminated differences and similarities between AID/APOBEC enzymes that can catalyse DNA mutation, RNA editing or both, suggesting that the strict functional classification of members of this family should be reconsidered. As many of these enzymes have been employed for targeted genome (or transcriptome) editing, a more holistic understanding will help improve the design of therapeutically relevant programmable base editors.

In this Perspective, Pecori et al. provide an overview of the AID/APOBEC cytidine deaminase family, discussing key structural features, how they contribute to viral and tumour evolution and how they can be harnessed for (potentially therapeutic) base-editing purposes.

Retroviruses of Bats: a Threat Waiting in the Wings?

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

Factors Regulating the Activity of LINE1 Retrotransposons

LINE-1 (L1) is a class of autonomous mobile genetic elements that form somatic mosaicisms in various tissues of the organism. The activity of L1 retrotransposons is strictly controlled by many factors in somatic and germ cells at all stages of ontogenesis. Alteration of L1 activity was noted in a number of diseases: in neuropsychiatric and autoimmune diseases, as well as in various forms of cancer. Altered activity of L1 retrotransposons for some pathologies is associated with epigenetic changes and defects in the genes involved in their repression. This review discusses the molecular genetic mechanisms of the retrotransposition and regulation of the activity of L1 elements. The contribution of various factors controlling the expression and distribution of L1 elements in the genome occurs at all stages of the retrotransposition. The regulation of L1 elements at the transcriptional, post-transcriptional and integration into the genome stages is described in detail. Finally, this review also focuses on the evolutionary aspects of L1 accumulation and their interplay with the host regulation system.

Variability in Codon Usage in Coronaviruses Is Mainly Driven by Mutational Bias and Selective Constraints on CpG Dinucleotide

The Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the third human-emerged virus of the 21st century from the Coronaviridae family, causing the ongoing coronavirus disease 2019 (COVID-19) pandemic. Due to the high zoonotic potential of coronaviruses, it is critical to unravel their evolutionary history of host species breadth, host-switch potential, adaptation and emergence, to identify viruses posing a pandemic risk in humans. We present here a comprehensive analysis of the composition and codon usage bias of the 82 Orthocoronavirinae members, infecting 47 different avian and mammalian hosts. Our results clearly establish that synonymous codon usage varies widely among viruses, is only weakly dependent on their primary host, and is dominated by mutational bias towards AU-enrichment and by CpG avoidance. Indeed, variation in GC3 explains around 34%, while variation in CpG frequency explains around 14% of total variation in codon usage bias. Further insight on the mutational equilibrium within Orthocoronavirinae revealed that most coronavirus genomes are close to their neutral equilibrium, the exception being the three recently infecting human coronaviruses, which lie further away from the mutational equilibrium than their endemic human coronavirus counterparts. Finally, our results suggest that, while replicating in humans, SARS-CoV-2 is slowly becoming AU-richer, likely until attaining a new mutational equilibrium.

Large-scale genome sampling reveals unique immunity and metabolic adaptations in bats

Comprising more than 1,400 species, bats possess adaptations unique among mammals including powered flight, unexpected longevity, and extraordinary immunity. Some of the molecular mechanisms underlying these unique adaptations includes DNA repair, metabolism and immunity. However, analyses have been limited to a few divergent lineages, reducing the scope of inferences on gene family evolution across the Order Chiroptera. We conducted an exhaustive comparative genomic study of 37 bat species, one generated in this study, encompassing a large number of lineages, with a particular emphasis on multi-gene family evolution across immune and metabolic genes. In agreement with previous analyses, we found lineage-specific expansions of the APOBEC3 and MHC-I gene families, and loss of the proinflammatory PYHIN gene family. We inferred more than 1,000 gene losses unique to bats, including genes involved in the regulation of inflammasome pathways such as epithelial defence receptors, the natural killer gene complex and the interferon-gamma induced pathway. Gene set enrichment analyses revealed genes lost in bats are involved in defence response against pathogen-associated molecular patterns and damage-associated molecular patterns. Gene family evolution and selection analyses indicate bats have evolved fundamental functional differences compared to other mammals in both innate and adaptive immune system, with the potential to enhance antiviral immune response while dampening inflammatory signalling. In addition, metabolic genes have experienced repeated expansions related to convergent shifts to plant-based diets. Our analyses support the hypothesis that, in tandem with flight, ancestral bats had evolved a unique set of immune adaptations whose functional implications remain to be explored.

Examination of the APOBEC3 Barrier to Cross Species Transmission of Primate Lentiviruses

The transmission of viruses from animal hosts into humans have led to the emergence of several diseases. Usually these cross-species transmissions are blocked by host restriction factors, which are proteins that can block virus replication at a specific step. In the natural virus host, the restriction factor activity is usually suppressed by a viral antagonist protein, but this is not the case for restriction factors from an unnatural host. However, due to ongoing viral evolution, sometimes the viral antagonist can evolve to suppress restriction factors in a new host, enabling cross-species transmission. Here we examine the classical case of this paradigm by reviewing research on APOBEC3 restriction factors and how they can suppress human immunodeficiency virus (HIV) and simian immunodeficiency virus (SIV). APOBEC3 enzymes are single-stranded DNA cytidine deaminases that can induce mutagenesis of proviral DNA by catalyzing the conversion of cytidine to promutagenic uridine on single-stranded viral (-)DNA if they escape the HIV/SIV antagonist protein, Vif. APOBEC3 degradation is induced by Vif through the proteasome pathway. SIV has been transmitted between Old World Monkeys and to hominids. Here we examine the adaptations that enabled such events and the ongoing impact of the APOBEC3-Vif interface on HIV in humans.

Foamy Viruses, Bet, and APOBEC3 Restriction

Non-human primates (NHP) are an important source of viruses that can spillover to humans and, after adaptation, spread through the host population. Whereas HIV-1 and HTLV-1 emerged as retroviral pathogens in humans, a unique class of retroviruses called foamy viruses (FV) with zoonotic potential are occasionally detected in bushmeat hunters or zookeepers. Various FVs are endemic in numerous mammalian natural hosts, such as primates, felines, bovines, and equines, and other animals, but not in humans. They are apathogenic, and significant differences exist between the viral life cycles of FV and other retroviruses. Importantly, FVs replicate in the presence of many well-defined retroviral restriction factors such as TRIM5α, BST2 (Tetherin), MX2, and APOBEC3 (A3). While the interaction of A3s with HIV-1 is well studied, the escape mechanisms of FVs from restriction by A3 is much less explored. Here we review the current knowledge of FV biology, host restriction factors, and FV-host interactions with an emphasis on the consequences of FV regulatory protein Bet binding to A3s and outline crucial open questions for future studies.

The Battle between Retroviruses and APOBEC3 Genes: Its Past and Present

The APOBEC3 family of proteins in mammals consists of cellular cytosine deaminases and well-known restriction factors against retroviruses, including lentiviruses. APOBEC3 genes are highly amplified and diversified in mammals, suggesting that their evolution and diversification have been driven by conflicts with ancient viruses. At present, lentiviruses, including HIV, the causative agent of AIDS, are known to encode a viral protein called Vif to overcome the antiviral effects of the APOBEC3 proteins of their hosts. Recent studies have revealed that the acquisition of an anti-APOBEC3 ability by lentiviruses is a key step in achieving successful cross-species transmission. Here, we summarize the current knowledge of the interplay between mammalian APOBEC3 proteins and viral infections and introduce a scenario of the coevolution of mammalian APOBEC3 genes and viruses.

Potential APOBEC-mediated RNA editing of the genomes of SARS-CoV-2 and other coronaviruses and its impact on their longer term evolution

Members of the APOBEC family of cytidine deaminases show antiviral activities in mammalian cells through lethal editing in the genomes of small DNA viruses, herpesviruses and retroviruses, and potentially those of RNA viruses such as coronaviruses. Consistent with the latter, APOBEC-like directional C→U transitions of genomic plus-strand RNA are greatly overrepresented in SARS-CoV-2 genome sequences of variants emerging during the COVID-19 pandemic. A C→U mutational process may leave evolutionary imprints on coronavirus genomes, including extensive homoplasy from editing and reversion at targeted sites and the occurrence of driven amino acid sequence changes in viral proteins. If sustained over longer periods, this process may account for the previously reported marked global depletion of C and excess of U bases in human seasonal coronavirus genomes. This review synthesizes the current knowledge on APOBEC evolution and function and the evidence of their role in APOBEC-mediated genome editing of SARS-CoV-2 and other coronaviruses.

-

SARS-CoV-2 sequence variants contain an overabundance of C- > U transitions

-

C- > U transitions are the hallmark of the activity of APOBEC cytosine deaminases

-

Further work is needed to determine APOBEC's role in coronavirus evolution

Rapid Evolution of HERC6 and Duplication of a Chimeric HERC5/6 Gene in Rodents and Bats Suggest an Overlooked Role of HERCs in Mammalian Immunity

Studying the evolutionary diversification of mammalian antiviral defenses is of main importance to better understand our innate immune repertoire. The small HERC proteins are part of a multigene family, including HERC5 and HERC6, which have probably diversified through complex evolutionary history in mammals. Here, we performed mammalian-wide phylogenetic and genomic analyses of HERC5 and HERC6, using 83 orthologous sequences from bats, rodents, primates, artiodactyls, and carnivores—the top five representative groups of mammalian evolution. We found that HERC5 has been under weak and differential positive selection in mammals, with only primate HERC5 showing evidences of pathogen-driven selection. In contrast, HERC6 has been under strong and recurrent adaptive evolution in mammals, suggesting past and widespread genetic arms-races with viral pathogens. Importantly, the rapid evolution of mammalian HERC6 spacer domain suggests that it might be a host-pathogen interface, targeting viral proteins and/or being the target of virus antagonists. Finally, we identified a HERC5/6 chimeric gene that arose from independent duplication in rodent and bat lineages and encodes for a conserved HERC5 N-terminal domain and divergent HERC6 spacer and HECT domains. This duplicated chimeric gene highlights adaptations that potentially contribute to rodent and bat immunity. Our findings open new research avenues on the functions of HERC6 and HERC5/6 in mammals, and on their implication in antiviral innate immunity.

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 genomes carry the signatures of their habitats

Coronaviruses such as SARS-CoV-2 regularly infect host tissues that express antiviral proteins (AVPs) in abundance. Understanding how they evolve to adapt or evade host immune responses is important in the effort to control the spread of infection. Two AVPs that may shape viral genomes are the zinc finger antiviral protein (ZAP) and the apolipoprotein B mRNA editing enzyme-catalytic polypeptide-like 3 (APOBEC3). The former binds to CpG dinucleotides to facilitate the degradation of viral transcripts while the latter frequently deaminates C into U residues which could generate notable viral sequence variations. We tested the hypothesis that both APOBEC3 and ZAP impose selective pressures that shape the genome of an infecting coronavirus. Our investigation considered a comprehensive number of publicly available genomes for seven coronaviruses (SARS-CoV-2, SARS-CoV, and MERS infecting Homo sapiens, Bovine CoV infecting Bos taurus, MHV infecting Mus musculus, HEV infecting Sus scrofa, and CRCoV infecting Canis lupus familiaris). We show that coronaviruses that regularly infect tissues with abundant AVPs have CpG-deficient and U-rich genomes; whereas those that do not infect tissues with abundant AVPs do not share these sequence hallmarks. Among the coronaviruses surveyed herein, CpG is most deficient in SARS-CoV-2 and a temporal analysis showed a marked increase in C to U mutations over four months of SARS-CoV-2 genome evolution. Furthermore, the preferred motifs in which these C to U mutations occur are the same as those subjected to APOBEC3 editing in HIV-1. These results suggest that both ZAP and APOBEC3 shape the SARS-CoV-2 genome: ZAP imposes a strong CpG avoidance, and APOBEC3 constantly edits C to U. Evolutionary pressures exerted by host immune systems onto viral genomes may motivate novel strategies for SARS-CoV-2 vaccine development.

Deaminase-Independent Mode of Antiretroviral Action in Human and Mouse APOBEC3 Proteins

Apolipoprotein B mRNA editing enzyme, catalytic polypeptide-like 3 (APOBEC3) proteins (APOBEC3s) are deaminases that convert cytosines to uracils predominantly on a single-stranded DNA, and function as intrinsic restriction factors in the innate immune system to suppress replication of viruses (including retroviruses) and movement of retrotransposons. Enzymatic activity is supposed to be essential for the APOBEC3 antiviral function. However, it is not the only way that APOBEC3s exert their biological function. Since the discovery of human APOBEC3G as a restriction factor for HIV-1, the deaminase-independent mode of action has been observed. At present, it is apparent that both the deaminase-dependent and -independent pathways are tightly involved not only in combating viruses but also in human tumorigenesis. Although the deaminase-dependent pathway has been extensively characterized so far, understanding of the deaminase-independent pathway remains immature. Here, we review existing knowledge regarding the deaminase-independent antiretroviral functions of APOBEC3s and their molecular mechanisms. We also discuss the possible unidentified molecular mechanism for the deaminase-independent antiretroviral function mediated by mouse APOBEC3.

A Potent Postentry Restriction to Primate Lentiviruses in a Yinpterochiropteran Bat

The COVID-19 pandemic suggests that bat innate immune systems are insufficiently characterized relative to the medical importance of these animals. Retroviruses, e.g., HIV-1, can be severe pathogens when they cross species barriers, and bat restrictions corresponding to retroviruses are comparatively unstudied. Here, we compared the abilities of retroviruses from three genera (Lentivirus, Gammaretrovirus, and Spumavirus) to infect cells of the large fruit-eating bat P. alecto and other mammals. We identified a major, specific postentry restriction to primate lentiviruses. HIV-1 and SIVmac are potently blocked at early life cycle steps, but nonprimate lentiviruses and foamy retroviruses are entirely unrestricted. Despite acting postentry and in a CypA-dependent manner with features reminiscent of antiretroviral factors from other mammals, this restriction was not saturable with virus-like particles and was independent of P. alecto TRIM5, TRIM21, TRIM22, TRIM34, and MX2. These results identify a novel restriction and highlight cyclophilin-capsid interactions as ancient species-specific determinants of retroviral infection.

Bats are primary reservoirs for multiple lethal human viruses, such as Ebola, Nipah, Hendra, rabies, severe acute respiratory syndrome coronavirus (SARS-CoV), Middle East respiratory syndrome-related coronavirus (MERS-CoV), and, most recently, SARS-CoV-2. The innate immune systems of these immensely abundant, anciently diverged mammals remain insufficiently characterized. While bat genomes contain many endogenous retroviral elements indicative of past exogenous infections, little is known about restrictions to extant retroviruses. Here, we describe a major postentry restriction in cells of the yinpterochiropteran bat Pteropus alecto. Primate lentiviruses (HIV-1, SIVmac) were potently blocked at early life cycle steps, with up to 1,000-fold decreases in infectivity. The block was specific, because nonprimate lentiviruses such as equine infectious anemia virus and feline immunodeficiency virus were unimpaired, as were foamy retroviruses. Interspecies heterokaryons demonstrated a dominant block consistent with restriction of incoming viruses. Several features suggested potential TRIM5 (tripartite motif 5) or myxovirus resistance protein 2 (MX2) protein restriction, including postentry action, cyclosporine sensitivity, and reversal by capsid cyclophilin A (CypA) binding loop mutations. Viral nuclear import was significantly reduced, and this deficit was substantially rescued by cyclosporine treatment. However, saturation with HIV-1 virus-like particles did not relieve the restriction at all. P. alecto TRIM5 was inactive against HIV-1 although it blocked the gammaretrovirus N-tropic murine leukemia virus. Despite major divergence in a critical N-terminal motif required for human MX2 activity, P. alecto MX2 had anti-HIV activity. However, this did not quantitatively account for the restriction and was independent of and synergistic with an additional CypA-dependent restriction. These results reveal a novel, specific restriction to primate lentiviruses in the Pteropodidae and advance understanding of bat innate immunity.

Species-Specific Evolution of Ebola Virus during Replication in Human and Bat Cells

Ebola virus (EBOV) causes a severe, often fatal disease in humans and nonhuman primates. Within the past decade, EBOV has caused two large and difficult-to-control outbreaks, one of which recently ended in the Democratic Republic of the Congo. Bats are the likely reservoir of EBOV, but little is known of their relationship with the virus. We perform serial passages of EBOV in human and bat cells and use circular sequencing to compare the short-term evolution of the virus. Virus populations passaged in bat cells have sequence markers indicative of host RNA editing enzyme activity, including evidence for ADAR editing of the EBOV glycoprotein. Multiple regions in the EBOV genome appear to have undergone adaptive evolution when passaged in bat and human cells. Individual mutated viruses are rescued and characterized. Our results provide insight into the host species-specific evolution of EBOV and highlight the adaptive flexibility of the virus.

Footprint of the host restriction factors APOBEC3 on the genome of human viruses

APOBEC3 enzymes are innate immune effectors that introduce mutations into viral genomes. These enzymes are cytidine deaminases which transform cytosine into uracil. They preferentially mutate cytidine preceded by thymidine making the 5'TC motif their favored target. Viruses have evolved different strategies to evade APOBEC3 restriction. Certain viruses actively encode viral proteins antagonizing the APOBEC3s, others passively face the APOBEC3 selection pressure thanks to a depleted genome for APOBEC3-targeted motifs. Hence, the APOBEC3s left on the genome of certain viruses an evolutionary footprint. The aim of our study is the identification of these viruses having a genome shaped by the APOBEC3s. We analyzed the genome of 33,400 human viruses for the depletion of APOBEC3-favored motifs. We demonstrate that the APOBEC3 selection pressure impacts at least 22% of all currently annotated human viral species. The papillomaviridae and polyomaviridae are the most intensively footprinted families; evidencing a selection pressure acting genome-wide and on both strands. Members of the parvoviridae family are differentially targeted in term of both magnitude and localization of the footprint. Interestingly, a massive APOBEC3 footprint is present on both strands of the B19 erythroparvovirus; making this viral genome one of the most cleaned sequences for APOBEC3-favored motifs. We also identified the endemic coronaviridae as significantly footprinted. Interestingly, no such footprint has been detected on the zoonotic MERS-CoV, SARS-CoV-1 and SARS-CoV-2 coronaviruses. In addition to viruses that are footprinted genome-wide, certain viruses are footprinted only on very short sections of their genome. That is the case for the gamma-herpesviridae and adenoviridae where the footprint is localized on the lytic origins of replication. A mild footprint can also be detected on the negative strand of the reverse transcribing HIV-1, HIV-2, HTLV-1 and HBV viruses. Together, our data illustrate the extent of the APOBEC3 selection pressure on the human viruses and identify new putatively APOBEC3-targeted viruses.

Six reference-quality genomes reveal evolution of bat adaptations

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

Retrocopying expands the functional repertoire of APOBEC3 antiviral proteins in primates

Host-virus arms races are inherently asymmetric; viruses evolve much more rapidly than host genomes. Thus, there is high interest in discovering mechanisms by which host genomes keep pace with rapidly evolving viruses. One family of restriction factors, the APOBEC3 (A3) cytidine deaminases, has undergone positive selection and expansion via segmental gene duplication and recombination. Here, we show that new copies of A3 genes have also been created in primates by reverse transcriptase-encoding elements like LINE-1 or endogenous retroviruses via a process termed retrocopying. First, we discovered that all simian primate genomes retain the remnants of an ancient A3 retrocopy: A3I. Furthermore, we found that some New World monkeys encode up to ten additional APOBEC3G (A3G) retrocopies. Some of these A3G retrocopies are transcribed in a variety of tissues and able to restrict retroviruses. Our findings suggest that host genomes co-opt retroelement activity in the germline to create new host restriction factors as another means to keep pace with the rapid evolution of viruses. (163).

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

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

Machine learning using intrinsic genomic signatures for rapid classification of novel pathogens: COVID-19 case study

The 2019 novel coronavirus (renamed SARS-CoV-2, and generally referred to as the COVID-19 virus) has spread to 184 countries with over 1.5 million confirmed cases. Such major viral outbreaks demand early elucidation of taxonomic classification and origin of the virus genomic sequence, for strategic planning, containment, and treatment. This paper identifies an intrinsic COVID-19 virus genomic signature and uses it together with a machine learning-based alignment-free approach for an ultra-fast, scalable, and highly accurate classification of whole COVID-19 virus genomes. The proposed method combines supervised machine learning with digital signal processing (MLDSP) for genome analyses, augmented by a decision tree approach to the machine learning component, and a Spearman's rank correlation coefficient analysis for result validation. These tools are used to analyze a large dataset of over 5000 unique viral genomic sequences, totalling 61.8 million bp, including the 29 COVID-19 virus sequences available on January 27, 2020. Our results support a hypothesis of a bat origin and classify the COVID-19 virus as Sarbecovirus, within Betacoronavirus. Our method achieves 100% accurate classification of the COVID-19 virus sequences, and discovers the most relevant relationships among over 5000 viral genomes within a few minutes, ab initio, using raw DNA sequence data alone, and without any specialized biological knowledge, training, gene or genome annotations. This suggests that, for novel viral and pathogen genome sequences, this alignment-free whole-genome machine-learning approach can provide a reliable real-time option for taxonomic classification.

Infectious KoRV-related retroviruses circulating in Australian bats

Bats are reservoirs of emerging viruses that are highly pathogenic to other mammals, including humans. Despite the diversity and abundance of bat viruses, to date they have not been shown to harbor exogenous retroviruses. Here we report the discovery and characterization of a group of koala retrovirus-related (KoRV-related) gammaretroviruses in Australian and Asian bats. These include the Hervey pteropid gammaretrovirus (HPG), identified in the scat of the Australian black flying fox (Pteropus alecto), which is the first reproduction-competent retrovirus found in bats. HPG is a close relative of KoRV and the gibbon ape leukemia virus (GALV), with virion morphology and Mn²⁺-dependent virion-associated reverse transcriptase activity typical of a gammaretrovirus. In vitro, HPG is capable of infecting bat and human cells, but not mouse cells, and displays a similar pattern of cell tropism as KoRV-A and GALV. Population studies reveal the presence of HPG and KoRV-related sequences in several locations across northeast Australia, as well as serologic evidence for HPG in multiple pteropid bat species, while phylogenetic analysis places these bat viruses as the basal group within the KoRV-related retroviruses. Taken together, these results reveal bats to be important reservoirs of exogenous KoRV-related gammaretroviruses.

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

This work represents a large, order-wide evolutionary analysis of the order Chiroptera (bats). Our pipeline for assembling sequence data and curating orthologous multiple sequence alignments includes methods for improving results when combining genomic and transcriptomic data sources. The resulting phylogenetic tree divides the order Chiroptera into Yinpterochiroptera and Yangochiroptera, in disagreement with the previous division into Megachiroptera and Microchiroptera and in agreement with some other recent molecular studies, and also provides evidence for other contested branch placements. We also performed a genome-wide analysis of positive selection and found 181 genes with signatures of positive selection. Enrichment analysis shows these positively selected genes to be primarily related to immune responses but also, surprisingly, collagen formation.

Historically, the evolution of bats has been analyzed using a small number of genetic loci for many species or many genetic loci for a few species. Here we present a phylogeny of 18 bat species, each of which is represented in 1,107 orthologous gene alignments used to build the tree. We generated a transcriptome sequence of Hypsignathus monstrosus, the African hammer-headed bat, and additional transcriptome sequence for Rousettus aegyptiacus, the Egyptian fruit bat. We then combined these data with existing genomic and transcriptomic data from 16 other bat species. In the analysis of such datasets, there is no clear consensus on the most reliable computational methods for the curation of quality multiple sequence alignments since these public datasets represent multiple investigators and methods, including different source materials (chromosomal DNA or expressed RNA). Here we lay out a systematic analysis of parameters and produce an advanced pipeline for curating orthologous gene alignments from combined transcriptomic and genomic data, including a software package: the Mismatching Isoform eXon Remover (MIXR). Using this method, we created alignments of 11,677 bat genes, 1,107 of which contain orthologs from all 18 species. Using the orthologous gene alignments created, we assessed bat phylogeny and also performed a holistic analysis of positive selection acting in bat genomes. We found that 181 genes have been subject to positive natural selection. This list is dominated by genes involved in immune responses and genes involved in the production of collagens.

Genetic and mechanistic basis for APOBEC3H alternative splicing, retrovirus restriction, and counteraction by HIV-1 protease

Human APOBEC3H (A3H) is a single-stranded DNA cytosine deaminase that inhibits HIV-1. Seven haplotypes (I–VII) and four splice variants (SV154/182/183/200) with differing antiviral activities and geographic distributions have been described, but the genetic and mechanistic basis for variant expression and function remains unclear. Using a combined bioinformatic/experimental analysis, we find that SV200 expression is specific to haplotype II, which is primarily found in sub-Saharan Africa. The underlying genetic mechanism for differential mRNA splicing is an ancient intronic deletion [del(ctc)] within A3H haplotype II sequence. We show that SV200 is at least fourfold more HIV-1 restrictive than other A3H splice variants. To counteract this elevated antiviral activity, HIV-1 protease cleaves SV200 into a shorter, less restrictive isoform. Our analyses indicate that, in addition to Vif-mediated degradation, HIV-1 may use protease as a  counter-defense mechanism against A3H in >80% of sub-Saharan African populations.

Human APOBEC3H has several haplotypes and splice variants with distinct anti-HIV-1 activities, but the genetics underlying the expression of these variants are unclear. Here, the authors identify an intronic deletion in A3H haplotype II resulting in production of the most active splice variant, which is counteracted by HIV-1 protease.

Going to Bat(s) for Studies of Disease Tolerance

A majority of viruses that have caused recent epidemics with high lethality rates in people, are zoonoses originating from wildlife. Among them are filoviruses (e.g., Marburg, Ebola), coronaviruses (e.g., SARS, MERS), henipaviruses (e.g., Hendra, Nipah) which share the common features that they are all RNA viruses, and that a dysregulated immune response is an important contributor to the tissue damage and hence pathogenicity that results from infection in humans. Intriguingly, these viruses also all originate from bat reservoirs. Bats have been shown to have a greater mean viral richness than predicted by their phylogenetic distance from humans, their geographic range, or their presence in urban areas, suggesting other traits must explain why bats harbor a greater number of zoonotic viruses than other mammals. Bats are highly unusual among mammals in other ways as well. Not only are they the only mammals capable of powered flight, they have extraordinarily long life spans, with little detectable increases in mortality or senescence until high ages. Their physiology likely impacted their history of pathogen exposure and necessitated adaptations that may have also affected immune signaling pathways. Do our life history traits make us susceptible to generating damaging immune responses to RNA viruses or does the physiology of bats make them particularly tolerant or resistant? Understanding what immune mechanisms enable bats to coexist with RNA viruses may provide critical fundamental insights into how to achieve greater resilience in humans.