Exploring the genome and transcriptome of the cave nectar bat Eonycteris spelaea with PacBio long-read sequencing

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

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

Genome assembly and isoform analysis of a highly heterozygous New Zealand fisheries species, the tarakihi (Nemadactylus macropterus)

Although being some of the most valuable and heavily exploited wild organisms, few fisheries species have been studied at the whole-genome level. This is especially the case in New Zealand, where genomics resources are urgently needed to assist fisheries management. Here, we generated 55 Gb of short Illumina reads (92× coverage) and 73 Gb of long Nanopore reads (122×) to produce the first genome assembly of the marine teleost tarakihi [Nemadactylus macropterus (Forster, 1801)], a highly valuable fisheries species in New Zealand. An additional 300 Mb of Iso-Seq reads were obtained to assist in gene annotation. The final genome assembly was 568 Mb long with an N50 of 3.37 Mb. The genome completeness was high, with 97.8% of complete Actinopterygii Benchmarking Universal Single-Copy Orthologs. Heterozygosity values estimated through k-mer counting (1.00%) and bi-allelic SNPs (0.64%) were high compared with the same values reported for other fishes. Iso-Seq analysis recovered 91,313 unique transcripts from 15,515 genes (mean ratio of 5.89 transcripts per gene), and the most common alternative splicing event was intron retention. This highly contiguous genome assembly and the isoform-resolved transcriptome will provide a useful resource to assist the study of population genomics and comparative eco-evolutionary studies in teleosts and related organisms.

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.

Perspectives and Benefits of High-Throughput Long-Read Sequencing in Microbial Ecology

Short-read, high-throughput sequencing (HTS) methods have yielded numerous important insights into microbial ecology and function. Yet, in many instances short-read HTS techniques are suboptimal, for example, by providing insufficient phylogenetic resolution or low integrity of assembled genomes. Single-molecule and synthetic long-read (SLR) HTS methods have successfully ameliorated these limitations. In addition, nanopore sequencing has generated a number of unique analysis opportunities, such as rapid molecular diagnostics and direct RNA sequencing, and both Pacific Biosciences (PacBio) and nanopore sequencing support detection of epigenetic modifications. Although initially suffering from relatively low sequence quality, recent advances have greatly improved the accuracy of long-read sequencing technologies. In spite of great technological progress in recent years, the long-read HTS methods (PacBio and nanopore sequencing) are still relatively costly, require large amounts of high-quality starting material, and commonly need specific solutions in various analysis steps. Despite these challenges, long-read sequencing technologies offer high-quality, cutting-edge alternatives for testing hypotheses about microbiome structure and functioning as well as assembly of eukaryote genomes from complex environmental DNA samples.

Novel de Novo Genome of Cynopterus brachyotis Reveals Evolutionarily Abrupt Shifts in Gene Family Composition across Fruit Bats

Major novel physiological or phenotypic adaptations often require accompanying modifications at the genic level. Conversely, the detection of considerable contractions and/or expansions of gene families can be an indicator of fundamental but unrecognized physiological change. We sequenced a novel fruit bat genome (Cynopterus brachyotis) and adopted a comparative approach to reconstruct the evolution of fruit bats, mapping contractions and expansions of gene families along their evolutionary history. Despite a radical change in life history as compared with other bats (e.g., loss of echolocation, large size, and frugivory), fruit bats have undergone surprisingly limited change in their genic composition, perhaps apart from a potentially novel gene family expansion relating to telomere protection and longevity. In sharp contrast, within fruit bats, the new Cynopterus genome bears the signal of unusual gene loss and gene family contraction, despite its similar morphology and lifestyle to two other major fruit bat lineages. Most missing genes are regulatory, immune-related, and olfactory in nature, illustrating the diversity of genomic strategies employed by bats to contend with responses to viral infection and olfactory requirements. Our results underscore that significant fluctuations in gene family composition are not always associated with obvious examples of novel physiological and phenotypic adaptations but may often relate to less-obvious shifts in immune strategies.

Loss of sweet taste despite the conservation of sweet receptor genes in insectivorous bats

The sense of taste provides key information on diet, but evolution of taste receptor genes in vertebrates is sometimes unable to predict their feeding ecology. Here we use behavioral experiments and functional assays to demonstrate the loss of sweet taste despite the conservation of sweet receptor genes in insectivorous bats. Although sweet taste receptor genes were highly conserved between frugivorous and insectivorous bats at the sequence level, our behavioral experiments revealed dramatic divergence in two bat species with distinct diets: the insectivorous bat showed no preference for natural sugars, whereas the frugivorous bat showed strong preferences for sucrose and fructose. Our cell-based assays from multiple representative bat species across the phylogeny further supported the behavioral preference tests.

The evolution of taste perception is usually associated with the ecology and dietary changes of organisms. However, the association between feeding ecology and taste receptor evolution is unclear in some lineages of vertebrate animals. One example is the sweet taste receptor gene Tas1r2. Previous analysis of partial sequences has revealed that Tas1r2 has undergone equally strong purifying selection between insectivorous and frugivorous bats. To test whether the sweet taste function is also important in bats with contrasting diets, we examined the complete coding sequences of both sweet taste receptor genes (Tas1r2 and Tas1r3) in 34 representative bat species. Although these two genes are highly conserved between frugivorous and insectivorous bats at the sequence level, our behavioral experiments revealed that an insectivorous bat (Myotis ricketti) showed no preference for natural sugars, whereas the frugivorous species (Rousettus leschenaultii) showed strong preferences for sucrose and fructose. Furthermore, while both sweet taste receptor genes are expressed in the taste tissue of insectivorous and frugivorous bats, our cell-based assays revealed striking functional divergence: the sweet taste receptors of frugivorous bats are able to respond to natural sugars whereas those of insectivorous bats are not, which is consistent with the behavioral preference tests, suggesting that functional evolution of sweet taste receptors is closely related to diet. This comprehensive study suggests that using sequence conservation alone could be misleading in inferring protein and physiological function and highlights the power of combining behavioral experiments, expression analysis, and functional assays in molecular evolutionary studies.

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

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

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.

Immunophenotyping monocytes, macrophages and granulocytes in the Pteropodid bat Eonycteris spelaea

Bats are asymptomatic reservoir hosts for several highly pathogenic viruses. Understanding this enigmatic relationship between bats and emerging zoonotic viruses requires tools and approaches which enable the comparative study of bat immune cell populations and their functions. We show that bat genomes have a conservation of immune marker genes which delineate phagocyte populations in humans, while lacking key mouse surface markers such as Ly6C and Ly6G. Cross-reactive antibodies against CD44, CD11b, CD14, MHC II, and CD206 were multiplexed to characterize circulating monocytes, granulocytes, bone-marrow derived macrophages (BMDMs) and lung alveolar macrophages (AMs) in the cave nectar bat Eonycteris spelaea. Transcriptional profiling of bat monocytes and BMDMs identified additional markers – including MARCO, CD68, CD163, CD172α, and CD88 – which can be used to further characterize bat myeloid populations. Bat cells often resembled their human counterparts when comparing immune parameters that are divergent between humans and mice, such as the expression patterns of certain immune cell markers. A genome-wide comparison of immune-related genes also revealed a much closer phylogenetic relationship between bats and humans compared to rodents. Taken together, this study provides a set of tools and a comparative framework which will be important for unravelling viral disease tolerance mechanisms in bats.

The temporal RNA virome patterns of a lesser dawn bat (Eonycteris spelaea) colony revealed by deep sequencing

The virosphere is largely unexplored and the majority of viruses are yet to be represented in public sequence databases. Bats are rich reservoirs of viruses, including several zoonoses. In this study, high throughput sequencing (HTS) of viral RNA extracted from swabs of four body sites per bat per timepoint is used to characterize the virome through a longitudinal study of a captive colony of fruit nectar bats, species Eonycteris spelaea in Singapore. Through unbiased shotgun and target enrichment sequencing, we identify both known and previously unknown viruses of zoonotic relevance and define the population persistence and temporal patterns of viruses from families that have the capacity to jump the species barrier. To our knowledge, this is the first study that combines probe-based viral enrichment with HTS to create a viral profile from multiple swab sites on individual bats and their cohort. This work demonstrates temporal patterns of the lesser dawn bat virome, including several novel viruses. Given the known risk for bat-human zoonoses, a more complete understanding of the viral dynamics in South-eastern Asian bats has significant implications for disease prevention and control. The findings of this study will be of interest to U.S. Department of Defense personnel stationed in the Asia-Pacific region and regional public health laboratories engaged in emerging infectious disease surveillance efforts.

A comprehensive annotation and differential expression analysis of short and long non-coding RNAs in 16 bat genomes

Although bats are increasingly becoming the focus of scientific studies due to their unique properties, these exceptional animals are still among the least studied mammals. Assembly quality and completeness of bat genomes vary a lot and especially non-coding RNA (ncRNA) annotations are incomplete or simply missing. Accordingly, standard bioinformatics pipelines for gene expression analysis often ignore ncRNAs such as microRNAs or long antisense RNAs. The main cause of this problem is the use of incomplete genome annotations. We present a complete screening for ncRNAs within 16 bat genomes. NcRNAs affect a remarkable variety of vital biological functions, including gene expression regulation, RNA processing, RNA interference and, as recently described, regulatory processes in viral infections. Within all investigated bat assemblies, we annotated 667 ncRNA families including 162 snoRNAs and 193 miRNAs as well as rRNAs, tRNAs, several snRNAs and lncRNAs, and other structural ncRNA elements. We validated our ncRNA candidates by six RNA-Seq data sets and show significant expression patterns that have never been described before in a bat species on such a large scale. Our annotations will be usable as a resource (rna.uni-jena.de/supplements/bats) for deeper studying of bat evolution, ncRNAs repertoire, gene expression and regulation, ecology and important host–virus interactions.

Multiple Long-Read Sequencing Survey of Herpes Simplex Virus Dynamic Transcriptome

Long-read sequencing (LRS) has become increasingly important in RNA research due to its strength in resolving complex transcriptomic architectures. In this regard, currently two LRS platforms have demonstrated adequate performance: the Single Molecule Real-Time Sequencing by Pacific Biosciences (PacBio) and the nanopore sequencing by Oxford Nanopore Technologies (ONT). Even though these techniques produce lower coverage and are more error prone than short-read sequencing, they continue to be more successful in identifying polycistronic RNAs, transcript isoforms including splice and transcript end variants, as well as transcript overlaps. Recent reports have successfully applied LRS for the investigation of the transcriptome of viruses belonging to various families. These studies have substantially increased the number of previously known viral RNA molecules. In this work, we used the Sequel and MinION technique from PacBio and ONT, respectively, to characterize the lytic transcriptome of the herpes simplex virus type 1 (HSV-1). In most samples, we analyzed the poly(A) fraction of the transcriptome, but we also performed random oligonucleotide-based sequencing. Besides cDNA sequencing, we also carried out native RNA sequencing. Our investigations identified more than 2,300 previously undetected transcripts, including coding, and non-coding RNAs, multi-splice transcripts, as well as polycistronic and complex transcripts. Furthermore, we found previously unsubstantiated transcriptional start sites, polyadenylation sites, and splice sites. A large number of novel transcriptional overlaps were also detected. Random-primed sequencing revealed that each convergent gene pair produces non-polyadenylated read-through RNAs overlapping the partner genes. Furthermore, we identified novel replication-associated transcripts overlapping the HSV-1 replication origins, and novel LAT variants with very long 5' regions, which are co-terminal with the LAT-0.7kb transcript. Overall, our results demonstrated that the HSV-1 transcripts form an extremely complex pattern of overlaps, and that entire viral genome is transcriptionally active. In most viral genes, if not in all, both DNA strands are expressed.

The Impact of cDNA Normalization on Long-Read Sequencing of a Complex Transcriptome

Normalization of cDNA is widely used to improve the coverage of rare transcripts in analysis of transcriptomes employing next-generation sequencing. Recently, long-read technology has been emerging as a powerful tool for sequencing and construction of transcriptomes, especially for complex genomes containing highly similar transcripts and transcript-spliced isoforms. Here, we analyzed the transcriptome of sugarcane, a highly polyploidy plant genome, by PacBio isoform sequencing (Iso-Seq) of two different cDNA library preparations, with and without a normalization step. The results demonstrated that, while the two libraries included many of the same transcripts, many longer transcripts were removed, and many new generally shorter transcripts were detected by normalization. For the same input cDNA and data yield, the normalized library recovered more total transcript isoforms and number of predicted gene families and orthologous groups, resulting in a higher representation for the sugarcane transcriptome, compared to the non-normalized library. The non-normalized library, on the other hand, included a wider transcript length range with more longer transcripts above ∼1.25 kb and more transcript isoforms per gene family and gene ontology terms per transcript. A large proportion of the unique transcripts comprising ∼52% of the normalized library were expressed at a lower level than the unique transcripts from the non-normalized library, across three tissue types tested including leaf, stalk, and root. About 83% of the total 5,348 predicted long noncoding transcripts was derived from the normalized library, of which ∼80% was derived from the lowly expressed fraction. Functional annotation of the unique transcripts suggested that each library enriched different functional transcript fractions. This demonstrated the complementation of the two approaches in obtaining a complete transcriptome of a complex genome at the sequencing depth used in this study.

Diversity and Evolution of Viral Pathogen Community in Cave Nectar Bats (Eonycteris spelaea)

Bats are unique mammals, exhibit distinctive life history traits and have unique immunological approaches to suppression of viral diseases upon infection. High-throughput next-generation sequencing has been used in characterizing the virome of different bat species. The cave nectar bat, Eonycteris spelaea, has a broad geographical range across Southeast Asia, India and southern China, however, little is known about their involvement in virus transmission. Here we investigate the diversity and abundance of viral communities from a colony of Eonycteris spelaea residing in Singapore. Our results detected 47 and 22 different virus families from bat fecal and urine samples, respectively. Among these, we identify a large number of virus families including Adenoviridae, Flaviviridae, Reoviridae, Papillomaviridae, Paramyxoviridae, Parvoviridae, Picornaviridae, and Polyomaviridae. In most cases, viral sequences from Eonycteris spelaea are genetically related to a group of bat viruses from other bat genera (e.g., Eidolon, Miniopterus, Rhinolophus and Rousettus). The results of this study improve our knowledge of the host range, spread and evolution of several important viral pathogens. More significantly, our findings provide a baseline to study the temporal patterns of virus shedding and how they correlate with bat phenological trends.

Multifactorial processes underlie parallel opsin loss in neotropical bats

The loss of previously adaptive traits is typically linked to relaxation in selection, yet the molecular steps leading to such repeated losses are rarely known. Molecular studies of loss have tended to focus on gene sequences alone, but overlooking other aspects of protein expression might underestimate phenotypic diversity. Insights based almost solely on opsin gene evolution, for instance, have made mammalian color vision a textbook example of phenotypic loss. We address this gap by investigating retention and loss of opsin genes, transcripts, and proteins across ecologically diverse noctilionoid bats. We find multiple, independent losses of short-wave-sensitive opsins. Mismatches between putatively functional DNA sequences, mRNA transcripts, and proteins implicate transcriptional and post-transcriptional processes in the ongoing loss of S-opsins in some noctilionoid bats. Our results provide a snapshot of evolution in progress during phenotypic trait loss, and suggest vertebrate visual phenotypes cannot always be predicted from genotypes alone.

Bats are famous for using their hearing to explore their environments, yet fewer people are aware that these flying mammals have both good night and daylight vision. Some bats can even see in color thanks to two light-sensitive proteins at the back of their eyes: S-opsin which detects blue and ultraviolet light and L-opsin which detects green and red light. Many species of bat, however, are missing one of these proteins and cannot distinguish any colors; in other words, they are completely color-blind.

Some bat species found in Central and South America have independently lost their ability to see blue-ultraviolet light and have thus also lost their color vision. These bats have diverse diets – ranging from insects to fruits and even blood – and being able to distinguish color may offer an advantage in many of their activities, including hunting or foraging. The vision genes in these bats, therefore, give scientists an opportunity to explore how a seemingly important trait can be lost at the molecular level.

Sadier, Davies et al. now report that S-opsin has been lost more than a dozen times during the evolutionary history of these Central and South American bats. The analysis used samples from 55 species, including animals caught from the wild and specimens from museums. As with other proteins, the instructions encoded in the gene sequence for S opsin need to be copied into a molecule of RNA before they can be translated into protein. As expected, S-opsin was lost several times because of changes in the gene sequence that disrupted the formation of the protein. However, at several points in these bats’ evolutionary history, additional changes have taken place that affected the production of the RNA or the protein, without an obvious change to the gene itself. This finding suggests that other studies that rely purely on DNA to understand evolution may underestimate how often traits may be lost. By capturing ‘evolution in action’, these results also provide a more complete picture of the molecular targets of evolution in a diverse set of bats.

Exploring the genome and transcriptome of the cave nectar bat Eonycteris spelaea with PacBio long-read sequencing

Background

In the past two decades, bats have emerged as an important model system to study host-pathogen interactions. More recently, it has been shown that bats may also serve as a new and excellent model to study aging, inflammation, and cancer, among other important biological processes. The cave nectar bat or lesser dawn bat (Eonycteris spelaea) is known to be a reservoir for several viruses and intracellular bacteria. It is widely distributed throughout the tropics and subtropics from India to Southeast Asia and pollinates several plant species, including the culturally and economically important durian in the region. Here, we report the whole-genome and transcriptome sequencing, followed by subsequent de novo assembly, of the E. spelaea genome solely using the Pacific Biosciences (PacBio) long-read sequencing platform.

Findings

The newly assembled E. spelaea genome is 1.97 Gb in length and consists of 4,470 sequences with a contig N50 of 8.0 Mb. Identified repeat elements covered 34.65% of the genome, and 20,640 unique protein-coding genes with 39,526 transcripts were annotated.

Conclusions

We demonstrated that the PacBio long-read sequencing platform alone is sufficient to generate a comprehensive de novo assembled genome and transcriptome of an important bat species. These results will provide useful insights and act as a resource to expand our understanding of bat evolution, ecology, physiology, immunology, viral infection, and transmission dynamics.