Non-lethal scute transcriptomics reveals biomarker candidates for environmental stress in Mexican Morelet's...

Aquatic ecosystems in southeastern Mexico face increasing chemical pollution, yet non-lethal molecular tools to monitor native reptile populations remain limited. We performed RNA sequencing (RNA-seq) on non-lethal caudal scute biopsies of 43 Morelet's crocodiles (Crocodylus moreletii) across six wild populations and a captive reference group. After mapping to the Alligator mississippiensis genome, 18,107 expressed genes were analyzed using multivariate, differential-expression, co-expression network, and single-sample pathway scoring approaches. Wild crocodiles showed a consistent shift along a dominant transcriptional axis contrasting dermal/extracellular matrix remodeling with epidermal-metabolic/barrier programs. A conserved core of differentially expressed genes (DEGs) and wild-associated modules were enriched for extracellular matrix organization, secretory stress handling, vesicle trafficking and adhesion pathways, while union DEGs captured site-specific alterations in lipid metabolism, ion transport, neuronal signaling and developmental processes. Integrating network connectivity, differential expression strength, statistical support, and receiver operating characteristic performance, we prioritized 10 up-regulated genes as a preliminary scute biomarker panel that discriminates wild from captive individuals. This landscape-scale scute transcriptomic framework supports non-lethal omics-based biomonitoring in crocodilians and other non-model wildlife.