Phylogenomic analyses provide insights into primate evolution

Proximate mechanisms underlying the coevolution of diet quality and relative brain size in primates

Multiple primate species, including humans, have evolved brains that are surprisingly large relative to their body sizes. Studies of this variation have focused on either proximate (how) or ultimate (why) explanations by correlating species-average brain sizes with, e.g. the rate of genetic changes or certain socioecological variables, respectively. Here, we combined proximate and ultimate perspectives to identify genes that modulated the coevolutionary relationship between diet quality and relative brain size in primates. For n = 50 species, we estimated selection pressure (i.e. root-to-tip dN/dS) for approximately 8K genes and collected brain size, body size and diet quality data. We first used this novel dataset to build on previous studies and bolster findings that neurogenesis-related genes facilitate evolutionary changes in brain size. We then applied phylogenetic partial correlation analysis (to identify genes correlated with both brain size and diet quality) and phylogenetic path analysis (to compare different causal models). We found dozens of genes that may have facilitated the coevolution of diet quality and brain size in primates and show that these genes are involved in neurodevelopment and energy metabolism. This is likely to reflect that higher-quality diets provide more energy to grow and maintain metabolically expensive brains. Our novel approach provides new insight into the drivers of primate brain size evolution.This article is part of the Theo Murphy meeting issue 'Selection shapes diverse animal minds'.

Chromosome fusions shaped karyotype evolution and evolutionary relationships in the model family Brassicaceae

The ancestral crucifer karyotype and 22 conserved genomic blocks (CGBs) facilitate phylogenomic analyses in the Brassicaceae. Chromosomal rearrangements reshuffled CGBs of ancestral chromosomes during karyotype evolution. Here, we identify eight protochromosomes representing the common ancestral karyotype (ACBK) of the two Brassicoideae supertribes: Camelinodae (Lineage I) and Brassicodae (Lineage II). The characterization of multiple cascading fusion events allows us to infer evolutionary relationships based on these events. In the Camelinodae, the ACBK first evolved into the AKI genome, which remained conserved in the Cardamineae, whereas it was altered to tAKI by a reciprocal translocation that preceded the diversification of most Camelinodae tribes. The identified fusion breakpoints largely overlap with CGB boundaries, suggesting that CGBs are mainly disrupted by chromosome fusions. Our results demonstrate the stable inheritance of chromosome fusions and their importance for reconstructing evolutionary relationships. The chromosomal breakpoint approach provides a basis for ancestral state reconstruction based on chromosome-level genome assemblies.

Unveiling cell-type-specific mode of evolution in comparative single-cell expression data

While methodology for determining the mode of evolution in coding sequences has been well established, evaluation of adaptation events in emerging types of phenotype data needs further development. Here we propose an analysis framework (expression variance decomposition, EVaDe) for comparative single-cell expression data based on phenotypic evolution theory. After decomposing the gene expression variance into separate components, we use two strategies to identify genes exhibiting large between-taxon expression divergence and small within-cell-type expression noise in certain cell types, attributing this pattern to putative adaptive evolution. In a dataset of primate prefrontal cortex, we find that such human-specific key genes enrich with neurodevelopment-related functions, while most other genes exhibit neutral evolution patterns. Specific neuron types are found to harbor more of these key genes than other cell types, thus likely to have experienced more extensive adaptation. Reassuringly, at molecular sequence level, the key genes are significantly associated with the rapidly evolving conserved non-coding elements. An additional case analysis comparing the naked mole-rat (NMR) with the mouse suggests that innate-immunity-related genes and cell types have undergone putative expression adaptation in NMR. Overall, the EVaDe framework may effectively probe adaptive evolution mode in single-cell expression data.

Restoring flowcell type and basecaller configuration from FASTQ files of nanopore sequencing data

As nanopore sequencing has been widely adopted, data accumulation has surged, resulting in over 700,000 public datasets. While these data hold immense potential for advancing genomic research, their utility is compromised by the absence of flowcell type and basecaller configuration in about 85% of the data and associated publications. These parameters are essential for many analysis algorithms, and their misapplication can lead to significant drops in performance. To address this issue, we present LongBow, designed to infer flowcell type and basecaller configuration directly from the base quality value patterns of FASTQ files. LongBow has been tested on 66 in-house basecalled FAST5/POD5 datasets and 1989 public FASTQ datasets, achieving accuracies of 95.33% and 91.45%, respectively. We demonstrate its utility by reanalyzing nanopore sequencing data from the COVID-19 Genomics UK (COG-UK) project. The results show that LongBow is essential for reproducing reported genomic variants and, through a LongBow-based analysis pipeline, we discovered substantially more functionally important variants while improving accuracy in lineage assignment. Overall, LongBow is poised to play a critical role in maximizing the utility of public nanopore sequencing data, while significantly enhancing the reproducibility of related research.

Mitogenomes of museum specimens provide new insight into species classification and recently reduced diversity of highly endangered Nomascus gibbons

Our findings reveal that the western black crested gibbon (Nomascus concolor) did not divide into different subspecies, and the relatively low level of genetic diversity emphasizes the importance of monitoring this indicator for vulnerable wildlife. Meanwhile, phylogeographic analysis of the Nomascus genus shows a north-to-south trend of ancestral geographic distribution.

Complete sequencing of ape genomes

The most dynamic and repetitive regions of great ape genomes have traditionally been excluded from comparative studies1-3. Consequently, our understanding of the evolution of our species is incomplete. Here we present haplotype-resolved reference genomes and comparative analyses of six ape species: chimpanzee, bonobo, gorilla, Bornean orangutan, Sumatran orangutan and siamang. We achieve chromosome-level contiguity with substantial sequence accuracy (<1 error in 2.7 megabases) and completely sequence 215 gapless chromosomes telomere-to-telomere. We resolve challenging regions, such as the major histocompatibility complex and immunoglobulin loci, to provide in-depth evolutionary insights. Comparative analyses enabled investigations of the evolution and diversity of regions previously uncharacterized or incompletely studied without bias from mapping to the human reference genome. Such regions include newly minted gene families in lineage-specific segmental duplications, centromeric DNA, acrocentric chromosomes and subterminal heterochromatin. This resource serves as a comprehensive baseline for future evolutionary studies of humans and our closest living ape relatives.

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.

Host ZAP activity correlates with the levels of CpG suppression in primate lentiviruses

Zinc-finger antiviral protein (ZAP) is thought to drive the suppression of CpG dinucleotides in many viruses to mimic the composition of their host genomes. However, in vivo evidence is sparse. Here, we investigated the reasons for unusually high CpG levels in SIVmus and SIVmon from mustached and mona monkeys, descendants of one of the precursors of HIV-1. We show that SIVmus is not resistant to ZAP inhibition. Instead, these Cercopithecus monkey hosts differ from other primate species by a splice site mutation and express the poorly active extralarge XL rather than the highly active L isoform of ZAP. Similarly, higher CpG levels in endogenous prosimian lentiviruses were associated with low activity of the corresponding host lemur ZAPs. In addition, lemur genes also show lower CpG suppression than other primates. Thus, the antiviral activity of ZAP not only affects suppression of CpG dinucleotides in viral transcripts but possibly also host genomes.

Genetic Diversity, Genetic Structure, and Demographic History of Black Snub-Nosed Monkey (Rhinopithecus strykeri) in the Gaoligong Mountains, Southwestern China

The Gaoligong Mountains, located in the southeastern Tibetan Plateau, is one of the world's biodiversity hotspots and provides a refugium for many endangered endemic animals. In this study, we reported a population genetic study on black snub-nosed monkey (Rhinopithecus strykeri), a critically endangered primate endemic to the Gaoligong Mountains, yet their large-scale population genetic study remains underexplored. Here, we performed population genetic analyses from two geographical populations (Pianma and Luoma) based on targeted genomic single-nucleotide polymorphism (SNP) data (37.7 K) and mitochondrial DNA (mtDNA) control region (842 bp). Both nuclear SNP data and mtDNA revealed relatively low levels of genetic variation in both populations compared to other reported primates, which is most likely to be explained by loss of historical genetic diversity due to inbreeding and long-term small effective population size, thus potentially aggravating the effects of inbreeding and genetic depression. Phylogenetic and population structure analyses for mtDNA revealed two deep lineages (approximately 0.69 million years ago), but limited genetic differentiation in nuclear data, which might have been caused by the Late Cenozoic uplift of the Tibetan Plateau and glacial refuge, and subsequent secondary contact as a result of historically high and bidirectional gene flow between populations. Ecological niche modeling and landscape connectivity analyses also showed historical and recent connectivity between two geographical populations. The demographic history inferred from both mtDNA and nuclear data revealed at least two continuous declines in the effective population size occurring around 43 Kya and 8-10 Kya, respectively, probably due to Pleistocene glaciations and subsequent human activities. Our results provide the first detailed and comprehensive genetic insights into the genetic diversity, population structure, and demographic history of a critically endangered species, and provide essential baseline information to guide conservation efforts.

Genomic and phenotypic evidence support visual and olfactory shifts in primate evolution

Sensory trade-offs between vision and olfaction in the evolution and radiation of primates have long been debated. However, insights have been limited by a lack of sensory gene sequences and accompanying functional predictions. Here we conduct large-scale functional analyses of visual and olfactory receptors and related brain regions across extant primates. Our results reveal a visual shift from ultraviolet to violet colour sensitivity in early haplorrhine primates, followed by acceleration in the rhodopsin retinal release rates at the origin of anthropoids, both of which are expected to greatly enhance visual acuity under brighter light conditions. Additionally, we find that the sensitivity of olfactory receptors shifted from narrowly to broadly tuned early in anthropoid evolution. In contrast, strepsirrhines appear to have retained sensitive dim-light vision and underwent functional enhancement of narrowly tuned olfactory receptors. Our models indicate that this would have enhanced odorant discrimination and facilitated olfaction-mediated physiology and behaviour. These differences in tuning patterns of olfactory receptors between major primate lineages mirror well-established morphological differences in external anatomy and brain structures, revealing new mechanisms of olfactory adaptation and evolutionary plasticity. Our multisystem analyses reveal patterns of co-evolution in genomic, molecular and neuroanatomical traits that are consistent with a sensory 'reallocation' rather than strict trade-offs.

An atlas of early human mandibular endochondral and osteogenic paracrine signaling regions of Meckel's cartilage

The mandible, also known as the lower jaw, is the only bone in the skull that can move and is essential for speaking and chewing. Meckel's cartilage (MC) is a temporary structure that supports the formation of the mandible, but how MC is involved in the ossification of the mandible is poorly understood. Through the use of single-cell RNA sequencing and single-cell spatial transcriptomics analyses, a spatiotemporal atlas of MC in human fetuses from 7 to 15 wk postconception was established, highlighting the role of MC in the ossification of the mandible. Importantly, we revealed that two populations of MC contributed to mandibular ossification through different mechanisms. The anterior MC can differentiate into osteolineage cells, as shown in an in vivo lineage tracing mouse model. The intermediate MC facilitates intramembranous ossification through cell-cell communications, possibly through signaling ligands like BMP5, BMP7, SEMA3A, PDGFC, and FGF7. This study suggests that MC plays a crucial role in mediating mandibular ossification through distinct mechanisms, providing valuable insights for understanding oral and craniofacial diseases and disorders in the future.

Generation and characterization of induced pluripotent stem cells of small apes

Small apes (family Hylobatidae), encompassing gibbons and siamangs, occupy a pivotal evolutionary position within the hominoid lineage, bridging the gap between great apes and catarrhine monkeys. Although they possess distinctive genomic and phenotypic features-such as rapid chromosomal rearrangements and adaptations for brachiation-functional genomic studies on small apes have been hindered by the limited availability of biological samples and developmental models. Here, we address this gap by successfully reprogramming primary skin fibroblasts from three small ape species: lar gibbons (Hylobates lar), Abbott's gray gibbons (Hylobates abbotti), and siamangs (Symphalangus syndactylus). Using Sendai virus-based stealth RNA vectors, we generated 31 reprogrammed cell lines, five of which were developed into transgene-free induced pluripotent stem cells. These iPSCs displayed canonical features of primed pluripotency, both morphologically and molecularly, consistent with other primate iPSCs. Directed differentiation experiments confirmed the capacity of the small ape iPSCs to generate cells representing all three germ layers. In particular, their successful differentiation into limb bud mesoderm cells underscores their utility in investigating the molecular and developmental mechanisms unique to small ape forelimb evolution. Transcriptomic profiling of small ape iPSCs revealed significant upregulation of pluripotency-associated genes, alongside elevated expression of transposable elements. Remarkably, LAVA retrotransposons-a class of elements specific to small apes-exhibited particularly high expression levels in these cells. Comparative transcriptomic analyses with iPSCs from humans, great apes, and macaques identified evolutionary trends and clade-specific gene expression signatures. These signatures highlighted processes linked to genomic stability and cell death, providing insights into small ape-specific adaptations. This study positions small ape iPSCs as a transformative tool for advancing functional genomics and evolutionary developmental biology. By facilitating detailed investigations into hominoid genome evolution and phenotypic diversification, this system bridges critical gaps in comparative research, enabling deeper exploration of the genetic and cellular underpinnings of small ape-specific traits.

The Principle of Cortical Development and Evolution

Human's robust cognitive abilities, including creativity and language, are made possible, at least in large part, by evolutionary changes made to the cerebral cortex. This paper reviews the biology and evolution of mammalian cortical radial glial cells (primary neural stem cells) and introduces the concept that a genetically step wise process, based on a core molecular pathway already in use, is the evolutionary process that has molded cortical neurogenesis. The core mechanism, which has been identified in our recent studies, is the extracellular signal-regulated kinase (ERK)-bone morphogenic protein 7 (BMP7)-GLI3 repressor form (GLI3R)-sonic hedgehog (SHH) positive feedback loop. Additionally, I propose that the molecular basis for cortical evolutionary dwarfism, exemplified by the lissencephalic mouse which originated from a larger gyrencephalic ancestor, is an increase in SHH signaling in radial glia, that antagonizes ERK-BMP7 signaling. Finally, I propose that: (1) SHH signaling is not a key regulator of primate cortical expansion and folding; (2) human cortical radial glial cells do not generate neocortical interneurons; (3) human-specific genes may not be essential for most cortical expansion. I hope this review assists colleagues in the field, guiding research to address gaps in our understanding of cortical development and evolution.

Mitochondrial DNA for Phylogeny Building: Assessing Individual and Grouped mtGenes as Proxies for the mtGenome in Platyrrhines

Phylogenetic trees are analytic tools used in primate studies to elucidate evolutionary relationships. Because of its relative ease to sequence and rapid evolution compared to nuclear genomes, mitochondrial DNA is frequently used for phylogeny building. This project evaluated the effectiveness of using individual or grouped mitochondrial genes (mtGenes) as a proxy for the mitochondrial genome (mtGenome) in phylogeny building within two nested primate datasets, Cebidae and Platyrrhini, with differing divergence dates. mtGene utility rankings were determined based on congruence values to the mtGenome tree. mtGenes trees were also assessed on tree resolution and ability to sort nested clades. We found that most individual mtGenes, including ribosomal genes (12S and 16S), COX genes, most ND genes, and d-Loop are not appropriate for use as proxies for the mtGenome when tree building in either the Cebidae or Platyrrhini set. On average, grouped mtGenes outperformed individual mtGenes in both sets, and mtGene and grouped mtGene rankings varied between sets. Pairing CYB and COX3 together or pairing ND2 and CYB worked well in both the Cebidae set and the Platyrrhini set. We also found that nucleotide diversity is not a predictor of mtGene performance. Instead, it may be that unique mtGene or mtGene system evolutionary history impacts mtGene performance.

Functional, Immunogenetic, and Structural Convergence in Influenza Immunity between Humans and Macaques

Human B cell immunity to the influenza hemagglutinin (HA) stem region, a universal influenza vaccine target, is often stereotyped and immunogenetically restricted, posing challenges for study outside humans. Here, we show that macaques vaccinated with a HA stem immunogen elicit human-like public B cell lineages targeting two major conserved sites of vulnerability, the central stem and anchor epitopes. Central stem antibodies were predominantly derived from V H 1-138, the macaque homolog of human V H 1-69, a V H -gene preferentially used in human central stem broadly neutralizing antibodies (bnAbs). Similarly, macaques produced anchor bnAbs with the human-like NWP motif. Both bnAb lineages were functionally and structurally analogous to their human counterparts, with recognition mediated largely by germline-encoded motifs. Thus the macaque immunoglobulin repertoire supports human-like public bnAb responses to influenza HA. Moreover, this underscores the utility of homologous germline-encoded immunity, suggesting that immune repertoires of macaques and humans may have been similarly shaped during evolution.

HIGHLIGHTS

Functional human-like public antibody lineages can be elicited to HA stem supersites in macaques. Macaque central stem bnAbs are predominantly derived from V H 1-138, a V H -gene homologous to human V H 1-69. The human-like CDR L3 NWP anchor epitope-targeting lineage can be elicited in macaques.Central stem and anchor bnAbs from humans and macaques engage their respective epitopes with atomic level similarity.

Long-read sequencing and genome assembly of natural history collection samples and challenging specimens

Museum collections harbor millions of samples, largely unutilized for long-read sequencing. Here, we use ethanol-preserved samples containing kilobase-sized DNA to show that amplification-free protocols can yield contiguous genome assemblies. Additionally, using a modified amplification-based protocol, employing an alternative polymerase to overcome PCR bias, we assemble the 3.1 Gb maned sloth genome, surpassing the previous 500 Mb protocol size limit. Our protocol also improves assemblies of other difficult-to-sequence molluscs and arthropods, including millimeter-sized organisms. By highlighting collections as valuable sample resources and facilitating genome assembly of tiny and challenging organisms, our study advances efforts to obtain reference genomes of all eukaryotes.

Phylogenetic Signal in Primate Tooth Enamel Proteins and its Relevance for Paleoproteomics

Ancient tooth enamel, and to some extent dentin and bone, contain characteristic peptides that persist for long periods of time. In particular, peptides from the enamel proteome (enamelome) have been used to reconstruct the phylogenetic relationships of fossil taxa. However, the enamelome is based on only about 10 genes, whose protein products undergo fragmentation in vivo and post mortem. This raises the question as to whether the enamelome alone provides enough information for reliable phylogenetic inference. We address these considerations on a selection of enamel-associated proteins that has been computationally predicted from genomic data from 232 primate species. We created multiple sequence alignments for each protein and estimated the evolutionary rate for each site. We examined which sites overlap with the parts of the protein sequences that are typically isolated from fossils. Based on this, we simulated ancient data with different degrees of sequence fragmentation, followed by phylogenetic analysis. We compared these trees to a reference species tree. Up to a degree of fragmentation that is similar to that of fossil samples from 1 to 2 million years ago, the phylogenetic placements of most nodes at family level are consistent with the reference species tree. We tested phylogenetic analysis on combinations of different enamel proteins and found that the composition of the proteome can influence deep splits in the phylogeny. With our methods, we provide guidance for researchers on how to evaluate the potential of paleoproteomics for phylogenetic studies before sampling valuable ancient specimens.

Comparative single-cell multiome identifies evolutionary changes in neural progenitor cells during primate brain development

Understanding the cellular and genetic mechanisms driving human-specific features of cortical development remains a challenge. We generated a cell-type resolved atlas of transcriptome and chromatin accessibility in the developing macaque and mouse prefrontal cortex (PFC). Comparing with published human data, our findings demonstrate that although the cortex cellular composition is overall conserved across species, progenitor cells show significant evolutionary divergence in cellular properties. Specifically, human neural progenitors exhibit extensive transcriptional rewiring in growth factor and extracellular matrix (ECM) pathways. Expression of the human-specific progenitor marker ITGA2 in the fetal mouse cortex increases the progenitor proliferation and the proportion of upper-layer neurons. These transcriptional divergences are primarily driven by altered activity in the distal regulatory elements. The chromatin regions with human-gained accessibility are enriched with human-specific sequence changes and polymorphisms linked to intelligence and neuropsychiatric disorders. Our results identify evolutionary changes in neural progenitors and putative gene regulatory mechanisms shaping primate brain evolution.

Perspectives and opportunities in forensic human, animal, and plant integrative genomics in the Pangenome era

The Human Pangenome Reference Consortium, the Chinese Pangenome Consortium, and other plant and animal pangenome projects have announced the completion of pilot work aimed at constructing high-quality, haplotype-resolved reference graph genomes representative of global ethno-linguistically different populations or different plant and animal species. These graph-based, gapless pangenome references, which are enriched in terms of genomic diversity, completeness, and contiguity, have the potential for enhancing long-read sequencing (LRS)-based genomic research, as well as improving mappability and variant genotyping on traditional short-read sequencing platforms. We comprehensively discuss the advancements in pangenome-based genomic integrative genomic discoveries across forensic-related species (humans, animals, and plants) and summarize their applications in variant identification and forensic genomics, epigenetics, transcriptomics, and microbiome research. Recent developments in multiplexed array sequencing have introduced a highly efficient and programmable technique to overcome the limitations of short forensic marker lengths in LRS platforms. This technique enables the concatenation of short RNA transcripts and DNA fragments into LRS-optimal molecules for sequencing, assembly, and genotyping. The integration of new pangenome reference coordinates and corresponding computational algorithms will benefit forensic integrative genomics by facilitating new marker identification, accurate genotyping, high-resolution panel development, and the updating of statistical algorithms. This review highlights the necessity of integrating LRS-based platforms, pangenome-based study designs, and graph-based pangenome references in short-read mapping and LRS-based innovations to achieve precision forensic science.

Mutation Rate Variation and Other Challenges in 2-LTR Dating of Primate Endogenous Retrovirus Integrations

The time of integration of germline-targeting Long Terminal Repeat (LTR) retroposons, such as endogenous retroviruses (ERVs), can be estimated by assessing the nucleotide divergence between the LTR sequences flanking the viral genes. Due to the viral replication mechanism, both LTRs are identical at the moment of integration, when the provirus becomes part of the host genome. After that time, proviral sequences evolve within the host DNA. When the mutation rate is known, nucleotide divergence between the LTRs would then be a measure of time elapsed since integration. Though frequently used, the approach has been complicated by the choice of host mutation rate and, to a lesser extent, by the method selected to estimate nucleotide divergence. As a result, outcomes can be incompatible with, for instance, speciation events identified from the fossil record. The review will give an overview of research reporting LTR-retroposon dating, and a summary of important factors to consider, including the quality, assembly, and alignment of sequences, the mutation rate of foreign DNA in host genomes, and the choice of a distance estimation method. Primates will here be the focus of the analysis because their genomes, ERVs, and fossil record have been extensively studied. However, most of the factors discussed have a wide applicability in the vertebrate field.

Unprecedented female mutation bias in the aye-aye, a highly unusual lemur from Madagascar

Every mammal studied to date has been found to have a male mutation bias: male parents transmit more de novo mutations to offspring than female parents, contributing increasingly more mutations with age. Although male-biased mutation has been studied for more than 75 years, its causes are still debated. One obstacle to understanding this pattern is its near universality-without variation in mutation bias, it is difficult to find an underlying cause. Here, we present new data on multiple pedigrees from two primate species: aye-ayes (Daubentonia madagascariensis), a member of the strepsirrhine primates, and olive baboons (Papio anubis). In stark contrast to the pattern found across mammals, we find a much larger effect of maternal age than paternal age on mutation rates in the aye-aye. In addition, older aye-aye mothers transmit substantially more mutations than older fathers. We carry out both computational and experimental validation of our results, contrasting them with results from baboons and other primates using the same methodologies. Further, we analyze a set of DNA repair and replication genes to identify candidate mutations that may be responsible for the change in mutation bias observed in aye-ayes. Our results demonstrate that mutation bias is not an immutable trait, but rather one that can evolve between closely related species. Further work on aye-ayes (and possibly other lemuriform primates) should help to explain the molecular basis for sex-biased mutation.

Alu-Sc-mediated exonization generated a mitochondrial LKB1 gene variant found only in higher order primates

The tumor suppressor LKB1/STK11 plays important roles in regulating cellular metabolism and stress responses and its mutations are associated with various cancers. We recently identified a novel exon 1b within intron 1 of human LKB1/STK11, which generates an alternatively spliced, mitochondria-targeting LKB1 isoform important for regulating mitochondrial oxidative stress. Here we examined the formation of this novel exon 1b and uncovered its relatively late emergence during evolution. Analyses of putative exon 1b genomic sequences within the primate superfamily indicated that the exonization of LKB1/STK11 exon 1b was mediated by the conserved retrotransposable element Alu-Sc. While putative exon 1b sequences are recognizable in most members of the primate family from New World Monkeys onwards, characteristically functional LKB1/STK11 exon 1b, with translation start and 5' and 3' splice sites, could only be found in greater apes and human, and interestingly, correlates with their increased body mass and longevity development.

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.

Genomic insights into ecological adaptation of oaks revealed by phylogenomic analysis of multiple species

Understanding the ecological adaptation of tree species can not only reveal the evolutionary potential but also benefit biodiversity conservation under global climate change. Quercus is a keystone genus in Northern Hemisphere forests, and its wide distribution in diverse ecosystems and long evolutionary history make it an ideal model for studying the genomic basis of ecological adaptations. Here we used a newly sequenced genome of Quercus gilva, an evergreen oak species from East Asia, with 18 published Fagales genomes to determine how Fagaceae genomes have evolved, identify genomic footprints of ecological adaptability in oaks in general, as well as between evergreen and deciduous oaks. We found that oak species exhibited a higher degree of genomic conservation and stability, as indicated by the absence of large-scale chromosomal structural variations or additional whole-genome duplication events. In addition, we identified expansion and tandem repetitions within gene families that contribute to plant physical and chemical defense (e.g., cuticle biosynthesis and oxidosqualene cyclase genes), which may represent the foundation for the ecological adaptation of oak species. Circadian rhythm and hormone-related genes may regulate the habits of evergreen and deciduous oaks. This study provides a comprehensive perspective on the ecological adaptations of tree species based on phylogenetic, genome evolutionary, and functional genomic analyses.

Ecological and anthropogenic effects on the genomic diversity of lemurs in Madagascar

Ecological variation and anthropogenic landscape modification have had key roles in the diversification and extinction of mammals in Madagascar. Lemurs represent a radiation with more than 100 species, constituting roughly one-fifth of the primate order. Almost all species of lemurs are threatened with extinction, but little is known about their genetic diversity and demographic history. Here, we analyse high-coverage genome-wide resequencing data from 162 unique individuals comprising 50 species of Lemuriformes, including multiple individuals from most species. Genomic diversity varies widely across the infraorder and yet is broadly consistent among individuals within species. We show widespread introgression in multiple genera and generally high levels of genomic diversity likely resulting from allele sharing that occurred during periods of connectivity and fragmentation during climatic shifts. We find distinct patterns of demographic history in lemurs across the ecogeographic regions of Madagascar within the last million years. Within the past 2,000 years, lemurs underwent major declines in effective population size that corresponded to the timing of human population expansion in Madagascar. In multiple regions of the island, we identified chronological trajectories of inbreeding that are consistent across genera and species, suggesting localized effects of human activity. Our results show how the extraordinary diversity of these long-neglected, endangered primates has been influenced by ecological and anthropogenic factors.

Cis-Regulatory Evolution of CCNB1IP1 Driving Gradual Increase of Cortical Size and Folding in primates

Neocortex expansion has a concerted relationship with folding, underlying evolution of human cognitive functions. However, molecular mechanisms underlying this significant evolutionary process remains unknown. Here, using tree shrew as an outgroup of primates, we identify a new regulator CCNB1IP1, which acquired its expression before the emergence of primates. Following the evolution of cis-regulatory elements, the CCNB1IP1 expression has steadily increased over the course of primate brain evolution, mirroring the gradual increase of neocortex. Mechanistically, we elucidated that CCNB1IP1 expression can cause an increase in neural progenitors through shortening G1 phase. Consistently, the CCNB1IP1 knock-in mouse model exhibited traits associated with enhanced learning and memory abilities. Together, our study reveals how changes in CCNB1IP1 expression may have contributed to the gradual evolution in primate brain.

If you please, draw me a neuron - linking evolutionary tinkering with human neuron evolution

Animal speciation often involves novel behavioral features that rely on nervous system evolution. Human-specific brain features have been proposed to underlie specialized cognitive functions and to be linked, at least in part, to the evolution of synapses, neurons, and circuits of the cerebral cortex. Here, we review recent results showing that, while the human cortex is composed of a repertoire of cells that appears to be largely similar to the one found in other mammals, human cortical neurons do display specialized features at many levels, from gene expression to intrinsic physiological properties. The molecular mechanisms underlying human species-specific neuronal features remain largely unknown but implicate hominid-specific gene duplicates that encode novel molecular modifiers of neuronal function. The identification of human-specific genetic modifiers of neuronal function brings novel insights on brain evolution and function and, could also provide new insights on human species-specific vulnerabilities to brain disorders.

Reconstructing human-specific regulatory functions in model systems

Uniquely human physical traits, such as an expanded cerebral cortex and changes in limb morphology that allow us to use tools and walk upright, are in part due to human-specific genetic changes that altered when, where, and how genes are expressed during development. Over 20 000 putative regulatory elements with potential human-specific functions have been discovered. Understanding how these elements contributed to human evolution requires identifying candidates most likely to have shaped human traits, then studying them in genetically modified animal models. Here, we review the progress and challenges in generating and studying such models and propose a pathway for advancing the field. Finally, we highlight that large-scale collaborations across multiple research domains are essential to decipher what makes us human.

Haplotype-phased genome assemblies and annotation of the northern white-cheeked gibbon (Nomascus leucogenys)

Nomascus leucogenys is a critically endangered species of small apes. Here, we sequenced and assembled the male genome of N. leucogenys, using PacBio and Hi-C datasets, with a particular focus on its Y-chromosome. The resulting high-quality haplotype-phased assemblies are at chromosome-scale, with scaffold/contig N50 values of 124.2/102.2 Mb for Haplotype 1 and 121.2/85.67 Mb for Haplotype 2. The assembled Y-chromosome spans 16.06 Mb. BUSCO assessment indicated completeness scores exceeding 95%. We predicted 18,925 protein-coding genes (23,783 mRNAs), including 58 genes on the Y-chromosome. Approximately 50% of the genome comprises repetitive elements. These comprehensive genome datasets will serve as a valuable resource for future studies on the genetics and protection of gibbons and improve our understanding on the evolution of Y-chromosome-related genes in primates.

A whole-genome scan for evidence of recent positive and balancing selection in aye-ayes (Daubentonia madagascariensis) utilizing a well-fit evolutionary baseline model

The aye-aye (Daubentonia madagascariensis) is one of the 25 most endangered primate species in the world, maintaining amongst the lowest genetic diversity of any primate measured to date. Characterizing patterns of genetic variation within aye-aye populations, and the relative influences of neutral and selective processes in shaping that variation, is thus important for future conservation efforts. In this study, we performed the first whole-genome scans for recent positive and balancing selection in the species, utilizing high-coverage population genomic data from newly sequenced individuals. We generated null thresholds for our genomic scans by creating an evolutionarily appropriate baseline model that incorporates the demographic history of this aye-aye population, and identified a small number of candidate genes. Most notably, a suite of genes involved in olfaction - a key trait in these nocturnal primates - were identified as experiencing long-term balancing selection. We also conducted analyses to quantify the expected statistical power to detect positive and balancing selection in this population using site frequency spectrum-based inference methods, once accounting for the potentially confounding contributions of population history, recombination and mutation rate variation, and purifying and background selection. This work, presenting the first high-quality, genome-wide polymorphism data across the functional regions of the aye-aye genome, thus provides important insights into the landscape of episodic selective forces in this highly endangered species.

The evolution of ageing: classic theories and emerging ideas

Ageing is generally regarded as a non-adaptive by-product of evolution. Based on this premise three classic evolutionary theories of ageing have been proposed. These theories have dominated the literature for several decades. Despite their individual nuances, the common thread which unites them is that they posit that ageing results from a decline in the intensity of natural selection with chronological age. Empirical evidence has been identified which supports each theory. However, a consensus remains to be fully established as to which theory best accounts for the evolution of ageing. A consequence of this uncertainty are counter arguments which advocate for alternative theoretical frameworks, such as those which propose an adaptive origin for ageing, senescence, or death. Given this backdrop, this review has several aims. Firstly, to briefly discuss the classic evolutionary theories. Secondly, to evaluate how evolutionary forces beyond a monotonic decrease in natural selection can affect the evolution of ageing. Thirdly, to examine alternatives to the classic theories. Finally, to introduce a pluralistic interpretation of the evolution of ageing. The basis of this pluralistic theoretical framework is the recognition that certain evolutionary ideas will be more appropriate depending on the organism, its ecological context, and its life history.

Environmental Driving of Adaptation Mechanism on Rumen Microorganisms of Sheep Based on Metagenomics and Metabolomics Data Analysis

Natural or artificial selection causes animals to adapt to their environment. The adaptive changes generated by the rumen population and metabolism form the basis of ruminant evolution. In particular, the adaptive drive for environmental adaptation reflects the high-quality traits of sheep that have migrated from other places or have been distant from their origins for a long time. The Hu sheep is the most representative sheep breed in the humid and low-altitude environments (Tai Lake region) in East Asia and has been widely introduced into the arid and high-altitude environments (Tibetan Plateau and Hotan region), resulting in environmental adaptive changes in the Hu sheep. In this study, a joint analysis of the rumen microbial metagenome and metabolome was conducted on Hu sheep from different regions (area of origin and area of introduction) with the objective of investigating the quality traits of Hu sheep and identifying microorganisms that influence the adaptive drive of ruminants. The results demonstrated that the growth performance of Hu sheep was altered due to changes in rumen tissue and metabolism following their introduction to the arid area at relatively high altitude. Metagenomic and metabolomic analyses (five ramsper area) revealed that 3580 different microorganisms and 732 different metabolites were identified in the rumen fluid of arid sheep. Among these, the representative upregulated metabolites were 4,6-isocanedione, methanesulfonic acid and N2-succinyl-L-arginine, while the dominant microorganism was Prevotella ruminicola. The downregulated metabolites were identified as campesterol, teprenone and dihydroclavaminic acid, while the disadvantaged microorganisms were Dialister_succinatiphilus, Prevotella_sp._AGR2160, Prevotella_multisaccharivorax and Selenomonas_bovis. The results of the Pearson analysis indicated that the rumen microbiota and metabolite content of sheep were significantly altered and highly correlated following their relocation from a humid lowland to an arid upland. In particular, the observed changes in rumen microorganisms led to an acceleration of body metabolism, rendering sheep highly adaptable to environmental stress. Prevotella_ruminicola was identified as playing an important role in this process. These findings provide insights into the environmental adaptation mechanisms of sheep.

Chromosome-level genome assembly of the cashmere goat

The goat, an early domesticated ruminant, is a reliable source of cashmere, meat and milk in global agricultural production. Despite this, the genome of cashmere-rich goats has yet to be characterized. Here, we assembled the nearly complete genome of a cashmere goat from a highly economically valuable Inner Mongolian Cashmere buck, utilizing a combination of PacBio HiFi, ONT ultra-long reads, and Hi-C technologies. The size of this genome is 2.76 Gb, with a contig N50 of 95.22 Mb. All assembled sequences were anchored onto 29 autosomes and both sex chromosomes, with only two gaps present on the X chromosome. We identified 1,333.29 Mb (48.26%) of repetitive sequences and predicted 22,480 protein-coding genes. Assembly quality assessment of the genome demonstrated that our assembled cashmere goat genome surpasses the continuity, completeness, and accuracy of other published goat genomes. Taken together, we provided the first cashmere goat assembly, bridging the gap in the genome of important economic breeds of domestic goats, and providing a valuable reference resource for goat genetics and genome research.

Whole genomes of Amazonian uakari monkeys reveal complex connectivity and fast differentiation driven by high environmental dynamism

Despite showing the greatest primate diversity on the planet, genomic studies on Amazonian primates show very little representation in the literature. With 48 geolocalized high coverage whole genomes from wild uakari monkeys, we present the first population-level study on platyrrhines using whole genome data. In a very restricted range of the Amazon rainforest, eight uakari species (Cacajao genus) have been described and categorized into the bald and black uakari groups, based on phenotypic and ecological differences. Despite a slight habitat overlap, we show that posterior to their split 0.92 Mya, bald and black uakaris have remained independent, without gene flow. Nowadays, these two groups present distinct genetic diversity and group-specific variation linked to pathogens. We propose differing hydrology patterns and effectiveness of geographic barriers have modulated the intra-group connectivity and structure of bald and black uakari populations. With this work we have explored the effects of the Amazon rainforest's dynamism on wild primates' genetics and increased the representation of platyrrhine genomes, thus opening the door to future research on the complexity and diversity of primate genomics.

A hybrid genome assembly of the endangered aye-aye (Daubentonia madagascariensis)

The aye-aye (Daubentonia madagascariensis) is the only extant member of the Daubentoniidae primate family. Although several reference genomes exist for this endangered strepsirrhine primate, the predominant usage of short-read sequencing has resulted in limited assembly contiguity and completeness, and no protein-coding gene annotations have yet been released. Here, we present a novel, fully annotated, chromosome-level hybrid de novo assembly for the species based on a combination of Oxford Nanopore Technologies long reads and Illumina short reads and scaffolded using genome-wide chromatin interaction data-a community resource that will improve future conservation efforts as well as primate comparative analyses.

Complete sequencing of ape genomes

We present haplotype-resolved reference genomes and comparative analyses of six ape species, namely: chimpanzee, bonobo, gorilla, Bornean orangutan, Sumatran orangutan, and siamang. We achieve chromosome-level contiguity with unparalleled sequence accuracy (<1 error in 500,000 base pairs), completely sequencing 215 gapless chromosomes telomere-to-telomere. We resolve challenging regions, such as the major histocompatibility complex and immunoglobulin loci, providing more in-depth evolutionary insights. Comparative analyses, including human, allow us to investigate the evolution and diversity of regions previously uncharacterized or incompletely studied without bias from mapping to the human reference. This includes newly minted gene families within lineage-specific segmental duplications, centromeric DNA, acrocentric chromosomes, and subterminal heterochromatin. This resource should serve as a definitive baseline for all future evolutionary studies of humans and our closest living ape relatives.

Deciphering the role of structural variation in human evolution: a functional perspective

Advances in sequencing technologies have enabled the comparison of high-quality genomes of diverse primate species, revealing vast amounts of divergence due to structural variation. Given their large size, structural variants (SVs) can simultaneously alter the function and regulation of multiple genes. Studies estimate that collectively more than 3.5% of the genome is divergent in humans versus other great apes, impacting thousands of genes. Functional genomics and gene-editing tools in various model systems recently emerged as an exciting frontier - investigating the wide-ranging impacts of SVs on molecular, cellular, and systems-level phenotypes. This review examines existing research and identifies future directions to broaden our understanding of the functional roles of SVs on phenotypic innovations and diversity impacting uniquely human features, ranging from cognition to metabolic adaptations.

Chromosome-level genome and population genomics of the intermediate horseshoe bat ( Rhinolophus affinis) reveal the molecular basis of virus tolerance in Rhinolophus and echolocation call frequency variation

Horseshoe bats (genus Rhinolophus, family Rhinolophidae) represent an important group within chiropteran phylogeny due to their distinctive traits, including constant high-frequency echolocation, rapid karyotype evolution, and unique immune system. Advances in evolutionary biology, supported by high-quality reference genomes and comprehensive whole-genome data, have significantly enhanced our understanding of species origins, speciation mechanisms, adaptive evolutionary processes, and phenotypic diversity. However, genomic research and understanding of the evolutionary patterns of Rhinolophus are severely constrained by limited data, with only a single published genome of R. ferrumequinum currently available. In this study, we constructed a high-quality chromosome-level reference genome for the intermediate horseshoe bat ( R. affinis). Comparative genomic analyses revealed potential genetic characteristics associated with virus tolerance in Rhinolophidae. Notably, we observed expansions in several immune-related gene families and identified various genes functionally associated with the SARS-CoV-2 signaling pathway, DNA repair, and apoptosis, which displayed signs of rapid evolution. In addition, we observed an expansion of the major histocompatibility complex class II (MHC-II) region and a higher copy number of the HLA- DQB2 gene in horseshoe bats compared to other chiropteran species. Based on whole-genome resequencing and population genomic analyses, we identified multiple candidate loci (e.g., GLI3) associated with variations in echolocation call frequency across R. affinis subspecies. This research not only expands our understanding of the genetic characteristics of the Rhinolophus genus but also establishes a valuable foundation for future research.

Interaction between the TBC1D24 TLDc domain and the KIBRA C2 domain is disrupted by two epilepsy-associated TBC1D24 missense variants

Mutations of human TBC1D24 are associated with deafness, epilepsy, or DOORS syndrome (deafness, onychodystrophy, osteodystrophy, cognitive disability, and seizures). The causal relationships between TBC1D24 variants and the different clinical phenotypes are not understood. Our hypothesis is that phenotypic heterogeneity of missense mutations of TBC1D24 results, in part, from perturbed binding of different protein partners. To discover novel protein partners of TBC1D24, we conducted yeast two-hybrid (Y2H) screen using mouse full-length TBC1D24 as bait. Kidney and brain protein (KIBRA), a scaffold protein encoded by Wwc1, was identified as a partner of TBC1D24. KIBRA functions in the Hippo signaling pathway and is important for human cognition and memory. The TBC1D24 TLDc domain binds to KIBRA full-length and to its C2 domain, confirmed by Y2H assays. No interaction was detected with Y2H assays between the KIBRA C2 domain and TLDc domains of NCOA7, MEAK7, and OXR1. Moreover, the C2 domains of other WWC family proteins do not interact with the TLDc domain of TBC1D24, demonstrating specificity. The mRNAs encoding TBC1D24 and KIBRA proteins in mouse are coexpressed at least in a subset of hippocampal cells indicating availability to interact in vivo. As two epilepsy-associated recessive variants (Gly511Arg and Ala515Val) in the TLDc domain of human TBC1D24 disrupt the interaction with the human KIBRA C2 domain, this study reveals a pathogenic mechanism of TBC1D24-associated epilepsy, linking the TBC1D24 and KIBRA pathways. The interaction of TBC1D24-KIBRA is physiologically meaningful and necessary to reduce the risk of epilepsy.

Vocal labeling of others by nonhuman primates

Humans, dolphins, and elephants are the only known species that vocally label their conspecifics. It remains unclear whether nonhuman primates share this ability. We recorded spontaneous "phee-call" dialogues between pairs of marmoset monkeys. We discovered that marmosets use these calls to vocally label their conspecifics. Moreover, they respond more consistently and correctly to calls that are specifically directed at them. Analysis of calls from multiple monkeys revealed that family members use similar calls and acoustic features to label others and perform vocal learning. These findings shed light on the complexities of social vocalizations among nonhuman primates and suggest that marmoset vocalizations may provide a model for understanding aspects of human language, thereby offering new insights into the evolution of social communication.

Spatiotemporal omics for biology and medicine

The completion of the Human Genome Project has provided a foundational blueprint for understanding human life. Nonetheless, understanding the intricate mechanisms through which our genetic blueprint is involved in disease or orchestrates development across temporal and spatial dimensions remains a profound scientific challenge. Recent breakthroughs in cellular omics technologies have paved new pathways for understanding the regulation of genomic elements and the relationship between gene expression, cellular functions, and cell fate determination. The advent of spatial omics technologies, encompassing both imaging and sequencing-based methodologies, has enabled a comprehensive understanding of biological processes from a cellular ecosystem perspective. This review offers an updated overview of how spatial omics has advanced our understanding of the translation of genetic information into cellular heterogeneity and tissue structural organization and their dynamic changes over time. It emphasizes the discovery of various biological phenomena, related to organ functionality, embryogenesis, species evolution, and the pathogenesis of diseases.

Enhancing recognition and interpretation of functional phenotypic sequences through fine-tuning pre-trained genomic models

BACKGROUND

Decoding human genomic sequences requires comprehensive analysis of DNA sequence functionality. Through computational and experimental approaches, researchers have studied the genotype-phenotype relationship and generate important datasets that help unravel complicated genetic blueprints. Thus, the recently developed artificial intelligence methods can be used to interpret the functions of those DNA sequences.

METHODS

This study explores the use of deep learning, particularly pre-trained genomic models like DNA_bert_6 and human_gpt2-v1, in interpreting and representing human genome sequences. Initially, we meticulously constructed multiple datasets linking genotypes and phenotypes to fine-tune those models for precise DNA sequence classification. Additionally, we evaluate the influence of sequence length on classification results and analyze the impact of feature extraction in the hidden layers of our model using the HERV dataset. To enhance our understanding of phenotype-specific patterns recognized by the model, we perform enrichment, pathogenicity and conservation analyzes of specific motifs in the human endogenous retrovirus (HERV) sequence with high average local representation weight (ALRW) scores.

RESULTS

We have constructed multiple genotype-phenotype datasets displaying commendable classification performance in comparison with random genomic sequences, particularly in the HERV dataset, which achieved binary and multi-classification accuracies and F1 values exceeding 0.935 and 0.888, respectively. Notably, the fine-tuning of the HERV dataset not only improved our ability to identify and distinguish diverse information types within DNA sequences but also successfully identified specific motifs associated with neurological disorders and cancers in regions with high ALRW scores. Subsequent analysis of these motifs shed light on the adaptive responses of species to environmental pressures and their co-evolution with pathogens.

CONCLUSIONS

These findings highlight the potential of pre-trained genomic models in learning DNA sequence representations, particularly when utilizing the HERV dataset, and provide valuable insights for future research endeavors. This study represents an innovative strategy that combines pre-trained genomic model representations with classical methods for analyzing the functionality of genome sequences, thereby promoting cross-fertilization between genomics and artificial intelligence.

Structural variation in humans and our primate kin in the era of telomere-to-telomere genomes and pangenomics

Structural variants (SVs) account for the majority of base pair differences both within and between primate species. However, our understanding of inter- and intra-species SV has been historically hampered by the quality of draft primate genomes and the absence of genome resources for key taxa. Recently, advances in long-read sequencing and genome assembly have begun to radically reshape our understanding of SVs. Two landmark achievements include the publication of a human telomere-to-telomere (T2T) genome as well as the development of the first human pangenome reference. In this review, we first look back to the major works laying the foundation for these projects. We then examine the ways in which T2T genome assemblies and pangenomes are transforming our understanding of and approach to primate SV. Finally, we discuss what the future of primate SV research may look like in the era of T2T genomes and pangenomics.

PhyloAln: A Convenient Reference-Based Tool to Align Sequences and High-Throughput Reads for Phylogeny and Evolution in the Omic Era

The current trend in phylogenetic and evolutionary analyses predominantly relies on omic data. However, prior to core analyses, traditional methods typically involve intricate and time-consuming procedures, including assembly from high-throughput reads, decontamination, gene prediction, homology search, orthology assignment, multiple sequence alignment, and matrix trimming. Such processes significantly impede the efficiency of research when dealing with extensive data sets. In this study, we develop PhyloAln, a convenient reference-based tool capable of directly aligning high-throughput reads or complete sequences with existing alignments as a reference for phylogenetic and evolutionary analyses. Through testing with simulated data sets of species spanning the tree of life, PhyloAln demonstrates consistently robust performance compared with other reference-based tools across different data types, sequencing technologies, coverages, and species, with percent completeness and identity at least 50 percentage points higher in the alignments. Additionally, we validate the efficacy of PhyloAln in removing a minimum of 90% foreign and 70% cross-contamination issues, which are prevalent in sequencing data but often overlooked by other tools. Moreover, we showcase the broad applicability of PhyloAln by generating alignments (completeness mostly larger than 80%, identity larger than 90%) and reconstructing robust phylogenies using real data sets of transcriptomes of ladybird beetles, plastid genes of peppers, or ultraconserved elements of turtles. With these advantages, PhyloAln is expected to facilitate phylogenetic and evolutionary analyses in the omic era. The tool is accessible at https://github.com/huangyh45/PhyloAln.

Expression of Truncated Products at the 5'-Terminal Region of RIPK2 and Evolutive Aspects that Support Their Biological Importance

Alternative splicing is the process of generating different mRNAs from the same primary transcript, which contributes to increase the transcriptome and proteome diversity. Abnormal splicing has been associated with the development of several diseases including cancer. Given that mutations and abnormal levels of the RIPK2 transcript and RIP-2 protein are frequent in tumors, and that RIP-2 modulates immune and inflammatory responses, we investigated alternative splicing events that result in partial deletions of the kinase domain at the N-terminus of RIP-2. We also investigated the structure and expression of the RIPK2 truncated variants and isoforms in different environments. In addition, we searched data throughout Supraprimates evolution that could support the biological importance of RIPK2 alternatively spliced products. We observed that human variants and isoforms were differentially regulated following temperature stress, and that the truncated transcript was more expressed than the long transcript in tumor samples. The inverse was found for the longer protein isoform. The truncated variant was also detected in chimpanzee, gorilla, hare, pika, mouse, rat, and tree shrew. The fact that the same variant has been preserved in mammals with divergence times up to 70 million years raises the hypothesis that it may have a functional significance.

Discovering Fragile Clades and Causal Sequences in Phylogenomics by Evolutionary Sparse Learning

Phylogenomic analyses of long sequences, consisting of many genes and genomic segments, reconstruct organismal relationships with high statistical confidence. But, inferred relationships can be sensitive to excluding just a few sequences. Currently, there is no direct way to identify fragile relationships and the associated individual gene sequences in species. Here, we introduce novel metrics for gene-species sequence concordance and clade probability derived from evolutionary sparse learning models. We validated these metrics using fungi, plant, and animal phylogenomic datasets, highlighting the ability of the new metrics to pinpoint fragile clades and the sequences responsible. The new approach does not necessitate the investigation of alternative phylogenetic hypotheses, substitution models, or repeated data subset analyses. Our methodology offers a streamlined approach to evaluating major inferred clades and identifying sequences that may distort reconstructed phylogenies using large datasets.

Reduction of bitter taste receptor gene family in folivorous colobine primates relative to omnivorous cercopithecine primates

Bitter taste perception is important in preventing animals from ingesting potentially toxic compounds. Whole-genome assembly (WGA) data have revealed that bitter taste receptor genes (TAS2Rs) comprise a multigene family with dozens of intact and disrupted genes in primates. However, publicly available WGA data are often incomplete, especially for multigene families. In this study, we employed a targeted capture (TC) approach specifically probing TAS2Rs for ten species of cercopithecid primates with diverse diets, including eight omnivorous cercopithecine species and two folivorous colobine species. We designed RNA probes for all TAS2Rs that we modeled to be intact in the common ancestor of cercopithecids ("ancestral-cercopithecid TAS2R gene set"). The TC was followed by short-read and high-depth massive-parallel sequencing. TC retrieved more intact TAS2R genes than found in WGA databases. We confirmed a large number of gene "births" at the common ancestor of cercopithecids and found that the colobine common ancestor and the cercopithecine common ancestor had contrasting trajectories: four gene "deaths" and three gene births, respectively. The number of intact TAS2R genes was markedly reduced in colobines (25-28 detected via TC and 20-26 detected via WGA analysis) as compared with cercopithecines (27-36 via TC and 19-30 via WGA). Birth or death events occurred at almost every phylogenetic-tree branch, making the composition of intact genes variable among species. These results show that evolutionary change in intact TAS2R genes is a complex process, refute a simple general prediction that herbivory favors more TAS2R genes, and have implications for understanding dietary adaptations and the evolution of detoxification abilities.

Adaptation of SIVmac to baboon primary cells results in complete absence of in vivo baboon infectivity

While simian immunodeficiency virus (SIV) infection is non-pathogenic in naturally infected African nonhuman primate hosts, experimental or accidental infection in rhesus macaques often leads to AIDS. Baboons, widely distributed throughout Africa, do not naturally harbor SIV, and experimental infection of baboons with SIVmac results in transient low-level viral replication. Elucidation of mechanisms of natural immunity in baboons could uncover new targets of antiviral intervention. We tested the hypothesis that an SIVmac adapted to replicate in baboon primary cells will gain the capacity to establish chronic infections in vivo. Here, we generated SIVmac variants in baboon cells through serial passage in PBMC from different donors (SIVbn-PBMC s1), in PBMC from the same donors (SIVbn-PBMC s2), or in isolated CD4 cells from the same donors used for series 2 (SIVbn-CD4). While SIVbn-PBMC s1 and SIVbn-CD4 demonstrated increased replication capacity, SIVbn-PBMC s2 did not. Pharmacological blockade of CCR5 revealed SIVbn-PBMC s1 could more efficiently use available CCR5 than SIVmac, a trait we hypothesize arose to circumvent receptor occupation by chemokines. Sequencing analysis showed that all three viruses accumulated different types of mutations, and that more non-synonymous mutations became fixed in SIVbn-PBMC s1 than SIVbn-PBMC s2 and SIVbn-CD4, supporting the notion of stronger fitness pressure in PBMC from different genetic backgrounds. Testing the individual contribution of several newly fixed SIV mutations suggested that is the additive effect of these mutations in SIVbn-PBMC s1 that contributed to its enhanced fitness, as recombinant single mutant viruses showed no difference in replication capacity over the parental SIVmac239 strain. The replicative capacity of SIVbn-PBMC passage 4 (P4) s1 was tested in vivo by infecting baboons intravenously with SIVbn-PBMC P4 s1 or SIVmac251. While animals infected with SIVmac251 showed the known pattern of transient low-level viremia, animals infected with SIVbn-PBMC P4 s1 had undetectable viremia or viral DNA in lymphoid tissue. These studies suggest that adaptation of SIV to grow in baboon primary cells results in mutations that confer increased replicative capacity in the artificial environment of cell culture but make the virus unable to avoid the restrictive factors generated by a complex multicellular organism.

Effects of nonsynonymous single nucleotide polymorphisms of the KIAA1217, SNTA1 and LTBP1 genes on the growth traits of Ujumqin sheep

Sheep body size can directly reflect the growth rates and fattening rates of sheep and is also an important index for measuring the growth performance of meat sheep. In this study, high-resolution resequencing data from four sheep breeds (Dorper sheep, Suffolk sheep, Ouessant sheep, and Shetland sheep) were analyzed. The nonsynonymous single nucleotide polymorphisms of three candidate genes (KIAA1217, SNTA1, and LTBP1) were also genotyped in 642 healthy Ujumqin sheep using MALDI-TOFMS and the genotyping results were associated with growth traits. The results showed that different genotypes of the KIAA1217 g.24429511T>C locus had significant effects on the chest circumferences of Ujumqin sheep. The SNTA1 g.62222626C>A locus had different effects on the chest depths, shoulder widths and rump widths of Ujumqin sheep. This study showed that these two sites can be used for marker-assisted selection, which will be beneficial for future precision molecular breeding.

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.