Size Evolution Research Articles

A theoretical and experimental investigation of continuous oil–water gravity separation

In this study, we have developed a mathematical model for a three-phase separator. The model consists of two sections: the inlet section and the separation section, separated by a perforated calming baffle. In the inlet section, two dispersion layers undergo droplet size evolution due to turbulent breakage and coalescence, described by a spatially homogeneous PBE. In the separation section, the two dispersion layers flow alongside each other and interact at an interface. The volumetric flow and velocity profiles are influenced by interfacial coalescence, with considerations for plug and laminar flow assumptions. The model incorporates droplet gravity-driven transport using the Kumar and Hartland model, binary and interfacial coalescence employing a film drainage model, and an effective diffusion term to account for the formation of the dense packed layer which ensures a physical volume fraction range of 0–1. Steady-state and transient numerical solvers are developed to solve the resulting advection–diffusion equations. Additionally, a series of experiments were conducted using a lab-scale multi-parallel pipes separator to investigate the impact of varying volume fractions and flow rates on the separation efficiency of the equipment. The model results are compared with the experimental data which shows relatively good agreement.

Viviparity imparts a macroevolutionary signature of ecological opportunity in the body size of female Liolaemus lizards

Viviparity evolved ~115 times across squamate reptiles, facilitating the colonization of cold habitats, where oviparous species are scarce or absent. Whether the ecological opportunity furnished by such colonization reconfigures phenotypic diversity and accelerates evolution is unclear. We investigated the association between viviparity and patterns and rates of body size evolution in female Liolaemus lizards, the most species-rich tetrapod genus from temperate regions. Here, we discover that viviparous species evolve ~20% larger optimal body sizes than their oviparous relatives, but exhibit similar rates of body size evolution. Through a causal modeling approach, we find that viviparity indirectly influences body size evolution through shifts in thermal environment. Accordingly, the colonization of cold habitats favors larger body sizes in viviparous species, reconfiguring body size diversity in Liolaemus. The catalyzing influence of viviparity on phenotypic evolution arises because it unlocks access to otherwise inaccessible sources of ecological opportunity, an outcome potentially repeated across the tree of life.

In vitro gastrointestinal digestion of cow’s and sheep’s dairy products: Impact of species and structure

Sheep’s milk (SM) is known to differ from cow’s milk (CM) in nutritional composition and physicochemical properties, which may lead to different digestion behaviours. This work aimed to investigate the impact of the species (cow vs sheep) and the structure (milk vs yogurt) on the digestion of dairy products. Using an in vitro static gastrointestinal digestion model, CM, SM, cow’s milk yogurt (CY) and sheep’s milk yogurt (SY) were compared on particle size evolution, microscopic observations, degree of lipolysis, degree of proteolysis, specific protein degradation and calcium bioaccessibility. Species and structure affected particle size evolution during the gastric phase resulting in smaller particles for yogurts compared to milks as well as for CM products compared to SM products. Species impacted lipid composition and lipolysis, with SM products presenting higher short/medium-chain fatty acids content and higher intestinal degree of lipolysis. Proteolysis was influenced by structure, with milks showing higher intestinal degree of proteolysis compared to yogurts. Caseins were digested faster in CM, ⍺-lactalbumin was digested faster in SM despite its higher concentration, and during gastric digestion β-lactoglobulin was more degraded in CM products compared to SM products and more in yogurts compared to milks. Lastly, SM products released more bioaccessible calcium than CM products. In conclusion, species (cow vs sheep) impacted more the digestion compared to the structure (milk vs yogurt). In fact, SM was different from CM mainly due to a denser protein network that might slow down the accessibility of the enzyme to its substrate which induce a delay of gastric disaggregation and thus lead to slower the digestion of the nutrients.

Balancing Polarization and Breakdown for High Capacitive Energy Storage by Microstructure Design.

The compromise of contradictive parameters, polarization, and breakdown strength, is necessary to achieve a high energy storage performance. The two can be tuned, regardless of material types, by controlling microstructures: amorphous states possess higher breakdown strength, while crystalline states have larger polarization. However, how to achieve a balance of amorphous and crystalline phases requires systematic and quantitative investigations. Herein, the trade-off between polarization and breakdown field is comprehensively evaluated with the evolution of microstructure, i.e., grain size and crystallinity, by phase-field simulations. The results indicate small grain size (≈10-35nm) with moderate crystallinity (≈60-80%) is more beneficial to maintain relatively high polarization and breakdown field simultaneously, consequently contributing to a high overall energy storage performance. Experimentally, therefore an ultrahigh energy density of 131 J cm-3 is achieved with a high efficiency of 81.6% in the microcrystal-amorphous dual-phase Bi3NdTi4O12 films. This work provides a guidance to substantially enhance dielectric energy storage by a simple and effective microstructure design.

Opportunity begets opportunity to drive macroevolutionary dynamics of a diverse lizard radiation.

Evolution proceeds unevenly across the tree of life, with some lineages accumulating diversity more rapidly than others. Explaining this disparity is challenging as similar evolutionary triggers often do not result in analogous shifts across the tree, and similar shifts may reflect different evolutionary triggers. We used a combination of approaches to directly consider such context-dependency and untangle the complex network of processes that shape macroevolutionary dynamics, focusing on Pleurodonta, a diverse radiation of lizards. Our approach shows that some lineage-wide signatures are lost when conditioned on sublineages: while viviparity appears to accelerate diversification, its effect size is overestimated by its association with the Andean mountains. Conversely, some signals that erode at broader phylogenetic scales emerge at shallower ones. Mountains, in general, do not affect speciation rates; rather, the occurrence in the Andean mountains specifically promotes diversification. Likewise, the evolution of larger sizes catalyzes diversification rates, but only within certain ecological and geographical settings. We caution that conventional methods of fitting models to entire trees may mistakenly assign diversification heterogeneity to specific factors despite evidence against their plausibility. Our study takes a significant stride toward disentangling confounding factors and identifying plausible sources of ecological opportunities in the diversification of large evolutionary radiations.

Evolutionary patterns of variations in chromosome counts and genome sizes show positive correlations with taxonomic diversity in tropical gingers.

Cytogenetic traits such as an organism's chromosome number and genome size are taxonomically critical as they are instrumental in defining angiosperm diversity. Variations in these traits can be traced to evolutionary processes such as polyploidization, although geographic variations across cytogenetic traits remain underexplored. In the pantropical monocot family Zingiberaceae (~1500 species), cytogenetic traits have been well documented; however, the role of these traits in shaping taxonomic diversity and biogeographic patterns of gingers is not known. A time-calibrated Bayesian phylogenetic tree was constructed for 290 taxa covering three of the four subfamilies in Zingiberaceae. We tested models of chromosome number and genome size evolution within the family and whether lineage age, taxonomic diversity, and distributional range explain the variations in the cytogenetic traits. Tests were carried out at two taxonomic ranks: within Zingiberaceae and within genus Hedychium using correlations, generalized linear models and phylogenetic least square models. The most frequent changes in chromosome number within Zingiberaceae were noted to be demi-polyploidization and polyploidization (~57% of the time), followed by ascending dysploidy (~27%). The subfamily Zingiberoideae showed descending dysploidy at its base, while Alpinioideae showed polyploidization at its internal nodes. Although chromosome counts and genome sizes did not corroborate with each other, suggesting that they are not equivalent; higher chromosome number variations and higher genome size variations were associated with higher taxonomic diversity and wider biogeographic distribution. Within Zingiberaceae, multiple incidences of polyploidization were discovered, and cytogenetic events appear to have reduced the genome sizes and increased taxonomic diversity, distributional ranges and invasiveness.

Prospective validation of in silico clinical trials prediction using an EGFR-mutated NSCLC mechanistic computational disease model.

8614 Background: Advanced lung adenocarcinoma (aLUAD) is divided into multiple molecularly defined subsets, each with specific biological characteristics and responses to targeted therapy. Consequently, conducting practice-changing and cost-effective randomized trials is challenging. In silico clinical trials, utilizing mechanistic computational models, offer a potential solution, yet their ability to predict clinical trial results has not been proven. Methods: We updated a previously published version of the computational mechanistic model of EGFR-mutated aLUAD, notably by incorporating metastases. Disease model was integrated with physiologically based pharmacokinetic (PBPK) models of cisplatin, pemetrexed, and osimertinib, offering a more comprehensive understanding of the disease's dynamics. These models are founded on causal connections between biological entities and phenomena, establishing a link between the specific biological behaviors of EGFR-mutated aLUAD and the evolution of tumor size and progression over time. We have independently and prospectively simulated the outcomes of the FLAURA2 and MARIPOSA randomized trials dedicated to EGFR-mutated aLUAD. This was achieved by replicating the experimental protocols and generating cohorts of virtual twins of patients using publicly available data. Results: In silico results were released before the actual outcomes of the FLAURA2 and MARIPOSA trials were publicly available. The simulations yielded predictions with overlapping confidence intervals, similar hazard ratios, median survivals, and shapes of curves compared to the actual trial data. Bootstrapped weighted log-rank tests were used to compare the simulated and observed Kaplan-Meier curves for the two arms of FLAURA2 and for the osimertinib arm of MARIPOSA. Between 94 and 98% of the tests were statistically non-significant (α=0.05), confirming the accuracy of the predictions. Conclusions: This first-of-its-kind prediction of clinical trial outcomes demonstrates that in silico trials, when based on robust mechanistic models, can be a reliable tool for enhancing the design of actual trials, particularly for EGFR-mutated aLUAD. They offer a promising avenue for overcoming the challenges of comparator arms in single-arm trials and may serve as a new standard for formulating statistical hypotheses in future studies.

Grain size sensitive modelling of the nonlinear behaviour and fatigue damage of Inconel 718 superalloy

Fatigue tests conducted on various types of materials often exhibit significant scatter, particularly in High Cycle Fatigue tests, attributed to imperfections located at the scale of the microstructure. In metallic materials, one contributor to this dispersion, but not the only one, is the variation in grain size, which cannot be uniform in a structure that has undergone a complex thermal history during its manufacturing process. This paper introduces a new model for predicting the evolution of the grain size in the nickel-based superalloy Inconel® 718 based on the applied thermal load and explores its impact on cyclic elasto-viscoplastic behaviour and fatigue life predictions. The proposed approach is applied to a complete numerical life analysis of an aircraft turbine disk.

Magnesium‐based nanocomposites developed through multi‐pass friction stir processing and strengthening mechanisms exploration

AbstractIn the current experimental work, using secondary phase hard nano titanium carbide (TiC) particles as reinforcement, two different magnesium metal matrices i. e., AZ31B/TiC and AZ61 A/TiC composite materials were synthesized by friction stir processing. Using the traditional testing approach for the developed materials, the simultaneous gain in metallurgical, mechanical, electrical, and tribological characteristics compared to the base substrate was examined. The microstructure study results for AZ31B/TiC and AZ61 A/TiC composites showed a uniform distribution of reinforced particles as well as an evolution in grain size, from 82 μm to 4.2 μm and from 74 μm to 3.7 μm, respectively, which consequently contribute in a significant gain in the microhardness of both composites i. e., around 2.2 times and 2.67 times respectively, greater than the base metal. When compared to monolithic alloys, the synthesized AZ31B/TiC and AZ61 A/TiC composites showed improvements in the areas of tensile strength, compressive strength, and coefficient of friction up to 1.81 times, and 1.64 times, 1.74 times and 1.58 times, and 57.92 % and 58.47 %, respectively. Furthermore, these improvements in characteristics also increase the final strengthening of the nanocomposite and reduce electrical conductivity.

Life history changes associated with over 400 generations of artificial selection on body size in Drosophila.

Body size is a trait that shapes many aspects of a species' development and evolution. Larger body size is often beneficial in animals, but it can also be associated with life history costs in natural systems. Similarly, miniaturization, the evolution of extremely small adult body size, is found in every major animal group, yet carries its own life history trade-offs. Given that these effects can depend on an animal's environment and life stage and have mainly been studied in species that are already specialized for their size, the life history changes associated with evolutionary shifts in body size warrant additional investigation. Here, we used Drosophila melanogaster populations that had undergone over 400 generations of artificial selection on body size to investigate the changes in life history traits associated with the evolution of extremely large and extremely small body sizes. Populations selected for small body size experienced strong trade-offs in multiple life history traits, including reduced female fecundity and lower juvenile viability. Although we found positively correlated changes in egg size associated with selection for both large and small body size, after adjusting for female body size, females from populations selected for large size had the lowest relative investment per egg and females from populations selected for small size had the highest relative investment per egg. Taken together, our results suggest that egg size may be a key constraint on the evolution of body size in D. melanogaster, providing insight into the broader phenomenon of body size evolution in insects.

Response to 223RaCl2 or 177Lu-iPSMA treatment in bone metastases due to castration-resistant prostate cancer as a function of absorbed dose

A new therapy for the treatment of bone metastases due to castration-resistant prostate cancer (mCRPC) is the administration of 223Ra-dichloride (223RaCl2) or 177Lu-iPSMA. Both radionuclides have a gamma emission, but the response to the treatment is not obtained by means of gammagraphic images. This evaluation is made from positron images that require a second radiopharmaceutical application to the patient as well as the appropriate installation and instrumentation. ObjectiveTo evaluate the response to treatment in bone metastases due to castration-resistant prostate cancer in patients treated with 223RaCl2 or 177Lu-iPSMA, by the uptake percentage and radiation absorbed dose in the metastatic lesion derived from gammagraphic images. MethodsTwo male patients with mCRCP received 4 and 3 cycles of 223RaCl2 and 177Lu-iPSMA, respectively. Lesion size evolution was assessed by absorbed dose calculated by MIRD methodology. Activity % was determined from gammagraphic images. ResultsTreatment with 223RaCl2 did not significantly reduce lesion size. However, the patient reported a significant reduction in pain. The patient treated with 177Lu-iPSMA showed a sustained reduction in lesion size. About 30% regression was achieved. ConclusionsThe proposed method avoids the additional use of positron radiopharmaceuticals and allows to correlate the response to treatment with the absorbed dose in the lesions.

Small body size is associated with increased evolutionary lability of wing skeleton proportions in birds

Birds are represented by 11,000 species and a great variety of body masses. Modular organisation of trait evolution across birds has facilitated simultaneous adaptation of different body regions to divergent ecological requirements. However, the role modularity has played in avian body size evolution, especially small-bodied, rapidly evolving and diverse avian subclades, such as hummingbirds and songbirds, is unknown. Modularity is influenced by the intersection of biomechanical restrictions, adaptation, and developmental controls, making it difficult to uncover the contributions of single factors such as body mass to skeletal organisation. We develop a novel framework to decompose this complexity, assessing factors underlying the modularity of skeletal proportions in fore-limb propelled birds distributed across a range of body masses. We demonstrate that differences in body size across birds triggers a modular reorganisation of flight apparatus proportions consistent with biomechanical expectations. We suggest weakened integration within the wing facilitates radiation in small birds. Our framework is generalisable to other groups and has the capacity to untangle the multi-layered complexity intrinsic to modular evolution.

A highly resolved nuclear phylogeny uncovers strong phylogenetic conservatism and correlated evolution of fruit color and size in Solanum L.

Mutualisms between plants and fruit-eating animals were key to the radiation of angiosperms. Still, phylogenetic uncertainties limit our understanding of fleshy-fruit evolution, as in the case of Solanum, a genus with remarkable fleshy-fruit diversity, but with unresolved phylogenetic relationships. We used 1786 nuclear genes from 247 species, including 122 newly generated transcriptomes/genomes, to reconstruct the Solanum phylogeny and examine the tempo and mode of the evolution of fruit color and size. Our analysis resolved the backbone phylogeny of Solanum, providing high support for its clades. Our results pushed back the origin of Solanum to 53.1 million years ago (Ma), with most major clades diverging between 35 and 27 Ma. Evolution of Solanum fruit color and size revealed high levels of trait conservatism, where medium-sized berries that remain green when ripe are the likely ancestral form. Our analyses revealed that fruit size and color are evolutionary correlated, where dull-colored fruits are two times larger than black/purple and red fruits. We conclude that the strong phylogenetic conservatism shown in the color and size of Solanum fruits could limit the influences of fruit-eating animals on fleshy-fruit evolution. Our findings highlight the importance of phylogenetic constraints on the diversification of fleshy-fruit functional traits.

Evolution of Litter Size: Proximate and Ultimate Mechanisms.

Relative reproductive success and failure are the ultimate determinants of Darwinian fitness. As such, reproductive traits and variation therein have an immediate and considerable impact on the evolutionary trajectory of lineages. Historically, significant attention has been paid to the ecological and evolutionary processes (ultimate factors) that shape the diversity and canalization of reproductive traits within groups to better our understanding of organismal diversity and population or species resilience. In contrast, the physiological systems that mediate variation within and among species (i.e., the proximate factors) in reproductive traits remain a significant black box. To-date, there is comparatively little information about how proximate mechanisms constrain or promote evolutionary potential in reproductive traits. In this mini-review, we focus on litter size in Eutherian mammals as a trait with relatively well-defined diversity (litter sizes are well-described both within and across species) and for which some genetic determinants have been identified. We discuss both the ultimate and potential proximate determinants of litter size with special attention to the breadth of physiological traits that may act as "toggle" switches for evolution of litter size. We close with a brief discussion of the role that physiological plasticity may play in the evolution of litter size and lay out several forward-looking areas for future research.

Evolutionary–developmental (evo-devo) dynamics of hominin brain size

Brain size tripled in the human lineage over four million years, but why this occurred remains uncertain. Here, to study what caused this brain expansion, I mathematically model the evolutionary and developmental (evo-devo) dynamics of hominin brain size. The model recovers (1) the evolution of brain and body sizes of seven hominin species starting from brain and body sizes of the australopithecine scale, (2) the evolution of the hominin brain–body allometry and (3) major patterns of human development and evolution. I show that the brain expansion recovered is not caused by direct selection for brain size but by its genetic correlation with developmentally late preovulatory ovarian follicles. This correlation is generated over development if individuals experience a challenging ecology and seemingly cumulative culture, among other conditions. These findings show that the evolution of exceptionally adaptive traits may not be primarily caused by selection for them but by developmental constraints that divert selection.

Diversity and evolution of transposable elements in the plant-parasitic nematodes

BackgroundTransposable elements (TEs) are mobile DNA sequences that propagate within genomes, occupying a significant portion of eukaryotic genomes and serving as a source of genetic variation and innovation. TEs can impact genome dynamics through their repetitive nature and mobility. Nematodes are incredibly versatile organisms, capable of thriving in a wide range of environments. The plant-parasitic nematodes are able to infect nearly all vascular plants, leading to significant crop losses and management expenses worldwide. It is worth noting that plant parasitism has evolved independently at least three times within this nematode group. Furthermore, the genome size of plant-parasitic nematodes can vary substantially, spanning from 41.5 Mbp to 235 Mbp. To investigate genome size variation and evolution in plant-parasitic nematodes, TE composition, diversity, and evolution were analysed in 26 plant-parasitic nematodes from 9 distinct genera in Clade IV.ResultsInterestingly, despite certain species lacking specific types of DNA transposons or retrotransposon superfamilies, they still exhibit a diverse range of TE content. Identification of species-specific TE repertoire in nematode genomes provides a deeper understanding of genome evolution in plant-parasitic nematodes. An intriguing observation is that plant-parasitic nematodes possess extensive DNA transposons and retrotransposon insertions, including recent sightings of LTR/Gypsy and LTR/Pao superfamilies. Among them, the Gypsy superfamilies were found to encode Aspartic proteases in the plant-parasitic nematodes.ConclusionsThe study of the transposable element (TE) composition in plant-parasitic nematodes has yielded insightful discoveries. The findings revealed that certain species exhibit lineage-specific variations in their TE makeup. Discovering the species-specific TE repertoire in nematode genomes is a crucial element in understanding the evolution of genomes in plant-parasitic nematodes. It allows us to gain a deeper insight into the intricate workings of these organisms and their genetic makeup. With this knowledge, we are gaining a fundamental piece in the puzzle of understanding the evolution of these parasites. Moreover, recent transpositions have led to the acquisition of new TE superfamilies, especially Gypsy and Pao retrotransposons, further expanding the diversity of TEs in these nematodes. Significantly, the widely distributed Gypsy superfamily possesses proteases that are exclusively associated with parasitism during nematode-host interactions.These discoveries provide a deeper understanding of the TE landscape within plant-parasitic nematodes.

High-Grade REE accumulation in regolith: Insights from supergene alteration of an apatite-rich vein at the Kapunda Cu mine, South Australia

AbstractThe stratiform and vein-hosted Kapunda Cu deposit in South Australia contains a saprolitized hydrothermal vein with 12.37 wt.% total rare earth oxide (TREO). The vein was analyzed by X-ray diffraction, scanning electron microscopy, synchrotron-based X-ray fluorescence microscopy and electron backscatter diffraction to understand the controls that govern high-grade REE accumulation during periods of intense weathering. Petrological assessments indicate the transformation of an apatite-calcite-aluminosilicate-bearing protolith to a supergene assemblage of Fe-oxides, kaolinite and REE-phosphate minerals that include rhabdophane-(Ce), monazite-(Ce) and florencite-(Ce). This transformation was facilitated by progressive acidification of the weathering fluid, which is indicated by: 1) the increasing crystallinity of authigenic Fe-oxides and kaolinite, which led to REE desorption; 2) the textural evolution and increase in grain size of authigenic REE-phosphates from nanoscopic crystallites, to acicular needles, to micro-scale hexagonal prisms; 3) the late dissolution of REE-phosphates; and 4) the replacement of goethite by jarosite, whose sulfate component originated from the oxidation and weathering of proximal sulfide minerals. Alongside the depletion of pH-buffering carbonate minerals that are indicated by the preservation of calcite menisci, this sulfide dissolution also facilitated acid generation. Results illustrate how highly acidic weathering fluids might facilitate either REE mobilization or REE accumulation in regolith. High-grade REE accumulation under acidic supergene conditions is prioritized when the host-rock contains a significant source of depositional ligands (i.e., phosphate in the form of apatite) that can be readily leached during intense weathering. Exploration companies should therefore assay routinely for REEs in any heavily weathered phosphatic rock, due to the observed efficiency of phosphate minerals as geochemical traps for REE accumulation.

Genome Assembly of the Dyeing Poison Frog Provides Insights into the Dynamics of Transposable Element and Genome-Size Evolution.

Genome size varies greatly across the tree of life and transposable elements are an important contributor to this variation. Among vertebrates, amphibians display the greatest variation in genome size, making them ideal models to explore the causes and consequences of genome size variation. However, high-quality genome assemblies for amphibians have, until recently, been rare. Here, we generate a high-quality genome assembly for the dyeing poison frog, Dendrobates tinctorius. We compare this assembly to publicly available frog genomes and find evidence for both large-scale conserved synteny and widespread rearrangements between frog lineages. Comparing conserved orthologs annotated in these genomes revealed a strong correlation between genome size and gene size. To explore the cause of gene-size variation, we quantified the location of transposable elements relative to gene features and find that the accumulation of transposable elements in introns has played an important role in the evolution of gene size in D. tinctorius, while estimates of insertion times suggest that many insertion events are recent and species-specific. Finally, we carry out population-scale mobile-element sequencing and show that the diversity and abundance of transposable elements in poison frog genomes can complicate genotyping from repetitive element sequence anchors. Our results show that transposable elements have clearly played an important role in the evolution of large genome size in D. tinctorius. Future studies are needed to fully understand the dynamics of transposable element evolution and to optimize primer or bait design for cost-effective population-level genotyping in species with large, repetitive genomes.

Time-Dependent Multi-Particle Model Describing the Hydrogen Absorption of Nanocrystalline Magnesium Powders: A Case Study

Classical kinetic models describing the hydrogen absorption of nanocrystalline metallic hydrides generally do not involve any parameter related to the change in the crystallite size during the hydrogenation at constant temperature. In the present investigation, ball-milled nanocrystalline Mg powders exhibiting lognormal crystallite size distribution have been subjected to hydrogen absorption in a Sievert-type apparatus. Partially absorbed states were achieved by interrupting the hydrogenation cycle at different hydrogen content, i.e., when 15%, 50%, and 90% of Mg powder transformed to MgH2. The evolution of the characteristic size of the nucleating MgH2 phase was determined from X-ray diffraction analysis. Considering the crystallite size distribution of the as-milled powder agglomerate as well as the growth during the isothermal hydrogenation process, a time-dependent multi-particle reaction function ∝CV¯t;R(t) was developed. It was shown unambiguously for this case study that the measured hydrogen absorption curve of the ball-milled Mg powder shows the best correlation with this model when it is compared to classical kinetic functions or the previously developed multi-particle reaction function excluding the change in the average crystallite size during hydrogenation.

Mitochondrial Genomic Evidence of Selective Constraints in Small-Bodied Terrestrial Cetartiodactyla.

Body size may drive the molecular evolution of mitochondrial genes in response to changes in energy requirements across species of different sizes. In this study, we perform selection pressure analysis and phylogenetic independent contrasts (PIC) to investigate the association between molecular evolution of mitochondrial genome protein-coding genes (mtDNA PCGs) and body size in terrestrial Cetartiodactyla. Employing selection pressure analysis, we observe that the average non-synonymous/synonymous substitution rate ratio (ω) of mtDNA PCGs is significantly reduced in small-bodied species relative to their medium and large counterparts. PIC analysis further confirms that ω values are positively correlated with body size (R2 = 0.162, p = 0.0016). Our results suggest that mtDNA PCGs of small-bodied species experience much stronger purifying selection as they need to maintain a heightened metabolic rate. On the other hand, larger-bodied species may face less stringent selective pressures on their mtDNA PCGs, potentially due to reduced relative energy expenditure per unit mass. Furthermore, we identify several genes that undergo positive selection, possibly linked to species adaptation to specific environments. Therefore, despite purifying selection being the predominant force in the evolution of mtDNA PCGs, positive selection can also occur during the process of adaptive evolution.