Tree of life reveals clock-like speciation and diversification

Complete Mitochondrial Genomes of Endemic Charrs of Salvelinus: Implications for Taxonomy and Species Identification

Gene fragments on which databases usually rely are not always suited for identification of endemic species of disputable origin. Clarification of the unresolved issues requires more comprehensive mitochondrial genome-based phylogenetics to enhance understanding of the relationship and delineation of species as a basis for their conservation. This study focuses on endemic charrs (Salvelinus, Salmonidae). The complete mitochondrial genomes (mitogenomes) of Longhead charr Salvelinus kronocius and Schmidt's charr Salvelinus schmidti from Lake Kronotskoye (Kamchatka Peninsula, Russia) were sequenced and compared with the mitogenomes of a putative ancestor species, Northern Dolly Varden Salvelinus malma malma. We analysed 36 mitogenomes, including six ones sequenced in this study, 26 charr sequences published in GenBank, and four outgroup sequences. We obtained additional evidence that the lacustrine charrs belong to the Beringian lineage with very low level of mitogenome divergence. Our findings demonstrate that the mitogenome sequencing serves as a valuable tool in clarifying the taxonomic status of narrow-range endemics or rare charrs and may be useful for future genetic investigations of charrs and other groups with high level of diversity toward proper assessment of the status of the controversial taxa. Expanding genetic resources is necessary for species identification and biodiversity monitoring, as well as for conservation measures aimed at controlling adverse natural and anthoropogenic impacts on fish population systems.

Urea Amidolyase as an Enzyme for Urea Utilisation in Phytoplankton: Functional Display in Chlamydomonas reinhardtii

Urea is an important source of nitrogen for many phytoplankton with the potential to stimulate harmful algal blooms, but the molecular machinery underpinning urea uptake and assimilation by algae is not fully understood. Urease (URE) is commonly regarded as the responsible enzyme, but urea amidolyase (UAL), albeit known to exist, has hardly been studied. Here, the species distribution, expression patterns and functional roles of UAL are examined. We found a widespread occurrence of UAL across six major phytoplankton lineages, along with evidence of a potential URE-independent evolutionary trajectory and lineage-specific losses. Quantitative analyses based on marine planktonic metagenomes and metatranscriptomes revealed that UAL is as prevalent as URE, but exhibits higher expression levels in phytoplankton than in bacteria, suggesting that UAL plays a crucial role in nitrogen nutrition in marine phytoplankton. Furthermore, using the CRISPR/Cas9 genome editing method and Chlamydomonas reinhardtii as the algal model, we showed that DUR2 in UAL is essential for urea utilisation, as its knockout completely abolishes the ability of algae to grow under urea as the sole nitrogen source. This study unveils an unappreciated mechanism in algae for utilising urea as a nutrient, underscores the need to consider both URE and UAL enzyme systems to model urea utilisation by algae and provides a crucial gene (DUR2) as a potential genetic marker for detecting the contribution of UAL to urea utilisation in phytoplankton.

Migration Strategy and Diet Affect the Metabolism of Passerine Birds During Endurance Flight

Bird migration varies greatly in overall distance and length of single flight bouts. Therefore, we expect that metabolic adaptations may also differ widely among migrants. Endurance flight is mainly fuelled by fat and complemented by protein. The proportions of lipids and protein accumulated before flights, and spent during flight, depend on food type. A fruit diet facilitates fattening more than a proteinous arthropod or seed diet. Adaptations to maximize lipid use during flight vary with the length of flight bouts. We expect that the type of diet and migration strategy (length of flight bouts, overall migration distance) affects flight metabolism. On a Swiss Alpine pass, we caught 30 species of nocturnal and diurnal migrant passerines out of natural migratory flight and compared them with conspecifics kept inactive. We examined the effects of migration strategy, primary diet, and body fat stores on plasma concentrations of six metabolites of the fat, protein, and carbohydrate metabolism, used as indicators of relative fuel use. During migratory flight, immediate migration strategy (short hops during day or long bouts during night), general migration strategy (long- and short-distance migration) and diet affected metabolite levels, while fat stores had an additional effect. Triglyceride and free fatty acid levels were high in night-migrants and frugivores. Uric acid and glucose levels were high in insectivores. Glucose, ß-hydroxy-butyrate, and glycerol were less dependent on day/night-migration or diet. The metabolic profiles indicate that migrant passerines vary in the degree of fat use depending on migration strategy, diet, and current fat stores. Day-migrating short-distance migrant insectivores and granivores used protein or glycogen to a higher degree than night migrants. Frugivores maximized fat use. Long-distance migrants favored fat use even further. Hence, long-distance night-migrant frugivores with high fat loads appear best adapted for fat use during migratory flight.

Genomic, Phenotypic and Environmental Correlates of Speciation in the Midwife Toads (Alytes)

Speciation, i.e., the formation of new species, implies that diverging populations evolve genetic, phenotypic or ecological factors that promote reproductive isolation (RI), but the relative contributions of these factors remain elusive. Here we test which of genomic, bioacoustic, morphological, and environmental differences best predicts RI across a continuum of divergence in the midwife toads (genus Alytes), a group of Western Mediterranean amphibians, using a total evidence approach. We found that, without strong geographic barriers to dispersal, the extent of introgression across hybrid zones between phylogeographic lineages, which should reflect the strength of RI, predominantly covaries with genomic divergence. Overall phenotypic differentiation becomes substantial only between well established, fully isolated species. These results suggest that speciation in midwife toads initially involve cryptic lineages, which probably evolve RI through intrinsic (genetic) hybrid incompatibilities. As they continue to diverge, these nascent species eventually differentiate externally, which potentially enforces pre-mating barriers and facilitates sympatry. This speciation scenario has practical implications for species delimitation, notably when using hybrid zones and divergence thresholds as proxies for reproductive isolation.

The Candida Genome Database: annotation and visualization updates

The Candida Genome Database (CGD; www.candidagenome.org) is unique in being both a model organism database and a fungal pathogen database. As a fungal pathogen database, CGD hosts locus pages for 5 species of the best-studied pathogenic fungi in the Candida group. As a model organism database, the species Candida albicans serves as a model both for other Candida spp. and for non-Candida fungi that form biofilms and undergo routine morphogenic switching from the planktonic form to the filamentous form, which is not done by other model yeasts. As pathogenic Candida species have become increasingly drug resistant, the high lethality of invasive candidiasis in immunocompromised people is increasingly alarming. There is a pressing need for additional research into basic Candida biology, epidemiology and phylogeny, and potential new antifungals. CGD serves the needs of this diverse research community by curating the entire gene-based Candida experimental literature as it is published, extracting, organizing, and standardizing gene annotations. Gene pages were added for the species Candida auris, recent outbreaks of which have been labeled an "urgent" threat. Most recently, we have begun linking clinical data on disease to relevant Literature Topics to improve searchability for clinical researchers. Because CGD curates for multiple species and most research focuses on aspects related to pathogenicity, we focus our curation efforts on assigning Literature Topic tags, collecting detailed mutant phenotype data, and assigning controlled Gene Ontology terms with accompanying evidence codes. Our Summary pages for each feature include the primary name and all aliases for that locus, a description of the gene and/or gene product, detailed ortholog information with links, a JBrowse window with a visual view of the gene on its chromosome, summarized phenotype, Gene Ontology, and sequence information, references cited on the summary page itself, and any locus notes. The database serves as a community hub, where we link to various types of reference material of relevance to Candida researchers, including colleague information, news, and notice of upcoming meetings. We routinely survey the community to learn how the field is evolving and how needs may have changed. For example, we asked our users which species we should next add to CGD, and the clear answer was Candida tropicalis. A key future challenge is management of the flood of high-throughput expression data to make it as useful as possible to as many researchers as possible. The central challenge for any community database is to turn data into knowledge, which the community can access, use, and build upon.

Fundamental constraints on vertebrate life history are shaped by aquatic-terrestrial transitions and reproductive mode

Vertebrate life histories evolve in response to selection imposed by abiotic and biotic environmental conditions while being limited by genetic, developmental, physiological, demographic and phylogenetic processes that constrain adaptation. Despite the well-recognized shifts in selective pressures accompanying transitions among environments, the conditions driving innovation and the consequences for life-history evolution remain outstanding questions. Here we compare the traits of vertebrates that occupy aquatic or terrestrial environments as juveniles to infer shifts in evolutionary constraints that explain differences in their life-history traits and thus their fundamental demographic rates. Our results emphasize the reduced potential for life-history diversification on land, especially that of reproductive strategies, which limits the scope of viable life-history strategies. Moreover, our study reveals differences between the evolution of viviparity in aquatic and terrestrial realms. Transitions from egg laying to live birth represent a major shift across life-history space for aquatic organisms, whereas terrestrial egg-laying organisms evolve live birth without drastic changes in life-history strategy. Whilst trade-offs in the allocation of resources place fundamental constraints on the way life histories can vary, ecological setting influences the position of species within the viable phenotypic space available for adaptive evolution.

Constraint of accessible chromatins maps regulatory loci involved in maize speciation and domestication

Comparative genomic studies can identify genes under evolutionary constraint or specialized for trait innovation. Growing evidence suggests that evolutionary constraint also acts on non-coding regulatory sequences, exerting significant impacts on fitness-related traits, although it has yet to be thoroughly explored in plants. Using the assay for transposase-accessible chromatin by sequencing (ATAC-seq), we profile over 80,000 maize accessible chromatin regions (ACRs), revealing that ACRs evolve faster than coding genes, with about one-third being maize-specific and regulating genes associated with speciation. We highlight the role of transposable elements (TEs) in driving intraspecific innovation of ACRs and identify hundreds of candidate ACRs potentially involved in transcriptional rewiring during maize domestication. Additionally, we demonstrate the importance of accessible chromatin in maintaining subgenome dominance and controlling complex trait variations. This study establishes a framework for analyzing the evolutionary trajectory of plant regulatory sequences and offers candidate loci for downstream exploration and application in maize breeding.

Phylogenetics, biogeography, and life history evolution in the broadly distributed treefrog genus Dendropsophus (Anura: Hylidae: Hylinae)

Dendropsophusis one of the most species-rich genera of hylid treefrogs. Recent studies integrating Sanger-generated mitochondrial and nuclear loci with phenomic characters (SP) have advanced understanding of this clade, but questions about its internal relationships and biogeographic history persist. To address these questions, we used anchored hybrid enrichment (AHE) to combine 432 nuclear loci for 78 taxa (72 % of species) with published data. Quantitatively, the impact of the AHE data was modest, with compositional differences in only three recognized clades and more than 80 % of the clades in the AHE + SP analyses also supported in the SP-only analyses. Nevertheless, the impact of AHE was crucial for resolving and increasing support for multiple nodes. We transferred one species of the formerD. ruschiigroup to theD. decipiensgroup and redefined theD. leucophyllatusgroup to avoid paraphyly. We estimated divergence times to reconstruct the clade's biogeographic history. We also examined evolution of oviposition sites and assessed its effect on lineage accumulation. Dendropsophuslikely originated ∼ 57 mya, predating the Andean uplift, with some taxa showing dispersal patterns less constrained by ecological changes than previously thought.

Comparative display behaviour of the native Iguana delicatissima with the non-native Iguana in the Guadeloupe Archipelago (Lesser Antilles)

Territorial and mating displays are a key mechanism affecting sexual selection, species recognition, and may affect success of non-native species in the presence of closely related ones. The green Iguana (Iguana iguana) has invaded Caribbean islands from their mainland range, affecting and possibly inter-breeding with native Iguana. However, important display behaviours have not been studied, even though it may be a rare opportunity to study adaptation and evolutionary change. Here, we describe the display-action-patterns (DAP) characterising the head-bob display of the two Iguana species of the Caribbean islands Iguana delicatissima (endemic Lesser Antillean Iguana) and non-native Iguana (closely related non-native species) in the Guadeloupe Archipelago. First, the DAP sequence of non-native Iguana is completely different from mainland animals, with reversed positions of the long and short head bobs. In the non-native Iguana, two long head-bobs (>1.0 s) separated by a highly variable pause (0.1 - 2.8 s) are followed by a series of 3-7 (up to 9) rapid head bobs (<0.5 s), is followed by one phase always with two longer. In native I. delicatissima, the DAP is always initiated with a vibratory gular movement of the extended dewlap in I. delicatissima, but not in the non-native Iguana. Duration of the DAP sequence is significantly shorter in non-native Iguana, always beginning with one singular long head bob (>1.5 s) followed by 3-6 more rapid head bobs. In addition, we observed significant among individual variation in all aspects of the DAP studied, indicating the potential for natural selection to further act on head-bob displays in the non-native and native species. This study now sets the stage to for future studies to determine if behaviour is a cause or consequence of invasion success.

Tracing the path from conservation to expansion evolutionary insights into NLR genes in oleaceae

The Oleaceae family, encompassing key genera such as Fraxinus (ash trees), Olea (olives), Jasminum (jasmine), Syringa (lilac), and Forsythia, plays a crucial ecological and economic role. Despite their importance, the evolutionary dynamics and immune system adaptations of their NLR (Nucleotide binding leucine-rich repeats) gene family remain largely unexplored. This study employs high-throughput comparative genomics to investigate NLR gene evolution across the Oleaceae family. The genus Fraxinus is widely distributed across both the New and Old Worlds, with 23 distinct species analyzed in this study. Our results reveal a predominant strategy of gene conservation in the evolution of the NLR gene family across these species. Geographical adaptation has played a significant role, particularly in Old World ash tree species, which exhibit dynamic patterns of gene expansion and contraction within the last 50 million years. Notably, genes acquired from an ancient whole genome duplication event (~ 35 Mya) have been retained across Fraxinus lineages. In contrast, the genus Olea (olives) has undergone extensive gene expansion driven by recent duplications and significant birth of novel NLR gene families. These differences in NLR gene evolution likely enhance Olea's ability to recognize diverse pathogens through recent expansions, while Fraxinus maintains specialized immune responses through conserved genes, with potential trade-offs in pathogen adaptation and energy efficiency. In terms of NLR distribution, all species of the Oleaceae family show an enhanced pseudogenization of TIR-NLRs and expansion in CCG10-NLR. However, the comparative RNA-seq expression analysis in olive suggests that partial NLR genes, despite their incomplete structure, have significant expression and may play important roles in plant immune responses. This study provides a comprehensive analysis of NLR gene evolution within the Oleaceae family, offering insights into the adaptive mechanisms of immune response evolution across diverse genera.

Five-factor theory of aging and death due to aging

This new theory of aging explains that aging and death due to aging are due to five factors, and also explains how these factors are interconnected and jointly lead to aging and death of the organism, pointing to many facts that strongly support it. The first factor is the harmful changes that occur in cellular structures. The second factor is the cessation of cell division in adult organisms, which leads to the inability to restore cellular structures. The third factor is the feature that cells do not die due to the accumulation of harmful changes that occur in the cells during the life of the organism. The fourth factor is the inability of stem cells to regenerate tissue by replacing such cells with new ones, because somatic cells do not die and there are no signals that stimulate the proliferation of stem cells and their differentiation into new ones that would replace dead cells. The fifth factor is that all cells die suddenly, due to the cessation of one of the vital functions of the organism, and not gradually during life, due to a decrease in the functionality of cells caused by the introduction of harmful changes in cellular structures, which would allow stem cells to regenerate tissues and keep the body young. Also, to show that this aging theory is valid, the theory gives its view of the evolution of five factors, which according to this theory lead to aging, which gives strong support to this theory.

Plant interactions associated with a directional shift in the richness range size relationship during the Glacial-Holocene transition in the Arctic

A nearly ubiquitous negative relationship between taxonomic richness and mean range-size (average area of taxa) is observed across space. However, the complexity of the mechanism limits its applicability for conservation or range prediction. We explore whether the relationship holds over time, and whether plant speciation, environmental heterogeneity, or plant interactions are major factors of the relationship within northeast Siberia and Alaska. By analysing sedimentary ancient DNA from seven lakes, we reconstruct plant richness, biotic environmental heterogeneity, and mean range-size over the last 30,000 years. We find positive richness to range-size relationships during the glacial period, shifting to negative during the interglacial period. Our results indicate neither speciation nor environmental heterogeneity is the principal driver. Network analyses show more positive interactions during the glacial period, which may contribute to positive richness to range-size relationships. Conversely, in the interglacial environment, negative interactions may result in negative relationships. Our findings suggest potential susceptibility to invasion but conservation advantages in far northern tundra given their positive interactions.

Chloroplast Genome Diversity and Molecular Evolution in Hypericaceae: New Insights from Three Hypericum Species

The Hypericaceae family, comprising nine genera and over seven hundred species, includes Hypericum plants traditionally used for medicinal purposes. In this study, we performed high-throughput sequencing on three Hypericum species: Hypericum acmosepalum, Hypericum addingtonii, and Hypericum beanii, and conducted comparative genomic analyses with related species. The chloroplast genome sizes were 152,654 bp, 122,570 bp, and 137,652 bp, respectively, with an average GC content of 37.9%. All genomes showed a quadripartite structure, with significant variations in IR regions (3231-26,846 bp). The total number of genes ranged from 91 to 129. SSRs were predominantly located in the LSC region, with mononucleotide repeats being dominant. Comparative analysis identified several hotspot regions, including accD, rpoC2, rpoB, and rpl22 in the LSC region and matK, rpl32, rpl33, and rps4 in the SSC region. Nucleotide polymorphism analysis revealed eight highly variable regions and eleven gene loci, providing potential molecular markers for species identification. Phylogenetic analysis indicated that Triadenum and Cratoxylum are closely related to Hypericum, with H. acmosepalum and H. beanii being closest relatives and Hypericum hookerianum as their sister species. These findings provide molecular tools for species identification and insights for conservation strategies of medicinal Hypericum species.

A new forest lizard fauna from Caribbean islands (Squamata, Diploglossidae, Celestinae)

The taxonomy of Neotropical forest lizards (Diploglossidae) has been the subject of controversy because of a paucity of diagnostic characters and genetic data. Recent molecular studies have produced phylogenies that are highly supported but have few individuals represented for each species. These studies have corrected generic names and defined new genera and subfamilies in Diploglossidae. However, they have shown that multiple species are not monophyletic or have high levels of genetic divergence, indicating the need to define new species. Three subfamilies, 12 genera, and 56 species of diploglossid lizards are currently recognized; 25 of these are in the subfamily Celestinae. We conducted a systematic revision of Caribbean celestine forest lizards (from the Cayman Islands, Jamaica, and Hispaniola) using DNA sequence data from 372 individuals, supplemented by both conventional and unconventional morphological characters from 958 preserved specimens. In some cases, we obtained DNA sequence data from museum specimens, including types, nearly 200 years old. We propose and use a new species delimitation method based on time of divergence. We define 17 new species, elevate 17 subspecies, and elevate one species from synonymy, resulting in 35 newly recognized species. Additionally, we synonymize two pairs of previously recognized subspecies and one pair of species. This increases the number of celestine species from 25 to 59 and raises the total number of diploglossids to 90 species. Of those, 63 occur on Caribbean islands and all are endemic to those islands. Fourteen Caribbean celestine species (24%) are Critically Endangered, 17 species (29%) are Endangered, and 1 species (2%) is Vulnerable, resulting in a proportion of threatened species (54%) more than twice as high as the average for reptiles, based on IUCN Redlist criteria. Three of the Critically Endangered species are possibly extinct because of human activities during the last two centuries. Several of the surviving species are near extinction and in need of immediate protection. Extensive forest loss on Caribbean islands has led to the decline of Caribbean forest lizards, which rely on forests as their primary habitat. In addition to deforestation, the introduction of the Small Indian Mongoose is in part responsible for the decline of Caribbean diploglossid lizards. That invasive predator was introduced as a biological control of rats in sugar cane fields in the late 19th Century (1872-1900), immediately resulting in a mass extinction of reptiles. The ground-dwelling and diurnal habits of diploglossids have made them particularly susceptible to mongoose predation.

Completing a molecular timetree of primates

Primates, consisting of apes, monkeys, tarsiers, and lemurs, are among the most charismatic and well-studied animals on Earth, yet there is no taxonomically complete molecular timetree for the group. Combining the latest large-scale genomic primate phylogeny of 205 recognized species with the 400-species literature consensus tree available from TimeTree.org yields a phylogeny of just 405 primates, with 50 species still missing despite having molecular sequence data in the NCBI GenBank. In this study, we assemble a timetree of 455 primates, incorporating every species for which molecular data are available. We use a synthetic approach consisting of a literature review for published timetrees, de novo dating of untimed trees, and assembly of timetrees from novel alignments. The resulting near-complete molecular timetree of primates allows testing of two long-standing alternate hypotheses for the origins of primate biodiversity: whether species richness arises at a constant rate, in which case older clades have more species, or whether some clades exhibit faster rates of speciation than others, in which case, these fast clades would be more species-rich. Consistent with other large-scale macroevolutionary analyses, we found that the speciation rate is similar across the primate tree of life, albeit with some variation in smaller clades.

The last of their kind: Is the genus Scutiger (Anura: Megophryidae) a relict element of the paleo-Transhimalaya biota?

The orographic evolution of the Himalaya-Tibet Mountain system continues to be a subject of controversy, leading to considerable uncertainty regarding the environment and surface elevation of the Tibetan Plateau during the Cenozoic era. As many geoscientific (but not paleontological) studies suggest, elevations close to modern heights exist in vast areas of Tibet since at least the late Paleogene, implicating the presence of large-scale alpine environments for more than 30 million years. To explore a recently proposed alternative model that assumes a warm temperate environment across paleo-Tibet, we carried out a phylogeographic survey using genomic analyses of samples covering the range of endemic lazy toads (Scutiger) across the Himalaya-Tibet orogen. We identified two main clades, with several, geographically distinct subclades. The long temporal gap between the stem and crown age of Scutiger may suggest high extinction rates. Diversification within the crown group, depending on the calibration, occurred either from the Mid-Miocene or Late-Miocene and continued until the Holocene. The present-day Himalayan Scutiger fauna could have evolved from lineages that existed on the southern edges of the paleo-Tibetan area (the Transhimalaya = Gangdese Shan), while extant species living on the eastern edge of the Plateau originated probably from the eastern edges of northern parts of the ancestral Tibetan area (Hoh Xil, Tanggula Shan). Based on the Mid-Miocene divergence time estimation and ancestral area reconstruction, we propose that uplift-associated aridification of a warm temperate Miocene-Tibet, coupled with high extirpation rates of ancestral populations, and species range shifts along drainage systems and epigenetic transverse valleys of the rising mountains, is a plausible scenario explaining the phylogenetic structure of Scutiger. This hypothesis aligns with the fossil record but conflicts with geoscientific concepts of high elevated Tibetan Plateau since the late Paleogene. Considering a Late-Miocene/Pliocene divergence time, an alternative scenario of dispersal from SE Asia into the East, Central, and West Himalaya cannot be excluded, although essential evolutionary and biogeographic aspects remain unresolved within this model.

The early transcriptional and post-transcriptional responses to fluconazole in sensitive and resistant Candida albicans

Candida albicans is a leading cause of fungal infections in immunocompromised patients. Management of candidemia relies on a few antifungal agents, with fluconazole being first line therapy. The emergence of fluconazole-resistant strains highlights the pressing need to improve our molecular understanding of the drug response mechanisms. By sequencing the 5'P mRNA degradation intermediates, we establish that co-translational mRNA decay occurs in C. albicans and characterize how in vivo 5´-3´ exonuclease degradation trails the last translating ribosome. Thus, the study of the 5' Phosphorylated mRNA degradome (5PSeq) offers a simple and affordable way to measure ribosome dynamics and identify codon specific ribosome stalls in response to drugs and amino acid deprivation. Building upon this, we combine RNA-Seq and 5PSeq to study the early response of sensitive and resistant C. albicans isolates to fluconazole. Our results highlight that transcriptional responses, rather than changes in ribosome dynamics, are the main driver of Candida resistance to fluconazole.

Shell Constraints on Evolutionary Body Size-Limb Size Allometry Can Explain Morphological Conservatism in the Turtle Body Plan

Turtles are a small clade of vertebrates despite having existed since the Late Triassic. Turtles have a conservative body plan relative to other amniotes, characterized by the presence of a shell and quadrupedality. This morphology is even retained in strong ecological specialists, such as sea turtles, which are secondarily adapted to marine locomotion by strong allometric scaling in their hands. It is possible that the body plan of turtles is strongly influenced by the presence of the shell, acting as a constraint to achieving greater diversity of body forms. Here, we explore the evolutionary allometric relationships of fore- and hindlimb stylopodia (i.e., humerus and femur) with one another as well as their relationship with shell size (carapace length) to assess evidence of constraint. All turtles, including Triassic shelled stem turtles, have near-isometric relationships that do not vary strongly between clades, and evolve at slow evolutionary rates. This indeed indicates that body proportions of turtles are constrained to a narrow range of possibilities. Minor allometric deviations are seen in highly aquatic sea turtles and softshell turtles, which modified their shells by bone losses. Our allometric regressions allow accurate body size estimations for fossils. Several independent sea turtle lineages converged on maximum sizes of 2.2 m of shell length, which may be a biological maximum for the group.

Inferring the locomotor ecology of two of the oldest fossil squirrels: influence of operationalization, trait, body size and machine learning method

Correlations between morphology and lifestyle of extant taxa are useful for predicting lifestyles of extinct relatives. Here, we infer the locomotor behaviour of Palaeosciurus goti from the middle Oligocene and Palaeosciurus feignouxi from the lower Miocene of France using their femoral morphology and different machine learning methods. We used two ways to operationalize morphology, in the form of a geometric morphometric shape dataset and a multivariate dataset of 11 femoral traits. The predictive models were built and tested using more than half (180) of the extant species of squirrel relatives. Both traditional models such as linear discriminant analysis and more sophisticated models like neural networks had the greatest predictive power. However, the predictive power also depended on the operationalization and the femoral traits used to build the model. We also found that predictive power tended to improve with increasing body size. Contrary to previous suggestions, the older species, P. goti, was most likely arboreal, whereas P. feignouxi was more likely terrestrial. This provides further evidence that arboreality was already the most common locomotor ecology among the earliest squirrels, while a predominantly terrestrial locomotor behaviour evolved shortly afterwards, before the vast establishment of grasslands in Europe.

Image-Based Quantitative Analysis of Epidermal Morphology in Wild Potato Leaves

The epidermal leaf patterns of plants exhibit remarkable diversity in cell shapes, sizes, and arrangements, driven by environmental interactions that lead to significant adaptive changes even among closely related species. The Solanaceae family, known for its high diversity of adaptive epidermal structures, has traditionally been studied using qualitative phenotypic descriptions. To advance this, we developed a workflow combining multi-scale computer vision, image processing, and data analysis to extract digital descriptors for leaf epidermal cell morphology. Applied to nine wild potato species, this workflow quantified key morphological parameters, identifying descriptors for trichomes, stomata, and pavement cells, and revealing interdependencies among these traits. Principal component analysis (PCA) highlighted two main axes, accounting for 45% and 21% of variance, corresponding to features such as guard cell shape, trichome length, stomatal density, and trichome density. These axes aligned well with the historical and geographical origins of the species, separating southern from Central American species, and forming distinct clusters for monophyletic groups. This workflow thus establishes a quantitative foundation for investigating leaf epidermal cell morphology within phylogenetic and geographic contexts.

Intrinsic structural disorder on proteins is involved in the interactome evolution

New mathematical tools help understand cell functions, adaptability, and evolvability to discover hidden variables to predict phenotypes that could be tested in the future in wet labs. Different models have been successfully used to discover the properties of the protein-protein interaction networks or interactomes. I found that in the hyperbolic Popularity-Similarity model, cellular proteins with the highest contents of structural intrinsic disorder cluster together in many different eukaryotic interactomes and this is not the case for the prokaryotic E. coli, where proteins with high degree of intrinsic disorder are scarce. I also found that the normalized theta variable from the Popularity-Similarity model for orthologues proteins correlate to the complexity of the organisms in analysis.

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.

Functional and regulatory diversification of Period genes responsible for circadian rhythm in vertebrates

The Period genes (Per) play essential roles in modulating the molecular circadian clock timing in a broad range of species, which regulates the physiological and cellular rhythms through the transcription-translation feedback loop. While the Period gene paralogs are widely observed among vertebrates, the evolutionary history and the functional diversification of Per genes across vertebrates are not well known. In this study, we comprehensively investigated the evolution of Per genes at the copy number and sequence levels, including de novo binding motif discovery by comparative genomics. We also determined the lineage-specific transcriptome landscape across tissues and developmental stages and phenotypic effects in public RNA-seq data sets of model species. We observed multiple lineage-specific gain and loss events Per genes, though no simple association was observed between ecological factors and Per gene numbers in each species. Among salmonid fish species, the per3 gene has been lost in the majority, whereas those retaining the per3 gene exhibit not a signature of relaxed selective constraint but rather a signature of intensified selection. We also determined the signature of adaptive diversification of the CRY-binding region in Per1 and Per3, which modulates the circadian rhythm. We also discovered putative regulatory sequences, which are lineage-specific, suggesting that these cis-regulatory elements may have evolved rapidly and divergently across different lineages. Collectively, our findings revealed the evolution of Per genes and their fine-tuned contribution to the plastic and precise regulation of circadian rhythms in various vertebrate taxa.

A time-calibrated phylogeny of the diversification of Holoadeninae frogs

The phylogeny of the major lineages of Amphibia has received significant attention in recent years, although evolutionary relationships within families remain largely neglected. One such overlooked group is the subfamily Holoadeninae, comprising 73 species across nine genera and characterized by a disjunct geographical distribution. The lack of a fossil record for this subfamily hampers the formulation of a comprehensive evolutionary hypothesis for their diversification. Aiming to fill this gap, we inferred the phylogenetic relationships and divergence times for Holoadeninae using molecular data and calibration information derived from the fossil record of Neobatrachia. Our inferred phylogeny confirmed most genus-level associations, and molecular dating analysis placed the origin of Holoadeninae in the Eocene, with subsequent splits also occurring during this period. The climatic and geological events that occurred during the Oligocene-Miocene transition were crucial to the dynamic biogeographical history of the subfamily. However, the wide highest posterior density intervals in our divergence time estimates are primarily attributed to the absence of Holoadeninae fossil information and, secondarily, to the limited number of sampled nucleotide sites.

Retrospect and prospect of Nicotiana tabacum genome sequencing

Investigating plant genomes offers crucial foundational resources for exploring various aspects of plant biology and applications, such as functional genomics and breeding practices. With the development in sequencing and assembly technology, several Nicotiana tabacum genomes have been published. In this paper, we reviewed the progress on N. tabacum genome assembly and quality, from the initial draft genomes to the recent high-quality chromosome-level assemblies. The application of long-read sequencing, optical mapping, and Hi-C technologies has significantly improved the contiguity and completeness of N. tabacum genome assemblies, with the latest assemblies having a contig N50 size over 50 Mb. Despite these advancements, further improvements are still required and possible, particularly on the development of pan-genome and telomere-to-telomere (T2T) genomes. These new genomes will capture the genomic diversity and variations among different N. tabacum cultivars and species, and provide a comprehensive view of the N. tabacum genome structure and gene content, so to deepen our understanding of the N. tabacum genome and facilitate precise breeding and functional genomics.

Constraint on boric acid resistance and tolerance evolvability in Candida albicans

Boric acid is a broad-spectrum antimicrobial used to treat vulvovaginal candidiasis when patients relapse on the primary azole drug fluconazole. Candida albicans is the most common cause of vulvovaginal candidiasis, colloquially referred to as a "vaginal yeast infection". Little is known about the propensity of C. albicans to develop BA resistance or tolerance (the ability of a subpopulation to grow slowly in high levels of drug). We evolved 96 replicates from eight diverse C. albicans strains to increasing BA concentrations to test the evolvability of BA resistance and tolerance. Replicate growth was individually assessed daily, with replicates passaged when they had reached an optical density consistent with exponential growth. Many replicates went extinct quickly. Although some replicates could grow in much higher levels of BA than the ancestral strains, evolved populations isolated from the highest terminal BA levels (after 11 weeks of passages) surprisingly showed only modest growth improvements and only at low levels of BA. No large increases in resistance or tolerance were observed in the evolved replicates. Overall, our findings illustrate that there may be evolutionary constraints limiting the emergence of BA resistance and tolerance, which could explain why it remains an effective treatment for recurrent yeast infections.

Chromosome-scale genome assembly reveals insights into the evolution and ecology of the harmful algal bloom species Phaeocystis globosa Scherffel

The phytoplankton Phaeocystis globosa plays an important role in sulfur cycling and climate control, and can develop harmful algal blooms (HABs). Here we report a chromosome-scale reference genome assembly of P. globosa, which enable in-depth analysis of molecular underpinnings of important ecological characteristics. Comparative genomic analyses detected two-rounds of genome duplications that may have fueled evolutionary innovations. The genome duplication may have resulted in the formation of dual HiDP and LoDP dimethylsulphoniopropionate (DMSP) biosynthesis pathways in P. globosa. Selective gene family expansions may have strengthened biological pathways critical for colonial formation that is often associated with the development of algal blooms. The copy numbers of rhodopsin genes are variable in different strains, suggesting that rhodopsin genes may play a role in strain-specific adaptation to ecological factors. The successful reconstruction of the P. globosa genome sets up an excellent platform that facilitates in-depth research on bloom development and DMSP metabolism.

The making of a nucleic acid sensor at the dawn of jawed vertebrate evolution

Self and nonself discrimination is fundamental to immunity. However, it remains largely enigmatic how the mechanisms of distinguishing nonself from self originated. As an intracellular nucleic acid sensor, protein kinase R (PKR) recognizes double-stranded RNA (dsRNA) and represents a crucial component of antiviral innate immunity. Here, we combine phylogenomic and functional analyses to show that PKR proteins probably originated from a preexisting kinase protein through acquiring dsRNA binding domains at least before the last common ancestor of jawed vertebrates during or before the Silurian period. The function of PKR appears to be conserved across jawed vertebrates. Moreover, we repurpose a protein closely related to PKR proteins into a putative dsRNA sensor, recapturing the making of PKR. Our study illustrates how a nucleic acid sensor might have originated via molecular tinkering with preexisting proteins and provides insights into the origins of innate immunity.

Exploring Paleogene Tibet's warm temperate environments through target enrichment and phylogenetic niche modelling of Himalayan spiny frogs (Paini, Dicroglossidae)

The Cenozoic topographic development of the Himalaya-Tibet orogen (HTO) substantially affected the paleoenvironment and biodiversity patterns of High Asia. However, concepts on the evolution and paleoenvironmental history of the HTO differ massively in timing, elevational increase and sequence of surface uplift of the different elements of the orogen. Using target enrichment of a large set of transcriptome-derived markers, ancestral range estimation and paleoclimatic niche modelling, we assess a recently proposed concept of a warm temperate paleo-Tibet in Asian spiny frogs of the tribe Paini and reconstruct their historical biogeography. That concept was previously developed in invertebrates. Because of their early evolutionary origin, low dispersal capacity, high degree of local endemism, and strict dependence on temperature and humidity, the cladogenesis of spiny frogs may echo the evolution of the HTO paleoenvironment. We show that diversification of main lineages occurred during the early to Mid-Miocene, while the evolution of alpine taxa started during the late Miocene/early Pliocene. Our distribution and niche modelling results indicate range shifts and niche stability that may explain the modern disjunct distributions of spiny frogs. They probably maintained their (sub)tropical or (warm)temperate preferences and moved out of the ancestral paleo-Tibetan area into the Himalaya as the climate shifted, as opposed to adapting in situ. Based on ancestral range estimation, we assume the existence of low-elevation, climatically suitable corridors across paleo-Tibet during the Miocene along the Kunlun, Qiangtang and/or Gangdese Shan. Our results contribute to a deeper understanding of the mechanisms and processes of faunal evolution in the HTO.

Freshwater fish and the Cretaceous/Palaeogene boundary: a critical assessment of survivorship patterns

Mass extinctions are major influences on both the phylogenetic structure of the modern biota and our ability to reconstruct broad-based patterns of evolutionary history. The most recent mass extinction is also the most famous-that which implicates a bolide impact in defining the Cretaceous/Palaeogene boundary (K/Pg). Although the biotic effects of this event receive intensive scrutiny, certain ecologically important and diverse groups remain woefully understudied. One such group is the freshwater ray-finned fishes (Actinopterygii). These fish represent 25% of modern vertebrate diversity, yet the isolated and fragmentary nature of their K/Pg fossil record limits our understanding of their diversity dynamics across this event. Here, we address this problem using diversification analysis of molecular-based phylogenies alongside a morphotype analysis of fossils recovered from a unique site in the Denver Basin of western North America that provides unprecedented K/Pg resolution. Our results reveal previously unrecognized signals of post-K/Pg diversification in freshwater clades and suggest that the change was driven by localized and sporadic patterns of extinction. Supported inferences regarding the effects of the K/Pg event on freshwater fish also inform our expectations of how freshwater faunas might recover from the current biodiversity crisis.

DateLife: Leveraging Databases and Analytical Tools to Reveal the Dated Tree of Life

Chronograms-phylogenies with branch lengths proportional to time-represent key data on timing of evolutionary events, allowing us to study natural processes in many areas of biological research. Chronograms also provide valuable information that can be used for education, science communication, and conservation policy decisions. Yet, achieving a high-quality reconstruction of a chronogram is a difficult and resource-consuming task. Here we present DateLife, a phylogenetic software implemented as an R package and an R Shiny web application available at www.datelife.org, that provides services for efficient and easy discovery, summary, reuse, and reanalysis of node age data mined from a curated database of expert, peer-reviewed, and openly available chronograms. The main DateLife workflow starts with one or more scientific taxon names provided by a user. Names are processed and standardized to a unified taxonomy, allowing DateLife to run a name match across its local chronogram database that is curated from Open Tree of Life's phylogenetic repository, and extract all chronograms that contain at least two queried taxon names, along with their metadata. Finally, node ages from matching chronograms are mapped using the congruification algorithm to corresponding nodes on a tree topology, either extracted from Open Tree of Life's synthetic phylogeny or one provided by the user. Congruified node ages are used as secondary calibrations to date the chosen topology, with or without initial branch lengths, using different phylogenetic dating methods such as BLADJ, treePL, PATHd8, and MrBayes. We performed a cross-validation test to compare node ages resulting from a DateLife analysis (i.e, phylogenetic dating using secondary calibrations) to those from the original chronograms (i.e, obtained with primary calibrations), and found that DateLife's node age estimates are consistent with the age estimates from the original chronograms, with the largest variation in ages occurring around topologically deeper nodes. Because the results from any software for scientific analysis can only be as good as the data used as input, we highlight the importance of considering the results of a DateLife analysis in the context of the input chronograms. DateLife can help to increase awareness of the existing disparities among alternative hypotheses of dates for the same diversification events, and to support exploration of the effect of alternative chronogram hypotheses on downstream analyses, providing a framework for a more informed interpretation of evolutionary results.

Tetrapod sperm length evolution in relation to body mass is shaped by multiple trade-offs

Sperm length is highly variable across species and many questions about its variation remain open. Although variation in body mass may affect sperm length evolution through its influence on multiple factors, the extent to which sperm length variation is linked to body mass remains elusive. Here, we use the Pareto multi-task evolution framework to investigate the relationship between sperm length and body mass across tetrapods. We find that tetrapods occupy a triangular Pareto front, indicating that trade-offs shape the evolution of sperm length in relation to body mass. By exploring the factors predicted to influence sperm length evolution, we find that sperm length evolution is mainly driven by sperm competition and clutch size, rather than by genome size. Moreover, the triangular Pareto front is maintained within endotherms, internal fertilizers, mammals and birds, suggesting similar evolutionary trade-offs within tetrapods. Finally, we demonstrate that the Pareto front is robust to phylogenetic dependencies and finite sampling bias. Our findings provide insights into the evolutionary mechanisms driving interspecific sperm length variation and highlight the importance of considering multiple trade-offs in optimizing reproductive traits.

Integrating Phylogenies with Chronology to Assemble the Tree of Life

Reconstructing the global Tree of Life necessitates computational approaches to integrate numerous molecular phylogenies with limited species overlap into a comprehensive supertree. Our survey of published literature shows that individual phylogenies are frequently restricted to specific taxonomic groups due to the expertise of investigators and molecular evolutionary considerations, resulting in any given species present in a minuscule fraction of phylogenies. We present a novel approach, called the chronological supertree algorithm (Chrono-STA), that can build a supertree of species from such data by using node ages in published molecular phylogenies scaled to time. Chrono-STA builds a supertree of organisms by integrating chronological data from molecular timetrees. It fundamentally differs from existing approaches that generate consensus phylogenies from gene trees with missing taxa, as Chrono-STA does not impute nodal distances, use a guide tree as a backbone, or reduce phylogenies to quartets. Analyses of simulated and empirical datasets show that Chrono-STA can combine taxonomically restricted timetrees with extremely limited species overlap. For such data, approaches that impute missing distances or assemble phylogenetic quartets did not perform well. We conclude that integrating phylogenies via temporal dimension enhances the accuracy of reconstructed supertrees that are also scaled to time.

Convergent evolution toward a slow pace of life predisposes insular endotherms to anthropogenic extinctions

Island vertebrates have evolved a number of morphological, physiological, and life history characteristics that set them apart from their mainland relatives. However, to date, the evolution of metabolism and its impact on the vulnerability to extinction of insular vertebrates remains poorly understood. This study used metabolic data from 2813 species of tetrapod vertebrates, including 695 ectothermic and 2118 endothermic species, to reveal that island mammals and birds evolved convergent metabolic strategies toward a slow pace of life. Insularity was associated with shifts toward slower metabolic rates and greater generation lengths in endotherms, while insularity just drove the evolution of longer generation lengths in ectotherms. Notably, a slow pace of life has exacerbated the extinction of insular endemic species in the face of anthropogenic threats. These findings have important implications for understanding physiological adaptations associated with the island syndrome and formulating conservation strategies across taxonomic groups with different metabolic modes.

Compound osteoderms preserved in amber reveal the oldest known skink

Scincidae is one of the most species-rich and cosmopolitan clades of squamate reptiles. Abundant disarticulated fossil material has also been attributed to this group, however, no complete pre-Cenozoic crown-scincid specimens have been found. A specimen in Burmite (99 MYA) is the first fossil that can be unambiguously referred to this clade. Our analyses place it as nested within extant skinks, supported by the presence of compound osteoderms formed by articulated small ostedermites. The specimen has a combination of dorsal and ventral compound osteoderms and overlapping cycloid scales that is limited to skinks. We propose that this type of osteoderm evolved as a response to an increased overlap of scales, and to reduced stiffness of the dermal armour. Compound osteoderms could be a key innovation that facilitated diversification in this megadiverse family.

Identifying genes within pathways in unannotated genomes with PaGeSearch

In biological research, the identification and comparison of genes within specific pathways across the genomes of various species are invaluable. However, annotating the entire genome is resource intensive, and sequence similarity searches often yield results that are not actually genes. To address these limitations, we introduce Pathway Gene Search (PaGeSearch), a tool designed to identify genes from predefined lists, especially those in specific pathways, within genomes. The tool uses an initial sequence similarity search to identify relevant genomic regions, followed by targeted gene prediction and neural network-based result filtering. PaGeSearch suggests the regions that are most likely the orthologs of the genes in the query and is designed to be applicable for species within five classes: mammals, fish, birds, eudicotyledons, and Liliopsida. Compared with GeMoMa and miniprot, PaGeSearch generally outperforms in terms of sensitivity and positive predictive value, as well as negative predictive value. Also, the exon coverage of gene models from PaGeSearch is higher compared with those in GeMoMa and miniprot. Although its performance shows increased variability when applied to actual biological pathways, it nonetheless maintains an acceptable level of accuracy. Evaluating PaGeSearch across different assembly levels, chromosome, scaffold, and contig shows minimal variation in outcomes, indicating that PaGeSearch is resilient to variations in assembly quality.

Evolution and adaptation of terrestrial plant-associated Plantibacter species into remote marine environments

Microbes are thought to be distributed and circulated around the world, but the connection between marine and terrestrial microbiomes remains largely unknown. We use Plantibacter, a representative genus associated with plants, as our research model to investigate the global distribution and adaptation of plant-related bacteria in plant-free environments, particularly in the remote Southern Ocean and the deep Atlantic Ocean. The marine isolates and their plant-associated relatives shared over 98% whole-genome average nucleotide identity (ANI), indicating recent divergence and ongoing speciation from plant-related niches to marine environments. Comparative genomics revealed that the marine strains acquired new genes via horizontal gene transfer from non-Plantibacter species and refined existing genes through positive selection to improve adaptation to new habitats. Meanwhile, marine strains retained the ability to interact with plants, such as modifying root system architecture and promoting germination. Furthermore, Plantibacter species were found to be widely distributed in marine environments, revealing an unrecognized phenomenon that plant-associated microbiomes have colonized the ocean, which could serve as a reservoir for plant growth-promoting microbes. This study demonstrates the presence of an active reservoir of terrestrial plant growth-promoting bacteria in remote marine systems and advances our understanding of the microbial connections between plant-associated and plant-free environments at the genome level.

Geographical patterns and determinants of insect biodiversity in China

Insects play important roles in the maintenance of ecosystem functioning and the provision of livelihoods for millions of people. However, compared with terrestrial vertebrates and angiosperms, such as the giant panda, crested ibis, and the metasequoia, insect conservation has not attracted enough attention, and a basic understanding of the geographical biodiversity patterns for major components of insects in China is lacking. Herein, we investigated the geographical distribution of insect biodiversity across multiple dimensions (taxonomic, genetic, and phylogenetic diversity) based on the spatial distribution and molecular DNA sequencing data of insects. Our analysis included 18 orders, 360 families, 5,275 genera, and 14,115 species of insects. The results revealed that Southwestern and Southeastern China harbored higher insect biodiversity and numerous older lineages, representing a museum, whereas regions located in Northwestern China harbored lower insect biodiversity and younger lineages, serving as an evolutionary cradle. We also observed that mean annual temperature and precipitation had significantly positive effects, whereas altitude had significantly negative effects on insect biodiversity in most cases. Moreover, cultivated vegetation harbored the highest insect taxonomic and phylogenetic diversity, and needleleaf and broadleaf mixed forests harbored the highest insect genetic diversity. These results indicated that human activities may positively contribute to insect spatial diversity on a regional scale. Our study fills a knowledge gap in insect spatial diversity in China. These findings could help guide national-level conservation plans and the post-2020 biodiversity conservation framework.

Bacterial sepsis triggers stronger transcriptomic immune responses in larger primates

Empirical data relating body mass to immune defence against infections remain limited. Although the metabolic theory of ecology predicts that larger organisms would have weaker immune responses, recent studies have suggested that the opposite may be true. These discoveries have led to the safety factor hypothesis, which proposes that larger organisms have evolved stronger immune defences because they carry greater risks of exposure to pathogens and parasites. In this study, we simulated sepsis by exposing blood from nine primate species to a bacterial lipopolysaccharide (LPS), measured the relative expression of immune and other genes using RNAseq, and fitted phylogenetic models to determine how gene expression was related to body mass. In contrast to non-immune-annotated genes, we discovered hypermetric scaling in the LPS-induced expression of innate immune genes, such that large primates had a disproportionately greater increase in gene expression of immune genes compared to small primates. Hypermetric immune gene expression appears to support the safety factor hypothesis, though this pattern may represent a balanced evolutionary mechanism to compensate for lower per-transcript immunological effectiveness. This study contributes to the growing body of immune allometry research, highlighting its importance in understanding the complex interplay between body size and immunity over evolutionary timescales.

Complete Chloroplast Genome of Alternanthera sessilis and Comparative Analysis with Its Congeneric Invasive Weed Alternanthera philoxeroides

Alternanthera sessilis is considered the closest relative to the invasive weed Alternanthera philoxeroides in China, making it an important native species for studying the invasive mechanisms and adaptations of A. philoxeroides. Chloroplasts play a crucial role in a plant's environmental adaptation, with their genomes being pivotal in the evolution and adaptation of both invasive and related species. However, the chloroplast genome of A. sessilis has remained unknown until now. In this study, we sequenced and assembled the complete chloroplast genome of A. sessilis using high-throughput sequencing. The A. sessilis chloroplast genome is 151,935 base pairs long, comprising two inverted repeat regions, a large single copy region, and a small single copy region. This chloroplast genome contains 128 genes, including 8 rRNA-coding genes, 37 tRNA-coding genes, 4 pseudogenes, and 83 protein-coding genes. When compared to the chloroplast genome of the invasive weed A. philoxeroides and other Amaranthaceae species, we observed significant variations in the ccsA, ycf1, and ycf2 regions in the A. sessilis chloroplast genome. Moreover, two genes, ccsA and accD, were found to be undergoing rapid evolution due to positive selection pressure. The phylogenetic trees were constructed for the Amaranthaceae family, estimating the time of independent species formation between A. philoxeroides and A. sessilis to be approximately 3.5186-8.8242 million years ago. These findings provide a foundation for understanding the population variation within invasive species among the Alternanthera genus.

Human gene regulatory evolution is driven by the divergence of regulatory element function in both cis and trans

Gene regulatory divergence between species can result from cis-acting local changes to regulatory element DNA sequences or global trans-acting changes to the regulatory environment. Understanding how these mechanisms drive regulatory evolution has been limited by challenges in identifying trans-acting changes. We present a comprehensive approach to directly identify cis- and trans-divergent regulatory elements between human and rhesus macaque lymphoblastoid cells using assay for transposase-accessible chromatin coupled to self-transcribing active regulatory region (ATAC-STARR) sequencing. In addition to thousands of cis changes, we discover an unexpected number (∼10,000) of trans changes and show that cis and trans elements exhibit distinct patterns of sequence divergence and function. We further identify differentially expressed transcription factors that underlie ∼37% of trans differences and trace how cis changes can produce cascades of trans changes. Overall, we find that most divergent elements (67%) experienced changes in both cis and trans, revealing a substantial role for trans divergence-alone and together with cis changes-in regulatory differences between species.

Clogmia albipunctata (Williston, 1893) midgut physiology: pH control and functional relationship with Lower Diptera (nematoceran) especially with hematophagous species

Clogmia albipunctata (Williston, 1893) is a non-hematophagous insect belonging to the order Diptera, suborder Nematocera (Lower Diptera) and family Psychodidae. In the present work, we investigated how C. albipunctata control their midgut pH under different physiological conditions, comparing their midgut physiology with some nematoceran hematophagous species. The C. albipunctata midgut pH was measured after ingestion of sugar, protein and under the effect of the alkalinizing hormone released in the hemolymph of the hematophagous sand fly Lutzomyia longipalpis obtained just after a blood meal. The midgut pH of unfed or sugar-fed C. albipunctata is 5.5-6, and its midgut underwent alkalinization after protein ingestion or under treatment with hemolymph collected from blood fed L. longipalpis. These results suggested that in nematocerans, mechanisms for pH control seem shared between hematophagous and non-hematophagous species. This kind of pH control is convenient for successful blood digestion. The independent evolution of many hematophagous groups from the Lower Diptera suggests that characteristics involved in midgut pH control were already present in non-hematophagous species and represent a readiness for adaptation to this feeding mode.

Cross-species transmission and host range genes in poxviruses

The persistent epidemic of human mpox, caused by mpox virus (MPXV), raises concerns about the future spread of MPXV and other poxviruses. MPXV is a typical zoonotic virus which can infect human and cause smallpox-like symptoms. MPXV belongs to the Poxviridae family, which has a relatively broad host range from arthropods to vertebrates. Cross-species transmission of poxviruses among different hosts has been frequently reported and resulted in numerous epidemics. Poxviruses have a complex linear double-strand DNA genome that encodes hundreds of proteins. Genes related to the host range of poxvirus are called host range genes (HRGs). This review briefly introduces the taxonomy, phylogeny and hosts of poxviruses, and then comprehensively summarizes the current knowledge about the cross-species transmission of poxviruses. In particular, the HRGs of poxvirus are described and their impacts on viral host range are discussed in depth. We hope that this review will provide a comprehensive perspective about the current progress of researches on cross-species transmission and HRG variation of poxviruses, serving as a valuable reference for academic studies and disease control in the future.

A cyst-forming coccidian with large geographical range infecting forest and commensal rodents: Sarcocystis muricoelognathis sp. nov

BACKGROUND

The geographic distribution and host-parasite interaction networks of Sarcocystis spp. in small mammals in eastern Asia remain incompletely known.

METHODS

Experimental infections, morphological and molecular characterizations were used for discrimination of a new Sarcocystis species isolated from colubrid snakes and small mammals collected in Thailand, Borneo and China.

RESULTS

We identified a new species, Sarcocystis muricoelognathis sp. nov., that features a relatively wide geographic distribution and infects both commensal and forest-inhabiting intermediate hosts. Sarcocystis sporocysts collected from rat snakes (Coelognathus radiatus, C. flavolineatus) in Thailand induced development of sarcocysts in experimental SD rats showing a type 10a cyst wall ultrastructure that was identical with those found in Rattus norvegicus from China and the forest rat Maxomys whiteheadi in Borneo. Its cystozoites had equal sizes in all intermediate hosts and locations, while sporocysts and cystozoites were distinct from other Sarcocystis species. Partial 28S rRNA sequences of S. muricoelognathis from M. whiteheadi were largely identical to those from R. norvegicus in China but distinct from newly sequenced Sarcocystis zuoi. The phylogeny of the nuclear 18S rRNA gene placed S. muricoelognathis within the so-called S. zuoi complex, including Sarcocystis attenuati, S. kani, S. scandentiborneensis and S. zuoi, while the latter clustered with the new species. However, the phylogeny of the ITS1-region confirmed the distinction between S. muricoelognathis and S. zuoi. Moreover, all three gene trees suggested that an isolate previously addressed as S. zuoi from Thailand (KU341120) is conspecific with S. muricoelognathis. Partial mitochondrial cox1 sequences of S. muricoelognathis were almost identical with those from other members of the group suggesting a shared, recent ancestry. Additionally, we isolated two partial 28S rRNA Sarcocystis sequences from Low's squirrel Sundasciurus lowii that clustered with those of S. scandentiborneensis from treeshews.

CONCLUSIONS

Our results provide strong evidence of broad geographic distributions of rodent-associated Sarcocystis and host shifts between commensal and forest small mammal species, even if the known host associations remain likely only snapshots of the true associations.

Association of ARRDC3 and NFIA variants with bovine congestive heart failure in feedlot cattle

BACKGROUND

Bovine congestive heart failure (BCHF) has become increasingly prevalent among feedlot cattle in the Western Great Plains of North America with up to 7% mortality in affected herds. BCHF is an untreatable complex condition involving pulmonary hypertension that culminates in right ventricular failure and death. Genes associated with BCHF in feedlot cattle have not been previously identified. Our aim was to search for genomic regions associated with this disease.

METHODS

A retrospective, matched case-control design with 102 clinical BCHF cases and their unaffected pen mates was used in a genome-wide association study. Paired nominal data from approximately 560,000 filtered single nucleotide polymorphisms (SNPs) were analyzed with McNemar's test.

RESULTS

Two independent genomic regions were identified as having the most significant association with BCHF: the arrestin domain-containing protein 3 gene ( ARRDC3), and the nuclear factor IA gene ( NFIA, mid- p-values, 1x10 -8 and 2x10 -7, respectively). Animals with two copies of risk alleles at either gene were approximately eight-fold more likely to have BCHF than their matched pen mates with either one or zero risk alleles at both genes (CI 95 = 3-17). Further, animals with two copies of risk alleles at both genes were 28-fold more likely to have BCHF than all others ( p-value = 1×10 -7, CI 95 = 4-206). A missense variant in ARRDC3 (C182Y) represents a potential functional variant since the C182 codon is conserved among all other jawed vertebrate species observed. A two-SNP test with markers in both genes showed 29% of 273 BCHF cases had homozygous risk genotypes in both genes, compared to 2.5% in 198 similar unaffected feedlot cattle. This and other DNA tests may be useful for identifying feedlot animals with the highest risk for BCHF in the environments described here.

CONCLUSIONS

Although pathogenic roles for variants in the ARRDC3 and NFIA genes are unknown, their discovery facilitates classifying animals by genetic risk and allows cattle producers to make informed decisions for selective breeding and animal health management.

Association of ARRDC3 and NFIA variants with bovine congestive heart failure in feedlot cattle

Background

Bovine congestive heart failure (BCHF) has become increasingly prevalent among feedlot cattle in the Western Great Plains of North America with up to 7% mortality in affected herds. BCHF is an untreatable complex condition involving pulmonary hypertension that culminates in right ventricular failure and death. Genes associated with BCHF in feedlot cattle have not been previously identified. Our aim was to search for genomic regions associated with this disease.

Methods

A retrospective, matched case-control design with 102 clinical BCHF cases and their unaffected pen mates was used in a genome-wide association study. Paired nominal data from approximately 560,000 filtered single nucleotide polymorphisms (SNPs) were analyzed with McNemar's test.

Results

Two independent genomic regions were identified as having the most significant association with BCHF: the arrestin domain-containing protein 3 gene ( ARRDC3), and the nuclear factor IA gene ( NFIA, mid- p-values, 1x10 -8 and 2x10 -7, respectively). Animals with two copies of risk alleles at either gene were approximately eight-fold more likely to have BCHF than their matched pen mates with either one or zero risk alleles at both genes (CI 95 = 3-17). Further, animals with two copies of risk alleles at both genes were 28-fold more likely to have BCHF than all others ( p-value = 1×10 -7, CI 95 = 4-206). A missense variant in ARRDC3 (C182Y) represents a potential functional variant since the C182 codon is conserved among all other jawed vertebrate species observed. A two-SNP test with markers in both genes showed 29% of 273 BCHF cases had homozygous risk genotypes in both genes, compared to 2.5% in 198 similar unaffected feedlot cattle. This and other DNA tests may be useful for identifying feedlot animals with the highest risk for BCHF in the environments described here.

Conclusions

Although pathogenic roles for variants in the ARRDC3 and NFIA genes are unknown, their discovery facilitates classifying animals by genetic risk and allows cattle producers to make informed decisions for selective breeding and animal health management.

Ancient and Modern Genomes Reveal Microsatellites Maintain a Dynamic Equilibrium Through Deep Time

Microsatellites are widely used in population genetics, but their evolutionary dynamics remain poorly understood. It is unclear whether microsatellite loci drift in length over time. This is important because the mutation processes that underlie these important genetic markers are central to the evolutionary models that employ microsatellites. We identify more than 27 million microsatellites using a novel and unique dataset of modern and ancient Adélie penguin genomes along with data from 63 published chordate genomes. We investigate microsatellite evolutionary dynamics over 2 timescales: one based on Adélie penguin samples dating to ∼46.5 ka and the other dating to the diversification of chordates aged more than 500 Ma. We show that the process of microsatellite allele length evolution is at dynamic equilibrium; while there is length polymorphism among individuals, the length distribution for a given locus remains stable. Many microsatellites persist over very long timescales, particularly in exons and regulatory sequences. These often retain length variability, suggesting that they may play a role in maintaining phenotypic variation within populations.

Global conservation status of the jawed vertebrate Tree of Life

Human-driven extinction threatens entire lineages across the Tree of Life. Here we assess the conservation status of jawed vertebrate evolutionary history, using three policy-relevant approaches. First, we calculate an index of threat to overall evolutionary history, showing that we expect to lose 86-150 billion years (11-19%) of jawed vertebrate evolutionary history over the next 50-500 years. Second, we rank jawed vertebrate species by their EDGE scores to identify the highest priorities for species-focused conservation of evolutionary history, finding that chondrichthyans, ray-finned fish and testudines rank highest of all jawed vertebrates. Third, we assess the conservation status of jawed vertebrate families. We found that species within monotypic families are more likely to be threatened and more likely to be in decline than other species. We provide a baseline for the status of families at risk of extinction to catalyse conservation action. This work continues a trend of highlighting neglected groups-such as testudines, crocodylians, amphibians and chondrichthyans-as conservation priorities from a phylogenetic perspective.

Haplotype-resolved genome assembly provides insights into evolutionary history of the Actinidia arguta tetraploid

Actinidia arguta, known as hardy kiwifruit, is a widely cultivated species with distinct botanical characteristics such as small and smooth-fruited, rich in beneficial nutrients, rapid softening and tolerant to extremely low temperatures. It contains the most diverse ploidy types, including diploid, tetraploid, hexaploid, octoploid, and decaploid. Here we report a haplotype-resolved tetraploid genome (A. arguta cv. 'Longcheng No.2') containing four haplotypes, each with 40,859, 41,377, 39,833 and 39,222 protein-coding genes. We described the phased genome structure, synteny, and evolutionary analyses to identify and date possible WGD events. Ks calculations for both allelic and paralogous genes pairs throughout the assembled haplotypic individuals showed its tetraploidization is estimated to have formed ~ 1.03 Mya following Ad-α event occurred ~ 18.7 Mya. Detailed annotations of NBS-LRRs or CBFs highlight the importance of genetic variations coming about after polyploidization in underpinning ability of immune responses or environmental adaptability. WGCNA analysis of postharvest quality indicators in combination with transcriptome revealed several transcription factors were involved in regulating ripening kiwi berry texture. Taking together, the assembly of an A. arguta tetraploid genome provides valuable resources in deciphering complex genome structure and facilitating functional genomics studies and genetic improvement for kiwifruit and other crops.

Moderating the neutralist-selectionist debate: exactly which propositions are we debating, and which arguments are valid?

Half a century after its foundation, the neutral theory of molecular evolution continues to attract controversy. The debate has been hampered by the coexistence of different interpretations of the core proposition of the neutral theory, the 'neutral mutation-random drift' hypothesis. In this review, we trace the origins of these ambiguities and suggest potential solutions. We highlight the difference between the original, the revised and the nearly neutral hypothesis, and re-emphasise that none of them equates to the null hypothesis of strict neutrality. We distinguish the neutral hypothesis of protein evolution, the main focus of the ongoing debate, from the neutral hypotheses of genomic and functional DNA evolution, which for many species are generally accepted. We advocate a further distinction between a narrow and an extended neutral hypothesis (of which the latter posits that random non-conservative amino acid substitutions can cause non-ecological phenotypic divergence), and we discuss the implications for evolutionary biology beyond the domain of molecular evolution. We furthermore point out that the debate has widened from its initial focus on point mutations, and also concerns the fitness effects of large-scale mutations, which can alter the dosage of genes and regulatory sequences. We evaluate the validity of neutralist and selectionist arguments and find that the tested predictions, apart from being sensitive to violation of underlying assumptions, are often derived from the null hypothesis of strict neutrality, or equally consistent with the opposing selectionist hypothesis, except when assuming molecular panselectionism. Our review aims to facilitate a constructive neutralist-selectionist debate, and thereby to contribute to answering a key question of evolutionary biology: what proportions of amino acid and nucleotide substitutions and polymorphisms are adaptive?