Next Article in Journal
Effect of MSCs and MSC-Derived Extracellular Vesicles on Human Blood Coagulation
Next Article in Special Issue
Genetic Diversity and Differentiation at Structurally Varying MHC Haplotypes and Microsatellites in Bottlenecked Populations of Endangered Crested Ibis
Previous Article in Journal
3,5 Diiodo-l-Thyronine (T2) Promotes the Browning of White Adipose Tissue in High-Fat Diet-Induced Overweight Male Rats Housed at Thermoneutrality
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Cells
Volume 8
Issue 3
10.3390/cells8030257
cells-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Reproductive Strategy Inferred from Major Histocompatibility Complex-Based Inter-Individual, Sperm-Egg, and Mother-Fetus Recognitions in Giant Pandas (Ailuropoda melanoleuca)

by
Ying Zhu
1,†,
Qiu-Hong Wan
1,†,
He-Min Zhang
2 and
Sheng-Guo Fang
1,*
1
MOE Key Laboratory of Biosystems Homeostasis & Protection, State Conservation Centre for Gene Resources of Endangered Wildlife, College of Life Sciences, Zhejiang University, Hangzhou 310058, China
2
China Conservation and Research Center for the Giant Panda, No. 98 Tongjiang Road, Dujiangyan 611800, China
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Cells 2019, 8(3), 257; https://doi.org/10.3390/cells8030257
Submission received: 13 December 2018 / Revised: 18 February 2019 / Accepted: 13 March 2019 / Published: 19 March 2019
(This article belongs to the Special Issue Major Histocompatibility Complex (MHC) in Health and Disease)
Browse Figures
Versions Notes

Abstract

:
Few major histocompatibility complex (MHC)-based mate choice studies include all MHC genes at the inter-individual, sperm-egg, and mother-fetus recognition levels. We tested three hypotheses of female mate choice in a 17-year study of the giant panda (Ailuropoda melanoleuca) while using ten functional MHC loci (four MHC class I loci: Aime-C, Aime-F, Aime-I, and Aime-L; six MHC class II loci: Aime-DRA, Aime-DRB3, Aime-DQA1, Aime-DQA2, Aime-DQB1, and Aime-DQB2); five super haplotypes (SuHa, SuHaI, SuHaII, DQ, and DR); and, seven microsatellites. We found female choice for heterozygosity at Aime-C, Aime-I, and DQ and for disassortative mate choice at Aime-C, DQ, and DR at the inter-individual recognition level. High mating success occurred in MHC-dissimilar mating pairs. No significant results were found based on any microsatellite parameters, suggesting that MHCs were the mate choice target and there were no signs of inbreeding avoidance. Our results indicate Aime-DQA1- and Aime-DQA2-associated disassortative selection at the sperm-egg recognition level and a possible Aime-C- and Aime-I-associated assortative maternal immune tolerance mechanism. The MHC genes were of differential importance at the different recognition levels, so all of the functional MHC genes should be included when studying MHC-dependent reproductive mechanisms.

Graphical Abstract

1. Introduction

Understanding the genetic basis and the driving forces of mate choice in animals has always been a major goal of evolutionary ecologists [1,2]. It has been proposed that females prefer males who can maximize their reproductive success and increase offspring quality/fitness [3,4]. Immunocompetence is undoubtedly an essential index of an individual’s fitness [5,6] and major histocompatibility complex (MHC) genes are suitable candidates in investigating the genetic basis underlying mate choice decisions, as MHC molecules can recognize and present antigens to T cells and trigger immune reactions [7,8,9,10]. A growing body of evidence shows the association between pathogen resistance and MHC haplotypes or alleles [11,12,13,14].
In recent years, an increasing number of studies have focused on MHC-associated mate choice, including three non-exclusive hypotheses that could explain the MHC-based mate choice. Firstly, heterozygous advantage: according to this hypothesis, the choosy sex could obtain additive benefits from mating with MHC-heterozygous mates whose disease resistance might be inherited by their offspring [15,16]. For example, in tuco-tucos (Ctenomys spp.), females prefer MHC-heterozygous males [17]. In the scarlet rosefinch (Carpodacus erythrinus), social males with low MHC heterozygosity are cheated on by their females more frequently than highly MHC-heterozygous males [18]. Secondly, the genetic compatibility hypothesis: the choosy sex is assumed to select MHC-dissimilar partners, resulting in the production of offspring with diverse genotypes that can recognize a broad array of pathogens and hence increase their fitness [15,16,19]. In the grey mouse lemur (Microcebus murinus), the fathers have lower allele sharing and a greater amino acid distance to the mother than the randomly assigned males [20]. A study on blue petrels (Halobaena caerulea) revealed that females mated more with functionally (not evolutionary) MHC-dissimilar males than with random males [21]. Similar patterns have been observed in great frigatebirds (Fregata minor) [22], the Chinese rose bitterling (Rhodeus ocellatus) [23], the pot-bellied seahorse (Hippocampus abdominalis) [24], the fat-tailed dwarf lemur (Cheirogaleus medius) [25], and mandrills (Mandrillus spp.) [26]. Furthermore, to acquire the best immunogenetic composition for their offspring, females would choose mates with the most appropriate MHC diversity (maximum or intermediate), which is also referred to as the optimal hypothesis, and it is an extension of the genetic compatibility hypothesis [10]. Studies in sticklebacks (family Gasterosteidae) have found that females with many alleles prefer males with few alleles, and vice versa, in order to obtain an optimal level of MHC diversity in their offspring for resistance against parasites and pathogens [27,28]. Thirdly, the inbreeding avoidance hypothesis: the choosy sex is expected to seek dissimilar partners, not only with regard to MHC, but also genome-wide, in order to gain fitness benefits [29,30].
Mate choice results from non-random reproductive investment, which could happen at the precopulatory stage or/and at the postcopulatory stage. In the precopulatory stage (or individual recognition level), the choosy sex uses visual, acoustic, or odor cues to choose mates [15]. In the postcopulatory stage, the females increase offspring quality/fitness by their eggs differentiating between sperm during fertilization (sperm-egg recognition), and by following a differential allocation strategy during embryo implantation (mother-fetus recognition) [16]. Human leukocyte antigen (HLA) class I and class II molecules are expressed on sperm cell surfaces [31,32,33], making the complementary cryptic female choice for MHC genotypes possible at the gamete level. Sperm selection that targets different levels of MHC diversity has been reported in several species. In some rodent species, females refuse to accept MHC-similar sperm [34], and in humans, females that mate with a male who has the same HLA haplotype tend to have a greater chance of spontaneous abortion [35,36,37]. In red junglefowl (Gallus gallus), the eggs favor sperm that is from MHC-dissimilar males [38]. In contrast, a fertilization advantage for MHC-similar mates has been observed in Atlantic salmon (Salmo salar) [39], Chinook salmon (Oncorhynchus tshawytscha) [40], and guppies (Poecilia reticulata) [41]. A recent study of the three-spined stickleback (Gasterosteus aculeatus) revealed that sperm selection resulted in offspring with an intermediate level of MHC diversity [42].
The giant panda (Ailuropoda melanoleuca) is an endangered species and its captive breeding has always been focused in China. Giant pandas usually have polyandrous/polygynous multiple mating systems [43]. Females are very choosy, and they have significantly higher copulation and birth rates when paired with preferred males [44]. In breeding programs, natural mating as well as artificial insemination are adopted to increase reproductive success [45]. Although the number of female giant pandas that are available to breed has increased, the fertilization rate is still low (~50% for nearly 20 years). Furthermore, female giant pandas have spontaneous abortions when they are inseminated by sperm from a male that they do not like [46]. Therefore, there may be an MHC-based mechanism that determines mate choice and fertilization in giant pandas.
It is best to investigate all MHC genes concerning mate choice, as MHC class I and class II molecules mainly present intracellular and extracellular pathogen-derived antigens, respectively [9], and the selection of MHC genes might influence mate choice results [47]. Few studies have focused on a large region of the MHC due to a lack of structural knowledge and an effective genotyping method [48], including a giant panda mate choice study that only used three MHC class II genes [49]. In previous studies, we characterized six functional MHC class II genes (Aime-DRA, Aime-DRB3, Aime-DQA1, Aime-DQA2, Aime-DQB1, and Aime-DQB2) [50,51] and four classical MHC class I genes (Aime-C, Aime-F, Aime-I, and Aime-L) [52], and we developed their genotyping protocols in the giant panda. Therefore, our previous studies provide a good foundation to investigate the relationship between a large number of MHC genes and female choice at the inter-individual, sperm-egg, and mother-fetus recognition levels in giant pandas.
In the present study, we took advantage of multiple years of observations and the possession of genetic data of a captive population in the Wolong Chinese Research and Conservation Center for the Giant Panda. We aimed to: (1) test three mate choice hypotheses at the individual recognition level (MHC-heterozygous choice, MHC-compatibility choice, and inbreeding avoidance) while using 10 MHC genes and seven microsatellites. If the heterozygote advantage is the main driving force of female mate choice, more diverse males should be preferred than less diverse ones, regardless of the females’ MHC genotypes. If female choice favors the production of offspring with high fitness that is based on MHC compatibility, we expect that females would choose MHC-dissimilar partners that provide the most appropriate level of MHC dissimilarity (maximum or intermediate) in the offspring. If mate choice aims to avoid inbreeding, then we expect genome-wide dissimilarity between partners (MHC and microsatellites); (2) test whether cryptic female choice for MHC compatibility occurs at the gamete level; and, (3) characterize successful embryo implantations by comparing the MHC genotypes of mothers and their offspring.

2. Materials and Methods

2.1. Study Species and Behavioral Observations of Inter-Individual Recognition

Giant pandas were housed in the Wolong Chinese Research and Conservation Center for the Giant Panda. Female giant pandas come into estrus from February to June [43], when their appetite decreases, males’ urination frequency increases, and they rub their genitalia against the ground or walls [43,46]. Two rutting females and males are usually put into two adjacent cages. Neighboring pandas have access to each other through a cage fence with full sight, hearing, and smell, but limited touch. Female pandas face multiple males that are consecutively presented. Veterinarians use pre-mating behaviors to determine when the males are sent to females for mating. If a female is interested in a male, then they usually respond by sniffing and pushing the fence, before the male is sent to the female’s cage to copulate. If a female is uninterested, then she does not respond or behaves aggressively. Female mate choice at the inter-individual recognition level was determined by naturally analyzing mating pairs that were recorded in a studbook [53]. There were 33 females and 21 males in the breeding program from 1991 and 2008, except for 1994, but we failed to obtain samples from five females and four males. Therefore, for the inter-individual recognition study, we included 182 natural mating events that involved 28 females and 17 males.

2.2. Sperm-Egg and Mother-Fetus Recognition

Female giant pandas have an annual estrus cycle with spontaneous ovulation [54,55]. They naturally mate with multiple males and require artificial insemination with males that they mate with and with males that they do not mate with but have good-quality sperm. Sperm from all males have a chance to access the egg, but only one of them is successful through a mechanism at the gamete stage, and sperm-egg recognition may occur at the gamete stage. Over 17 years, a total of 80 offspring (zygotes observed) were produced. Other zygotes were combinations of egg-to-sperm haplotypes, except for offspring haplotypes. The information that is required for natural mating and artificial insemination was acquired from the studbook and SPARKS 1.5 [53].
In addition, we compared the difference between the zygotes observed and other zygotes (combinations of egg-to-sperm haplotypes, except for offspring haplotypes) to mothers in successful embryo implantation events, and called it the “mother-fetus recognition level”.

2.3. DNA Extraction, MHC and Microsatellite Genotyping, and the Definition of a Super Haplotype

We collected 110 blood and the fecal samples. Blood samples were obtained during a routine medical examination and preserved in liquid nitrogen. Fecal samples were less than two days old and stored in 95% ethanol. Genomic DNA extraction from the blood and fecal samples was conducted, as described by Wan et al. [56].
We used seven microsatellite loci (Aim-3, Aim-5, Aim-10, Aim-11, Aim-13, Aim-14, and Aim-16) that performed well in a previous paternity test [57] to conduct paternity analysis for 113 individuals. Polymerase chain reaction (PCR) amplification and genotyping mirrored that described by Zhang et al. [58].
We used 10 functional Aime-MHC loci, including four class I loci (Aime-C, Aime-F, Aime-I, and Aime-L) and six class II loci (Aime-DRA, Aime-DRB3, Aime-DQA1, Aime-DQA2, Aime-DQB1, and Aime-DQB2). We genotyped all individuals at the polymorphic MHC class I exon 2–3 regions and MHC class II exon 2. The primer sets, PCR amplification, and genotyping were as described in two previous studies [52,59].
A physical MHC map revealed that as well as four MHC class I genes, there were also six MHC class II genes linked together [52,60]. In addition to using the above-mentioned independent MHC loci, we analyzed the allele linkage relationships of four MHC class I genes, six MHC class II genes, all MHC genes, genes in the DQ region, and genes in the DR region. For example, we identified homozygotes in four MHC class I genes of the offspring and then inferred the linkage relationship between the mother and father according to Mendel’s law. We named the linkage of four MHC class I genes as SuperHaplotypeI (SuHaI). Using the same procedure, we named the linkage relationships between six MHC class II genes and all MHC genes SuHaII and SuHa, respectively. Genes in the DQ and DR regions were simply called DQ and DR, respectively.

2.4. Paternity Test

We performed a paternity analysis with seven microsatellite loci while using an exclusive method that was based on Mendel’s law, and MHC genotype data confirmed the results. Information on mothers was obtained from the studbook [53].

2.5. Data Analysis

2.5.1. Female Mate Choice at the Inter-Individual Recognition Level

Three mate choice hypotheses, heterozygote advantage, genetic compatibility, and inbreeding avoidance were tested at the inter-individual recognition level.
We tested whether females prefer heterozygous males by utilizing three parameters: (1) The number of heterozygotes (Hobs) at MHC loci in males. We used 0 and 1 to represent homozygote and heterozygote, respectively; (2) Multilocus heterozygosity (MLH) in males, i.e., the proportion of heterozygous microsatellite loci that accounted for the total number of microsatellite loci; (3) d2 value (microsatellites), i.e., the genetic distance between the two alleles. We calculated the d2 for each individual according to the formula d2 = 1/n Σn (ai − aj)2, where ai and aj refer to the repeat length of two individual alleles [61].
We tested whether female giant pandas prefer MHC-dissimilar or -similar males (choice for genetic compatibility) with three parameters: (1) The number of MHC alleles that is shared between females and males as “Nas”, which is twice the number of alleles that is shared by females and males divided by the total number of females and males [62]; (2) The pairwise functional amino acid distance between MHC alleles of females and males, calculated as “Faadis” = Dab + DaB + DAB + Dab, where A, a, B, and b are four alleles in two mates [63]. Each amino acid was characterized by five physicochemical variables: z1 (hydrophobicity), z2 (steric bulk), z3 (polarity), and z4 and z5 (electronic effects) [64]. The functional amino acid distance was the Euclidean metric of the two vectors, which consists of two alleles from a female and male [65]. We not only considered the distance at all sites, but also at the antigen binding sites (ABSs), as it is ABSs that determine pathogen binding and recognition and they are considered a functional region [14,66,67]. We defined ABSs according to human sequences [68].
We tested the inbreeding avoidance hypothesis using Queller and Goodnight’s relatedness (microsatellites) [69], as well as Nas and Faadis. Queller and Goodnight’s relatedness was calculated in SPAGeDi v 1.4 [70] to estimate the genetic similarity between female and male giant pandas. If mate choice aims to avoid inbreeding, we would expect to see choice that is based on compatibility at the genome level (MHC and microsatellites), or female mate choice was targeted to MHCs.
We adopted three approaches to test each hypothesis. Firstly, we performed a randomization test, which is a nonparametric approach that is based on Monte Carlo sampling, to test whether female giant pandas randomly choose their partners for mating. We simulated the natural mating scenario (same accessible males for each female in the respective year), with the exception that females randomly chose males. We generated a null distribution by allowing each female to randomly choose 10,000 times between all males in the respective year. We then compared the mean of the values that were obtained from 182 naturally mated pairs or males involved in natural mating with a set of randomly matched pairs. We calculated the exact p values as twice the proportion of the simulations, which gave higher values than those observed. If females preferred heterozygous males, we expected that the mean values (Hobs, MLH, and d2) of males that were naturally mated would be significantly higher than those of the randomly assigned males. If females preferred MHC-dissimilar males, we expected that the mean Nas of naturally mated pairs would be significantly lower than that of the randomly assigned pairs, or that the Faadis of naturally mated pairs would be significantly higher than that of the randomly assigned pairs. If females aim to avoid inbreeding, then we expected that the mean values (Faadis and relatedness) would be higher than those of randomly assigned pairs. The randomization tests were conducted in ResamplingStats v 4.0 (ResamplingStats Inc., Arlington, TX, USA).
Secondly, we used a paired Student’s t-test (or Wilcoxon test if the model assumptions were not met) to test whether the natural mating group exhibited higher heterozygosity and/or greater MHC divergence than non-mating group by comparing the mean values of males that are involved in natural mating and those of males that were not involved in mating. Hobs, MLH, and d2 were used to test the heterozygote advantage hypothesis, and Nas, Faadis, and relatedness were used to test the compatibility and inbreeding avoidance hypotheses.
We employed a model-based approach to explore the association between genetic variables and the probability of natural mating success. We performed a generalized linear mixed model (GLMM) with a binomial distribution and logic link function for 956 paring events, and the dependent variable was coded as 1 for a natural mating event and 0 for those that did not. The MHC genetic variables (Hobs, Nas, and Faadis) were included as the fixed effects and female ID (to avoid female pseudoreplication), male ID, and year (to account for differences among breeding years) were included as random effects.

2.5.2. Sperm-Egg Recognition Level

We used a randomization test, a paired Student’s t-test, and a GLMM to test the MHC compatibility hypothesis at the sperm-egg recognition level.
Randomization was based on the data obtained (same males for fertilization or natural mating in the respective year), with the exception that the eggs randomly chose sperm. We generated a null distribution by allowing each egg to randomly choose 10,000 times between all of the sperm of the respective year. Subsequently, we compared the mean MHC divergence from the zygotes that were observed (Faadis = DAB, where A and B are two haplotypes in two gametes, N = 65) with a distribution of egg-to-sperm divergence from randomly assigned zygotes. The exact p values were calculated, as described above.
Paired Student t-tests (or Wilcoxon tests if the model assumptions were not met) were used to test whether the observed combination of eggs and sperm indicated greater MHC divergence by comparing the mean values of the zygotes observed (offspring) with other zygotes.
We used a GLMM with a binomial distribution and logic link to test whether MHC divergence between eggs and sperm influenced breeding success (successful births). The dependent variable was coded as 1 for a breeding event and 0 for a non-breeding event. Egg-to-sperm MHC divergence (Faadis) was included as a fixed effect and female ID, male ID, and year were included as the random effects.
The paired Student t-tests, Wilcoxon tests, and GLMMs were conducted in SPSS v. 20.0 (SPSS Inc., Chicago, IL, USA).

2.5.3. Mother-Fetus Recognition Level

We used paired Student t-tests or Wilcoxon tests, as appropriate, to test whether the zygotes observed (offspring) had lower MHC divergence (Faadis) from the mother than other zygotes (combinations of egg-to-sperm haplotypes, except for offspring haplotypes) by comparing the mean Faadis of zygotes observed to the mother with other zygotes to the mother.
A GLMM with a binomial distribution and logic link was used to explore the association between genetic variables and the probability of breeding success. The dependent variable was coded as 1 for a breeding event and 0 for a non-breeding event. The MHC divergence of zygotes from the mother was included as a fixed effect, and the random effects mirrored those described above.
We adjusted the p values with a false discovery rate of 20%, following the Benjamini–Hochberg procedure for multiple testing [71,72].

2.6. Declarations Ethics Statement

All blood and fecal samples were collected from captive pandas that were housed in the Wolong Chinese Research and Conservation Center for the Giant Panda. Blood samples were collected during routine examinations with permission from the China Giant Panda Protection and Management Office. We obtained specific permission from the China Research and Conservation Center for the Giant Panda to take fecal samples from captive individuals during the non-breeding season. Permission to use the samples was given by the State Conservation Center for Gene Resources of Endangered Wildlife of China, where they were deposited.

3. Results

3.1. Microsatellite and MHC Diversity

The male MLH ranged from 0.286 to 1 (median = 0.714), d2 ranged from 0.714 to 41.571 (median = 25.857), and Queller and Goodnight’s relatedness ranged from −0.484 to 0.829 (median = −0.070).
A total of 47 MHC sequences were isolated, which are the same as the published sequences, and included 22 MHC class I alleles and 25 MHC class II alleles: seven Aime-C (Aime-C*01–08, JX987000–JX987005, and JX987007), one Aime-F (Aime-F*01 and JX9870008), seven Aime-I (Aime-I*01–07 and JX987009–JX987015), seven Aime-L (Aime-L*01–06 and JX987016–JX987021), seven DQA1 (DQA1*01–07), three DQA2 (DQA2*01–03), 6 DQB1 (DQB1*01–06), seven DRB3 (DRB3*01–05, 07–08), one DRA, and one DQB2. We obtained only one allele at Aime-F, Aime-DRA, and Aime-DQB2, so we excluded these three loci from the analysis.
Furthermore, we obtained 36 SuHa, 20 SuHaI, 18 SuHaII, 18 DQ, and seven DR. The linkage relationships of SuHaII, SuHaI, SuHa, and DQ are presented in Figures S1–S4, respectively. The number of super haplotypes in the DR region was the same as the allele number at DRB3, as there were only two loci in the DR region, with DRA being a homozygote.
The mean number of variable amino acid sites in MHC class II genes was higher than that in MHC class I genes with respect to the whole exon and all of the ABS sites (12.8% vs. 11.5%, respectively, and 41.0% vs. 25.7%, respectively). This was found by comparing the variable amino acid sites in SuHaI and SuHaII (Table S1). However, the mean difference between pairwise alleles in MHC class I genes was greater than that in MHC class II genes (Table S1), as was the case for SuHaI and SuHaII. Moreover, the number of variable amino acid sites and the difference between the pairwise alleles in the DR region were both greater than those in the DQ region.

3.2. Inter-Individual Recognition

The number of males that were accessible to females ranged from three to 11 during the 17-year period, and females naturally mated with 1–4 males.

3.2.1. Heterozygosity Advantage

Males Involved in Natural Mating Versus Randomly Assigned Males

Concerning super haplotypes, males that were involved in natural mating had significantly more heterozygotes at SuHa, SuHaII, and DQ than the randomly assigned males (PSuHa = 0.009, PSuHaII = 0.011, and PDQ = 0.001; Figure 1a). Regarding individual loci, males that are involved in natural mating had more heterozygotes at Aime-C, Aime-I, and Aime-DQB1 than randomly assigned males (PC = 0.020, PI = 0.000, and PDQB1 = 0.000; Figure 1b). However, there were no significant differences in heterozygosity at any other super haplotype or locus between the males that naturally mated and randomly assigned males (Figure 1). Furthermore, we found no significant difference in MLH or d2 between males that naturally mated and those that were randomly assigned (Table 1).

Males Involved in Natural Mating Versus Natural Non-Mating Males

Males that were involved in natural mating had a higher proportion of heterozygotes at SuHa, SuHaI, SuHaII, and SuHaDQ than those natural non-mating males (males not involved in mating, PSuHa = PSuHaI = PSuHaII = PDQ = 0.000; Figure 2a). This pattern was found at all loci, except Aime-L and Aime-DRB (PC = 0.000, PI = 0.000, PDQA1 = 0.001, PDQA2 = 0.007, and PDQB1 = 0.000; Figure 2b). We found no significant differences in MLH or d2 between males that were involved in natural mating and other males (Table 1).

Relationship between Male MHC Heterozygosity and Natural Mating Success

Being encouraged by the results above, we combined seven polymorphic MHC loci (DQA1, DQA2, DQB1, DRB3, Aime-C, Aime-I, and Aime-L) to model the relationship between male MHC heterozygosity and natural mating success, and found a nonsignificant relationship between overall male MHC heterozygosity and natural mating success (F = 2.963, df = 954, P = 0.86). Female ID, male ID, and year did not have any effect on natural mating success.

3.2.2. Genetic Compatibility and Inbreeding Avoidance

Naturally Mated Pairs versus Randomly Assigned Pairs

There was a significant difference in Nas at SuHa, SuHaII, DQ, and DR between naturally mated pairs and randomly assigned pairs (PSuHa = 0.006, PSuHaII = 0.000, PDQ = 0.000, and PDR = 0.000; Figure 3a). Males that were involved in natural mating had significantly lower Nas at Aime-DQA1, DQA2, DQB1, and DRB3 than the randomly assigned pairs (PDQA1 = 0.038, PDQA2 = 0.001, PDQB1 = 0.000, and PDRB3 = 0.000; Figure 3b). However, we did not find any significant differences in Nas at SuHaI or other loci (Figure 3).
Naturally mated pairs had significantly higher Faadis at SuHa, SuHaII, DQ, DR, Aime-C, DQA1, DQA2, DQB1, and DRB3 than randomly assigned pairs with respect to all sites and the ABS sites (Table 2). However, we did not find a significant difference in Faadis at SuHaI, Aime-I, or Aime-L between the pairs with respect to all sites and ABS sites (Table 2). In addition, there was no significant difference between the naturally mated pairs and randomly assigned pairs in genetic relatedness at microsatellites (Table 1).

Naturally Mated Pairs versus Non-Mating Pairs

The comparison between naturally mated pairs and non-mating pairs that are based on allele sharing, functional amino acid distances, and genetic relatedness was the same as between naturally mated pairs and randomly assigned pairs (Figure 4 and Figure 5 and Table 1 and Table 3). We found no evidence to support the inbreeding avoidance hypothesis when comparing the results of the MHC genes (significant) and microsatelites (nonsignificant).

Relationship between MHC Divergence between Females and Males and Natural Mating Success

The model-based analysis confirmed the randimization results. We found a negative relationship between allele sharing females and males and natural mating success at SuHa, SuhaII, DQ, and DR. Separately examining the MHC loci revealed that the overall negative trend might have been caused by lower allele sharing at DQA1, DQA2, DQB1, and DRB3 between the mated pairs (Figure 6). A positive relationship was found between the functional amino acid distance of females to males and natural mating success at SuHa, SuhaII, and DQ, and at three individual loci (DQA1, DQA2, and DQB1; Table 4). The model-based analysis indicated that higher natural mating success resulted from more MHC-dissimilar mating pairs. These results support the MHC genetic compatibility hypothesis.

3.3. Sperm-Egg Recognition

We successfully assigned 16 fathers to 80 offspring, including 30 twins. In 65 cases, adult females accepted sperm from 1–4 adult males, including males that are involved in natural mating and artificial insemination.

3.3.1. Observed Zygotes (Offspring) Versus Randomly Assigned Zygotes

We found no significant difference between the observed zygotes and the randomly assigned zygotes in Faadis (egg-to-sperm MHC divergence) at all types of super haplotype and all MHC loci, except for DQA1 and DQA2 (Table S2). The Faadis of the observed zygotes at DQA1 and DQA2 was outside the 97.5% hypothesis distribution with respect to all sites and the ABS sites (Figure 7, Table S2).

3.3.2. Observed Zygotes (Offspring) Versus Other Zygotes

The results between the observed zygotes and other zygotes (combinations of egg-to-sperm haplotypes, except for offspring haplotypes)) were similar to those between the observed zygotes and randomly assigned zygotes (Table 5, Figure 8). Observed zygotes had higher Faadis values at DQA1 and DQA2 than other zygotes at all sites and ABS sites (PDQA1-ABS = 0.008, PDQA1-ALL = 0.012, PDQA2-ABS = 0.000, and PDQA2-ALL = 0.000; Figure 8).

3.3.3. Relationship between Breeding Success and MHC Divergence of Zygotes

The model-based analysis revealed that a large MHC functional amino acid distance between eggs and sperm at DQ (mainly at DQA1 and DQA2) and Aime-C resulted in higher breeding success, indicating the preference for maximum MHC divergence between eggs and sperm (Table 6).

3.4. Mother-Fetus Recognition

There was no significant difference in Faadis between the observed zygotes and their mothers and other zygotes and their mothers at all super haplotypes and all loci, except for Aime-C and Aime-I (Figure 9, Table S3). The observed zygotes had significantly lower Faadis to mothers at Aime-C and Aime-I than other zygotes (p < 0.05, Figure 9).
We did not find a significant relationship between breeding success and the functional amino acid distance of mothers to a combination of egg and sperm haplotypes at any MHC loci (Table S4).

4. Discussion

4.1. Female Choice at the Inter-Individual Recognition Level

We found that female giant pandas preferred males with MHC-heterozygous and MHC-dissimilar genotypes, which supports the heterozygosity advantage and disassortative choice of the compatibility hypotheses, respectively. It has been proposed that females usually use male ornaments as a cue for choosing MHC-heterozygous males, and olfaction as a cue for choosing MHC-dissimilar males [15]. Body color varies among giant pandas, and it may reflect condition [73,74,75]. Therefore, we propose that female giant pandas assess males’ physical condition in order to choose MHC-heterozygous mates by evaluating the depth of body color. During the mating season, giant pandas secrete musk from their anal glands and in their urine, which is highly odorous [43], thus it is possible that females use odor to choose among male MHC genotypes. Animals can simultaneously use two cues, and the two strategies are not mutually exclusive, particularly in species with both good vision and olfaction [19,76,77]. A recent study on giant pandas reported multimodal signal behavior between females and males, including olfaction, vision, and hearing [55], supporting the use of two strategies. It has been reported that female pandas use odor cues prior to face-to-face meetings, so we propose that females choose MHC-dissimilar mating partners based on odor cues alone, and they use body color depth to visually choose MHC-heterozygous males. A “multiple-cue strategy” was also found in a study of female mice that changed strategy when the males’ urinary scent-marking rate changed [78], and in lizards that use coloration as a mate choice cue at long distances, but use odor at short distances [79]. Using multiple cues to choose a mate may be more common than expected, as females gain more information and reduce their selection costs [80].
The preference for MHC-dissimilar mates might be due to inbreeding avoidance in giant panda populations. In general, if mate choice aims to avoid inbreeding, then we expected to see significantly lower relatedness in the observed pairs than in randomly assigned pairs [20,63,81]. However, our results revealed no significant difference between observed pairs and randomly assigned pairs (Table 1), suggesting that female giant pandas do not attempt to avoid inbreeding. Avoiding inbreeding may not be as important to females as maximizing the number of MHC polymorphisms in their offspring in order to resist pathogens. In addition, we did not find any evidence for the heterozygote advantage hypothesis that is based on the MLH and d2 results, suggesting that overall genetic diversity has no effect on mate choice. These findings indicate that female mate choice targets functional MHC genes rather than other regions, or it is a byproduct of inbreeding avoidance.
Some studies have reported that giant pandas are susceptible to parasites and viruses [82,83,84,85], suggesting that the MHC genes are important in female mate choice at the inter-individual recognition level. Furthermore, our results show that females favored partners that were the most MHC-dissimilar to themselves, resulting in offspring with high heterozygosity. Whether offspring have high immunocompetence should be addressed in future studies regarding the relationship between MHC heterozygosity and immunocompetence in the giant panda.

4.2. Cryptic Female Choice at the Sperm-Egg Recognition Level

Female giant pandas usually mate with multiple males, possibly to increase fertilization success and offspring genetic quality, as has been found in other species [3,86,87]. Our results revealed that sperm and eggs do not randomly combine, and that sperm from zygotes observed was more dissimilar to eggs at DQA1 and DQA2 than sperm from other zygotes (Figure 7 and Figure 8 and Table 5). Two mechanisms may explain this: cryptic female choice or sperm competition [15,16,20], but we could not identify which was the most important. Nevertheless, combined sperm and eggs had the maximum functional amino acid distance, which supports the female MHC-disassortative choice of compatibility hypothesis that is described above.

4.3. Mother-Fetus Recognition Level

The observed zygotes were more similar to mothers at Aime-C and Aime-I than the randomly assigned zygotes and other zygotes (Figure 9). These findings suggest that the observed zygotes had higher compatibility at Aime-C and Aime-I, which could decrease or block graft rejections from the mother’s immune system. MHC class I molecules play an important role in the immune reaction between the mother and fetus, e.g., HLA-G molecules are present at the mother-child interface of trophoblastic cells and protect the fetus from the lytic activity of maternal uterine natural killer cells [88,89]. The identification of an HLA-G ortholog in the giant panda would elucidate the effects of mother-fetus immunity in this species; however, MHC class I genes have similar loci, making HLA-G orthologs difficult to identify in the giant panda.

4.4. Hierarchical and Cooperative Effects

From the super haplotype level to individual loci, female mate choice was hierarchical in nature. For example, we found that SuHa, which represents all of the functional MHC genes in giant pandas, predicted MHC-heterozygous female choice. However, when separately analyzing SuHaI and SuHaII, the effect could only be observed in SuHaII. When we then excluded DR from SuHaII and analyzed DQ alone, the heterozygosity advantage could still be detected. In contrast, the heterozygosity advantage could not be detected when only analyzing DR. Finally, when we focused on the individual locus level, only one out of three DQ genes (DQB1) still predicted female choice. A similar pattern was found for the genetic compatibility hypothesis.
Alternatively, multiple MHC genes that act cooperatively can explain the above example. The effect of DQ seemed to be larger than that of DQB1 alongside DQA1 and DQA2. The integration of the loci or super haplotypes had a greater effect than the sum of the loci or super haplotypes.
The SuHa results were inconsistent with those of SuHaII with respect to inter-individual recognition, suggesting that SuHaII is more important than SuHaI. This may have been caused by greater variation in SuHaII than in SuHaI.

4.5. Relative Importance of MHC-I, MHC-II, DQ, and DR

MHC class I, MHC class II, DQ, and DR were differentially important at the inter-individual, gamete, and mother-fetus recognition levels. More MHC II genes predicted compatible female mate choice than the MHC I genes, e.g., naturally mated males were more dissimilar to their partners at SuHaII but not at SuHaI. Furthermore, all four of the MHC class II genes (DQA1, DQA2, DQB1, and DRB3) predicted compatible female mate choice, while only one out of three MHC class I genes had a significant result. At the gamete level, the sperm were more dissimilar to eggs at DQ than at DR. At the mother-fetus recognition level, the observed zygotes (offspring) were more similar to their mothers at MHC class I genes, but not at MHC class II genes.
Genes that predict female mate choice or gamete selection may be more important than other MHC loci that are involved in pathogen resistance, as revealed by many studies in which distinct molecules that are coded by different alleles recognize specific pathogens, and their ability to resist pathogens differs [11,12,13,14], such as the DQ region. The DQ region in the giant panda differs to that in other mammals, because it contains more genes and alleles than the DR region [59]. DRB3 exhibited the most variation among the seven individual polymorphic loci. MHC polymorphisms may be driven by sexual selection [10], which is in line with our results at DQ and DRB3.
Our findings suggest that MHC loci do not play equal roles in female mate choice at the inter-individual recognition level or at other levels, and that targeted MHC loci may be key for female mate choice at the individual recognition level. Therefore, more MHC loci should be isolated and large MHC regions surveyed.

5. Conclusions

Our results for Aime-C, Aime-I, and DQ support the heterozygosity hypothesis, while the results for Aime-C, DQ, and DR support the genetic compatibility hypothesis. Comparisons of combined sperm and other sperm-to-eggs revealed that sperm competition or Aime-DQA1- and DQA2-associated gamete selection occurred. The comparison of zygotes observed (offspring) and other zygotes revealed the possible Aime-C- and Aime-I-associated maternal immune tolerance mechanisms. We suggest that captive breeding programs should consider the MHC constitution. Our study provides a good foundation for studying the relationships between the MHC constitution, individual fitness (lifetime reproductive success), and mate choice cues in the giant panda.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4409/8/3/257/s1. Figure S1, Allele linkage relationships of six major histocompatibility complex (MHC) class II genes for SuHaII. Asterisk indicates recombinational SuHaII. Figure S2, Allele linkage relationships of four major histocompatibility complex (MHC) class I genes for SuHaI. Figure S3, Allele linkage relationships between SuHaI and SuHaII for SuHa. Figure S4, Allele linkage relationships of DQ genes. Asterisk indicates recombinational DQ super haplotype. Table S1, Genetic variation in major histocompatibility complex (MHC) class I and class II molecules and super haplotypes. Table S2, Functional amino acids in observed zygotes and randomly assigned zygotes. Table S3, Mean functional amino acid distances at five super haplotypes between mothers and observed zygotes and between mothers and other zygotes. Table S4, Association between breeding success and major histocompatibility complex (MHC) functional amino acid distances of the mother and a combination of egg and sperm haplotypes.

Author Contributions

Conceptualization, S.-G.F. and Q.-H.W.; Methodology, Q.-H.W. and Y.Z.; Experiment and Analysis, Y.Z, Q.-H.W. and S.-G.F.; Resources, H.-M.Z.; Writing–Original Draft Preparation, Y.Z, Q.-H.W. and S.-G.F.; Writing–Review & Editing, Y.Z., Q.-H.W. and S.-G.F.; Supervision, S.-G.F. and Q.-H.W.; Project Administration and Funding Acquisition, S.-G.F.

Funding

This work was supported by a National Key Program (2016YFC0503200) from the Ministry of Science and Technology of China, a special grant for the giant panda from the State Forestry Administration of the People’s Republic of China, and the Fundamental Research Funds for the Central Universities of the People’s Republic of China.

Acknowledgments

We thank China Research and Conservation Center for the giant panda (Wolong) for precious samples of giant panda.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Andersson, M.B. Sexual Selection; Princeton University Press: Princeton, NJ, USA, 1994. [Google Scholar]
  2. Andersson, M.; Simmons, L.W. Sexual selection and mate choice. Trends Ecol. Evol. 2006, 21, 296–302. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Jennions, M.D.; Petrie, M. Why do females mate multiply? A review of the genetic benefits. Biol. Rev. 2000, 75, 21–64. [Google Scholar] [CrossRef]
  4. Searcy, W.A. The evolutionary effects of mate selection. Annu. Rev. Ecol. Syst. 1982, 13, 57–85. [Google Scholar] [CrossRef]
  5. Folstad, I.; Karter, A.J. Parasites, bright males, and the immunocompetence handicap. Am. Nat. 1992, 603–622. [Google Scholar] [CrossRef]
  6. Hamilton, W.D.; Zuk, M. Heritable true fitness and bright birds: A role for parasites? Science 1982, 218, 384–387. [Google Scholar] [CrossRef]
  7. Klein, J. The Natural History of the Major Histocompatibility Complex; Wiley & Sons: New York, NY, USA, 1986. [Google Scholar]
  8. Sommer, S. The importance of immune gene variability (MHC) in evolutionary ecology and conservation. Front. Zool. 2005, 2, 16. [Google Scholar] [CrossRef] [PubMed]
  9. Piertney, S.B.; Oliver, M.K. The evolutionary ecology of the major histocompatibility complex. Heredity 2006, 96, 7–21. [Google Scholar] [CrossRef] [PubMed]
  10. Milinski, M. The major histocompatibility complex, sexual selection, and mate choice. Annu. Rev. Ecol. Evol. Syst. 2006, 37, 159–186. [Google Scholar] [CrossRef]
  11. de Campos-Lima, P.O.; Gavioli, R.; Zhang, Q.J.; Wallace, L.E.; Dolcetti, R.; Rowe, M.; Rickinson, A.B.; Masucci, M.G. HLA-A11 epitope loss isolates of Epstein-Barr virus from a highly A11+ population. Science 1993, 260, 98–100. [Google Scholar] [CrossRef] [PubMed]
  12. Grimholt, U.; Larsen, S.; Nordmo, R.; Midtlyng, P.; Kjoeglum, S.; Storset, A.; Saebø, S.; Stet, R.J.M. MHC polymorphism and disease resistance in Atlantic salmon (Salmo salar); facing pathogens with single expressed major histocompatibility class I and class II loci. Immunogenetics 2003, 55, 210–219. [Google Scholar] [CrossRef]
  13. Harf, R.; Sommer, S. Association between major histocompatibility complex class II DRB alleles and parasite load in the hairy-footed gerbil, Gerbillurus paeba, in the southern Kalahari. Mol. Ecol. 2004, 14, 85–91. [Google Scholar] [CrossRef] [PubMed]
  14. Langefors, Å.; Lohm, J.; Grahn, M.; Andersen, Ø.; von Schantz, T. Association between major histocompatibility complex class IIB alleles and resistance to Aeromonas salmonicida in Atlantic salmon. Proc. R. Soc. Lond. B Biol. Sci. 2001, 268, 479–485. [Google Scholar] [CrossRef] [PubMed]
  15. Mays, H.L.; Hill, G.E. Choosing mates: Good genes versus genes that are a good fit. Trends Ecol. Evol. 2004, 19, 554–559. [Google Scholar] [CrossRef] [PubMed]
  16. Neff, B.D.; Pitcher, T.E. Genetic quality and sexual selection: An integrated framework for good genes and compatible genes. Mol. Ecol. 2005, 14, 19–38. [Google Scholar] [CrossRef] [PubMed]
  17. Cutrera, A.P.; Fanjul, M.S.; Zenuto, R.R. Females prefer good genes: MHC-associated mate choice in wild and captive tuco-tucos. Anim. Behav. 2012, 83, 847–856. [Google Scholar] [CrossRef]
  18. Promerová, M.; Vinkler, M.; Bryja, J.; Poláková, R.; Schnitzer, J.; Munclinger, P.; Albrecht, T. Occurrence of extra-pair paternity is connected to social male’s MHC-variability in the scarlet rosefinch Carpodacus erythrinus. J. Avian Biol. 2011, 42, 5–10. [Google Scholar] [CrossRef]
  19. Tregenza, T.; Wedell, N. Genetic compatibility, mate choice and patterns of parentage: Invited review. Mol. Ecol. 2000, 9, 1013–1027. [Google Scholar] [CrossRef] [PubMed]
  20. Schwensow, N.; Eberle, M.; Sommer, S. Compatibility counts: MHC-associated mate choice in a wild promiscuous primate. Proc. R. Soc. B Biol. Sci. 2008, 275, 555–564. [Google Scholar] [CrossRef]
  21. Strandh, M.; Westerdahl, H.; Pontarp, M.; Canback, B.; Dubois, M.P.; Miquel, C.; Taberlet, P.; Bonadonna, F. Major histocompatibility complex class II compatibility, but not class I, predicts mate choice in a bird with highly developed olfaction. Proc. R. Soc. B Biol. Sci. 2012, 279, 4457–4463. [Google Scholar] [CrossRef]
  22. Juola, F.A.; Dearborn, D.C. Sequence-based evidence for major histocompatibility complex-disassortative mating in a colonial seabird. Proc. R. Soc. B Biol. Sci. 2012, 279, 153–162. [Google Scholar] [CrossRef]
  23. Reichard, M.; Spence, R.; Bryjova, A.; Bryja, J.; Smith, C. Female Rose Bitterling Prefer MHC-Dissimilar Males: Experimental Evidence. PLoS ONE 2012, 7. [Google Scholar] [CrossRef] [PubMed]
  24. Bahr, A.; Sommer, S.; Mattle, B.; Wilson, A.B. Mutual mate choice in the potbellied seahorse (Hippocampus abdominalis). Behav. Ecol. 2012, 23, 869–878. [Google Scholar] [CrossRef]
  25. Schwensow, N.; Fietz, J.; Dausmann, K.; Sommer, S. MHC-associated mating strategies and the importance of overall genetic diversity in an obligate pair-living primate. Evol. Ecol. 2008, 22, 617–636. [Google Scholar] [CrossRef]
  26. Setchell, J.; Charpentier, M.; Abbott, K.; Wickings, E.; Knapp, L. Opposites attract: MHC-associated mate choice in a polygynous primate. J. Evol. Biol. 2010, 23, 136–148. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Reusch, T.B.H.; Haberli, M.; Aeschlimann, P.B.; Milinski, M. Female sticklebacks count alleles in a strategy of sexual selection explaining MHC polymorphism. Nature 2001, 414, 300–302. [Google Scholar] [CrossRef] [PubMed]
  28. Aeschlimann, P.B.; Haberli, M.A.; Reusch, T.B.H.; Boehm, T.; Milinski, M. Female sticklebacks Gasterosteus aculeatus use self-reference to optimize MHC allele number during mate selection. Behav. Ecol. Sociobiol. 2003, 54, 119–126. [Google Scholar] [CrossRef]
  29. Arkush, K.D.; Giese, A.R.; Mendonca, H.L.; McBride, A.M.; Marty, G.D.; Hedrick, P.W. Resistance to three pathogens in the endangered winter-run chinook salmon (Oncorhynchus tshawytscha): Effects of inbreeding and major histocompatibility complex genotypes. Can. J. Fish. Aquat. Sci. 2002, 59, 966–975. [Google Scholar] [CrossRef]
  30. Meagher, S.; Penn, D.J.; Potts, W.K. Male-Male Competition Magnifies Inbreeding Depression in Wild House Mice. Proc. Natl. Acad. Sci. USA 2000, 97, 3324–3329. [Google Scholar] [CrossRef]
  31. MartinVilla, J.M.; Luque, I.; MartinezQuiles, N.; Corell, A.; Regueiro, J.R.; Timon, M.; ArnaizVillena, A. Diploid expression of human leukocyte antigen class I and class II molecules on spermatozoa and their cyclic inverse correlation with inhibin concentration. Biol. Reprod. 1996, 55, 620–629. [Google Scholar] [CrossRef] [Green Version]
  32. Martin-Villa, J.M.; Longas, J.; Arnaiz-Villena, A. Cyclic expression of HLA class I and II molecules on the surface of purified human spermatozoa and their control by serum inhibin B levels. Biol. Reprod. 1999, 61, 1381–1386. [Google Scholar] [CrossRef]
  33. Paradisi, R.; Neri, S.; Pession, A.; Magrini, E.; Bellavia, E.; Ceccardi, S.; Flamigni, C. Human leukocyte antigen II expression in sperm cells: Comparison between fertile and infertile men. Arch. Androl. 2000, 45, 203–213. [Google Scholar] [PubMed]
  34. Rulicke, T.; Chapuisat, M.; Homberger, F.R.; Macas, E.; Wedekind, C. MHC-genotype of progeny influenced by parental infection. Proc. R. Soc. B Biol. Sci. 1998, 265, 711–716. [Google Scholar] [CrossRef] [PubMed]
  35. Knapp, L.A.; Ha, J.C.; Sackett, G.P. Parental MHC antigen sharing and pregnancy wastage in captive pigtailed macaques. J. Reprod. Immunol. 1996, 32, 73–88. [Google Scholar] [CrossRef]
  36. Ober, C.; Weitkamp, L.R.; Cox, N.; Dytch, H.; Kostyu, D.; Elias, S. HLA and mate choice in humans. Am. J. Hum. Genet. 1997, 61, 497–504. [Google Scholar] [CrossRef]
  37. Schacter, B.; Weitkamp, L.R.; Johnson, W.E. Parental HLA compatibility, fetal wastage and neural-tube defects-evidence for a T/T-linke locus in humans. Am. J. Hum. Genet. 1984, 36, 1082–1091. [Google Scholar]
  38. Lovlie, H.; Gillingham, M.A.F.; Worley, K.; Pizzari, T.; Richardson, D.S. Cryptic female choice favours sperm from major histocompatibility complex-dissimilar males. Proc. R. Soc. B Biol. Sci. 2013, 280, 20131296. [Google Scholar] [CrossRef]
  39. Yeates, S.; Einum, S.; Fleming, I.; Megens, H.-J.; Hindar, K.; Holt, W.; Van Look, K.; Gage, M. Atlantic salmon eggs favour sperm in competition that have similar major histocompatibility alleles. Proc. R. Soc. Lond. B Biol. Sci. 2009, 276, 559–566. [Google Scholar] [CrossRef]
  40. Gessner, C.; Nakagawa, S.; Zavodna, M.; Gemmell, N.J. Sexual selection for genetic compatibility: The role of the major histocompatibility complex on cryptic female choice in Chinook salmon (Oncorhynchus tshawytscha). Heredity 2017, 118. [Google Scholar] [CrossRef]
  41. Gasparini, C.; Congiu, L.; Pilastro, A. Major histocompatibility complex similarity and sexual selection: Different does not always mean attractive. Mol. Ecol. 2015, 24, 4286–4295. [Google Scholar] [CrossRef]
  42. Lenz, T.L.; Hafer, N.; Samonte, I.E.; Yeates, S.E.; Milinski, M. Cryptic haplotype-specific gamete selection yields offspring with optimal MHC immune genes. Evolution 2018, 72, 2478–2490. [Google Scholar] [CrossRef]
  43. Hu, J.C. Research on the Giant Panda; Shanghai Publishing House of Science and Technology: Shanghai, China, 2001. [Google Scholar]
  44. Martin-Wintle, M.S.; Shepherdson, D.; Zhang, G.; Zhang, H.; Li, D.; Zhou, X.; Li, R.; Swaisgood, R.R. Free mate choice enhances conservation breeding in the endangered giant panda. Nat. Commun. 2015, 6. [Google Scholar] [CrossRef] [PubMed]
  45. Zhang, H.M.; Wang, P.Y. The Study on Reproduction in Giant Panda; China FORESTRY publishing House: Bei Jing, China, 2003. [Google Scholar]
  46. Feng, W.H.; Zhang, A.J. The Research on Breeding and Disease of Giant Panda; Sichuan Publishing House of Science and Technology: Cheng Du, China, 1991. [Google Scholar]
  47. Huchard, E.; Baniel, A.; Schliehe-Diecks, S.; Kappeler, P.M. MHC-disassortative mate choice and inbreeding avoidance in a solitary primate. Mol. Ecol. 2013, 22, 4071–4086. [Google Scholar] [CrossRef] [PubMed]
  48. Babik, W. Methods for MHC genotyping in non-model vertebrates. Mol. Ecol. Resour. 2010, 10, 237–251. [Google Scholar] [CrossRef] [PubMed]
  49. Yu, L.; Nie, Y.; Yan, L.; Hu, Y.; Wei, F. No evidence for MHC-based mate choice in wild giant pandas. Ecol. Evol. 2018, 8, 8642–8651. [Google Scholar] [CrossRef] [PubMed]
  50. Wan, Q.H.; Zeng, C.J.; Ni, X.W.; Pan, H.J.; Fang, S.G. Giant panda genomic data provide insight into the birth-and-death process of mammalian major histocompatibility complex class II genes. PLoS ONE 2009, 4, e4147. [Google Scholar] [CrossRef] [PubMed]
  51. Wan, Q.H.; Zhang, P.; Ni, X.W.; Wu, H.L.; Chen, Y.Y.; Kuang, Y.Y.; Ge, Y.F.; Fang, S.G. A novel HURRAH protocol reveals high numbers of monomorphic MHC class II loci and two asymmetric multi-locus haplotypes in the Père David’s deer. PLoS ONE 2011, 6, e14518. [Google Scholar] [CrossRef]
  52. Zhu, Y.; Sun, D.D.; Ge, Y.F.; Yu, B.; Chen, Y.Y.; Wan, Q.H. Isolation and characterization of class I MHC genes in the giant panda (Ailuropoda melanoleuca). Chin. Sci. Bull. 2012, 57, 1–8. [Google Scholar] [CrossRef]
  53. Xie, Z.; Gipps, J. The 2009 International Studbood for Giant Panda (Ailuropoda melanoleuca); Chinese Association of Zoological Garden: Beijing, China, 2009. [Google Scholar]
  54. Durrant, B.S.; Olson, M.A.; Amodeo, D.; Anderson, A.; Russ, K.D.; Campos-Morales, R.; Gual-Sill, F.; Garza, J.R. Vaginal cytology and vulvar swelling as indicators of impending estrus and ovulation in the giant panda (Ailuropoda melanoleuca). Zoo Biol. 2003, 22, 313–321. [Google Scholar] [CrossRef]
  55. Owen, M.A.; Swaisgood, R.R.; McGeehan, L.; Zhou, X.P.; Lindburg, D.G. Dynamics of Male-Female Multimodal Signaling Behavior across the Estrous Cycle in Giant Pandas (Ailuropoda melanoleuca). Ethology 2013, 119, 869–880. [Google Scholar] [CrossRef]
  56. Wan, Q.H.; Zhu, L.; Wu, H.; Fang, S.G. Major histocompatibility complex class II variation in the giant panda (Ailuropoda melanoleuca). Mol. Ecol. 2006, 15, 2441–2450. [Google Scholar] [CrossRef]
  57. Li, D.; Cui, H.; Wang, C.; Ling, S.; Huang, Z.; Zhang, H. A fast and effective method to perform paternity testing for Wolong giant pandas. Chin. Sci. Bull. 2011, 56, 2559–2564. [Google Scholar] [CrossRef] [Green Version]
  58. Zhang, H.M.; Guo, Y.; Li, D.S.; Wang, P.Y.; Fang, S.G. Sixteen novel microsatellite loci developed for the giant panda (Ailuropoda melanoleuca). Conserv. Genet. 2009, 10, 589–592. [Google Scholar] [CrossRef]
  59. Chen, Y.Y.; Zhu, Y.; Wan, Q.H.; Lou, J.K.; Li, W.J.; Ge, Y.F.; Fang, S.G. Patterns of adaptive and neutral diversity identify the Xiaoxiangling Mountains as a refuge for the giant panda. PLoS ONE 2013. [Google Scholar] [CrossRef]
  60. Zeng, C.J.; Pan, H.J.; Gong, S.B.; Yu, J.Q.; Wan, Q.H.; Fang, S.G. Giant panda BAC library construction and assembly of a 650-kb contig spanning major histocompatibility complex class II region. BMC Genom. 2007, 8, 315. [Google Scholar] [CrossRef]
  61. Coltman, D.; Slate, J. Microsatellite measures of inbreeding: A meta-analysis. Evolution 2003, 57, 971–983. [Google Scholar] [CrossRef]
  62. Wetton, J.H.; Carter, R.E.; Parkin, D.T.; Walters, D. Demographic study of a wild house sparrow population by DNA fingerprinting. Nature 1987, 327, 147–149. [Google Scholar] [CrossRef]
  63. Landry, C.; Garant, D.; Duchesne, P.; Bernatchez, L. ‘Good genes as heterozygosity’: The major histocompatibility complex and mate choice in Atlantic salmon (Salmo salar). Proc. R. Soc. Lond. B Biol. Sci. 2001, 268, 1279–1285. [Google Scholar] [CrossRef]
  64. Sandberg, M.; Eriksson, L.; Jonsson, J.; Sjöström, M.; Wold, S. New chemical descriptors relevant for the design of biologically active peptides. A multivariate characterization of 87 amino acids. J. Med. Chem. 1998, 41, 2481–2491. [Google Scholar] [CrossRef]
  65. Huchard, E.; Weill, M.; Cowlishaw, G.; Raymond, M.; Knapp, L.A. Polymorphism, haplotype composition, and selection in the Mhc-DRB of wild baboons. Immunogenetics 2008, 60, 585–598. [Google Scholar] [CrossRef]
  66. Hill, G.E. Plumage coloration is a sexually selected indicator of male quality. Nature 1991, 350, 337–339. [Google Scholar] [CrossRef]
  67. Kurtz, J.; Wegner, K.M.; Kalbe, M.; Reusch, T.B.; Schaschl, H.; Hasselquist, D.; Milinski, M. MHC genes and oxidative stress in sticklebacks: An immuno-ecological approach. Proc. Rroc. Soc. B Biol. Sci. 2006, 273, 1407–1414. [Google Scholar] [CrossRef] [PubMed]
  68. Bjorkman, P.J.; Saper, M.A.; Samraoui, B.; Bennett, W.S.; Strominger, J.L.; Wiley, D.C. The foreigh antigen bing site and T cell recognition regions of class I histocompatibility antigens. Nat. Immunol. 1987, 329, 512–518. [Google Scholar]
  69. Queller, D.C.; Goodnight, K.F. Estimating relatedness using genetic markers. Evolution 1989, 258–275. [Google Scholar] [CrossRef] [PubMed]
  70. Hardy, O.J.; Vekemans, X. SPAGeDi: A versatile computer program to analyse spatial genetic structure at the individual or population levels. Mol. Ecol. Notes 2002, 2, 618–620. [Google Scholar] [CrossRef]
  71. Benjamini, Y.; Hochberg, Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. J. R. Stat. Soc. Ser. B Stat. Methodol. 1995, 57, 289–300. [Google Scholar] [CrossRef]
  72. Yekutieli, D.; Benjamini, Y. Resampling-based false discovery rate controlling multiple test procedures for correlated test statistics. J. Stat. Plan. Infer. 1999, 82, 171–196. [Google Scholar] [CrossRef]
  73. McGraw, K.J.; Hill, G.E. Differential effects of endoparasitism on the expression of carotenoid-and melanin-based ornamental coloration. Proc. R. Soc. Lond. B Biol. Sci. 2000, 267, 1525–1531. [Google Scholar] [CrossRef] [PubMed]
  74. Fitze, P.S.; Richner, H. Differential effects of a parasite on ornamental structures based on melanins and carotenoids. Behav. Ecol. 2002, 13, 401–407. [Google Scholar] [CrossRef]
  75. Siefferman, L.; Hill, G.E. Structural and melanin coloration indicate parental effort and reproductive success in male eastern bluebirds. Behav. Ecol. 2003, 14, 855–861. [Google Scholar] [CrossRef] [Green Version]
  76. Evans, J.; Magurran, A. Multiple benefits of multiple mating in guppies. Proc. Natl. Acad. Sci. USA 2000, 97, 10074–10076. [Google Scholar] [CrossRef] [Green Version]
  77. Colegrave, N.; Kotiaho, J.S.; Tomkins, J.L. Mate choice or polyandry: Reconciling genetic compatibility and good genes sexual selection. Evol. Ecol. Res. 2002, 4, 911–917. [Google Scholar]
  78. Roberts, S.C.; Gosling, L.M. Genetic similarity and quality interact in mate choice decisions by female mice. Nat. Genet. 2003, 35, 103–106. [Google Scholar] [CrossRef] [PubMed]
  79. Lopez, P.; Martin, J. Pheromonal recognition of females takes precedence over the chromatic cue in male Iberian wall lizards Podarcis hispanica. Ethology 2001, 107, 901–912. [Google Scholar] [CrossRef]
  80. Candolin, U. The use of multiple cues in mate choice. Biol. Rev. 2003, 78, 575–595. [Google Scholar] [CrossRef] [Green Version]
  81. Huchard, E.; Knapp, L.A.; Wang, J.; Raymond, M.; Cowlishaw, G. MHC, mate choice and heterozygote advantage in a wild social primate. Mol. Ecol. 2010, 19, 2545–2561. [Google Scholar] [CrossRef] [PubMed]
  82. Feng, W.H.; Wang, R.L.; Zhong, S.M.; Ye, Z.Y.; Cui, X.Z.; Zeng, J.H. Analysis on the Dead Cause of the Anatomical Carcass of Giant Panda (Ailuropoda Melanoleuca); Sichuan Scientific & Technical Publishers: Chengdu, China, 1991; pp. 244–248. [Google Scholar]
  83. Ye, Z.Y. The Control of the Diseases of Giant Panda in Field: Report of 50 Cases; Sichuan Scientific & Technical Publishers: Chengdu, China, 1991. [Google Scholar]
  84. Mainka, S.A.; Qiu, X.M.; He, T.M.; Appel, M.J. Serologic survey of giant pandas (Ailuropoda melanoleuca), and domestic dogs and cats in the Wolong Reserve, China. J. Wildl. Dis. 1994, 30, 86–89. [Google Scholar] [CrossRef]
  85. Qin, Q.; Li, D.S.; Zhang, H.M.; Hou, R.; Zhang, Z.H.; Zhang, C.L.; Zhang, J.G.; Wei, F.W. Serosurvey of selected viruses in captive giant pandas ( Ailuropoda melanoleuca) in China. Vet. Microbiol. 2010, 142, 199–204. [Google Scholar] [CrossRef]
  86. Fedorka, K.M.; Mousseau, T.A. Material and genetic benefits of female multiple mating and polyandry. Anim. Behav. 2002, 64, 361–367. [Google Scholar] [CrossRef]
  87. Yasui, Y. A" good-sperm" model can explain the evolution of costly multiple mating by females. Am. Nat. 1997, 149, 573–584. [Google Scholar] [CrossRef]
  88. Schmidt, C.M.; Orr, H.T. Maternal/fetal interactions: The role of the MHC class I molecule HLA-G. Crit. Rev. Immunol. 1992, 13, 207–224. [Google Scholar]
  89. Rouas-Freiss, N.; GoncAlves, R.M.-B.; Menier, C.; Dausset, J.; Carosella, E.D. Direct evidence to support the role of HLA-G in protecting the fetus from maternal uterine natural killer cytolysis. Proc. Natl. Acad. Sci. USA 1997, 94, 11520–11525. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Cells 08 00257 g001 550
Figure 1. Mean individual major histocompatibility complex heterozygosity (Hobs) of males involved in natural mating and randomly assigned males. (a) five super haplotypes; and, (b) seven individual loci. Two-tailed 95% confidence intervals are indicated by black lines.
Figure 1. Mean individual major histocompatibility complex heterozygosity (Hobs) of males involved in natural mating and randomly assigned males. (a) five super haplotypes; and, (b) seven individual loci. Two-tailed 95% confidence intervals are indicated by black lines.
Cells 08 00257 g001
Cells 08 00257 g002 550
Figure 2. Proportions of heterozygotes and homozygotes in males involved in natural mating and natural non-mating males. (a) five super haplotypes; and, (b) seven individual loci. Asterisk shows results with p values that are smaller than 0.05.
Figure 2. Proportions of heterozygotes and homozygotes in males involved in natural mating and natural non-mating males. (a) five super haplotypes; and, (b) seven individual loci. Asterisk shows results with p values that are smaller than 0.05.
Cells 08 00257 g002
Cells 08 00257 g003 550
Figure 3. Mean allele sharing values (Nas) in females and males involved in natural mating and in females and randomly assigned males. (a) five super haplotypes; and, (b) seven individual loci. Two-tailed 95% confidence intervals are indicated by black lines.
Figure 3. Mean allele sharing values (Nas) in females and males involved in natural mating and in females and randomly assigned males. (a) five super haplotypes; and, (b) seven individual loci. Two-tailed 95% confidence intervals are indicated by black lines.
Cells 08 00257 g003
Cells 08 00257 g004 550
Figure 4. Proportion of allele sharing (Nas) in naturally mated pairs and non-mating pairs (a) five super haplotypes, and (b) seven individual loci. Asterisk shows results with p values that are smaller than 0.05.
Figure 4. Proportion of allele sharing (Nas) in naturally mated pairs and non-mating pairs (a) five super haplotypes, and (b) seven individual loci. Asterisk shows results with p values that are smaller than 0.05.
Cells 08 00257 g004
Cells 08 00257 g005 550
Figure 5. Functional amino acid distances (Faadis) in naturally mated pairs and non-mating pairs. (a) functional amino acid distance of all sites; and, (b) functional amino acid distance of antigen binding sites (ABSs). Asterisk shows results with p values that are smaller than 0.05.
Figure 5. Functional amino acid distances (Faadis) in naturally mated pairs and non-mating pairs. (a) functional amino acid distance of all sites; and, (b) functional amino acid distance of antigen binding sites (ABSs). Asterisk shows results with p values that are smaller than 0.05.
Cells 08 00257 g005
Cells 08 00257 g006 550
Figure 6. Relationship between shared major histocompatibility complex (MHC) alleles and frequency of successful natural mating. (a) SuHa, including four MHC I loci (Aime-C, Aime-F, Aime-I, and Aime-L) and six MHC II loci; (b) SuHaII, including six MHC II loci (DQA1, DQA2, DQB1, DQB2, DRB3, and DRA); (c) DQ, including DQB1, DQB2, DQA1, and DQA2; and, (d) DR, including DRB3 and DRA; (e–h) Four individual loci.
Figure 6. Relationship between shared major histocompatibility complex (MHC) alleles and frequency of successful natural mating. (a) SuHa, including four MHC I loci (Aime-C, Aime-F, Aime-I, and Aime-L) and six MHC II loci; (b) SuHaII, including six MHC II loci (DQA1, DQA2, DQB1, DQB2, DRB3, and DRA); (c) DQ, including DQB1, DQB2, DQA1, and DQA2; and, (d) DR, including DRB3 and DRA; (e–h) Four individual loci.
Cells 08 00257 g006
Cells 08 00257 g007 550
Figure 7. Frequency distributions of functional amino acid distances (Faadis) between eggs and randomly assigned sperm after 10,000 simulations. (a) whole exon2 of Aime-DQA1; (b) antigen binding site (ABS) within exon2 of Aime-DQA1; (c) whole exon2 of Aime-DQA2; and, (d) ABS within exon2 of Aime-DQA2. Two-tailed 95% confidence intervals are indicated by black dashed lines. MHC, major histocompatibility complex.
Figure 7. Frequency distributions of functional amino acid distances (Faadis) between eggs and randomly assigned sperm after 10,000 simulations. (a) whole exon2 of Aime-DQA1; (b) antigen binding site (ABS) within exon2 of Aime-DQA1; (c) whole exon2 of Aime-DQA2; and, (d) ABS within exon2 of Aime-DQA2. Two-tailed 95% confidence intervals are indicated by black dashed lines. MHC, major histocompatibility complex.
Cells 08 00257 g007
Cells 08 00257 g008 550
Figure 8. Functional amino acid distance (Faadis) for genetic compatibility at seven major histocompatibility complex loci between observed zygotes and other zygotes. (a) functional amino acid distance of all sites; and, (b) functional amino acid distance of antigen binding sites (ABSs). Asterisk shows results with p values that are smaller than 0.05.
Figure 8. Functional amino acid distance (Faadis) for genetic compatibility at seven major histocompatibility complex loci between observed zygotes and other zygotes. (a) functional amino acid distance of all sites; and, (b) functional amino acid distance of antigen binding sites (ABSs). Asterisk shows results with p values that are smaller than 0.05.
Cells 08 00257 g008
Cells 08 00257 g009 550
Figure 9. Functional amino acid distances (Faadis) at seven major histocompatibility complex loci between mothers to observed zygotes (offspring) and between mothers to other zygotes. (a) functional amino acid distances at all sites; and, (b) functional amino acid distances at antigen binding sites (ABSs). Asterisk shows results with p values that are smaller than 0.05.
Figure 9. Functional amino acid distances (Faadis) at seven major histocompatibility complex loci between mothers to observed zygotes (offspring) and between mothers to other zygotes. (a) functional amino acid distances at all sites; and, (b) functional amino acid distances at antigen binding sites (ABSs). Asterisk shows results with p values that are smaller than 0.05.
Cells 08 00257 g009
Table
Table 1. Microsatellite results from randomization and paired Student t-tests.
Table 1. Microsatellite results from randomization and paired Student t-tests.
Compatibility TestHeterozygote Advantage Test
RelatednessMLHd2
RandomizationObserved−0.0870.75125.961
Test aSimulated mean−0.0650.74427.491
95% CI[−0.090, −0.040][0.713, 0.775][24.319, 30.648]
P0.0790.6510.345
Paired test bPotential fathers−0.0930.78225.157
Other males−0.0610.74127.480
Statisticst = 1.592Z = 0.556t = 0.913
P0.1150.5780.364
a Comparison between the natural mating group and the randomly assigned group. b Comparison between the natural mating group and the non-mating group.
Table
Table 2. Functional amino acid distances for female genetic compatibility between the observed pairs and randomly assigned pairs.
Table 2. Functional amino acid distances for female genetic compatibility between the observed pairs and randomly assigned pairs.
LocusRegionSimulated Mean [95% CI]Observe MeanP
SuHaABS123.873 [121.759, 125.899]127.8320.001
ALL136.412 [134.088, 138.677]140.7280.000
SuHaIABS93.169 [90.939, 95.359]94.3260.297
ALL101.310 [98.883, 103.635]102.2850.416
SuHaIIABS73.374 [71.216, 75.579]78.0330.000
ALL82.600 [80.143, 85.028]87.4620.000
DQABS53.815 [51.715, 55.876]58.0000.000
ALL57.319 [55.059, 59.577]61.4330.001
DRABS45.051 [43.434, 46.662]47.4980.002
ALL53.531 [51.588, 55.516]56.5100.003
CABS60.340 [58.005, 62.674]62.8400.036
ALL63.839 [61.359, 66.265]66.4310.038
IABS46.787 [44.406, 49.218]45.2180.196
ALL50.388 [47.956, 52.819]48.6810.172
LABS33.315 [31.085, 35.588]32.6240.548
ALL39.862 [37.267, 42.449]38.5840.340
DQA1ABS32.307 [30.740, 33.938]35.5870.000
ALL35.121 [30.740,33.938]38.4000.000
DQA2ABS11.596 [10.674, 12.542]13.0760.001
ALL11.603 [10.674, 12.542]13.0760.001
DQB1ABS32.821 [30.624, 35.005]36.0020.005
ALL34.555 [32.231, 36.817]37.9240.003
DRB3ABS45.051 [43.434, 46.662]47.4980.002
ALL53.531 [51.588, 55.516]56.5100.003
Note: Bolded figures indicate a significant difference after false discovery rate correction. ABS, antigen binding site.
Table
Table 3. Mean functional amino acid distances at five super haplotypes in naturally mated pairs and non-mating pairs.
Table 3. Mean functional amino acid distances at five super haplotypes in naturally mated pairs and non-mating pairs.
Locus Natural Mated PairsNon-Mating PairsStatisticsP
SuHaABS126.209119.677t = 3.3100.001
ALL138.990132.018t = 3.2980.001
SuHaIABS92.31590.284t = 1.0540.294
ALL100.21598.262t = 0.9320.353
SuHaIIABS77.54870.363t = 3.9600.000
ALL86.82879.532t = 3.6300.000
DQABS57.96150.322Z = −4.6040.000
ALL61.37553.769Z = −4.4350.000
DRABS46.87944.284t = 2.2980.024
ALL55.73252.649t = 2.2350.027
Note: Bolded figures indicate a significant difference after false discovery rate correction. ABS, antigen binding site.
Table
Table 4. Association between natural mating success and major histocompatibility complex functional amino acid distances of females and males.
Table 4. Association between natural mating success and major histocompatibility complex functional amino acid distances of females and males.
ABSALL
FPFP
SuHa9.0040.0038.5180.004
SuHaI1.2410.2660.9190.338
SuHaII10.9330.0018.6340.003
DQ11.6490.00110.110.002
DR3.330.0682.8610.091
C3.8490.0503.4980.062
I0.4060.5240.4750.491
L0.1220.7270.4690.494
DQA110.9150.0019,4650.002
DQA24.6880.0314.6880.031
DQB16.9530.0096.880.009
DRB33.330.0682.8610.091
Note: There were 954 degrees of freedom for each genetic variable tested. Bolded figures indicate a significant difference after false discovery rate correction. ABS, antigen binding site.
Table
Table 5. Mean functional amino acid distances at five super haplotypes in observed zygotes and other zygotes.
Table 5. Mean functional amino acid distances at five super haplotypes in observed zygotes and other zygotes.
Locus Observed ZygotesOther ZygotesStatisticsP
SuHaABS31.11931.153Z = −1.1840.237
ALL34.05634.383t = −0.2280.820
SuHaIABS22.29023.649Z = −0.0960.924
ALL24.08225.684Z = −0.0270.978
SuHaIIABS19.95918.785Z = −1.6960.090
ALL22.13621.261Z = −1.1080.268
DQABS15.16413.717Z = −2.6240.009
ALL15.99714.677Z = −2.5260.012
DRABS11.85711.633Z = −0.3040.761
ALL13.92413.943Z = −0.0230.982
Note: Bolded figures indicate a significant difference after false discovery rate correction. ABS, antigen binding site.
Table
Table 6. Association between breeding success and major histocompatibility complex functional amino acid distances between egg and sperm haplotypes.
Table 6. Association between breeding success and major histocompatibility complex functional amino acid distances between egg and sperm haplotypes.
ABSALL
dfFPFP
SuHa4720.6710.4130.9890.320
SuHaI4671.8270.1772.0720.151
SuHaII4530.5310.4660.1200.730
DQ4530.0010.98212.320.000
DR3440.0070.9350.1010.751
C2830.2050.6510.1350.714
I2565.4400.0205.7490.017
L2520.9020.3431.1360.288
DQA13997.0580.0086.3250.012
DQA214616.0900.00016.090.000
DQB12420.1070.7430.0720.788
DRB33440.0070.9350.1010.751
Note: Bolded figures indicate a significant difference after false discovery rate correction. ABS, antigen binding site.

Share and Cite

MDPI and ACS Style

Zhu, Y.; Wan, Q.-H.; Zhang, H.-M.; Fang, S.-G. Reproductive Strategy Inferred from Major Histocompatibility Complex-Based Inter-Individual, Sperm-Egg, and Mother-Fetus Recognitions in Giant Pandas (Ailuropoda melanoleuca). Cells 2019, 8, 257. https://doi.org/10.3390/cells8030257

AMA Style

Zhu Y, Wan Q-H, Zhang H-M, Fang S-G. Reproductive Strategy Inferred from Major Histocompatibility Complex-Based Inter-Individual, Sperm-Egg, and Mother-Fetus Recognitions in Giant Pandas (Ailuropoda melanoleuca). Cells. 2019; 8(3):257. https://doi.org/10.3390/cells8030257

Chicago/Turabian Style

Zhu, Ying, Qiu-Hong Wan, He-Min Zhang, and Sheng-Guo Fang. 2019. "Reproductive Strategy Inferred from Major Histocompatibility Complex-Based Inter-Individual, Sperm-Egg, and Mother-Fetus Recognitions in Giant Pandas (Ailuropoda melanoleuca)" Cells 8, no. 3: 257. https://doi.org/10.3390/cells8030257

APA Style

Zhu, Y., Wan, Q. -H., Zhang, H. -M., & Fang, S. -G. (2019). Reproductive Strategy Inferred from Major Histocompatibility Complex-Based Inter-Individual, Sperm-Egg, and Mother-Fetus Recognitions in Giant Pandas (Ailuropoda melanoleuca). Cells, 8(3), 257. https://doi.org/10.3390/cells8030257

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop