The Egyptian Rousette Genome Reveals Unexpected Features of Bat Antiviral Immunity.

Bats harbor many viruses asymptomatically, including several notorious for causing extreme virulence in humans. To identify differences between antiviral mechanisms in humans and bats, we sequenced, assembled, and analyzed the genome of Rousettus aegyptiacus, a natural reservoir of Marburg virus and the only known reservoir for any filovirus. We found an expanded and diversified KLRC/KLRD family of natural killer cell receptors, MHC class I genes, and type I interferons, which dramatically differ from their functional counterparts in other mammals. Such concerted evolution of key components of bat immunity is strongly suggestive of novel modes of antiviral defense. An evaluation of the theoretical function of these genes suggests that an inhibitory immune state may exist in bats. Based on our findings, we hypothesize that tolerance of viral infection, rather than enhanced potency of antiviral defenses, may be a key mechanism by which bats asymptomatically host viruses that are pathogenic in humans.

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

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

Universal features of Nsp1-mediated translational shutdown by coronaviruses.

Nonstructural protein 1 (Nsp1) produced by coronaviruses inhibits host protein synthesis. The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Nsp1 C-terminal domain was shown to bind the ribosomal mRNA channel to inhibit translation, but it is unclear whether this mechanism is broadly used by coronaviruses, whether the Nsp1 N-terminal domain binds the ribosome, or how Nsp1 allows viral RNAs to be translated. Here, we investigated Nsp1 from SARS-CoV-2, Middle East respiratory syndrome coronavirus (MERS-CoV), and Bat-Hp-CoV coronaviruses using structural, biophysical, and biochemical experiments, revealing a conserved role for the C-terminal domain. Additionally, the N-terminal domain of Bat-Hp-CoV Nsp1 binds to the decoding center of the 40S subunit, where it would prevent mRNA and eIF1A accommodation. Structure-based experiments demonstrated the importance of decoding center interactions in all three coronaviruses and showed that the same regions of Nsp1 are necessary for the selective translation of viral RNAs. Our results provide a mechanistic framework to understand how Nsp1 controls preferential translation of viral RNAs.

The large bat Helitron DNA transposase forms a compact monomeric assembly that buries and protects its covalently bound 5'-transposon end.

Helitrons are widespread eukaryotic DNA transposons that have significantly contributed to genome variability and evolution, in part because of their distinctive, replicative rolling-circle mechanism, which often mobilizes adjacent genes. Although most eukaryotic transposases form oligomers and use RNase H-like domains to break and rejoin double-stranded DNA (dsDNA), Helitron transposases contain a single-stranded DNA (ssDNA)-specific HUH endonuclease domain. Here, we report the cryo-electron microscopy structure of a Helitron transposase bound to the 5'-transposon end, providing insight into its multidomain architecture and function. The monomeric transposase forms a tightly packed assembly that buries the covalently attached cleaved end, protecting it until the second end becomes available. The structure reveals unexpected architectural similarity to TraI, a bacterial relaxase that also catalyzes ssDNA movement. The HUH active site suggests how two juxtaposed tyrosines, a feature of many replication initiators that use HUH nucleases, couple the conformational shift of an α-helix to control strand cleavage and ligation reactions.

Pleiotropy of positive selection in ancient ACE2 suggests an alternative hypothesis for bat-specific adaptations to host coronaviruses.

Angiotensin-convertingenzyme 2 (ACE2) has dual functions, regulating cardiovascular physiology and serving as the receptor for coronaviruses. Bats, the only true flying mammals and natural viral reservoirs, have evolved positive alterations in traits related to both functions of ACE2. This suggests significant evolutionary changes in ACE2 during bat evolution. To test this hypothesis, we examine the selection pressure in ACE2 along the ancestral branch of all bats (AncBat-ACE2), where powered flight and bat-coronavirus coevolution occurred, and detect a positive selection signature. To assess the functional effects of positive selection, we resurrect AncBat-ACE2 and its mutant (AncBat-ACE2-mut) created by replacing the positively selected sites. Compared to AncBat-ACE2-mut, AncBat-ACE2 exhibits stronger enzymatic activity, enhances mice's performance in exercise fatigue, and shows lower affinity to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Our findings indicate the functional pleiotropy of positive selection in the ancient ACE2 of bats, providing an alternative hypothesis for the evolutionary origin of bats' defense against coronaviruses.

Unconventional IFNω-like Genes Dominate the Type I IFN Locus and the Constitutive Antiviral Responses in Bats.

Bats are the natural reservoir hosts of some viruses, some of which may spill over to humans and cause global-scale pandemics. Different from humans, bats may coexist with high pathogenic viruses without showing symptoms of diseases. As one of the most important first defenses, bat type I IFNs (IFN-Is) were thought to play a role during this virus coexistence and thus were studied in recent years. However, there are arguments about whether bats have a contracted genome locus or constitutively expressed IFNs, mainly due to species-specific findings. We hypothesized that because of the lack of pan-bat analysis, the common characteristics of bat IFN-Is have not been revealed yet. In this study, we characterized the IFN-I locus for nine Yangochiroptera bats and three Yinpterochiroptera bats on the basis of their high-quality bat genomes. We also compared the basal expression in six bats and compared the antiviral and antiproliferative activity and the thermostability of representative Rhinolophus bat IFNs. We found a dominance of unconventional IFNω-like responses in the IFN-I system, which is unique to bats. In contrast to IFNα-dominated IFN-I loci in the majority of other mammals, bats generally have shorter IFN-I loci with more unconventional IFNω-like genes (IFNω or related IFNαω), but with fewer or even no IFNα genes. In addition, bats generally have constitutively expressed IFNs, the highest expressed of which is more likely an IFNω-like gene. Likewise, the highly expressed IFNω-like protein also demonstrated the best antiviral activity, antiproliferative activity, or thermostability, as shown in a representative Rhinolophus bat species. Overall, we revealed pan-bat unique, to our knowledge, characteristics in the IFN-I system, which provide insights into our understanding of the innate immunity that contributes to a special coexistence between bats and 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.

Distinct Genes with Similar Functions Underlie Convergent Evolution in Myotis Bat Ecomorphs.

Convergence offers an opportunity to explore to what extent evolution can be predictable when genomic composition and environmental triggers are similar. Here, we present an emergent model system to study convergent evolution in nature in a mammalian group, the bat genus Myotis. Three foraging strategies-gleaning, trawling, and aerial hawking, each characterized by different sets of phenotypic features-have evolved independently multiple times in different biogeographic regions in isolation for millions of years. To investigate the genomic basis of convergence and explore the functional genomic changes linked to ecomorphological convergence, we sequenced and annotated 17 new genomes and screened 16,426 genes for positive selection and associations between relative evolutionary rates and foraging strategies across 30 bat species representing all Myotis ecomorphs across geographic regions as well as among sister groups. We identify genomic changes that describe both phylogenetic and ecomorphological trends. We infer that colonization of new environments may have first required changes in genes linked to hearing sensory perception, followed by changes linked to fecundity and development, metabolism of carbohydrates, and heme degradation. These changes may be linked to prey acquisition and digestion and match phylogenetic trends. Our findings also suggest that the repeated evolution of ecomorphs does not always involve changes in the same genes but rather in genes with the same molecular functions such as developmental and cellular processes.

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

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

MHC class II proteins mediate cross-species entry of bat influenza viruses.

Zoonotic influenza A viruses of avian origin can cause severe disease in individuals, or even global pandemics, and thus pose a threat to human populations. Waterfowl and shorebirds are believed to be the reservoir for all influenza A viruses, but this has recently been challenged by the identification of novel influenza A viruses in bats1,2. The major bat influenza A virus envelope glycoprotein, haemagglutinin, does not bind the canonical influenza A virus receptor, sialic acid or any other glycan1,3,4, despite its high sequence and structural homology with conventional haemagglutinins. This functionally uncharacterized plasticity of the bat influenza A virus haemagglutinin means the tropism and zoonotic potential of these viruses has not been fully determined. Here we show, using transcriptomic profiling of susceptible versus non-susceptible cells in combination with genome-wide CRISPR-Cas9 screening, that the major histocompatibility complex class II (MHC-II) human leukocyte antigen DR isotype (HLA-DR) is an essential entry determinant for bat influenza A viruses. Genetic ablation of the HLA-DR α-chain rendered cells resistant to infection by bat influenza A virus, whereas ectopic expression of the HLA-DR complex in non-susceptible cells conferred susceptibility. Expression of MHC-II from different bat species, pigs, mice or chickens also conferred susceptibility to infection. Notably, the infection of mice with bat influenza A virus resulted in robust virus replication in the upper respiratory tract, whereas mice deficient for MHC-II were resistant. Collectively, our data identify MHC-II as a crucial entry mediator for bat influenza A viruses in multiple species, which permits a broad vertebrate tropism.

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

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

Signals of positive selection in genomes of palearctic Myotis-bats coexisting with a fungal pathogen.

Disease can act as a driving force in shaping genetic makeup across populations, even species, if the impacts influence a particularly sensitive part of their life cycles. White-nose disease is caused by a fungal pathogen infecting bats during hibernation. The mycosis has caused massive population declines of susceptible species in North America, particularly in the genus Myotis. However, Myotis bats appear to tolerate infection in Eurasia, where the fungal pathogen has co-evolved with its bat hosts for an extended period of time. Therefore, with susceptible and tolerant populations, the fungal disease provides a unique opportunity to tease apart factors contributing to tolerance at a genomic level to and gain an understanding of the evolution of non-harmful in host-parasite interactions. To investigate if the fungal disease has caused adaptation on a genomic level in Eurasian bat species, we adopted both whole-genome sequencing approaches and a literature search to compile a set of 300 genes from which to investigate signals of positive selection in genomes of 11 Eurasian bats at the codon-level. Our results indicate significant positive selection in 38 genes, many of which have a marked role in responses to infection. Our findings suggest that white-nose syndrome may have applied a significant selective pressure on Eurasian Myotis-bats in the past, which can contribute their survival in co-existence with the pathogen. Our findings provide an insight on the selective pressure pathogens afflict on their hosts using methodology that can be adapted to other host-pathogen study systems.

Increased genome size is caused by heterochromatin addition in two non-related bat species, Hesperoptenus doriae and Philetor brachypterus (Vespertilionidae, Chiroptera, Mammalia).

The average genome size (GS) of bats, which are the only mammals capable of powered flight, is approximately 18% smaller than that of closely related mammalian orders. The low nuclear DNA content of Chiroptera is comparable to that of birds, which are also characterized by a high metabolic rate. Only a few chiropteran taxa possess notable amounts of constitutive heterochromatin. Here, we studied the karyotypes of two non-related vesper bat species with unusually high amounts of constitutive heterochromatin: Hesperoptenus doriae and Philetor brachypterus. Conventional staining methods and whole-chromosome painting with probes derived from Myotis myotis (2n = 44), showing a karyotype close to that of the presumed ancestor of Vespertilionidae, revealed Robertsonian fusions as the main type of rearrangement leading to the exceptionally reduced diploid chromosome number of 2n = 26 in both species. Moreover, both karyotypes are characterized by large blocks of pericentromeric heterochromatin composed of CMA-positive and DA-DAPI-positive segments. In H. doriae, the heterochromatin accumulation has resulted in a genome size of 3.22 pg (1C), which is 40% greater than the mean genome size for the family. For P. brachypterus, a genome size of 2.94 pg was determined, representing an increase of about 28%. Most notably, in H. doriae, the presence of additional constitutive heterochromatin correlates with an extended mitotic cell cycle duration in vitro. A reduction in diploid chromosome number to 30 or lower is discussed as a possible cause of the accumulation of pericentromeric heterochromatin in Vespertilionidae.

Chiropterans Are a Hotspot for Horizontal Transfer of DNA Transposons in Mammalia.

Horizontal transfer of transposable elements (TEs) is an important mechanism contributing to genetic diversity and innovation. Bats (order Chiroptera) have repeatedly been shown to experience horizontal transfer of TEs at what appears to be a high rate compared with other mammals. We investigated the occurrence of horizontally transferred (HT) DNA transposons involving bats. We found over 200 putative HT elements within bats; 16 transposons were shared across distantly related mammalian clades, and 2 other elements were shared with a fish and two lizard species. Our results indicate that bats are a hotspot for horizontal transfer of DNA transposons. These events broadly coincide with the diversification of several bat clades, supporting the hypothesis that DNA transposon invasions have contributed to genetic diversification of bats.

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

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

Musculoskeletal morphogenesis supports the convergent evolution of bat laryngeal echolocation.

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

Comparative genomics reveals convergent signals associated with the high metabolism and longevity in birds and bats.

Birds and bats have long lifespans relative to their body size compared with non-flying animals. However, the genomic basis associated with longer lifespan of flying species despite their higher metabolism was unclear. In this study, we hypothesized that genes involved in the regulation of metabolism and lifespan changed with the acquisition of flight and searched for genes that show specific evolutionary patterns in flying species. As a result, we identified several genes that show different evolutionary rates in bird and bat lineages. Genes in pathways involved in lifespan regulation were conserved in birds, while they evolved at an accelerated rate in bats. We also searched for genes in which convergent amino acid substitutions occurred in birds and bats and found such substitutions in genes involved in cancer, reactive oxygen species control and immunity. Our study revealed genomic changes associated with the acquisition of flight in birds and bats and suggested that multiple genes involved in the regulation of lifespan and metabolism support both high metabolism and longevity in flying species.

Cryptic hybridization between the ancient lineages of Natterer's bat (Myotis nattereri).

Studying hybrid zones that form between morphologically cryptic taxa offers valuable insights into the mechanisms of cryptic speciation and the evolution of reproductive barriers. Although hybrid zones have long been the focus of evolutionary studies, the awareness of cryptic hybrid zones increased recently due to rapidly growing evidence of biological diversity lacking obvious phenotypic differentiation. The characterization of cryptic hybrid zones with genome-wide analysis is in its early stages and offers new perspectives for studying population admixture and thus the impact of gene flow. In this study, we investigate the population genomics of the Myotis nattereri complex in one of its secondary contact zones, where a putative hybrid zone is formed between two of its cryptic lineages. By utilizing a whole-genome shotgun sequencing approach, we aim to characterize this cryptic hybrid zone in detail. Demographic analysis suggests that the cryptic lineages diverged during the Pliocene, c. 3.6 million years ago. Despite this ancient separation, the populations in the contact zone exhibit mitochondrial introgression and a considerable amount of mixing in nuclear genomes. The genomic structure of the populations corresponds to geographic locations and the genomic admixture changes along a geographic gradient. These findings suggest that there is no effective hybridization barrier between both lineages, nevertheless, their population structure is shaped by dispersal barriers. Our findings highlight how such deeply diverged cryptic lineages can still readily hybridize in secondary contact.

Phylogeography of horseshoe bat sarbecoviruses in Vietnam and neighbouring countries. Implications for the origins of SARS-CoV and SARS-CoV-2.

Previous studies on horseshoe bats (Rhinolophus spp.) have described many coronaviruses related to SARS-CoV (SARSCoVr) in China and only a few coronaviruses related to SARS-CoV-2 (SARSCoV2r) in Yunnan (southern China), Cambodia, Laos and Thailand. Here, we report the results of several field missions carried out in 2017, 2021 and 2022 across Vietnam during which 1218 horseshoe bats were sampled from 19 locations. Sarbecoviruses were detected in 11% of faecal RNA extracts, with much more positives among Rhinolophus thomasi (46%). We assembled 38 Sarbecovirus genomes, including 32 SARSCoVr, four SARSCoV2r, and two recombinants of SARSCoVr and SARSCoV2r (RecSar), one showing a Spike protein very similar to SARS-CoV-2. We detected a bat co-infected with four coronaviruses, including two sarbecoviruses. Our analyses revealed that Sarbecovirus genomes evolve in Vietnam under strong geographical and host constraints. First, we found evidence for a deep separation between viruses from northern Vietnam and those from central and southern Vietnam. Second, we detected only SARSCoVr in Rhinolophus thomasi, both SARSCoVr and SARSCoV2r in Rhinolophus affinis, and only RecSar in Rhinolophus pusillus captured close to the border with China. Third, the bias in favour of Uracil in synonymous third codon positions of SARSCoVr extracted from R. thomasi showed a negative correlation with latitudes. Our results also provided support for an emergence of SARS-CoV in horseshoe bats from northern Yunnan and emergence of SARS-CoV-2 in horseshoe bats from northern Indochina subtropical forests (southern Yunnan, northern Laos and north-western Vietnam).

RNA editing sites and triplet usage in exomes of bat RNA virus genomes of the family Paramyxoviridae.

Bats contain a diverse spectrum of viral species in their bodies. The RNA virus family Paramyxoviridae tends to infect several vertebrate species, which are accountable for a variety of devastating infections in both humans and animals. Viruses of this kind include measles, mumps, and Hendra. Some synonymous codons are favoured over others in mRNAs during gene-to-protein synthesis process. Such phenomenon is termed as codon usage bias (CUB). Our research emphasized many aspects that shape the CUB of genes in the Paramyxoviridae family found in bats. Here, the nitrogenous base A occurred the most. AT was found to be abundant in the coding sequences of the Paramyxoviridae family. RSCU data revealed that A or T ending codons occurred more frequently than predicted. Furthermore, 3 overrepresented codons (CAT, AGA, and GCA) and 7 underrepresented codons (CCG, TCG, CGC, CGG, CGT, GCG and ACG) were detected in the viral genomes. Correspondence analysis, neutrality plot, and parity plots highlight the combined impact of mutational pressure and natural selection on CUB. The neutrality plot of GC12 against GC3 yielded a regression coefficient value of 0.366, indicating that natural selection had a significant (63.4 %) impact. Moreover, RNA editing analysis was done, which revealed the highest frequency of C to T mutations. The results of our research revealed the pattern of codon usage and RNA editing sites in Paramyxoviridae genomes.