Oestrous cycles and the breeding season of the Père David's deer hind (Elaphurus davidianus)

Reproductive cycles were studied in a group of tame Père David's deer hinds. The non-pregnant hind is seasonally polyoestrous and, in animals studied over 2 years, the breeding season began in early August (2 August +/- 3.3 days; s.e.m., N = 9) and ended in mid-December (18 December +/- 5.7 days; N = 8) and early January (6 January +/- 3.2 days; N = 11) in consecutive years. During the anoestrous period, plasma progesterone concentrations were low (0.2 +/- 0.01 ng/ml) or non-detectable. There was a small, transient increase in progesterone values before the onset of the first cycle of the breeding season. In daily samples taken during an oestrous cycle in which hinds were mated by a marked vasectomized stag, progesterone concentrations remained low (less than 0.5 ng/ml) for a period of about 6 days around the time of oestrus, showed a significant increase above oestrous levels by Day 4 (Day 0 = day of oestrus) and then continued to increase for 18 +/- 2.8 days to reach mean maximum luteal levels of 3.5 +/- 0.6 ng/ml. The plasma progesterone profiles from a number of animals indicated that marking of the hinds by the vasectomized stag did not occur at each ovulation during the breeding season and therefore an estimate of the cycle length could not be determined by this method. In the following year, detection of oestrus in 5 hinds was based on behavioural observations made in the absence of the stag. A total of 19 oestrous cycles with a mean length of 19.5 +/- 0.6 days was observed.(ABSTRACT TRUNCATED AT 250 WORDS)     

Pulsatile secretion of LH during the periovulatory and luteal phases of the oestrous cycle in the Père David's deer hind (Elaphurus davidianus)

Changes in the secretion of LH during the oestrous cycle were studied in 5 tame Père David's deer in which ovulation was synchronized with progesterone implants and prostaglandin injections. Plasma LH concentrations were measured in samples collected at 15-min intervals for a 36-h period, starting 16 h after the removal of the progesterone implants (follicular phase), and for a further 10-h period 10 days after the removal of the progesterone implants (luteal phase). In all animals, there was a preovulatory surge of LH and behavioural oestrus which occurred at a mean time of 59.6 h (+/- 3.25) and 69 h respectively following implant removal. LH pulse frequency was significantly higher during the follicular phase (0.59 +/- 0.03 pulses/h) than the luteal phase (0.24 +/- 0.2 pulses/h), thus confirming in deer findings from research on domesticated ruminants. There were no significant differences between the follicular and luteal phases in mean plasma LH concentrations (0.57 +/- 0.09 and 0.74 +/- 0.13 ng/ml) or mean pulse amplitude (0.99 +/- 0.14 and 1.05 +/- 0.21 ng/ml) for the follicular and luteal phase respectively. The long interval from the removal of progesterone to the onset of the LH surge and the absence of a significant difference in mean LH concentration or pulse amplitude in the follicular and luteal phases resemble published data for cattle but differ from sheep in which there is a short interval from luteal regression to the onset of the surge and a marked increase in LH pulse amplitude during the luteal phase.     

LH secretion and response to GnRH during seasonal anoestrus of the Pére David's deer hind (Elaphurus davidianus)

The pattern of LH secretion and response to exogenous GnRH was determined on 5 occasions during seasonal anoestrus of the Père David's deer hind. LH pulse frequency was low (3.3 +/- 0.6 pulses/18 h) in early anoestrus (February), increased significantly in mid-anoestrus (April; 8.4 +/- 1.4 pulses/18 h) and thereafter declined slightly in late anoestrus (June; 6.3 +/- 0.25 pulses/18 h). Mean LH concentrations also showed significant changes during anoestrus with higher levels in mid-anoestrus (April; 0.85 +/- 0.12 ng/ml) when compared with other times (0.53 +/- 0.05, 0.60 +/- 0.10, 0.68 +/- 0.06 and 0.71 +/- 0.05 ng/ml for February, March, May and June, respectively). LH pulse amplitude showed no significant changes during the study. The LH response to intravenous injections of 2 micrograms GnRH was lowest in early anoestrus (February), increased significantly in mid-anoestrus (April) and remained high through late anoestrus. The response during the luteal phase was similar to that seen during late anoestrus. These results indicate that seasonal anoestrus in the Père David's deer hind is not a uniform state but is characterized by an early period of 'deep' anoestrus.     

Efficacy of intermittent or continuous administration of GnRH in inducing ovulation in early and late seasonal anoestrus in the Père David's deer hind (Elaphurus davidianus)

Père David's deer hinds were treated with GnRH, administered as intermittent i.v. injections (2.0 micrograms/injection at 2-h intervals) for 4 days, or as a continuous s.c. infusion (1.0 micrograms/h) for 14 days. These treatments were given early (February-March) and late (May-June) in the period of seasonal anoestrus. The administration of repeated injections of GnRH increased mean LH concentrations from pretreatment values of 0.54 +/- 0.09 to 2.10 +/- 0.25 ng/ml over the first 8 h of treatment in early anoestrus, and from 0.62 +/- 0.11 to 2.73 +/- 0.49 ng/ml in late anoestrus. The mean amplitude of GnRH-induced LH episodes was greater (P less than 0.01) in late (4.03 +/- 0.28 ng/ml) than in early (3.12 +/- 0.26 ng/ml) anoestrus, but within each replicate (early or late anoestrus), neither mean LH episode amplitude nor mean plasma LH concentrations differed significantly between the four periods of intensive blood sampling. On the basis of their progesterone profiles, 6/12 hinds had ovulated in response to treatment with injections of GnRH (1/6 in early anoestrus and 5/6 in late anoestrus), and oestrus and a preovulatory LH surge were recorded in all of these animals. Oestrus and a preovulatory LH surge were also recorded in one other animal treated in early anoestrus in which progesterone concentrations remained low. The mean times of onset of oestrus (91.0 +/- 1.00 and 62.4 +/- 0.98 h) and of the preovulatory LH surge (85.8 +/- 3.76 and 59.4 +/- 0.25 h) both occurred significantly earlier (P less than 0.001) in animals treated in late anoestrus. Continuous infusion of GnRH to seasonally anoestrous hinds resulted in an increase in mean plasma LH concentrations, but this response did not differ significantly between early (2.15 +/- 0.28 ng/ml) and late (2.48 +/- 0.26 ng/ml) anoestrus. Ovulation, based on progesterone profiles, occurred in 2/7 hinds in early anoestrus and in 4/6 hinds in late anoestrus. Oestrus was detected in all except one of these hinds. The mean time of onset of oestrus occurred earlier in animals treated in late anoestrus (66.2 +/- 0.32 h and 46.7 +/- 0.67 h, P less than 0.01). The administration of GnRH, given either intermittently or continuously, will induce ovulation in a proportion of seasonally anoestrous Père David's deer. However, more animals ovulate in response to this treatment in late than in early anoestrus (75% compared with 23%).     

QTL for live weight traits in Père David's x red deer interspecies hybrids

Interspecies hybrids between Père David's deer (Elaphurus davidianus) and red deer (Cervus elaphus) have proved to be a powerful resource in the search for quantitative trait loci (QTL) in deer. Several regions of the genome with significant effects on live weight and growth rates in backcross hybrids were detected. These include putative QTL for 6-month live weight (LOD 3.90) on linkage group 12, for 14-month live weight (LOD 3.19) on linkage group 1, three putative QTL for growth rate from 3 to 6 months (LOD 4.19 on linkage group 12, LOD 3.92 on linkage group 12, and LOD 3.34 on linkage group 5). In addition, linkage groups 20 and 1 appear to be associated with live weight traits between 9 and 16 months. The variance in traits explained by these QTL ranged between 5.3% and 11.2%. Allele substitution with Père David's alleles at different loci had both positive and negative effects on live weights and growth rates.     

QTL for pubertal and seasonality traits in male Père David's x red deer interspecies hybrids

The unique Pere David's (Elaphurus davidianus) x red deer (Cervus elaphus) backcross hybrid has been used to search for evidence of quantitative trait loci (QTL) for antler pubertal (date and live weight at pedicle initiation) and antler seasonality (date of antler cleaning and casting) traits in temperate species of deer. Analyses using marker information revealed evidence for a QTL for date at pedicle initiation (LOD = 3.7) and live weight at pedicle initiation (LOD = 3.1). These QTL explained 13% and 11% of the phenotypic variance in these traits, respectively.     

Seasonal home range patterns of the reintroduced and rewild female Père David’s deer Elaphurus davidianus

To learn the seasonal home range pattern and its role in the ecological adaptation of rewild Père David’s deer, we took the investigation from spring to winter in three-collared females that released into the coastal wetlands of the Yellow Sea, China. Using GPS collars, minimum convex polygons and fixed kernel estimation, we recorded the movement and estimated the home ranges of the deer. The results showed that the collared deer inhabited only 42.85% of the coastal wetland habitat that had been designated for them. The home ranges were larger during winter and spring than in autumn and summer. In summer and autumn, the collared deer remained in the reserve, which was dominated by spartina grass. In the spring and winter, they were lured by the green leaves of wheat in the farmlands outside the reserve but were then driven back by people. Home range patterns reflected the seasonality of the availability of food resources and human activity. In addition to the frequently utilized habitats, much larger wetland areas are required for such deer. As the population density increases gradually, more available space will be needed to meet the fundamental needs of the seasonal dispersal of rewild Père David’s deer.     

Seasonal changes of reproductive behavior and fecal steroid concentrations in Père David's deer.

From July 1997 to September 1997 and from March 1998 to July 1998, we studied reproductive behaviors of Père David's deer in Dafeng, China. During the field behavioral observations, we collected fresh voided fecal samples from the Père David's stags and hinds periodically and kept those samples under -20 degrees C until laboratory analysis. We analyzed the fecal testosterone, estradiol, and progesterone concentrations in those samples using radioimmunoassay. During this study, we also recorded 17 types of male reproductive behaviors and nine types of female reproductive behaviors. Reproductive behaviors and the fecal steroid concentrations showed overt seasonal fluctuations. There were statistically significant correlations between some male reproductive behaviors, such as anogenital sniffing, urine sniffing, urine spraying, wallowing, bellowing, antler adorning, Antler swags mud, chasing, herding hinds, chin resting, mounting and copulating, with the fecal testosterone concentrations. These results suggested that seasonal reproductive behaviors in stags are strongly associated with circulating testosterone. We also found that some female reproductive behaviors fluctuate corresponding with changes in fecal estrogen concentrations. Although there was no direct evidence to confirm the correlations between female reproductive behaviors and fecal estrogen in our experiment, we could not rule out that reproductive activities of hinds were largely related to ovarian estrogen secretion, and estrogen is necessary for inducing female reproductive behaviors.     

cDNA cloning, recombinant expression and bioactivity of Père David's deer BAFF

B cell activating factor (BAFF), a member of the TNF family, is a critical cytokine for the survival, proliferation, maturation, and differentiation of B cells. In the present study, Père David's deer BAFF (miBAFF) was amplified from Elaphurus davidianus using RT-PCR. This is the first BAFF cloned from a member of Cervidae family. The open reading frame (ORF) of the miBAFF cDNA consists of 843 bases that encode a 280-amino acid protein bearing typical TNF homology domain. Sequence alignment shows that miBAFF shares 39.3%-97% sequence homology with the BAFF sequences of other mammals. Comparative protein modeling predicted that the 3D structure of the soluble mature portion of miBAFF (misBAFF) is very similar to that of human BAFF (hsBAFF). Recombinant misBAFF fused to a SUMO-tag was efficiently expressed in Escherichia coli BL21 (DE3) cells. The protein molecular weight of ~36 KDa was determined using sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) and Western blotting. In vitro, purified misBAFF was shown to promote the survival and proliferation of Père David's deer peripheral blood lymphocytes and mouse B cells. These results indicate that miBAFF plays an important role in the survival/proliferation of mouse B cells and, shows highly conserved evolutionarily, leading to functional cross-reactivity that exists between mouse and Père David's deer BAFF.     

Wherever I may roam—Human activity alters movements of red deer (Cervus elaphus) and elk (Cervus canadensis) across two continents

Human activity and associated landscape modifications alter the movements of animals with consequences for populations and ecosystems worldwide. Species performing long-distance movements are thought to be particularly sensitive to human impact. Despite the increasing anthropogenic pressure, it remains challenging to understand and predict animals' responses to human activity. Here we address this knowledge gap using 1206 Global Positioning System movement trajectories of 815 individuals from 14 red deer (Cervus elaphus) and 14 elk (Cervus canadensis) populations spanning wide environmental gradients, namely the latitudinal range from the Alps to Scandinavia in Europe, and the Greater Yellowstone Ecosystem in North America. We measured individual-level movements relative to the environmental context, or movement expression, using the standardized metric Intensity of Use, reflecting both the directionality and extent of movements. We expected movement expression to be affected by resource (Normalized Difference Vegetation Index, NDVI) predictability and topography, but those factors to be superseded by human impact. Red deer and elk movement expression varied along a continuum, from highly segmented trajectories over relatively small areas (high intensity of use), to directed transitions through restricted corridors (low intensity of use). Human activity (Human Footprint Index, HFI) was the strongest driver of movement expression, with a steep increase in Intensity of Use as HFI increased, but only until a threshold was reached. After exceeding this level of impact, the Intensity of Use remained unchanged. These results indicate the overall sensitivity of Cervus movement expression to human activity and suggest a limitation of plastic responses under high human pressure, despite the species also occurring in human-dominated landscapes. Our work represents the first comparison of metric-based movement expression across widely distributed populations of a deer genus, contributing to the understanding and prediction of animals' responses to human activity.     

Crossing the border? Structure of the red deer (Cervus elaphus) population from the Bavarian–Bohemian forest ecosystem

Anthropogenic impact such as overhunting and habitat fragmentation has reduced the total red deer population (Cervus elaphus) across Europe. In Germany remaining subpopulations are even confined to designated areas with limited or no gene flow among them. Red deer populations inhabiting the Bavarian–Bohemian forest ecosystem had been divided by a fortified State border between Germany and former Czechoslovakia. To assess red deer genetic diversity more than two decades after the removal of the fortifications, we analysed a population from the Bavarian Forest National Park, Germany, and one from the National Park ?umava, Czech Republic, using 11 microsatellite loci and a 910 bp long section of the mitochondrial control region (mtDNA). Bayesian analyses of microsatellite allele frequencies favoured the presence of a single population in the Bavarian-Bohemian forest ecosystem over other population genetic structures. This admixture was supported by a lack of population pairwise differentiation between German and Czech red deer microsatellite genotypes in the analysis of molecular variance (AMOVA, F_ST = 0.009, p = 0.383). Contrastingly, AMOVA revealed a highly significant matrilinear differentiation of mtDNA between the two samples (Φ_ST = 0.285, p = 0.002), whereby German red deer belonged predominantly to haplogroup A (western Europe) and Czech red deer predominantly to haplogroup C (eastern Europe). In combination, these findings indicated a high degree of philopatry by does and extensive gene flow across the former border mediated by stags. They also identified the Bavarian–Bohemian forest ecosystem as part of a suture zone between western and eastern European red deer matrilines.     

Diet composition and selection of Père David's deer in Hubei Shishou Milu National Nature Reserve,China

Hubei Shishou Milu National Nature Reserve is an ideal place to restore the wild population of Père David's deer (Elaphurus davidianus). Understanding foraging ecology and diet composition is essential for assessing population development or establishing long-term effective conservation measures for endangered species. However, little is known about the diet composition of Père David's deer and its diet selection mechanism. In this study, we used stable isotope technology to investigate the diet composition of Père David's deer according to various tissues (i.e., fur, muscle, liver, heart, and feces) and seasons, and evaluated the correlation between the nutrient composition of plants and diet composition. Bayesian isotope analysis showed that the autumn and winter diet estimated by fur and fecal samples indicated a diet dominated by C₃ grasses (42.7%–57.2%, mean), while the summer diet estimated by muscle and liver samples was dominated by C₃ forbs (30.9%–41.6%, mean). The Pearson correlation test indicated that the contribution of winter diet composition reflected by fur and fecal samples was associated with correlations with crude protein (r = .666, p < .01) and soluble sugars (r = .695, p < .01). The results indicated that crude protein and soluble sugars were important factors influencing the winter diet selection of Père David's deer. In the context of the current reintroduction facing many challenges, such as habitat fragmentation, wetland degradation, and human disturbance, comprehensively evaluating the diet selection mechanism of Père David's deer under different resource specificities and temporal changes should be considered in the future.     

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

The Père David's deer is a highly inbred, but recovered, species, making it interesting to consider their adaptive molecular evolution from an immunological perspective. Prior to this study, genomic sequencing was the only method for isolating all functional MHC genes within a certain species. Here, we report a novel protocol for isolating MHC class II loci from a species, and its use to investigate the adaptive evolution of this endangered deer at the level of multi-locus haplotypes. This protocol was designated "HURRAH" based on its various steps and used to estimate the total number of MHC class II loci. We confirmed the validity of this novel protocol in the giant panda and then used it to examine the Père David's deer. Our results revealed that the Père David's deer possesses nine MHC class II loci and therefore has more functional MHC class II loci than the eight genome-sequenced mammals for which full MHC data are currently available. This could potentially account at least in part for the strong survival ability of this species in the face of severe bottlenecking. The results from the HURRAH protocol also revealed that: (1) All of the identified MHC class II loci were monomorphic at their antigen-binding regions, although DRA was dimorphic at its cytoplasmic tail; and (2) these genes constituted two asymmetric functional MHC class II multi-locus haplotypes: DRA1*01 ～ DRB1 ～ DRB3 ～ DQA1 ～ DQB2 (H1) and DRA1*02 ～ DRB2 ～ DRB4 ～ DQA2 ～ DQB1 (H2). The latter finding indicates that the current members of the deer species have lost the powerful ancestral MHC class II haplotypes of nine or more loci, and have instead fixed two relatively weak haplotypes containing five genes. As a result, the Père David's deer are currently at risk for increased susceptibility to infectious pathogens.     

Eight novel microsatellite markers from the Père David’s deer (<Emphasis Type=Italic>Elaphurus davidianus</Emphasis>)

Père David’s deer underwent known bottlenecks but has recovered to more than 2000 individuals in China, making it interesting to assess its genetic variability from a microsatellite perspective. For this purpose, we developed eight novel microsatellite loci using magnetic-bead enrichment methods. These microsatellite markers revealed a low level of genetic variation in the David’s deer, showing two alleles each locus, 0.31 and 0.38 for average observed and expected heterozyogosities, respectively. Combined with previous cross-species primers, this set of markers would play an important role in future genetic investigation of the David’s deer.     

Selective culling of Iberian red deer stags (<Emphasis Type="Italic">Cervus elaphus hispanicus</Emphasis>) by selective montería in Spain

Hunting can affect animal populations not only by increasing mortality but also by introducing selection components associated with particular features of individuals. In addition to the most widespread hunting system in Spain for Iberian red deer stags (Cervus elaphus hispanicus) called montería, there are also selective monterías aimed at culling poor-trophy males in order to improve the average quality of the trophies for commercial hunt. This way of removing poor-trophy males contrasts with the most common procedure of shooting individual males by selective stalking that is used in other areas of Europe. Also, due to the hunting procedure by which most deer are shot while running chased by dogs, it is doubtful whether hunters are actually producing a selective impact on deer populations. In this paper, we compare data of males shot in commercial montería and in selective montería in Southern Spain. We found that males in selective montería were smaller in body size and in antler size than in commercial hunts, even correcting by age, although the selective effect was stronger at some ages. We discuss the implications of this practice for sustainable use and conservation.     

The Père David's deer MHC class I genes show unexpected diversity patterns, with monomorphic classical genes but polymorphic nonclassical genes and pseudogenes

Père David's deer (Elaphurus davidianus) is a highly inbred species that arose from 11 founders but now comprises a population of about 3,000 individuals, making it interesting to investigate the adaptive variation of this species from the major histocompatibility complex (MHC) perspective. In this study, we isolated Elda-MHC class I loci using magnetic bead-based cDNA hybridization, and examined the molecular variations of these loci using single-strand conformation polymorphism (SSCP) and sequence analysis. We obtained seven MHC class I genes, which we designated F1, F12, G2, I7, AF, I8, and C1. Our analyses of stop codons, phylogenetic trees, amino acid conservation, and G+C content revealed that F1, F12, G2, and I7 were classical genes, AF was a nonclassical gene, and I8 and C1 were pseudogenes. Our subsequent molecular examinations showed that the diversity pattern in the Père David's deer was unusual. Most mammals have more polymorphic classical class I loci vs. the nonclassical and neutral genes. In contrast, the Père David's deer was found to be monomorphic at classical genes F1, F12, G2, and I7, dimorphic at the nonclassical AF gene, dimorphic at pseudogene I8, and tetramorphic at pseudogene C1. The adverse polymorphism patterns of Elda-I genes might provide evidence for selection too faster deplete MHC variation than drift in the bottlenecked populations, while the postbottleneck tetramorphism of the C1 pseudogene appears to be evidence of strong historical balancing selection.     

Eighteen novel polymorphic microsatellite loci developed from the Père David’s deer (<Emphasis Type="Italic">Elaphurus davidianus</Emphasis>)

Père David’s deer is a severely bottlenecked species but without showing inbreeding depression, making it essential to develop molecular markers to explore her genetic mechanism of population recovery. In this study, we isolated 18 novel polymorphic microsatellite loci from a dinucleotide-enriched library. This suit of markers presented 2–3 alleles for each locus and their observed and expected heterozygosities were 0.057–0.610 and 0.056–0.598, respectively. These new microsatellite loci had an average of 2.12 alleles and thus contributed to relatively low exclusion probabilities of parentage and paternity testing (0.768 and 0.921). However, when these loci were examined in combination with previous microsatellite markers, overall probabilities of parentage and paternity exclusion went up to 0.905 and 0.990, respectively, showing that these 26 microsatellite loci should be adopted together in future genetic analyses for this highly inbred species.