The "Vampirome": Transcriptome and proteome analysis of the principal and accessory submaxillary glands of the vampire bat Desmodus rotundus, a vector of human rabies

Integrative multi-omics analysis reveals the contribution of neoVTX genes to venom diversity of Synanceia verrucosa

BACKGROUND

Animal venom systems are considered as valuable model for investigating the molecular mechanisms underlying phenotypic evolution. Stonefish are the most venomous and dangerous fish because of severe human envenomation and occasionally fatalities, whereas the genomic background of their venom has not been fully explored compared with that in other venomous animals.

RESULTS

In this study, we followed modern venomic pipelines to decode the Synanceia verrucosa venom components. A catalog of 478 toxin genes was annotated based on our assembled chromosome-level genome. Integrative analysis of the high-quality genome, the transcriptome of the venom gland, and the proteome of crude venom revealed mechanisms underlying the venom complexity in S. verrucosa. Six tandem-duplicated neoVTX subunit genes were identified as the major source for the neoVTX protein production. Further isoform sequencing revealed massive alternative splicing events with a total of 411 isoforms demonstrated by the six genes, which further contributed to the venom diversity. We then characterized 12 dominantly expressed toxin genes in the venom gland, and 11 of which were evidenced to produce the venom protein components, with the neoVTX proteins as the most abundant. Other major venom proteins included a presumed CRVP, Kuntiz-type serine protease inhibitor, calglandulin protein, and hyaluronidase. Besides, a few of highly abundant non-toxin proteins were also characterized and they were hypothesized to function in housekeeping or hemostasis maintaining roles in the venom gland. Notably, gastrotropin like non-toxin proteins were the second highest abundant proteins in the venom, which have not been reported in other venomous animals and contribute to the unique venom properties of S. verrucosa.

CONCLUSIONS

The results identified the major venom composition of S. verrucosa, and highlighted the contribution of neoVTX genes to the diversity of venom composition through tandem-duplication and alternative splicing. The diverse neoVTX proteins in the venom as lethal particles are important for understanding the adaptive evolution of S. verrucosa. Further functional studies are encouraged to exploit the venom components of S. verrucosa for pharmaceutical innovation.

Different species of Chiroptera: Immune cells and molecules

The distinct composition and immune response characteristics of bats' innate and adaptive immune systems, which enable them to serve as host of numerous serious zoonotic viruses without falling ill, differ substantially from those of other mammals, it have garnered significant attention. In this article, we offer a systematic review of the names, attributes, and functions of innate and adaptive immune cells & molecules across different bat species. This includes descriptions of the differences shown by research between 71 bat species in 10 families, as well as comparisons between bats and other mammals. Studies of the immune cells & molecules of different bat species are necessary to understand the unique antiviral immunity of bats. By providing comprehensive information on these unique immune responses, it is hoped that new insights will be provided for the study of co-evolutionary dynamics between viruses and the bat immune system, as well as human antiviral immunity.

Gonadal Transcriptomic Analysis Reveals Novel Sex-Related Genes in Bactrocera dorsalis

Bactrocera dorsalis (Hendel) (Diptera: Tephritidae) is one of the most devastating agricultural pests worldwide due to its high reproductive and invasive abilities. The elucidation of its gonadal developmental characteristics and the identification of sex-related genes will provide a useful genetic basis for reproductive-based pest control. Here, the gonadal transcriptome of B. dorsalis was sequenced, and novel gonad-specific expressed genes were analyzed. A total of 1338, 336, 35, and 479 differentially expressed genes (DEGs) were found in the testis (TE), ovary (OV), female accessory gland (FAG), and male accessory gland (MAG), respectively. Furthermore, 463 highly expressed gonad-specific genes were identified, with the TE having the highest number of specific highly expressed genes, at 402, followed by 51 in the OV, 9 in the MAG, and only 1 in the FAG. Strikingly, approximately half of highly expressed gonad-specific genes were uncharacterized. Then, it was found that 35, 17, 3, 2, and 1 of 202 uncharacterized highly expressed TE-specific genes encoded proteins that contained transmembrane domains, signal peptides, high-mobility group boxes, the zinc finger domain, and the BTB/POZ domain, respectively. Interestingly, approximately 40% of uncharacterized highly expressed gonad-specific genes encoding proteins were not predicted to possess functional motifs or domains. Finally, the spatiotemporal expression and sequence characterization of six novel highly expressed gonad-specific genes were analyzed. Altogether, our findings provide a valuable dataset for future functional analyses of sex-related genes and potential target sites for pest control.

The evolution of sanguivory in vampire bats: origins and convergences

Blood-feeding (sanguivory) has evolved more than two dozen times among birds, fishes, insects, arachnids, molluscs, crustaceans, and annelids; however, among mammals, it is restricted to the vampire bats. Here, the authors revisit the question of how it evolved in that group. Evidence to date suggests that the ancestors of phyllostomids were insectivorous, and that carnivory, omnivory, and nectarivory evolved among phyllostomids after vampire bats diverged. Frugivory likely also evolved after vampire bats diverged, but the phylogeny is ambiguous on that point. However, vampire bats lack any genetic evidence of a frugivorous past, and the behavioural progression from frugivory to sanguivory is difficult to envision. Thus, the most parsimonious scenario is that sanguivory evolved in an insectivorous ancestor to vampire bats via ectoparasite-eating, wound-feeding, or some combination of the two—all feeding habits found among blood-feeding birds today. Comparing vampire bats with other sanguivores, the authors find several remarkable examples of convergence. Further, it was found that blood-feeding has been ca. 50 times more likely to evolve in a vertebrate lineage than in an invertebrate one. The authors hypothesize that this difference exists because vertebrates are more likely than invertebrates to have the biochemical necessities required to assimilate the components of vertebrate blood.

The Fast and the Furriest: Investigating the Rate of Selection on Mammalian Toxins

The evolution of venom and the selection pressures that act on toxins have been increasingly researched within toxinology in the last two decades, in part due to the exceptionally high rates of diversifying selection observed in animal toxins. In 2015, Sungar and Moran proposed the 'two-speed' model of toxin evolution linking evolutionary age of a group to the rates of selection acting on toxins but due to a lack of data, mammals were not included as less than 30 species of venomous mammal have been recorded, represented by elusive species which produce small amounts of venom. Due to advances in genomics and transcriptomics, the availability of toxin sequences from venomous mammals has been increasing. Using branch- and site-specific selection models, we present the rates of both episodic and pervasive selection acting upon venomous mammal toxins as a group for the first time. We identified seven toxin groups present within venomous mammals, representing Chiroptera, Eulipotyphla and Monotremata: KLK1, Plasminogen Activator, Desmallipins, PACAP, CRiSP, Kunitz Domain One and Kunitz Domain Two. All but one group (KLK1) was identified by our results to be evolving under both episodic and pervasive diversifying selection with four toxin groups having sites that were implicated in the fitness of the animal by TreeSAAP (Selection on Amino Acid Properties). Our results suggest that venomous mammal ecology, behaviour or genomic evolution are the main drivers of selection, although evolutionary age may still be a factor. Our conclusion from these results indicates that mammalian toxins are following the two-speed model of selection, evolving predominately under diversifying selection, fitting in with other younger venomous taxa like snakes and cone snails-with high amounts of accumulating mutations, leading to more novel adaptions in their toxins.

An inhibitor of leukotriene-A4 hydrolase from bat salivary glands facilitates virus infection

SignificanceAn immunosuppressant protein (MTX), which facilitates virus infection by inhibiting leukotriene A4 hydrolase (LTA4H) to produce the lipid chemoattractant leukotriene B4 (LTB4), was identified and characterized from the submandibular salivary glands of the bat Myotis pilosus. To the best of our knowledge, this is a report of an endogenous LTA4H inhibitor in animals. MTX was highly concentrated in the bat salivary glands, suggesting a mechanism for the generation of immunological privilege and immune tolerance and providing evidence of viral shedding through oral secretions. Moreover, given that the immunosuppressant MTX selectively inhibited the proinflammatory activity of LTA4H, without affecting its antiinflammatory activity, MTX might be a potential candidate for the development of antiinflammatory drugs by targeting the LTA4-LTA4H-LTB4 inflammatory axis.

Nectar-feeding bats and birds show parallel molecular adaptations in sugar metabolism enzymes

In most vertebrates, the demand for glucose as the primary substrate for cellular respiration is met by the breakdown of complex carbohydrates, or energy is obtained by protein and lipid catabolism. In contrast, a few bat and bird species have convergently evolved to subsist on nectar, a sugar-rich mixture of glucose, fructose, and sucrose.^1-4 How these nectar-feeders have adapted to cope with life-long high sugar intake while avoiding the onset of metabolic syndrome and diabetes^5-7 is not understood. We analyzed gene sequences obtained from 127 taxa, including 22 nectar-feeding bat and bird genera that collectively encompass four independent origins of nectarivory. We show these divergent taxa have undergone pervasive molecular adaptation in sugar catabolism pathways, including parallel selection in key glycolytic and fructolytic enzymes. We also uncover convergent amino acid substitutions in the otherwise evolutionarily conserved aldolase B (ALDOB), which catalyzes rate-limiting steps in fructolysis and glycolysis, and the mitochondrial gatekeeper pyruvate dehydrogenase (PDH), which links glycolysis and the tricarboxylic acid cycle. Metabolomic profile and enzyme functional assays are consistent with increased respiratory flux in nectar-feeding bats and help explain how these taxa can both sustain hovering flight and efficiently clear simple sugars. Taken together, our results indicate that nectar-feeding bats and birds have undergone metabolic adaptations that have enabled them to exploit a unique energy-rich dietary niche among vertebrates.

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

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

Surveying the Vampire Bat (Desmodus rotundus) Serum Proteome: A Resource for Identifying Immunological Proteins and Detecting Pathogens

Bats are increasingly studied as model systems for longevity and as natural hosts for some virulent viruses. Yet the ability to characterize immune mechanisms of viral tolerance and to quantify infection dynamics in wild bats is often limited by small sample volumes and few species-specific reagents. Here, we demonstrate how proteomics can overcome these limitations by using data-independent acquisition-based shotgun proteomics to survey the serum proteome of 17 vampire bats (Desmodus rotundus) from Belize. Using just 2 μL of sample and relatively short separations of undepleted serum digests, we identified 361 proteins across 5 orders of magnitude. Levels of immunological proteins in vampire bat serum were then compared to human plasma via published databases. Of particular interest were antiviral and antibacterial components, circulating 20S proteasome complex and proteins involved in redox activity. Lastly, we used known virus proteomes to putatively identify Rh186 from Macacine herpesvirus 3 and ORF1a from Middle East respiratory syndrome-related coronavirus, indicating that mass spectrometry-based techniques show promise for pathogen detection. Overall, these results can be used to design targeted mass-spectrometry assays to quantify immunological markers and detect pathogens. More broadly, our findings also highlight the application of proteomics in advancing wildlife immunology and pathogen surveillance.

Identification of Novel Toxin Genes from the Stinging Nettle Caterpillar Parasa lepida (Cramer, 1799): Insights into the Evolution of Lepidoptera Toxins

Simple Summary

Many caterpillar species can produce toxins that cause harmful reactions to humans, varying from mild irritation to death. Currently, there is very limited knowledge about caterpillar toxin diversity, because only a few species have been investigated. We used the transcriptome technique to identify candidate toxin genes from the nettle caterpillar Parasa lepida (Cramer, 1799). It is a common pest of oil palm, coconut, and mango in South and South-East Asia, which can cause severe pain and allergic responses to those in contact with them. We reported 168 candidate toxin genes. Most of them are members of the toxin genes families commonly recruited in animal venoms such as serine protease and serine protease inhibitors. However, we identified 21 novel genes encoding knottin-like peptides expressed at a high level in the transcriptome. Their predicted 3D structures are similar to neurotoxins in scorpion and tarantula. Our study suggests that P. lepida venom contains diverse toxin proteins that potentially cause allergic reactions and pain. This study sheds light on the hidden diversity of toxin proteins in caterpillar lineage, which could be future fruitful new drug sources.

Abstract

Many animal species can produce venom for defense, predation, and competition. The venom usually contains diverse peptide and protein toxins, including neurotoxins, proteolytic enzymes, protease inhibitors, and allergens. Some drugs for cancer, neurological disorders, and analgesics were developed based on animal toxin structures and functions. Several caterpillar species possess venoms that cause varying effects on humans both locally and systemically. However, toxins from only a few species have been investigated, limiting the full understanding of the Lepidoptera toxin diversity and evolution. We used the RNA-seq technique to identify toxin genes from the stinging nettle caterpillar, Parasa lepida (Cramer, 1799). We constructed a transcriptome from caterpillar urticating hairs and reported 34,968 unique transcripts. Using our toxin gene annotation pipeline, we identified 168 candidate toxin genes, including protease inhibitors, proteolytic enzymes, and allergens. The 21 P. lepida novel Knottin-like peptides, which do not show sequence similarity to any known peptide, have predicted 3D structures similar to tarantula, scorpion, and cone snail neurotoxins. We highlighted the importance of convergent evolution in the Lepidoptera toxin evolution and the possible mechanisms. This study opens a new path to understanding the hidden diversity of Lepidoptera toxins, which could be a fruitful source for developing new drugs.

A Transcriptomic Approach to the Recruitment of Venom Proteins in a Marine Annelid

The growing number of known venomous marine invertebrates indicates that chemical warfare plays an important role in adapting to diversified ecological niches, even though it remains unclear how toxins fit into the evolutionary history of these animals. Our case study, the Polychaeta Eulalia sp., is an intertidal predator that secretes toxins. Whole-transcriptome sequencing revealed proteinaceous toxins secreted by cells in the proboscis and delivered by mucus. Toxins and accompanying enzymes promote permeabilization, coagulation impairment and the blocking of the neuromuscular activity of prey upon which the worm feeds by sucking pieces of live flesh. The main neurotoxins ("phyllotoxins") were found to be cysteine-rich proteins, a class of substances ubiquitous among venomous animals. Some toxins were phylogenetically related to Polychaeta, Mollusca or more ancient groups, such as Cnidaria. Some toxins may have evolved from non-toxin homologs that were recruited without the reduction in molecular mass and increased specificity of other invertebrate toxins. By analyzing the phylogeny of toxin mixtures, we show that Polychaeta is uniquely positioned in the evolution of animal venoms. Indeed, the phylogenetic models of mixed or individual toxins do not follow the expected eumetazoan tree-of-life and highlight that the recruitment of gene products for a role in venom systems is complex.

Developing critical thinking in STEM education through inquiry-based writing in the laboratory classroom

Laboratory pedagogy is moving away from step-by-step instructions and toward inquiry-based learning, but only now developing methods for integrating inquiry-based writing (IBW) practices into the laboratory course. Based on an earlier proposal (Science 2011;332:919), we designed and implemented an IBW sequence in a university bioinformatics course. We automatically generated unique, double-blinded, biologically plausible DNA sequences for each student. After guided instruction, students investigated sequences independently and responded through IBW writing assignments. IBW assignments were structured as condensed versions of a scientific research article, and because the sequences were double blinded, they were also assessed as authentic science and evaluated on clarity and persuasiveness. We piloted the approach in a seven-day workshop (35 students) at Perdana University in Malaysia. We observed dramatically improved student engagement and indirect evidence of improved learning outcomes over a similar workshop without IBW. Based on student feedback, initial discomfort with the writing component abated in favor of an overall positive response and increasing comfort with the high demands of student writing. Similarly, encouraging results were found in a semester length undergraduate module at the National University of Singapore (155 students).

Adaptive Evolution of C-Type Lysozyme in Vampire Bats

In mammals, chicken-type (c-type) lysozymes are part of the innate immune system, killing bacteria by degrading peptidoglycan in their cell walls. Many of the studies on the evolution of c-type lysozymes have focused on its new digestive function, including the duplicated stomach lysozymes in ruminants. Similarly, in bats, gene duplications and subsequent adaptive evolution of c-type lysozyme have been reported in a clade of insectivorous species, which might have been driven by the need to digest chitin. However, no studies on the evolution of c-type lysozyme have been carried out in the second largest and dietary diverse bat family Phyllostomidae, which includes insectivorous, frugivorous, nectarivorous and sanguivorous species. Here, we sequenced and analyzed c-type lysozyme genes from four phyllostomid bats, the common vampire bat, the white-winged vampire bat, the lesser long-nosed bat and the big fruit-eating bat. Only a single lysozyme gene was identified in each of these species. Evidence for positive selection on mature lysozyme was found on lineages leading to vampire bats, but not other bats with single copy lysozyme genes. Moreover, several amino acid substitutions found in mature lysozymes from the sanguivorous clade are predicted to have functional impacts, adding further evidence for the adaptive evolution of lysozyme in vampire bats. Functional adaptation of vampire bat lysozymes could be associated with anti-microbial defense, possibly driven by the specialized sanguivory-related habits of vampire bats.

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.

De Novo Transcriptome Assembly and Functional Annotation in Five Species of Bats

High-throughput RNA sequencing is a powerful tool that allows us to perform gene prediction and analyze tissue-specific overexpression of genes, but also at species level comparisons can be performed, although in a more restricted manner. In the present study complete liver transcriptomes of five tropical bat species were De novo assembled and annotated. Highly expressed genes in the five species were involved in glycolysis and lipid metabolism pathways. Cross-species differential expression analysis was conducted using single copy orthologues shared across the five species. Between 22 and 29 orthologs were upregulated for each species. We detected upregulated expression in Artibeus jamaicensis genes related to fructose metabolism pathway. Such findings can be correlated with A. jamaicensis dietary habits, as it was the unique frugivorous species included. This is the first report of transcriptome assembly by RNA-seq in these species, except for A. jamaicensis and as far as our knowledge is the first cross-species comparisons of transcriptomes and gene expression in tropical bats.

Chitinase genes (CHIAs) provide genomic footprints of a post-Cretaceous dietary radiation in placental mammals

The end-Cretaceous extinction led to a massive faunal turnover, with placental mammals radiating in the wake of nonavian dinosaurs. Fossils indicate that Cretaceous stem placentals were generally insectivorous, whereas their earliest Cenozoic descendants occupied a variety of dietary niches. It is hypothesized that this dietary radiation resulted from the opening of niche space, following the extinction of dinosaurian carnivores and herbivores. We provide the first genomic evidence for the occurrence and timing of this dietary radiation in placental mammals. By comparing the genomes of 107 placental mammals, we robustly infer that chitinase genes (CHIAs), encoding enzymes capable of digesting insect exoskeletal chitin, were present as five functional copies in the ancestor of all placental mammals, and the number of functional CHIAs in the genomes of extant species positively correlates with the percentage of invertebrates in their diets. The diverse repertoire of CHIAs in early placental mammals corroborates fossil evidence of insectivory in Cretaceous eutherians, with descendant lineages repeatedly losing CHIAs beginning at the Cretaceous/Paleogene (K/Pg) boundary as they radiated into noninsectivorous niches. Furthermore, the timing of gene loss suggests that interordinal diversification of placental mammals in the Cretaceous predates the dietary radiation in the early Cenozoic, helping to reconcile a long-standing debate between molecular timetrees and the fossil record. Our results demonstrate that placental mammal genomes, including humans, retain a molecular record of the post-K/Pg placental adaptive radiation in the form of numerous chitinase pseudogenes.

Hologenomic adaptations underlying the evolution of sanguivory in the common vampire bat

Adaptation to specialized diets often requires modifications at both genomic and microbiome levels. We applied a hologenomic approach to the common vampire bat (Desmodus rotundus), one of the only three obligate blood-feeding (sanguivorous) mammals, to study the evolution of its complex dietary adaptation. Specifically, we assembled its high-quality reference genome (scaffold N50 = 26.9 Mb, contig N50 = 36.6 kb) and gut metagenome, and compared them against those of insectivorous, frugivorous and carnivorous bats. Our analyses showed a particular common vampire bat genomic landscape regarding integrated viral elements, a dietary and phylogenetic influence on gut microbiome taxonomic and functional profiles, and that both genetic elements harbour key traits related to the nutritional (for example, vitamin and lipid shortage) and non-nutritional (for example, nitrogen waste and osmotic homeostasis) challenges of sanguivory. These findings highlight the value of a holistic study of both the host and its microbiota when attempting to decipher adaptations underlying radical dietary lifestyles.

Are Fireworms Venomous? Evidence for the Convergent Evolution of Toxin Homologs in Three Species of Fireworms (Annelida, Amphinomidae)

Amphinomids, more commonly known as fireworms, are a basal lineage of marine annelids characterized by the presence of defensive dorsal calcareous chaetae, which break off upon contact. It has long been hypothesized that amphinomids are venomous and use the chaetae to inject a toxic substance. However, studies investigating fireworm venom from a morphological or molecular perspective are scarce and no venom gland has been identified to date, nor any toxin characterized at the molecular level. To investigate this question, we analyzed the transcriptomes of three species of fireworms-Eurythoe complanata, Hermodice carunculata, and Paramphinome jeffreysii-following a venomics approach to identify putative venom compounds. Our venomics pipeline involved de novo transcriptome assembly, open reading frame, and signal sequence prediction, followed by three different homology search strategies: BLAST, HMMER sequence, and HMMER domain. Following this pipeline, we identified 34 clusters of orthologous genes, representing 13 known toxin classes that have been repeatedly recruited into animal venoms. Specifically, the three species share a similar toxin profile with C-type lectins, peptidases, metalloproteinases, spider toxins, and CAP proteins found among the most highly expressed toxin homologs. Despite their great diversity, the putative toxins identified are predominantly involved in three major biological processes: hemostasis, inflammatory response, and allergic reactions, all of which are commonly disrupted after fireworm stings. Although the putative fireworm toxins identified here need to be further validated, our results strongly suggest that fireworms are venomous animals that use a complex mixture of toxins for defense against predators.

De novo transcriptome assembly for the spiny mouse (Acomys cahirinus)

Spiny mice of the genus Acomys display several unique physiological traits, including menstruation and scar-free wound healing; characteristics that are exceedingly rare in mammals, and of considerable interest to the scientific community. These unique attributes, and the potential for spiny mice to accurately model human diseases, are driving increased use of this genus in biomedical research, however little genetic information is accessible for this species. This project aimed to generate a draft transcriptome for the Common spiny mouse (Acomys cahirinus). Illumina sequencing of RNA from 15 organ types (male and female) produced 451 million, 150 bp paired-end reads (92.4Gbp). An extensive survey of de novo transcriptome assembly approaches using Trinity, SOAPdenovo-Trans, and Oases at multiple kmer lengths was conducted, producing 50 single-kmer assemblies from this dataset. Non-redundant transcripts from all assemblies were merged into a meta-assembly using the EvidentialGene tr2aacds pipeline, producing the largest gene catalogue to date for Acomys cahirinus. This study provides the first detailed characterization of the spiny mouse transcriptome. It validates use of the EvidentialGene tr2aacds pipeline in mammals to augment conventional de novo assembly approaches, and provides a valuable scientific resource for further investigation into the unique physiological characteristics inherent in the genus Acomys.

Comparative genomics reveals contraction in olfactory receptor genes in bats

Gene loss and gain during genome evolution are thought to play important roles in adaptive phenotypic diversification. Among mammals, bats possess the smallest genomes and have evolved the unique abilities of powered flight and laryngeal echolocation. To investigate whether gene family evolution has contributed to the genome downsizing and phenotypic diversification in this group, we performed comparative evolutionary analyses of complete proteome data for eight bat species, including echolocating and non-echolocating forms, together with the proteomes of 12 other laurasiatherian mammals. Our analyses revealed extensive gene loss in the most recent ancestor of bats, and also of carnivores (both >1,000 genes), although this gene contraction did not appear to correlate with the reduction in genome size in bats. Comparisons of highly dynamic families suggested that expansion and contraction affected genes with similar functions (immunity, response to stimulus) in all laurasiatherian lineages. However, the magnitude and direction of these changes varied greatly among groups. In particular, our results showed contraction of the Olfactory Receptor (OR) gene repertoire in the last common ancestor of all bats, as well as that of the echolocating species studied. In contrast, non-echolocating fruit bats showed evidence of expansion in ORs, supporting a "trade-off" between sensory modalities.

Systems Biomedicine of Rabies Delineates the Affected Signaling Pathways

The prototypical neurotropic virus, rabies, is a member of the Rhabdoviridae family that causes lethal encephalomyelitis. Although there have been a plethora of studies investigating the etiological mechanism of the rabies virus and many precautionary methods have been implemented to avert the disease outbreak over the last century, the disease has surprisingly no definite remedy at its late stages. The psychological symptoms and the underlying etiology, as well as the rare survival rate from rabies encephalitis, has still remained a mystery. We, therefore, undertook a systems biomedicine approach to identify the network of gene products implicated in rabies. This was done by meta-analyzing whole-transcriptome microarray datasets of the CNS infected by strain CVS-11, and integrating them with interactome data using computational and statistical methods. We first determined the differentially expressed genes (DEGs) in each study and horizontally integrated the results at the mRNA and microRNA levels separately. A total of 61 seed genes involved in signal propagation system were obtained by means of unifying mRNA and microRNA detected integrated DEGs. We then reconstructed a refined protein–protein interaction network (PPIN) of infected cells to elucidate the rabies-implicated signal transduction network (RISN). To validate our findings, we confirmed differential expression of randomly selected genes in the network using Real-time PCR. In conclusion, the identification of seed genes and their network neighborhood within the refined PPIN can be useful for demonstrating signaling pathways including interferon circumvent, toward proliferation and survival, and neuropathological clue, explaining the intricate underlying molecular neuropathology of rabies infection and thus rendered a molecular framework for predicting potential drug targets.

Vasoconstriction induced by salivary gland extracts from ixodid ticks

In their quest for blood, most haematophagous parasites secrete vasodilators in their saliva to counter the host haemostatic response of vasoconstriction. Surprisingly, salivary gland extracts from adult female Dermacentor reticulatus and Rhipicephalus appendiculatus ticks induced constriction in a rat femoral artery model; males induced vasoconstriction or vasodilation depending on the time of feeding. Based on comparative HPLC fractionation, the active compounds inducing vasoconstriction do not appear to be prostaglandins (which ticks normally use as vasodilators). Vasoconstriction may be unique to ixodid ticks, helping them control blood flow during their prolonged blood-feeding of up to 10 days or more.

Arachnids of medical importance in Brazil: main active compounds present in scorpion and spider venoms and tick saliva

Arachnida is the largest class among the arthropods, constituting over 60,000 described species (spiders, mites, ticks, scorpions, palpigrades, pseudoscorpions, solpugids and harvestmen). Many accidents are caused by arachnids, especially spiders and scorpions, while some diseases can be transmitted by mites and ticks. These animals are widely dispersed in urban centers due to the large availability of shelter and food, increasing the incidence of accidents. Several protein and non-protein compounds present in the venom and saliva of these animals are responsible for symptoms observed in envenoming, exhibiting neurotoxic, dermonecrotic and hemorrhagic activities. The phylogenomic analysis from the complementary DNA of single-copy nuclear protein-coding genes shows that these animals share some common protein families known as neurotoxins, defensins, hyaluronidase, antimicrobial peptides, phospholipases and proteinases. This indicates that the venoms from these animals may present components with functional and structural similarities. Therefore, we described in this review the main components present in spider and scorpion venom as well as in tick saliva, since they have similar components. These three arachnids are responsible for many accidents of medical relevance in Brazil. Additionally, this study shows potential biotechnological applications of some components with important biological activities, which may motivate the conducting of further research studies on their action mechanisms.

Evidence for a lectin specific for sulfated glycans in the salivary gland of the malaria vector, Anopheles gambiae

Salivary gland homogenate (SGH) from the female mosquitoes Anopheles gambiae, An. stephensi, An. freeborni, An. dirus and An. albimanus were found to exhibit hemagglutinating (lectin) activity. Lectin activity was not found for male An. gambiae, or female Ae aegypti, Culex quinquefasciatus, Phlebotomus duboscqi, and Lutzomyia longipalpis. With respect to species-specificity, An. gambiae SGH agglutinates red blood cells (RBC) from humans, horse, sheep, goat, pig, and cow; it is less active for rats RBC, and not detectable for guinea-pigs or chicken RBC. Notably, lectin activity was inhibited by low concentrations of dextran sulfate 50–500 K, fucoidan, heparin, laminin, heparin sulfate proteoglycan, sialyl-containing glycans (e.g. 3′-sialyl Lewis X, and 6′-sialyl lactose), and gangliosides (e.g. GM3, GD1, GD1b, GTB1, GM1, GQ1B), but not by simple sugars. These results imply that molecule(s) in the salivary gland target sulfated glycans. SGH from An. gambiae was also found to promote agglutination of HL-60 cells which are rich in sialyl Lewis X, a glycan that decorates PSGL-1, the neutrophils receptor that interacts with endothelial cell P-selectin. Accordingly, SGH interferes with HL-60 cells adhesion to immobilized P-selectin. Because An. gambiae SGH expresses galectins, one member of this family (herein named Agalectin) was expressed in E. coli. Recombinant Agalectin behaves as a non-covalent homodimer. It does not display lectin activity, and does not interact with 500 candidates tested in a Glycan microarray. Gel-filtration chromatography of the SGH of An. gambiae identified a fraction with hemagglutinating activity, which was analyzed by 1D PAGE followed by in-gel tryptic digestion, and nano-LC MS/MS. This approach identified several genes which emerge as candidates for a lectin targeting sulfated glycans, the first with this selectivity to be reported in the SGH of a blood-sucking arthropod. The role of salivary molecules (sialogenins) with lectin activity is discussed in the context of inflammation, and parasite-vector-host interactions.

A Polychaete's powerful punch: venom gland transcriptomics of Glycera reveals a complex cocktail of toxin homologs

Glycerids are marine annelids commonly known as bloodworms. Bloodworms have an eversible proboscis adorned with jaws connected to venom glands. Bloodworms prey on invertebrates, and it is known that the venom glands produce compounds that can induce toxic effects in animals. Yet, none of these putative toxins has been characterized on a molecular basis. Here we present the transcriptomic profiles of the venom glands of three species of bloodworm, Glycera dibranchiata, Glycera fallax and Glycera tridactyla, as well as the body tissue of G. tridactyla. The venom glands express a complex mixture of transcripts coding for putative toxin precursors. These transcripts represent 20 known toxin classes that have been convergently recruited into animal venoms, as well as transcripts potentially coding for Glycera-specific toxins. The toxins represent five functional categories: Pore-forming and membrane-disrupting toxins, neurotoxins, protease inhibitors, other enzymes, and CAP domain toxins. Many of the transcripts coding for putative Glycera toxins belong to classes that have been widely recruited into venoms, but some are homologs of toxins previously only known from the venoms of scorpaeniform fish and monotremes (stonustoxin-like toxin), turrid gastropods (turripeptide-like peptides), and sea anemones (gigantoxin I-like neurotoxin). This complex mixture of toxin homologs suggests that bloodworms employ venom while predating on macroscopic prey, casting doubt on the previously widespread opinion that G. dibranchiata is a detritivore. Our results further show that researchers should be aware that different assembly methods, as well as different methods of homology prediction, can influence the transcriptomic profiling of venom glands.

Are ticks venomous animals?

INTRODUCTION

As an ecological adaptation venoms have evolved independently in several species of Metazoa. As haematophagous arthropods ticks are mainly considered as ectoparasites due to directly feeding on the skin of animal hosts. Ticks are of major importance since they serve as vectors for several diseases affecting humans and livestock animals. Ticks are rarely considered as venomous animals despite that tick saliva contains several protein families present in venomous taxa and that many Ixodida genera can induce paralysis and other types of toxicoses. Tick saliva was previously proposed as a special kind of venom since tick venom is used for blood feeding that counteracts host defense mechanisms. As a result, the present study provides evidence to reconsider the venomous properties of tick saliva.

RESULTS

Based on our extensive literature mining and in silico research, we demonstrate that ticks share several similarities with other venomous taxa. Many tick salivary protein families and their previously described functions are homologous to proteins found in scorpion, spider, snake, platypus and bee venoms. This infers that there is a structural and functional convergence between several molecular components in tick saliva and the venoms from other recognized venomous taxa. We also highlight the fact that the immune response against tick saliva and venoms (from recognized venomous taxa) are both dominated by an allergic immunity background. Furthermore, by comparing the major molecular components of human saliva, as an example of a non-venomous animal, with that of ticks we find evidence that ticks resemble more venomous than non-venomous animals. Finally, we introduce our considerations regarding the evolution of venoms in Arachnida.

CONCLUSIONS

Taking into account the composition of tick saliva, the venomous functions that ticks have while interacting with their hosts, and the distinguishable differences between human (non-venomous) and tick salivary proteins, we consider that ticks should be referred to as venomous ectoparasites.

Desmolaris, a novel factor XIa anticoagulant from the salivary gland of the vampire bat (Desmodus rotundus) inhibits inflammation and thrombosis in vivo

The identity of vampire bat saliva anticoagulant remained elusive for almost a century. Sequencing the salivary gland genes from the vampire bat Desmodus rotundus identified Desmolaris as a novel 21.5-kDa naturally deleted (Kunitz 1-domainless) form of tissue factor pathway inhibitor. Recombinant Desmolaris was expressed in HEK293 cells and characterized as a slow, tight, and noncompetitive inhibitor of factor (F) XIa by a mechanism modulated by heparin. Desmolaris also inhibits FXa with lower affinity, independently of protein S. In addition, Desmolaris binds kallikrein and reduces bradykinin generation in plasma activated with kaolin. Truncated and mutated forms of Desmolaris determined that Arg32 in the Kunitz-1 domain is critical for protease inhibition. Moreover, Kunitz-2 and the carboxyl-terminus domains mediate interaction of Desmolaris with heparin and are required for optimal inhibition of FXIa and FXa. Notably, Desmolaris (100 μg/kg) inhibited FeCl3-induced carotid artery thrombus without impairing hemostasis. These results imply that FXIa is the primary in vivo target for Desmolaris at antithrombotic concentrations. Desmolaris also reduces the polyphosphate-induced increase in vascular permeability and collagen- and epinephrine-mediated thromboembolism in mice. Desmolaris emerges as a novel anticoagulant targeting FXIa under conditions in which the coagulation activation, particularly the contact pathway, plays a major pathological role.

Dracula's children: molecular evolution of vampire bat venom

While vampire bat oral secretions have been the subject of intense research, efforts have concentrated only on two components: DSPA (Desmodus rotundus salivary plasminogen activator) and Draculin. The molecular evolutionary history of DSPA has been elucidated, while conversely draculin has long been known from only a very small fragment and thus even the basic protein class was not even established. Despite the fact that vampire bat venom has a multitude of effects unaccounted by the documented bioactivities of DSPA and draculin, efforts have not been made to establish what other bioactive proteins are secreted by their submaxillary gland. In addition, it has remained unclear whether the anatomically distinct anterior and posterior lobes of the submaxillary gland are evolving on separate gene expression trajectories or if they remain under the shared genetic control. Using a combined proteomic and transcriptomic approach, we show that identical proteins are simultaneously expressed in both lobes. In addition to recovering the known structural classes of DSPA, we recovered a novel DSPA isoform as well as obtained a very large sequence stretch of draculin and thus established that it is a mutated version of the lactotransferrin scaffold. This study reveals a much more complex secretion profile than previously recognised. In addition to obtaining novel versions of scaffolds convergently recruited into other venoms (allergen-like, CRiSP, kallikrein, Kunitz, lysozyme), we also documented novel expression of small peptides related to calcitonin, PACAP, and statherin. Other overexpressed protein types included BPI-fold, lacritin, and secretoglobin. Further, we investigate the molecular evolution of various vampire bat venom-components and highlight the dominant role of positive selection in the evolution of these proteins. Conspicuously many of the proteins identified in the proteome were found to be homologous to proteins with known activities affecting vasodilation and platelet aggregation. We show that vampire bat venom proteins possibly evade host immune response by the mutation of the surface chemistry through focal mutagenesis under the guidance of positive Darwinian selection. These results not only contribute to the body of knowledge regarding haematophagous venoms but also provide a rich resource for novel lead compounds for use in drug design and development.

BIOLOGICAL SIGNIFICANCE

These results have direct implications in understanding the molecular evolutionary history of vampire bat venom. The unusual peptides discovered reinforce the value of studying such neglected taxon for biodiscovery.

Salivary antigen-5/CAP family members are Cu2+-dependent antioxidant enzymes that scavenge O₂₋. and inhibit collagen-induced platelet aggregation and neutrophil oxidative burst

The function of the antigen-5/CAP family of proteins found in the salivary gland of bloodsucking animals has remained elusive for decades. Antigen-5 members from the hematophagous insects Dipetalogaster maxima (DMAV) and Triatoma infestans (TIAV) were expressed and discovered to attenuate platelet aggregation, ATP secretion, and thromboxane A2 generation by low doses of collagen (<1 μg/ml) but no other agonists. DMAV did not interact with collagen, glycoprotein VI, or integrin α2β1. This inhibitory profile resembles the effects of antioxidants Cu,Zn-superoxide dismutase (Cu,Zn-SOD) in platelet function. Accordingly, DMAV was found to inhibit cytochrome c reduction by O2[Symbol: see text] generated by the xanthine/xanthine oxidase, implying that it exhibits antioxidant activity. Moreover, our results demonstrate that DMAV blunts the luminescence signal of O2[Symbol: see text] generated by phorbol 12-myristate 13-acetate-stimulated neutrophils. Mechanistically, inductively coupled plasma mass spectrometry and fluorescence spectroscopy revealed that DMAV, like Cu,Zn-SOD, interacts with Cu(2+), which provides redox potential for catalytic removal of O2[Symbol: see text]. Notably, surface plasmon resonance experiments (BIAcore) determined that DMAV binds sulfated glycosaminoglycans (e.g. heparin, KD ~100 nmol/liter), as reported for extracellular SOD. Finally, fractions of the salivary gland of D. maxima with native DMAV contain Cu(2+) and display metal-dependent antioxidant properties. Antigen-5/CAP emerges as novel family of Cu(2+)-dependent antioxidant enzymes that inhibit neutrophil oxidative burst and negatively modulate platelet aggregation by a unique salivary mechanism.