Adaptive evolution of the enigmatic Takakia now facing climate change in Tibet.

The most extreme environments are the most vulnerable to transformation under a rapidly changing climate. These ecosystems harbor some of the most specialized species, which will likely suffer the highest extinction rates. We document the steepest temperature increase (2010-2021) on record at altitudes of above 4,000 m, triggering a decline of the relictual and highly adapted moss Takakia lepidozioides. Its de-novo-sequenced genome with 27,467 protein-coding genes includes distinct adaptations to abiotic stresses and comprises the largest number of fast-evolving genes under positive selection. The uplift of the study site in the last 65 million years has resulted in life-threatening UV-B radiation and drastically reduced temperatures, and we detected several of the molecular adaptations of Takakia to these environmental changes. Surprisingly, specific morphological features likely occurred earlier than 165 mya in much warmer environments. Following nearly 400 million years of evolution and resilience, this species is now facing extinction.

Modeling the mosaic structure of bacterial genomes to infer their evolutionary history.

The chronology and phylogeny of bacterial evolution are difficult to reconstruct due to a scarce fossil record. The analysis of bacterial genomes remains challenging because of large sequence divergence, the plasticity of bacterial genomes due to frequent gene loss, horizontal gene transfer, and differences in selective pressure from one locus to another. Therefore, taking advantage of the rich and rapidly accumulating genomic data requires accurate modeling of genome evolution. An important technical consideration is that loci with high effective mutation rates may diverge beyond the detection limit of the alignment algorithms used, biasing the genome-wide divergence estimates toward smaller divergences. In this article, we propose a novel method to gain insight into bacterial evolution based on statistical properties of genome comparisons. We find that the length distribution of sequence matches is shaped by the effective mutation rates of different loci, by the horizontal transfers, and by the aligner sensitivity. Based on these inputs, we build a model and show that it accounts for the empirically observed distributions, taking the Enterobacteriaceae family as an example. Our method allows to distinguish segments of vertical and horizontal origins and to estimate the time divergence and exchange rate between any pair of taxa from genome-wide alignments. Based on the estimated time divergences, we construct a time-calibrated phylogenetic tree to demonstrate the accuracy of the method.

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

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

Dating Ammonia-Oxidizing Bacteria with Abundant Eukaryotic Fossils.

Evolution of a complete nitrogen (N) cycle relies on the onset of ammonia oxidation, which aerobically converts ammonia to nitrogen oxides. However, accurate estimation of the antiquity of ammonia-oxidizing bacteria (AOB) remains challenging because AOB-specific fossils are absent and bacterial fossils amenable to calibrate molecular clocks are rare. Leveraging the ancient endosymbiosis of mitochondria and plastid, as well as using state-of-the-art Bayesian sequential dating approach, we obtained a timeline of AOB evolution calibrated largely by eukaryotic fossils. We show that the first AOB evolved in marine Gammaproteobacteria (Gamma-AOB) and emerged between 2.1 and 1.9 billion years ago (Ga), thus postdating the Great Oxidation Event (GOE; 2.4 to 2.32 Ga). To reconcile the sedimentary N isotopic signatures of ammonia oxidation occurring near the GOE, we propose that ammonia oxidation likely occurred at the common ancestor of Gamma-AOB and Gammaproteobacterial methanotrophs, or the actinobacterial/verrucomicrobial methanotrophs which are known to have ammonia oxidation activities. It is also likely that nitrite was transported from the terrestrial habitats where ammonia oxidation by archaea took place. Further, we show that the Gamma-AOB predated the anaerobic ammonia-oxidizing (anammox) bacteria, implying that the emergence of anammox was constrained by the availability of dedicated ammonia oxidizers which produce nitrite to fuel anammox. Our work supports a new hypothesis that N redox cycle involving nitrogen oxides evolved rather late in the ocean.

Functional Optimization in Distinct Tissues and Conditions Constrains the Rate of Protein Evolution.

Understanding the main determinants of protein evolution is a fundamental challenge in biology. Despite many decades of active research, the molecular and cellular mechanisms underlying the substantial variability of evolutionary rates across cellular proteins are not currently well understood. It also remains unclear how protein molecular function is optimized in the context of multicellular species and why many proteins, such as enzymes, are only moderately efficient on average. Our analysis of genomics and functional datasets reveals in multiple organisms a strong inverse relationship between the optimality of protein molecular function and the rate of protein evolution. Furthermore, we find that highly expressed proteins tend to be substantially more functionally optimized. These results suggest that cellular expression costs lead to more pronounced functional optimization of abundant proteins and that the purifying selection to maintain high levels of functional optimality significantly slows protein evolution. We observe that in multicellular species both the rate of protein evolution and the degree of protein functional efficiency are primarily affected by expression in several distinct cell types and tissues, specifically, in developed neurons with upregulated synaptic processes in animals and in young and fast-growing tissues in plants. Overall, our analysis reveals how various constraints from the molecular, cellular, and species' levels of biological organization jointly affect the rate of protein evolution and the level of protein functional adaptation.

Loss of endogenous estrogen alters mitochondrial metabolism and muscle clock-related protein Rbm20 in female mdx mice.

Female carriers of a Duchenne muscular dystrophy (DMD) gene mutation manifest exercise intolerance and metabolic anomalies that may be exacerbated following menopause due to the loss of estrogen, a known regulator of skeletal muscle function and metabolism. Here, we studied the impact of estrogen depletion (via ovariectomy) on exercise tolerance and muscle mitochondrial metabolism in female mdx mice and the potential of estrogen replacement therapy (using estradiol) to protect against functional and metabolic perturbations. We also investigated the effect of estrogen depletion, and replacement, on the skeletal muscle proteome through an untargeted proteomic approach with TMT-labelling. Our study confirms that loss of estrogen in female mdx mice reduces exercise capacity, tricarboxylic acid cycle intermediates, and citrate synthase activity but that these deficits are offset through estrogen replacement therapy. Furthermore, ovariectomy downregulated protein expression of RNA-binding motif factor 20 (Rbm20), a critical regulator of sarcomeric and muscle homeostasis gene splicing, which impacted pathways involving ribosomal and mitochondrial translation. Estrogen replacement modulated Rbm20 protein expression and promoted metabolic processes and the upregulation of proteins involved in mitochondrial dynamics and metabolism. Our data suggest that estrogen mitigates dystrophinopathic features in female mdx mice and that estrogen replacement may be a potential therapy for post-menopausal DMD carriers.

Clockor2: Inferring Global and Local Strict Molecular Clocks Using Root-to-Tip Regression.

Molecular sequence data from rapidly evolving organisms are often sampled at different points in time. Sampling times can then be used for molecular clock calibration. The root-to-tip (RTT) regression is an essential tool to assess the degree to which the data behave in a clock-like fashion. Here, we introduce Clockor2, a client-side web application for conducting RTT regression. Clockor2 allows users to quickly fit local and global molecular clocks, thus handling the increasing complexity of genomic datasets that sample beyond the assumption of homogeneous host populations. Clockor2 is efficient, handling trees of up to the order of 104 tips, with significant speed increases compared with other RTT regression applications. Although clockor2 is written as a web application, all data processing happens on the client-side, meaning that data never leave the user's computer. Clockor2 is freely available at https://clockor2.github.io/.

LSM12 and ME31B/DDX6 Define Distinct Modes of Posttranscriptional Regulation by ATAXIN-2 Protein Complex in Drosophila Circadian Pacemaker Neurons.

ATAXIN-2 (ATX2) has been implicated in human neurodegenerative diseases, yet it remains elusive how ATX2 assembles specific protein complexes to execute its physiological roles. Here we employ the posttranscriptional co-activator function of Drosophila ATX2 to demonstrate that LSM12 and ME31B/DDX6 are two ATX2-associating factors crucial for sustaining circadian rhythms. LSM12 acts as a molecular adaptor for the recruitment of TWENTY-FOUR (TYF) to ATX2. The ATX2-LSM12-TYF complex thereby stimulates TYF-dependent translation of the rate-limiting clock gene period (per) to maintain 24 hr periodicity in circadian behaviors. In contrast, ATX2 contributes to NOT1-mediated gene silencing and associates with NOT1 in a ME31B/DDX6-dependent manner. The ME31B/DDX6-NOT1 complex does not affect PER translation but supports high-amplitude behavioral rhythms along with ATX2, indicating a PER-independent clock function of ATX2. Taken together, these data suggest that the ATX2 complex may switch distinct modes of posttranscriptional regulation through its associating factors to control circadian clocks and ATX2-related physiology.

Two light sensors decode moonlight versus sunlight to adjust a plastic circadian/circalunidian clock to moon phase.

Many species synchronize their physiology and behavior to specific hours. It is commonly assumed that sunlight acts as the main entrainment signal for ∼24-h clocks. However, the moon provides similarly regular time information. Consistently, a growing number of studies have reported correlations between diel behavior and lunidian cycles. Yet, mechanistic insight into the possible influences of the moon on ∼24-h timers remains scarce. We have explored the marine bristleworm Platynereis dumerilii to investigate the role of moonlight in the timing of daily behavior. We uncover that moonlight, besides its role in monthly timing, also schedules the exact hour of nocturnal swarming onset to the nights' darkest times. Our work reveals that extended moonlight impacts on a plastic clock that exhibits <24 h (moonlit) or >24 h (no moon) periodicity. Abundance, light sensitivity, and genetic requirement indicate that the Platynereis light receptor molecule r-Opsin1 serves as a receptor that senses moonrise, whereas the cryptochrome protein L-Cry is required to discriminate the proper valence of nocturnal light as either moonlight or sunlight. Comparative experiments in Drosophila suggest that cryptochrome's principle requirement for light valence interpretation is conserved. Its exact biochemical properties differ, however, between species with dissimilar timing ecology. Our work advances the molecular understanding of lunar impact on fundamental rhythmic processes, including those of marine mass spawners endangered by anthropogenic change.

Fossil-calibrated molecular clock data enable reconstruction of steps leading to differentiated multicellularity and anisogamy in the Volvocine algae.

BACKGROUND: Throughout its nearly four-billion-year history, life has undergone evolutionary transitions in which simpler subunits have become integrated to form a more complex whole. Many of these transitions opened the door to innovations that resulted in increased biodiversity and/or organismal efficiency. The evolution of multicellularity from unicellular forms represents one such transition, one that paved the way for cellular differentiation, including differentiation of male and female gametes. A useful model for studying the evolution of multicellularity and cellular differentiation is the volvocine algae, a clade of freshwater green algae whose members range from unicellular to colonial, from undifferentiated to completely differentiated, and whose gamete types can be isogamous, anisogamous, or oogamous. To better understand how multicellularity, differentiation, and gametes evolved in this group, we used comparative genomics and fossil data to establish a geologically calibrated roadmap of when these innovations occurred. RESULTS: Our ancestral-state reconstructions, show that multicellularity arose independently twice in the volvocine algae. Our chronograms indicate multicellularity evolved during the Carboniferous-Triassic periods in Goniaceae + Volvocaceae, and possibly as early as the Cretaceous in Tetrabaenaceae. Using divergence time estimates we inferred when, and in what order, specific developmental changes occurred that led to differentiated multicellularity and oogamy. We find that in the volvocine algae the temporal sequence of developmental changes leading to differentiated multicellularity is much as proposed by David Kirk, and that multicellularity is correlated with the acquisition of anisogamy and oogamy. Lastly, morphological, molecular, and divergence time data suggest the possibility of cryptic species in Tetrabaenaceae. CONCLUSIONS: Large molecular datasets and robust phylogenetic methods are bringing the evolutionary history of the volvocine algae more sharply into focus. Mounting evidence suggests that extant species in this group are the result of two independent origins of multicellularity and multiple independent origins of cell differentiation. Also, the origin of the Tetrabaenaceae-Goniaceae-Volvocaceae clade may be much older than previously thought. Finally, the possibility of cryptic species in the Tetrabaenaceae provides an exciting opportunity to study the recent divergence of lineages adapted to live in very different thermal environments.

Timing without coding: How do long non-coding RNAs regulate circadian rhythms?

Long non-coding RNAs (lncRNAs) are a new class of regulatory RNAs that play important roles in disease development and a variety of biological processes. Recent studies have underscored the importance of lncRNAs in the circadian clock system and demonstrated that lncRNAs regulate core clock genes and the core clock machinery in mammals. In this review, we provide an overview of our current understanding of how lncRNAs regulate the circadian clock without coding a protein. We also offer additional insights into the challenges in understanding the functions of lncRNAs and other unresolved questions in the field. We do not cover other regulatory ncRNAs even though they also play important roles; readers are highly encouraged to refer to other excellent reviews on this topic.

More than the clock: distinct regulation of muscle function and metabolism by PER2 and RORα.

Circadian rhythms, governed by the dominant central clock, in addition to various peripheral clocks, regulate almost all biological processes, including sleep-wake cycles, hormone secretion and metabolism. In certain contexts, the regulation and function of the peripheral oscillations can be decoupled from the central clock. However, the specific mechanisms underlying muscle-intrinsic clock-dependent modulation of muscle function and metabolism remain unclear. We investigated the outcome of perturbations of the primary and secondary feedback loops of the molecular clock in skeletal muscle by specific gene ablation of Period circadian regulator 2 (Per2) and RAR-related orphan receptor alpha (Rorα), respectively. In both models, a dampening of core clock gene oscillation was observed, while the phase was preserved. Moreover, both loops seem to be involved in the homeostasis of amine groups. Highly divergent outcomes were seen for overall muscle gene expression, primarily affecting circadian rhythmicity in the PER2 knockouts and non-oscillating genes in the RORα knockouts, leading to distinct outcomes in terms of metabolome and phenotype. These results highlight the entanglement of the molecular clock and muscle plasticity and allude to specific functions of different clock components, i.e. the primary and secondary feedback loops, in this context. The reciprocal interaction between muscle contractility and circadian clocks might therefore be instrumental to determining a finely tuned adaptation of muscle tissue to perturbations in health and disease. KEY POINTS: Specific perturbations of the primary and secondary feedback loop of the molecular clock result in specific outcomes on muscle metabolism and function. Ablation of Per2 (primary loop) or Rorα (secondary loop) blunts the amplitude of core clock genes, in absence of a shift in phase. Perturbation of the primary feedback loop by deletion of PER2 primarily affects muscle gene oscillation. Knockout of RORα and the ensuing modulation of the secondary loop results in the aberrant expression of a large number of non-clock genes and proteins. The deletion of PER2 and RORα affects muscle metabolism and contractile function in a circadian manner, highlighting the central role of the molecular clock in modulating muscle plasticity.

Detecting Episodic Evolution through Bayesian Inference of Molecular Clock Models.

Molecular evolutionary rate variation is a key aspect of the evolution of many organisms that can be modeled using molecular clock models. For example, fixed local clocks revealed the role of episodic evolution in the emergence of SARS-CoV-2 variants of concern. Like all statistical models, however, the reliability of such inferences is contingent on an assessment of statistical evidence. We present a novel Bayesian phylogenetic approach for detecting episodic evolution. It consists of computing Bayes factors, as the ratio of posterior and prior odds of evolutionary rate increases, effectively quantifying support for the effect size. We conducted an extensive simulation study to illustrate the power of this method and benchmarked it to formal model comparison of a range of molecular clock models using (log) marginal likelihood estimation, and to inference under a random local clock model. Quantifying support for the effect size has higher sensitivity than formal model testing and is straight-forward to compute, because it only needs samples from the posterior and prior distribution. However, formal model testing has the advantage of accommodating a wide range molecular clock models. We also assessed the ability of an automated approach, known as the random local clock, where branches under episodic evolution may be detected without their a priori definition. In an empirical analysis of a data set of SARS-CoV-2 genomes, we find "very strong" evidence for episodic evolution. Our results provide guidelines and practical methods for Bayesian detection of episodic evolution, as well as avenues for further research into this phenomenon.

Epimutations Define a Fast-Ticking Molecular Clock in Plants.

Stochastic gains and losses of DNA methylation at CG dinucleotides are a frequent occurrence in plants. These spontaneous 'epimutations' occur at a rate that is 100 000 times higher than the genetic mutation rate, are effectively neutral at the genome-wide scale, and are stably inherited across mitotic and meiotic cell divisions. Mathematical models have been extraordinarily successful at describing how epimutations accumulate in plant genomes over time, making this process one of the most predictable epigenetic phenomena to date. Here, we propose that their high rate and effective neutrality make epimutations a powerful new molecular clock for timing evolutionary events of the recent past and for age dating of long-lived perennials such as trees.

Diversity and biogeography of South American mud turtles elucidated by multilocus DNA sequencing (Testudines: Kinosternidae).

Kinosternon is the most speciose genus of extant turtles, with 22 currently recognized species, distributed across large parts of the Americas. Most species have small distributions, but K. leucostomum and K. scorpioides range from Mexico to South America. Previous studies have found discordance between mitochondrial and nuclear phylogenies in some kinosternid groups, with the current taxonomy following the nuclear-based results. Herein, based on extended molecular, geographic, and taxonomic sampling, we explore the phylogeographic structure and taxonomic limits for K. leucostomum and the K. scorpioides group and present a fossil-calibrated nuclear time tree for Kinosternon. Our results reveal contrasting differentiation patterns for the K. scorpioides group and K. leucostomum, despite overlapping distributions. Kinosternon leucostomum shows only shallow geographic divergence, whereas the K. scorpioides group is polyphyletic with up to 10 distinct taxa, some of them undescribed. We support the elevation of K. s. albogulare and K. s. cruentatum to species level. Given the deep divergence within the genus Kinosternon, we propose the recognition of three subgenera, Kinosternon, Cryptochelys and Thyrosternum, and the abandonment of the group-based classification, at least for the K. leucostomum and K. scorpioides groups. Our results show an initial split in Kinosternon that gave rise to two main radiations, one Nearctic and one mainly Neotropical. Most speciation events in Kinosternon occurred during the Quaternary and we hypothesize that they were mediated by both climatic and geological events. Additionally, our data imply that at least three South American colonizations occurred, two in the K. leucostomum group, and one in the K. scorpioides group. Additionally, we hypothesize that discordance between mitochondrial and nuclear phylogenetic signal is due to mitochondrial capture from an extinct kinosternine lineage.

Phylogeny and evolution of dissimilatory sulfite reduction in prokaryotes.

Sulfate is the second most common nonmetallic ion in modern oceans, as its concentration dramatically increased alongside tectonic activity and atmospheric oxidation in the Proterozoic. Microbial sulfate/sulfite metabolism, involving organic carbon or hydrogen oxidation, is linked to sulfur and carbon biogeochemical cycles. However, the coevolution of microbial sulfate/sulfite metabolism and Earth's history remains unclear. Here, we conducted a comprehensive phylogenetic analysis to explore the evolutionary history of the dissimilatory sulfite reduction (Dsr) pathway. The phylogenies of the Dsr-related genes presented similar branching patterns but also some incongruencies, indicating the complex origin and evolution of Dsr. Among these genes, dsrAB is the hallmark of sulfur-metabolizing prokaryotes. Our detailed analyses suggested that the evolution of dsrAB was shaped by vertical inheritance and multiple horizontal gene transfer events and that selection pressure varied across distinct lineages. Dated phylogenetic trees indicated that key evolutionary events of dissimilatory sulfur-metabolizing prokaryotes were related to the Great Oxygenation Event (2.4-2.0 Ga) and several geological events in the "Boring Billion" (1.8-0.8 Ga), including the fragmentation of the Columbia supercontinent (approximately 1.6 Ga), the rapid increase in marine sulfate (1.3-1.2 Ga), and the Neoproterozoic glaciation event (approximately 1.0 Ga). We also proposed that the voluminous iron formations (approximately 1.88 Ga) might have induced the metabolic innovation of iron reduction. In summary, our study provides new insights into Dsr evolution and a systematic view of the coevolution of dissimilatory sulfur-metabolizing prokaryotes and the Earth's environment.

Phylogeny and biogeography of harmochirine jumping spiders (Araneae: Salticidae).

We use ultraconserved elements (UCE) and Sanger data to study the phylogeny, age, and biogeographical history of harmochirine jumping spiders, a group that includes the species-rich genus Habronattus, whose remarkable courtship has made it the focus of studies of behaviour, sexual selection, and diversification. We recovered 1947 UCE loci from 43 harmochirine taxa and 4 outgroups, yielding a core dataset of 193 UCEs with at least 50 % occupancy. Concatenated likelihood and ASTRAL analyses confirmed the separation of harmochirines into two major clades, here designated the infratribes Harmochirita and Pellenita. Most are African or Eurasian with the notable exception of a clade of pellenites containing Habronattus and Pellenattus of the Americas and Havaika and Hivanua of the Pacific Islands. Biogeographical analysis using the DEC model favours a dispersal of the clade's ancestor from Eurasia to the Americas, from which Havaika's ancestor dispersed to Hawaii and Hivanua's ancestor to the Marquesas Islands. Divergence time analysis on 32 loci with 85 % occupancy, calibrated by fossils and island age, dates the dispersal to the Americas at approximately 4 to 6 million years ago. The explosive radiation of Habronattus perhaps began only about 4 mya. The phylogeny clarifies both the evolution of sexual traits (e.g., the terminal apophyses was enlarged in Pellenes and not subsequently lost) and the taxonomy. Habronattus is confirmed as monophyletic. Pellenattus is raised to the status of genus, and 13 species moved into it as new combinations. Bianor stepposus Logunov, 1991 is transferred to Sibianor, and Pellenes bulawayoensis Wesolowska, 1999 is transferred to Neaetha. A molecular clock rate estimate for spider UCEs is presented and its utility to inform prior distributions is discussed.

Molecular phylogenetic and estimation of evolutionary divergence and biogeography of the family Schizoparmaceae and allied families (Diaporthales, Ascomycota).

The Diaporthales includes 32 families, many of which are important plant pathogens, endophytes and saprobes, e.g., members of the families Pseudoplagiostomataceae, Pyrisporaceae and Schizoparmaceae. Nucleotide sequences derived from five genetic loci including: ITS, LSU, TEF1-α, TUB2 and RPB2 were used for Bayesian evolutionary analysis to determine divergence times and evolutionary relationships within the Schizoparmaceae. Molecular clock analyses revealed that the ancestor of Schizoparmaceae split during the Upper Cretaceous period approximately 75.7 Mya (95 % highest posterior density of 60.3-91.3 Mya). Reconstructing ancestral state in phylogenies (RASP) with using the Bayesian Binary Markov chain Monte Carlo (BBM) Method to reconstruct the historical biogeography for the family Schizoparmaceae indicated its most likely origin in Africa. Based on taxonomic and phylogenetic analyses, the Pseudoplagiostomataceae and Pyrisporaceae relationships were clarified and a total of four species described herein. For Pseudoplagiostomataceae, three new species and one known species that include, Pseudoplagiostoma fafuense sp. nov., Ps. ilicis sp. nov., Ps. sanmingense sp. nov. and Ps. bambusae are described and a key of Pseudoplagiostomataceae is provided. With respect to Pyrisporaceae, we considered Pseudoplagiostoma castaneae to be a synonym of Pyrispora castaneae. In addition, a new species of Schizoparmaceae, Coniella fujianensis sp. nov. is described and illustrated.

Evaluating the Accuracy of Methods for Detecting Correlated Rates of Molecular and Morphological Evolution.

Determining the link between genomic and phenotypic change is a fundamental goal in evolutionary biology. Insights into this link can be gained by using a phylogenetic approach to test for correlations between rates of molecular and morphological evolution. However, there has been persistent uncertainty about the relationship between these rates, partly because conflicting results have been obtained using various methods that have not been examined in detail. We carried out a simulation study to evaluate the performance of 5 statistical methods for detecting correlated rates of evolution. Our simulations explored the evolution of molecular sequences and morphological characters under a range of conditions. Of the methods tested, Bayesian relaxed-clock estimation of branch rates was able to detect correlated rates of evolution correctly in the largest number of cases. This was followed by correlations of root-to-tip distances, Bayesian model selection, independent sister-pairs contrasts, and likelihood-based model selection. As expected, the power to detect correlated rates increased with the amount of data, both in terms of tree size and number of morphological characters. Likewise, greater among-lineage rate variation in the data led to improved performance of all 5 methods, particularly for Bayesian relaxed-clock analysis when the rate model was mismatched. We then applied these methods to a data set from flowering plants and did not find evidence of a correlation in evolutionary rates between genomic data and morphological characters. The results of our study have practical implications for phylogenetic analyses of combined molecular and morphological data sets, and highlight the conditions under which the links between genomic and phenotypic rates of evolution can be evaluated quantitatively.

The Cretaceous Diatom Database: A tool for investigating early diatom evolution.

The Cretaceous period is the time of the first appearance of the diatoms in the fossil record. These fossils give us direct evidence of the age and early evolution of the diatom lineage. The fossil record, however, is incomplete and therefore often extrapolated through time-calibrated phylogenies. These two approaches offer different perspectives on the early evolution of diatoms, which is still poorly understood. We compiled the first comprehensive Cretaceous Diatom Database, a tool to investigate the taxonomy, diversity, and occurrence of the earliest known diatom lineages. To further aid the integration and use of the oldest diatom fossils in molecular clock analyses, we present a set of well-documented Cretaceous fossils that can be placed onto molecular phylogenetic trees of extant and extinct species, making them ideal candidates for the calibration of molecular clocks. The analysis of the fossil record and the Cretaceous Diatom Database revealed Cretaceous diversity is substantially greater than previously thought, yet considerable taxonomic work is still needed. The Cretaceous Diatom Database and the list of Cretaceous fossils for calibrating molecular clocks represent valuable resources for future evolutionary and taxonomic studies of modern and fossil diatoms.