Junpeng Zhang

Virome Analysis for Identification of Novel Mammalian Viruses in Bat Species from Chinese Provinces

ABSTRACT Bats are natural hosts for a large variety of zoonotic viruses. This study aimed to describe the range of bat viromes, including viruses from mammals, insects, fungi, plants, and phages, in 11 insectivorous bat species (216 bats in total) common in six provinces of China. To analyze viromes, we used sequence-independent PCR amplification and next-generation sequencing technology (Solexa Genome Analyzer II; Illumina). The viromes were identified by sequence similarity comparisons to known viruses. The mammalian viruses included those of the Adenoviridae, Herpesviridae, Papillomaviridae, Retroviridae, Circoviridae, Rhabdoviridae, Astroviridae, Flaviridae, Coronaviridae, Picornaviridae, and Parvovirinae; insect viruses included those of the Baculoviridae, Iflaviridae, Dicistroviridae, Tetraviridae, and Densovirinae; fungal viruses included those of the Chrysoviridae, Hypoviridae, Partitiviridae, and Totiviridae; and phages included those of the Caudovirales, Inoviridae, and Microviridae and unclassified phages. In addition to the viruses and phages associated with the insects, plants, and bacterial flora related to the diet and habitation of bats, we identified the complete or partial genome sequences of 13 novel mammalian viruses. These included herpesviruses, papillomaviruses, a circovirus, a bocavirus, picornaviruses, a pestivirus, and a foamy virus. Pairwise alignments and phylogenetic analyses indicated that these novel viruses showed little genetic similarity with previously reported viruses. This study also revealed a high prevalence and diversity of bat astroviruses and coronaviruses in some provinces. These findings have expanded our understanding of the viromes of bats in China and hinted at the presence of a large variety of unknown mammalian viruses in many common bat species of mainland China.

Deciphering the bat virome catalog to better understand the ecological diversity of bat viruses and the bat origin of emerging infectious diseases

Studies have demonstrated that ~60%–80% of emerging infectious diseases (EIDs) in humans originated from wild life. Bats are natural reservoirs of a large variety of viruses, including many important zoonotic viruses that cause severe diseases in humans and domestic animals. However, the understanding of the viral population and the ecological diversity residing in bat populations is unclear, which complicates the determination of the origins of certain EIDs. Here, using bats as a typical wildlife reservoir model, virome analysis was conducted based on pharyngeal and anal swab samples of 4440 bat individuals of 40 major bat species throughout China. The purpose of this study was to survey the ecological and biological diversities of viruses residing in these bat species, to investigate the presence of potential bat-borne zoonotic viruses and to evaluate the impacts of these viruses on public health. The data obtained in this study revealed an overview of the viral community present in these bat samples. Many novel bat viruses were reported for the first time and some bat viruses closely related to known human or animal pathogens were identified. This genetic evidence provides new clues in the search for the origin or evolution pattern of certain viruses, such as coronaviruses and noroviruses. These data offer meaningful ecological information for predicting and tracing wildlife-originated EIDs.

Historical male‐mediated introgression in horseshoe bats revealed by multilocus DNA sequence data

Instances of hybridization between mammalian taxa in the wild are rarely documented. To test for introgression between sibling species of horseshoe bat (Rhinolophus yunanensis and R. pearsoni) and two subspecies of the latter (R. p. pearsoni and R. p. chinensis), we sequenced two mtDNA and two ncDNA markers in individuals sampled from multiple localities within their overlapping ranges. The interspecific mtDNA gene tree corresponded to the expected taxonomic divisions, and coalescent‐based analyses suggested divergence occurred around 4 MYA. However, these relationships strongly conflicted with those recovered from two independent nuclear gene trees, in which R. yunanensis clustered with R. p. pearsoni to the exclusion of R. p. chinensis. This geographically widespread discordance is best explained by large‐scale historical introgression of ncDNA from R. yunanensis to R. pearsoni by male‐mediated exchange in mixed species colonies during Pleistocene glacial periods, when ranges may have contracted and overlapped more than at present. Further species tree–gene tree conflicts were detected between R. p. pearsoni and R. p. chinensis, also indicating past and/or current introgression in their overlapping regions. However, here the patterns point to asymmetric mtDNA introgression without ncDNA introgression. Analyses of coalescence times indicate this exchange has occurred subsequent to the divergence of these subspecies from their common ancestor. Our work highlights the importance of using multiple data sets for reconstructing phylogeographic histories and resolving taxonomic relationships.

Novel SARS-like betacoronaviruses in bats, China, 2011.

To clarify the evolutionary relationships among betavoronaviruses that infect bats, we analyzed samples collected during 2010–2011 from 14 insectivorous bat species in China. We identified complete genomes of 2 novel betacoronaviruses in Rhinolophus pusillus and Chaerephon plicata bats, which showed close genetic relationships with severe acute respiratory syndrome coronaviruses.

Unique expression patterns of multiple key genes associated with the evolution of mammalian flight.

Bats are the only mammals capable of true flight. Critical adaptations for flight include a pair of dramatically elongated hands with broad wing membranes. To study the molecular mechanisms of bat wing evolution, we perform genomewide mRNA sequencing and in situ hybridization for embryonic bat limbs. We identify seven key genes that display unique expression patterns in embryonic bat wings and feet, compared with mouse fore- and hindlimbs. The expression of all 5′HoxD genes (Hoxd9–13) and Tbx3, six known crucial transcription factors for limb and digit development, is extremely high and prolonged in the elongating wing area. The expression of Fam5c, a tumour suppressor, in bat limbs is bat-specific and significantly high in all short digit regions (the thumb and foot digits). These results suggest multiple genetic changes occurred independently during the evolution of bat wings to elongate the hand digits, promote membrane growth and keep other digits short. Our findings also indicate that the evolution of limb morphology depends on the complex integration of multiple gene regulatory networks and biological processes that control digit formation and identity, chondrogenesis, and interdigital regression or retention.

Adaptive Evolution in the Glucose Transporter 4 Gene Slc2a4 in Old World Fruit Bats (Family: Pteropodidae)

Frugivorous and nectarivorous bats are able to ingest large quantities of sugar in a short time span while avoiding the potentially adverse side-effects of elevated blood glucose. The glucose transporter 4 protein (GLUT4) encoded by the Slc2a4 gene plays a critical role in transmembrane skeletal muscle glucose uptake and thus glucose homeostasis. To test whether the Slc2a4 gene has undergone adaptive evolution in bats with carbohydrate-rich diets in relation to their insect-eating sister taxa, we sequenced the coding region of the Slc2a4 gene in a number of bat species, including four Old World fruit bats (Pteropodidae) and three New World fruit bats (Phyllostomidae). Our molecular evolutionary analyses revealed evidence that Slc2a4 has undergone a change in selection pressure in Old World fruit bats with 11 amino acid substitutions detected on the ancestral branch, whereas, no positive selection was detected in the New World fruit bats. We noted that in the former group, amino acid replacements were biased towards either Serine or Isoleucine, and, of the 11 changes, six were specific to Old World fruit bats (A133S, A164S, V377F, V386I, V441I and G459S). Our study presents preliminary evidence that the Slc2a4 gene has undergone adaptive changes in Old World fruit bats in relation to their ability to meet the demands of a high sugar diet.

ORF8-Related Genetic Evidence for Chinese Horseshoe Bats as the Source of Human Severe Acute Respiratory Syndrome Coronavirus

Abstract Several lineage B betacoronaviruses termed severe acute respiratory syndrome (SARS)–like CoVs (SL-CoVs) were identified from Rhinolophus bats in China. These viruses are characterized by a set of unique accessory open reading frames (ORFs) that are located between the M and N genes. Among unique accessory ORFs, ORF8 is most hypervariable. In this study, the ORF8s of all SL-CoVs were classified into 3 types, and, for the first time, it was found that very few SL-CoVs from Rhinolophus sinicus have ORF8s that are identical to that of human SARS-CoV. This finding provides new genetic evidence for Chinese horseshoe bats as the source of human SARS-CoV.

Lineage divergence and historical gene flow in the Chinese horseshoe bat (Rhinolophus sinicus).

Closely related taxa living in sympatry provide good opportunities to investigate the origin of barriers to gene flow as well as the extent of reproductive isolation. The only two recognized subspecies of the Chinese rufous horseshoe bat Rhinolophus sinicus are characterized by unusual relative distributions in which R. s. septentrionalis is restricted to a small area within the much wider range of its sister taxon R. s. sinicus. To determine the history of lineage divergence and gene flow between these taxa, we applied phylogenetic, demographic and coalescent analyses to multi-locus datasets. MtDNA gene genealogies and microsatellite-based clustering together revealed three divergent lineages of sinicus, corresponding to Central China, East China and the offshore Hainan Island. However, the central lineage of sinicus showed a closer relationship with septentrionalis than with other lineages of R. s. sinicus, in contrary to morphological data. Paraphyly of sinicus could result from either past asymmetric mtDNA introgression between these two taxa, or could suggest septentrionalis evolved in situ from its more widespread sister subspecies. To test between these hypotheses, we applied coalescent-based phylogenetic reconstruction and Approximate Bayesian Computation (ABC). We found that septentrionalis is likely to be the ancestral taxon and therefore a recent origin of this subspecies can be ruled out. On the other hand, we found a clear signature of asymmetric mtDNA gene flow from septentrionalis into central populations of sinicus yet no nuclear gene flow, thus strongly pointing to historical mtDNA introgression. We suggest that the observed deeply divergent lineages within R. sinicus probably evolved in isolation in separate Pleistocene refugia, although their close phylogeographic correspondence with distinct eco-environmental zones suggests that divergent selection might also have promoted broad patterns of population genetic structure.

First Knockdown Gene Expression in Bat (Hipposideros armiger) Brain Mediated by Lentivirus

Lentivirus-mediated RNA interference (RNAi) is a potent experimental tool for investigating gene functions in vitro and in vivo. It has advantages that transgenic technology lacks. However, in vivo applications are difficult to apply in the central nervous system of non-model organisms due to the lack of a standard brain atlas and genetic information. Here, we report the development of an in vivo gene delivery system used in bat brain tissue for the first time, based on lentivirus (LV) vectors expressing short hairpin RNA (shRNA) targeting Hipposideros armiger forkhead box P2 (FoxP2). In vitro transfection into HEK 293T cell with the vector bearing the cassettes encoding FoxP2 shRNA verified the knockdown efficiency. Pseudovirus particles were administered via stereotactic intracerebral microinjection into the anterior cingulate cortex of H. armiger. FoxP2 is of major interest because of its role in sensorimotor coordination and probably in echolocation. Subsequent in situ hybridization validated the in vivo silencing of the target gene. This report demonstrates that LV-mediated expression of RNAi could achieve effective gene silencing in bats, a non-model organism, and will assist in elucidating the functions of bat genes.

Prenatal development supports a single origin of laryngeal echolocation in bats

Bat laryngeal echolocation is considered as one of the most complex and diverse modes of auditory sensory perception in animals and its evolutionary history has been the cause of many scientific controversies in the past two decades. To date, the majority of scientific evidence supports that bats (Chiroptera) are divided into two subordinal groups: Yinpterochiroptera, containing the laryngeal echolocating superfamily Rhinolophidae as sister taxa to the non-laryngeal echolocating family Pteropodidae; and Yangochiroptera, containing all other laryngeal echolocating lineages. This topology has led to an unanswered question in mammalian biology: was laryngeal echolocation lost in the ancestral pteropodids or gained convergently in the echolocating bat lineages? To date, there is insufficient and conflicting evidence from fossil, genomic, morphological and phylogenomic data to resolve this question. We detail an ontogenetic study of fetal cochlear development from seven species of bats and five outgroup mammals and show that in early fetal development, all bats including the non-laryngeal echolocating pteropodids have a similarly large cochlea typically associated with laryngeal echolocation abilities. The subsequent cochlear growth rate in the pteropodids is the slowest of all mammals and leads to the pteropodids and the non-echolocating lineages eventually sharing a similar cochlear morphospace as adults. The results suggest that pteropodids maintain a vestigial developmental stage indicative of past echolocation capabilities and thus support a single origin of laryngeal echolocation in bats.