Homology comparison of UCP1 in Eothenomys miletus with other species
The nucleotide sequences of UCP1 in
Eothenomys miletus showed the highest similarity with those of
Myodes glareolus,
Microtus oregoni and
Arvicola amphibius, while lower similarity was observed with
Ochotona curzoniae (81.696%) and
Bos mutus (80.022%) (Table 3). Phylogenetic tree analysis aligned with sequence comparisons (Fig 1). Consequently, it is hypothesized that the
UCP1 gene is relatively conserved in rodents, yet it also exhibits specificity among different species
(Liu et al., 2016 and
Fu et al., 2006). This finding is analogous to the results of
Gao, (2014) analysis of partial sequences of the
UCP1 gene in the Chinese-Myanmar tree shrew and also supports the view of
Beáta et al., (2024) that the
UCP1 gene exhibits both conservation and species specificity between humans and rodents.
Bioinformatics analysis of UCP1 gene sequence
The
UCP1 gene in
Eothenomys miletus encodes 307 amino acids (Table 4), with a molecular weight of 33 kDa, theoretical isoelectric point (pI) of 9.23, molecular formula C
1498H
2402N
392O
434S
16, total atom count of 4,742, instability index of 32.61, extinction coefficient of 24,785, grand average of hydropathicity (GRAVY) of 0.157 and a half-life of 30 hours. The protein is alkaline, with threonine (Thr, T) being the most abundant amino acid (11.1%) (Fig 2). It contains 19 negatively charged and 28 positively charged residues. Post-translational modifications (PTMs) include 33 phosphorylation sites (15 serine [S], 17 threonine [T] and 1 tyrosine [Y]) (Fig 3), 1 potential N-glycosylation site at position 188 (Fig 4) and 5 O-glycosylation sites at positions 3, 5, 6, 7 and 12.
Eothenomys miletus is an endemic species of the Hengduan Mountains in China
(Luo et al., 2000). In order to adapt to environmental changes at different altitudes in this region, the species enhances its thermogenic capacity by upregulating UCP1 expression in brown adipose tissue, thereby better coping with cold environments
(Han et al., 2022). Furthermore, post-translational modifications (PTMs), including phosphorylation and glycosylation, have been demonstrated to play a pivotal regulatory role in the process of mRNA translation into proteins and are implicated in a variety of biological processes (
Wang and Zhang, 2019). Therefore, the present study conducted a bioinformatics analysis of the UCP1 protein in
Eothenomys miletus, systematically elucidating its basic characteristics, including amino acid composition, molecular weight, hydrophilic/hydrophobic amino acid ratio and phosphorylation and glycosylation sites. This study provides important data for further exploration of the functional role of the
UCP1 gene.
Subcellular localization and structural features
UCP1 is localized to mitochondria, lacks a GPI anchor or transmembrane domains and contains ~11.20 amino acids embedded in transmembrane regions (Fig 5). The first 60 amino acids include 5.52 transmembrane helices, with a 21.93% probability of cytoplasmic orientation. The protein exhibits strong hydrophobicity (peak value: 2.311) and weak hydrophilicity (minimum: -2.033), with dominant hydrophobic regions at positions 4-7 (Fig 6). Signal peptide analysis (D=0.214; cutoff=0.450) confirmed UCP1 as a non-secretory protein (Fig 7). Secondary structure comprises 51.79% α-helices, 37.13% random coils and 11.07% extended strands (Fig 8). Three Mito_carr domains span residues 10–107, 109–206 and 209–300 (Fig 9). Tertiary modeling (SWISS-MODEL) yielded a high-quality structure (GMQE=0.87, 82.27% sequence identity) (Fig 10). The first identification of UCP1 was in the fat mitochondria of hamsters, rats and guinea pigs in the mid-1970s, with particularly high activity in BAT (
Rodríguez-Cuenca et al., 2010). Subsequent studies revealed that UCP1 is also expressed in the thymus of mice
(Adams et al., 2008). As a pivotal mitochondrial inner membrane protein, UCP1 exerts a pivotal role in temperature regulation, energy metabolism and obesity control in animals
(Zhou et al., 2022). This study shows that the UCP1 protein in
E. miletus is localized in mitochondria. The secondary and tertiary structure prediction results of the UCP1 protein are consistent, indicating significant structural stability and functionality. This provides a foundation of data for a better under-standing of UCP1.
Selection pressure analysis
In this study, the ML method in MEGA X was used to construct the phylogenetic tree of the
UCP1 gene (Fig 11). The topology of the phylogenetic tree constructed by MEGA X was basically consistent with the traditional classification and the bootstrap values were also basically consistent, with most above 90 and a few below 70. It is inferred that the evolutionary rate of the
UCP1 gene is moderate and has basically formed a relatively stable function.
Likelihood ratio tests (2ΔlnL) revealed significant differences between models M0 and M3 (2ΔlnL = 92.3228, p<0.05) but not between M7 and M8 (p = 0.99985601). M3 indicated three site classes (p0=0.02029, p1=0.02341, p2=0.58008; ω=0, 0, 0.05405), while M1 showed two classes (p0=0.82334, p1=0.17666; ω=0.11383, 1.00000). No positive selection signals were detected. During the course of a long-term evolutionary process, selective pressure does not invariably act on the entirety of a gene; rather, it acts on specific functional sites, such as particular sites within the coding region of a gene (
Kamath and Getz, 2011). In this study, no positive selection sites were identified in the
UCP1 gene across the selected species, indicating that UCP1 is a relatively conserved gene.
A comparison of the high-latitude, low-altitude (foreground branch) model (M0 vs. M2: 2ΔlnL = 0.23402, p>0.05) with the low-latitude, high-altitude model (2ΔlnL = 7.17748, p<0.05, dN/dS = 0.15496) indicates that UCP1 in low-latitude, high-altitude rodents tends toward purifying selection (Table 5). This evolutionary phenomenon can be inferred from the environmental adaptation mechanisms of
Eothenomys miletus, a species endemic to the Hengduan Mountains. The low-latitude, high-altitude regions of the Hengduan Mountains possess unique environmental characteristics, including small annual temperature fluctuations and relatively abundant food resources
(Gong et al., 2001). This environmental stability may have reduced the selective pressure on the UCP1 gene.
Mu, (2015) observed that this species employs metabolic, regulatory and hormonal mechanisms to mitigate the effects of environmental fluctuations on the organism. Research has demonstrated that regulatory processes, encompassing metabolic and hormonal regulation, have been shown to provide adaptive contributions in response to environmental stress
(Scherbarth et al., 2010; Kuzmenko et al., 2024; Kuti et al., 2022). Consequently, these regulatory adaptive mechanisms can be employed to mitigate selection pressure on the UCP1 gene. In relatively stable environments, natural selection tends to preserve the existing functional conservatism of the UCP1 gene rather than promote new adaptive mutations. Animals have been shown to optimize the thermogenic efficiency of brown adipose tissue by maintaining the stability of the UCP1 gene. This gene may confer evolutionary advantages more effectively than developing new metabolic pathways (
Kamath and Getz, 2011).
From a functional perspective, the UCP1 gene primarily participates in energy metabolism and body temperature regulation processes, which are crucial for the survival and reproduction of the species. During the course of long-term evolution, it is imperative for this gene to preserve relatively stable functions, thereby ensuring the normal physiological activities of the organism. For instance, in cold environments, the UCP1 protein converts chemical energy into thermal energy through uncoupling oxidative phosphorylation to maintain body temperature
(Dieckmann et al., 2022). This fundamental and critical function subjects the gene to strong purifying selection, continuously eliminating harmful mutations and making it difficult for adaptive mutations represented by positive selection sites to accumulate. Consequently, this process maintains functional conservation. Furthermore, UCP1 demonstrates adaptive differentiation in thermal environments: high-latitude mammals exhibit higher evolutionary rates
(Su et al., 2025), while low-latitude species, such as
Heteroce- phalus glaber, exhibit loss-of-function mutations
(Kim et al., 2011). The UCP1 gene, located in low-latitude, high-altitude regions, functions as an evolutionary "hub," balancing ancestral thermogenesis (cold adaptation) with regulatory flexibility (environmental variability). The coexistence of thermal stress and hypoxia in high-altitude regions may drive a trade-off between thermogenesis efficiency and energy conservation, providing a natural model for studying metabolic adaptation thresholds in cross-latitude dispersal. However, further research is needed in subsequent studies to elucidate the precise function of the
UCP1 gene and its regulatory mechanisms.