Model adequacy tests for probabilistic models of chromosome-number evolution

ChromEvol v.3: modeling rate heterogeneity in chromosome number evolution

Changes in chromosome numbers are a prominent driver of plant evolution, impacting ecological diversification, stress tolerance, and phenotypes. ChromEvol is a widely used software tool for deciphering patterns of chromosome-number change along a phylogeny of interest. It evaluates the fit of alternative models to the data, estimates transition rates of different types of events, and infers the expected number of events along each branch of the phylogeny. We introduce ChromEvol v.3, featuring multiple novel methodological advancements that capture variation in the transition rates along a phylogeny. This version better allows researchers to identify how dysploidy and polyploidy rates change based on the number of chromosomes in the genome, with respect to a discrete trait, or at certain subclades of the phylogeny. We demonstrate the applicability of the new models on the Solanaceae phylogeny. Our analyses identify four chromosome-number transition regimes that characterize distinct Solanaceae clades and demonstrate an association between self-compatibility and altered dynamics of chromosome-number evolution. ChromEvol v.3, available at https://github.com/anatshafir1/chromevol, offers researchers a more flexible, comprehensive, and accurate tool to investigate the evolution of chromosome numbers and the various processes affecting it.

Shifts in Chromosome Evolution Rates Shape the Karyotype Patterns of Leafcutting Ants

Trait evolution has become a central focus in evolutionary biology, with phylogenetic comparative methods offering a framework to study how and why traits vary among species. Identifying variations in trait evolution rates within phylogenies is important for uncovering the mechanisms behind these differences. Karyotype variation, which is substantial across eukaryotic organisms, plays an essential role in species diversification. This study investigates karyotype variation within the leafcutting ant clade, focusing on chromosome number and morphology. We aim to determine whether karyotypic traits are phylogenetically dependent and how different evolutionary models explain karyotype diversity. Previous models have been insufficient in explaining these variations. To address these gaps, we employ modern phylogenetic methods to assess the impact of chromosomal fissions and fusions on karyotype evolution. By evaluating various evolutionary models-particularly the Brownian motion model, which suggests neutral chromosomal changes-we pursue for the further understanding the mode and tempo of karyotype evolution in ants. Our research examines how shifts in chromosomal change rates contribute to divergence among leafcutting ant species and assesses the role of chromosomal changes in the clade's evolutionary trajectory. Comparative analysis of leafcutting ant ideograms suggests that shared karyotype traits are strongly related to species relationships. This implies that karyotype diversification in leafcutting ants follows a phylogenetic trajectory at varying rates, with differences in karyotype traits reflecting the evolutionary distance between lineages. Particularly, the increase in the chromosome number of Acromyrmex is likely due to fission rearrangements rather than demi or polyploidization. We discuss and provide insights into the mechanisms driving karyotype variation and its implications for leafcutting ant diversification.

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.

Complex patterns of ploidy in a holocentric plant clade (Schoenus, Cyperaceae) in the Cape biodiversity hotspot

BACKGROUND AND AIMS

It is unclear how widespread polyploidy is throughout the largest holocentric plant family - the Cyperaceae. Because of the prevalence of chromosomal fusions and fissions, which affect chromosome number but not genome size, it can be impossible to distinguish if individual plants are polyploids in holocentric lineages based on chromosome count data alone. Furthermore, it is unclear how differences in genome size and ploidy levels relate to environmental correlates within holocentric lineages, such as the Cyperaceae.

METHODS

We focus our analyses on tribe Schoeneae, and more specifically the southern African clade of Schoenus. We examine broad-scale patterns of genome size evolution in tribe Schoeneae and focus more intensely on determining the prevalence of polyploidy across the southern African Schoenus by inferring ploidy level with the program ChromEvol, as well as interpreting chromosome number and genome size data. We further investigate whether there are relationships between genome size/ploidy level and environmental variables across the nutrient-poor and summer-arid Cape biodiversity hotspot.

KEY RESULTS

Our results show a large increase in genome size, but not chromosome number, within Schoenus compared to other species in tribe Schoeneae. Across Schoenus, there is a positive relationship between chromosome number and genome size, and our results suggest that polyploidy is a relatively common process throughout the southern African Schoenus. At the regional scale of the Cape, we show that polyploids are more often associated with drier locations that have more variation in precipitation between dry and wet months, but these results are sensitive to the classification of ploidy level.

CONCLUSIONS

Polyploidy is relatively common in the southern African Schoenus, where a positive relationship is observed between chromosome number and genome size. Thus, there may be a high incidence of polyploidy in holocentric plants, whose cell division properties differ from monocentrics.

Inferring Chromosome Number Changes Along a Phylogeny Using chromEvol

Chromosome numbers have long been used for the identification of key genomic events such as polyploidy and dysploidy. These inferences are often challenging, particularly when applied to large phylogenies, or clades in which more than a few chromosome number transitions had occurred. Here we describe the chromEvol computational framework that infers shifts in chromosome numbers along a phylogeny using probabilistic models of chromosome number change. Given chromosome count data and an associated phylogeny, chromEvol identifies such patterns by fitting probabilistic models of chromosome number evolution to the data. We describe the chromEvol workflow using available online tools, including the specification of the desired models, the examination of model fit to the data, and the inference of ploidy levels. The pipeline can be used by the wide scientific community and requires no previous computational or programming skills.

Relationship between fruit phenotypes and domestication in hexaploid populations of biribá (Annona mucosa) in Brazilian Amazonia

Background

Biribá (Annona mucosa Jacq.) is a fruit tree domesticated in Amazonia and has polyploid populations. The species presents ample phenotypic variation in fruit characteristics, including weight (100-4,000 g) and differences in carpel protrusions. Two cytotypes are recorded in the literature (2n = 28, 42) and genome size records are divergent (2C = 4.77, 5.42 and 6.00 pg). To decipher the role of polyploidy in the domestication of A. mucosa, we examined the relationships among phenotypic variation, chromosome number and genome size, and which came first, polyploidization or domestication.

Methodology

We performed chromosome counts of A. mucosa from central and western Brazilian Amazonia, and estimated genome size by flow cytometry. We performed phylogenetic reconstruction with publicly available data using a Bayesian framework, time divergence analysis and reconstructed the ancestral chromosome number for the genus Annona and for A. mucosa.

Results

We observed that variation in fruit phenotypes is not associated with variation in chromosome number and genome size. The most recent common ancestor of A. mucosa is inferred to be polyploid and diverged before domestication.

Conclusions

We conclude that, when domesticated, A. mucosa was already polyploid and we suggest that human selection is the main evolutionary force behind fruit size and fruit morphological variation in Annona mucosa.

The Chromosome Counts Database (CCDB)

For decades, plant biologists have been interested in the determination and documentation of chromosome numbers for extant taxa. This central cytological character has been used as an important phylogenetic marker and as an indicator for major genomic events such as polyploidy and dysploidy. Due to their significance and the relative ease by which chromosome numbers can be obtained, chromosome numbers have been extensively recorded across the plant kingdom and documented in a wide variety of resources. This makes the collection process a wearing task, often leading to partial data retrieval. In 2015, the Chromosome Counts Database (CCDB) was assembled, being an online unified community resource. This database compiles dozens of different chromosome counts sources, of which a significant portion had been unavailable before in a digitized, searchable format. The vast amount of data assembled in CCDB has already enabled a large number of analyses to examine the evolution of different plant hierarchies, as well as the application of various follow-up analyses, such as ploidy-level inference using chromEvol. CCDB ( http://ccdb.tau.ac.il/ ) encourages data sharing among the botanical community and is expected to continue expanding as additional chromosome numbers are recorded.

Simulation Tests of Methods in Evolution, Ecology, and Systematics: Pitfalls, Progress, and Principles

Complex statistical methods are continuously developed across the fields of ecology, evolution, and systematics (EES). These fields, however, lack standardized principles for evaluating methods, which has led to high variability in the rigor with which methods are tested, a lack of clarity regarding their limitations, and the potential for misapplication. In this review, we illustrate the common pitfalls of method evaluations in EES, the advantages of testing methods with simulated data, and best practices for method evaluations. We highlight the difference between method evaluation and validation and review how simulations, when appropriately designed, can refine the domain in which a method can be reliably applied. We also discuss the strengths and limitations of different evaluation metrics. The potential for misapplication of methods would be greatly reduced if funding agencies, reviewers, and journals required principled method evaluation.

Chromosome Number and Genome Size Evolution in Brasolia and Sobralia (Sobralieae, Orchidaceae)

Despite the clear circumscription of tribe Sobralieae (Orchidaceae), its internal relationships are still dubious. The recently delimited genus Brasolia, based on previous Sobralia species, is now assumed to be paraphyletic, with a third genus, Elleanthus, nested in it. The morphology of these three genera is significantly different, indicating the necessity of new data for a better genera delimitation. Though morphology and molecular data are available, cytogenetics data for Sobralieae is restricted to two Sobralia and one Elleanthus species. Aiming to evaluate the potential of cytogenetic data for Brasolia-Elleanthus-Sobralia genera delimitation, we present chromosome number and genome size data for 21 and 20 species, respectively, and used such data to infer the pattern of karyotype evolution in these genera. The analysis allowed us to infer x = 24 as the base chromosome number and genome size of average 1C-value of 5.0 pg for the common ancestor of Brasolia-Elleanthus-Sobralia. The recurrent descending dysploidy in Sobralieae and the punctual genome upsize suggest a recent diversification in Sobralieae but did not allow differing between Brasolia and Sobralia. However, the basal position of tribe Sobralieae in the subfamily Epidendroideae makes this tribe of interest to further studies clarifying the internal delimitation and pattern of karyotype evolution.

The Evolution of Cytogenetic Traits in Cuscuta (Convolvulaceae), the Genus With the Most Diverse Chromosomes in Angiosperms

Karyotypes are characterized by traits such as chromosome number, which can change through whole-genome duplication and dysploidy. In the parasitic plant genus Cuscuta (Convolvulaceae), chromosome numbers vary more than 18-fold. In addition, species of this group show the highest diversity in terms of genome size among angiosperms, as well as a wide variation in the number and distribution of 5S and 35S ribosomal DNA (rDNA) sites. To understand its karyotypic evolution, ancestral character state reconstructions were performed for chromosome number, genome size, and position of 5S and 35S rDNA sites. Previous cytogenetic data were reviewed and complemented with original chromosome counts, genome size estimates, and rDNA distribution assessed via fluorescence in situ hybridization (FISH), for two, seven, and 10 species, respectively. Starting from an ancestral chromosome number of x = 15, duplications were inferred as the prevalent evolutionary process. However, in holocentric clade (subgenus Cuscuta), dysploidy was identified as the main evolutionary mechanism, typical of holocentric karyotypes. The ancestral genome size of Cuscuta was inferred as approximately 1C = 12 Gbp, with an average genome size of 1C = 2.8 Gbp. This indicates an expansion of the genome size relative to other Convolvulaceae, which may be linked to the parasitic lifestyle of Cuscuta. Finally, the position of rDNA sites varied mostly in species with multiple sites in the same karyotype. This feature may be related to the amplification of rDNA sites in association to other repeats present in the heterochromatin. The data suggest that different mechanisms acted in different subgenera, generating the exceptional diversity of karyotypes in Cuscuta.