The Impact of Reconstruction Methods, Phylogenetic Uncertainty and Branch Lengths on Inference of Chromosome Number Evolution in American Daisies (Melampodium, Asteraceae)

Drift drives the evolution of chromosome number II: The impact of range size on genome evolution in Carnivora

Chromosome number is a fundamental genomic trait that is often the first recorded characteristic of a genome. Across large clades, a common pattern emerges: many or even most lineages exhibit relative stasis, while a handful of lineages or species exhibit striking variation. Despite recent developments in comparative methods, most of this heterogeneity is still poorly understood. It is essential to understand why some lineages have rapid rates of chromosome number evolution, as it can impact a variety of other traits. Previous research suggests that biased female meiotic drive may shape rates of karyotype evolution in some mammals. However, Carnivora exhibits variation that this female meiotic drive model cannot explain. We hypothesize that variation in effective population size may underlie rate variation in Carnivora. To test this hypothesis, we estimated rates of fusions and fissions while accounting for range size, which we use as a proxy for effective population size. We reason fusions and fissions are deleterious or underdominant and that only in lineages with small range sizes will these changes be able to fix due to genetic drift. In this study, we find that the rates of fusions and fissions are elevated in taxa with small range sizes relative to those with large range sizes. Based on these findings, we conclude that 1) naturally occurring structural mutations that change chromosome number are underdominant or mildly deleterious, and 2) when population sizes are small, structural rearrangements may play an important role in speciation and reduction in gene flow among populations.

A non-homogeneous model of chromosome-number evolution to reveal shifts in the transition patterns across the phylogeny

Changes in chromosome numbers, including polyploidy and dysploidy events, play a key role in eukaryote evolution as they could expediate reproductive isolation and have the potential to foster phenotypic diversification. Deciphering the pattern of chromosome-number change within a phylogeny currently relies on probabilistic evolutionary models. All currently available models assume time homogeneity, such that the transition rates are identical throughout the phylogeny. Here, we develop heterogeneous models of chromosome-number evolution that allow multiple transition regimes to operate in distinct parts of the phylogeny. The partition of the phylogeny to distinct transition regimes may be specified by the researcher or, alternatively, identified using a sequential testing approach. Once the number and locations of shifts in the transition pattern are determined, a second search phase identifies regimes with similar transition dynamics, which could indicate on convergent evolution. Using simulations, we study the performance of the developed model to detect shifts in patterns of chromosome-number evolution and demonstrate its applicability by analyzing the evolution of chromosome numbers within the Cyperaceae plant family. The developed model extends the capabilities of probabilistic models of chromosome-number evolution and should be particularly helpful for the analyses of large phylogenies that include multiple distinct subclades.

An Overview to the Index to Chromosome Numbers in Asteraceae Database: Revisiting Base Chromosome Numbers, Polyploidy, Descending Dysploidy, and Hybridization

A brief overview to the Index to Chromosome Numbers in Asteraceae database is provided. The database contains karyological information on Asteraceae and has been repeatedly improved and updated and is now hosted at the National Bioscience Database center. Also, we take the opportunity to revisit the evolution of base chromosome numbers in Asteraceae, emphasizing the phenomena of polyploidy, descending dysploidy, and hybridization, common in the family. Chromosome numbers for species included in one of the most recent phylogenetic treatments of the Asteraceae were obtained from the Index to Chromosome Numbers in Asteraceae database were mapped on to the modified phylogeny diagram, and base chromosome numbers were determined for each branch of the phylogeny. Results for tribal base numbers were the same as those hypothesized in our previous work with additional base numbers added for tribes not previously recognized but supported by newer phylogenetic methods. The Asteraceae show an ancestral base chromosome number of x = 9 and originated in the Antarctica (Gondowanaland) in Cretaceous (80 Mys ago). The x = 9 number has been retained through successive South American lineages of the Barnadesieeae, Gochnatieae, Stiffieae, Wunderlichieae, Astereae, and Senecioneae following northward migration. Northward migration to Africa was accompanied with x = 10 becoming the dominant base chromosome number as the family evolved multiple additional tribes. Northward migration to Australasia with x = 9 was in Astereae and the families Goodeneaseae, Menyanthaceae, and Stylydiaceae. The evolution of the North American Heliantheae alliance began with the appearance of x2 = 19 which persisted in multiple additional new tribes. Frequent dysploidy decreases, polyploidy and hybridization occurred throughout the history of the family.

Genome sequencing reveals chromosome fusion and extensive expansion of genes related to secondary metabolism in Artemisia argyi

Artemisia argyi, as famous as Artemisia annua, is a medicinal plant with huge economic value in the genus of Artemisia and has been widely used in the world for about 3000 years. However, a lack of the reference genome severely hinders the understanding of genetic basis for the active ingredient synthesis of A. argyi. Here, we firstly report a complex chromosome-level genome assembly of A. argyi with a large size of 8.03 Gb, with features of high heterozygosity (2.36%), high repetitive sequences (73.59%) and a huge number of protein-coding genes (279 294 in total). The assembly reveals at least three rounds of whole-genome duplication (WGD) events, including a recent WGD event in the A. argyi genome, and a recent burst of transposable element, which may contribute to its large genome size. The genomic data and karyotype analyses confirmed that A. argyi is an allotetraploid with 34 chromosomes. Intragenome synteny analysis revealed that chromosomes fusion event occurred in the A. argyi genome, which elucidates the changes in basic chromosome numbers in Artemisia genus. Significant expansion of genes related to photosynthesis, DNA replication, stress responses and secondary metabolism were identified in A. argyi, explaining the extensive environmental adaptability and rapid growth characteristics. In addition, we analysed genes involved in the biosynthesis pathways of flavonoids and terpenoids, and found that extensive gene amplification and tandem duplication contributed to the high contents of metabolites in A. argyi. Overall, the reference genome assembly provides scientific support for evolutionary biology, functional genomics and breeding in A. argyi and other Artemisia species.

The diverse terminology of reptile eggshell microstructure and its effect on phylogenetic comparative analyses

Reptile eggshell ensures water and gas exchange during incubation and plays a key role in reproductive success. The diversity of reptilian incubation and life history strategies has led to many clade-specific structural adaptations of their eggshell, which have been studied in extant taxa (i.e. birds, crocodilians, turtles, and lepidosaurs). Most studies on non-avian eggshells were performed over 30 years ago and categorized reptile eggshells into two main types: "hard" and "soft" - sometimes with a third intermediate category, "semi-rigid." In recent years, however, debate over the evolution of eggshell structure of major reptile clades has revealed how definitions of hard and soft eggshells influence inferred deep-time evolutionary patterns. Here, we review the diversity of extant and fossil eggshell with a focus on major reptile clades, and the criteria that have been used to define hard, soft, and semi-rigid eggshells. We show that all scoring approaches that retain these categories discretize continuous quantitative traits (e.g. eggshell thickness) and do not consider independent variation of other functionally important microstructural traits (e.g. degree of calcification, shell unit inner structure). We demonstrate the effect of three published approaches to discretizing eggshell type into hard, semi-rigid, and soft on ancestral state reconstructions using 200+ species representing all major extant and extinct reptile clades. These approaches result in different ancestral states for all major clades including Archosauria and Dinosauria, despite a difference in scoring for only 1-4% of the sample. Proposed scenarios of reptile eggshell evolution are highly conditioned by sampling, tree calibration, and lack of congruence between definitions of eggshell type. We conclude that the traditional "soft/hard/semi-rigid" classification of reptilian eggshells should be abandoned and provide guidelines for future descriptions focusing on specific functionally relevant characteristics (e.g. inner structures of shell units, pores, and membrane elements), analyses of these traits in a phylogenetic context, and sampling of previously undescribed taxa, including fossil eggs.

Chromosome restructuring and number change during the evolution of Morus notabilis and Morus alba

Mulberry (Morus spp.) is an economically important plant as the main food plant used for rearing domesticated silkworm and it has multiple uses in traditional Chinese medicine. Two basic chromosome numbers (Morus notabilis, n = 7, and Morus alba, n = 14) have been reported in the genus Morus, but the evolutionary history and relationship between them remain unclear. In the present study, a 335-Mb high-quality chromosome-scale genome was assembled for the wild mulberry species M. notabilis. Comparative genomic analyses indicated high chromosomal synteny between the 14 chromosomes of cultivated M. alba and the six chromosomes of wild M. notabilis. These results were successfully verified by fluorescence in situ hybridization. Chromosomal fission/fusion events played crucial roles in the chromosome restructuring process between M. notabilis and M. alba. The activity of the centromere was another key factor that ensured the stable inheritance of chromosomes. Our results also revealed that long terminal repeat retrotransposons were a major driver of the genome divergence and evolution of the mulberry genomes after they diverged from each other. This study provides important insights and a solid foundation for studying the evolution of mulberry, allowing the accelerated genetic improvement of cultivated mulberry species.

Do Meristic Characters Used in Phylogenetic Analysis Evolve in an Ordered Manner?

The use of ordered characters in phylogenetic analysis has been inconsistent throughout the history of phylogenetic inference. It has become more widespread in recent years, and some have advocated that all characters representing continuous or meristic traits should be ordered as a matter of course. Here, using the example of dental evolution, we examine two factors that may impact on whether meristic characters actually evolve in an ordered manner: the regulatory hierarchy governing the development of teeth that allows large sections of the entire tooth row to be suppressed in a single transition and regionalization of the tooth row where different modules have a degree of independence in their evolution. These are studied using both empirical and simulated data. Models of evolution of such characters are examined over molecular phylogenies to see if ordered or unordered models fit best. Simulations of tooth-row evolution are designed to incorporate changes in region size and multiple levels of developmental control to suppress individual regions or the entire row. The empirical analyses show that in a clade with largely homodont dentition the characters evolve in an ordered manner, but if dentition is heterodont with distinct regionalization their evolution better fits an unordered model. In the simulations, even if teeth are added and removed from the tooth row in an ordered manner, dividing the row into independently evolving modules can lead to characters covering multiple modules better fitting an unordered model of evolution. Adding the ability to suppress regions or the entire tooth row has a variable effect depending on the rates of suppression relative to the rates of addition and subtraction of individual teeth. We therefore advise not following a single policy when deciding whether to order meristic traits but to base the decision on a priori knowledge of the focal clade's evolution and developmental biology. [Discrete characters; ordered characters; phylogeny; teeth.].

Lineage-specific patterns of chromosome evolution are the rule not the exception in Polyneoptera insects

The structure of a genome can be described at its simplest by the number of chromosomes and the sex chromosome system it contains. Despite over a century of study, the evolution of genome structure on this scale remains recalcitrant to broad generalizations that can be applied across clades. To address this issue, we have assembled a dataset of 823 karyotypes from the insect group Polyneoptera. This group contains orders with a range of variations in chromosome number, and offer the opportunity to explore the possible causes of these differences. We have analysed these data using both phylogenetic and taxonomic approaches. Our analysis allows us to assess the importance of rates of evolution, phylogenetic history, sex chromosome systems, parthenogenesis and genome size on variation in chromosome number within clades. We find that fusions play a key role in the origin of new sex chromosomes, and that orders exhibit striking differences in rates of fusions, fissions and polyploidy. Our results suggest that the difficulty in finding consistent rules that govern evolution at this scale may be due to the presence of many interacting forces that can lead to variation among groups.

Differential Genome Size and Repetitive DNA Evolution in Diploid Species of Melampodium sect. Melampodium (Asteraceae)

Plant genomes vary greatly in composition and size mainly due to the diversity of repetitive DNAs and the inherent propensity for their amplification and removal from the host genome. Most studies addressing repeatome dynamics focus on model organisms, whereas few provide comprehensive investigations across the genomes of related taxa. Herein, we analyze the evolution of repeats of the 13 species in Melampodium sect. Melampodium, representing all but two of its diploid taxa, in a phylogenetic context. The investigated genomes range in size from 0.49 to 2.27 pg/1C (ca. 4.5-fold variation), despite having the same base chromosome number (x = 10) and very strong phylogenetic affinities. Phylogenetic analysis performed in BEAST and ancestral genome size reconstruction revealed mixed patterns of genome size increases and decreases across the group. High-throughput genome skimming and the RepeatExplorer pipeline were utilized to determine the repeat families responsible for the differences in observed genome sizes. Patterns of repeat evolution were found to be highly correlated with phylogenetic position, namely taxonomic series circumscription. Major differences found were in the abundances of the SIRE (Ty1-copia), Athila (Ty3-gypsy), and CACTA (DNA transposon) lineages. Additionally, several satellite DNA families were found to be highly group-specific, although their overall contribution to genome size variation was relatively small. Evolutionary changes in repetitive DNA composition and genome size were complex, with independent patterns of genome up- and downsizing throughout the evolution of the analyzed diploids. A model-based analysis of genome size and repetitive DNA composition revealed evidence for strong phylogenetic signal and differential evolutionary rates of major lineages of repeats in the diploid genomes.

Phylogeny, character evolution and spatiotemporal diversification of the species-rich and world-wide distributed tribe Rubieae (Rubiaceae)

The Rubiaceae tribe Rubieae has a world-wide distribution with up to 1,000 species. These collectively exhibit an enormous ecological and morphological diversity, making Rubieae an excellent group for macro- and microevolutionary studies. Previous molecular phylogenetic analyses used only a limited sampling within the tribe or missed lineages crucial for understanding character evolution in this group. Here, we analyze sequences from two plastid spacer regions as well as morphological and biogeographic data from an extensive and evenly distributed sampling to establish a sound phylogenetic framework. This framework serves as a basis for our investigation of the evolution of important morphological characters and the biogeographic history of the Rubieae. The tribe includes three major clades, the Kelloggiinae Clade (Kelloggia), the Rubiinae Clade (Didymaea, Rubia) and the most species-rich Galiinae Clade (Asperula, Callipeltis, Crucianella, Cruciata, Galium, Mericarpaea, Phuopsis, Sherardia, Valantia). Within the Galiinae Clade, the largest genera Galium and Asperula are para- and polyphyletic, respectively. Smaller clades, however, usually correspond to currently recognized taxa (small genera or sections within genera), which may be used as starting points for a refined classification in this clade. Life-form (perennial versus annual), flower shape (long versus short corolla tube) and fruit characters (dry versus fleshy, with or without uncinate hairs) are highly homoplasious and have changed multiple times independently. Inference on the evolution of leaf whorls, a characteristic feature of the tribe, is sensitive to model choice. Multi-parted leaf whorls appear to have originated from opposite leaves with two small interpetiolar stipules that are subsequently enlarged and increased in number. Early diversification of Rubieae probably started during the Miocene in western Eurasia. Disjunctions between the Old and the New World possibly are due to connections via a North Atlantic land bridge. Diversification of the Galiineae Clade started later in the Miocene, probably in the Mediterranean, from where lineages reached, often multiple times, Africa, eastern Asia and further on the Americas and Australia.

Dating the Species Network: Allopolyploidy and Repetitive DNA Evolution in American Daisies (Melampodium sect. Melampodium, Asteraceae)

Allopolyploidy has played an important role in the evolution of the flowering plants. Genome mergers are often accompanied by significant and rapid alterations of genome size and structure via chromosomal rearrangements and altered dynamics of tandem and dispersed repetitive DNA families. Recent developments in sequencing technologies and bioinformatic methods allow for a comprehensive investigation of the repetitive component of plant genomes. Interpretation of evolutionary dynamics following allopolyploidization requires both the knowledge of parentage and the age of origin of an allopolyploid. Whereas parentage is typically inferred from cytogenetic and phylogenetic data, age inference is hampered by the reticulate nature of the phylogenetic relationships. Treating subgenomes of allopolyploids as if they belonged to different species (i.e., no recombination among subgenomes) and applying cross-bracing (i.e., putting a constraint on the age difference of nodes pertaining to the same event), we can infer the age of allopolyploids within the framework of the multispecies coalescent within BEAST2. Together with a comprehensive characterization of the repetitive DNA fraction using the RepeatExplorer pipeline, we apply the dating approach in a group of closely related allopolyploids and their progenitor species in the plant genus Melampodium (Asteraceae). We dated the origin of both the allotetraploid, Melampodium strigosum, and its two allohexaploid derivatives, Melampodium pringlei and Melampodium sericeum, which share both parentage and the direction of the cross, to the Pleistocene ($<$1.4 Ma). Thus, Pleistocene climatic fluctuations may have triggered formation of allopolyploids possibly in short intervals, contributing to difficulties in inferring the precise temporal order of allopolyploid species divergence of M. sericeum and M. pringlei. The relatively recent origin of the allopolyploids likely played a role in the near-absence of major changes in the repetitive fraction of the polyploids' genomes. The repetitive elements most affected by the postpolyploidization changes represented retrotransposons of the Ty1-copia lineage Maximus and, to a lesser extent, also Athila elements of Ty3-gypsy family.