Conservation genetics focuses on the effects of contemporary genetic structuring on long-term survival of a species. It helps wildlife managers protect biodiversity by identifying a series of conservation units, which include species, evolutionarily significant units (ESUs), management units (MUs), action units (AUs), and family nets (FNs). Although mitochondrial DNA (mtDNA) evolves 5-10 times faster than single-copy nuclear DNA (scnDNA), it records few traces of contemporary events. Thus, mtDNA can be used to resolve taxonomic uncertainties and ESUs. Variable number of tandem repeats (VNTRs) evolve 100-1000 times faster than scnDNA and provide a powerful tool for analyzing recent and contemporary events. VNTR analysis techniques include polymerase chain reaction (PCR)-based microsatellite assays and oligonucleotide probing. Size homoplasy problems in PCR-based microsatellite assays can strongly affect the inference of recent population history. The high homozygosity in endangered species is reflected in a relatively low number and level of variability in microsatellite loci. This combined with "allelic dropout" and "misprinting" errors contributes to the generation of highly biased genetic data following analyses of natural populations. Thus, in conservation genetics, microsatellites are of limited use for identifying ESUs, MUs, and AUs. In contrast to PCR-based microsatellite analysis, oligonucleotide probing avoids errors resulting from PCR amplification. It is particularly suitable for inferring recent population history and contemporary gene flow between fragmented subpopulations. Oligonucleotide fingerprinting generates individual-specific DNA banding patterns and thus provides a highly precise tool for monitoring demography of natural populations. Hence, DNA fingerprinting is powerful for distinguishing ESUs, MUs, AUs, and FNs. The use of oligonucleotide fingerprinting and fecal DNA is opening new areas for conservation genetics.

Crocodilians are diving reptiles that can hold their breath under water for long periods of time and are crepuscular animals with excellent sensory abilities. They comprise a sister lineage of birds and have no sex chromosome. Here we report the genome sequence of the endangered Chinese alligator (Alligator sinensis) and describe its unique features. The next-generation sequencing generated 314 Gb of raw sequence, yielding a genome size of 2.3 Gb. A total of 22 200 genes were predicted in Alligator sinensis using a de novo, homology- and RNA-based combined model. The genetic basis of long-diving behavior includes duplication of the bicarbonate-binding hemoglobin gene, co-functioning of routine phosphate-binding and special bicarbonate-binding oxygen transport, and positively selected energy metabolism, ammonium bicarbonate excretion and cardiac muscle contraction. Further, we elucidated the robust Alligator sinensis sensory system, including a significantly expanded olfactory receptor repertoire, rapidly evolving nerve-related cellular components and visual perception, and positive selection of the night vision-related opsin and sound detection-associated otopetrin. We also discovered a well-developed immune system with a considerable number of lineage-specific antigen-presentation genes for adaptive immunity as well as expansion of the tripartite motif-containing C-type lectin and butyrophilin genes for innate immunity and expression of antibacterial peptides. Multifluorescence in situ hybridization showed that alligator chromosome 3, which encodes DMRT1, exhibits significant synteny with chicken chromosome Z. Finally, population history analysis indicated population admixture 0.60-1.05 million years ago, when the Qinghai-Tibetan Plateau was uplifted.

The major histocompatibility complex (MHC) plays an important role in immune response. Avian MHCs are not well characterized, only reporting highly compact Galliformes MHCs and extensively fragmented zebra finch MHC. We report the first genomic structure of an endangered Pelecaniformes (crested ibis) MHC containing 54 genes in three regions spanning ~500 kb. In contrast to the loose BG (26 loci within 265 kb) and Class I (11 within 150) genomic structures, the Core Region is condensed (17 within 85). Furthermore, this Region exhibits a COL11A2 gene, followed by four tandem MHC class II αβ dyads retaining two suites of anciently duplicated “αβ” lineages. Thus, the crested ibis MHC structure is entirely different from the known avian MHC architectures but similar to that of mammalian MHCs, suggesting that the fundamental structure of ancestral avian class II MHCs should be “COL11A2-IIαβ1-IIαβ2.” The gene structures, residue characteristics, and expression levels of the five class I genes reveal inter-locus functional divergence. However, phylogenetic analysis indicates that these five genes generate a well-supported intra-species clade, showing evidence for recent duplications. Our analyses suggest dramatic structural variation among avian MHC lineages, help elucidate avian MHC evolution, and provide a foundation for future conservation studies.

Habitat destruction and human activity have greatly impacted the natural history of the giant panda (Ailuropoda melanoleuca). Although the genetic diversity of neutral markers has been examined in this endangered species, no previous work has examined adaptive molecular polymorphisms in the giant panda. Here, the major histocompatibility complex (MHC) class II DRB locus was investigated in the giant panda, using single-strand conformation polymorphism (SSCP) and sequence analysis. Comparisons of DNA samples extracted from faecal and blood samples from the same individual revealed that the two materials yielded similar quantities and qualities of DNA, as well as identical SSCP patterns and allelic sequences, demonstrating the reliability of DNA isolation from panda faeces. Analysis of faecal samples from 60 giant pandas revealed relatively low number of alleles: seven alleles. However, the alleles were quite divergent, varying from each other by a range of 7-47 nucleotide substitutions (4-25 amino acid substitutions). Construction of a neighbour-joining tree and comparisons among DRB alleles from other species revealed that both similar and highly divergent alleles survived in the bottlenecked panda populations. Despite species-specific primers used and excellent faecal DNA isolated, a lower level of heterozygosity than expected was still observed due to inbreeding. There were three types of evidence supporting the presence of balancing selection in the giant panda: (i) an obvious excess of nonsynonymous substitutions over synonymous at the antigen-binding positions; (ii) trans-species evolution of two alleles between the giant panda and other felids; and (iii) a more even distribution of alleles than expected from neutrality.

Over the last 100 years giant pandas (Ailuropoda melanoleuca) have been separated into six completely isolated mountain ranges. DNA fingerprinting revealed different differentiation patterns in giant pandas, including early-stage, late-stage, recent divergence and recent genetic depression. A separation around 10 000 years ago resulted in highly significant differences in DNA fingerprints and morphological characteristics between Qinling and Sichuan populations. Supported by morphological differences, the genetic data were used to classify the Qinling population as a new subspecies, A. m. qinlingensis, while the Sichuan populations were classified into the original subspecies, A. m. melanoleuca. Thus, the Qinling population deserves management as a separate unit. In the Sichuan populations, two management units were defined, including Qionglai-Minshan and Daxiangling-Xiaoxiangling-Liangshan. Our data suggest urgent measures are needed to establish green corridors between subpopulations in each mountain range to increase gene flow and genetic variation to ensure long-term survival.

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

As a natural hibernator, the Chinese alligator (Alligator sinensis) is an ideal and intriguing model to investigate changes in microbial community structure and function caused by hibernation. In this study, we used 16S rRNA profiling and metagenomic analysis to compare the composition, diversity, and functional capacity in the gut microbiome of hibernating vs. active Chinese alligators. Our results show that gut microbial communities undergo seasonal restructuring in response to seasonal cycles of feeding and fasting in the Chinese alligator, but this animal harbors a core gut microbial community primarily dominated by Proteobacteria, Fusobacteria, Bacteroidetes, and Firmicutes across the gut regions. During hibernation, there is an increase in the abundance of bacterial taxa (e.g., the genus Bacteroides) that can degrade host mucin glycans, which allows adaptation to winter fasting. This is accompanied by the enrichment of mucin oligosaccharide-degrading enzyme and carbohydrate-active enzyme families. In contrast, during the active phase (feeding), active Chinese alligators exhibit a carnivore gut microbiome dominated by Fusobacteria, and there is an increase in the relative abundance of bacteria (e.g., Cetobacterium somerae) with known proteolytic and amino acids-fermentating functions that improve host protein-rich food digestion efficiency. In addition, seasonal variations in the expression of β-defensins play a protective role in intestinal immunity. These findings provide insights into the functional adaptations of host–gut microbe symbioses to seasonal dietary shifts to maintain gut homeostasis and health, especially in extreme physiological states.

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.