Non-neutral evolution of the major histocompatibility complex class II gene DRB1 in the sac-winged bat Saccopteryx bilineata

Henipavirus Immune Evasion and Pathogenesis Mechanisms: Lessons Learnt from Natural Infection and Animal Models

Nipah henipavirus (NiV) and Hendra henipavirus (HeV) are zoonotic emerging paramyxoviruses causing severe disease outbreaks in humans and livestock, mostly in Australia, India, Malaysia, Singapore and Bangladesh. Both are bat-borne viruses and in humans, their mortality rates can reach 60% in the case of HeV and 92% for NiV, thus being two of the deadliest viruses known for humans. Several factors, including a large cellular tropism and a wide zoonotic potential, con-tribute to their high pathogenicity. This review provides an overview of HeV and NiV pathogenicity mechanisms and provides a summary of their interactions with the immune systems of their different host species, including their natural hosts bats, spillover-hosts pigs, horses, and humans, as well as in experimental animal models. A better understanding of the interactions between henipaviruses and their hosts could facilitate the development of new therapeutic strategies and vaccine measures against these re-emerging viruses.

Major histocompatibility complex variation is similar in little brown bats before and after white-nose syndrome outbreak

White-nose syndrome (WNS), caused by the fungal pathogen Pseudogymnoascus destructans (Pd), has driven alarming declines in North American hibernating bats, such as little brown bat (Myotis lucifugus). During hibernation, infected little brown bats are able to initiate anti-Pd immune responses, indicating pathogen-mediated selection on the major histocompatibility complex (MHC) genes. However, such immune responses may not be protective as they interrupt torpor, elevate energy costs, and potentially lead to higher mortality rates. To assess whether WNS drives selection on MHC genes, we compared the MHC DRB gene in little brown bats pre- (Wisconsin) and post- (Michigan, New York, Vermont, and Pennsylvania) WNS (detection spanning 2014-2015). We genotyped 131 individuals and found 45 nucleotide alleles (27 amino acid alleles) indicating a maximum of 3 loci (1-5 alleles per individual). We observed high allelic admixture and a lack of genetic differentiation both among sampling sites and between pre- and post-WNS populations, indicating no signal of selection on MHC genes. However, post-WNS populations exhibited decreased allelic richness, reflecting effects from bottleneck and drift following rapid population declines. We propose that mechanisms other than adaptive immunity are more likely driving current persistence of little brown bats in affected regions.

Can extreme MHC class I diversity be a feature of a wide geographic range? The example of Seba's short-tailed bat (Carollia perspicillata)

The major histocompatibility complex (MHC) is one of the most diverse genetic regions under pathogen-driven selection because of its central role in antigen binding and immunity. The highest MHC variability, both in terms of the number of individual alleles and gene copies, has so far been found in passerine birds; this is probably attributable to passerine adaptation to both a wide geographic range and a diverse array of habitats. If extraordinary high MHC variation and duplication rates are adaptive features under selection during the evolution of ecologically and taxonomically diverse species, then similarly diverse MHC architectures should be found in bats. Bats are an extremely species-rich mammalian group that is globally widely distributed. Many bat species roost in multitudinous groups and have high contact rates with pathogens, conspecifics, and allospecifics. We have characterized the MHC class I diversity in 116 Panamanian Seba's short-tailed bats (Carollia perspicillata), a widely distributed, generalist, neotropical species. We have detected a remarkable individual and population-level diversity of MHC class I genes, with between seven and 22 alleles and a unique genotype in each individual. This diversity is comparable with that reported in passerine birds and, in both taxonomic groups, further variability has evolved through length polymorphisms. Our findings support the hypothesis that, for species with a geographically broader range, high MHC class I variability is particularly adaptive. Investigation of the details of the underlying adaptive processes and the role of the high MHC diversity in pathogen resistance are important next steps for a better understanding of the role of bats in viral evolution and as carriers of several deadly zoonotic viruses.

Prelude to a panzootic: Gene flow and immunogenetic variation in northern little brown myotis vulnerable to bat white-nose syndrome

The fungus that causes bat white-nose syndrome (WNS) recently leaped from eastern North America to the Pacific Coast. The pathogen’s spread is associated with the genetic population structure of a host ( Myotis lucifugus). To understand the fine-scale neutral and immunogenetic variation among northern populations of M. lucifugus, we sampled 1142 individuals across the species’ northern range. We used genotypes at 11 microsatellite loci to reveal the genetic structure of, and directional gene flow among, populations to predict the likely future spread of the pathogen in the northwest and to estimate effective population size ( Ne). We also pyrosequenced the DRB1-like exon 2 of the class II major histocompatibility complex (MHC) in 160 individuals to explore immunogenetic selection by WNS. We identified three major neutral genetic clusters: Eastern, Montane Cordillera (and adjacent sampling areas), and Haida Gwaii, with admixture at intermediate areas and significant substructure west of the prairies. Estimates of Ne were unexpectedly low (289–16 000). Haida Gwaii may provide temporary refuge from WNS, but the western mountain ranges are not barriers to its dispersal in M. lucifugus and are unlikely to slow its spread. Our major histocompatibility complex (MHC) data suggest potential selection by WNS on the MHC, but gene duplication limited the immunogenetic analyses.

Insights into the ancestral organisation of the mammalian MHC class II region from the genome of the pteropid bat, Pteropus alecto

BACKGROUND

Bats are an extremely successful group of mammals and possess a variety of unique characteristics, including their ability to co-exist with a diverse range of pathogens. The major histocompatibility complex (MHC) is the most gene dense and polymorphic region of the genome and MHC class II (MHC-II) molecules play a vital role in the presentation of antigens derived from extracellular pathogens and activation of the adaptive immune response. Characterisation of the MHC-II region of bats is crucial for understanding the evolution of the MHC and of the role of pathogens in shaping the immune system.

RESULTS

Here we describe the relatively contracted MHC-II region of the Australian black flying-fox (Pteropus alecto), providing the first detailed insight into the MHC-II region of any species of bat. Twelve MHC-II genes, including one locus (DRB2) located outside the class II region, were identified on a single scaffold in the bat genome. The presence of a class II locus outside the MHC-II region is atypical and provides evidence for an ancient class II duplication block. Two non-classical loci, DO and DM and two classical, DQ and DR loci, were identified in P. alecto. A putative classical, DPB pseudogene was also identified. The bat's antigen processing cluster, though contracted, remains highly conserved, thus supporting its importance in antigen presentation and disease resistance.

CONCLUSIONS

This detailed characterisation of the bat MHC-II region helps to fill a phylogenetic gap in the evolution of the mammalian class II region and is a stepping stone towards better understanding of the immune responses in bats to viral, bacterial, fungal and parasitic infections.

Spatial pattern of genetic diversity and selection in the MHC class II DRB of three Neotropical bat species

Background

Although bats are natural reservoirs of many pathogens, few studies have been conducted on the genetic variation and detection of selection in major histocompatibility complex (MHC) genes. These genes are critical for resistance and susceptibility to diseases, and host–pathogen interactions are major determinants of their extensive polymorphism. Here we examined spatial patterns of diversity of the expressed MHC class II DRB gene of three sympatric Neotropical bats, Carollia perspicillata and Desmodus rotundus (Phyllostomidae), and Molossus molossus (Molossidae), all of which use the same environments (e.g., forests, edge habitats, urban areas). Comparison with neutral marker (mtDNA D-loop) diversity was performed at the same time.

Results

Twenty-three DRB alleles were identified in 19 C. perspicillata, 30 alleles in 35 D. rotundus and 20 alleles in 28 M. molossus. The occurrence of multiple DRB loci was found for the two Phyllostomidae species. The DRB polymorphism was high in all sampling sites and different signatures of positive selection were detected depending on the environment. The patterns of DRB diversity were similar to those of neutral markers for C. perspicillata and M. molossus. In contrast, these patterns were different for D. rotundus for which a geographical structure was highlighted. A heterozygote advantage was also identified for this species. No recombination or gene conversion event was found and phylogenetic relationships showed a trans-species mode of evolution in the Phyllostomids.

Conclusions

This study of MHC diversity demonstrated the strength of the environment and contrasting pathogen pressures in shaping DRB diversity. Differences between positively selected sites identified in bat species highlighted the potential role of gut microbiota in shaping immune responses. Furthermore, multiple geographic origins and/or population admixtures observed in C. perspicillata and M. molossus populations acted as an additional force in shaping DRB diversity. In contrast, DRB diversity of D. rotundus was shaped by environment rather than demographic history.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0802-1) contains supplementary material, which is available to authorized users.

Evolution and comparative analysis of the bat MHC-I region

Bats are natural hosts to numerous viruses and have ancient origins, having diverged from other eutherian mammals early in evolution. These characteristics place them in an important position to provide insights into the evolution of the mammalian immune system and antiviral immunity. We describe the first detailed partial map of a bat (Pteropus alecto) MHC-I region with comparative analysis of the MHC-I region and genes. The bat MHC-I region is highly condensed, yet relatively conserved in organisation, and is unusual in that MHC-I genes are present within only one of the three highly conserved class I duplication blocks. We hypothesise that MHC-I genes first originated in the β duplication block, and subsequently duplicated in a step-wise manner across the MHC-I region during mammalian evolution. Furthermore, bat MHC-I genes contain unique insertions within their peptide-binding grooves potentially affecting the peptide repertoire presented to T cells, which may have implications for the ability of bats to control infection without overt disease.

Antiviral immune responses of bats: a review

Despite being the second most species-rich and abundant group of mammals, bats are also among the least studied, with a particular paucity of information in the area of bat immunology. Although bats have a long history of association with rabies, the emergence and re-emergence of a number of viruses from bats that impact human and animal health has resulted in a resurgence of interest in bat immunology. Understanding how bats coexist with viruses in the absence of disease is essential if we are to begin to develop therapeutics to target viruses in humans and susceptible livestock and companion animals. Here, we review the current status of knowledge in the field of bat antiviral immunology including both adaptive and innate mechanisms of immune defence and highlight the need for further investigations in this area. Because data in this field are so limited, our discussion is based on both scientific discoveries and theoretical predictions. It is hoped that by provoking original, speculative or even controversial ideas or theories, this review may stimulate further research in this important field. Efforts to understand the immune systems of bats have been greatly facilitated in recent years by the availability of partial genome sequences from two species of bats, a megabat, Pteropus vampyrus, and a microbat, Myotis lucifugus, allowing the rapid identification of immune genes. Although bats appear to share most features of the immune system with other mammals, several studies have reported qualitative and quantitative differences in the immune responses of bats. These observations warrant further investigation to determine whether such differences are associated with the asymptomatic nature of viral infections in bats.

The immune gene repertoire of an important viral reservoir, the Australian black flying fox

Background

Bats are the natural reservoir host for a range of emerging and re-emerging viruses, including SARS-like coronaviruses, Ebola viruses, henipaviruses and Rabies viruses. However, the mechanisms responsible for the control of viral replication in bats are not understood and there is little information available on any aspect of antiviral immunity in bats. Massively parallel sequencing of the bat transcriptome provides the opportunity for rapid gene discovery. Although the genomes of one megabat and one microbat have now been sequenced to low coverage, no transcriptomic datasets have been reported from any bat species. In this study, we describe the immune transcriptome of the Australian flying fox, Pteropus alecto, providing an important resource for identification of genes involved in a range of activities including antiviral immunity.

Results

Towards understanding the adaptations that have allowed bats to coexist with viruses, we have de novo assembled transcriptome sequence from immune tissues and stimulated cells from P. alecto. We identified about 18,600 genes involved in a broad range of activities with the most highly expressed genes involved in cell growth and maintenance, enzyme activity, cellular components and metabolism and energy pathways. 3.5% of the bat transcribed genes corresponded to immune genes and a total of about 500 immune genes were identified, providing an overview of both innate and adaptive immunity. A small proportion of transcripts found no match with annotated sequences in any of the public databases and may represent bat-specific transcripts.

Conclusions

This study represents the first reported bat transcriptome dataset and provides a survey of expressed bat genes that complement existing bat genomic data. In addition, these data provide insight into genes relevant to the antiviral responses of bats, and form a basis for examining the roles of these molecules in immune response to viral infection.

Independent evolution of functional MHC class II DRB genes in New World bat species

Genes of the major histocompatibility complex (MHC) play a pivotal role in the vertebrate immune system and are attractive markers for functional, fitness-related, genetic variation. Although bats (Chiroptera) represent the second largest mammalian order and are prone to various emerging infectious diseases, little is known about MHC evolution in bats. In the present study, we examined expressed MHC class II DRB sequences (exons 1 to 4) of New World bat species, Saccopteryx bilineata, Carollia perspicillata, Noctilio albiventris and Noctilio leporinus (only exon 2). We found a wide range of copy number variation of DRB loci with one locus detected in the genus Noctilio and up to ten functional loci observed in S. bilineata. Sequence variation between alleles of the same taxa was high with evidence for positive selection. We found statistical support for recombination or gene conversion events among sequences within the same but not between bat species. Phylogenetic relationships among DRB alleles provided strong evidence for independent evolution of the functional MHC class II DRB genes in the three investigated species, either by recent gene duplication, or homogenization of duplicated loci by frequent gene conversion events. Phylogenetic analysis of all available chiropteran DRB exon 2 sequences confirmed their monophyletic origin within families, but revealed a possible trans-species mode of evolution pattern in congeneric bat species, e.g. within the genera Noctilio and Myotis. This is the first study investigating phylogenetic relationships of MHC genes within bats and therefore contributes to a better understanding of MHC evolution in one of the most dominant mammalian order.

Bats and lyssaviruses

Numerous bat species have been identified as important reservoirs of zoonotic viral pathogens. Rabies and rabies-related viruses constitute one of the most important viral zoonoses and pose a significant threat to public health across the globe. Whereas rabies virus (RABV) appears to be restricted to bats of the New World, related lyssavirus species have not been detected in the Americas and have only been detected in bat populations across Africa, Eurasia, and Australia. Currently, 11 distinct species of lyssavirus have been identified, 10 of which have been isolated from bat species and all of which appear to be able to cause encephalitis consistent with that seen with RABV infection of humans. In contrast, whereas lyssaviruses are apparently able to cause clinical disease in bats, it appears that these lyssaviruses may also be able to circulate within bat populations in the absence of clinical disease. This feature of these highly encephalitic viruses, alongside many other aspects of lyssavirus infection in bats, is poorly understood. Here, we review what is known of the complex relationship between bats and lyssaviruses, detailing both natural and experimental infections of these viruses in both chiropteran and nonchiropteran models. We also discuss potential mechanisms of virus excretion, transmission both to conspecifics and spill-over of virus into nonvolant species, and mechanisms of maintenance within bat populations. Importantly, we review the significance of neutralizing antibodies reported within bat populations and discuss the potential mechanisms by which highly neurovirulent viruses such as the lyssaviruses are able to infect bat species in the absence of clinical disease.

MHC class II DRB diversity, selection pattern and population structure in a neotropical bat species, Noctilio albiventris

Genes of the major histocompatibility complex (MHC) have a crucial role in the immune response of vertebrates, alter the individual odour and are involved in shaping mating preferences. Pathogen-mediated selection, sexual selection and maternal-fetal interactions have been proposed as the main drivers of frequently observed high levels of polymorphism in functionally important parts of the MHC. Bats constitute the second largest mammalian order and have recently emerged as important vectors of infectious diseases. In addition, Chiroptera are interesting study subjects in evolutionary ecology in the context of olfactory communication, mate choice and associated fitness benefits. Thus, it is surprising that they belong to the least studied mammalian taxa in terms of their MHC diversity. In this study, we investigated the variability in the functionally important MHC class II gene DRB, evidence for selection and population structure in the group-living lesser bulldog bat, Noctilio albiventris, in Panama. We found a single expressed, polymorphic Noal-DRB gene. The substitution pattern of the nucleotide sequences of the 18 detected alleles provided evidence for positive selection acting above the evolutionary history of the species in shaping MHC diversity. Roosting colonies were not genetically differentiated but females showed lower levels of heterozygosity than males, which might be a sign that the sexes differ in the selection pressures acting on the MHC. This study provides the prerequisites for further investigations of the role of the individual MHC constitution in parasite resistance, olfactory communication and mate choice in N. albiventris and other bats.

Natural selection coupled with intragenic recombination shapes diversity patterns in the major histocompatibility complex class II genes of the giant panda

Ample variations of the major histocompatibility complex (MHC) genes are essential for vertebrates to adapt to various environmental conditions. In this study, we investigated the genetic variations and evolutionary patterns of seven functional MHC class II genes (one DRA, two DRB, two DQA, and two DQB) of the giant panda. The results showed the presence of two monomorphic loci (DRA and DQB2) and five polymorphic loci with different numbers of alleles (seven at DRB1, six at DRB3, seven at DQA1, four at DQA2, six at DQB1). The presence of balancing selection in the giant panda was supported by the following pieces of evidence: (1) The observed heterozygosity was higher than expected. (2) Amino acid heterozygosity was significantly higher at antigen-binding sites (ABS) compared with non-ABS sequences. (3) The selection parameter omega (d(N)/d(S)) was significantly higher at ABS compared with non-ABS sequences. (4) Approximately 95.45% of the positively selected codons (P>0.95) were located at or adjacent to an ABS. Furthermore, this study showed that (1) The Qinling subspecies exhibited high omega values across each locus (all >1), supporting its extensive positive selection. (2) The Sichuan subspecies displayed small omega at DRB1 (omega<0.72) and DQA2 (omega<0.48), suggesting that these sites underwent strong purifying selection. (3) Intragenic recombination was detected in DRB1, DQA1, and DQB1. The molecular diversity in classic Aime-MHC class II genes implies that the giant panda had evolved relatively abundant variations in its adaptive immunity along the history of host-pathogen co-evolution. Collectively, these findings indicate that natural selection accompanied by recombination drives the contrasting diversity patterns of the MHC class II genes between the two studied subspecies of giant panda.

Class II DRB polymorphism and sequence diversity in two vesper bats in the genus Myotis

Almost no studies have been done with respect to major histocompatibility complex (MHC) polymorphism and sequence diversity in bats, although they account for one in five living mammalian species. We analysed MHC Class II DRB polymorphism and sequence diversity in two Mexican verpertilionid bat species, the widespread continental species Myotis velifer and the narrowly distributed (and endangered) island endemic Myotis vivesi. We find extensive DRB polymorphism in the widespread M. velifer, similar to that commonly reported in other mammals. The geographically restricted M. vivesi by contrast shows only very limited polymorphism. We conclude that M. vivesi has undergone a dramatic loss of MHC polymorphism. The significance of this inference in light of other information on population structure and genetic diversity in this species is discussed.

Fellatio by fruit bats prolongs copulation time

Oral sex is widely used in human foreplay, but rarely documented in other animals. Fellatio has been recorded in bonobos Pan paniscus, but even then functions largely as play behaviour among juvenile males. The short-nosed fruit bat Cynopterus sphinx exhibits resource defence polygyny and one sexually active male often roosts with groups of females in tents made from leaves. Female bats often lick their mate's penis during dorsoventral copulation. The female lowers her head to lick the shaft or the base of the male's penis but does not lick the glans penis which has already penetrated the vagina. Males never withdrew their penis when it was licked by the mating partner. A positive relationship exists between the length of time that the female licked the male's penis during copulation and the duration of copulation. Furthermore, mating pairs spent significantly more time in copulation if the female licked her mate's penis than if fellatio was absent. Males also show postcopulatory genital grooming after intromission. At present, we do not know why genital licking occurs, and we present four non-mutually exclusive hypotheses that may explain the function of fellatio in C. sphinx.