Ancient hybridization and admixture in macaques (genus Macaca) inferred from whole genome sequences.

The evolutionary history of the stump-tailed macaque (Macaca arctoides) and its genetic relationship to other macaques is a subject of continuing controversy. Here, we have reported the first genome sequences of two stump-tailed macaques and one Assamese macaque (M. assamensis). Additionally, we have investigated the genetic diversity between macaque species and analyzed ancient hybridization events. Genome-wide analyses demonstrated that the stump-tailed macaque is more closely related to sinica species than to fascicularis/mulatta species. This topology contradicts the mitochondrial sequence-based phylogeny that places the stump-tailed macaque into the fascicularis/mulatta group. However, our results further show that stump-tailed macaques have genetic backgrounds distinct from sinica species, and present evidence of gene flows with rhesus macaques. We suggest that an ancient introgression occurred after stump-tailed macaques diverged from sinica species. The distinct gene flow between proto-arctoides and proto-mulatta resulted in the transfer of rhesus macaque-type mitochondria into proto-arctoides. The rhesus macaque-type mitochondria remained in populations because of genetic drift during the bottleneck. The PSMC results and morphological and geographic evidence are consistent with the mitochondria capture pattern in the stump-tailed macaque. The molecular clock estimates suggest that the mitochondrial transference into stump-tailed macaques occurred 0.4-1.4 million years ago. Furthermore, we detected extensive admixtures between different macaque species, indicating that gene flow has played an important role in the evolutionary history of the genus Macaca.

Prevalence of gastrointestinal parasites in captive non-human primates of twenty-four zoological gardens in China.

Captive primates are susceptible to gastrointestinal (GIT) parasitic infections, which are often zoonotic and can contribute to morbidity and mortality. Fecal samples were examined by the means of direct smear, fecal flotation, fecal sedimentation, and fecal cultures. Of 26.51% (317/1196) of the captive primates were diagnosed gastrointestinal parasitic infections. Trichuris spp. were the most predominant in the primates, while Entamoeba spp. were the most prevalent in Old World monkeys (P < 0.05). These preliminary data will improve the management of captive primates and the safety of animal keepers and visitors.

Multilocus genotype and subtype analysis of Cryptosporidium andersoni derived from a Bactrian camel (Camelus bactrianus) in China.

Fecal specimens from two Bactrian camels were collected in the Ya'an city zoo of China and were examined for Cryptosporidium by centrifugal flotation. One specimen was found to be parasitized by Cryptosporidium via microscopy, and the oocysts were measured to have an average size of 7.03 × 5.50 µm (n > 50). The isolate was genotyped by polymerase chain reaction (PCR) amplification and DNA sequence analysis of the partial 18S rRNA, COWP, and A135 genes, and was confirmed to be Cryptosporidium andersoni with minor nucleotide differences. Multilocus sequence typing (MLST) analysis indicated that the subtype of the camel-derived C. andersoni isolate was A4, A4, A4, and A1 at the four minisatellite loci (MS1, MS2, MS3, and MS16, respectively). Therefore, this isolate belongs to the most common MLST subtype reported in cattle in China and is distinct from two other known camel C. andersoni MLST subtypes (A6, A4, A2, A1 and A6, A5, A2, A1). Animal transmission experiments demonstrated that the C. andersoni isolate was not infectious to immunosuppressed or immunocompetent Kun-ming mice, Sprague-Dawley rats, and hamsters but was biologically similar to most bovine C. andersoni isolates characterized so far. Therefore, transmission of this camel-derived C. andersoni isolate is very likely to occur between camels and bovine.

Mebendazole in the treatment of Hymenolepis nana infections in the captive ring-tailed lemur (Lemur catta), China.

The present study was conducted to evaluate the effectiveness of mebendazole in the treatment of Hymenolepis nana infection in ring-tailed lemurs (Lemur catta). Ten (L. catta) from the Chengdu Zoological Garden in China, which were naturally infected with H. nana, were treated with mebendazole (10 mg/kg for 5 days). A posttreatment fecal examination was conducted 10 and 20 days after the start of treatment. All treatments resulted in a decrease in the number of eggs per gram in the posttreatment sample compared with the pretreatment sample. Reduction of mean egg count was 97.6% and 100% on days 10 and 20, respectively. The results indicated that mebendazole has marked efficacy against H. nana infections in L. catta.

Comparative efficacy of ivermectin and levamisole for reduction of migrating and encapsulated larvae of Baylisascaris transfuga in mice.

The comparative efficacy of 2 anthelmintics (ivermectin and levamisole) against Baylisascaris transfuga migrating and encapsulated larvae was studied in mice. A total of 60 BALB/c mice inoculated each with about 1,000 embryonated B. transfuga eggs were equally divided into 6 groups (A-F) randomly. Mice of groups A and B were treated with ivermectin and levamisole, respectively, on day 3 post-infection (PI). Mice of groups A-C were killed on day 13 PI. Similarly, groups D and E were treated with ivermectin and levamisole, respectively, on day 14 PI, and all mice of groups D-F were treated on day 24 PI. The groups C and F were controls. Microexamination was conducted to count the larvae recovering from each mouse. The percentages of reduction in the number of migrating larvae recovered from group A (ivermectin) and B (levamisole) were 88.3% and 81.1%, respectively. In addition, the reduction in encapsulated larvae counts achieved by ivermectin (group D) and levamisole (group E) was 75.0% and 49.2%, respectively. The results suggested that, to a certain extent, both anthelmintics appeared to be more effective against migrating larvae than encapsulated larvae. However, in the incipient stage of infection, ivermectin may be more competent than levamisole as a larvicidal drug for B. transfuga.

Comparative Microbiome Analysis Reveals the Ecological Relationships Between Rumen Methanogens, Acetogens, and Their Hosts.

Ruminant methane, which is generated by methanogens through the consumption of hydrogen and supports the normal function of the rumen ecosystem, is a major source of greenhouse gases. Reductive acetogenesis by acetogens is a possible alternative sink that can dispose of hydrogen for acetate production. However, the distribution of rumen methanogens and acetogens along with the relationships among methanogens, acetogens, and their host are poorly understood. Therefore, we investigated the rumen methanogen and acetogen communities of 97 individual animals representing 14 ruminant species within three ruminant families Cervidae (deer), Bovidae (bovid), and Moschidae (musk deer). The results showed that the Methanobrevibacter spp. and acetogens associated with Eubacteriaceae were the most widespread methanogens and acetogens, respectively. However, other methanogens and acetogens exhibited host specificity in the rumen of reindeer and Chinese muntjac deer. Acetogen and methanogen communities were not correlated in these species, and the phylosymbiosis signature between host phylogeny and the composition of both communities was lacking. The abundance of Methanobrevibacter gottschalkii was negatively correlated with the degree of papillation of the rumen wall. Finally, co-occurrence analysis showed that the variation of the predicted methane yields was characterized by the interactive patterns between methanogens, acetogens, and concentrations of rumen metabolites. Our results show that rumen methanogen and acetogen communities have low compositional interdependence and do not exhibit parallel host evolution, which suggests that the strategies for mitigating methane production should be based on a species-specific rumen microbiota analysis.

Can-SINE dynamics in the giant panda and three other Caniformia genomes.

BACKGROUND: Although repeat sequences constitute about 37% of carnivore genomes, the characteristics and distribution of repeat sequences among carnivore genomes have not been fully investigated. Based on the updated Repbase library, we re-annotated transposable elements (TEs) in four Caniformia genomes (giant panda, polar bear, domestic dog, and domestic ferret) and performed a systematic, genome-wide comparison focusing on the Carnivora-specific SINE family, Can-SINEs. RESULTS: We found the majority of young recently integrated transposable elements are LINEs and SINEs in carnivore genomes. In particular, SINEC1_AMe, SINEC1B_AMe and SINEC_C1 are the top three most abundant Can-SINE subfamilies in the panda and polar bear genomes. Transposition in transposition analysis indicates that SINEC1_AMe and SINEC1B_AMe are the most active subfamilies in the panda and the polar bear genomes. SINEC2A1_CF and SINEC1A_CF subfamilies show a higher retrotransposition activity in the dog genome, and MVB2 subfamily is the most active Can-SINE in the ferret genome. As the giant panda is an endangered icon species, we then focused on the identification of panda specific Can-SINEs. With the panda-associated two-way genome alignments, we identified 250 putative panda-specific (PPS) elements (139 SINEC1_AMes and 111 SINEC1B_AMes) that inserted in the panda genome but were absent at the orthologous regions of the other three genomes. Further investigation of these PPS elements allowed us to identify a new Can-SINE subfamily, the SINEC1_AMe2, which was distinguishable from the current SINEC1_AMe consensus by four non-CpG sites. SINEC1_AMe2 has a high copy number (> 100,000) in the panda and polar bear genomes and the vast majority (> 96%) of the SINEC1_AMe2 elements have divergence rates less than 10% in both genomes. CONCLUSIONS: Our results suggest that Can-SINEs show lineage-specific retransposition activity in the four genomes and have an important impact on the genomic landscape of different Caniformia lineages. Combining these observations with results from the COSEG, Network, and target site duplication analysis, we suggest that SINEC1_AMe2 is a young mobile element subfamily and currently active in both the panda and polar bear genomes.

Microbial Biogeography Along the Gastrointestinal Tract of a Red Panda.

The red panda (Ailurus fulgens) is a herbivorous carnivore that is protected worldwide. The gastrointestinal tract (GIT) microbial community has widely acknowledged its vital role in host health, especially in diet digestion; However, no study to date has revealed the GIT microbiota in the red panda. Here, we characterized the microbial biogeographical characteristics in the GIT of a red panda using high-throughput sequencing technology. Significant differences were observed among GIT segments by beta diversity of microbiota, which were divided into four distinct groups: the stomach, small intestine, large intestine, and feces. The stomach and duodenum showed less bacterial diversity, but contained higher bacterial abundance and the most unclassified tags. The number of species in the stomach and small intestine samples was higher than that of the large intestine and fecal samples. A total of 133 core operational taxonomic units were obtained from the GIT samples with 97% sequence identity. Proteobacteria (52.16%), Firmicutes (10.09%), and Bacteroidetes (7.90%) were the predominant phyla in the GIT of the red panda. Interestingly, Escherichia-Shigella were largely abundant in the stomach, small intestine, and feces whereas the abundance of Bacteroides in the large intestine was high. Overall, our study provides a deeper understanding of the gut biogeography of the red panda microbial population. Future research will be important to investigate the microbial culture, metagenomics and metabolism of red panda GIT, especially in Escherichia-Shigella.

Cellulose-dependent expression and antibacterial characteristics of surfactin from Bacillus subtilis HH2 isolated from the giant panda.

Surfactin secreted by Bacillus subtilis can confer strong, diverse antipathogenic effects, thereby benefitting the host. Carbon source is an important factor for surfactin production. However, the mechanism that bacteria utilize cellulose, the most abundant substance in the intestines of herbivores, to produce surfactin remains unclear. Here, we used B. subtilis HH2, isolated from the feces of a giant panda, as a model to determine changes in surfactin expression in the presence of different concentrations of cellulose by quantitative polymerase chain reaction and high-performance liquid chromatography. We further investigated the antimicrobial effects of surfactin against three common intestinal pathogens (Escherichia coli, Staphylococcus aureus, and Salmonella enterica) and its resistance to high temperature (60-121°C), pH (1-12), trypsin (100-300 µg/mL, pH 8), and pepsin (100-300 µg/mL, pH 2). The results showed that the surfactin expressed lowest in bacteria cultured in the presence of 1% glucose medium as the carbon source, whereas increased in an appropriate cellulose concentration (0.67% glucose and 0.33% cellulose). The surfactin could inhibit E. coli and Staphylococcus aureus, but did not affect efficiently for Salmonella enterica. The antibacterial ability of surfactin did not differ according to temperature (60-100°C), pH (2-11), trypsin (100-300 µg/mL), and pepsin (100-300 µg/mL; P > 0.05), but decreased significantly at extreme environments (121°C, pH 1 or 12; P < 0.05) compared with that in the control group (37°C, pH = 7, without any protease). In conclusion, our findings indicated that B. subtilis HH2 could increase surfactin expression in an appropriate cellulose environment and thus provide benefits to improve the intestinal health of herbivores.

Complete mitochondrial genome sequence of the Asian golden cat, Catopuma temminckii.

In this study, the mitochondrial genome of Asian golden cat (Catopuma temminckii) is sequenced. The mitochondrial genome was 16,985 bp long, including 13 protein-coding genes, 22 transfer RNA (tRNA) genes, 2 ribosomal RNA (rRNA) genes, 1 control region and 1 origin of light-strand replication. The overall base composition of the mitochondrial genome was 32.76% A, 27.49 % T, 25.75 % C, and 13.99 % G. The complete mitochondrial genome of Catopuma temminckii could contribute to understanding taxonomic status and phylogenetic relationship of genus Catopuma.

The complete mitochondrial genome of Stump-tailed Macaques (Macaca arctoides).

Stump-tailed macaque (Macaca arctoides) has been enlisted as the Near Threatened species in the IUCN Red List. In this study, the complete mitochondial genome of M. arctoides was determined. The mitogenome was 16,559 bp in length with an A + T content of 56.8%, containing 2 ribosomal RNA genes (rRNA), 22 transfer RNA genes (tRNA), 13 protein-coding genes (PCGs) and a control region (D-loop). Ten protein-coding genes (ND1, COX1, COX2, ND4, COX3, ND4L, ATP8, ATP6, CYTB, ND6) started with ATG codons, while ND2, ND3 and ND5 initiated with ATT, ATC and ATA, respectively. Eight PCGs ended with complete termination codons except for COX3, ND3, ND4, and CYTB terminated with incompleted codon T.

Mitochondrial Genome and Nuclear Markers Provide New Insight into the Evolutionary History of Macaques.

The evolutionary history of macaques, genus Macaca, has been under debate due to the short times of divergence. In this study, maternal, paternal, and biparental genetic systems were applied to infer phylogenetic relationships among macaques and to trace ancient hybridization events in their evolutionary history. Using a PCR display method, 17 newly phylogenetically informative Alu insertions were identified from M. assamensis. We combined presence/absence analysis of 84 Alu elements with mitochondrial genomes as well as nuclear sequences (five autosomal genes, two Y chromosomal genes, and one X chromosomal fragment) to reconstruct a robust macaque phylogeny. Topologies generated from different inherited markers were similar supporting six well defined species groups and a close relationship of M. assamensis and M. thibetana, but differed in the placing of M. arctoides. Both Alu elements and nuclear genes supported that M. arctoides was close to the sinica group, whereas the mitochondrial data clustered it into the fascicularis/mulatta lineage. Our results reveal that a sex-biased hybridization most likely occurred in the evolutionary history of M. arctoides, and suggest an introgressive pattern of male-mediated gene flow from the ancestors of M. arctoides to the M. mulatta population followed by nuclear swamping. According to the estimation of divergence dates, the hybridization occurred around 0.88~1.77 mya (nuclear data) or 1.38~2.56 mya (mitochondrial data). In general, our study indicates that a combination of various molecular markers could help explain complicated evolutionary relationships. Our results have provided new insights into the evolutionary history of macaques and emphasize that hybridization might play an important role in macaque evolution.

Cloning and characterization of giant panda (Ailuropoda melanoleuca) IL-18 binding protein.

The giant panda (Ailuropoda melanoleuca) is an endangered species. Interleukin-18 (IL-18) plays an important role in the innate and adaptive immune responses by inducing IFN-γ. IL-18 has been implicated in the pathogenesis of various diseases. IL-18 binding protein (IL-18BP) is an intrinsic inhibitor of IL-18 that possesses higher affinity to IL-18. In this study, we cloned and characterized IL-18BP in giant panda (AmIL-18BP) from the spleen. The amino acid sequence of giant panda IL-18BP ORF shared about 65% identities with other species. To evaluate the effects of AmIL-18BP on the immune responses, we expressed the recombinant AmIL-18BP in Escherichia coli BL21 (DE3).The fusing protein PET-AmIL-18BP was purified by nickel affinity column chromatography. The biological function of purified PET-AmIL-18BP was determined on mice splenocyte by qRT-PCR. The results showed that AmIL-18BP was functional and could significantly reduce IFN-γ production in murine splenocytes. These results will facilitate the study of protecting giant panda on etiology and immunology.

Multilocus genotypes and broad host-range of Enterocytozoon bieneusi in captive wildlife at zoological gardens in China.

BACKGROUND: Enterocytozoon bieneusi is a common opportunistic pathogen that is widely detected in humans, domestic animals and wildlife, and poses a challenge to public health. The present study was performed to evaluate the prevalence, genotypic diversity and zoonotic potential of E. bieneusi among wildlife at Chengdu and Bifengxia zoological gardens in Sichuan Province, China. RESULTS: Of the 272 fresh fecal samples harvested from 70 captive wildlife species at Chengdu Zoo (n = 198) and Bifengxia Zoo (n = 74), 21 (10.6 %) and 22 (29.7 %) tested positive for E. bieneusi by internal transcribed spacer (ITS) sequencing analysis, respectively. Specifically, genotypes D, Peru 6, CHB1, BEB6, CHS9, SC02 and SC03, and genotypes D, CHB1, SC01 and SC02 were detected in the Chengdu and Bifengxia Zoo samples, respectively. Five known genotypes (D, Peru 6, BEB6, CHS9 and CHB1) and three novel genotypes (SC01, SC02 and SC03) were clustered into the zoonotic group (group 1) and host-adapted group (group 2). Multilocus sequence typing (MLST) analysis targeting three microsatellites (MS1, MS3 and MS7) and one minisatellite (MS4) were successfully sequenced for 37, 33, 35 and 37 specimens, generating 8, 3, 11 and 15 distinct locus types, respectively. Altogether, we identified 27 multilocus genotypes (MLGs) among the E. bieneusi isolates by MLST. These data highlight the high genetic diversity of E. bieneusi among zoo wildlife. CONCLUSIONS: To our knowledge, this is the first report on the prevalence and genotypic diversity of E. bieneusi infections among captive wildlife in zoos in southwest China. Notably, we identified three novel E. bieneusi genotypes, as well as six new mammalian hosts (Asian golden cats, Tibetian blue bears, blackbucks, hog deer, Malayan sun bears and brown bears) for this organism. Moreover, the occurrence of zoonotic genotypes suggests that wildlife may act as reservoirs of E. bieneusi that can serve as a source of human microsporidiosis. The findings presented here should contribute to the control of zoonotic disease in China.

Adjuvant effects of recombinant giant panda (Ailuropoda melanoleuca) IL-18 on the canine distemper disease vaccine in mice.

Canine distemper virus (CDV) is a morbillivirus known to cause morbidity and mortality in a broad range of animals. Giant pandas (Ailuropoda melanoleuca), especially captive ones, are susceptible to natural infection with CDV. Interleukin-18 (IL-18) is a powerful adjuvant molecule that can enhance the development of antigen-specific immunity and vaccine efficacy. In this study, a giant panda IL-18 gene eukaryotic expression plasmid (pcAmIL-18) was constructed. Female BALB/c mice were muscularly inoculated with the plasmids pcAmIL-18, pcDNA3.1 and PBS, respectively. They were subsequently injected with an attenuated CDV vaccine for dogs, and the induced humoral and cellular responses were evaluated. The results showed that pcAmIL-18 remarkably improved the level of specific antibody, IFN-γ and IL-2 in mice sera, the T lymphocyte proliferation index and the percentage of CD4(+) and CD8(+) cells. These data indicated that pcAmIL-18 is a potential adjuvant that promotes specific immunity.

New subtype of Cryptosporidium cuniculus isolated from rabbits by sequencing the Gp60 gene.

Given the paucity of literature available on rabbits infected with Cryptosporidium in Sichuan Province (China), 290 fecal samples were collected from rabbits in the animal house of Sichuan Agricultural University, China and examined for Cryptosporidium oocysts using the Sheather's sucrose flotation technique and a modified acid-fast staining method. Three samples tested positive (prevalence = 1.03%). The positive isolates were genotyped by sequence analysis of the 18S rRNA, HSP70, COWP, and Cp135 genes and characterized by PCR-restriction fragment length polymorphism (RFLP) analysis of the 18S rRNA gene. Phylogenetic analysis was established using the neighbor-joining (NJ) method. All the isolates were identified as Cryptosporidium cuniculus. Further subtyping of the positive isolates was performed by DNA sequencing of the 60-kDa glycoprotein (gp60) gene. Only 1 subtype family was detected, Va, which was proposed to be a new subtype, VaA31. This study is the first report about the prevalence, genetic identification, and Cp135 gene of C. cuniculus in rabbits in Sichuan Province, China. The obtained results indicate that the C. cuniculus subtype in rabbits in Sichuan Province is unique.

A new genotype of Cryptosporidium from giant panda (Ailuropoda melanoleuca) in China.

Fifty-seven fecal samples were collected from giant pandas (Ailuropoda melanoleuca) in the China Conservation and Research Centre for the Giant Panda (CCRCGP) in Sichuan and examined for Cryptosporidium oocysts by Sheather's sugar flotation technique. An 18-year-old male giant panda was Cryptosporidium positive, with oocysts of an average size of 4.60×3.99 µm (n=50). The isolate was genetically analyzed using the partial 18S rRNA, 70 kDa heat shock protein (HSP70), Cryptosporidium oocyst wall protein (COWP) and actin genes. Multi-locus genetic characterization indicated that the present isolate was different from known Cryptosporidium species and genotypes. The closest relative was the Cryptosporidium bear genotype, with 11, 10, and 6 nucleotide differences in the 18S rRNA, HSP70, and actin genes, respectively. Significant differences were also observed in the COWP gene compared to Cryptosporidium mongoose genotype. The homology to the bear genotype at the 18S rRNA locus was 98.6%, which is comparable to that between Cryptosporidium parvum and Cryptosporidium hominis (99.2%), or between Cryptosporidium muris and Cryptosporidium andersoni (99.4%). Therefore, the Cryptosporidium in giant pandas in this study is considered as a new genotype: the Cryptosporidium giant panda genotype.

The mitochondrial genome of Baylisascaris procyonis.

BACKGROUND: Baylisascaris procyonis (Nematoda: Ascaridida), an intestinal nematode of raccoons, is emerging as an important helminthic zoonosis due to serious or fatal larval migrans in animals and humans. Despite its significant veterinary and public health impact, the epidemiology, molecular ecology and population genetics of this parasite remain largely unexplored. Mitochondrial (mt) genomes can provide a foundation for investigations in these areas and assist in the diagnosis and control of B. procyonis. In this study, the first complete mt genome sequence of B. procyonis was determined using a polymerase chain reaction (PCR)-based primer-walking strategy. METHODOLOGY/PRINCIPAL FINDINGS: The circular mt genome (14781 bp) of B. procyonis contained 12 protein-coding, 22 transfer RNA and 2 ribosomal RNA genes congruent with other chromadorean nematodes. Interestingly, the B. procyonis mtDNA featured an extremely long AT-rich region (1375 bp) and a high number of intergenic spacers (17), making it unique compared with other secernentean nematodes characterized to date. Additionally, the entire genome displayed notable levels of AT skew and GC skew. Based on pairwise comparisons and sliding window analysis of mt genes among the available 11 Ascaridida mtDNAs, new primer pairs were designed to amplify specific short fragments of the genes cytb (548 bp fragment) and rrnL (200 bp fragment) in the B. procyonis mtDNA, and tested as possible alternatives to existing mt molecular beacons for Ascaridida. Finally, phylogenetic analysis of mtDNAs provided novel estimates of the interrelationships of Baylisasaris and Ascaridida. CONCLUSIONS/SIGNIFICANCE: The complete mt genome sequence of B. procyonis sequenced here should contribute to molecular diagnostic methods, epidemiological investigations and ecological studies of B. procyonis and other related ascaridoids. The information will be important in refining the phylogenetic relationships within the order Ascaridida and enriching the resource of markers for systematic, population genetic and evolutionary biological studies of parasitic nematodes of socio-economic importance.