Genetic diversity and population differentiation of the Chinese soft-shelled turtle (Pelodiscus sinensis) in three geographical populations

Single-cell transcriptomes from turtle livers reveal sensitivity of hepatic immune cells to bacteria-infection

The liver is important in the synthesis, metabolism and storage of nutrients, detoxification and immune response of the body, and the liver immune response against exogenous pathogens from the intestinal tract plays a key role in the immune activities. However, the cellular composition of the liver immune atlas remains sparsely studied in reptiles. We used single-cell RNA sequencing to identify the cellular profile of the liver of the Chinese soft-shelled turtle (Pelodiscus sinensis). We obtained the transcriptional landscape based on 9938 cells from the fractionation of fresh hepatic tissues from two individuals, uninfected and infected with bacteria (Aeromonas hydrophila). We identified seven hepatic immune cell subsets, including plasma, erythroid, T/NK, B, endothelial, dendritic and Kupffer cells. Bacteria-infection altered the number of liver immune cells, as revealed by the fact that the infected turtle had more plasma, endothelial and Kupffer cells and fewer T/NK, dendritic and erythroid cells than did the uninfected turtle. Our study is the first to provide a comprehensive view of the hepatic immune landscape of P. sinensis at the single-cell resolution that outlines the characteristics of immune cells in the turtle liver and provides a liver transcriptome baseline for turtle immunology.

Population genetics of wild and captive Trachemysvenusta (Gray, 1856) (Reptilia: Emydidae) in the Usumacinta river basin in Mexico

The Mesoamerican slider Trachemysvenusta is endemic to Central America and Southern Mexico. Several human-mediated disturbances, including habitat degradation and illegal hunting for food, have impacted its populations along the Usumacinta river basin. The extent to which these disturbances have affected the genetic diversity and population structure of T. venusta inhabiting the basin remains unresolved. To this end, we analyzed eight microsatellite markers in five wild populations of T. venusta from the middle and lower reaches of the basin as well as one captive population. Our results show high levels of genetic diversity for all analyzed populations, low F _ST values, high gene flow and no genetic structure, indicating an absence of genetic differentiation across sites and, thus, a single panmictic population for the basin. Evidence of a genetic bottleneck was observed in two of the wild populations (and the captive one), indicating some impact from disturbances, whether from poaching or habitat fragmentation, despite the seemingly high connectivity of most populations. Results are discussed in terms of the relative importance of genetic parameters for the conservation of T. venusta, particularly in light of the importance of demographic stochasticity in local conditions undergoing rapid changes.

Millennium-old farm breeding of Chinese softshell turtles (Pelodiscus spp.) results in massive erosion of biodiversity

Chinese softshell turtles (Pelodiscus spp.) are widely distributed, ranging from the Amur and Ussuri Rivers in the Russian Far East through the Korean Peninsula, Japan, and eastern, central, and southern China to southern Vietnam. In East and Southeast Asia, Chinese softshell turtles are traditionally exploited for food and have been farm-bred in China since the Spring and Autumn Period, more than 2400 years ago. Currently, the annual production of Pelodiscus amounts to 340,000 t in China alone. Using mitochondrial DNA (2428 bp) and five nuclear loci (3704 bp), we examined broad sampling of wild and farm-bred Pelodiscus to infer genetic and taxonomic differentiation. We discovered four previously unknown mitochondrial lineages, all from China. One lineage from Jiangxi is deeply divergent and sister to the mitochondrial lineage of Pelodiscus axenaria. The nuclear loci supported species status for P. axenaria and the new lineage from Jiangxi. Pelodiscus maackii and P. parviformis, both harboring distinct mitochondrial lineages, were not differentiated from P. sinensis in the studied nuclear markers. The same is true for two new mitochondrial lineages from Zhejiang, China, represented by only one individual each, and another new lineage from Anhui, Guangdong, Jiangxi and Zhejiang, China. However, Vietnamese turtles yielding a mitochondrial lineage clustering within P. sinensis were distinct in nuclear markers, suggesting that these populations could represent another unknown species with introgressed mitochondria. Its species status is also supported by the syntopic occurrence with P. sinensis in northern Vietnam and by morphology. In addition, we confirmed sympatry of P. axenaria and P. parviformis in Guangxi, China, and found evidence for sympatry of P. sinensis and the new putative species from Jiangxi, China. We also discovered evidence for hybridization in turtle farms and for the occurrence of alien lineages in the wild (Zhejiang, China), highlighting the risk of genetic pollution of native stock. In the face of the large-scale breeding of Pelodiscus, we claim that the long-term survival of distinct genetic lineages and species can only be assured when an upscale market segment for pure-bred softshell turtles is established, making the breeding of pure lineages lucrative for turtle farms. Our findings underline that the diversity of Pelodiscus is currently underestimated and threatened by anthropogenic admixture. We recommend mass screening of genetic and morphological variation of Chinese softshell turtles as a first step to understand and preserve their diversity.

Yangtze River, an insignificant genetic boundary in tufted deer (Elaphodus cephalophus): the evidence from a first population genetics study

Great rivers were generally looked at as the geographical barrier to gene flow for many taxonomic groups. The Yangtze River is the third largest river in the world, and flows across South China and into the East China Sea. Up until now, few studies have been carried out to evaluate its effect as a geographical barrier. In this study, we attempted to determine the barrier effect of the Yangtze River on the tufted deer (Elaphodus cephalophus) using the molecular ecology approach. Using mitochondrial DNA control region (CR) sequences and 13 nuclear microsatellite loci, we explored the genetic structure and gene flow in two adjacent tufted deer populations (Dabashan and Wulingshan populations), which are separated by the Yangtze River. Results indicated that there are high genetic diversity levels in the two populations, but no distinguishable haplotype group or potential genetic cluster was detected which corresponded to specific geographical population. At the same time, high gene flow was observed between Wulingshan and Dabashan populations. The tufted deer populations experienced population decrease from 0.3 to 0.09 Ma BP, then followed by a distinct population increase. A strong signal of recent population decline (T = 4,396 years) was detected in the Wulingshan population by a Markov-Switching Vector Autoregressions(MSVAR) process population demography analysis. The results indicated that the Yangtze River may not act as an effective barrier to gene flow in the tufted deer. Finally, we surmised that the population demography of the tufted deer was likely affected by Pleistocene climate fluctuations and ancient human activities.

Genetic variation in Whitmania pigra, Hirudo nipponica and Poecilobdella manillensis, three endemic and endangered species in China using SSR and TRAP markers

Leeches are not only important medicinal animals worldwide but also are endangered. We aimed to (i) explore the level of genetic diversity within/among populations of three leeches, (ii) assess genetic differentiation among these three leeches, and (iii) discuss an appropriate strategy for conserving leech germplasm. A total of 315 individuals of Whitmania pigra, Hirudo nipponica and Poecilobdella manillensis from 21 populations were collected in China and Vietnam. The genetic structure and genetic diversity among and within the 21 populations were evaluated using target region amplified polymorphism (TRAP) and simple sequence repeat (SSR) markers. Sixteen pairs of TRAP primers generated a total of 398 fragments, of which 396 (99.50%) were polymorphic; fourteen pairs of SSR primers generated a total of 60 fragments, of which 59 (98.33%) were polymorphic. Shannon's index (I) and Nei's gene diversity index (H) for the three leeches were high at the species level (I=0.4980 and H=0.3323 for TRAPs, I=0.4487 and H=0.2969 for SSRs in W. pigra; I=0.4147/0.3769, H=0.2788/0.2566 for H. nipponica; and I=0.4616/0.4717, H=0.3099/0.3203 for P. manillensis). However, low genetic diversity was determined at the population level; the average genetic diversity measures within populations were H=0.1767/0.1376, I=0.2589/0.2043 for W. pigra, H=0.2149/0.2021, I=0.3184/0.3000 for H. nipponica and H=0.2850/0.2724, I=0.4152/0.3967 for P. manillensis. We conclude that there was limited gene exchange within/among populations and species, as the gene flow number (Nm) was 0.5493/0.5807. However, for all three species, the genetic diversity was different at the population level. Gene differentiation (Gst) and Nm were 0.4682 /0.5364 and 0.5678/0.4321 for W. pigra, 0.2294/0.2127 and 1.6797/1.8512 for H. nipponica and 0.1214/0.1496 and 3.6202/2.8412 for P. manillensis. STRUCTURE analysis, Unweighted Pair-Group Method with Arithmetic means (UPGMA) cluster analysis and Principal Coordinates Analysis (PCOA) all yielded similar results. The isolation-by-distance pattern was not significant for any of the three species by the Mantel test. These data emphasize the need for management, conservation, and rehabilitation of this animal species. Finally, an appropriate strategy for conserving leech is proposed.