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

Introduction

Genes of the major histocompatibility complex (MHC) play an essential role in the adaptive immune response of vertebrates and are crucial for the understanding of the influence of natural selection on genetic diversity in wild populations (Bernatchez and Landry, 2003). MHC-encoded cell surface glycoproteins present antigens of intracellular (MHC class I genes) and extracellular (MHC class II genes) origin to T cells that trigger the appropriate immune response (Klein, 1986; Janeway and Travers, 2002; Kumanovics et al., 2003). Class II-derived molecules are heterodimers of two associated polypeptide chains, an _α_- and a _β_-chain. The _β_-chain consists of an antigen-binding β _1_-domain, an immunoglobulin-like _β2_-domain and a transmembrane domain. Certain amino-acid positions of the β _1_-region, the so-called antigen-binding sites (ABS, encoded by exon 2), show high levels of genetic variation (Brown et al., 1993) with a higher rate of non-synonymous (_d_N, amino-acid altering) over synonymous (_d_S, silent) nucleotide substitutions (Hughes and Nei, 1988, 1989). According to neutrality theory (Kimura, 1977), the rate of synonymous nucleotide substitution is predicted to be larger than the rate of non-synonymous substitutions because a change of the amino-acid composition is more likely to be deleterious. The pattern of an elevated rate of non-synonymous substitutions at ABS is considered as clear evidence for ‘positive Darwinian selection’ (Hughes and Nei 1988, 1989) shaping genetic variation (reviewed in Hughes, 2007). Others used the term ‘historical positive selection’ (Hedrick, 1999) or ‘positive selection over evolutionary time scale’ (Bernatchez and Landry, 2003) to describe this observation. This polymorphism enables the immune system to recognise an extensive range of pathogens and is therefore crucial for the immunological fitness of individuals and, thus, animal populations (Edwards and Potts, 1996; Bernatchez and Landry, 2003). The high levels of polymorphism at MHC class II loci found in most vertebrate species are thought to be maintained by different forms of balancing selection, including heterozygosity advantage (Doherty and Zinkernagel, 1975) and rare allele advantage (also called negative frequency-dependent selection; Takahata and Nei, 1990; Slade and McCallum, 1992) as well as selection that varies in space and time (Hill, 1991; Hedrick, 2002). Distinguishing between the different forms of balancing selection in natural populations is difficult, as they are not mutually exclusive and may operate together with other selective and neutral forces (Sommer, 2005; Piertney and Oliver, 2006; Spurgin and Richardson, 2010).

Genes of the MHC are known to alter the individual body odour (for example, Penn and Potts, 1998; Milinski, 2006; Kwak et al., 2009) and may therefore be involved in mate choice decisions (for example, Wedekind et al., 1995; Chaix et al., 2008; Ilmonen et al., 2009), thus promoting genetic structure within and among social groups (Chesser, 1990; Kundu and Faulkes, 2004; Matocq and Lacy, 2004; Cutrera and Lacey, 2006). Sexual selection processes may not only be supported in direct mate choice situations but also through postcopulatory mechanisms including maternal–fetal interactions (for example, Wedekind et al., 1996; Ober and Van der Ven, 1997). Several studies have already highlighted that reproduction among MHC-dissimilar mates is favoured (for example, primates: Schwensow et al., 2008; and humans: Chaix et al., 2008; Ilmonen et al., 2009); mostly because a genetically diverse offspring genotype is advantageous in the defence against pathogens (Parham and Otha, 1996; Penn and Potts, 1999). Yet, the debate is still controversial. In some populations, MHC-based mate choice was not affirmed (humans: Hedrick and Black, 1997; Chaix et al., 2008; and ruminants: Patterson and Pemberton, 1997). Some studies have suggested that an intermediate rather than the highest level of MHC diversity is optimal (Wegner et al., 2003; Woelfig et al., 2009).

Bats are not only the second largest mammalian order but also the most gregarious of all mammals. Some bat colonies harbour several million individuals. Thus, social communication among bats can be complex, involving not only visual and acoustical but also olfactory signals (McCracken and Wilkinson, 2000; Altringham and Fenton, 2003). Many bat species produce distinct odours from a variety of glands (for example, Quay, 1970; Dapson et al., 1977; Voigt and von Helversen, 1999; Scully et al., 2000; Caspers et al., 2009). Pilot studies have shown that bats use odours for kin and individual recognition (Gustin and McCracken, 1987; Bouchard, 2001; Safi and Kerth, 2003) and during male–female interactions (Voigt and von Helversen, 1999; Bouchard, 2001). Recently, it has been discovered that bats are reservoir hosts and potentially important vectors of many infectious diseases (Calisher et al., 2006; Wong et al., 2007; Field, 2009; Wibbelt et al., 2010). Given the combination of these facts, it becomes apparent that the study of bat immunogenetics is of fundamental relevance not only for understanding the mechanisms underlying the zoonotic potential of bats but also to shed light on the link between MHC, pathogen resistance and olfaction-based mate choice in mammals.

In the present study, we examined the MHC class II DRB region of a neotropical bat species, the lesser bulldog bat, Noctilio albiventris, in Panama. N. albiventris lives in social groups of so far unknown relatedness. Besides acoustical signals, group members also communicate among each other via volatile compounds from the so-called oily spots in the subaxillary region. In addition, males possess inguinal pockets next to the scrotum that provide a distinct odour, particularly during the reproductive period (Studier and Lavoie, 1984; Brook and Decker, 1996). In summary, N. albiventris harbours morphologically distinct scent-producing organs, suggesting that volatiles are of key importance not only for individual communication but also for mate choice.

Our specific aims were: (1) to design reliable primers to characterise MHC class II DRB exon 2 diversity in N. albiventris, using a genome walking approach that enables the amplification of specific DNA fragments in situations where the sequence of only one primer is known (modified from Ko et al., 2003, 2) to examine evidence for selection acting on MHC in N. albiventris, that is, a higher rate of non-synonymous over synonymous substitutions in ABS would indicate positive selection acting above the evolutionary history of the species and (3) to investigate whether subpopulations, roosting colonies and the sexes differ in their genetic population structure. Differences in the MHC constitution of subpopulations and roosting colonies could indicate local adaptation as a result of recent selection events, whereas differences between sexes could indicate MHC-related sexual selection mechanisms.

Materials and methods

Study site and sampling

Our main study site was the village Gamboa (09.07° N and 079.41° W) in Panama. Bats were caught with mist nets (see Dechmann et al., 2009) as they emerged at dusk from their roost during 2006–2008. Individuals caught from the same daytime roost were assumed to be from the same colony. A second study site was Barro Colorado Island (BCI; 09.10° N and 079.51° W) situated 15 km away from Gamboa, where animals foraging above water were caught in mist nets set up along the boat docks of a marina (Figure 1). Bats from the two sites were considered separate subpopulations as their foraging areas did not overlap (Dechmann et al., 2009). Bats were sexed, aged and body measurements taken. From all 215 bats, we collected a 4-mm skin sample from the wing membrane using a sterile biopsy punch (Worthington Wilmer and Barrat, 1996). Skin samples were stored in 96% ethanol until DNA isolation.

Figure 1

figure 1

Map of the study sites Gamboa and BCI in Panama (modified from the Smithonian Tropical Research Institute mapserver).

Full size image

For optimal primer design and to prove expression of MHC alleles, a liver sample was collected from an euthanised male N. albiventris. The liver sample was preserved in RNA later (Sigma-Aldrich, Steinheim, Germany) and stored at −20 °C until subsequent analysis.

Genomic DNA and RNA isolation and complementary DNA library construction

Genomic DNA (gDNA) was extracted using the DNeasy Tissue Kit (Qiagen, Hilden, Germany) following the manufacturer's protocol. RNA was isolated from the liver using the RNeasy Kit (Qiagen) that includes a DNase treatment according to the manufacturer's instructions. The quantity was determined by measuring the absorbance at 260 nm and purity (that is presence of contaminants) was assessed by the ratio at 260/280 nm, which should be ∼1.8 for DNA and ∼2.0 for RNA (Nanotrop, Peqlab, Erlangen, Germany). Up to 5 μg of total RNA were used to construct a first strand complementary DNA (cDNA) library. For reverse transcription, we included 5 μg Oligo(dT)12–18 Primer (Invitrogen, Karlsruhe, Germany) and 200 U of Revert Aid H-Minus M-MulV reverse transcriptase (Fermentas, St Leon-Rot, Germany) together with 20 U of a ribonuclease inhibitor (RiboLock, Fermentas) as instructed in the manufacturer's protocol. cDNA was checked spectrophotometrically (Nanotrop) and by electrophoresis on a 1.5% agarose gel. Reproducibility was ensured by two independent repeats of the whole procedure.

Primer design using cDNA and vectorette PCR

An MHC class II DRB cDNA alignment of different mammalian species was constructed by a BLAST search (http://www.ncbi.nlm.nih.gov) of the GenBank database. Primers complementary to conserved parts of exon 1 (Ex1) and exon 3 (Ex3) were designed (Figure 2). They were used together with primer L729 situated in exon 4 (Bowen et al., 2004, Figure 2) in different combinations to amplify DRB sequences from the cDNA of N. albiventris. PCR amplification was carried out in a total reaction volume of 20 μl, including 20 ng cDNA, 1 U of proofreading polymerase (HotStar HiFidelity Polymerase, Qiagen), 1 × HotStar HiFidelity PCR buffer, 0.3 mM deoxynucleotid triphosphates (dNTPs) and 0.4 μM of each primer (Sigma-Aldrich, Steinheim, Germany). Thermocycling started with an initial denaturation step at 96 °C for 10 min, followed by 35 cycles of denaturation at 96 °C for 45 s, annealing at 58 °C for 45 s and elongation at 72 °C for 2 min. An additional 15 min elongation step followed the last cycle. PCR products were purified (Cycle pure Kit, Peqlab) and cloned using the pCR4-TOPO TA cloning Kit (Invitrogen). At least 23 recombinant clones per primer combination were amplified using the vector primers T7for and M13rev, sequenced using BigDye Terminator v.3.1 Cycle Sequencing Kit (Applied Biosystems, Foster City, CA, USA) and analysed on ABI PRISM 310 Automated Genetic Analyzer (Applied Biosystems). All primer combinations amplified the same two MHC class II DRB alleles in N. albiventris. These Noal sequences were used to design the species-specific primers Ex2a and Ex2b (Figure 2) binding to conserved sites of exon 2 (inferred from the established cDNA alignment and available sequences of the bat Saccopteryx bilineata (GenBank accession numbers: EF533888-EF533900)), which were used for the subsequent first vectorette PCR assay.

Figure 2

figure 2

Positions and sequences of PCR primers to amplify the indicated fragments of the MHC class II DRB gene in N. albiventris based on gDNA (a) and cDNA (b). Boxes symbolise exons, dark lines mark introns and dotted lines mark not sequenced parts of the introns. Open circles indicate positions of ligated vectorettes (see Materials and Methods for details). The primer _DRB_-L729 is taken from Bowen et al. (2004).

Full size image

Vectorette PCR is a method to amplify DNA fragments of interest, where the sequence information is only available on one side, and is described in detail by Ko et al. (2003). We constructed vectorette libraries with gDNA of five individuals. Up to 5 μg of gDNA, because of available template of each animal, were digested with restriction enzymes _Eco_RI, _Xap_I, Fsp, and _Csp_I (Fast Digest, Fermentas). In total, 1–5 μg gDNA were digested by 1–5 U restriction enzyme in a total volume of 50 or 100 μl depending on the amount of gDNA. Double-stranded vectorettes (1 μM) consisting of vect53 and vect57TTAA (_Eco_RI and _Xap_I) or vect53 and vect57AT (Fsp and _Csp_I) (Ko et al., 2003) were ligated to the sticky ends of the digested gDNA using 2 U of T4DNA ligase per 1 μg digested gDNA (Rapid DNA Ligation Kit, Fermentas).

Vectorette PCR assays consisted of two PCRs. The first PCR was followed by a second nested PCR to avoid false positives. A step-down scheme was always applied, using HotStar Taq Master Mix (Qiagen) according to the user's manual in a total reaction volume of 25 μl with the following conditions: initial activation of HotStar Taq and denaturation at 95 °C for 14 min, followed by 5 cycles of denaturation at 95 °C for 60 s, annealing at 67 °C for 60 s and elongation at 72 °C for 2 min, five cycles with annealing at 63 °C, followed by 15 cycles denaturation at 95 °C for 45 s, annealing at 60 °C for 45 s and extension at 72 °C for 120 s, 15 cycles with annealing at 57 °C and a final extension at 72 °C for 10 min. Amplification products were checked on a 1.5% agarose gel and purified with cycle Pure Kit (Peqlab).

The first vectorette PCR assay was performed with the primer Ex2a and the vectorette primer C20 (Ko et al., 2003), and the nested vectorette PCR included the nested primer Ex2b in combination with the nested vectorette primer B21 (Ko et al., 2003). The resulting intron 2 sequences allowed the design of _Noctilio_-specific intron 2 primers (Int2a and Int2b) (Figure 2). These were used to amplify intron 1 sequences in another vectorette PCR assay. Finally, according to the intron 1 sequences, the primers Int1a and Int1b were developed. Int1a and Int2a were used to amplify the whole 270 bp MHC DRB class II exon 2 and partial introns (45 bp in intron1, 48 bp in intron2, Figure 2) and were applied to genotype all 215 N. albiventris individuals.

Genotyping and allele identification

Amplification was performed in 25 μl volumes containing 0.4 μM of each primer (Sigma-Aldrich), 1 × reaction buffer (10 mM Tris HCl, 50 mM KCl, 0.1% Triton x-100, 0.2 mg ml−1 bovine serum albumin, MP Biomedicals, Heidelberg, Germany). 0.2 mM of each dNTP and 1 U Taq polymerase (MP Biomedicals). Step-down PCR was performed as follows: initial incubation at 95 °C for 5 min, followed by 5 cycles of denaturation at 95 °C for 10 s, annealing at 67 °C for 10 s and elongation at 72 °C for 30 s, five cycles with annealing at 64 °C and 25 cycles denaturation at 95 °C for 10 s, annealing at 60 °C for 10 s and elongation at 72 °C for 30 s, final extension was at 72 °C for 3 min. Amplicons were genotyped by single strand confirmation polymorphism on a polyacrylamide gel as described in Schad et al. (2004). Allele identification was done by excising the single strands from the gel matrix and diluting them in distilled water. Thereafter a reamplification of diluted single strands with primers Int1a and Int2a was performed before sequence analyses as described above. An autonomous amplicon with primer Int1b and Int2b of each individual was directly sequenced to confirm the individual single strand confirmation polymorphism pattern.

Statistical analyses

We edited and aligned nucleotide sequences manually using MEGA 4.0 (Tamura et al., 2007). We also used this programme for calculating the _P_-distance of amino-acid sequences as a measurement of functional MHC class II DRB divergence (Nei and Kumar, 2000) and for analysing the relative rates of non-synonymous (_d_N) and synonymous (_d_S) base pair substitutions according to Nei and Gojobori (1986), applying the Jukes–Cantor correction for multiple hits (Jukes and Cantor, 1969). Calculations were applied for all sites and separately for putative ABS and non-ABS, assuming functional congruence to human ABS of the HLA-DR1 (DRA/DRB1*0101) molecule (Brown et al., 1993). The _d_N/_d_S ratios of all sites as well as for ABS and non-ABS separately were compared with an implemented _Z_-test (Nei and Kumar, 2000) to test for positive selection.

We identified species-specific positively selected sites (PSS) with maximum-likelihood analysis using CODEML (included in PAML version 3.15 software package; Yang, 1997; Yang, 2007) and compared these PSS with the human ABS. First, we fitted models with different assumptions of selection patterns to the sequence data. We used the models M7 (β) and M8 (β and ω) as described in Yang et al. (2007). M7 served as a null model, where the ω-ratio varies according to the β- distribution and does not allow positive selected sites (0<ω<1). M8 adds a class of sites to account for the possible occurrence of PSS (ω>1). The models were compared using a likelihood ratio test by calculating the likelihood difference 2Δ_l_=2(_l_1−_l_0) and then compared with a _χ_2-distribution, with the degrees of freedom equal to the difference in the number of estimated parameters (Yang and Bielawski, 2000). In the next step, after likelihood ratio test provided evidence for positive selection, the Bayes empirical Bayes method (Yang et al., 2005) integrated in CODEML was used to identify the sites under positive selection with the cutoff posterior probability set at _P_b=95%.

Allele frequencies, observed and expected heterozygosity, and deviation from Hardy–Weinberg expectations were estimated using the software ARLEQUIN version 3.0 (Excoffier et al., 2005). Pairwise FST based on haplotype frequencies were calculated to infer population subdivision (10 000 permutations, Wright, 1951, 1965). Allelic richness (R) as a measure of the number of alleles independent of sample size was estimated using the rarefaction method as implemented in FSTAT version2.9.3 ([Goudet, 2001](/articles/hdy2010173#ref-CR32 "Goudet (2001). FSTAT version, a program to estimate and test gene diversities and fixation. Available from http://www.unil.ch/izea/softwares/fstat.html

            .")).

The _χ_2-tests were used to compare the number of heterozygote individuals between the sexes. Differences in the mean individual amino-acid distance of males and females were investigated by analysis of variance (ANOVA). Calculations were always two-tailed with significance level at α=0.05 and performed using SPSS version 16.0 (SPSS Inc., Chicago, IL, USA). Bonferroni corrected significant levels were applied for multiple comparisons (Rice, 1989).

Results

Capture and sampling

In total, we caught 215 individuals (91 males and 123 females), 185 in Gamboa (59 males and 108 females) and 29 in BCI (22 males and 7 females). All of them were genotyped for the MHC class II DRB gene. In Gamboa, we caught bats from seven roosts. We included only colonies with more than 10 individuals in statistical analyses (colony A: _N_=52, 22 males and 30 females; colony B: _N_=74, 20 males and 54 females; colony C: _N_=27, 13 males and 14 females; and colony D: _N_=14, 4 males and 10 females).

MHC class II DRB diversity

Using cDNA of one male N. albiventris, all primer combinations amplified the same two MHC class II DRB alleles differing only in fragment length (628, 527 and 376 bp, respectively, Figure 2). We designed species-specific intron primers (Int1a and Int2a) for amplification of the whole MHC class II DRB exon 2 (270 bp) of N. albiventris flanked by short intron sequences with a total length of 363 bp (Figure 2) by means of the vectorette PCR approach based on gDNA. No indels or stop codons were found and no more than two alleles per individual were amplified, suggesting that a single DRB locus was expressed. We named this locus _MhcNoal_-DRB according to the established MHC nomenclature (Klein et al., 1990; Ellis et al., 2006). The nucleotide and deduced amino-acid sequence showed high homology with human (84%), canine (83%), equine (83%) and other mammalian class II DRB genes.

We detected 18 different alleles of _Noal_-DRB in the 215 individuals of N. albiventris (GenBank accession numbers: HM347941-HM347958). In the nucleotide sequences, we observed 71 (26.3%) variable positions and the alleles differed by 2 to 40 (average 26.3±2.9) nucleotide positions. All alleles had a unique amino-acid sequence, whereas 38 (43%) out of 89 amino acids were polymorphic. They differed by 1 to 26 (average 17.6±2.5) amino-acid positions.

The most common allele Noal-DRB*10 occurred at a frequency of 0.255 followed by Noal-DRB*02 (0.175) and Noal-DRB*04 (0.117). The remaining alleles occurred in 2 to 36 individuals with a frequency ranging between 0.005 and 0.089 (Table 1). Observed heterozygosity (0.902) was higher than expected (0.871) and did not deviate from Hardy–Weinberg expectations (Table 2). The individual Noal-DRB exon 2 distance of an individual based on amino-acid sequence ranged from 0 (homozygote) to 0.287 with an average of 0.177±0.078 (Table 2).

Table 1 Allele frequencies of MHC class II DRB alleles in the whole population and in subpopulations, roosting colonies and sexes

Full size table

Table 2 MHC class II DRB exon 2 variability in N. albiventris in the whole population and in subpopulations, roosting colonies and sexes

Full size table

Evidence for historical selection

Two approaches were used to test for historical positive selection acting on the examined exon 2 of Noal-DRB locus. First, the averaged rates of non-synonymous (_d_N) and synonymous (_d_S) base pair substitutions of all sequences were calculated for all sites and separately for putative ABS and non-ABS, assuming functional congruence to human ABS of the HLA-DR1 molecule (Brown et al., 1993). Non-synonymous substitutions occurred at a significantly higher rate than synonymous ones (_d_N/_d_S=2.90, _Z_-test, _P_=0.004) especially in the regions that code for ABS (_d_N/_d_S=3.14, _Z_-test, _P_=0.002; Table 3). Second, we observed a significantly higher log likelihood estimate (2Δ_l_=59.67, d.f.=2, P<0.0001) for model M8 (positive selection) than for its corresponding null model M7 (no positive selection). The Bayes empirical Bayes approach under model M8 inferred 17 significant sites (PSS) to be under positive selection with the cutoff posterior probability set at 95%. In total, 14 of these sites were congruent with predicted ABS of the human HLA-DR1 _β_-chain (Brown et al., 1993). The other three sites were located in close proximity within one to two amino-acid positions to the human ABS (Table 4, Figure 3). In the Noal-DRB sequences, all but one (61 W) of certain amino-acid positions, which are conserved in human HLA-DR1 molecules (Brown et al., 1993; Stern et al., 1994), were also conserved presenting the identical amino acids. At position 61, tryptophan was replaced by leucine in six Noal alleles (33.3%) and identified as PSS (Table 4, Figure 3).

Table 3 Non-synonymous (_d_N) and synonymous (_d_S) substitutions (±s.e.) as well as their ratio in ABS and non-ABS assuming concordance with the human HLA-DR1 molecule (Brown et al., 1993)

Full size table

Table 4 Identification of species-specific PSS by likelihood analysis in MHC class II DRB exon 2 amino-acid sequences of N. albiventris

Full size table

Figure 3

figure 3

Amino-acid sequence variation of 18 MHC class II DRB exon 2 alleles of N. albiventris and the human _HLA_-DRB1*0101 sequence. Dots mark identity with the top sequence. Numbers indicate the amino-acid positions of the _β1_-domain, * signify the amino-acid positions of ABS and conserved sites (ConsSites) of the human HLA-DR1 _β_-chain (Brown et al., 1993; Stern et al., 1994). PSS indicates species-specific positive selected sites identified by CODEML (Yang et al., 2005).

Full size image

Population structure and MHC

We observed no genetic differentiation between the two subpopulations of Gamboa and BCI (FST=0.003, _P_=0.23). Allele frequencies of the two subpopulations are shown in Table 1. Numbers of alleles were similar after correction for differences in sample size. Observed heterozygosity was high in both subpopulations without a significant deviation from Hardy–Weinberg expectations (Table 2).

All four roosting colonies with more than 10 individuals captured in the village of Gamboa showed high levels of heterozygosity and allelic richness (R) ranging from 9 to 11.1. Observed heterozygosity exceeded the expected value in all but one colony (colony B) and did not deviate significantly from Hardy–Weinberg equilibrium (Table 2). Pairwise FST statistics showed a slight differentiation between colonies A, C and D but this significance was lost after Bonferroni correction (α′ ⩽0.008; Table 5).

Table 5 Pairwise differentiation between roosting colonies using conventional F-statistic based on haplotype frequencies (Wright, 1965)

Full size table

We found no significant differentiation within males and females when comparing colonies (males: range of FST: 0.001–0.026, _P_-values not significant; females: FST: <0.001–0.039, _P_-values not significant; Bonferroni corrected significance level α′⩽0.008). We found no significant differentiation between males and females within colonies either (FSTA=0.009, _P_=0.87; FSTB=0.001, _P_=0.44; FSTC=0.007, _P_=0.54; FSTD=0.030, _P_=0.82; and FSTAll=0.002, _P_=0.87).

However, males and females differed significantly in their individual amino-acid distance (ANOVA: F=8.48, _P_=0.004, d.f.=1, Figure 4, Table 2). Furthermore, the observed heterozygosity of males was significantly higher than in females (overall: _χ_2=7.73, d.f.=1, _P_=0.005). The analyses were not significant when colonies were analysed separately (A: _χ_2=0.81, _P_=0.36; B: _χ_2=2.54, _P_=0.11; C: _χ_2=3.13, _P_=0.07; and D: _χ_2=0.93, _P_=0.33). In males, the observed heterozygosity exceeded the expected value in the overall sample as well as in all colonies analysed separately but did not deviate significantly from Hardy–Weinberg equilibrium (overall: Hobs=0.97, Hex_P_=0.89; A: Hobs=1.00, Hex_P_=0.91; B: Hobs=0.96, Hex_P_=0.89; C: Hobs=1.00, Hex_P_=0.88; D: Hobs=1.00, Hex_P_=0.89; and BCI: Hobs=1.00, Hex_P_=0.89). In females, the observed heterozygosity was always lower than in males but close to the expected value in the overall sample as well as in two out of four colonies (overall: Hobs=0.85, Hex_P_=0.86; A: Hobs=0.96, Hex_P_=0.90; B: Hobs=0.83, Hex_P_=0.86; C: Hobs=0.78, Hex_P_=0.79; D: Hobs=0.75, Hex_P_=0.75; and BCI: Hobs=1.00, Hex_P_=0.92; Table 2).

Figure 4

figure 4

Mean (±s.e.) amino-acid distance between individual MHC class II DRB exon 2 alleles of N. albiventris in males and females. ANOVA: F=8.48, d.f.=1, **_P_=0.004; N, sample size.

Full size image

Discussion

To date, not much is known about the genetic structure and polymorphism of the MHC in bats (Meyer and Brunner, 2007) mainly because of a lack of sequence data for primer design. This is astonishing given the size of the order, the importance of olfactory signals for social communication in bats and their zoonotic relevance as reservoir hosts for many pathogens. Therefore, our primary aim was to determine MHC class II DRB exon 2 diversity in the lesser bulldog bat, N. albiventris. Further, we tested for evidence of positive selection acting above the evolutionary history of the species. And third, we were interested whether local adaptation and sexual selection may shape the contemporary genetic structure of N. albiventris in a Central American population.

MHC class II DRB diversity

Well-designed primers are essential for population genetic analyses. In this study, cDNA analyses followed by a vectorette PCR approach offered the possibility of intron mapping. Species-specific primers binding to conserved segments of the flanking introns of N. albiventris were designed to amplify the whole DRB exon 2. The reliability of the resulting species-specific sequence patterns of the locus of interest is very high because the incidence of non-amplifying alleles can be neglected. This might turn into a problem when exon primers are used to amplify parts of polymorphic genes. The occurrence of non-amplifying alleles would preclude the use of data for many population genetic purposes because they can seriously bias population genetic analyses (Dakin and Avise, 2004; Cummings et al., 2010).

We found evidence for a single expressed MHC class II DRB locus in N. albiventris with moderate allelic variability of 18 alleles detected in 215 individuals. This allelic variability is within the range of MHC class II DRB polymorphism of other mammalian species (for example, Sommer, 2005). The only other DRB gene studied in a bat, Saccopteryx bilineata, revealed a rather low allelic variability (11 alleles in 85 individuals; Meyer and Brunner, 2007). However, this might be an underestimate of the actual variability because preliminary studies on RNA and DNA with species-specific developed primers revealed evidence for at least five DRB loci in Saccopteryx bilineata (Schad et al., unpublished data). In general, the occurrence of different MHC class II genes (DP, DO, DM, DQ and DR) is conserved in mammals. However, the number of functional _α_- and _β_-genes is highly variable because of species-specific local duplication events. Yet, it can even vary between individuals of the same species (see Kumanovics et al., 2003). The DRB is the most widely studied and usually the most diverse class II gene, not only with respect to high allelic richness but also in terms of gene duplications (for example, Doxiadis et al., 2000; Bowen et al., 2004; Babik et al., 2005; Schwensow et al., 2007). We cannot rule out that other MHC class II genes may provide a higher variability in bats as it is described, for instance, in the cetacean MHC class II DQB gene (Baker et al., 2006). As a future task, it will be necessary to investigate other loci to evaluate MHC diversity in bats more comprehensively.

Evidence for historical selection

We found a higher rate of non-synonymous versus synonymous nucleotide substitutions over the entire sequence and especially in the putative ABS but not in non-ABS (Brown et al., 1993). Such elevated _d_N/_d_S ratios are a widely accepted sign for historical positive selection in polymorphic MHC genes (Hughes and Nei 1988, 1989; Nei and Kumar, 2000; Nielsen, 2001; Hughes, 2007). They indicate positive selection on ABS acting above the evolutionary history of the species (Hedrick, 1999; Bernatchez and Landry, 2003; Sommer, 2005). We subsequently confirmed these results in a species-specific analysis where different codon evolution models were compared by maximum-likelihood analyses (Yang and Bielawski, 2000; Yang, 2007). This method has been suggested as a powerful tool with a high likelihood of detecting effects of positive selection (Wong et al., 2004). Our data fitted best with the model incorporating positive selection and 17 species-specific PSS were identified by the Bayes empirical Bayes analysis. A total of 14 were congruent with human ABS of the HLA-DR1 _β_-chain. Three PSS were situated outside the human ABS and some of the human ABS were not identified as PSS in Noal-DRB. Comparable results have been reported from other species (Kundu and Faulkes, 2004; Schwensow et al., 2007; Babik et al., 2008; Meyer-Lucht et al., 2008). In general, high congruence of PSS with human ABS is assumed to demonstrate homologous functionality of the molecule. Contrarily, human ABS that are not identified as PSS might be not involved in the antigen recognition and binding in the respective species. Furthermore, in most species investigated so far additional PSS have been reported. These findings suggest species-specific selection pressures acting on MHC genes because of a different pathogen exposure. In addition, certain amino-acid residuals of HLA-DR1 molecules are highly conserved and involved in universal hydrogen bond of antigens (61 W, 81 H and 82 N) or are responsible for the stability of DR1 heterodimers in building salt bridges between the dimers (52 E and 55 R) (Brown et al., 1993; Stern et al., 1994). Also in the Noal-DRB sequences, all of these positions, except position 61, were conserved indicating similar conserved functionality of the molecules.

Population structure and MHC

Mating behaviour as a correlate and driver of social structure has been suggested in addition to pathogens as another main subject of selection on MHC loci in natural populations (Hambuch and Lacey, 2002; Kundu and Faulkes, 2004; Cutrera and Lacey, 2006). Animals form social groups and colonies in response to cooperative interactions as well as mating tactics. Furthermore, patterns of genetic subdivision are also shaped by the extent and nature of philopatric behaviour (Travis et al., 1995; Sommer et al., 2002; Solomon, 2003; Cutrera and Lacey, 2006). We examined the genetic structure of three different levels of social formations or units (subpopulation, roosting colonies and sexes) in our population of N. albiventris. We found no genetic differentiation between the two subpopulations (Gamboa and BCI, separated by 15 km), indicating the presence of gene flow at a larger spatial scale. Both subpopulations showed similar levels of polymorphism at the Noal-DRB locus. Observed heterozygosity was high in the two subpopulations and the allele frequencies and the number of alleles were similar after correcting for differences in sample sizes. This implies an equivalent selection pressure maintaining diversity at the Noal-DRB locus in both subpopulations in the recent past.

To investigate social structure based on roosting habits, we compared four colonies in the village of Gamboa, all located in the range of 1.5 km2. All colonies showed similar levels of genetic variation (heterozygosity, allelic richness and allele frequencies). F-statistics revealed only limited effects of subdivisions. The colony D showed a slight differentiation compared with the colonies A and C (not significant after Bonferroni correction), which might be rather the result of missing rare alleles because of small sample size than an effect of population structure based on roosting habits. Dechmann et al. (2009) distinguished social groups of females (2–5 individuals) by the fact that they emerged simultaneously from a roost. We could not find genetic differences between male and female members of different colonies. The composition of colonies might vary in time, indicating a fission–fusion society rather than stable associations. Neutral markers would offer the possibility to gain insights into population dynamic processes like kin relationships of social groups of females and roosting colonies. Thus, further ecological studies on demographic structure as well as genetic analyses adding neutral markers will help to fully understand the social system of this species.

We tested for gender-specific differences in the MHC constitution to assess the occurrence of sexual selection. We found no population differentiation between males and females. However, in the overall sample, males showed a significantly higher heterozygosity rate and also a higher individual amino-acid distance than females. In males, the observed heterozygosity exceeded the expected value in the overall sample as well as in all within-colony comparisons. In females, the observed heterozygosities were always lower than in males but were almost identical with the expected ones in the overall sample and in two out of four colonies. As in both sexes observed and expected heterozygosities did not significantly deviate from Hardy–Weinberg expectations, it remains unclear whether the difference between sexes is because of a higher heterozygosity rate in males or a lower heterozygosity rate in females. An increased heterozygosity rate in males could indicate balancing selection in form of a heterozygote advantage and a deficit in heterozygote females could suggest a reduced selection pressure to maintain diversity in females. It might be that the selection intensity in the investigated Noctilio population is not strong enough to detect significant deviations from Hardy–Weinberg expectations. It is well known that levels of allelic diversity in relation to sample size have an effect on the statistical power to detect significant deviations from Hardy–Weinberg expectations (Seddon and Ellegren, 2004). We did not find small-scale population structure that would offer a likely explanation for this gender-specific difference. Male-based long-distance dispersal is unlikely as it would promote the occurrence of new alleles in the male population, raising the overall polymorphism at the MHC locus that we did not observe. A methodological error because of DNA quality differences between the sexes seems also to be unlikely because collection and treatment of samples have been the same throughout the study. Although direct female choice for heterozygote males would lead to heterozygote offspring in general including daughters, it cannot explain the sex-specific bias in heterozygosity. Currently, the differences in heterozygosity in males and females is most likely due to MHC-mediated postcopulatory mechanisms (for example, caused by myotic drive, gametic selection and maternal–fetal interactions) or by sex-specific survival differences that result in an increased rate of MHC heterozygote males or a deficit in heterozygote females.

The possibility of the existence of postcopulatory mechanisms resulting in sex-specific differences in MHC heterozygosity was reported only in a few studies so far. Dorak et al. (2002) found an increased heterozygosity for MHC class II DRB lineages in newborn male babies and suggested that negative selection of homozygotes might be restricted to male offspring only. A deficit in MHC homozygosity in newborn males was also observed in mice (Hamilton and Hellstrom 1978) and rats (Palm 1969, 1970, 1974). Some studies in humans investigated the compatibility at different HLA loci between mothers and infants and observed differences in the sex ratio assuming a different fetal loss in males and females. The results are heterogeneous and differ between HLA loci (Ober et al., 1987; Astolfi et al., 1990, 1996). In all of these studies, the underlying mechanisms have not been investigated. The ongoing discussion has been reviewed by Fernandez et al. (1999) and Ziegler et al. (2005). The overall conclusion is that the MHC is critical for numerous aspects of mammalian reproduction concerning spermatogenesis (Ziegler et al., 2002, 2005), a sperm-selective egg-cumulus complex (Wedekind et al., 1996; Rülicke et al., 1998; Eisenach and Giojalas, 2006), and viability and development of the foetus (Gill, 1992; Ober et al., 1987; Astolfi et al., 1990; Wedekind et al., 1996; Ziegler et al., 2005). However, in all these processes the selective forces might act because of MHC-linked genes, for example, olfactory receptor genes, transcription factors and others rather than by the MHC itself (Ho et al., 1990; Gill, 1992; Ziegler et al., 2002, 2005; Eisenach and Giojalas, 2006). To the best of our knowledge, sex-specific differences in offspring survival based on MHC heterozygosity have not been reported yet. But they have been investigated using a microsatellite-based measure of outbreeding (mean d 2) in a few species. In the bat Rhinolophus ferrumequinum, outbreeding was positively associated with significantly increased survival in male offspring only (Rossiter et al., 2001). The authors postulate that outbreeding at the microsatellite markers reflects immunocompetence, which in turn influences mortality. They also suggest that characterisation of MHC loci may provide a suitable test for their hypothesis. Even so, in the red deer (Cervus elaphus) a contrary pattern was found as male offspring survival was negatively associated with outbreeding at microsatellite markers (Coulson et al., 1999). In mice, a reduced survivorship was reported for inbred adult males but not for females, most likely as a consequence of aggressive interactions of males in the defence of territories (Meagher et al., 2000). Ongoing studies might reveal the mechanisms that have contributed to the sex-biased diversity pattern in the investigated population of N. albiventris. In addition, neutral markers would provide more detailed information about social structure, dispersal behaviour and gene flow of both sexes.

In this first study on the MHC class II variability of the lesser bulldog bat, we detected high genetic variation and evidence for historical positive selection acting on a single expressed Noal-DRB locus. The polymorphism at the antigen-binding region of the molecule is considered as the precondition to cope with a variety of pathogens. No population differentiation between subpopulations, roosting colonies and sexes was observed, but males revealed a significantly higher heterozygosity rate and genetic variability in terms of the genetic distance between the individual MHC alleles than females. We are aware that at this state of the investigation no conclusions on the underlying mechanisms can be made but our data can provide the basis for further research on the role of the MHC constitution in host–pathogen interactions, individual body odours and sexual selection in a highly interesting bat species.

Accession codes

Accessions

GenBank/EMBL/DDBJ

References

Download references

Acknowledgements

We thank Silke L Voigt-Heucke, Antje Kretzschmar, Megan Lupek, Martin Wikelski and Apanie Wood for help with catching bats and sample collection. We thank Wieslaw Babik for initial advice on the vectorette PCR approach. We are grateful to Charles Gunn for linguistic advice. We thank three anonymous reviewers for helpful comments on former versions of the manuscript. The German Science Foundation funded C. Voigt, S. Speck and D. Dechmann (VO 890/11-1). All capture and handling of animals, as well as collection and export of samples were done in concordance with Panamanian laws under permits from the Panamanian authorities (ANAM) and under an IUCAC protocol. We are indebted to the Smithsonian Tropical Research Institute and their staff for help with logistics and interaction with the Panamanian authorities.

Author information

Authors and Affiliations

  1. Leibniz Institute for Zoo and Wildlife Research, Berlin, Germany
    J Schad, D K N Dechmann, C C Voigt & S Sommer
  2. Max Planck Institute of Ornithology, Vogelwarte Radolfzell, Radolfzell, Germany
    D K N Dechmann
  3. Smithsonian Tropical Research Institute, Balboa, Ancon, Republic of Panama
    D K N Dechmann

Authors

  1. J Schad
    You can also search for this author inPubMed Google Scholar
  2. D K N Dechmann
    You can also search for this author inPubMed Google Scholar
  3. C C Voigt
    You can also search for this author inPubMed Google Scholar
  4. S Sommer
    You can also search for this author inPubMed Google Scholar

Corresponding author

Correspondence toS Sommer.

Ethics declarations

Competing interests

The authors declare no conflict of interest.

Rights and permissions

About this article

Cite this article

Schad, J., Dechmann, D., Voigt, C. et al. MHC class II DRB diversity, selection pattern and population structure in a neotropical bat species, Noctilio albiventris.Heredity 107, 115–126 (2011). https://doi.org/10.1038/hdy.2010.173

Download citation

Keywords