A genomic timescale for placental mammal evolution

CloseRead: a tool for assessing assembly errors in immunoglobulin loci applied to vertebrate long-read genome assemblies

Despite tremendous advances in long-read sequencing, some structurally complex and repeat-rich genomic regions remain challenging to assemble. Furthermore, we lack tools to assess local assembly quality, making it hard to identify problems and assess progress. Here we develop a new approach "CloseRead" for visualizing local assembly quality and diagnosing errors using multiple metrics. We apply CloseRead to evaluate how well immunoglobulin loci, paradigmatic cases of structurally complex regions, are assembled in 74 state-of-the-art vertebrate genomes. We then show that targeted, local re-assembly can correct the specific errors identified by CloseRead, highlighting the value of an iterative approach to genome assembly.

Host taxonomy and environment shapes insectivore viromes and viral spillover risks in Southwestern China

BACKGROUND

Zoonotic viruses originating from small mammals pose significant challenges to public health on a global scale. Insectivores, serving as natural reservoirs for a diverse array of zoonotic viruses, are known to carry a multitude of viral species. However, compared to the extensive research conducted on rodents (Rodentia) and bats (Chiroptera), the role of insectivores in harboring and transmitting unknown pathogens remains underexplored, which may lead to a severe underestimation of their contributions and impact to global public health.

RESULTS

This study employed a meta-transcriptomic approach to profile the viromes of 214 individual insectivores, encompassing 13 species from the families Soricidae, Erinaceidae, and Talpidae, collected across 12 counties in Yunnan Province, a recognized zoonotic hotspot. Based on virus reads, the analysis identified 42 viral families associated with vertebrates, highlighting significant virome diversity and host-specific viral tropisms among shrews, hedgehogs, and moles, along with notable geographic and environmental specificity of the viruses. Shrews exhibited greater viral richness and abundance compared to hedgehogs and moles, with variations influenced predominantly by host taxonomy, altitude, and geographic location. A total of 114 RNA-dependent RNA polymerase sequences were obtained, leading to the identification of 68 viruses, including 57 novel species. Instances of host jumping were observed in 11 viruses, with potential pathogenic viruses related to Mojiang paramyxovirus and members of the Hantaviridae family. Cross-species transmission was predominantly observed in viruses carried by shrews, while moles may play a pivotal role in facilitating viral transmission among insectivores.

CONCLUSIONS

This study enhances the understanding of the high diversity of mammalian viruses among insectivores in a relatively confined region and underscores the associations between virome composition and related zoonotic risks, providing a foundation for proactive measures to prevent and control the spillover of emerging zoonotic pathogens and potential future outbreaks. Video Abstract.

Accurate, scalable, and fully automated inference of species trees from raw genome assemblies using ROADIES

Current genome sequencing initiatives across a wide range of life forms offer significant potential to enhance our understanding of evolutionary relationships and support transformative biological and medical applications. Species trees play a central role in many of these applications; however, despite the widespread availability of genome assemblies, accurate inference of species trees remains challenging due to the limited automation, substantial domain expertise, and computational resources required by conventional methods. To address this limitation, we present ROADIES, a fully automated pipeline to infer species trees starting from raw genome assemblies. In contrast to the prominent approach, ROADIES incorporates a unique strategy of randomly sampling segments of the input genomes to generate gene trees. This eliminates the need for predefining a set of loci, limiting the analyses to a fixed number of genes, and performing the cumbersome gene annotation and/or whole genome alignment steps. ROADIES also eliminates the need to infer orthology by leveraging existing discordance-aware methods that allow multicopy genes. Using the genomic datasets from large-scale sequencing efforts across four diverse life forms (placental mammals, pomace flies, birds, and budding yeasts), we show that ROADIES infers species trees that are comparable in quality to the state-of-the-art studies but in a fraction of the time and effort, including on challenging datasets with rampant gene tree discordance and complex polyploidy. With its speed, accuracy, and automation, ROADIES has the potential to vastly simplify species tree inference, making it accessible to a broader range of scientists and applications.

Phylogenomics establishes an Early Miocene reconstruction of reef vertebrate diversity

Oceans blanket more than two-thirds of Earth's surface, yet marine biodiversity is disproportionately concentrated in coral reefs. Investigating the origins of this exceptional diversity is crucial for predicting how reefs will respond to anthropogenic disturbances. Here, we use a genome-scale dataset to reconstruct the evolutionary history of the wrasses and parrotfishes (Labridae), which rank among the most species-rich and ecologically diverse lineages of reef fishes. We show that major labrid clades experienced pulses of evolutionary innovation and accelerated diversification during the Miocene approximately 20 to 15 million years ago that the origin of no single phenotypic trait can explain. These results draw parallels to the evolutionary histories of many clades after mass extinctions and corroborate recent fossil evidence for an Early Miocene extinction event in oceanic vertebrates and changes in coral reef faunal composition. Our data provide genomic evidence for a major Early Miocene reassembly of reef faunas.

The many origins of extremophile fishes

Extremophiles survive in environments that are considered uninhabitable for most living things. The evolution of extremophiles is of great interest because of how they may have contributed to the assembly of ecosystems, yet the evolutionary dynamics that drive extremophile evolution remain obscure. Here, we investigate the evolution of extremophiles in Zoarcoidea, a lineage of over 300 species of fishes that have colonized both poles, the deep sea, and hydrothermal vents. We show that a pulse of habitat invasion occurred across over 20 different zoarcoid lineages within the last 8 million years, far after the origin of their prototypical innovation for surviving in cold water: type III antifreeze protein. Instead, a secondary burst of anatomical, physiological and life history traits and a handful of founder events in extreme ecosystems appear to have propelled zoarcoid diversification. These results decentralize the role of prototypical changes to organismal biology in shaping extremophile radiations and provide a clear example of how a combination of ancient adaptations and recent contingency shapes the origination of lineages in challenging habitats.

Complete omission of exon 21 from Slc12a2 transcripts in mice results in hearing loss

Hereditary hearing loss is highly heterogeneous. SLC12A2 is linked to autosomal dominant nonsyndromic hearing loss, DFNA78, with all the pathogenic variants affecting the exon 21. The gene encodes a cotransporter NKCC1 crucial for regulating intracellular osmotic pressure and producing endolymph in the cochlea. We generated two mouse strains with heterologous Slc12a2 variants in the splice site of the exon 21 (Em1: NM_009194.3:c.2912-2 A > G and Em2: c.2912-4_2913del). Slc12a2Em2/Em2 mice with complete skip of the exon 21 showed reduced endolymph on postnatal day 1 (P1), reduced stria vascularis (StV) and no auditory brainstem responses at 4 weeks. Reduced StV size was considered to be due to rebalance osmotic pressure, and upregulation of Cldn9 revealed by RNA-seq was considered as tissue response to repair the gaps from reduced cell sizes in the Slc12a2Em2/Em2 cochlea. Female Slc12a2Em2/+ mice also exhibited mild elevation of ABR thresholds in several sound frequencies. Slc12a2Em1/Em1 mice showed normal hearing, presumably due to sufficient cotransporter activity from the 9 bases shorter transcript by cryptic splicing. Minigene assays indicated that a single nucleotide difference between humans and mice at the 5' end of the exon 21 affects exon 21 splicing. Slc12a2Em2 mouse is proposed as a model for studying DFNA78 pathology.

Sexual Size Dimorphism Correlates With the Number of Androgen Response Elements in Mammals, But Only in Small-Bodied Species

Sexual size dimorphism is common throughout the animal kingdom, but its evolution and development remain difficult to explain given most of the genome is shared between males and females. Sex-biased regulation of genes via sex hormone signaling offers an intuitive mechanism by which males and females could develop different body sizes. One prediction of this hypothesis is that the magnitude of sexual size dimorphism scales with the number of androgen response elements or estrogen response elements, the DNA motifs to which sex hormone receptors bind. Here, we test this hypothesis using 268 mammalian species with full genome assemblies and annotations. We find that in the two smallest-bodied lineages (Chiroptera and Rodentia), sexual size dimorphism increases (male-larger) as the number of androgen response elements in a genome increases. In fact, myomorph rodents-which are especially small-bodied with high sexual size dimorphism-show an explosion of androgen receptor elements in their genomes. In contrast, the three large-bodied lineages (orders Carnivora, Cetartiodactyla, and Primates) do not show this relationship, instead following Rensch's Rule, or the observation that sexual size dimorphism increases with overall body size. One hypothesis to unify these observations is that small-bodied organisms like bats and rodents tend to reach peak reproductive fitness quickly and are more reliant on hormonal signaling to achieve sexual size dimorphism over relatively short time periods. Our study uncovers a previously unappreciated relationship between sexual size dimorphism, body size, and hormone signaling that likely varies in ways related to life history.

wQFM-TREE: highly accurate and scalable quartet-based species tree inference from gene trees

Motivation

methods are becoming increasingly popular for species tree estimation from multi-locus data in the presence of gene tree discordance. Accurate Species TRee Algorithm (ASTRAL), a leading method in this class, solves the Maximum Quartet Support Species Tree problem within a constrained solution space, while heuristics like Weighted Quartet Fiduccia-Mattheyses (wQFM) and Weighted Quartet MaxCut (wQMC) use weighted quartets and a divide-and-conquer strategy. Recent studies showed wQFM to be more accurate than ASTRAL and wQMC, though its scalability is hindered by the computational demands of explicitly generating and weighting Θ ( n 4 ) quartets. Here, we introduce wQFM-TREE, a novel summary method that enhances wQFM by avoiding explicit quartet generation and weighting, enabling its application to large datasets.

Results

Extensive simulations under diverse and challenging model conditions, with hundreds or thousands of taxa and genes, consistently demonstrate that wQFM-TREE matches or improves upon the accuracy of ASTRAL. It outperformed ASTRAL in 25 of 27 model conditions (statistically significant in 20) involving 200-1000 taxa. Moreover, applying wQFM-TREE to re-analyze the green plant dataset from the One Thousand Plant Transcriptomes Initiative produced a tree highly congruent with established evolutionary relationships of plants. wQFM-TREE's remarkable accuracy and scalability make it a strong competitor to leading methods. Its algorithmic and combinatorial innovations also enhance quartet-based computations, advancing phylogenetic estimation.

Availability and implementation

wQFM-TREE is freely available in open source form at https://github.com/abdur-rafi/wQFM-TREE.

New remarkably complete skeleton of Mixodectes reveals arboreality in a large Paleocene primatomorphan mammal following the Cretaceous-Paleogene mass extinction

Mixodectids are poorly understood placental mammals from the Paleocene of western North America that have variably been considered close relatives of euarchontan mammals (primates, dermopterans, and scandentians) with hypothesized relationships to colugos, extinct plagiomenids, and/or microsyopid plesiadapiforms. Here we describe the most complete dentally associated skeleton yet recovered for a mixodectid, specifically Mixodectes pungens from the early Paleocene of the San Juan Basin, New Mexico. A partial skull with all the teeth erupted and associated axial skeleton, forelimbs, and hind limbs, with epiphyses fused, indicate that it was a mature adult. Results from cladistic analyses incorporating new data robustly support primatomorphan (Primates + Dermoptera) affinities of Mixodectidae, but relationships within Euarchonta are less clear, with Mixodectes recovered as a stem primatomorphan, stem dermopteran, or stem primate. Analyses of postcrania suggest that M. pungens was a relatively large (~ 1.3 kg), claw-climbing arborealist capable of frequent clinging on large diameter vertical supports. With teeth suggesting an omnivorous diet that included leaves, M. pungens occupied a unique ecological niche in the early Paleocene of North America that differed from contemporary, arboreal plesiadapiforms that were smaller and more frugivorous. Euarchontans were thus a more diverse radiation in the early Cenozoic than previously appreciated.

CASTER: Direct species tree inference from whole-genome alignments

Genomes contain mosaics of discordant evolutionary histories, challenging the accurate inference of the tree of life. Although genome-wide data are routinely used for discordance-aware phylogenomic analyses, because of modeling and scalability limitations, the current practice leaves out large chunks of genomes. As more high-quality genomes become available, we urgently need discordance-aware methods to infer the tree directly from a multiple genome alignment. In this study, we introduce Coalescence-Aware Alignment-Based Species Tree Estimator (CASTER), a theoretically justified site-based method that eliminates the need to predefine recombination-free loci. CASTER is scalable to hundreds of mammalian whole genomes. We demonstrate the accuracy and scalability of CASTER in simulations that include recombination and apply CASTER to several biological datasets, showing that its per-site scores can reveal both biological and artifactual patterns of discordance across the genome.

Interpreting mammalian synonymous site conservation in light of the unwanted transcript hypothesis

Mammalian genomes are biased towards GC bases at third codon positions, likely due to a GC-biased ancestral genome and the selectively neutral recombination-related process of GC-biased gene conversion. The unwanted transcript hypothesis posits that this high GC content at synonymous sites may be beneficial for protecting against spurious transcripts, particularly in species with low effective population sizes. Utilising a 240 placental mammal genome alignment and single-base resolution conservation scores, we interpret sequence conservation at mammalian four-fold degenerate sites in this context and find evidence in support of the unwanted transcript hypothesis, including a strong GC bias, high conservation at sites relating to exon splicing, less human genetic variation at conserved four-fold degenerate sites, and conservation of sites important for epigenetic regulation of developmental genes. Additionally, we show that high conservation of four-fold degenerate sites in essential developmental genes, including homeobox genes, likely relates to the low mutation rates experienced by these genes.

Identification of the shortest species-specific oligonucleotide sequences

Despite the exponential increase in sequencing information driven by massively parallel DNA sequencing technologies, universal and succinct genomic fingerprints for each organism are still missing. Identifying the shortest species-specific nucleotide sequences offers insights into species evolution and holds potential practical applications in agriculture, wildlife conservation, and healthcare. We propose a new method for sequence analysis termed nucleic "quasi-primes," the shortest occurring sequences in each of 45,076 organismal reference genomes, present in one genome and absent from every other examined genome. In the human genome, we find that the genomic loci of nucleic quasi-primes are most enriched for genes associated with brain development and cognitive function. In a single-cell case study focusing on the human primary motor cortex, nucleic quasi-prime genes account for a significantly larger proportion of the variation based on average gene expression. Nonneuronal cell types, including astrocytes, endothelial cells, microglia perivascular-macrophages, oligodendrocytes, and vascular and leptomeningeal cells, exhibit significant activation of quasi-prime-containing gene associations related to cancer, whereas simultaneously suppressing quasi-prime-containing genes are associated with cognitive, mental, and developmental disorders. We also show that human disease-causing variants, eQTLs, mQTLs, and sQTLs are 4.43-fold, 4.34-fold, 4.29-fold, and 4.21-fold enriched at human quasi-prime loci, respectively. These findings indicate that nucleic quasi-primes are genomic loci linked to the evolution of species-specific traits, and in humans, they provide insights in the development of cognitive traits and human diseases, including neurodevelopmental disorders.

The human genome encodes a multitude of novel miRNAs

Human cells generate a vast complexity of noncoding RNAs, the "RNA dark matter," which includes a vast small RNA (sRNA) transcriptome. The biogenesis, biological relevance, and mechanisms of action of most of these transcripts remain unknown, and they are widely assumed to represent degradation products. Here, we aimed to functionally characterize human sRNA transcriptome by attempting to answer the following question-can a significant number of novel sRNAs correspond to novel members of known classes, specifically, microRNAs (miRNAs)? By developing and validating a miRNA discovery pipeline, we show that at least 2726 novel canonical miRNAs, majority of which represent novel miRNA families, exist in just one human cell line compared to just 1914 known miRNA loci. Moreover, potentially tens of thousands of miRNAs remain to be discovered. Strikingly, many novel miRNAs map to exons of protein-coding genes emphasizing a complex and interleaved architecture of the genome. The existence of so many novel members of a functional class of sRNAs suggest that the human sRNA transcriptome harbors a multitude of novel regulatory molecules. Overall, these results suggest that we are at the very beginning of understanding the true functional complexity of the sRNA component of the "RNA dark matter."

Simplicity within biological complexity

Motivation

Heterogeneous, interconnected, systems-level, molecular (multi-omic) data have become increasingly available and key in precision medicine. We need to utilize them to better stratify patients into risk groups, discover new biomarkers and targets, repurpose known and discover new drugs to personalize medical treatment. Existing methodologies are limited and a paradigm shift is needed to achieve quantitative and qualitative breakthroughs.

Results

In this perspective paper, we survey the literature and argue for the development of a comprehensive, general framework for embedding of multi-scale molecular network data that would enable their explainable exploitation in precision medicine in linear time. Network embedding methods (also called graph representation learning) map nodes to points in low-dimensional space, so that proximity in the learned space reflects the network's topology-function relationships. They have recently achieved unprecedented performance on hard problems of utilizing few omic data in various biomedical applications. However, research thus far has been limited to special variants of the problems and data, with the performance depending on the underlying topology-function network biology hypotheses, the biomedical applications, and evaluation metrics. The availability of multi-omic data, modern graph embedding paradigms and compute power call for a creation and training of efficient, explainable and controllable models, having no potentially dangerous, unexpected behaviour, that make a qualitative breakthrough. We propose to develop a general, comprehensive embedding framework for multi-omic network data, from models to efficient and scalable software implementation, and to apply it to biomedical informatics, focusing on precision medicine and personalized drug discovery. It will lead to a paradigm shift in the computational and biomedical understanding of data and diseases that will open up ways to solve some of the major bottlenecks in precision medicine and other domains.

The Genomic Landscape, Causes, and Consequences of Extensive Phylogenomic Discordance in Murine Rodents

A species tree is a central concept in evolutionary biology whereby a single branching phylogeny reflects relationships among species. However, the phylogenies of different genomic regions often differ from the species tree. Although tree discordance is widespread in phylogenomic studies, we still lack a clear understanding of how variation in phylogenetic patterns is shaped by genome biology or the extent to which discordance may compromise comparative studies. We characterized patterns of phylogenomic discordance across the murine rodents-a large and ecologically diverse group that gave rise to the laboratory mouse and rat model systems. Combining recently published linked-read genome assemblies for seven murine species with other available rodent genomes, we first used ultraconserved elements (UCEs) to infer a robust time-calibrated species tree. We then used whole genomes to examine finer-scale patterns of discordance across ∼12 million years of divergence. We found that proximate chromosomal regions tended to have more similar phylogenetic histories. There was no clear relationship between local tree similarity and recombination rates in house mice, but we did observe a correlation between recombination rates and average similarity to the species tree. We also detected a strong influence of linked selection whereby purifying selection at UCEs led to appreciably less discordance. Finally, we show that assuming a single species tree can result in substantial deviation from the results with gene trees when testing for positive selection under different models. Collectively, our results highlight the complex relationship between phylogenetic inference and genome biology and underscore how failure to account for this complexity can mislead comparative genomic studies.

Evolution of 3D Chromatin Folding

Studies examining the evolution of genomes have focused mainly on sequence conservation. However, the inner working of a cell implies tightly regulated crosstalk between complex gene networks controlled by small dispersed regulatory elements of physically contacting DNA regions. How these different levels of chromatin organization crosstalk in different species underpins the potential for genome evolutionary plasticity. We review the evolution of chromatin organization across the Animal Tree of Life. We introduce general aspects of the mode and tempo of genome evolution to later explore the multiple layers of genome organization. We argue that both genome and chromosome size modulate patterns of chromatin folding and that chromatin interactions facilitate the formation of lineage-specific chromosomal reorganizations, especially in germ cells. Overall, analyzing the mechanistic forces involved in the maintenance of chromatin structure and function of the germ line is critical for understanding genome evolution, maintenance, and inheritance.

WASTER: Practical de novo phylogenomics from low-coverage short reads

The advent of affordable whole-genome sequencing has spurred numerous large-scale projects aimed at inferring the tree of life, yet achieving a complete species-level phylogeny remains a distant goal due to significant costs and computational demands. Traditional species tree inference methods, though effective, are hampered by the need for high-coverage sequencing, high-quality genomic alignments, and extensive computational resources. To address these challenges, this study introduces WASTER, a novel de novo tool for inferring species trees directly from short-read sequences. WASTER employs a k-mer based approach for identifying variable sites, circumventing the need for genome assembly and alignment. Using simulations, we demonstrate that WASTER achieves accuracy comparable to that of traditional alignment-based methods, even for low sequencing depth, and has substantially higher accuracy than other alignment-free methods. We validate WASTER's efficacy on real data, where it accurately reconstructs phylogenies of eukaryotic species with as low depth as 1.5X. WASTER provides a fast and efficient solution for phylogeny estimation in cases where genome assembly and/or alignment may bias analyses or is challenging, for example due to low sequencing depth. It also provides a method for generating guide trees for tree-based alignment algorithms. WASTER's ability to accurately estimate trees from low-coverage sequencing data without relying on assembly and alignment will lead to substantially reduced sequencing and computational costs in phylogenomic projects.

Functional differences between rodent and human PD-1 linked to evolutionary divergence

Mechanistic understanding of the inhibitory immunoreceptor PD-1 is largely based on mouse models, but human and mouse PD-1 share only 59.6% amino acid identity. Here, we found that human PD-1 is more inhibitory than mouse PD-1, owing to stronger interactions with the ligands PD-L1 and PD-L2 and more efficient recruitment of the effector phosphatase Shp2. In a mouse melanoma model with adoptively transferred T cells, humanization of a PD-1 intracellular domain disrupted the antitumor activity of CD8+ T cells and increased the magnitude of anti-PD-1 response. We identified a motif highly conserved across vertebrate PD-1 orthologs, absent in rodents, as a key determinant for differential Shp2 recruitment. Evolutionary analysis suggested that PD-1 underwent a rodent lineage-specific functional attenuation during evolution. Together, our study uncovers species-specific features of the PD-1 pathway, with implications for PD-1 evolution and differential anti-PD-(L)1 responses in mouse models and human patients.

Detecting diversifying selection for a trait from within and between-species genotypes and phenotypes

To quantify selection acting on a trait, methods have been developed using either within or between-species variation. However, methods using within-species variation do not integrate the changes at the macro-evolutionary scale. Conversely, current methods using between-species variation usually discard within-species variation, thus not accounting for processes at the micro-evolutionary scale. The main goal of this study is to define a neutrality index for a quantitative trait, by combining within- and between-species variation. This neutrality index integrates nucleotide polymorphism and divergence for normalizing trait variation. As such, it does not require estimation of population size nor of time of speciation for normalization. Our index can be used to seek deviation from the null model of neutral evolution, and test for diversifying selection. Applied to brain mass and body mass at the mammalian scale, we show that brain mass is under diversifying selection. Finally, we show that our test is not sensitive to the assumption that population sizes, mutation rates and generation time are constant across the phylogeny, and automatically adjust for it.

The allometry of brain size in Euarchontoglires: clade-specific patterns and their impact on encephalization quotients

The timing and nature of evolutionary shifts in the relative brain size of Primates have been extensively studied. Less is known, however, about the scaling of the brain-to-body size in their closest living relatives, i.e., among other members of Euarchontoglires (Dermoptera, Scandentia, Lagomorpha, Rodentia). Ordinary least squares (OLS), reduced major axis (RMA), and phylogenetic generalized least squares (PGLS) regressions were fitted to the largest euarchontogliran data set of brain and body mass, comprising 715 species. Contrary to previous inferences, lagomorph brain sizes (PGLS slope = 0.465; OLS slope = 0.593) scale relative to body mass similarly to rodents (PGLS = 0.526; OLS = 0.638), and differently than primates (PGLS = 0.607; OLS = 0.794). There is a shift in the pattern of the scaling of the brain in Primates, with Strepsirrhini occupying an intermediate stage similar to Scandentia but different from Rodentia and Lagomorpha, while Haplorhini differ from all other groups in the OLS and RMA analyses. The unique brain-body scaling relationship of Primates among Euarchontoglires illustrates the need for clade-specific metrics for relative brain size (i.e., encephalization quotients; EQs) for more restricted taxonomic entities than Mammalia. We created clade-specific regular and phylogenetically adjusted EQ equations at superordinal, ordinal, and subordinal levels. When using fossils as test cases, our results show that generalized mammalian equations underestimate the encephalization of the stem lagomorph Megalagus turgidus in the context of lagomorphs, overestimate the encephalization of the stem primate Microsyops annectens and the early euprimate Necrolemur antiquus, but provide similar EQ values as our new strepsirrhine-specific EQ when applied to the early euprimate Adapis parisiensis.

Genomics, the diversification of mammals, and the evolution of placentation

When and why did variations in placental structure and function evolve? Such questions cannot be addressed without a reliable version of mammalian phylogeny. Twenty-five years ago, the mammalian tree was reshaped by molecular phylogenetics. Soon it was shown, in contrast to prevailing theories, that the common ancestor of placental mammals had invasive placentation. Subsequently, evolution of many other features of extraembryonic membranes was addressed. This endeavour stimulated research to fill gaps in our knowledge of placental morphology. Last year the mammalian tree was again revised based on a large set of genomic data. With that in mind, this review provides an update on placentation in the nineteen orders of placental mammals, incorporating much recent data. The principal features such as shape, interdigitation, the interhaemal barrier and the yolk sac are summarized in synoptic tables. The evolution of placental traits and its timing is then explored by reference to the revised mammalian tree. Examples are the early appearance of epitheliochorial placentation in the common ancestor of artiodactyls, perissodactyls, pangolins and carnivores (with reversion to invasive forms in the latter) and later refinements such as the binucleate trophoblast cells and placentomes of ruminants. In primates, the intervillous space gradually evolved from the more basic labyrinth whereas trophoblast invasion of the decidua was a late development in humans and great apes. Only seldom can we glimpse the "why" of placental evolution. The best examples concern placental hormones, including some striking examples of convergent evolution such as the chorionic gonadotropins of primates and equids. In concluding, I review current ideas about what drives placental evolution and identify significant gaps in our knowledge of placentation, including several relevant to the evolution of placentation in primates.

Estimating the proportion of beneficial mutations that are not adaptive in mammals

Mutations can be beneficial by bringing innovation to their bearer, allowing them to adapt to environmental change. These mutations are typically unpredictable since they respond to an unforeseen change in the environment. However, mutations can also be beneficial because they are simply restoring a state of higher fitness that was lost due to genetic drift in a stable environment. In contrast to adaptive mutations, these beneficial non-adaptive mutations can be predicted if the underlying fitness landscape is stable and known. The contribution of such non-adaptive mutations to molecular evolution has been widely neglected mainly because their detection is very challenging. We have here reconstructed protein-coding gene fitness landscapes shared between mammals, using mutation-selection models and a multi-species alignments across 87 mammals. These fitness landscapes have allowed us to predict the fitness effect of polymorphisms found in 28 mammalian populations. Using methods that quantify selection at the population level, we have confirmed that beneficial non-adaptive mutations are indeed positively selected in extant populations. Our work confirms that deleterious substitutions are accumulating in mammals and are being reverted, generating a balance in which genomes are damaged and restored simultaneously at different loci. We observe that beneficial non-adaptive mutations represent between 15% and 45% of all beneficial mutations in 24 of 28 populations analyzed, suggesting that a substantial part of ongoing positive selection is not driven solely by adaptation to environmental change in mammals.

Patterns of chromosome evolution in ruminants

Studying when and where gross genomic rearrangements occurred during evolution is key to understanding changes in genome structure with functional consequences that might eventually lead to speciation. Here we identified chromosome rearrangements in ruminants, a clade characterized by large chromosome differences. Using 26 genome assemblies, we reconstructed five ancestral karyotypes and classified the rearrangement events occurring in each lineage. With these reconstructions, we then identified evolutionary breakpoints regions (EBRs) and synteny fragments. Ruminant karyotype evolution is characterized by inversions, while interchromosomal rearrangements occurred preferentially in the oldest ancestor of ruminants. We found that EBRs are depleted of protein coding genes, including housekeeping genes. Similarly, EBRs are not enriched in high GC regions, suggesting that meiotic double strand breaks might not be their origin. Overall, our results characterize at fine detail the location of chromosome rearrangements in ruminant evolution and provide new insights into the formation of EBRs.

A versatile attention-based neural network for chemical perturbation analysis and its potential to aid surgical treatment: an experimental study

Deep learning models have emerged as rapid, accurate, and effective approaches for clinical decisions. Through a combination of drug screening and deep learning models, drugs that may benefit patients before and after surgery can be discovered to reduce the risk of complications or speed recovery. However, most existing drug prediction methods have high data requirements and lack interpretability, which has a limited role in adjuvant surgical treatment. To address these limitations, the authors propose the attention-based convolution transpositional interfusion network (ACTIN) for flexible and efficient drug discovery. ACTIN leverages the graph convolution and the transformer mechanism, utilizing drug and transcriptome data to assess the impact of chemical pharmacophores containing certain elements on gene expression. Remarkably, just with only 393 training instances, only one-tenth of the other models, ACTIN achieves state-of-the-art performance, demonstrating its effectiveness even with limited data. By incorporating chemical element embedding disparity and attention mechanism-based parameter analysis, it identifies the possible pharmacophore containing certain elements that could interfere with specific cell lines, which is particularly valuable for screening useful pharmacophores for new drugs tailored to adjuvant surgical treatment. To validate its reliability, the authors conducted comprehensive examinations by utilizing transcriptome data from the lung tissue of fatal COVID-19 patients as additional input for ACTIN, the authors generated novel lead chemicals that align with clinical evidence. In summary, ACTIN offers insights into the perturbation biases of elements within pharmacophore on gene expression, which holds the potential for guiding the development of new drugs that benefit surgical treatment.

Reticulate evolution: Detection and utility in the phylogenomics era

Phylogenomics has enriched our understanding that the Tree of Life can have network-like or reticulate structures among some taxa and genes. Two non-vertical modes of evolution - hybridization/introgression and horizontal gene transfer - deviate from a strictly bifurcating tree model, causing non-treelike patterns. However, these reticulate processes can produce similar patterns to incomplete lineage sorting or recombination, potentially leading to ambiguity. Here, we present a brief overview of a phylogenomic workflow for inferring organismal histories and compare methods for distinguishing modes of reticulate evolution. We discuss how the timing of coalescent events can help disentangle introgression from incomplete lineage sorting and how horizontal gene transfer events can help determine the relative timing of speciation events. In doing so, we identify pitfalls of certain methods and discuss how to extend their utility across the Tree of Life. Workflows, methods, and future directions discussed herein underscore the need to embrace reticulate evolutionary patterns for understanding the timing and rates of evolutionary events, providing a clearer view of life's history.

The origin and evolution of shrews (Soricidae, Mammalia)

Shrews are among the most speciose of mammalian clades, but their evolutionary history is poorly understood. Their fossil record is fragmentary and even the anatomy of living groups is not well documented. Here, we incorporate the oldest, most complete fossil shrew yet known into the first phylogenetic analysis of the group to include molecular, morphological and temporal data. Our study reveals previously unknown diversity among total- and crown-group soricids. This includes a novel element of the mammalian skeleton: a robust, needle-like sesamoid extending cranially from the second thoracic neural arch in Myosoricini, comparable in length to these species' humeri. Additionally, 'red-toothed' shrews have an unusually elongate basicranium, and 'white-toothed' shrews probably evolved from a common ancestor with dental pigmentation. The fossil Domnina and crown soricids have a double-jaw articulation and incomplete zygomatic arch, but unlike nearly all crown species, Domnina has open vomeronasal canals and a tympanic process of the basisphenoid. Domnina and other heterosoricids are phylogenetically outside crown Soricidae. The oldest, well-supported total-group soricoids are North American, not Asian, and Soricidae probably originated during the Palaeocene or early Eocene. The diverse mammalian genus Crocidura originated and began to diversify during the Miocene.

Extensive longevity and DNA virus-driven adaptation in nearctic Myotis bats

The genus Myotis is one of the largest clades of bats, and exhibits some of the most extreme variation in lifespans among mammals alongside unique adaptations to viral tolerance and immune defense. To study the evolution of longevity-associated traits and infectious disease, we generated near-complete genome assemblies and cell lines for 8 closely related species of Myotis. Using genome-wide screens of positive selection, analyses of structural variation, and functional experiments in primary cell lines, we identify new patterns of adaptation contributing to longevity, cancer resistance, and viral interactions in bats. We find that Myotis bats have some of the most significant variation in cancer risk across mammals and demonstrate a unique DNA damage response in primary cells of the long-lived M. lucifugus. We also find evidence of abundant adaptation in response to DNA viruses - but not RNA viruses - in Myotis and other bats in sharp contrast with other mammals, potentially contributing to the role of bats as reservoirs of zoonoses. Together, our results demonstrate how genomics and primary cells derived from diverse taxa uncover the molecular bases of extreme adaptations in non-model organisms.

Relaxed selection in evolution of genes regulating limb development gives clue to variation in forelimb morphology of cetaceans and other mammals

Cetaceans have evolved unique limb structures, such as flippers, due to genetic changes during their transition to aquatic life. However, the full understanding of the genetic and evolutionary mechanisms behind these changes is still developing. By examining 25 limb-related protein-coding genes across various mammalian species, we compared genetic changes between aquatic mammals, like whales, and other mammals with unique limb structures such as bats, rodents and elephants. Our findings revealed significant modifications in limb-related genes, including variations in the Hox, GDF5 and Evx genes. Notably, a relaxed selection in several key genes was observed, suggesting a lifting of developmental constraints, which might have facilitated the emergence of morphological innovations in cetacean limb morphology. We also uncovered non-synonymous changes, insertions and deletions in these genes, particularly in the polyalanine tract of HOXD13, which are distinctive to cetaceans or convergent with other aquatic mammals. These genetic variations correlated with the diverse and specialized limb structures observed in cetaceans, indicating a complex interplay of relaxed selection and specific mutations in mammalian limb evolution.

The Primate Major Histocompatibility Complex: An Illustrative Example of Gene Family Evolution

Gene families are groups of evolutionarily-related genes. One large gene family that has experienced rapid evolution is the Major Histocompatibility Complex (MHC), whose proteins serve critical roles in innate and adaptive immunity. Across the ~60 million year history of the primates, some MHC genes have turned over completely, some have changed function, some have converged in function, and others have remained essentially unchanged. Past work has typically focused on identifying MHC alleles within particular species or comparing gene content, but more work is needed to understand the overall evolution of the gene family across species. Thus, despite the immunologic importance of the MHC and its peculiar evolutionary history, we lack a complete picture of MHC evolution in the primates. We readdress this question using sequences from dozens of MHC genes and pseudogenes spanning the entire primate order, building a comprehensive set of gene and allele trees with modern methods. Overall, we find that the Class I gene subfamily is evolving much more quickly than the Class II gene subfamily, with the exception of the Class II MHC-DRB genes. We also pay special attention to the often-ignored pseudogenes, which we use to reconstruct different events in the evolution of the Class I region. We find that despite the shared function of the MHC across species, different species employ different genes, haplotypes, and patterns of variation to achieve a successful immune response. Our trees and extensive literature review represent the most comprehensive look into MHC evolution to date.

PTMoreR-enabled cross-species PTM mapping and comparative phosphoproteomics across mammals

To support PTM proteomic analysis and annotation in different species, we developed PTMoreR, a user-friendly tool that considers the surrounding amino acid sequences of PTM sites during BLAST, enabling a motif-centric analysis across species. By controlling sequence window similarity, PTMoreR can map phosphoproteomic results between any two species, perform site-level functional enrichment analysis, and generate kinase-substrate networks. We demonstrate that the majority of real P-sites in mice can be inferred from experimentally derived human P-sites with PTMoreR mapping. Furthermore, the compositions of 129 mammalian phosphoproteomes can also be predicted using PTMoreR. The method also identifies cross-species phosphorylation events that occur on proteins with an increased tendency to respond to the environmental factors. Moreover, the classic kinase motifs can be extracted across mammalian species, offering an evolutionary angle for refining current motifs. PTMoreR supports PTM proteomics in non-human species and facilitates quantitative phosphoproteomic analysis.

Distinct Genes with Similar Functions Underlie Convergent Evolution in Myotis Bat Ecomorphs

Convergence offers an opportunity to explore to what extent evolution can be predictable when genomic composition and environmental triggers are similar. Here, we present an emergent model system to study convergent evolution in nature in a mammalian group, the bat genus Myotis. Three foraging strategies-gleaning, trawling, and aerial hawking, each characterized by different sets of phenotypic features-have evolved independently multiple times in different biogeographic regions in isolation for millions of years. To investigate the genomic basis of convergence and explore the functional genomic changes linked to ecomorphological convergence, we sequenced and annotated 17 new genomes and screened 16,426 genes for positive selection and associations between relative evolutionary rates and foraging strategies across 30 bat species representing all Myotis ecomorphs across geographic regions as well as among sister groups. We identify genomic changes that describe both phylogenetic and ecomorphological trends. We infer that colonization of new environments may have first required changes in genes linked to hearing sensory perception, followed by changes linked to fecundity and development, metabolism of carbohydrates, and heme degradation. These changes may be linked to prey acquisition and digestion and match phylogenetic trends. Our findings also suggest that the repeated evolution of ecomorphs does not always involve changes in the same genes but rather in genes with the same molecular functions such as developmental and cellular processes.

A molecular phylogeny for all 21 families within Chiroptera (bats)

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.

Evolution across the adaptive landscape in a hyperdiverse beetle radiation

The extraordinary diversification of beetles on Earth is a textbook example of adaptive evolution. Yet, the tempo and drivers of this super-radiation remain largely unclear. Here, we address this problem by investigating macroevolutionary dynamics in darkling beetles (Coleoptera: Tenebrionidae), one of the most ecomorphologically diverse beetle families (with over 30,000 species). Using multiple genomic datasets and analytical approaches, we resolve the long-standing inconsistency over deep relationships in the family. In conjunction with a landmark-based dataset of body shape morphology, we show that the evolutionary history of darkling beetles is marked by ancient rapid radiations, frequent ecological transitions, and rapid bursts of morphological diversification. On a global scale, our analyses uncovered a significant pulse of phenotypic diversification proximal to the Cretaceous-Palaeogene (K/Pg) mass extinction and convergence of body shape associated with recurrent ecological specializations. On a regional scale, two major Australasian radiations, the Adeliini and the Heleine clade, exhibited contrasting patterns of ecomorphological diversification, representing phylogenetic niche conservatism versus adaptive radiation. Our findings align with the Simpsonian model of adaptive evolution across the macroevolutionary landscape and highlight a significant role of ecological opportunity in driving the immense ecomorphological diversity in a hyperdiverse beetle group.

Assessing Assembly Errors in Immunoglobulin Loci: A Comprehensive Evaluation of Long-read Genome Assemblies Across Vertebrates

Long-read sequencing technologies have revolutionized genome assembly producing near-complete chromosome assemblies for numerous organisms, which are invaluable to research in many fields. However, regions with complex repetitive structure continue to represent a challenge for genome assembly algorithms, particularly in areas with high heterozygosity. Robust and comprehensive solutions for the assessment of assembly accuracy and completeness in these regions do not exist. In this study we focus on the assembly of biomedically important antibody-encoding immunoglobulin (IG) loci, which are characterized by complex duplications and repeat structures. High-quality full-length assemblies for these loci are critical for resolving haplotype-level annotations of IG genes, without which, functional and evolutionary studies of antibody immunity across vertebrates are not tractable. To address these challenges, we developed a pipeline, "CloseRead", that generates multiple assembly verification metrics for analysis and visualization. These metrics expand upon those of existing quality assessment tools and specifically target complex and highly heterozygous regions. Using CloseRead, we systematically assessed the accuracy and completeness of IG loci in publicly available assemblies of 74 vertebrate species, identifying problematic regions. We also demonstrated that inspecting assembly graphs for problematic regions can both identify the root cause of assembly errors and illuminate solutions for improving erroneous assemblies. For a subset of species, we were able to correct assembly errors through targeted reassembly. Together, our analysis demonstrated the utility of assembly assessment in improving the completeness and accuracy of IG loci across species.

New member of Plasmodium (Vinckeia) and Plasmodium cyclopsi discovered in bats in Sierra Leone - nuclear sequence and complete mitochondrial genome analyses

Malaria remains the most important arthropod-borne infectious disease globally. The causative agent, Plasmodium, is a unicellular eukaryote that develops inside red blood cells. Identifying new Plasmodium parasite species that infect mammalian hosts can shed light on the complex evolution and diversity of malaria parasites. Bats feature a high diversity of microorganisms including seven separate genera of malarial parasites. Three species of Plasmodium have been reported so far, for which scarce reports exist. Here we present data from an investigation of Plasmodium infections in bats in the western Guinean lowland forest in Sierra Leone. We discovered a new Plasmodium parasite in the horseshoe bat Rhinolophus landeri. Plasmodium cyclopsi infections in a member of leaf-nosed bats, Doryrhina cyclops, exhibited a high prevalence of 100%. Phylogenetic analysis of complete mitochondrial genomes and nine nuclear markers recovered a close relationship between P. cyclopsi and the new Plasmodium parasite with the rodent species Plasmodium berghei, a widely used in vivo model to study malaria in humans. The data suggests that the "rodent/bat" Plasmodium (Vinckeia) clade represents a diverse group of malarial parasites that would likely expand with a systematic sampling of small mammals in tropical Africa. Identifying the bat Plasmodium repertoire is central to our understanding of the evolution of Plasmodium parasites in mammals.

African elephants address one another with individually specific name-like calls

Personal names are a universal feature of human language, yet few analogues exist in other species. While dolphins and parrots address conspecifics by imitating the calls of the addressee, human names are not imitations of the sounds typically made by the named individual. Labelling objects or individuals without relying on imitation of the sounds made by the referent radically expands the expressive power of language. Thus, if non-imitative name analogues were found in other species, this could have important implications for our understanding of language evolution. Here we present evidence that wild African elephants address one another with individually specific calls, probably without relying on imitation of the receiver. We used machine learning to demonstrate that the receiver of a call could be predicted from the call's acoustic structure, regardless of how similar the call was to the receiver's vocalizations. Moreover, elephants differentially responded to playbacks of calls originally addressed to them relative to calls addressed to a different individual. Our findings offer evidence for individual addressing of conspecifics in elephants. They further suggest that, unlike other non-human animals, elephants probably do not rely on imitation of the receiver's calls to address one another.

Torpor use in the wild by one of the world's largest bats

Torpor is widespread among bats presumably because most species are small, and torpor greatly reduces their high mass-specific resting energy expenditure, especially in the cold. Torpor has not been recorded in any bat species larger than 50 g, yet in theory could be beneficial even in the world's largest bats (flying-foxes; Pteropus spp.) that are exposed to adverse environmental conditions causing energy bottlenecks. We used temperature telemetry to measure body temperature in wild-living adult male grey-headed flying-foxes (P. poliocephalus; 799 g) during winter in southern Australia. We found that all individuals used torpor while day-roosting, with minimum body temperature reaching 27°C. Torpor was recorded following a period of cool, wet and windy weather, and on a day with the coldest maximum air temperature, suggesting it is an adaptation to reduce energy expenditure during periods of increased thermoregulatory costs and depleted body energy stores. A capacity for torpor among flying-foxes has implications for understanding their distribution, behavioural ecology and life history. Furthermore, our discovery increases the body mass of bats known to use torpor by more than tenfold and extends the documented use of this energy-saving strategy under wild conditions to all bat superfamilies, with implications for the evolutionary maintenance of torpor among bats and other mammals.

Priority effects transcend scales and disciplines in biology

Although primarily studied through the lens of community ecology, phenomena consistent with priority effects appear to be widespread across many different scenarios spanning a broad range of spatial, temporal, and biological scales. However, communication between these research fields is inconsistent and has resulted in a fragmented co-citation landscape, likely due to the diversity of terms used to refer to priority effects across these fields. We review these related terms, and the biological contexts in which they are used, to facilitate greater cross-disciplinary cohesion in research on priority effects. In breaking down these semantic barriers, we aim to provide a framework to better understand the conditions and mechanisms of priority effects, and their consequences across spatial and temporal scales.

The Chromosome-Scale Genome of Chitala ornata Illuminates the Evolution of Early Teleosts

Teleosts are the most prolific vertebrates, occupying the vast majority of aquatic environments, and their pectoral fins have undergone remarkable physiological transformations throughout their evolution. Studying early teleost fishes, such as those belonging to the Osteoglossiformes order, could offer crucial insights into the adaptive evolution of pectoral fins within this group. In this study, we have assembled a chromosomal-level genome for the Clown featherback (Chitala ornata), achieving the highest quality genome assembly for Osteoglossiformes to date, with a contig N50 of 32.78 Mb and a scaffold N50 of 40.73 Mb. By combining phylogenetic analysis, we determined that the Clown featherback diverged approximately 202 to 203 million years ago (Ma), aligning with continental separation events. Our analysis revealed the intriguing discovery that a unique deletion of regulatory elements is adjacent to the Gli3 gene, specifically in teleosts. This deletion might be tied to the specialized adaptation of their pectoral fins. Furthermore, our findings indicate that specific contractions and expansions of transposable elements (TEs) in teleosts, including the Clown featherback, could be connected to their adaptive evolution. In essence, this study not only provides a high-quality genomic resource for Osteoglossiformes but also sheds light on the evolutionary trajectory of early teleosts.

The Late Cretaceous eutherian Zalambdalestes reveals unique axis and complex evolution of the mammalian neck

The typical mammalian neck consisting of seven cervical vertebrae (C1-C7) was established by the Late Permian in the cynodont forerunners of modern mammals. This structure is precisely adapted to facilitate movements of the head during feeding, locomotion, predator evasion, and social interactions. Eutheria, the clade including crown placentals, has a fossil record extending back more than 125 million years revealing significant morphological diversification in the Mesozoic. Yet very little is known concerning the early evolution of eutherian cervical morphology and its functional adaptations. A specimen of Zalambdalestes lechei from the Late Cretaceous of Mongolia boasts exceptional preservation of an almost complete series of cervical vertebrae (C2-C7) revealing a highly modified axis (C2). The significance of this cervical morphology is explored utilizing an integrated approach combining comparative anatomical examination across mammals, muscle reconstruction, geometric morphometrics and virtual range of motion analysis. We compared the shape of the axis in Zalambdalestes to a dataset of 88 mammalian species (monotremes, marsupials, and placentals) using three-dimensional landmark analysis. The results indicate that the unique axis morphology of Zalambdalestes has no close analog among living mammals. Virtual range of motion analysis of the neck strongly implies Zalambdalestes was capable of exerting very forceful head movements and had a high degree of ventral flexion for an animal its size. These findings reveal unexpected complexity in the early evolution of the eutherian cervical morphology and suggest a feeding behavior similar to insectivores specialized in vermivory and defensive behaviors in Zalambdalestes akin to modern spiniferous mammals.

Accurate, scalable, and fully automated inference of species trees from raw genome assemblies using ROADIES

Inference of species trees plays a crucial role in advancing our understanding of evolutionary relationships and has immense significance for diverse biological and medical applications. Extensive genome sequencing efforts are currently in progress across a broad spectrum of life forms, holding the potential to unravel the intricate branching patterns within the tree of life. However, estimating species trees starting from raw genome sequences is quite challenging, and the current cutting-edge methodologies require a series of error-prone steps that are neither entirely automated nor standardized. In this paper, we present ROADIES, a novel pipeline for species tree inference from raw genome assemblies that is fully automated, easy to use, scalable, free from reference bias, and provides flexibility to adjust the tradeoff between accuracy and runtime. The ROADIES pipeline eliminates the need to align whole genomes, choose a single reference species, or pre-select loci such as functional genes found using cumbersome annotation steps. Moreover, it leverages recent advances in phylogenetic inference to allow multi-copy genes, eliminating the need to detect orthology. Using the genomic datasets released from large-scale sequencing consortia across three diverse life forms (placental mammals, pomace flies, and birds), we show that ROADIES infers species trees that are comparable in quality with the state-of-the-art approaches but in a fraction of the time. By incorporating optimal approaches and automating all steps from assembled genomes to species and gene trees, ROADIES is poised to improve the accuracy, scalability, and reproducibility of phylogenomic analyses.

Imbalanced speciation pulses sustain the radiation of mammals

The evolutionary histories of major clades, including mammals, often comprise changes in their diversification dynamics, but how these changes occur remains debated. We combined comprehensive phylogenetic and fossil information in a new "birth-death diffusion" model that provides a detailed characterization of variation in diversification rates in mammals. We found an early rising and sustained diversification scenario, wherein speciation rates increased before and during the Cretaceous-Paleogene (K-Pg) boundary. The K-Pg mass extinction event filtered out more slowly speciating lineages and was followed by a subsequent slowing in speciation rates rather than rebounds. These dynamics arose from an imbalanced speciation process, with separate lineages giving rise to many, less speciation-prone descendants. Diversity seems to have been brought about by these isolated, fast-speciating lineages, rather than by a few punctuated innovations.

A mammalian tripartite enhancer cluster controls hypothalamic Pomc expression, food intake, and body weight

Food intake and energy balance are tightly regulated by a group of hypothalamic arcuate neurons expressing the proopiomelanocortin (POMC) gene. In mammals, arcuate-specific POMC expression is driven by two cis-acting transcriptional enhancers known as nPE1 and nPE2. Because mutant mice lacking these two enhancers still showed hypothalamic Pomc mRNA, we searched for additional elements contributing to arcuate Pomc expression. By combining molecular evolution with reporter gene expression in transgenic zebrafish and mice, here, we identified a mammalian arcuate-specific Pomc enhancer that we named nPE3, carrying several binding sites also present in nPE1 and nPE2 for transcription factors known to activate neuronal Pomc expression, such as ISL1, NKX2.1, and ERα. We found that nPE3 originated in the lineage leading to placental mammals and remained under purifying selection in all mammalian orders, although it was lost in Simiiformes (monkeys, apes, and humans) following a unique segmental deletion event. Interestingly, ablation of nPE3 from the mouse genome led to a drastic reduction (>70%) in hypothalamic Pomc mRNA during development and only moderate (<33%) in adult mice. Comparison between double (nPE1 and nPE2) and triple (nPE1, nPE2, and nPE3) enhancer mutants revealed the relative contribution of nPE3 to hypothalamic Pomc expression and its importance in the control of food intake and adiposity in male and female mice. Altogether, these results demonstrate that nPE3 integrates a tripartite cluster of partially redundant enhancers that originated upon a triple convergent evolutionary process in mammals and that is critical for hypothalamic Pomc expression and body weight homeostasis.

Vocal learning-associated convergent evolution in mammalian proteins and regulatory elements

Vocal production learning ("vocal learning") is a convergently evolved trait in vertebrates. To identify brain genomic elements associated with mammalian vocal learning, we integrated genomic, anatomical, and neurophysiological data from the Egyptian fruit bat (Rousettus aegyptiacus) with analyses of the genomes of 215 placental mammals. First, we identified a set of proteins evolving more slowly in vocal learners. Then, we discovered a vocal motor cortical region in the Egyptian fruit bat, an emergent vocal learner, and leveraged that knowledge to identify active cis-regulatory elements in the motor cortex of vocal learners. Machine learning methods applied to motor cortex open chromatin revealed 50 enhancers robustly associated with vocal learning whose activity tended to be lower in vocal learners. Our research implicates convergent losses of motor cortex regulatory elements in mammalian vocal learning evolution.

The evolving world of milk oligosaccharides: Biochemical diversity understood by computational advances

Milk oligosaccharides, complex carbohydrates unique to mammalian milk, play crucial roles in infant nutrition and immune development. This review explores their biochemical diversity, tracing the evolutionary paths that have led to their variation across different species. We highlight the intersection of nutrition, biology, and chemistry in understanding these compounds. Additionally, we discuss the latest computational methods and analytical techniques that have revolutionized the study of milk oligosaccharides, offering insights into their structural complexity and functional roles. This brief but essential review not only aims to provide a deeper understanding of milk oligosaccharides but also discuss the road toward their potential applications.

Genomes, fossils, and the concurrent rise of modern birds and flowering plants in the Late Cretaceous

The phylogeny and divergence timing of the Neoavian radiation remain controversial despite recent progress. We analyzed the genomes of 124 species across all Neoavian orders, using data from 25,460 loci spanning four DNA classes, including 5,756 coding sequences, 12,449 conserved nonexonic elements, 4,871 introns, and 2,384 intergenic segments. We conducted a comprehensive sensitivity analysis to account for the heterogeneity across different DNA classes, leading to an optimal tree of Neoaves with high resolution. This phylogeny features a novel Neoavian dichotomy comprising two monophyletic clades: a previously recognized Telluraves (land birds) and a newly circumscribed Aquaterraves (waterbirds and relatives). Molecular dating analyses with 20 fossil calibrations indicate that the diversification of modern birds began in the Late Cretaceous and underwent a constant and steady radiation across the KPg boundary, concurrent with the rise of angiosperms as well as other major Cenozoic animal groups including placental and multituberculate mammals. The KPg catastrophe had a limited impact on avian evolution compared to the Paleocene-Eocene Thermal Maximum, which triggered a rapid diversification of seabirds. Our findings suggest that the evolution of modern birds followed a slow process of gradualism rather than a rapid process of punctuated equilibrium, with limited interruption by the KPg catastrophe. This study places bird evolution into a new context within vertebrates, with ramifications for the evolution of the Earth's biota.

Rapid adaptation of cellular metabolic rate to the MicroRNA complements of mammals and its relevance to the evolution of endothermy

The metabolic efficiency of mammalian cells depends on the attenuation of intrinsic translation noise by microRNAs. We devised a metric of cellular metabolic rate (cMR), rMR/Mexp optimally fit to the number of microRNA families (mirFam), that is robust to variation in mass and sensitive to body temperature (Tb), consistent with the heat dissipation limit theory of Speakman and Król (2010). Using mirFam as predictor, an Ornstein-Uhlenbeck process of stabilizing selection, with an adaptive shift at the divergence of Boreoeutheria, accounted for 95% of the variation in cMR across mammals. Branchwise rates of evolution of cMR, mirFam and Tb concurrently increased 6- to 7-fold at the divergence of Boreoeutheria, independent of mass. Cellular MR variation across placental mammals was also predicted by the sum of model conserved microRNA-target interactions, revealing an unexpected degree of integration of the microRNA-target apparatus into the energy economy of the mammalian cell.

Molecular Evidence for Relaxed Selection on the Enamel Genes of Toothed Whales (Odontoceti) with Degenerative Enamel Phenotypes

Different species of toothed whales (Odontoceti) exhibit a variety of tooth forms and enamel types. Some odontocetes have highly prismatic enamel with Hunter-Schreger bands, whereas enamel is vestigial or entirely lacking in other species. Different tooth forms and enamel types are associated with alternate feeding strategies that range from biting and grasping prey with teeth in most oceanic and river dolphins to the suction feeding of softer prey items without the use of teeth in many beaked whales. At the molecular level, previous studies have documented inactivating mutations in the enamel-specific genes of some odontocete species that lack complex enamel. At a broader scale, however, it is unclear whether enamel complexity across the full diversity of extant Odontoceti correlates with the relative strength of purifying selection on enamel-specific genes. Here, we employ sequence alignments for seven enamel-specific genes (ACP4, AMBN, AMELX, AMTN, ENAM, KLK4, MMP20) in 62 odontocete species that are representative of all extant families. The sequences for 33 odontocete species were obtained from databases, and sequences for the remaining 29 species were newly generated for this study. We screened these alignments for inactivating mutations (e.g., frameshift indels) and provide a comprehensive catalog of these mutations in species with one or more inactivated enamel genes. Inactivating mutations are rare in Delphinidae (oceanic dolphins) and Platanistidae/Inioidea (river dolphins) that have higher enamel complexity scores. By contrast, mutations are much more numerous in clades such as Monodontidae (narwhal, beluga), Ziphiidae (beaked whales), Physeteroidea (sperm whales), and Phocoenidae (porpoises) that are characterized by simpler enamel or even enamelless teeth. Further, several higher-level taxa (e.g., Hyperoodon, Kogiidae, Monodontidae) possess shared inactivating mutations in one or more enamel genes, which suggests loss of function of these genes in the common ancestor of each clade. We also performed selection (dN/dS) analyses on a concatenation of these genes and used linear regression and Spearman's rank-order correlation to test for correlations between enamel complexity and two different measures of selection intensity (# of inactivating mutations per million years, dN/dS values). Selection analyses revealed that relaxed purifying selection is especially prominent in physeteroids, monodontids, and phocoenids. Linear regressions and correlation analyses revealed a strong negative correlation between selective pressure (dN/dS values) and enamel complexity. Stronger purifying selection (low dN/dS) is found on branches with more complex enamel and weaker purifying selection (higher dN/dS) occurs on branches with less complex enamel or enamelless teeth. As odontocetes diversified into a variety of feeding modes, in particular, the suction capture of prey, a reduced reliance on the dentition for prey capture resulted in the relaxed selection of genes that are critical to enamel development.

Basicranial evidence suggests picrodontid mammals are not stem primates

The Picrodontidae from the middle Palaeocene of North America are enigmatic placental mammals that were allied with various mammalian groups but are generally now considered to have close affinities to paromomyid and palaechthonid plesiadapiforms based on proposed dental synapomorphies. The picrodontid fossil record consists entirely of dental and gnathic remains except for one partial cranium of Zanycteris paleocenus (AMNH 17180). Here, we use µCT technology to unveil previously undocumented morphology in AMNH 17180, describe and compare the basicranial morphology of a picrodontid for the first time, and incorporate these new data into cladistic analyses. The basicranial morphology of Z. paleocenus is distinct from plesiadapiforms and shares similarities with the Palaeogene Apatemyidae and Nyctitheriidae. Results of cladistic analyses incorporating these novel data suggest picrodontids are not stem primates nor euarchontan mammals and that the proposed dental synapomorphies uniting picrodontids with plesiadapiforms and, by extension, primates evolved independently. Results highlight the need to scrutinize proposed synapomorphies of highly autapomorphic taxa with limited fossil records.

Murine trophoblast organoids as a model for trophoblast development and CRISPR-Cas9 screening

The placenta becomes one of the most diversified organs during placental mammal radiation. The main in vitro model for studying mouse trophoblast development is the 2D differentiation model of trophoblast stem cells, which is highly skewed to certain lineages and thus hampers systematic screens. Here, we established culture conditions for the establishment, maintenance, and differentiation of murine trophoblast organoids. Murine trophoblast organoids under the maintenance condition contain stem cell-like populations, whereas differentiated organoids possess various trophoblasts resembling placental ones in vivo. Ablation of Nubpl or Gcm1 in trophoblast organoids recapitulated their deficiency phenotypes in vivo, suggesting that those organoids are valid in vitro models for trophoblast development. Importantly, we performed an efficient CRISPR-Cas9 screening in mouse trophoblast organoids using a focused sgRNA (single guide RNA) library targeting G protein-coupled receptors. Together, our results establish an organoid model to investigate mouse trophoblast development and a practicable approach to performing forward screening in trophoblast lineages.