Metabarcoding-based diet analysis is a valuable tool for understanding the feeding behavior of a wide range of species. However, many studies using these methods for wild animals assume accuracy and precision without experimental evaluation with known positive control food items. Here, we conducted a feeding trial experiment with a positive control community in pasture-raised chickens and assessed the efficacy of several commonly used DNA extraction kits and primer sets. We hand-fed 22 known food items, including insects and plants, to six backyard laying hens and collected their excreta for eight h. We evaluated the efficacy of three DNA extraction kits, three primer sets for plant identification (targeting rbcL, trnL, and internal transcribed spacer 2 [ITS2]), and three primer sets for arthropod identification (targeting cytochrome oxidase subunit I [COI]). The detection success rate of our positive control food items was highly variable, ranging from 2.04% to 93.88% for all kit/primer combinations and averaging 37.35% and 43.57% for the most effective kit/primer combination for plants and insects, respectively. Extraction kits using bead-based homogenization positively affected the recovery proportion of plant and insect DNA in excreta samples. The minimum time to detect known food items was 44 min post-feeding. Two COI primer sets significantly outperformed the third, and both recovery proportion and taxonomic resolution from ITS2 were significantly higher than those from rbcL and trnL. Taken together, these results display the potential variability that can be inherently present in DNA-based diet analyses and highlight the utility of experimental feeding trials in validating such approaches, particularly for omnivores with diverse diets.

Heterochely denotes the presence of dissimilarly sized chelipeds on opposite sides of the body, a prevalent occurrence in diverse crustaceans. Conversely, heterometry pertains to the quantifiable disparities in size between these chelipeds. Both chelipeds hold pivotal roles in activities such as foraging, mating, and defense. Consequently, individuals of both genders in heterochelic species exhibit this morphological pattern. Previous studies have identified sexual dimorphism in cheliped size, with males displaying larger major chelipeds compared to females, albeit solely relying on propodus length as a size proxy and focusing solely on the major cheliped. In our study, we meticulously examined 190 specimens of the common European hermit crab Pagurus bernhardus from two collections. We sought to elucidate allometric relationships and assess whether heterometry exhibited sex-based differences when adjusting for body size by using ratios. Our findings revealed that male chelipeds displayed hyperallometric growth relative to females, and all three calculated heterometry indices exhibited significant disparities between the sexes. Consequently, male specimens exhibited larger major and minor chelipeds, even when theoretically matched for body size with females. This phenomenon may be attributed, among other factors, to male–male contests. Should indirect mate selection favor males with larger chelipeds in proportion to their body size, this dynamic could potentiate sexual selection in their favor.

The Pseudogastromyzon genus, consisting of species predominantly distributed throughout southeastern China, has garnered increasing market attention in recent years due to its ornamental appeal. However, the overlapping diagnostic attributes render the commonly accepted criteria for interspecific identification unreliable, leaving the phylogenetic relationships among Pseudogastromyzon species unexplored. In the present study, we undertake molecular phylogenetic and morphological examinations of the Pseudogastromyzon genus. Our phylogenetic analysis of mitochondrial genes distinctly segregated Pseudogastromyzon species into two clades: the Pseudogastromyzon clade and the Labigastromyzon clade. A subsequent morphological assessment revealed that the primary dermal ridge (specifically, the second ridge) within the labial adhesive apparatus serves as an effective and precise interspecific diagnostic characteristic. Moreover, the distributional ranges of Pseudogastromyzon and Labigastromyzon are markedly distinct, exhibiting only a narrow area of overlap. Considering the morphological heterogeneity of the labial adhesive apparatus and the substantial division within the molecular phylogeny, we advocate for the elevation of the Labigastromyzon subgenus to the status of a separate genus. Consequently, we have ascertained the validity of the Pseudogastromyzon and Labigastromyzon species, yielding a total of six valid species. To facilitate future research, we present comprehensive descriptions of the redefined species and introduce novel identification keys.

The genus Pseudamnicola Paulucci, 1878, is commonly found throughout the Mediterranean region. The genus displays considerable levels of endemism, accompanied by notable systematic and taxonomic ambiguity. However, the application of molecular data has proven highly effective in clarifying taxonomy and unveiling the diversity of cryptic species within the genus. Therefore, we employed all cytochrome c oxidase subunit I sequence data available and generated new ones from Greece to infer the phylogeny of the genus throughout its Mediterranean range and estimate the divergence times as well as the ancestral area of diversification. Our phylogenetic and time-estimate analyses demonstrate that with 36 to 38 extant Pseudamnicola species and genetic divergences across species ranging from 0.5% to 11.9% on average, the genus underwent relatively recent diversification during late Miocene (6.53 Ma), and the primary speciation events occurred during Plio-Pleistocene. The Italian Peninsula and Islands and the Ionian Drainages as defined by the Freshwater Ecoregions of the World are the ancestral regions of the genus following two different dispersal routes. Our study contributes to deepening our understanding of Pseudamnicola phylogeny by using data from throughout its range for the first time. This phylogeny provides evidence and confirms previous studies that relatively recent habitat isolation, followed by founder and dispersal events, has been one of the primary reasons for the evolution of the genus Pseudamnicola in the Mediterranean basin.

Global climate change is expected to have a profound effect on species distribution. Due to the temperature constraints, some narrow niche species could shift their narrow range to higher altitudes or latitudes. In this study, we explored the correlation between species traits, genetic structure, and geographical range size. More specifically, we analyzed how these variables are affected by differences in fundamental niche breadth or dispersal ability in the members of two sympatrically distributed stream-dwelling amphibian species (frog, Quasipaa yei; salamander, Pachyhynobius shangchengensis), in Dabie Mountains, East China. Both species showed relatively high genetic diversity in most geographical populations and similar genetic diversity patterns (JTX, low; BYM, high) correlation with habitat changes and population demography. Multiple clustering analyses were used to disclose differentiation among the geographical populations of these two amphibian species. Q. yei disclosed the relatively shallow genetic differentiation, while P. shangchengensis showed an opposite pattern. Under different historical climatic conditions, all ecological niche modeling disclosed a larger suitable habitat area for Q. yei than for P. shangchengensis; these results indicated a wider environment tolerance or wider niche width of Q. yei than P. shangchengensis. Our findings suggest that the synergistic effects of environmental niche variation and dispersal ability may help shape genetic structure across geographical topology, particularly for species with extremely narrow distribution.

Collection specimens provide valuable and often overlooked biological material that enables addressing relevant, long-unanswered questions in conservation biology, historical biogeography, and other research fields. Here, we use preserved specimens to analyze the historical distribution of the black francolin (Francolinus francolinus, Phasianidae), a case that has recently aroused the interest of archeozoologists and evolutionary biologists. The black francolin currently ranges from the Eastern Mediterranean and the Middle East to the Indian subcontinent, but, at least since the Middle Ages, it also had a circum-Mediterranean distribution. The species could have persisted in Greece and the Maghreb until the 19th century, even though this possibility had been questioned due to the absence of museum specimens and scant literary evidence. Nevertheless, we identified four 200-year-old stuffed black francolins—presumably the only ones still existing—from these areas and sequenced their mitochondrial DNA control region. Based on the comparison with conspecifics (n = 396) spanning the entirety of the historic and current species range, we found that the new samples pertain to previously identified genetic groups from either the Near East or the Indian subcontinent. While disproving the former occurrence of an allegedly native westernmost subspecies, these results point toward the role of the Crown of Aragon in the circum-Mediterranean expansion of the black francolin, including the Maghreb and Greece. Genetic evidence hints at the long-distance transport of these birds along the Silk Road, probably to be traded in the commerce centers of the Eastern Mediterranean.

The marine environment presents challenges for wound healing in cetaceans, despite their remarkable recovery abilities with minimal infections or complications. However, the molecular mechanism underlying this efficient wound healing remains underexplored. To better understand the molecular mechanisms behind wound healing in cetaceans, we investigated the evolutionary patterns of 37 wound healing-related genes in representative mammals. We found wound healing-related genes experience adaptive evolution in cetaceans: (1) Three extrinsic coagulation pathway-related genes—tissue factor (F3), coagulation factor VII (F7), and coagulation factor X (F10)—are subject to positive selection in cetaceans, which might promote efficient hemostasis after injury; positive selection in transforming growth factor-beta 2 (TGF-β2), transforming growth factor-beta 3 (TGF-β3), and platelet-derived growth factor D (PDGFD), which play immunological roles in wound healing, may help cetaceans enhance inflammatory response and tissue debridement. (2) Coagulation factor XII (F12) is the initiation factor in the intrinsic coagulation pathway. It had a premature stop codon mutation and was subjected to selective stress relaxation in cetaceans, suggesting that the early termination of F12 may help cetaceans avoid the risk of vascular blockage during diving. (3) Fibrinogen alpha chain (FGA) and FIII, which were detected to contain the specific amino acid substitutions in marine mammals, indicating similar evolutionary mechanisms might exist among marine mammals to maintain strong wound-healing ability. Thus, our research provides further impetus to study the evolution of the wound healing system in cetaceans and other marine mammals, extending knowledge of preventing coagulation disorder and atherosclerosis in humans.

The evolutionary dynamics of the ecoregions of southern South America and the species that inhabit them have been poorly studied, and few biogeographic hypotheses have been proposed and tested. Quaternary climatic oscillations are among the most important processes that have led to the current distribution of genetic variation in different regions of the world. In this work, we studied the evolutionary history and distribution of the Córdoba vesper mouse (Calomys venustus), a characteristic rodent of the region of which little is known about its natural history. Since the population dynamics of this species are influenced by climatic factors, this rodent is a suitable model to study the effects of Quaternary climatic oscillations in central Argentina. The mitochondrial cytochrome b gene was sequenced to analyze the phylogeography of C. venustus, and ecological niche modeling tools were used to map its potential distributions. The results of these approaches were combined to provide additional spatially explicit information about this species’ past. Our results suggest that the Espinal was the area of origin of this species, which expanded demographically and spatially during the last glacial period. A close relationship was found between the Espinal and the Mountain Chaco. These results are consistent with previous studies and emphasize the role of the Espinal in the biogeographic history of southern South America as an area of origin of several species.

Blister beetles (Coleoptera: Meloidae) produce a natural defensive toxin cantharidin (CTD), which has been used for various cancer treatments and other diseases. Currently, the lack of chromosome-level reference genomes in Meloidae limits further understanding of the mechanism of CTD biosynthesis and environmental adaptation. In this study, the chromosome-level genome assembly of Mylabris phalerata was generated based on PacBio and Hi-C sequencing. This reference genome was about 136.68 Mb in size with contig N50 of 9.17 Mb and composed of 12 chromosomes. In comparison to six other Coleoptera insects, M. phalerata exhibited multiple expanded gene families enriched in juvenile hormone (JH) biosynthetic process pathway, farnesol dehydrogenase activity, and cytochrome P450, which may be related to CTD biosynthesis. Consistently, the transcriptomic analysis suggested the “terpenoid backbone biosynthesis” pathway and “the juvenile hormone” as putative core pathways of CTD biosynthesis and presented eight up-regulated differential expression genes in male adults as candidate genes. It is possible that the restricted feeding niche and lifestyle of M. phalerata were the cause of the gene family’s contraction of odorant binding proteins. The ABC transporters (ABCs) related to exporting bound toxins out of the cell and the resistance to the self-secreted toxins (e.g. CTD) were also contracted, possibly due to other self-protection strategies in M. phalerata. A foundation of understanding CTD biosynthesis and environmental adaptation of blister beetles will be established by our reference genome and discoveries.

Using a complex analysis of the molecular genetics, morphological, and ecological characteristics of Hippoboscidae flies, the phylogenetic structure and trends in the evolution of morphological characters that contribute to the ectoparasitic lifestyle of hippoboscid flies of the north of Eurasia were studied for the first time. The research was carried out on 26 Palearctic species from 10 genera. The analysis of molecular phylogeny revealed the levels of clustering of the family with the species predominantly parasitizing mammals or birds, the time of cluster formation, and the divergence of species in the Palearctic conditions. An independent adaptation to birds occurred in the genera Icosta, Pseudolynchia, Ornithoica, and others. Bird parasites are characterized by bifid tarsal claws, long hooks on pulvilli, and long empodium setae (except genus Ornithoica). Mammalian parasites are characterized by simple tarsal claws, short lobes of hooks on pulvilli, and zones on empodium with short setae. Specialization in empodium and pulvillus morphotypes and wing reduction are higher diverged in mammalian parasites than in bird parasites. The decrease of flight ability and wing reduction independently arose in different subfamilies of Hippoboscidae flies. Our results assume that the tribe Ornithomyini is a paraphyletic group, since, according to the complex of morphological features, the genus Ornithoica can be considered a separate lineage of evolution.

Musk secreted by male forest musk deer (Moschus berezovskii) musk glands is an invaluable component of medicine and perfume. Musk secretion depends on musk gland maturation; however, the mechanism of its development remains elusive. Herein, using single cell multiome ATAC + gene expression coupled with several bioinformatic analyses, a dynamic transcriptional cell atlas of musk gland development was revealed, and key genes and transcription factors affecting its development were determined. Twelve cell types, including two different types of acinar cells (Clusters 0 and 10) were identified. Single-nucleus RNA and single-nucleus ATAC sequencing analyses revealed that seven core target genes associated with musk secretion (Hsd17b2, Acacb, Lss, Vapa, Aldh16a1, Aldh7a1, and Sqle) were regulated by 12 core transcription factors (FOXO1, CUX2, RORA, RUNX1, KLF6, MGA, NFIC, FOXO3, ETV5, NR3C1, HSF4, and MITF) during the development of Cluster 0 acinar cells. Kyoto Encyclopedia of Genes and Genomes enrichment showed significant changes in the pathways associated with musk secretion during acinar cell development. Gene set variation analysis also revealed that certain pathways associated with musk secretion were enriched in 6-year-old acinar cells. A gene co-expression network was constructed during acinar cell development to provide a precise understanding of the connections between transcription factors, genes, and pathways. Finally, intercellular communication analysis showed that intercellular communication is involved in musk gland development. This study provides crucial insights into the changes and key factors underlying musk gland development, which serve as valuable resources for studying musk secretion mechanisms and promoting the protection of this endangered species.

Blister beetles (Coleoptera: Meloidae) are currently subdivided into three subfamilies: Eleticinae (a basal group), Nemognathinae, and Meloinae. These are all characterized by the endogenous production of the defensive terpene cantharidin (CA), whereas the two most derived subfamilies show a hypermetamorphic larval development. Here, we provide novel draft genome assemblies of five species sampled across the three blister beetle subfamilies (Iselma pallidipennis, Stenodera caucasica, Zonitis immaculata, Lydus trimaculatus, and Mylabris variabilis) and performed a comparative analysis with other available Meloidae genomes and the closely-related canthariphilous species (Pyrochroa serraticornis) to disclose adaptations at a molecular level. Our results highlighted the expansion and selection of genes potentially responsible for CA production and metabolism, as well as its mobilization and vesicular compartmentalization. Furthermore, we observed adaptive selection patterns and gain of genes devoted to epigenetic regulation, development, and morphogenesis, possibly related to hypermetamorphosis. We hypothesize that most genetic adaptations occurred to support both CA biosynthesis and hypermetamorphosis, two crucial aspects of Meloidae biology that likely contributed to their evolutionary success.

Bats, members of the Chiroptera order, rank as the second most diverse group among mammals. Recent molecular systematic studies on bats have successfully classified 21 families within two suborders: Yinpterochiroptera and Yangochiroptera. Nevertheless, the phylogeny within these 21 families has remained a subject of controversy. In this study, we have employed a balanced approach to establish a robust family-level phylogenetic hypothesis for bats, utilizing a more comprehensive molecular dataset. This dataset includes representative species from all 21 bat families, resulting in a reduced level of missing genetic information. The resulting phylogenetic tree comprises 21 lineages that are strongly supported, each corresponding to one of the bat families. Our findings support to place the Emballonuroidea superfamily as the basal lineage of Yangochiroptera, and that Myzopodidae should be situated as a basal lineage of Emballonuroidea, forming a sister relationship with the clade consisting of Nycteridae and Emballonuridae. Finally, we have conducted dating analyses on this newly resolved phylogenetic tree, providing divergence times for each bat family. Collectively, our study has employed a relatively comprehensive molecular dataset to establish a more robust phylogeny encompassing all 21 bat families. This improved phylogenetic framework will significantly contribute to our understanding of evolutionary processes, ecological roles, disease dynamics, and biodiversity conservation in the realm of bats.