Introduction {#Sec1}
============

The Chinese alligator (*Alligator sinensis*) is a very ancient species that had once been pushed to the fringe of extinction by massive hunting, territory deprivation and environment deterioration^[@CR1]^. In the 1960s, the wild Chinese alligators become very scarce, and the complicated environments of wetlands and swamps made them even harder to be spotted^[@CR2]^. The population of Chinese alligator didn't stop dropping until 1979 when Changxing Yinjiabian Chinese Alligator Nature Reserve (CYCANR) and Anhui Research Center for Chinese Alligator Reproduction (ARCCAR) were founded to conserve this species, starting with 7 and 212 founder alligators respectively. As of today, there are more than 8,000 captive Chinese alligators at both reserve sites, making the ARCCAR and CYCNAR populations not only the last but also the largest captive populations in China. On the other hand, the number of the Chinese alligators in the wild is dropping 20% every year and the current number is predicted to be lower than 120^[@CR2],\ [@CR3]^.

Judging from the population size alone, the Chinese alligators seem to be recovering quite well from the bottleneck event 30 years ago, and previous studies on the major histocompatibility complex (MHC) genes didn't raise alarms on any of the repercussions of the bottleneck effect. The MHC genes are not only the most functionally polymorphic genes in vertebrates, but also mediate immune responses, thus making them good adaptive markers to assess the evolutionary potential of fighting pathogens^[@CR4]^. Previous studies on the Chinese alligator MHC II genes had shown great polymorphism in exon 2, exon 3 and even intron 2 sequences^[@CR5]--[@CR8]^. After the publication of crocodile, alligator and gharial genomes^[@CR9]--[@CR11]^, the structure and polymorphism of the MHC genes in the crocodilians has been massively annotated and studied^[@CR12],\ [@CR13]^. However, all the above-mentioned studies have been using universal primers to investigate the Crocodylia MHC genes. Considering the evidence of MHC gene duplication^[@CR12],\ [@CR14]^, it's possible that the cross-locus amplification would cause sequences from different loci being regarded as alleles from a single locus, thus distorting the true intra-locus alignment and elevating the rate of non-synonymous (*d* ~N~) over synonymous (*d* ~S~) substitutions. Hence, the true level of genetic diversity at single MHC genes of the Crocodylia is still pending for exploration.

In our previous studies on the Chinese alligator, He *et al*. constructed a bacterial artificial chromosome (BAC) library that contains 6.8-fold genome equivalents^[@CR15]^. Ye *et al*. constructed several contigs containing MHC genes based on the BAC library^[@CR16]^. Wan *et al*. published the Chinese alligator genome as 2.3 Gb in size, and annotated over 22,200 genes^[@CR16]^. These works enable us to design locus-specific primers to characterize MHC genes, hoping to investigate whether bottleneck event continues to impact the Chinese alligator population at the single gene level even after an exponential population growth, to unveil the actual polymorphism of a single MHC gene, and to elaborate the evolutionary forces that influence MHC genes.

Results {#Sec2}
=======

BLAST results and Genotyping of MHC genes {#Sec3}
-----------------------------------------

By BLASTing the BAC-end and known MHC sequences to the *A. sinensis* genome and predicting potential MHC genes, we found 3 MHC loci with intact coding sequences and gene structures. The 3 MHC loci can be pinpointed to different *A. sinensis* genome scaffolds: scaffold1303_1 (NCBI accession number: NW_005842753), scaffold364_1 (NCBI accession number: NW_005842546) and scaffold184_1 (NCBI accession number: NW_005842983). The 3 MHC loci were also verified to be from 3 BAC clones that had no overlap sequences, suggesting that they were independent loci. The BAC clones 1327C2 and 20A2 contained 2 MHC I loci and the BAC clone 1085A9 contained 1 MHC II locus; these MHC loci were accordingly named I1327, I20 and Beta1085, respectively.

The single strand conformation polymorphism (SSCP) results of these loci were all di-morphic: 2 alleles at each locus, which were confirmed by the sequencing results. Nonetheless, the gene I1321 is actually mono-morphic at the functional level due to the same amino acid sequence of the 2 alleles (Fig. [1](#Fig1){ref-type="fig"}). No more than two sequences were present in each animal, demonstrating no cross-locus amplification in our study, and none of the sequences showed deletions, insertions, or stop codons, showing the accurate genotyping at single functional MHC genes. In total, we obtained 6 nucleotide acid alleles but 5 amino acid alleles from the 3 MHC loci using the locus-specific genotyping primers (Fig. [1](#Fig1){ref-type="fig"}), suggesting an extremely low level of adaptive genetic variation in the endangered Chinese alligator.Figure 1Amino acid sequence alignment of 3 MHC loci. Dots represent identical amino acids to the top variant and crosses under the alignment depict putative antigen binding sites.

Allelic distribution, heterozygosity and Hardy-Weinberg Equilibrium test {#Sec4}
------------------------------------------------------------------------

The sequence analysis exhibits a quite low level of nucleotide diversity (0.0025---0.0128) in all the MHC loci (Table [1](#Tab1){ref-type="table"}). These results were unlike previous MHC studies in *A. sinensis* ^[@CR5],\ [@CR6]^, likely because of their use of degenerate primers that are capable of amplifying more than one MHC locus. However, our results coincided with results from Wan *et al*., which revealed low SNP heterozygosity throughout the *A. sinensis* genome^[@CR11]^.Table 1Genetic diversity at the MHC loci.Locus*N* ~A~Nucleotide length (bp)*Hdπ*Exon 2Intron 2Exon 3ExonIntronExonIntronI13272261151352760.444 ± 0.0640.467 ± 0.0170.25%0.20%I202267914; 9252760.457 ± 0.0550.356 ± 0.0250.34%0.20%Beta10852282631; 632---1.000 ± 0.0520.182 ± 0.0210.94%0.06%Note: *N* ~A~, the number of alleles; *Hd*, haplotype diversity; *π*, nucleotide diversity.

Although varying with loci, the allelic distribution patterns of the same MHC gene are similar among populations (Table [2](#Tab2){ref-type="table"}). The Hardy-Weinberg Equilibrium (*HWE*) test results show a severe heterozygote deficiency at the Beta1085 locus in ZJ and AH populations. Differently, the Tajima's D test results are strange (Table [2](#Tab2){ref-type="table"}): below 0 value at MHC class II locus (significant in ZJ population), and above 0 value at both MHC class I loci (significant in ZJ and AH population). The different Tajima's D values between the class I and II MHC loci might indicate different stages of bottleneck or different selection pressures.Table 2Genetic diversity of MHC in 3 Chinese alligator populations.PopulationLocusAllele 01Allele 02*H* ~O~*H* ~E~*P-value of HWE*Tajima's D testZJI1327**0.73**0.270.4690.3960.3961.711AHI1327**0.69**0.310.3750.4370.4412.083\*USAI1327**0.70**0.300.2000.4420.1321.518ZJI200.36**0.64**0.5310.4680.6992.611\*AHI200.22**0.78**0.3750.3471.0001.383USAI200.10**0.90**0.2000.1891.0000.949ZJBeta1085**0.98**0.020.0310.0920.015\*−1.953\*AHBeta1085**0.91**0.090.1250.2310.008\*−1.524USABeta1085**1.00**0.00------------Note: the predominant alleles are in bold; \* indicates *P* \< 0.05; *H* ~O~ and *H* ~E~ are observed and expected heterozygosities, respectively.

The repercussions of bottleneck effects can also be observed at the microsatellite loci (Table [3](#Tab3){ref-type="table"}), where allelic diversity is very low (2.5 alleles per locus) in all 3 populations. Despite the higher *H* ~O~ than *H* ~E~ at several loci, this phenomenon should be attributed to the selection of microsatellite markers with high-resolution power in the study of Ma^[@CR17]^.Table 3Genetic diversity of microsatellite markers in 3 Chinese alligator populations.**LocusZJ (n = 32)AH (n = 32)USA (n = 10)N** ~**A**~**H** ~**O**~**H** ~**E**~**N** ~**A**~**H** ~**O**~**H** ~**E**~**N** ~**A**~**H** ~**O**~**H** ~**E**~CXA-14220.844\*0.50320.5360.39920.8000.505CXA-4120.688\*0.47620.3930.43120.1110.111CXA-620.5000.41120.6150.50720.7000.479CXA-930.250\*0.40530.6540.57821.000\*0.526CXA-3430.500\*0.60830.786\*0.59030.900\*0.679CXA-4330.531\*0.58740.5360.54320.7000.521Mean2.50.5520.4982.6670.5870.5082.1670.7020.470SD0.5480.2000.0860.8160.1320.0780.4080.3120.190Note: *N* ~A~, number of alleles; *H* ~O~, observed heterozygosity; *H* ~E~, expected heterozytosity; \* shows significant deviation from HWE (*P* \< 0.05).

**Calculation of*d***~**N**~**and*d***~**S**~**substitutions** {#Sec5}
--------------------------------------------------------------

The *d* ~N~/*d* ~S~ ratio can provide useful information on the degree of selective pressure acting on a protein-coding gene^[@CR18]^. The *d* ~N~/*d* ~S~ \> 1 implies positive selection; *d* ~N~/*d* ~S~ \< 1 implies purifying (negative) selection; and *d* ~N~/*d* ~S~ = 1 indicates neutral (i.e. no) selection. As we can see in Table [4](#Tab4){ref-type="table"}, the Chinese alligator MHC exons exhibit extremely low *d* ~N~ and *d* ~S~ values due to the low allelic and nucleotide diversity, and the Z-tests show no significant P value at any locus. The following trends, however, can be observed: at the MHC class I loci, the *d* ~N~/*d* ~S~ ratios all exceed 1 at antigen binding sites (ABS), showing the sign of positive selection. Nonetheless, at the MHC class II locus, there is no nucleotide substitution at ABS positions (Fig. [1](#Fig1){ref-type="fig"}), thus indicating purifying selection. When the *d* ~N~/*d* ~S~ values were compared between the ABS and non-ABS sites, the MHC class I loci produced a larger ratio at ABSs than at non-ABSs while the MHC class II locus showed *d* ~S~ at non-ABSs in spite of zero *d* ~N~ changes (Table [4](#Tab4){ref-type="table"}), indicating balancing selection at play in the Chinese alligator.Table 4The *d* ~N~/*d* ~S~ ratio of 3 MHC loci in 3 Chinese alligator populations.LocusPosition*N*ZJAHUSA*d* ~N~*d* ~S~*d* ~N~/*d* ~S~*d* ~N~*d* ~S~*d* ~N~/*d* ~S~*d* ~N~*d* ~S~*d* ~N~/*d* ~S~I1327ABS560.008 ± 0.0050.005 ± 0.0061.600.008 ± 0.0060.006 ± 0.0061.330.009 ± 0.0060.006 ± 0.0061.50non-ABS125000000000All1810.002 ± 0.0020.002 ± 0.0021.000.003 ± 0.0020.002 ± 0.0021.500.003 ± 0.0020.002 ± 0.0021.50ABS560.007 ± 0.0050∞0.005 ± 0.0040∞0.003 ± 0.0020∞I20non-ABS1250.002 ± 0.0020.005 ± 0.0060.400.001 ± 0.0010.004 ± 0.0040.250.001 ± 0.0010.002 ± 0.0020.50All1810.003 ± 0.0020.004 ± 0.0040.750.003 ± 0.0010.003 ± 0.0031.000.001 ± 0.0010.002 ± 0.0010.50ABS49000000---------Beta1085non-ABS13900.001 ± 0.0010.000.002 ± 0.0020.005 ± 0.0040.40---------All18800.001 ± 0.0010.000.001 ± 0.0010.004 ± 0.0030.25---------Note: the *d* ~N~/*d* ~S~ ratios are calculated separately for antigen-binding site (ABS), non-ABS and all sites. Codon numbers (*N*) of the sites are shown. In the USA population, Beta1085 locus has only one allele, thus no *d* ~N~ and *d* ~S~ values are calculated.

Comparison of nucleotide substitution between exon and intron {#Sec6}
-------------------------------------------------------------

All 3 Chinese alligator populations shared the same allele sequences in all 3 MHC loci, and the only differences were the allelic frequencies at each locus. Thus, we computed nucleotide substitution on the basis of these unique allele sequences and without consideration of allelic frequency and population identity. We calculated the mean *d* ~S~ in exons and the mean number of nucleotide substitutions per site (*d*) in introns of the MHC class I and II genes, and plotted them against the nucleotide position (Fig. [2](#Fig2){ref-type="fig"}). The results show a much higher substitution rate in the exon 2 and 3 region than in the intron 2 at all MHC loci, which means the introns are younger than exons, suggesting the exons being constantly maintained by balancing selection for ages.Figure 2A plot of *d* ~S~ in exon 2 (and exon 3 at class I genes) versus *d* in intron2.

Comparison of differentiation degree between the MHC and microsatellite loci {#Sec7}
----------------------------------------------------------------------------

Since the MHC genes are constantly maintained by balancing selection, even in different populations the MHC allelic frequencies will be pretty similar due to the influence of balancing selection. The microsatellite loci, however, will take more diverted evolutionary paths under the pressure of genetic drift. We tested pairwise *F* ~ST~ of the 3 Chinese alligator populations at both MHC and microsatellite loci, and found significant population differentiation among 3 populations at the microsatellite loci (*P* \< 0.05) with no differentiation at the MHC loci (Table [5](#Tab5){ref-type="table"}). This contrast indicates that even though severely damaged by the bottleneck effect, the Chinese alligator MHC loci are still influenced by balancing selection, i.e. MHC genes maintain more similar alleles than neutral markers do.Table 5Pairwise *F* ~ST~ test among 3 Chinese alligator populations.ZJAHUSAZJ---0.014490.04720AH0.08424\*---0.00670USA0.08990\*0.06212\*---Note: the numbers on the left side below the diagonal and on the right side above the diagonal are the *F* ~st~ values of microsatellite and MHC loci, respectively. \* indicates *P* \< 0.05.

Discussion {#Sec8}
==========

Normally, when a population encounters bottleneck events, rare alleles are more likely to be lost than the common alleles, and positive Tajima's D values are expected; when population expands, the segregating sites accumulate at the rare frequencies, thus leading to negative Tajima's D values^[@CR19]^. In the Chinese alligator's case, the bottleneck effect 30 years ago was too severe for the population to attain ample rare alleles with low frequencies, as evidently shown by the di-morphism MHC loci in this study. Therefore, the positive Tajima's D value at the MHC class I loci should be a normal case for the Chinese alligator population. The negative Tajima's D value at the MHC class II locus, however, could indicate the ongoing population expansion or purifying selection, whose trace can also be found in the sequence alignments (Fig. [1](#Fig1){ref-type="fig"}) and *d* ~N~, *d* ~S~ results (Table [4](#Tab4){ref-type="table"}); the MHC class I loci have non-synonymous nucleotide substitutions at and next to the ABS sites whereas the MHC class II locus has no nucleotide substitution at all at any ABS sites.

As our abovementioned results have shown, even after 30 years, the bottleneck effect still influences the MHC genes of the Chinese alligator. That's why conventional methods -- such as *d* ~N~/*d* ~S~ values -- function poorly in detecting balancing selection. The significant *d* ~N~/*d* ~S~ values in previous studies^[@CR5],\ [@CR7]^ are more likely caused by their usage of universal primers, as Miller and Lambert pointed out, the usage of universal primers could produce a lower than expected *d* ~S~ value, leading an elevated *d* ~N~/*d* ~S~ ratio and thus giving "anomalous" results^[@CR20]^.

Therefore, we need to compare the theoretically balancing-selection-influenced sites to those non-balancing-selection-influenced sites to reveal the presence of balancing selection. The comparison between exons and introns would be an effective tool. Hughes and Yeager^[@CR21],\ [@CR22]^ discovered that in most genes, the mean number of nucleotide substitutions per site (*d*) in introns and the mean *d* ~S~ in exons are usually about equal -- for their mutations are both selectively neutral to a gene -- except when it comes to the MHC genes, the mean *d* ~S~ in exons is always much higher than the mean *d* in introns. They explained that exons in the MHC genes are very ancient because they have been maintained by balancing selection for a very long time, while the mutations in introns are controlled by recombination and genetic drift. When an ancient intron polymorphism is lost due to recombination or drift, it will be selectively neutral to the MHC gene. As time goes by, introns in the MHC genes will become evolutionarily younger on average than are the exons constantly maintained by balancing selection. The much older exons will consequently possess a much higher *d* ~S~ than that of the younger introns, hence mean *d* ~S~ \> mean *d* in MHC genes^[@CR22]^. Our findings in this study support Hughes and Yeager's theory.

The environmental factors could be another cause of the similar MHC genes among populations. Both AH and ZJ populations are intensively managed as captive populations, and their habitats are modified to natural wetlands surrounded by farming countryside. The isolated habitats may provide an environment with low level exposure to pathogens, and cause the MHC genes lack motivations to change; similar situations happen to many other species^[@CR20],\ [@CR23]--[@CR25]^. Even if the Chinese alligator populations are sustainable at present, they would still be susceptible to new pathogens when reintroduced to the wild in the future. While the balancing selection is working on the MHC genes, proper conservation strategy should be devised to protect this endangered species.

Materials and Methods {#Sec9}
=====================

Animal experiment ethics statement {#Sec10}
----------------------------------

All experiments were carried out in accordance with the guidelines issued by the Ethical Committee of Laboratory Animal of Zhejiang University, and all experimental protocols were approved by the Ethical Committee of Laboratory Animal of Zhejiang University.

Sampling {#Sec11}
--------

The *A. sinensis* samples (see Supplementary Table [S1](#MOESM1){ref-type="media"}) came from 3 captive populations, kindly provided by CYCANR (the Zhejiang population, acronym ZJ), ARCCAR (the Anhui population, acronym AH) and the Rockefeller Wildlife Refuge (the USA population, acronym USA). Although there are more than 8,000 Chinese alligators in the ZJ and AH populations, most of them are still juvenile while several hundreds of Chinese alligators are breeding adults. In this study, all sample donors were adult male/female Chinese alligators that were randomly captured, and we believe they can represent the species' natural range. The blood samples were taken during routine medical examinations. DNA was extracted using a traditional phenol-chloroform method^[@CR26]^.

Genome BLAST and primer design {#Sec12}
------------------------------

We screened the BAC clones using the universal MHC primers of Ye *et al*.^[@CR16]^ and pinpointed each MHC gene to their corresponding genome scaffolds (requiring \>95% in both gene coverage and sequence similarity) by BLASTing the end sequences of the target BACs and the known MHC sequence against the *A. sinensis* genome^[@CR11]^ using default parameters of a local BLAST tool downloaded from NCBI (ftp://ftp.ncbi.nlm.nih.gov/blast/executables/blast+/LATEST/). Then, we used Augustus (<http://augustus.gobics.de/>) to predict exon positions for each MHC sequence, and manually translated each exon to amino acid sequences in order to search for MHC loci with intact coding sequences and gene structures. All the MHC genes from the BAC and the genome were manually checked for gene structural integrity to rule out pseudogenes. The MHC genes with a complete gene structure were selected for subsequent studies. We gave prefix '*Alsi*' to all the MHC genes.

We then designed locus-specific primers to amplify the antigen binding domains of the MHC loci (exons 2 and 3 for the MHC class I, and exon 2 for the MHC class II). Intron 2 sequences were also amplified and sequenced to ensure the correctness of genotyping as well as to be a useful comparison material (see Supplementary Table [S2](#MOESM1){ref-type="media"}). Primers were designed using Primer Premier 5 (<http://www.premierbiosoft.com/>), and the primers' binding positions are illustrated in Supplementary Fig. [S1](#MOESM1){ref-type="media"}. We performed pre-experiments to investigate the polymorphism of each locus using 10 randomly selected blood samples in each population. We chose the polymorphic loci to genotype all the samples in 3 populations. Six microsatellite markers from our former works were also used as a neutral contrast^[@CR17]^.

PCR amplification and genotyping {#Sec13}
--------------------------------

The PCR amplifications were performed in a 10 μl reaction system containing 1 L of template DNA (30--50 ng/µL), 0.5 U rTaq DNA polymerase (TaKaRa), 1 µL of 10 × ExTaq buffer (TaKaRa), 1 µL of 2.5 mmol/L dNTPs, 0.2 µL of each primer (diluted to 10 mol/µL), and 6.5 µL of double-distilled water. The PCR was programmed as follows: initial denaturation at 95 °C for 5 min, followed by 35 cycles of 95 °C for 30 s, 63 °C for 30 s, and extension at 72 °C for 30 s, with a final extension of 72 °C for 5 min.

We adopted Single-Strand Conformation Polymorphism - heteroduplex (SSCP-HD) techniques^[@CR27]^ to genotype each individual and screen 3 Chinese alligator populations in order to examine the level of adaptive genetic variation in this endangered reptile. In the SSCP-HD analysis, we added 5 μl of 2X loading buffer (95% formamide, 10 mm NaOH, 0.25% bromophenol blue, 0.25% xylene cyanol) for each 10 μl of PCR product. And the mixture was denatured at 95 °C for 5 min, and swiftly transferred onto ice for cooling down. Then, the PCR products were run in an acrylamide gel consisting of 12% 37.5:1 acrylamide/bisacrylamide with 2.5% crosslinking, and sequences were separated in 0.5 × TBE running buffer at 16 °C by 150 V for 6.5 h on the Decode System (Bio-Rad). Finally, the SSCP gel was fixed in 10% acetic acid for 30 min, washed with dH2O, and achieved silver staining pictures. We repeated the PCR-SSCP process for 3 times at each locus to make sure the banding pattern was stable and consistent, and obtained the allele sequences by sequencing the homozygous individuals with unique banding patterns in BGI, Shanghai. Microsatellite loci were amplified and genotyped as described by Ma^[@CR17]^.

Data analyses {#Sec14}
-------------

We used Lasergene 7 (DNASTAR) for nucleotide sequence editing and Mega 5^[@CR28]^ for sequence alignment. We used DnaSP^[@CR29]^ to calculate the haplotype diversity (*Hd*), the nucleotide diversity (*π*), and Tajima's D values. The non-synonymous (*d* ~N~) and synonymous (*d* ~S~) substitution rates of the MHC exons were calculated using Mega 5, and we used K-estimator^[@CR30]^ with Kimura-2p method to compute nucleotide substitution rate of intron 2 (*d*) as well as the *d* ~S~ of exon 2 and exon 3, which were then plotted against nucleotide position using a sliding window size of 15 base pairs and steps of 3 base pairs. We used Kaufman's study^[@CR31]^ as a reference to annotate antigen-binding sites.

Allelic frequency, observed heterozygosity (*H* ~O~), expected heterozygosity (*H* ~E~) and *HWE*, as well as pairwise *F* ~ST~ among 3 populations were calculated in Arlequin 3.5^[@CR32]^.

Data availability statement {#Sec15}
---------------------------

The genomic sequences of MHC genes are collected from the scaffolds 1303_1 (GenBank accession number: NW_005842753), 364_1 (NW_005842546) and 184_1 (NW_005842983) of the Chinese alligator genome (GCA_000455745.1). The new nucleotide acid sequences obtained in this study are available in the supplementary Data [S1](#MOESM1){ref-type="media"}.

Electronic supplementary material
=================================

 {#Sec16}

Supplementary Information

Teng Zhai and Hai-Qiong Yang contributed equally to this work.

**Electronic supplementary material**

**Supplementary information** accompanies this paper at doi:10.1038/s41598-017-05640-2

**Publisher\'s note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Special thanks go to the CYCANR and the ARCCAR for providing samples. This work was supported by a grant from the National Natural Science Foundation of China (No. 31470439), a special grant from the State Forestry Administration, and the Fundamental Research Funds for the Central Universities of the P. R. China.

T.Z. and H.Q.Y. conducted the experiments, interpreted data and wrote the manuscript. R.C.Z., L.M.F. and G.H.Z. collected samples and participated in experiments. S.G.F. conceived the project, designed the research, interpreted data and wrote the manuscript. All authors reviewed and approved the manuscript.

Competing Interests {#FPar1}
===================

The authors declare that they have no competing interests.
