Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Evolutionary trajectories”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 recordsLinked to original sources
PubMed2025

Distinct evolutionary trajectories of subgenomic centromeres in polyploid wheat.

BACKGROUND: Centromeres are crucial for precise chromosome segregation and maintaining genome stability during cell division. However, their evolutionary dynamics, particularly in polyploid organisms with complex genomic architectures, remain largely enigmatic. Allopolyploid wheat, with its well-defined hierarchical ploidy series and recent polyploidization history, serves as an excellent model to explore centromere evolution. RESULTS: In this study, we perform a systematic comparative analysis of centromeres in common wheat and its corresponding ancestral species, utilizing the latest comprehensive reference genome assembly available. Our findings reveal that wheat centromeres predominantly consist of five types of centromeric-specific retrotransposon elements (CRWs), with CRW1 and CRW2 being the most prevalent. We identify distinct evolutionary trajectories in the functional centromeres of each subgenome, characterized by variations in copy number, insertion age, and CRW composition. By utilizing CENH3-ChIP data across various ploidy levels, we uncover a series of CRW invasion events that have shaped the evolution of AA subgenome centromeres. Conversely, the evolutionary process of the DD subgenome centromeres involves their expansion from diploid to hexaploid wheat, facilitating adaptation to a larger genomic context. Integration of complete einkorn centromere assemblies and Aegilops tauschii pan-genomes further revealed subgenome-specific centromere evolutionary trajectories. By inclusion of synthetic hexaploid from S2-S3 generations, alongside 2x/6 × natural accessions, we demonstrate that DD subgenome centromere expansion represents a gradual evolutionary process rather than an immediate response to polyploidization. CONCLUSIONS: Our study provides a comprehensive landscape of centromere adaptation, evolution, and maturation, along with insights into how retrotransposon invasions drive centromere evolution in polyploid wheat.

Centromere↗
PubMed2026

Parallel evolutionary trajectories rewire enteropathogenic Escherichia coli adhesion to restore host attachment.

Enteropathogenic Escherichia coli (EPEC) causes disease in children, presenting as chronic diarrhea that can impair physical and cognitive development. The attachment of typical EPEC (tEPEC) to the gut epithelium via bundle-forming pili (BFP) is a key factor in its virulence. Yet, infections by atypical EPEC (aEPEC), which lack BFP, have become increasingly common. To investigate how aEPEC recover host-attachment in the absence of BFP, we performed experimental evolution using a non-adherent E. coli, constructed to mimic the ancestor of aEPEC, and selected adherent progeny. Highly adherent variants evolved through phase-variable activation of type I fimbriae (T1F), followed by two alternative trajectories: bacterial filamentation, which increases T1F avidity, or point mutations in the T1F adhesin FimH that enhance ligand affinity. Extending our analysis to the genomes of 327 aEPEC strains isolated from infected patients revealed that similar FimH mutations are common. We further demonstrated experimentally that these naturally occurring variants often increase epithelial-attachment. Our findings implicate T1F in aEPEC pathogenesis and suggest it may be clinically relevant for anti-adhesion therapy. More broadly, these results indicate that impaired host-attachment can be rapidly compensated by upregulating and optimizing an alternative adhesin, and that combining experimental evolution with comparative genomics can reveal evolutionary trajectories occurring in nature.

Bacterial Adhesion↗
PubMed2025

The evolutionary trajectories and gene regulatory roles of nuclear-integrated plastid DNA: clues for enhancing environmental adaptation in Caryophyllales.

Environmental stimuli can induce the transfer of chloroplast DNA to the nuclear genome, resulting in nuclear-integrated plastid DNAs (NUPTs). However, their role in plant adaptability remains unclear. Species within the Caryophyllales order, known for their adaptation to extreme environments, provide an ideal model for studying the evolutionary dynamics and functions of NUPTs. In this study, we analyzed NUPTs in 24 Caryophyllales species to investigate their evolution and regulatory roles in gene expression, particularly in response to environmental stimuli. We found significant interspecies variation in NUPT abundance, ranging from 566 insertions in Amaranthus cruentus to 3585 in Beta vulgaris, with sizes spanning from 100 bp to over 100 kb. Approximately 62% of NUPTs were inserted within the last 20 million years, while some species exhibit insertion peaks dating back 49 million years. NUPT presence/absence polymorphisms in six related species suggest that NUPT insertions and deletions are dynamic processes influenced by phylogeny. NUPTs predominantly integrate into intergenic regions but also insert into genes and promoters, with certain regions acting as hotspots. Notably, NUPTs introduce numerous environmental-responsive cis-acting elements in promoter regions. Genes with NUPT insertions in their promoters are significantly enriched for functions related to environmental response. Further luciferase assays in Spinacia oleracea demonstrated that NUPT insertions can regulate the expression of genes related to environmental responses, indicating their potential role in adaptive evolution. Overall, our study provides insights into NUPT evolution and their influence on gene function and plant adaptability to environmental stimuli.

Plastids↗
PubMed2025

Long-read sequencing of single cell-derived melanoma subclones reveals divergent and parallel genomic and epigenomic evolutionary trajectories.

Tumor evolution is driven by various mutational processes, ranging from single-nucleotide variants (SNVs) to large structural variants (SVs) to dynamic shifts in DNA methylation. Current short-read sequencing methods struggle to accurately capture the full spectrum of these genomic and epigenomic alterations due to inherent technical limitations. To overcome that, here we introduce an approach for long-read sequencing of single-cell derived subclones, and use it to profile 23 subclones of a mouse melanoma cell line, characterized with distinct growth phenotypes and treatment responses. We develop a computational framework for harmonization and joint analysis of different variant types in the evolutionary context. Uniquely, our framework enables detection of recurrent amplifications of putative driver genes, generated by independent SVs across different lineages, suggesting parallel evolution. In addition, our approach revealed gradual and lineage-specific methylation changes associated with aggressive clonal phenotypes. We also show our set of phylogeny-constrained variant calls along with openly released sequencing data can be a valuable resource for the development of new computational methods.

Journal Article↗
PubMed2026

Divergent trajectories of genome architecture and chromosome evolution in ferns and angiosperms.

Ferns and angiosperms represent the two largest vascular plant lineages but exhibit striking genomic and ecological contrasts. We investigated whether differences in genome size, chromosome architecture, GC content, and stomatal traits reveal divergent evolutionary trajectories between these lineages. We assembled the most comprehensive dataset to date, integrating genome size, chromosome number and size, GC content, and stomatal traits for over 1100 fern species and compared it with an extensive angiosperm dataset. Ferns exhibited markedly lower variability and c. 16-fold slower rates of chromosome size evolution than angiosperms. A persistent positive relationship between genome size and chromosome number in ferns suggests limited cytological post-polyploid diploidization. While ferns generally possess larger stomata, this difference disappears after accounting for genome size, indicating that nucleotypic constraints, rather than lineage-specific physiology, dictate stomatal dimensions. Both groups share a unimodal GC-genome size relationship peaking at c. 14 Gbp. Larger fern chromosomes imply lower genome-wide recombination rates, potentially limiting genetic reshuffling and adaptive potential. Our results highlight fundamentally divergent evolutionary trajectories, likely shaped by meiotic symmetry in ferns and meiotic asymmetry, possibly centromere drive, and post-polyploid diploidization in angiosperms, defining the functional and genomic landscapes of these lineages across deep evolutionary timescales.

Genome, Plant↗
PubMed2026

Evolution, Mechanisms, and Therapeutic Implications of Mobile Tetracycline Destructases.

Tetracycline destructases (TDases) pose an emerging global threat by enzymatically inactivating all generations of tetracycline (Tet) antibiotics, including last-resort agents such as tigecycline. Despite their recent identification, TDases have rapidly disseminated worldwide, largely driven by mobile genetic elements and environmental reservoirs. This review synthesizes current knowledge on TDase genomics, structural and catalytic mechanisms, ecological niches, and clinical impacts. We detail the mechanistic distinctions between type 1 and type 2 TDases, emphasizing their divergent structural configurations and substrate specificity profiles. Additionally, we examine strategies for therapeutic intervention, highlighting progress in structure-guided inhibitor development. Key gaps remain in understanding ancestral reservoirs, evolutionary trajectories, and effective surveillance strategies. Addressing these areas through integrative evolutionary, biochemical, and ecological studies is critical for mitigating the clinical spread and therapeutic impact of TDases globally.

Humans↗
PubMed2026

Competing subclones and fitness diversity shape tumor evolution across cancer types.

MOTIVATION: Intratumor heterogeneity arises from ongoing somatic evolution and complicates cancer diagnosis, prognosis, and treatment. Reconstructing evolutionary dynamics typically requires spatiotemporal samples, which are often unavailable in clinical settings. Computational approaches that can infer tumor evolutionary history from single-timepoint bulk sequencing data remain limited. RESULTS: We present estimating evolutionary events through single-timepoint sequencing (TEATIME), a novel computational framework that models tumors as mixtures of two competing cell populations: an ancestral clone with baseline fitness and a derived subclone with elevated fitness. Using cross-sectional bulk sequencing data, TEATIME estimates mutation rates, timing of subclone emergence, relative fitness, and number of generations of growth. To quantify intratumor fitness asymmetries, we introduce a novel metric-fitness diversity-which captures the imbalance between competing cell populations and serves as a measure of functional intratumor heterogeneity. Applying TEATIME to 33 tumor types from The Cancer Genome Atlas, we revealed divergent as well as convergent evolutionary patterns. Notably, we found that immune-hot microenvironments constraint subclonal expansion and limit fitness diversity. Moreover, we detected temporal dependencies in mutation acquisition, where early driver mutations in ancestral clones epistatically shape the fitness landscape, predisposing specific subclones to selective advantages. These findings underscore the importance of intratumor competition and tumor-microenvironment interactions in shaping evolutionary trajectories, driving intratumor heterogeneity. Lastly, we demonstrate that TEATIME-derived evolutionary parameters and fitness diversity offer novel prognostic insights across multiple cancer types. AVAILABILITY AND IMPLEMENTATION: R implementation of TEATIME is available on GitHub (https://github.com/liliulab/TEATIME) and Zenodo (https://zenodo.org/records/17422174).

Neoplasms↗
PubMed2026

Consistent and idiosyncratic pleiotropy in shaping genetic correlations.

Pleiotropy, the phenomenon where a single mutation influences multiple phenotypic traits, creates genetic correlations that can constrain evolutionary trajectories. Yet genetic correlations differ in their persistence: some remain stable over long evolutionary timescales, whereas others change rapidly across generations or environments. One explanation is that similar values of genetic correlation, rG, can arise from different pleiotropic architectures: broadly aligned effects across many loci, or disproportionate covariance contributions from a few large effect loci. Motivated by the distinction between vertical and horizontal pleiotropy, here, we develop a bivariate marker effect framework for recombinant mapping populations that separates candidate large covariance contributors from the polygenic background correlation, rD. We define rD as the correlation among marker effects after trimming markers with unusually large covariance contributions. rD is a trait-pair summary of how consistently small and moderate effect markers align across the genome; high rD is expected when many perturbations propagate through shared developmental, physiological, causal, or geometric structure. Applying this framework to high-dimensional yeast single-cell morphology, we show that trait pairs with similar rG can differ substantially in rD, and that a small number of candidate outlier regions can strongly influence some marker effect correlations. We then test whether rD predicts the environmental stability of genetic correlations under geldanamycin-mediated Hsp90 perturbation. Trait pairs with stronger rD show smaller absolute changes in rG. These results suggest that genetic correlations supported by a strong polygenic marker effect background are more environmentally stable than correlations shaped primarily by a few large covariance contributors.

Genetic Pleiotropy↗
PubMed2025

Compensatory evolution to DNA replication stress is robust to nutrient availability.

Evolutionary repair refers to the compensatory evolution that follows perturbations in cellular processes. While evolutionary trajectories are often reproducible, other studies suggest they are shaped by genotype-by-environment (GxE) interactions. Here, we test the predictability of evolutionary repair in response to DNA replication stress-a severe perturbation impairing the conserved mechanisms of DNA synthesis, resulting in genetic instability. We conducted high-throughput experimental evolution on Saccharomyces cerevisiae experiencing constitutive replication stress, grown under different glucose availability. We found that glucose levels impact the physiology and adaptation rate of replication stress mutants. However, the genetics of adaptation show remarkable robustness across environments. Recurrent mutations collectively recapitulated the fitness of evolved lines and are advantageous across macronutrient availability. We also identified a novel role of the mediator complex of RNA polymerase II in adaptation to replicative stress. Our results highlight the robustness and predictability of evolutionary repair mechanisms to DNA replication stress and provide new insights into the evolutionary aspects of genome stability, with potential implications for understanding cancer development.

DNA Replication↗
PubMed2026

Evolutionary diversification of invertase paralogs couples carbon metabolism and sexual reproduction in fission yeasts.

Dynamic patterns of gene gain and loss play a major role in the diversification of eukaryotes, reflecting adaptation to a broad range of ecological contexts. Reconstructing the evolutionary history of genes provides a powerful framework for understanding how functional innovation shapes life-history traits. Here we report a comprehensive analysis of gene gain and loss across the fission yeast clade, whose evolutionary trajectory remains elusive. Reductive evolution of metabolic genes is a major contributor to species diversification, as observed in other fungal taxa. Notably, we uncovered an evolutionary scenario in which an ancestral gene duplication was followed by lineage-specific loss of one or the other paralog, except in S. pombe, which retained both. We demonstrate that these paralogs encode catalytically-active invertases, named Inv1 and Inv2, with distinct enzymatic properties, localization, regulation, and physiological roles. Inv1 is a secreted enzyme subject to glucose catabolite repression and is the sole invertase required for sucrose assimilation, resembling canonical yeast invertases. In contrast, Inv2 is intracellular, constitutively expressed, and required for inducing sexual differentiation in response to nutrient availability. Overall, these findings reveal an unexpected role for carbon metabolism in modulating the haploid-diploid cycle of fission yeasts, suggesting that diversification of core metabolic functions may contribute to adaptation to environments with distinct sugar compositions.

Evolution↗
PubMed2026

Evolution and heterogeneity of lethal metastatic bladder cancer subtypes.

Histological variation is a prognostic feature of metastatic urothelial cancer1-3, but its evolutionary trajectory remains poorly defined. We developed a metastatic bladder cancer rapid autopsy programme enriched in histological subtypes4 to profile individuals with terminal disease. Here by reconstructing the evolutionary histories of patient tumours, we show that metastasis-to-metastasis seeding is the dominant pattern of cancer spread and that increased polyclonal migration predicts poor prognosis. The burden, heterogeneity and timing of genomic alterations differ markedly among histological subtypes. Plasmacytoid and neuroendocrine variants develop early driver alterations associated with shorter survival. Mutational signature analyses and experimental models demonstrated that plasmacytoid tumours uniquely use the Fanconi anaemia pathway to mitigate chemotherapy-induced genomic scarring. Single-nucleus profiling revealed mixed cell states in histological subtypes and an association between transcriptional heterogeneity and patient survival. Characterization of the tumour microenvironment uncovered distinct immune states across subtypes, with plasmacytoid tumours exhibiting immune-inflamed profiles, whereas squamous tumours are predominantly immunosuppressive. Last, we demonstrate that post-mortem cell-free DNA captures genomic and transcriptional heterogeneity of the subtypes, which provides a potential strategy for noninvasive assessment of tumour identity and aggressiveness. Our results provide new insights into how tumour heterogeneity shapes the evolutionary history of disease progression in bladder cancer histological subtypes.

Journal Article↗
PubMed2026

Phylogenetic Constraints and Environmental Filtering Jointly Drive Adaptive Evolution in Phragmites australis: From Genetic Structure to Trait Decoupling on the Mongolian Plateau.

The Mongolian Plateau, a typical arid and semi-arid zone in Eurasia, is characterized by highly heterogeneous and fragmented wetland habitats. Phragmites australis, a common wetland species in this region, exhibits remarkable adaptability. Unraveling the coordination between phylogenetic history and local environmental filtering is crucial for elucidating its adaptive mechanisms. Integrating landscape genomics and trait-based phylogenetic analyses, we analyzed transcriptome-wide SNPs, multidimensional functional traits, and environmental variables across 90 individuals from 30 natural P. australis populations. This study aims to reveal the genetic and phenotypic variation patterns underlying population genetic structure and trait variation, specifically distinguishing the roles of geographic isolation, environmental filtering, and phylogenetic history. Results reveal a significant drainage-dependent pattern in genetic structure. Populations in hydrologically connected basins show extensive admixture, whereas those in isolated endorheic basins form distinct lineages. While geographic isolation underpins genetic differentiation, environmental filtering independently explains ~33.84% of the genetic variation, driven primarily by moisture heterogeneity (precipitation seasonality and soil moisture). Crucially, we observed differentiated evolutionary trajectories across functional traits. Structural traits (e.g., plant height, leaf thickness) are phylogenetically conserved; in contrast, physiological traits (e.g., water use efficiency) are decoupled from phylogeny, showing patterns consistent with high plasticity regulated by local environments. This evolutionary decoupling strategy enables P. australis to flexibly adapt to heterogeneous habitats while maintaining structural stability. This study uncovers the synergistic mechanisms by which geographic isolation and environmental filtering jointly shape the genetic patterns of this cosmopolitan species at a regional scale, clarifies that its evolutionary responses may depend heavily on the differentiated plasticity of trait types, and provides valuable regional insights into how widespread wetland species adapt to heterogeneous environments under global change.

Mongolia Plateau↗
PubMed2025

Cross-Kingdom Genomic Conservation of Putative Human Sleep-Related Genes: Phylogenomic Evidence From Chlamydomonas reinhardtii.

Sleep is a widespread and evolutionarily conserved process observed in diverse organisms, from jellyfish to mammals, hinting at its origin as a life-supporting mechanism over 500 million years ago. Although its fundamental purpose and mechanisms remain unclear, sleep's evolution and adaptive significance continue to be debated. This study explores the evolutionary origins of sleep using Chlamydomonas reinhardtii as a model organism, identifying 112 putative sleep-related genes across species and highlighting the evolutionary conservation of sleep-regulatory pathways. Additionally, discovering uncharacterized proteins with high sequence similarity and significant e-values suggests unexplored roles in sleep regulation, underscoring the potential of C. reinhardtii to reveal new insights into the molecular basis of sleep. This study provides a foundation for identifying previously unknown sleep-associated proteins, particularly within single-celled organisms, which may offer novel perspectives on the biological role of sleep. The study demonstrates that phylogenomic analysis of diverse model organisms can expand our understanding of the evolutionary trajectory of sleep and its fundamental function, paving the way for further research in sleep biology and its health implications. Overall, the fundamental functions of sleep observed in higher animal phyla originated from its primordial activities, demonstrating an evolutionary continuum wherein more specialized tasks were integrated with sleep's essential restorative properties.

Chlamydomonas reinhardtii↗
PubMed2025

Origin and Evolution of Bacterial Periplasmic Force Transducers.

In double-membraned bacteria, non-equilibrium processes that occur at the outer membrane are typically coupled to the chemiosmotically energized inner membrane. TolA and TonB are homologous proteins which energetically couple inner membrane motor proteins to the essential processes of outer membrane stabilization and substrate import, respectively. The evolutionary trajectories of these proteins have been difficult to elucidate due to low-sequence conservation, yet they are thought to transduce force similarly. Here, this problem was addressed using structural prediction approaches to identify and annotate force transduction operons to trace their distribution and evolutionary origins. In the process, we identify a novel outer membrane-tethering system and a previously unknown family of monomeric force transducers. This approach revealed putative tolA genes, and thus the core organizational principles of the tol-pal operon throughout diverse bacterial taxa. We discovered that the α-helical structure of the periplasm-spanning domain II of TolA previously thought its hallmark, is anomalous amongst most Tol-Pal systems. This structure is mainly prevalent in γ-proteobacteria, likely in adaptation to their lifestyle. Comparison of Tol-Pal and Ton system distribution suggests that TolA emerged from a TonB paralogue and co-emerged with Pal, the outer membrane-tethering lipoprotein that functionalizes the Tol-Pal system. We also determined that TolB, the Pal-mobilizing protein, likely emerged from a family of outer membrane proteins; and CpoB, a periplasmic factor that coordinates peptidoglycan remodeling with cell division, was originally a lipoprotein present in the ancestral Tol-Pal system. The extensive conservation of the Tol-Pal system throughout Gracilicutes highlights its significance in bacterial cell biology.

Evolution, Molecular↗
PubMed2026

Divergent evolutionary strategies in spider venoms: A comparative proteomic profiling of four sympatric species from Yunnan.

Spider venoms comprise complex cocktails of bioactive molecules evolved for predation and defense, representing a valuable resource for biological research and pharmaceutical discovery. In this study, we performed a systematic analysis of venom gland extracts from four common spider species indigenous to Yunnan, China: Agelena limbata, Hippasa lycosina, Lycosa grahami, and Sinopoda pengi. Using an integrated transcriptomic and proteomic targeted profiling approach, we successfully annotated 141 distinct toxins. Comparative analysis revealed significant interspecific heterogeneity, suggesting distinct evolutionary trajectories and "weapon system economics." Both A. limbata and L. grahami exhibited a "peptide-dominant" profile anchored by neurotoxic peptides and isomerases, optimized for rapid chemical paralysis. In contrast, S. pengi displayed a distinct "protein-dominant" signature enriched with high-molecular-weight enzymes and CAP superfamily proteins, likely functioning to facilitate tissue degradation and toxin diffusion. Occupying an intermediate position, H. lycosina demonstrated a hybrid composition. These findings suggest that although these species share the same geographical range, their venom systems have undergone divergent evolutionary adaptations driven by specific ecological niches and hunting strategies. This study represents the first systematic proteomic characterization of these venom components, providing a valuable reservoir of molecular candidates while highlighting the bioinformatic nuances of analyzing whole-gland homogenates.

Animals↗
PubMed2026

Transposable element-driven expansion of enhancer RNA repertoires underlies regulatory innovation and polyploid adaptation in cereal crops.

Cereal genomes have undergone repeated polyploidization and transposable element (TE) proliferation, collectively generating complex regulatory landscapes. However, the evolutionary trajectories and functional implications of these landscapes remain largely unexplored. Using chromatin-bound RNA sequencing across seven cereal species, we systematically mapped 45,952 regulatory element transcripts (RETs), including 32,867 distal RETs corresponding to enhancer RNAs (eRNAs). Our analysis revealed that 56% of lineage-specific eRNAs originated from TE expansions, indicating that TEs serve as major reservoirs of species-specific regulatory innovation in cereals. Notably, we identified remarkable conservation in defense-related functions, root-specific expression, and TE-derived origins of eRNAs across both ancient and recent evolutionary layers of Triticeae, suggesting recurrent recruitment of TE-derived, root-associated regulatory elements throughout Triticeae evolution. Furthermore, we found that young eRNA pairs in hexaploid wheat with high sequence similarity, many originating from RLG_famc8.3 and DTC_famc4.3, exhibited pronounced root specificity and coordinated expression, suggesting targeted amplification and refinement of successful ancestral regulatory strategies established after Triticeae divergence. To facilitate community access, we developed Cereal-eRNAdb (http://bioinfo.cemps.ac.cn/Cereal-eRNAdb/), a comprehensive database integrating 69,426 eRNAs with functional annotations across 296 samples. Our findings suggest that TE-mediated innovation of root-specific eRNAs may contribute to Triticeae adaptation and provide a foundational resource for exploiting regulatory variation in cereal crop breeding.

Enhancer RNAs↗
PubMed2025

Recurrent Evolutionary Innovations in Rodent and Primate Schlafen Genes.

SCHLAFEN proteins are a large family of RNase-related enzymes carrying essential immune and developmental functions. Despite these important roles, Schlafen genes display varying degrees of evolutionary conservation in mammals. While this appears to influence their molecular activities, a detailed understanding of these evolutionary innovations is still lacking. Here, we used in-depth phylogenomic approaches to characterize the evolutionary trajectories and selective forces shaping mammalian Schlafen genes. We traced lineage-specific Schlafen amplifications and found that recent duplicates evolved under distinct selective forces, supporting repeated subfunctionalization cycles. Codon-level natural selection analyses in primates and rodents identified recurrent positive selection over Schlafen protein domains engaged in viral interactions. Combining known crystal structures and predicted protein structures, we discovered a novel class of rapidly evolving residues enriched at the contact interface of SCHLAFEN protein dimers. Our results suggest that inter-SCHLAFEN compatibilities are under strong selective pressures and are likely to impact their molecular functions. We posit that cycles of genetic conflicts with pathogens and between paralogs drove Schlafens' recurrent evolutionary innovations in mammals.

Animals↗