Global climate change drives variations in species distribution patterns and affects biodiversity, potentially increasing the risk of species extinction. Investigating the potential distribution range of species under future global climate change is crucial for biodiversity conservation and ecosystem management. In this study, we collected distributional data for 5282 reptile species to assess their conservation status based on distributional ranges using species distribution models. Our predictions indicate that the potential distribution ranges for over half of these species are projected to decrease under different scenarios. Under future scenarios with relatively low carbon emissions, the increase in the number of threatened reptiles is significantly lower, highlighting the importance of human efforts. Surprisingly, we identified some endangered species that are projected to expand their distribution ranges, underscoring the potential positive effects of climate change on some special species. Our findings emphasize the increased extinction risk faced by reptile species due to climate change and highlight the urgent need to mitigate the effects of habitat degradation and human activities on their potential distribution in the future.

Biological invasions are among the threats to global biodiversity and social sustainability, especially on islands. Identifying the threshold of area at which non-native species begin to increase abruptly is crucial for early prevention strategies. The small-island effect (SIE) was proposed to quantify the nonlinear relationship between native species richness and area but has not yet been applied to non-native species and thus to predict the key breakpoints at which established non-native species start to increase rapidly. Based on an extensive global dataset, including 769 species of non-native birds, mammals, amphibians and reptiles established on 4277 islands across 54 archipelagos, we detected a high prevalence of SIEs across 66.7% of archipelagos. Approximately 50% of islands have reached the threshold area and thus may be undergoing a rapid increase in biological invasions. SIEs were more likely to occur in those archipelagos with more non-native species introduction events, more established historical non-native species, lower habitat diversity and larger archipelago area range. Our findings may have important implications not only for targeted surveillance of biological invasions on global islands but also for predicting the responses of both non-native and native species to ongoing habitat fragmentation under sustained land-use modification and climate change.

Epoxiconazole (EPX) is a world widely used chiral triazole fungicide in the agriculture field. The excessive application of this triazole may cause damage to lizards. However, limited information is known about the toxicokinetics of EPX on lizards. Our study aimed to investigate the enantioselective absorption, distribution, metabolism, and elimination (ADME) of EPX in lizards following low and high dose exposure (10 and 100 mg kg-1 bodyweitht (bw)). The results demonstrated that (+)-EPX was easier absorbed than (-)-EPX in lizard plasma. Both (+)-EPX and (-)-EPX were detected in the liver, gonad, kidney, skin, brain, and intestine, with (+)-EPX preferentially distributed in these tissues. The elimination of (-)-EPX was faster than that of (+)-EPX in lizard liver and kidney in the high dose groups. Chiral conversion was found between EPX enantiomers in lizard skin. Simultaneously, five metabolites including M2, M4, M10, M18 and M19 were detected in lizard liver and kidney after EPX enantiomers exposure. The relative concentrations of M2, M4, and M10 were higher in the liver and kidney of (-)-EPX groups than those produced from (+)-EPX groups. The metabolic enzymes CYP3A4 and SULT1A1 primarily mediated enantioselective metabolism of EPX. The conclusions drawn from this study significantly enhance our understanding of the enantioselective behaviors of chiral triazole fungicides in reptiles, offering essential guidance for assessing the risks associated with different enantiomers of triazole fungicides.

Fossil feathers have transformed our understanding of integumentary evolution in vertebrates. The evolution of feathers is associated with novel skin ultrastructures, but the fossil record of these changes is poor and thus the critical transition from scaled to feathered skin is poorly understood. Here we shed light on this issue using preserved skin in the non-avian feathered dinosaur Psittacosaurus. Skin in the non-feathered, scaled torso is three-dimensionally replicated in silica and preserves epidermal layers, corneocytes and melanosomes. The morphology of the preserved stratum corneum is consistent with an original composition rich in corneous beta proteins, rather than (alpha-) keratins as in the feathered skin of birds. The stratum corneum is relatively thin in the ventral torso compared to extant quadrupedal reptiles, reflecting a reduced demand for mechanical protection in an elevated bipedal stance. The distribution of the melanosomes in the fossil skin is consistent with melanin-based colouration in extant crocodilians. Collectively, the fossil evidence supports partitioning of skin development in Psittacosaurus: a reptile-type condition in non-feathered regions and an avian-like condition in feathered regions. Retention of reptile-type skin in non-feathered regions would have ensured essential skin functions during the early, experimental stages of feather evolution.

The Stat (signal transducer and activator of transcription) gene family plays a vital role in regulating immunity and the processes of cellular proliferation, differentiation, and apoptosis across diverse organisms. Although the functions of Stat genes in immunity have been extensively documented in many mammals, limited data are available for reptiles. We used phylogenetic analysis to identify eight putative members of the Stat family (Stat1-1, Stat1-2, Stat2, Stat3, Stat4, Stat5b, Stat6-1, and Stat6-2) within the genome of M. reevesii, a freshwater turtle found in East Asia. Sequence analysis showed that the Stat genes contain four conserved structural domains protein interaction domain, coiled-coil domain, DNA-binding domain, and Src homology domain 2. In addition, Stat1, Stat2, and Stat6 contain TAZ2bind, Apolipo_F, and TALPID3 structural domains. The mRNA levels of Stat genes were upregulated in spleen tissues at 4, 8, 12, and 16 h after administration of lipopolysaccharide, a potent activator of the immune system. Stat5b expression at 12-h LPS post-injection exhibited the most substantial difference from the control. The expression of Stat5b in spleen tissue cellular was verified by immunofluorescence. These results suggest that Stat5b plays a role in the immune response of M. reevesii and may prove to be as a positive marker of an immune response in future studies.

Reptiles are an important, yet often understudied, taxon in nature conservation. They play a significant role in ecosystems1 and can serve as indicators of environmental health, often responding more rapidly to human pressures than other vertebrate groups.2 At least 21% of reptiles are currently assessed as threatened with extinction by the IUCN.3 However, due to the lack of comprehensive global assessments until recently, they have been omitted from spatial studies addressing conservation or spatial prioritization (e.g., Rosauer et al.,4,5,6,7,8 Fritz and Rahbek,4,5,6,7,8 Farooq et al.,4,5,6,7,8 Meyer et al., 4,5,6,7,8 and Farooq et al.4,5,6,7,8). One important knowledge gap in conservation is the lack of spatially explicit information on the main threats to biodiversity,9 which significantly hampers our ability to respond effectively to the current biodiversity crisis.10,11 In this study, we calculate the probability of a reptile species in a specific location being affected by one of seven biodiversity threats-agriculture, climate change, hunting, invasive species, logging, pollution, and urbanization. We conducted the analysis at a global scale, using a 50 km × 50 km grid, and evaluated the impact of these threats by studying their relationship with the risk of extinction. We find that climate change, logging, pollution, and invasive species are most linked to extinction risk. However, we also show that there is considerable geographical variation in these results. Our study highlights the importance of going beyond measuring the intensity of threats to measuring the impact of these separately for various biogeographical regions of the world, with different historical contingencies, as opposed to a single global analysis treating all regions the same.

The zoonotic potential of the protist parasites Cryptosporidium spp. and Giardia duodenalis in amphibians and reptiles raises public health concerns due to their growing popularity as pets. This review examines the prevalence and diversity of these parasites in wild and captive amphibians and reptiles to better understand the zoonotic risk. Research on Giardia in both groups is limited, and zoonotic forms of Cryptosporidium or Giardia have not been reported in amphibians. Host-adapted Cryptosporidium species dominate in reptiles, albeit some reptiles have been found to carry zoonotic (C. hominis and C. parvum) and rodent-associated (C. tyzzeri, C. muris and C. andersoni) species, primarily through mechanical carriage. Similarly, the limited reports of Giardia duodenalis (assemblages A, B and E) in reptiles may also be due to mechanical carriage. Thus, the available evidence indicates minimal zoonotic risk associated with these organisms in wild and captive frogs and reptiles. The exact transmission routes for these infections within reptile populations remain poorly understood, particularly regarding the importance of mechanical carriage. Although the risk appears minimal, continued research and surveillance efforts are necessary to gain a more comprehensive understanding of the transmission dynamics and ultimately improve our ability to safeguard human and animal health.

BACKGROUND: The family Rhabdiasidae (Nematoda: Rhabditida) is a globally distributed group of nematode parasites, with over 110 species parasitic mainly in amphibians and reptiles. However, the systematic position of the family Rhabdiasidae in the order Rhabditida remains unsolved, and the evolutionary relationships among its genera are still unclear. Moreover, the present knowledge of the mitochondrial genomes of rhabdiasids remains limited.
METHODS: Two rhabdiasid species: Rhabdias kafunata Sata, Takeuchi & Nakano, 2020 and R. bufonis (Schrank, 1788) collected from the Asiatic toad Bufo gargarizans Cantor (Amphibia: Anura) in China, were identified based on morphology (light and scanning electron microscopy) and molecular characterization (sequencing of the nuclear 28S and ITS regions and mitochondrial cox1 and 12S genes). The complete mitochondrial genomes of R. kafunata and R. bufonis were also sequenced and annotated for the first time. Moreover, phylogenetic analyses based on the amino acid sequences of 12 protein-coding genes (PCGs) of the mitochondrial genomes were performed to clarify the systematic position of the family Rhabdiasidae in the order Rhabditida using maximum likelihood (ML) and Bayesian inference (BI). The phylogenetic analyses based on the 28S + ITS sequences, were also inferred to assess the evolutionary relationships among the genera within Rhabdiasidae.
RESULTS: The detailed morphology of the cephalic structures, vulva and eggs in R. kafunata and R. bufonis was revealed using scanning electron microscopy (SEM) for the first time. The characterization of 28S and ITS regions of R. kafunata was reported for the first time. The mitogenomes of R. kafunata and R. bufonis are 15,437 bp and 15,128 bp long, respectively, and both contain 36 genes, including 12 PCGs (missing atp8). Comparative mitogenomics revealed that the gene arrangement of R. kafunata and R. bufonis is different from all of the currently available mitogenomes of nematodes. Phylogenetic analyses based on the ITS + 28S data showed Neoentomelas and Kurilonema as sister lineages, and supported the monophyly of Entomelas, Pneumonema, Serpentirhabdias and Rhabdias. Mitochondrial phylogenomic results supported Rhabdiasidae as a member of the superfamily Rhabditoidea in the suborder Rhabditina, and its occurrance as sister to the family Rhabditidae.
CONCLUSIONS: The complete mitochondrial genome of R. kafunata and R. bufonis were reported for the first time, and two new gene arrangements of mitogenomes in Nematoda were revealed. Mitogenomic phylogenetic results indicated that the family Rhabdiasidae is a member of Rhabditoidea in Rhabditina, and is closely related to Rhabditidae. Molecular phylogenies based on the ITS + 28S sequence data supported the validity of Kurilonema, and showed that Kurilonema is sister to Neoentomelas. The present phylogenetic results also indicated that the ancestors of rhabdiasids seem to have initially infected reptiles, then spreading to amphibians.

Trait datasets are increasingly being used in studies investigating eco-evolutionary theory and global conservation initiatives. Reptiles are emerging as a key group for studying these questions because their traits are crucial for understanding the ability of animals to cope with environmental changes and their contributions to ecosystem processes. We collected data from earlier databases, and the primary literature to create an up-to-date dataset of reptilian traits, encompassing 40 traits from 12060 species of reptiles (Archelosauria: Crocodylia and Testudines, Rhynchocephalia, and Squamata: Amphisbaenia, Sauria, and Serpentes). The data were gathered from 1288 sources published between 1820 and 2023. The dataset includes morphological, physiological, behavioral, and life history traits, as well as information on the availability of genetic data, IUCN Red List assessments, and population trends.

Male and female Chinese soft-shelled turtles (Pelodiscus sinensis) have sex-dimorphic growth patterns, and males have higher commercial value because of their larger size and thicker calipash. Thus, developing sex-specific markers is beneficial to studies on all-male breeding in P. sinensis. Here, we developed an accurate and efficient workflow for the screening of sex-specific sequences with ZW or XY sex determination systems. Based on this workflow, female and male P. sinensis reference genomes of 2.23 Gb and 2.26 Gb were obtained using de novo assembly. After aligning and filtering, 4.01 Mb female-specific sequences were finally identified. Subsequently, the seven developed sex-specific primer pairs were 100% accurate in preliminary, population, and embryonic validation. The presence and absence of bands for the primers of P44, P45, P66, P67, P68, and P69, as well as two and one bands for the PB1 primer, indicate that the embryos are genetically female and male, respectively. NR and functional annotations identified several sex-determining candidate genes and related pathways, including Ran, Eif4et, and Crkl genes, and the insulin signaling pathway and the cAMP signaling pathway, respectively. Collectively, our results reveal that a ZW-type sex-determination system is present in P. sinensis and provide novel insights for the screening of sex-specific markers, sex-control breeding, and the studies of the sex determination mechanism of P. sinensis.