Phylogenomic discovery of deleterious mutations facilitates hybrid potato breeding.

Hybrid potato breeding will transform the crop from a clonally propagated tetraploid to a seed-reproducing diploid. Historical accumulation of deleterious mutations in potato genomes has hindered the development of elite inbred lines and hybrids. Utilizing a whole-genome phylogeny of 92 Solanaceae and its sister clade species, we employ an evolutionary strategy to identify deleterious mutations. The deep phylogeny reveals the genome-wide landscape of highly constrained sites, comprising ∼2.4% of the genome. Based on a diploid potato diversity panel, we infer 367,499 deleterious variants, of which 50% occur at non-coding and 15% at synonymous sites. Counterintuitively, diploid lines with relatively high homozygous deleterious burden can be better starting material for inbred-line development, despite showing less vigorous growth. Inclusion of inferred deleterious mutations increases genomic-prediction accuracy for yield by 24.7%. Our study generates insights into the genome-wide incidence and properties of deleterious mutations and their far-reaching consequences for breeding.

Expansion of the global RNA virome reveals diverse clades of bacteriophages.

High-throughput RNA sequencing offers broad opportunities to explore the Earth RNA virome. Mining 5,150 diverse metatranscriptomes uncovered >2.5 million RNA virus contigs. Analysis of >330,000 RNA-dependent RNA polymerases (RdRPs) shows that this expansion corresponds to a 5-fold increase of the known RNA virus diversity. Gene content analysis revealed multiple protein domains previously not found in RNA viruses and implicated in virus-host interactions. Extended RdRP phylogeny supports the monophyly of the five established phyla and reveals two putative additional bacteriophage phyla and numerous putative additional classes and orders. The dramatically expanded phylum Lenarviricota, consisting of bacterial and related eukaryotic viruses, now accounts for a third of the RNA virome. Identification of CRISPR spacer matches and bacteriolytic proteins suggests that subsets of picobirnaviruses and partitiviruses, previously associated with eukaryotes, infect prokaryotic hosts.

Quantitative fate mapping: A general framework for analyzing progenitor state dynamics via retrospective lineage barcoding.

Natural and induced somatic mutations that accumulate in the genome during development record the phylogenetic relationships of cells; whether these lineage barcodes capture the complex dynamics of progenitor states remains unclear. We introduce quantitative fate mapping, an approach to reconstruct the hierarchy, commitment times, population sizes, and commitment biases of intermediate progenitor states during development based on a time-scaled phylogeny of their descendants. To reconstruct time-scaled phylogenies from lineage barcodes, we introduce Phylotime, a scalable maximum likelihood clustering approach based on a general barcoding mutagenesis model. We validate these approaches using realistic in silico and in vitro barcoding experiments. We further establish criteria for the number of cells that must be analyzed for robust quantitative fate mapping and a progenitor state coverage statistic to assess the robustness. This work demonstrates how lineage barcodes, natural or synthetic, enable analyzing progenitor fate and dynamics long after embryonic development in any organism.

Characterizing genetic intra-tumor heterogeneity across 2,658 human cancer genomes.

Intra-tumor heterogeneity (ITH) is a mechanism of therapeutic resistance and therefore an important clinical challenge. However, the extent, origin, and drivers of ITH across cancer types are poorly understood. To address this, we extensively characterize ITH across whole-genome sequences of 2,658 cancer samples spanning 38 cancer types. Nearly all informative samples (95.1%) contain evidence of distinct subclonal expansions with frequent branching relationships between subclones. We observe positive selection of subclonal driver mutations across most cancer types and identify cancer type-specific subclonal patterns of driver gene mutations, fusions, structural variants, and copy number alterations as well as dynamic changes in mutational processes between subclonal expansions. Our results underline the importance of ITH and its drivers in tumor evolution and provide a pan-cancer resource of comprehensively annotated subclonal events from whole-genome sequencing data.

Identification of novel bat coronaviruses sheds light on the evolutionary origins of SARS-CoV-2 and related viruses.

Despite the discovery of animal coronaviruses related to SARS-CoV-2, the evolutionary origins of this virus are elusive. We describe a meta-transcriptomic study of 411 bat samples collected from a small geographical region in Yunnan province, China, between May 2019 and November 2020. We identified 24 full-length coronavirus genomes, including four novel SARS-CoV-2-related and three SARS-CoV-related viruses. Rhinolophus pusillus virus RpYN06 was the closest relative of SARS-CoV-2 in most of the genome, although it possessed a more divergent spike gene. The other three SARS-CoV-2-related coronaviruses carried a genetically distinct spike gene that could weakly bind to the hACE2 receptor in vitro. Ecological modeling predicted the co-existence of up to 23 Rhinolophus bat species, with the largest contiguous hotspots extending from South Laos and Vietnam to southern China. Our study highlights the remarkable diversity of bat coronaviruses at the local scale, including close relatives of both SARS-CoV-2 and SARS-CoV.

What one genus of showy moths can say about migration, adaptation, and wing pattern.

The Ornate Moth, Utetheisa ornatrix, has served as a model species in chemical ecology studies for decades. Like in the widely publicized stories of the Monarch and other milkweed butterflies, the Ornate Moth and its relatives are tropical insects colonizing whole continents assisted by their chemical defenses. With the recent advances in genomic techniques and evo-devo research, it is becoming a model for studies in other areas, from wing pattern development to phylogeography, from toxicology to epigenetics. We used a genomic approach to learn about Utetheisa's evolution, detoxification, dispersal abilities, and wing pattern diversity. We present an evolutionary genomic analysis of the worldwide genus Utetheisa, then focusing on U. ornatrix. Our reference genome of U. ornatrix reveals gene duplications in the regions possibly associated with detoxification abilities, which allows them to feed on toxic food plants. Finally, comparative genomic analysis of over 100 U. ornatrix specimens from the museum with apparent differences in wing patterns suggest the potential roles of cortex and lim3 genes in wing pattern formation of Lepidoptera and the utility of museum-preserved collection specimens for wing pattern research.

Revisiting the four Hexapoda classes: Protura as the sister group to all other hexapods.

Insects represent the most diverse animal group, yet previous phylogenetic analyses based on morphological and molecular data have failed to agree on the evolutionary relationships of early insects and their six-legged relatives (together constituting the clade Hexapoda). In particular, the phylogenetic positions of the three early-diverging hexapod lineages-the coneheads (Protura), springtails (Collembola), and two-pronged bristletails (Diplura)-have been debated for over a century, with alternative topologies implying drastically different scenarios of the evolution of the insect body plan and hexapod terrestrialization. We addressed this issue by sampling all hexapod orders and experimenting with a broad range of across-site compositional heterogeneous models designed to tackle ancient divergences. Our analyses support Protura as the earliest-diverging hexapod lineage ("Protura-sister") and Collembola as a sister group to Diplura, a clade corresponding to the original composition of Entognatha, and characterized by the shared possession of internal muscles in the antennal flagellum. The previously recognized 'Elliplura' hypothesis is recovered only under the site-homogeneous substitution models with partial supermatrices. Our cross-validation analysis shows that the site-heterogeneous CAT-GTR model, which recovers "Protura-sister," fits significantly better than homogeneous models. Furthermore, the morphologically unusual Protura are also supported as the earliest-diverging hexapod lineage by other lines of evidence, such as mitogenomes, comparative embryology, and sperm morphology, which produced results similar to those in this study. Our backbone phylogeny of hexapods will facilitate the exploration of the underpinnings of hexapod terrestrialization and megadiversity.

A region of suppressed recombination misleads neoavian phylogenomics.

Genomes are typically mosaics of regions with different evolutionary histories. When speciation events are closely spaced in time, recombination makes the regions sharing the same history small, and the evolutionary history changes rapidly as we move along the genome. When examining rapid radiations such as the early diversification of Neoaves 66 Mya, typically no consistent history is observed across segments exceeding kilobases of the genome. Here, we report an exception. We found that a 21-Mb region in avian genomes, mapped to chicken chromosome 4, shows an extremely strong and discordance-free signal for a history different from that of the inferred species tree. Such a strong discordance-free signal, indicative of suppressed recombination across many millions of base pairs, is not observed elsewhere in the genome for any deep avian relationships. Although long regions with suppressed recombination have been documented in recently diverged species, our results pertain to relationships dating circa 65 Mya. We provide evidence that this strong signal may be due to an ancient rearrangement that blocked recombination and remained polymorphic for several million years prior to fixation. We show that the presence of this region has misled previous phylogenomic efforts with lower taxon sampling, showing the interplay between taxon and locus sampling. We predict that similar ancient rearrangements may confound phylogenetic analyses in other clades, pointing to a need for new analytical models that incorporate the possibility of such events.

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.

Helix-loop-helix proteins and the advent of cellular diversity: 30 years of discovery.

Helix-loop-helix (HLH) proteins are dimeric transcription factors that control lineage- and developmental-specific gene programs. Genes encoding for HLH proteins arose in unicellular organisms >600 million years ago and then duplicated and diversified from ancestral genes across the metazoan and plant kingdoms to establish multicellularity. Hundreds of HLH proteins have been identified with diverse functions in a wide variety of cell types. HLH proteins orchestrate lineage specification, commitment, self-renewal, proliferation, differentiation, and homing. HLH proteins also regulate circadian clocks, protect against hypoxic stress, promote antigen receptor locus assembly, and program transdifferentiation. HLH proteins deposit or erase epigenetic marks, activate noncoding transcription, and sequester chromatin remodelers across the chromatin landscape to dictate enhancer-promoter communication and somatic recombination. Here the evolution of HLH genes, the structures of HLH domains, and the elaborate activities of HLH proteins in multicellular life are discussed.

Hox, homology, and parsimony: An organismal perspective.

Hox genes are important regulators in animal development. They often show a mosaic of conserved (e.g., longitudinal axis patterning) and lineage-specific novel functions (e.g., development of skeletal, sensory, or locomotory systems). Despite extensive research over the past decades, it remains controversial at which node in the animal tree of life the Hox cluster evolved. Its presence already in the last common metazoan ancestor has been proposed, although the genomes of both putative earliest extant metazoan offshoots, the ctenophores and the poriferans, are devoid of Hox sequences. The lack of Hox genes in the supposedly "simple"-built poriferans and their low number in cnidarians and the basally branching bilaterians, the xenacoelomorphs, seems to support the classical notion that the number of Hox genes is correlated with the degree of animal complexity. However, the 4-fold increase of the Hox cluster in xiphosurans, a basally branching chelicerate clade, as well as the situation in some teleost fishes that show a multitude of Hox genes compared to, e.g., human, demonstrates, that there is no per se direct correlation between organismal complexity and Hox number. Traditional approaches have tried to base homology on the morphological level on shared expression profiles of individual genes, but recent data have shown that, in particular with respect to Hox and other regulatory genes, complex gene-gene interactions rather than expression signatures of individual genes alone are responsible for shaping morphological traits during ontogeny. Accordingly, for sound homology assessments and reconstructions of character evolution on organ system level, additional independent datasets (e.g., morphological, developmental) need to be included in any such analyses. If supported by solid data, proposed structural homology should be regarded as valid and not be rejected solely on the grounds of non-parsimonious distribution of the character over a given phylogenetic topology.

Molecular evolution of the Thrombospondin superfamily.

Thrombospondins (TSPs) are multidomain, calcium-binding glycoproteins that have wide-ranging roles in vertebrates in cell interactions, extracellular matrix (ECM) organisation, angiogenesis, tissue remodelling, synaptogenesis, and also in musculoskeletal and cardiovascular functions. Land animals encode five TSPs, which assembly co-translationally either as trimers (subgroup A) or pentamers (subgroup B). The vast majority of research has focused on this canonical TSP family, which evolved through the whole-genome duplications that took place early in the vertebrate lineage. With benefit of the growth in genome- and transcriptome-predicted proteomes of a much wider range of animal species, examination of TSPs throughout metazoan phyla has revealed extensive conservation of subgroup B-type TSPs in invertebrates. In addition, these searches established that canonical TSPs are, in fact, one branch within a TSP superfamily that includes other clades designated mega-TSPs, sushi-TSPs and poriferan-TSPs. Despite the apparent simplicity of poriferans and cnidarians as organisms, these phyla encode a greater diversity of TSP superfamily members than vertebrates. We discuss here the molecular characteristics of the TSP superfamily members, current knowledge of their expression profiles and functions in invertebrates, and models for the evolution of this complex ECM superfamily.

The rise of taxon-specific epitope predictors.

Computational predictors of immunogenic peptides, or epitopes, are traditionally built based on data from a broad range of pathogens without consideration for taxonomic information. While this approach may be reasonable if one aims to develop one-size-fits-all models, it may be counterproductive if the proteins for which the model is expected to generalize are known to come from a specific subset of phylogenetically related pathogens. There is mounting evidence that, for these cases, taxon-specific models can outperform generalist ones, even when trained with substantially smaller amounts of data. In this comment, we provide some perspective on the current state of taxon-specific modelling for the prediction of linear B-cell epitopes, and the challenges faced when building and deploying these predictors.

Expanding Archaeal Diversity and Phylogeny: Past, Present, and Future.

The discovery of the Archaea is a major scientific hallmark of the twentieth century. Since then, important features of their cell biology, physiology, ecology, and diversity have been revealed. Over the course of some 40 years, the diversity of known archaea has expanded from 2 to about 30 phyla comprising over 20,000 species. Most of this archaeal diversity has been revealed by environmental 16S rRNA gene amplicon sequencing surveys using a broad range of universal and targeted primers. Of the few primers that target a large fraction of known archaeal diversity, all display a bias against recently discovered lineages, which limits studies aiming to survey overall archaeal diversity. Induced by genomic exploration of archaeal diversity, and improved phylogenomics approaches, archaeal taxonomic classification has been frequently revised. Due to computational limitations and continued discovery of new lineages, a stable archaeal phylogeny is not yet within reach. Obtaining phylogenetic and taxonomic consensus of archaea should be a high priority for the archaeal research community.

Complementing model species with model clades.

Model species continue to underpin groundbreaking plant science research. At the same time, the phylogenetic resolution of the land plant Tree of Life continues to improve. The intersection of these two research paths creates a unique opportunity to further extend the usefulness of model species across larger taxonomic groups. Here we promote the utility of the Arabidopsis thaliana model species, especially the ability to connect its genetic and functional resources, to species across the entire Brassicales order. We focus on the utility of using genomics and phylogenomics to bridge the evolution and diversification of several traits across the Brassicales to the resources in Arabidopsis, thereby extending scope from a model species by establishing a "model clade". These Brassicales-wide traits are discussed in the context of both the model species Arabidopsis thaliana and the family Brassicaceae. We promote the utility of such a "model clade" and make suggestions for building global networks to support future studies in the model order Brassicales.

An eosimiid primate of South Asian affinities in the Paleogene of Western Amazonia and the origin of New World monkeys.

Recent fossil discoveries in Western Amazonia revealed that two distinct anthropoid primate clades of African origin colonized South America near the Eocene/Oligocene transition (ca. 34 Ma). Here, we describe a diminutive fossil primate from Brazilian Amazonia and suggest that, surprisingly, a third clade of anthropoids was involved in the Paleogene colonization of South America by primates. This new taxon, Ashaninkacebus simpsoni gen. et sp. nov., has strong dental affinities with Asian African stem anthropoids: the Eosimiiformes. Morphology-based phylogenetic analyses of early Old World anthropoids and extinct and extant New World monkeys (platyrrhines) support relationships of both Ashaninkacebus and Amamria (late middle Eocene, North Africa) to the South Asian Eosimiidae. Afro-Arabia, then a mega island, played the role of a biogeographic stopover between South Asia and South America for anthropoid primates and hystricognathous rodents. The earliest primates from South America bear little adaptive resemblance to later Oligocene-early Miocene platyrrhine monkeys, and the scarcity of available paleontological data precludes elucidating firmly their affinities with or within Platyrrhini. Nonetheless, these data shed light on some of their life history traits, revealing a particularly small body size and a diet consisting primarily of insects and possibly fruit, which would have increased their chances of survival on a natural floating island during this extraordinary over-water trip to South America from Africa. Divergence-time estimates between Old and New World taxa indicate that the transatlantic dispersal(s) could source in the intense flooding events associated with the late middle Eocene climatic optimum (ca. 40.5 Ma) in Western Africa.

Correlated substitutions reveal SARS-like coronaviruses recombine frequently with a diverse set of structured gene pools.

Quantifying SARS-like coronavirus (SL-CoV) evolution is critical to understanding the origins of SARS-CoV-2 and the molecular processes that could underlie future epidemic viruses. While genomic analyses suggest recombination was a factor in the emergence of SARS-CoV-2, few studies have quantified recombination rates among SL-CoVs. Here, we infer recombination rates of SL-CoVs from correlated substitutions in sequencing data using a coalescent model with recombination. Our computationally-efficient, non-phylogenetic method infers recombination parameters of both sampled sequences and the unsampled gene pools with which they recombine. We apply this approach to infer recombination parameters for a range of positive-sense RNA viruses. We then analyze a set of 191 SL-CoV sequences (including SARS-CoV-2) and find that ORF1ab and S genes frequently undergo recombination. We identify which SL-CoV sequence clusters have recombined with shared gene pools, and show that these pools have distinct structures and high recombination rates, with multiple recombination events occurring per synonymous substitution. We find that individual genes have recombined with different viral reservoirs. By decoupling contributions from mutation and recombination, we recover the phylogeny of non-recombined portions for many of these SL-CoVs, including the position of SARS-CoV-2 in this clonal phylogeny. Lastly, by analyzing >400,000 SARS-CoV-2 whole genome sequences, we show current diversity levels are insufficient to infer the within-population recombination rate of the virus since the pandemic began. Our work offers new methods for inferring recombination rates in RNA viruses with implications for understanding recombination in SARS-CoV-2 evolution and the structure of clonal relationships and gene pools shaping its origins.

Scaphopoda is the sister taxon to Bivalvia: Evidence of ancient incomplete lineage sorting.

The almost simultaneous emergence of major animal phyla during the early Cambrian shaped modern animal biodiversity. Reconstructing evolutionary relationships among such closely spaced branches in the animal tree of life has proven to be a major challenge, hindering understanding of early animal evolution and the fossil record. This is particularly true in the species-rich and highly varied Mollusca where dramatic inconsistency among paleontological, morphological, and molecular evidence has led to a long-standing debate about the group's phylogeny and the nature of dozens of enigmatic fossil taxa. A critical step needed to overcome this issue is to supplement available genomic data, which is plentiful for well-studied lineages, with genomes from rare but key lineages, such as Scaphopoda. Here, by presenting chromosome-level genomes from both extant scaphopod orders and leveraging complete genomes spanning Mollusca, we provide strong support for Scaphopoda as the sister taxon of Bivalvia, revitalizing the morphology-based Diasoma hypothesis originally proposed 50 years ago. Our molecular clock analysis confidently dates the split between Bivalvia and Scaphopoda at ~520 Ma, prompting a reinterpretation of controversial laterally compressed Early Cambrian fossils, including Anabarella, Watsonella, and Mellopegma, as stem diasomes. Moreover, we show that incongruence in the phylogenetic placement of Scaphopoda in previous phylogenomic studies was due to ancient incomplete lineage sorting (ILS) that occurred during the rapid radiation of Conchifera. Our findings highlight the need to consider ILS as a potential source of error in deep phylogeny reconstruction, especially in the context of the unique nature of the Cambrian Explosion.

Testing the large genome constraint hypothesis in tropical rhizomatous herbs: life strategies, plant traits and habitat preferences in gingers.

Plant species with large genomes tend to be excluded from climatically more extreme environments with a shorter growing season. Species that occupy such environments are assumed to be under natural selection for more rapid growth and smaller genome size (GS). However, evidence for this is available only for temperate organisms. Here, we study the evolution of GS in two subfamilies of the tropical family Zingiberaceae to find out whether species with larger genomes are confined to environments where the vegetative season is longer. We tested our hypothesis on 337 ginger species from regions with contrasting climates by correlating their GS with an array of plant traits and environmental variables. We revealed 16-fold variation in GS which was tightly related to shoot seasonality. Negative correlations of GS with latitude, temperature and precipitation emerged in the subfamily Zingiberoidae, demonstrating that species with larger GS are excluded from areas with a shorter growing season. In the subfamily Alpinioideae, GS turned out to be correlated with the type of stem and light requirements and its members cope with seasonality mainly by adaptation to shady and moist habitats. The Ornstein-Uhlenbeck models suggested that evolution in regions with humid climates favoured larger GS than in drier regions. Our results indicate that climate seasonality exerts an upper constraint on GS not only in temperate regions but also in the tropics, unless species with large genomes find alternative ways to escape from that constraint.

Insertions and Deletions: Computational Methods, Evolutionary Dynamics, and Biological Applications.

Insertions and deletions constitute the second most important source of natural genomic variation. Insertions and deletions make up to 25% of genomic variants in humans and are involved in complex evolutionary processes including genomic rearrangements, adaptation, and speciation. Recent advances in long-read sequencing technologies allow detailed inference of insertions and deletion variation in species and populations. Yet, despite their importance, evolutionary studies have traditionally ignored or mishandled insertions and deletions due to a lack of comprehensive methodologies and statistical models of insertions and deletion dynamics. Here, we discuss methods for describing insertions and deletion variation and modeling insertions and deletions over evolutionary time. We provide practical advice for tackling insertions and deletions in genomic sequences and illustrate our discussion with examples of insertions and deletion-induced effects in human and other natural populations and their contribution to evolutionary processes. We outline promising directions for future developments in statistical methodologies that would allow researchers to analyze insertions and deletion variation and their effects in large genomic data sets and to incorporate insertions and deletions in evolutionary inference.