Integrating gene annotation with orthology inference at scale

Establishing genome sequencing and assembly for non-model and emerging model organisms: a brief guide

Reference genome assemblies are the basis for comprehensive genomic analyses and comparisons. Due to declining sequencing costs and growing computational power, genome projects are now feasible in smaller labs. De novo genome sequencing for non-model or emerging model organisms requires knowledge about genome size and techniques for extracting high molecular weight DNA. Next to quality, the amount of DNA obtained from single individuals is crucial, especially, when dealing with small organisms. While long-read sequencing technologies are the methods of choice for creating high quality genome assemblies, pure short-read assemblies might bear most of the coding parts of a genome but are usually much more fragmented and do not well resolve repeat elements or structural variants. Several genome initiatives produce more and more non-model organism genomes and provide rules for standards in genome sequencing and assembly. However, sometimes the organism of choice is not part of such an initiative or does not meet its standards. Therefore, if the scientific question can be answered with a genome of low contiguity in intergenic parts, missing the high standards of chromosome scale assembly should not prevent publication. This review describes how to set up an animal genome sequencing project in the lab, how to estimate costs and resources, and how to deal with suboptimal conditions. Thus, we aim to suggest optimal strategies for genome sequencing that fulfil the needs according to specific research questions, e.g. "How are species related to each other based on whole genomes?" (phylogenomics), "How do genomes of populations within a species differ?" (population genomics), "Are differences between populations relevant for conservation?" (conservation genomics), "Which selection pressure is acting on certain genes?" (identification of genes under selection), "Did repeats expand or contract recently?" (repeat dynamics).

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.

Unbiased anchors for reliable genome-wide synteny detection

Orthology inference lies at the foundation of comparative genomics research. The correct identification of loci which descended from a common ancestral sequence is not only complicated by sequence divergence but also duplication and other genome rearrangements. The conservation of gene order, i.e. synteny, is used in conjunction with sequence similarity as an additional factor for orthology determination. Current approaches, however, rely on genome annotations and are therefore limited. Here we present an annotation-free approach and compare it to synteny analysis with annotations. We find that our approach works better in closely related genomes whereas there is a better performance with annotations for more distantly related genomes. Overall, the presented algorithm offers a useful alternative to annotation-based methods and can outperform them in many cases.

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.

Quantifying the influence of genetic context on duplicated mammalian genes

Gene duplication is a fundamental part of evolutionary innovation. While single-gene duplications frequently exhibit asymmetric evolutionary rates between paralogs, the extent to which this applies to multi-gene duplications remains unclear. In this study, we investigate the role of genetic context in shaping evolutionary divergence within multi-gene duplications, leveraging microsynteny to differentiate source and target copies. Using a dataset of 193 mammalian genome assemblies and a bird outgroup, we systematically analyze patterns of sequence divergence between duplicated genes and reference orthologs. We find that target copies, those relocated to new genomic environments, exhibit elevated evolutionary rates compared to source copies in the ancestral location. This asymmetry is influenced by the distance between copies and the size of the target copy. We also demonstrate that the polarization of rate asymmetry in paralogs, the "choice" of the slowly evolving copy, is biased towards collective, block-wise polarization in multi-gene duplications. Our findings highlight the importance of genetic context in modulating post-duplication divergence, where differences in cis-regulatory elements and co-expressed gene clusters between source and target copies may be responsible. This study presents a large-scale test of asymmetric evolution in multi-gene duplications, offering new insight into how genome architecture shapes functional diversification of paralogs.

Genomic insights into marine environment adaptation and conservation of the threatened olive ridley turtle (Lepidochelys olivacea)

Sea turtles are marine flagship species and most of them are currently in a threatened state. Long-term surviving in the ocean has driven significant morphological and physiological changes for this group, which makes them an ideal model for studying adaptive evolution of marine environments. Herein, we present a chromosome-scale genome of Lepidochelys olivacea with a genome size of 2.22 Gb and a contig N50 of 97.3 Mb. Comparative genomic analyses uncovered a suite of adaptive changes in genes related to olfaction, vision, virus defense, and longevity, which may help explain the genetic underpinnings of its marine environment adaptation. We also observed that the genome-wide heterozygosity of L. olivacea was low (6.45e-4), consistent with its prolonged population decline. Overall, our study provides valuable genetic resources for understanding evolutionary adaptations to aquatic environment and for the conservation of this threatened species.

Comparative population pangenomes reveal unexpected complexity and fitness effects of structural variants

Structural variants (SVs) are widespread in vertebrate genomes, yet their evolutionary dynamics remain poorly understood. Using 45 long-read de novo genome assemblies and pangenome tools, we analyze SVs within three closely related species of North American jays (Aphelocoma, scrub-jays) displaying a 60-fold range in effective population size. We find rapid evolution of genome architecture, including ~100 Mb variation in genome size driven by dynamic satellite landscapes with unexpectedly long (> 10 kb) repeat units and widespread variation in gene content, influencing gene expression. SVs exhibit slightly deleterious dynamics modulated by variant length and population size, with strong evidence of adaptive fixation only in large populations. Our results demonstrate how population size shapes the distribution of SVs and the importance of pangenomes to characterizing genomic diversity.

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.

Bat genomes illuminate adaptations to viral tolerance and disease resistance

Zoonoses are infectious diseases transmitted from animals to humans. Bats have been suggested to harbour more zoonotic viruses than any other mammalian order1. Infections in bats are largely asymptomatic2,3, indicating limited tissue-damaging inflammation and immunopathology. To investigate the genomic basis of disease resistance, the Bat1K project generated reference-quality genomes of ten bat species, including potential viral reservoirs. Here we describe a systematic analysis covering 115 mammalian genomes that revealed that signatures of selection in immune genes are more prevalent in bats than in other mammalian orders. We found an excess of immune gene adaptations in the ancestral chiropteran branch and in many descending bat lineages, highlighting viral entry and detection factors, and regulators of antiviral and inflammatory responses. ISG15, which is an antiviral gene contributing to hyperinflammation during COVID-19 (refs. 4,5), exhibits key residue changes in rhinolophid and hipposiderid bats. Cellular infection experiments show species-specific antiviral differences and an essential role of protein conjugation in antiviral function of bat ISG15, separate from its role in secretion and inflammation in humans. Furthermore, in contrast to humans, ISG15 in most rhinolophid and hipposiderid bats has strong anti-SARS-CoV-2 activity. Our work reveals molecular mechanisms that contribute to viral tolerance and disease resistance in bats.

Orthology inference at scale with FastOMA

The surge in genome data, with ongoing efforts aiming to sequence 1.5 M eukaryotes in a decade, could revolutionize genomics, revealing the origins, evolution and genetic innovations of biological processes. Yet, traditional genomics methods scale poorly with such large datasets. Here, addressing this, 'FastOMA' provides linear scalability for orthology inference, enabling the processing of thousands of eukaryotic genomes within a day. FastOMA maintains the high accuracy and resolution of the well-established Orthologous Matrix (OMA) approach in benchmarks. FastOMA is available via GitHub at https://github.com/DessimozLab/FastOMA/ .

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.

High-quality sika deer omics data and integrative analysis reveal genic and cellular regulation of antler regeneration

The antler is the only organ that can fully regenerate annually in mammals. However, the regulatory pattern and mechanism of gene expression and cell differentiation during this process remain largely unknown. Here, we obtain comprehensive assembly and gene annotation of the sika deer (Cervus nippon) genome. We construct, together with large-scale chromatin accessibility and gene expression data, gene regulatory networks involved in antler regeneration, identifying four transcription factors, MYC, KLF4, NFE2L2, and JDP2, with high regulatory activity across the whole regeneration process. Comparative studies and luciferase reporter assay suggest the MYC expression driven by a cervid-specific regulatory element might be important for antler regenerative ability. We further develop a model called combinatorial TF Oriented Program (cTOP), which integrates single-cell data with bulk regulatory networks and find PRDM1, FOSL1, BACH1, and NFATC1 as potential pivotal factors in antler stem cell activation and osteogenic differentiation. Additionally, we uncover interactions within and between cell programs and pathways during the regeneration process. These findings provide insights into the gene and cell regulatory mechanisms of antler regeneration, particularly in stem cell activation and differentiation.

Positive Selection on Mammalian Immune Genes-Effects of Gene Function and Selective Constraint

Genome-wide analyses of various taxa have repeatedly shown that immune genes are important targets of positive selection. However, little is known about what factors determine which immune genes are under positive selection. To address this question, we here focus on the mammalian immune system and investigate the importance of gene function and other factors such as gene expression, protein-protein interactions, and overall selective constraint as determinants of positive selection. We compiled a list of >1,100 immune genes that were divided into six functional categories and analyzed using data from rodents. Genes encoding proteins that are in direct interactions with pathogens, such as pattern recognition receptors (PRRs), are often expected to be key targets of positive selection. We found that categories containing cytokines, cytokine receptors, and other cell surface proteins involved in, for example, cell-cell interactions were at least as important targets as PRRs, with three times higher rate of positive selection than nonimmune genes. The higher rate of positive selection on cytokines and cell surface proteins was partly an effect of these categories having lower selective constraint. Nonetheless, cytokines had a higher rate of positive selection than nonimmune genes even at a given level of selective constraint, indicating that gene function per se can also be a determinant of positive selection. These results have broad implications for understanding the causes of positive selection on immune genes, specifically the relative importance of host-pathogen coevolution versus other processes.

Different orthology inference algorithms generate similar predicted orthogroups among Brassicaceae species

Premise

Orthology inference is crucial for comparative genomics, and multiple algorithms have been developed to identify putative orthologs for downstream analyses. Despite the abundance of proposed solutions, including publicly available benchmarks, it is difficult to assess which tool is most suitable for plant species, which commonly have complex genomic histories.

Methods

We explored the performance of four orthology inference algorithms-OrthoFinder, SonicParanoid, Broccoli, and OrthNet-on eight Brassicaceae genomes in two groups: one group comprising only diploids and another set comprising the diploids, two mesopolyploids, and one recent hexaploid genome.

Results

The composition of the orthogroups reflected the species' ploidy and genomic histories, with the diploid set having a higher proportion of identical orthogroups. While the diploid + higher ploidy set had a lower proportion of orthogroups with identical compositions, the average degree of similarity between the orthogroups was not different from the diploid set.

Discussion

Three algorithms-OrthoFinder, SonicParanoid, and Broccoli-are helpful for initial orthology predictions. Results produced using OrthNet were generally outliers but could still provide detailed information about gene colinearity. With our Brassicaceae dataset, slight discrepancies were found across the orthology inference algorithms, necessitating additional analyses such as tree inference to fine-tune results.

De Novo Genome Assembly for an Endangered Lemur Using Portable Nanopore Sequencing in Rural Madagascar

As one of the most threatened mammalian taxa, lemurs of Madagascar are facing unprecedented anthropogenic pressures. To address conservation imperatives such as this, researchers have increasingly relied on conservation genomics to identify populations of particular concern. However, many of these genomic approaches necessitate high-quality genomes. While the advent of next-generation sequencing technologies and the resulting reduction in associated costs have led to the proliferation of genomic data and high-quality reference genomes, global discrepancies in genomic sequencing capabilities often result in biological samples from biodiverse host countries being exported to facilities in the Global North, creating inequalities in access and training within genomic research. Here, we present the first published reference genome for the endangered red-fronted brown lemur (Eulemur rufifrons) from sequencing efforts conducted entirely within the host country using portable Oxford Nanopore sequencing. Using an archived E. rufifrons specimen, we conducted long-read, nanopore sequencing at the Centre ValBio Research Station near Ranomafana National Park, in rural Madagascar, generating over 750 Gb of sequencing data from 10 MinION flow cells. Exclusively using this long-read data, we assembled 2.157 gigabase, 2980-contig nuclear assembly with an N50 of 101.6 Mb and a 17,108 bp mitogenome. The nuclear assembly had 30× average coverage and was comparable in completeness to other primate reference genomes, with a 96.1% BUSCO completeness score for primate-specific genes. As the first published reference genome for E. rufifrons and the only annotated genome available for the speciose Eulemur genus, this resource will prove vital for conservation genomic studies while our efforts exhibit the potential of this protocol to address research inequalities and build genomic capacity.

Genetic landscape in undiagnosed patients with syndromic hearing loss revealed by whole exome sequencing and phenotype similarity search

There are hundreds of rare syndromic diseases involving hearing loss, many of which are not targeted for clinical genetic testing. We systematically explored the genetic causes of undiagnosed syndromic hearing loss using a combination of whole exome sequencing (WES) and a phenotype similarity search system called PubCaseFinder. Fifty-five families with syndromic hearing loss of unknown cause were analyzed using WES after prescreening of several deafness genes depending on patient clinical features. Causative genes were identified in 22 families, including both established genes associated with syndromic hearing loss (PTPN11, CHD7, KARS1, OPA1, DLX5, MITF, SOX10, MYO7A, and USH2A) and those associated with nonsyndromic hearing loss (STRC, EYA4, and KCNQ4). Association of a DLX5 variant with incomplete partition type I (IP-I) anomaly of the inner ear was identified in a patient with cleft lip and palate and acetabular dysplasia. The study identified COL1A1, CFAP52, and NSD1 as causative genes through phenotype similarity search or by analogy. ZBTB10 was proposed as a novel candidate gene for syndromic hearing loss with IP-I. A mouse model with homozygous Zbtb10 frameshift variant resulted in embryonic lethality, suggesting the importance of this gene for early embryonic development. Our data highlight a wide spectrum of rare causative genes in patients with syndromic hearing loss, and demonstrate that WES analysis combined with phenotype similarity search is a valuable approach for clinical genetic testing of undiagnosed disease.

Screening microorganisms with robust and stable protein expression and secretion capacity

Robust and stable protein secretion is crucial for efficient recombinant protein production. Here, a novel and powerful platform using split GFP activated droplet sorting (SGADS) has been developed to significantly boost the yields of the protein of interest (POI). The SGADS platform leverages solubilizing peptide P17 and secretory expression in Bacillus subtilis to optimize two split GFP sensors: the P17-GFP1-9/GFP10-POI-GFP11 sensor for assessing protease activity and the P17-GFP1-10/GFP11-POI sensor for measuring secretion capacity. This innovative platform has demonstrated its effectiveness by successfully screening high-performance mutant strains capable of producing collagen, amylase, and protein glutaminase across a range of host organisms, including Escherichia coli, Bacillus subtilis, and Pichia pastoris. The substantial increases in production achieved with the SGADS platform highlight its broad applicability and potential in enhancing recombinant protein production.

Convergent relaxation of molecular constraint in herbivores reveals the changing role of liver and kidney functions across mammalian diets

Mammalia comprises a great diversity of diet types and associated adaptations. An understanding of the genomic mechanisms underlying these adaptations may offer insights for improving human health. Comparative genomic studies of diet that employ taxonomically restricted analyses or simplified diet classifications may suffer reduced power to detect molecular convergence associated with diet evolution. Here, we use a quantitative carnivory score-indicative of the amount of animal protein in the diet-for 80 mammalian species to detect significant correlations between the relative evolutionary rates of genes and changes in diet. We have identified six genes-ACADSB, CLDN16, CPB1, PNLIP, SLC13A2, and SLC14A2-that experienced significant changes in evolutionary constraint alongside changes in carnivory score, becoming less constrained in lineages evolving more herbivorous diets. We further consider the biological functions associated with diet evolution and observe that pathways related to amino acid and lipid metabolism, biological oxidation, and small molecule transport experienced reduced purifying selection as lineages became more herbivorous. Liver and kidney functions show similar patterns of constraint with dietary change. Our results indicate that these functions are important for the consumption of animal matter and become less important with the evolution of increasing herbivory. So, genes expressed in these tissues experience a relaxation of evolutionary constraint in more herbivorous lineages.

Orthrus: Towards Evolutionary and Functional RNA Foundation Models

In the face of rapidly accumulating genomic data, our ability to accurately predict key mature RNA properties that underlie transcript function and regulation remains limited. Pre-trained genomic foundation models offer an avenue to adapt learned RNA representations to biological prediction tasks. However, existing genomic foundation models are trained using strategies borrowed from textual or visual domains that do not leverage biological domain knowledge. Here, we introduce Orthrus, a Mamba-based mature RNA foundation model pre-trained using a novel self-supervised contrastive learning objective with biological augmentations. Orthrus is trained by maximizing embedding similarity between curated pairs of RNA transcripts, where pairs are formed from splice isoforms of 10 model organisms and transcripts from orthologous genes in 400+ mammalian species from the Zoonomia Project. This training objective results in a latent representation that clusters RNA sequences with functional and evolutionary similarities. We find that the generalized mature RNA isoform representations learned by Orthrus significantly outperform existing genomic foundation models on five mRNA property prediction tasks, and requires only a fraction of fine-tuning data to do so. Finally, we show that Orthrus is capable of capturing divergent biological function of individual transcript isoforms.

Evolutionary diversity of CXCL16-CXCR6: Convergent substitutions and recurrent gene loss in sauropsids

The CXCL16-CXCR6 axis is crucial for regulating the persistence of CD8 tissue-resident memory T cells (TRM). CXCR6 deficiency lowers TRM cell numbers in the lungs and depletes ILC3s in the lamina propria, impairing mucosal defence. This axis is linked to diseases like HIV/SIV, cancer, and COVID-19. Together, these highlight that the CXCL16-CXCR6 axis is pivotal in host immunity. Previous studies of the CXCL16-CXCR6 axis found genetic variation among species but were limited to primates and rodents. To understand the evolution and diversity of CXCL16-CXCR6 across vertebrates, we compared approximately 400 1-to-1 CXCR6 orthologs spanning diverse vertebrates. The unique DRF motif of CXCR6 facilitates leukocyte adhesion by interacting with cell surface-expressed CXCL16 and plays a key role in G-protein selectivity during receptor signalling; however, our findings show that this motif is not universal. The DRF motif is restricted to mammals, turtles, and frogs, while the DRY motif, typical in other CKRs, is found in snakes and lizards. Most birds exhibit the DRL motif. These substitutions at the DRF motif affect the receptor-Gi/o protein interaction. We establish recurrent CXCR6 gene loss in 10 out of 36 bird orders, including Galliformes and Passeriformes, Crocodilia, and Elapidae, attributed to segmental deletions and/or frame-disrupting changes. Notably, single-cell RNA sequencing of the lung shows a drop in TRM cells in species with CXCR6 loss, suggesting a possible link. The concurrent loss of ITGAE, CXCL16, and CXCR6 in chickens may have altered CD8 TRM cell abundance, with implications for immunity against viral diseases and vaccines inducing CD8 TRM cells.

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.

Practical Guidance and Workflows for Identifying Fast Evolving Non-Coding Genomic Elements Using PhyloAcc

Comparative genomics provides ample ways to study genome evolution and its relationship to phenotypic traits. By developing and testing alternate models of evolution throughout a phylogeny, one can estimate rates of molecular evolution along different lineages in a phylogeny and link these rates with observations in extant species, such as convergent phenotypes. Pipelines for such work can help identify when and where genomic changes may be associated with, or possibly influence, phenotypic traits. We recently developed a set of models called PhyloAcc, using a Bayesian framework to estimate rates of nucleotide substitution on different branches of a phylogenetic tree and evaluate their association with pre-defined or estimated phenotypic traits. PhyloAcc-ST and PhyloAcc-GT both allow users to define a priori a set of target lineages and then compare different models to identify loci accelerating in one or more target lineages. Whereas ST considers only one species tree across all input loci, GT considers alternate topologies for every locus. PhyloAcc-C simultaneously models molecular rates and rates of continuous trait evolution, allowing the user to ask whether the two are associated. Here, we describe these models and provide tips and workflows on how to prepare the input data and run PhyloAcc.

Haplotype-resolved genome and population genomics of the threatened garden dormouse in Europe

Genomic resources are important for evaluating genetic diversity and supporting conservation efforts. The garden dormouse (Eliomys quercinus) is a small rodent that has experienced one of the most severe modern population declines in Europe. We present a high-quality haplotype-resolved reference genome for the garden dormouse, and combine comprehensive short and long-read transcriptomics data sets with homology-based methods to generate a highly complete gene annotation. Demographic history analysis of the genome reveal a sharp population decline since the last interglacial, indicating an association between colder climates and population declines before anthropogenic influence. Using our genome and genetic data from 100 individuals, largely sampled in a citizen-science project across the contemporary range, we conduct the first population genomic analysis for this species. We find clear evidence for population structure across the species' core Central European range. Notably, our data show that the Alpine population, characterized by strong differentiation and reduced genetic diversity, is reproductively isolated from other regions and likely represents a differentiated evolutionary significant unit (ESU). The predominantly declining Eastern European populations also show signs of recent isolation, a pattern consistent with a range expansion from Western to Eastern Europe during the Holocene, leaving relict populations now facing local extinction. Overall, our findings suggest that garden dormouse conservation may be enhanced in Europe through the designation of ESUs.

RERconverge Expansion: Using Relative Evolutionary Rates to Study Complex Categorical Trait Evolution

Comparative genomics approaches seek to associate molecular evolution with the evolution of phenotypes across a phylogeny. Many of these methods lack the ability to analyze non-ordinal categorical traits with more than two categories. To address this limitation, we introduce an expansion to RERconverge that associates shifts in evolutionary rates with the convergent evolution of categorical traits. The categorical RERconverge expansion includes methods for performing categorical ancestral state reconstruction, statistical tests for associating relative evolutionary rates with categorical variables, and a new method for performing phylogeny-aware permutations, "permulations", on categorical traits. We demonstrate our new method on a three-category diet phenotype, and we compare its performance to binary RERconverge analyses and two existing methods for comparative genomic analyses of categorical traits: phylogenetic simulations and a phylogenetic signal based method. We present an analysis of how the categorical permulations scale with the number of species and the number of categories included in the analysis. Our results show that our new categorical method outperforms phylogenetic simulations at identifying genes and enriched pathways significantly associated with the diet phenotypes and that the categorical ancestral state reconstruction drives an improvement in our ability to capture diet-related enriched pathways compared to binary RERconverge when implemented without user input on phenotype evolution. The categorical expansion to RERconverge will provide a strong foundation for applying the comparative method to categorical traits on larger data sets with more species and more complex trait evolution than have previously been analyzed.

Gene expression differences associated with intrinsic hindfoot muscle loss in the jerboa, Jaculus jaculus

Vertebrate animals that run or jump across sparsely vegetated habitats, such as horses and jerboas, have reduced the number of distal limb bones, and many have lost most or all distal limb muscle. We previously showed that nascent muscles are present in the jerboa hindfoot at birth and that these myofibers are rapidly and completely lost soon after by a process that shares features with pathological skeletal muscle atrophy. Here, we apply an intra- and interspecies differential RNA-Seq approach, comparing jerboa and mouse muscles, to identify gene expression differences associated with the initiation and progression of jerboa hindfoot muscle loss. We show evidence for reduced hepatocyte growth factor and fibroblast growth factor signaling and an imbalance in nitric oxide signaling; all are pathways that are necessary for skeletal muscle development and regeneration. We also find evidence for phagosome formation, which hints at how myofibers may be removed by autophagy or by nonprofessional phagocytes without evidence for cell death or immune cell activation. Last, we show significant overlap between genes associated with jerboa hindfoot muscle loss and genes that are differentially expressed in a variety of human muscle pathologies and rodent models of muscle loss disorders. All together, these data provide molecular insight into the process of evolutionary and developmental muscle loss in jerboa hindfeet.

Human-specific protein-coding and lncRNA genes cast sex-biased genes in the brain and their relationships with brain diseases

BACKGROUND

Gene expression shows sex bias in the brain as it does in other organs. Since female and male humans exhibit noticeable differences in emotions, logical thinking, movement, spatial orientation, and even the incidence of neurological disorders, sex biases in the brain are especially interesting, but how they are determined, whether they are conserved or lineage specific, and what the consequences of the biases are, remain poorly explored and understood.

METHODS

Based on RNA-seq datasets from 16  and 14 brain regions in humans and macaques across developmental periods and from patients with brain diseases, we used linear mixed models (LMMs) to differentiate variations in gene expression caused by factors of interest and confounding factors and identify four types of sex-biased genes. Effect size and confidence in each effect were measured upon the local false sign rate (LFSR). We utilized the biomaRt R package to acquire orthologous genes in humans and macaques from the BioMart Ensembl website. Transcriptional regulation of sex-biased genes by sex hormones and lncRNAs were analyzed using the CellOracle, GENIE3, and Longtarget programs. Sex-biased genes' functions were revealed by gene set enrichment analysis using multiple methods.

RESULTS

Lineage-specific sex-biased genes greatly determine the distinct sex biases in human and macaque brains. In humans, those encoding proteins contribute directly to immune-related functions, and those encoding lncRNAs intensively regulate the expression of other sex-biased genes, especially genes with immune-related functions. The identified sex-specific differentially expressed genes (ssDEGs) upon gene expression in disease and normal samples also indicate that protein-coding ssDEGs are conserved in humans and macaques but that lncRNA ssDEGs are not conserved. The results answer the above questions, reveal an intrinsic relationship between sex biases in the brain and sex-biased susceptibility to brain diseases, and will help researchers investigate human- and sex-specific ncRNA targets for brain diseases.

CONCLUSIONS

Human-specific genes greatly cast sex-biased genes in the brain and their relationships with brain diseases, with protein-coding genes contributing to immune response related functions and lncRNA genes critically regulating sex-biased genes. The high proportions of lineage-specific lncRNAs in mammalian genomes indicate that sex biases may have evolved rapidly in not only the brain but also other organs.

The DoGA consortium expression atlas of promoters and genes in 100 canine tissues

The dog, Canis lupus familiaris, is an important model for studying human diseases. Unlike many model organisms, the dog genome has a comparatively poor functional annotation, which hampers gene discovery for development, morphology, disease, and behavior. To fill this gap, we established a comprehensive tissue biobank for both the dog and wolf samples. The biobank consists of 5485 samples representing 132 tissues from 13 dogs, 12 dog embryos, and 24 wolves. In a subset of 100 tissues from nine dogs and 12 embryos, we characterized gene expression activity for each promoter, including alternative and novel, i.e., previously not annotated, promoter regions, using the 5' targeting RNA sequencing technology STRT2-seq. We identified over 100,000 promoter region candidates in the recent canine genome assembly, CanFam4, including over 45,000 highly reproducible sites with gene expression and respective tissue enrichment levels. We provide a promoter and gene expression atlas with interactive, open data resources, including a data coordination center and genome browser track hubs. We demonstrated the applicability of Dog Genome Annotation (DoGA) data and resources using multiple examples spanning canine embryonic development, morphology and behavior, and diseases across species.

Quest for Orthologs in the Era of Biodiversity Genomics

The era of biodiversity genomics is characterized by large-scale genome sequencing efforts that aim to represent each living taxon with an assembled genome. Generating knowledge from this wealth of data has not kept up with this pace. We here discuss major challenges to integrating these novel genomes into a comprehensive functional and evolutionary network spanning the tree of life. In summary, the expanding datasets create a need for scalable gene annotation methods. To trace gene function across species, new methods must seek to increase the resolution of ortholog analyses, e.g. by extending analyses to the protein domain level and by accounting for alternative splicing. Additionally, the scope of orthology prediction should be pushed beyond well-investigated proteomes. This demands the development of specialized methods for the identification of orthologs to short proteins and noncoding RNAs and for the functional characterization of novel gene families. Furthermore, protein structures predicted by machine learning are now readily available, but this new information is yet to be integrated with orthology-based analyses. Finally, an increasing focus should be placed on making orthology assignments adhere to the findable, accessible, interoperable, and reusable (FAIR) principles. This fosters green bioinformatics by avoiding redundant computations and helps integrating diverse scientific communities sharing the need for comparative genetics and genomics information. It should also help with communicating orthology-related concepts in a format that is accessible to the public, to counteract existing misinformation about evolution.

Challenges in Assembling the Dated Tree of Life

The assembly of a comprehensive and dated Tree of Life (ToL) remains one of the most formidable challenges in evolutionary biology. The complexity of life's history, involving both vertical and horizontal transmission of genetic information, defies its representation by a simple bifurcating phylogeny. With the advent of genome and metagenome sequencing, vast amounts of data have become available. However, employing this information for phylogeny and divergence time inference has introduced significant theoretical and computational hurdles. This perspective addresses some key methodological challenges in assembling the dated ToL, namely, the identification and classification of homologous genes, accounting for gene tree-species tree mismatch due to population-level processes along with duplication, loss, and horizontal gene transfer, and the accurate dating of evolutionary events. Ultimately, the success of this endeavor requires new approaches that integrate knowledge databases with optimized phylogenetic algorithms capable of managing complex evolutionary models.

Evolutionary Innovations in Conserved Regulatory Elements Associate With Developmental Genes in Mammals

Transcriptional enhancers orchestrate cell type- and time point-specific gene expression programs. Genetic variation within enhancer sequences is an important contributor to phenotypic variation including evolutionary adaptations and human disease. Certain genes and pathways may be more prone to regulatory evolution than others, with different patterns across diverse organisms, but whether such patterns exist has not been investigated at a sufficient scale. To address this question, we identified signatures of accelerated sequence evolution in conserved enhancer elements throughout the mammalian phylogeny at an unprecedented scale. While different genes and pathways were enriched for regulatory evolution in different parts of the tree, we found a striking overall pattern of pleiotropic genes involved in gene regulatory and developmental processes being enriched for accelerated enhancer evolution. These genes were connected to more enhancers than other genes, which was the basis for having an increased amount of sequence acceleration over all their enhancers combined. We provide evidence that sequence acceleration is associated with turnover of regulatory function. Detailed study of one acceleration event in an enhancer of HES1 revealed that sequence evolution led to a new activity domain in the developing limb that emerged concurrently with the evolution of digit reduction in hoofed mammals. Our results provide evidence that enhancer evolution has been a frequent contributor to regulatory innovation at conserved developmental signaling genes in mammals.

Chromosome-level genome assembly of the bay scallop Argopecten irradians

The bay scallop, Argopecten irradians, is a species of major commercial, cultural, and ecological importance. It is endemic to the eastern coast of the United States, but has also been introduced to China, where it supports a significant aquaculture industry. Here, we provide an annotated chromosome-level reference genome assembly for the bay scallop, assembled using PacBio and Hi-C data. The total genome size is 845.9 Mb, distributed over 1,503 scaffolds with a scaffold N50 of 44.3 Mb. The majority (92.9%) of the assembled genome is contained within the 16 largest scaffolds, corresponding to the 16 chromosomes confirmed by Hi-C analysis. The assembly also includes the complete mitochondrial genome. Approximately 36.2% of the genome consists of repetitive elements. The BUSCO analysis showed a completeness of 96.2%. We identified 33,772 protein-coding genes. This genome assembly will be a valuable resource for future research on evolutionary dynamics, adaptive mechanisms, and will support genome-assisted breeding, contributing to the conservation and management of this iconic species in the face of environmental and pathogenic challenges.

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 assembled 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.

Cell type and cell signaling innovations underlying mammalian pregnancy

How fetal and maternal cell types have co-evolved to enable mammalian placentation poses a unique evolutionary puzzle. Here, we present a multi-species atlas integrating single-cell transcriptomes from six species bracketing therian mammal diversity. We find that invasive trophoblasts share a gene-expression signature across eutherians, and evidence that endocrine decidual cells evolved stepwise from an immunomodulatory cell type retained in Tenrec with affinity to human decidua of menstruation. We recover evolutionary patterns in ligand-receptor signaling: fetal and maternal cells show a pronounced tendency towards disambiguation, but a predicted arms race dynamic between them is limited. We reconstruct cell communication networks of extinct mammalian ancestors, finding strong integration of fetal trophoblast into maternal networks. Together, our results reveal a dynamic history of cell type and signaling evolution.

Synopsis

The fetal-maternal interface is one of the most intense loci of cell-cell signaling in the human body. Invasion of cells from the fetal placenta into the uterus, and the corresponding transformation of maternal tissues called decidualization, first evolved in the stem lineage of eutherian mammals( 1 , 2 ). Single-cell studies of the human fetal-maternal interface have provided new insight into the cell type diversity and cell-cell interactions governing this chimeric organ( 3-5 ). However, the fetal-maternal interface is also one of the most rapidly evolving, and hence most diverse, characters among mammals( 6 ), and an evolutionary analysis is missing. Here, we present and compare single-cell data from the fetal-maternal interface of species bracketing key events in mammal phylogeny: a marsupial (opossum, Monodelphis domestica ), the afrotherian Tenrec ecaudatus, and four Euarchontoglires - guinea pig and mouse (Rodentia) together with recent macaque and human data (primates) ( 4 , 5 , 7 ). We infer cell type homologies, identify a gene-expression signature of eutherian invasive trophoblast conserved over 99 million years, and discover a predecidual cell in the tenrec which suggests stepwise evolution of the decidual stromal cell. We reconstruct ancestral cell signaling networks, revealing the integration of fetal cell types into the interface. Finally, we test two long-standing theoretical predictions, the disambiguation hypothesis( 8 ) and escalation hypothesis( 9 ), at transcriptome-wide scale, finding divergence between fetal and maternal signaling repertoires but arms race dynamics restricted to a small subset of ligand-receptor pairs. In so doing, we trace the co-evolutionary history of cell types and their signaling across mammalian viviparity.

A comparative analysis of planarian genomes reveals regulatory conservation in the face of rapid structural divergence

The planarian Schmidtea mediterranea is being studied as a model species for regeneration, but the assembly of planarian genomes remains challenging. Here, we report a high-quality haplotype-phased, chromosome-scale genome assembly of the sexual S2 strain of S. mediterranea and high-quality chromosome-scale assemblies of its three close relatives, S. polychroa, S. nova, and S. lugubris. Using hybrid gene annotations and optimized ATAC-seq and ChIP-seq protocols for regulatory element annotation, we provide valuable genome resources for the planarian research community and a first comparative perspective on planarian genome evolution. Our analyses reveal substantial divergence in protein-coding sequences and regulatory regions but considerable conservation within promoter and enhancer annotations. We also find frequent retrotransposon-associated chromosomal inversions and interchromosomal translocations within the genus Schmidtea and, remarkably, independent and nearly complete losses of ancestral metazoan synteny in Schmidtea and two other flatworm groups. Overall, our results suggest that platyhelminth genomes can evolve without syntenic constraints.

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.

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.

A haplotype-resolved reference genome of a long-distance migratory bat, Pipistrellus nathusii (Keyserling & Blasius, 1839)

We present a complete, chromosome-scale reference genome for the long-distance migratory bat Pipistrellus nathusii. The genome encompasses both haplotypic sets of autosomes and the separation of both sex chromosomes by utilizing highly accurate long-reads and preserving long-range phasing information through the use of three-dimensional chromatin conformation capture sequencing (Hi-C). This genome, accompanied by a comprehensive protein-coding sequence annotation, provides a valuable genomic resource for future investigations into the genomic bases of long-distance migratory flight in bats as well as uncovering the genetic architecture, population structure and evolutionary history of Pipistrellus nathusii. The reference-quality genome presented here gives a fundamental resource to further our understanding of bat genetics and evolution, adding to the growing number of high-quality genetic resources in this field. Here, we demonstrate its use in the phylogenetic reconstruction of the order Chiroptera, and in particular, we present the resources to allow detailed investigations into the genetic drivers and adaptations related to long-distance migration.

The genome of African manatee Trichechus senegalensis reveals secondary adaptation to the aquatic environment

Sirenians exhibit unique aquatic adaptations, showcasing both convergent adaptive features shared with cetaceans and unique characteristics such as cold sensitivity and dense bones. Here, we report a chromosome-level genome of the African manatee (Trichechus senegalensis) with high continuity, completeness, and accuracy. We found that genes associated with osteopetrosis have undergone positive selection (CSF1R and LRRK1) or pseudogenized (FAM111A and IGSF23) in the African manatee, potentially contributing to the dense bone formation. The loss of KCNK18 may have increased their sensitivity to cold water temperatures. Moreover, we identified convergent evolutionary signatures in 392 genes among fully aquatic mammals, primarily enriched in skin or skeletal system development and circadian rhythm, which contributed to the transition from terrestrial to fully aquatic lifestyles. The African manatee currently possesses a small effective population size and low genome-wide heterozygosity. Overall, our study provides genetic resources for understanding the evolutionary characteristics and conservation efforts of this species.

Three-dimensional genome architecture persists in a 52,000-year-old woolly mammoth skin sample

Analyses of ancient DNA typically involve sequencing the surviving short oligonucleotides and aligning to genome assemblies from related, modern species. Here, we report that skin from a female woolly mammoth (†Mammuthus primigenius) that died 52,000 years ago retained its ancient genome architecture. We use PaleoHi-C to map chromatin contacts and assemble its genome, yielding 28 chromosome-length scaffolds. Chromosome territories, compartments, loops, Barr bodies, and inactive X chromosome (Xi) superdomains persist. The active and inactive genome compartments in mammoth skin more closely resemble Asian elephant skin than other elephant tissues. Our analyses uncover new biology. Differences in compartmentalization reveal genes whose transcription was potentially altered in mammoths vs. elephants. Mammoth Xi has a tetradic architecture, not bipartite like human and mouse. We hypothesize that, shortly after this mammoth's death, the sample spontaneously freeze-dried in the Siberian cold, leading to a glass transition that preserved subfossils of ancient chromosomes at nanometer scale.

Unraveling genomic features and phylogenomics through the analysis of three Mexican endemic Myotis genomes

Background

Genomic resource development for non-model organisms is rapidly progressing, seeking to uncover molecular mechanisms and evolutionary adaptations enabling thriving in diverse environments. Limited genomic data for bat species hinder insights into their evolutionary processes, particularly within the diverse Myotis genus of the Vespertilionidae family. In Mexico, 15 Myotis species exist, with three-M. vivesi, M. findleyi, and M. planiceps-being endemic and of conservation concern.

Methods

We obtained samples of Myotis vivesi, M. findleyi, and M. planiceps for genomic analysis. Each of three genomic DNA was extracted, sequenced, and assembled. The scaffolding was carried out utilizing the M. yumanensis genome via a genome-referenced approach within the ntJoin program. GapCloser was employed to fill gaps. Repeat elements were characterized, and gene prediction was done via ab initio and homology methods with MAKER pipeline. Functional annotation involved InterproScan, BLASTp, and KEGG. Non-coding RNAs were annotated with INFERNAL, and tRNAscan-SE. Orthologous genes were clustered using Orthofinder, and a phylogenomic tree was reconstructed using IQ-TREE.

Results

We present genome assemblies of these endemic species using Illumina NovaSeq 6000, each exceeding 2.0 Gb, with over 90% representing single-copy genes according to BUSCO analyses. Transposable elements, including LINEs and SINEs, constitute over 30% of each genome. Helitrons, consistent with Vespertilionids, were identified. Values around 20,000 genes from each of the three assemblies were derived from gene annotation and their correlation with specific functions. Comparative analysis of orthologs among eight Myotis species revealed 20,820 groups, with 4,789 being single copy orthogroups. Non-coding RNA elements were annotated. Phylogenomic tree analysis supported evolutionary chiropterans' relationships. These resources contribute significantly to understanding gene evolution, diversification patterns, and aiding conservation efforts for these endangered bat species.

The evolution history of an allotetraploid mangrove tree analysed with a new tool Allo4D

Mangrove species are broadly classified as true mangroves and mangrove associates. The latter are amphibious plants that can survive in the intertidal zone and reproduce naturally in terrestrial environments. Their widespread distribution and extensive adaptability make them ideal research materials for exploring adaptive evolution. In this study, we de novo assembled two genomes of mangrove associates (the allotetraploid Barringtonia racemosa (2n = 4x = 52) and diploid Barringtonia asiatica (2n = 2x = 26)) to investigate the role of allopolyploidy in the evolutionary history of mangrove species. We developed a new allotetraploid-dividing tool Allo4D to distinguish between allotetraploid scaffold-scale subgenomes and verified its accuracy and reliability using real and simulated data. According to the two subgenomes of allotetraploid B. racemosa divided using Allo4D, the allopolyploidization event was estimated to have occurred approximately one million years ago (Mya). We found that B. racemosa, B. asiatica, and Diospyros lotus shared a whole genome duplication (WGD) event during the K-Pg (Cretaceous-Paleozoic) period. K-Pg WGD and recent allopolyploidization events contributed to the speciation of B. racemosa and its adaptation to coastal habitats. We found that genes in the glucosinolates (GSLs) pathway, an essential pathway in response to various biotic and abiotic stresses, expanded rapidly in B. racemosa during polyploidization. In summary, this study provides a typical example of the adaptation of allopolyploid plants to extreme environmental conditions. The newly developed tool, Allo4D, can effectively divide allotetraploid subgenomes and explore the evolutionary history of polyploid plants, especially for species whose ancestors are unknown or extinct.

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.

Genomic signatures associated with recurrent scale loss in cyprinid fish

Scale morphology represents a fundamental feature of fish and a key evolutionary trait underlying fish diversification. Despite frequent and recurrent scale loss throughout fish diversification, comprehensive genome-wide analyses of the genomic signatures associated with scale loss in divergent fish lineages remain scarce. In the current study, we investigated genome-wide signatures, specifically convergent protein-coding gene loss, amino acid substitutions, and cis-regulatory sequence changes, associated with recurrent scale loss in two divergent Cypriniformes lineages based on large-scale genomic, transcriptomic, and epigenetic data. Results demonstrated convergent changes in many genes related to scale formation in divergent scaleless fish lineages, including loss of P/Q-rich scpp genes (e.g. scpp6 and scpp7), accelerated evolution of non-coding elements adjacent to the fgf and fgfr genes, and convergent amino acid changes in genes (e.g. snap29) under relaxed selection. Collectively, these findings highlight the existence of a shared genetic architecture underlying recurrent scale loss in divergent fish lineages, suggesting that evolutionary outcomes may be genetically repeatable and predictable in the convergence of scale loss in fish.

A chromosome-level genome assembly of an avivorous bat species (Nyctalus aviator)

Currently, three carnivorous bat species, namely Ia io, Nyctalus lasiopterus, and Nyctalus aviator, are known to actively prey on seasonal migratory birds (hereinafter referred to as "avivorous bats"). However, the absence of reference genomes impedes a thorough comprehension of the molecular adaptations of avivorous bat species. Herein, we present the high-quality chromosome-scale reference genome of N. aviator based on PacBio subreads, DNBSEQ short-reads and Hi-C sequencing data. The genome assembly size of N. aviator is 1.77 Gb, with a scaffold N50 of 102 Mb, of which 99.8% assembly was anchored into 21 pseudo-chromosomes. After masking 635.1 Mb repetitive sequences, a total of 19,412 protein-coding genes were identified, of which 99.3% were functionally annotated. The genome assembly and gene prediction reached 96.1% and 96.1% completeness of Benchmarking Universal Single-Copy Orthologs (BUSCO), respectively. This chromosome-level reference genome of N. aviator fills a gap in the existing information on the genomes of carnivorous bats, especially avivorous ones, and will be valuable for mechanism of adaptations to dietary niche expansion in bat species.

The NCBI Comparative Genome Viewer (CGV) is an interactive visualization tool for the analysis of whole-genome eukaryotic alignments

We report a new visualization tool for analysis of whole-genome assembly-assembly alignments, the Comparative Genome Viewer (CGV) (https://ncbi.nlm.nih.gov/genome/cgv/). CGV visualizes pairwise same-species and cross-species alignments provided by National Center for Biotechnology Information (NCBI) using assembly alignment algorithms developed by us and others. Researchers can examine large structural differences spanning chromosomes, such as inversions or translocations. Users can also navigate to regions of interest, where they can detect and analyze smaller-scale deletions and rearrangements within specific chromosome or gene regions. RefSeq or user-provided gene annotation is displayed where available. CGV currently provides approximately 800 alignments from over 350 animal, plant, and fungal species. CGV and related NCBI viewers are undergoing active development to further meet needs of the research community in comparative genome visualization.

Evolution of a novel adrenal cell type that promotes parental care

Cell types with specialized functions fundamentally regulate animal behaviour, and yet the genetic mechanisms that underlie the emergence of novel cell types and their consequences for behaviour are not well understood1. Here we show that the monogamous oldfield mouse (Peromyscus polionotus) has recently evolved a novel cell type in the adrenal gland that expresses the enzyme AKR1C18, which converts progesterone into 20α-hydroxyprogesterone. We then demonstrate that 20α-hydroxyprogesterone is more abundant in oldfield mice, where it induces monogamous-typical parental behaviours, than in the closely related promiscuous deer mice (Peromyscus maniculatus). Using quantitative trait locus mapping in a cross between these species, we ultimately find interspecific genetic variation that drives expression of the nuclear protein GADD45A and the glycoprotein tenascin N, which contribute to the emergence and function of this cell type in oldfield mice. Our results provide an example by which the recent evolution of a new cell type in a gland outside the brain contributes to the evolution of social behaviour.

Pervasive loss of regulated necrotic cell death genes in elephants, hyraxes, and sea cows (Paenungualta)

Gene loss can promote phenotypic differences between species, for example, if a gene constrains phenotypic variation in a trait, its loss allows for the evolution of a greater range of variation or even new phenotypes. Here, we explore the contribution of gene loss to the evolution of large bodies and augmented cancer resistance in elephants. We used genomes from 17 Afrotherian and Xenarthran species to identify lost genes, i.e., genes that have pseudogenized or been completely lost, and Dollo parsimony to reconstruct the evolutionary history of gene loss across species. We unexpectedly discovered a burst of gene losses in the Afrotherian stem lineage and found that the loss of genes with functions in regulated necrotic cell death modes was pervasive in elephants, hyraxes, and sea cows (Paenungulata). Among the lost genes are MLKL and RIPK3, which mediate necroptosis, and sensors that activate inflammasomes to induce pyroptosis, including AIM2, MEFV, NLRC4, NLRP1, and NLRP6. These data suggest that the mechanisms that regulate necrosis and pyroptosis are either extremely derived or potentially lost in these lineages, which may contribute to the repeated evolution of large bodies and cancer resistance in Paenungulates as well as susceptibility to pathogen infection.

Jump-starting life: balancing transposable element co-option and genome integrity in the developing mammalian embryo

Remnants of transposable elements (TEs) are widely expressed throughout mammalian embryo development. Originally infesting our genomes as selfish elements and acting as a source of genome instability, several of these elements have been co-opted as part of a complex system of genome regulation. Many TEs have lost transposition ability and their transcriptional potential has been tampered as a result of interactions with the host throughout evolutionary time. It has been proposed that TEs have been ultimately repurposed to function as gene regulatory hubs scattered throughout our genomes. In the early embryo in particular, TEs find a perfect environment of naïve chromatin to escape transcriptional repression by the host. As a consequence, it is thought that hosts found ways to co-opt TE sequences to regulate large-scale changes in chromatin and transcription state of their genomes. In this review, we discuss several examples of TEs expressed during embryo development, their potential for co-option in genome regulation and the evolutionary pressures on TEs and on our genomes.

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.

A High-Quality Blue Whale Genome, Segmental Duplications, and Historical Demography

The blue whale, Balaenoptera musculus, is the largest animal known to have ever existed, making it an important case study in longevity and resistance to cancer. To further this and other blue whale-related research, we report a reference-quality, long-read-based genome assembly of this fascinating species. We assembled the genome from PacBio long reads and utilized Illumina/10×, optical maps, and Hi-C data for scaffolding, polishing, and manual curation. We also provided long read RNA-seq data to facilitate the annotation of the assembly by NCBI and Ensembl. Additionally, we annotated both haplotypes using TOGA and measured the genome size by flow cytometry. We then compared the blue whale genome with other cetaceans and artiodactyls, including vaquita (Phocoena sinus), the world's smallest cetacean, to investigate blue whale's unique biological traits. We found a dramatic amplification of several genes in the blue whale genome resulting from a recent burst in segmental duplications, though the possible connection between this amplification and giant body size requires further study. We also discovered sites in the insulin-like growth factor-1 gene correlated with body size in cetaceans. Finally, using our assembly to examine the heterozygosity and historical demography of Pacific and Atlantic blue whale populations, we found that the genomes of both populations are highly heterozygous and that their genetic isolation dates to the last interglacial period. Taken together, these results indicate how a high-quality, annotated blue whale genome will serve as an important resource for biology, evolution, and conservation research.