The nature of intraspecific and interspecific genome size variation in taxonomically complex eyebrights

Environmental impacts on intraspecific variation in Ambrosia artemisiifolia genome size in Slovakia, Central Europe

The quantity of DNA in angiosperms exhibits variation attributed to many external influences, such as environmental factors, geographical features, or stress factors, which exert constant selection pressure on organisms. Since invasive species possess adaptive capabilities to acclimate to novel environmental conditions, ragweed (Ambrosia artemisiifolia L.) was chosen as a subject for investigating their influence on genome size variation. Slovakia has diverse climatic conditions, suitable for testing the hypothesis that air temperature and precipitation, the main limiting factors of ragweed occurrence, would also have an impact on its genome size. Our results using flow cytometry confirmed this hypothesis and also found a significant association with geographical features such as latitude, altitude, and longitude. We can conclude that plants growing in colder environments farther from oceanic influences exhibit smaller DNA amounts, while optimal growth conditions result in a greater variability in genome size, reflecting the diminished effect of selection pressure.

Repeat turnover meets stable chromosomes: repetitive DNA sequences mark speciation and gene pool boundaries in sugar beet and wild beets

Sugar beet and its wild relatives share a base chromosome number of nine and similar chromosome morphologies. Yet, interspecific breeding is impeded by chromosome and sequence divergence that is still not fully understood. Since repetitive DNAs are among the fastest evolving parts of the genome, we investigated, if repeatome innovations and losses are linked to chromosomal differentiation and speciation. We traced genome and chromosome-wide evolution across 13 beet species comprising all sections of the genera Beta and Patellifolia. For this, we combined short and long read sequencing, flow cytometry, and cytogenetics to build a comprehensive framework that spans the complete scale from DNA to chromosome to genome. Genome sizes and repeat profiles reflect the separation into three gene pools with contrasting evolutionary patterns. Among all repeats, satellite DNAs harbor most genomic variability, leading to fundamentally different centromere architectures, ranging from chromosomal uniformity in Beta and Patellifolia to the formation of patchwork chromosomes in Corollinae/Nanae. We show that repetitive DNAs are causal for the genome expansions and contractions across the beet genera, providing insights into the genomic underpinnings of beet speciation. Satellite DNAs in particular vary considerably between beet genomes, leading to the evolution of distinct chromosomal setups in the three gene pools, likely contributing to the barriers in beet breeding. Thus, with their isokaryotypic chromosome sets, beet genomes present an ideal system for studying the link between repeats, genomic variability, and chromosomal differentiation and provide a theoretical fundament for understanding barriers in any crop breeding effort.

The global distribution of angiosperm genome size is shaped by climate

Angiosperms, which inhabit diverse environments across all continents, exhibit significant variation in genome sizes, making them an excellent model system for examining hypotheses about the global distribution of genome size. These include the previously proposed large genome constraint, mutational hazard, polyploidy-mediated, and climate-mediated hypotheses. We compiled the largest genome size dataset to date, encompassing 16 017 (> 5% of known) angiosperm species, and analyzed genome size distribution using a comprehensive geographic distribution dataset for all angiosperms. We observed that angiosperms with large range sizes generally had small genomes, supporting the large genome constraint hypothesis. Climate was shown to exert a strong influence on genome size distribution along the global latitudinal gradient, while the frequency of polyploidy and the type of growth form had negligible effects. In contrast to the unimodal patterns along the global latitudinal gradient shown by plant size traits and polyploid proportions, the increase in angiosperm genome size from the equator to 40-50°N/S is probably mediated by different (mostly climatic) mechanisms than the decrease in genome sizes observed from 40 to 50°N northward. Our analysis suggests that the global distribution of genome sizes in angiosperms is mainly shaped by climatically mediated purifying selection, genetic drift, relaxed selection, and environmental filtering.

LocoGSE, a sequence-based genome size estimator for plants

Extensive research has focused on exploring the range of genome sizes in eukaryotes, with a particular emphasis on land plants, where significant variability has been observed. Accurate estimation of genome size is essential for various research purposes, but existing sequence-based methods have limitations, particularly for low-coverage datasets. In this study, we introduce LocoGSE, a novel genome size estimator designed specifically for low-coverage datasets generated by genome skimming approaches. LocoGSE relies on mapping the reads on single copy consensus proteins without the need for a reference genome assembly. We calibrated LocoGSE using 430 low-coverage Angiosperm genome skimming datasets and compared its performance against other estimators. Our results demonstrate that LocoGSE accurately predicts monoploid genome size even at very low depth of coverage (<1X) and on highly heterozygous samples. Additionally, LocoGSE provides stable estimates across individuals with varying ploidy levels. LocoGSE fills a gap in sequence-based plant genome size estimation by offering a user-friendly and reliable tool that does not rely on high coverage or reference assemblies. We anticipate that LocoGSE will facilitate plant genome size analysis and contribute to evolutionary and ecological studies in the field. Furthermore, at the cost of an initial calibration, LocoGSE can be used in other lineages.

Comparative microscopic investigations of leaf epidermis in four Ajuga species from Korea

The genus Ajuga is widely distributed in temperate to subtropical regions, and four species are currently recognized in Korea (A. decumbens, A. multiflora, A. nipponensis, and A. spectabilis), but epidermal anatomical differences across these species have never been described. A comparative study of the leaf micromorphological characteristics of Korean Ajuga species was performed using light microscopy (LM) and scanning electron microscopy (SEM) to elucidate their taxonomic usefulness and to assess leaf micromorphological diversity. Considerable diversity in epidermal and stomatal anatomy was observed across Korean Ajuga species. Species had both hypostomatic or amphistomatic leaves, with anomocytic, anisocytic, diactyic, or actinocytic stomatal complexes. Guard cell length across species ranged from 17.66 ± 0.57 μm to 32.50 ± 2.38 μm and correlated with genome size. Abnormal stomata were frequently observed in three species (A. decumbens, A. multiflora, and A. nipponensis) but not in A. spectabilis. Three types of glandular trichomes were found: peltate in all species, short-stalked in all species, and long-stalked glandular trichomes in A. multiflora. Among the investigated leaf micromophological characters, trichome type, epidermal cell shape, and stomatal morphology were all taxonomically informative traits at a species level. RESEARCH HIGHLIGHTS: A comprehensive micromorphological description of the leaf surface is provided for Korean Ajuga species using scanning electron microscopic (SEM) and light microscopic (LM) analyses. The diverse range of stomatal development and the occurrence of polymorphic stomatal types are documented for the first time in Korean Ajuga species. The great diversity in stomatal and trichome morphology in Korean Ajuga species are taxonomically useful traits for species identification.

Is there hybridization between diploid and tetraploid Euphrasia in a secondary contact zone?

PREMISE

Strong postzygotic reproductive isolating barriers are usually expected to limit the extent of natural hybridization between species with contrasting ploidy. However, genomic sequencing has revealed previously overlooked examples of natural cross-ploidy hybridization in some flowering plant genera, suggesting that the phenomenon may be more common than once thought. We investigated potential cross-ploidy hybridization in British eyebrights (Euphrasia, Orobanchaceae), a group from which 13 putative cross-ploidy hybrid combinations have been reported based on morphology.

METHODS

We analyzed a contact zone between diploid Euphrasia rostkoviana and tetraploid E. arctica in Wales. We sequenced part of the internal transcribed spacer (ITS) of nuclear ribosomal DNA and used genotyping by sequencing (GBS) to look for evidence of cross-ploidy hybridization and introgression.

RESULTS

Common variant sites in the ITS region were fixed between diploids and tetraploids, indicating a strong barrier to hybridization. Clustering analyses of 356 single-nucleotide polymorphisms (SNPs) generated using GBS clearly separated samples by ploidy and revealed strong genetic structure (F_ST  = 0.44). However, the F_ST distribution across all SNPs was bimodal, indicating potential differential selection on loci between diploids and tetraploids. Demographic inference suggested potential gene flow, limited to around one or fewer migrants per generation.

CONCLUSIONS

Our results suggest that recent cross-ploidy hybridization is rare or absent in a site of secondary contact in Euphrasia. While a strong ploidy barrier prevents hybridization over ecological timescales, such hybrids may form in stable populations over evolutionary timescales, potentially allowing cross-ploidy introgression to take place.

Genetic Diversity Assessment of Iranian Kentucky Bluegrass Accessions: II. Nuclear DNA Content and Its Association with Morphological and Geographical Features

Poa pratensis L. is a perennial turfgrass with high regeneration and fertility, resistance to cold and drought, and quick colonization. By facultative apomixis, this plant can create a wide range of ploidy levels (2n = 22 to 2n = 154), resulting in a wide range of chromosomal numbers and sexual and apomictic reproductive diversity. The plant materials included fifty accessions from Iran's Center, South, North, North-East, North-West, and West ecoregions. UPOV standards were used to measure the qualities that were researched. The squash technique of chromosome counting revealed that Iranian Kentucky bluegrass accessions had chromosomal counts ranging from 24 to 87. The relative sizes of the 2C genomes were measured using laser flow cytometry. The range of DNA content was fairly wide, ranging from 4.92 to 11.52 pg. DNA content has a strong positive correlation with elevation, a moderately positive correlation with flag leaf length and leaf sheath width, and a negative correlation with inflorescence anthocyanin color and leaf anthocyanin color. The genotypes and ecological zones of this plant in Iran were distinguished based on morphological diversity and DNA content. The results from this study could be useful in identifying and studying wild Kentucky bluegrass genotypes. It aids in predicting the location of rare genotypes used as breeding materials. It can also increase the plant's variability for future generations by introducing new ecotypes, with particular genomic and morphological traits, to previously cultivated populations. We expect that the findings of this study will aid in understanding the evolution of this plant in the context of Iran's climatic variety.

Novel Insights into the Nature of Intraspecific Genome Size Diversity in Cannabis sativa L

Cannabis sativa has been used for millennia in traditional medicine for ritual purposes and for the production of food and fibres, thus, providing important and versatile services to humans. The species, which currently has a worldwide distribution, strikes out for displaying a huge morphological and chemical diversity. Differences in Cannabis genome size have also been found, suggesting it could be a useful character to differentiate between accessions. We used flow cytometry to investigate the extent of genome size diversity across 483 individuals belonging to 84 accessions, with a wide range of wild/feral, landrace, and cultivated accessions. We also carried out sex determination using the MADC2 marker and investigated the potential of flow cytometry as a method for early sex determination. All individuals were diploid, with genome sizes ranging from 1.810 up to 2.152 pg/2C (1.189-fold variation), apart from a triploid, with 2.884 pg/2C. Our results suggest that the geographical expansion of Cannabis and its domestication had little impact on its overall genome size. We found significant differences between the genome size of male and female individuals. Unfortunately, differences were, however, too small to be discriminated using flow cytometry through the direct processing of combined male and female individuals.

Integration of Repeatomic and Cytogenetic Data on Satellite DNA for the Genome Analysis in the Genus Salvia (Lamiaceae)

Within the complicated and controversial taxonomy of cosmopolitan genus Salvia L. (Lamiaceae) are valuable species Salvia officinalis L. and Salvia sclarea L., which are important for the pharmaceutical, ornamental horticulture, food, and perfume industries. Genome organization and chromosome structure of these essential oil species remain insufficiently studied. For the first time, the comparative repeatome analysis of S. officinalis and S. sclarea was performed using the obtained NGS data, RepeatExplorer/TAREAN pipelines and FISH-based chromosome mapping of the revealed satellite DNA families (satDNAs). In repeatomes of these species, LTR retrotransposons made up the majority of their repetitive DNA. Interspecific variations in genome abundance of Class I and Class II transposable elements, ribosomal DNA, and satellite DNA were revealed. Four (S. sclarea) and twelve (S. officinalis) putative satDNAs were identified. Based on patterns of chromosomal distribution of 45S rDNA; 5S rDNA and the revealed satDNAs, karyograms of S. officinalis and S. sclarea were constructed. Promising satDNAs which can be further used as chromosome markers to assess inter- and intraspecific chromosome variability in Salvia karyotypes were determined. The specific localization of homologous satDNA and 45S rDNA on chromosomes of the studied Salvia species confirmed their common origin, which is consistent with previously reported molecular phylogenetic data.

Measuring the Invisible: The Sequences Causal of Genome Size Differences in Eyebrights (Euphrasia) Revealed by k-mers

Genome size variation within plant taxa is due to presence/absence variation, which may affect low-copy sequences or genomic repeats of various frequency classes. However, identifying the sequences underpinning genome size variation is challenging because genome assemblies commonly contain collapsed representations of repetitive sequences and because genome skimming studies by design miss low-copy number sequences. Here, we take a novel approach based on k-mers, short sub-sequences of equal length k, generated from whole-genome sequencing data of diploid eyebrights (Euphrasia), a group of plants that have considerable genome size variation within a ploidy level. We compare k-mer inventories within and between closely related species, and quantify the contribution of different copy number classes to genome size differences. We further match high-copy number k-mers to specific repeat types as retrieved from the RepeatExplorer2 pipeline. We find genome size differences of up to 230Mbp, equivalent to more than 20% genome size variation. The largest contributions to these differences come from rDNA sequences, a 145-nt genomic satellite and a repeat associated with an Angela transposable element. We also find size differences in the low-copy number class (copy number ≤ 10×) of up to 27 Mbp, possibly indicating differences in gene space between our samples. We demonstrate that it is possible to pinpoint the sequences causing genome size variation within species without the use of a reference genome. Such sequences can serve as targets for future cytogenetic studies. We also show that studies of genome size variation should go beyond repeats if they aim to characterise the full range of genomic variants. To allow future work with other taxonomic groups, we share our k-mer analysis pipeline, which is straightforward to run, relying largely on standard GNU command line tools.

Altitudinal Patterns in Adaptive Evolution of Genome Size and Inter-Genome Hybridization Between Three Elymus Species From the Qinghai–Tibetan Plateau

Genome size variation and hybridization occur frequently within or between plant species under diverse environmental conditions, which enrich species diversification and drive the evolutionary process. Elymus L. is the largest genus in Triticeae with five recognized basic genomes (St, H, P, W, and Y). However, the data on population cytogenetics of Elymus species are sparse, especially whether genome hybridization and chromosomal structure can be affected by altitude are still unknown. In order to explore the relationship between genome sizes, we studied interspecific hybridization and altitude of Elymus species at population genetic and cytological levels. Twenty-seven populations at nine different altitudes (2,800–4,300 m) of three Elymus species, namely, hexaploid E. nutans (StHY, 2n = 6x = 42), tetraploid E. burchan-buddae (StY, 2n = 4x = 28), and E. sibiricus (StH, 2n = 4x = 28), were sampled from the Qinghai–Tibetan Plateau (QTP) to estimate whether intraspecific variation could affect the genomic relationships by genomic in situ hybridization (GISH), and quantify the genome size of Elymus among different altitude ecological groups by flow cytometry. The genome size of E. nutans, E. burchan-buddae, and E. sibiricus varied from 12.38 to 22.33, 8.81 to 18.93, and 11.46 to 20.96 pg/2C with the averages of 19.59, 12.39, and 16.85 pg/2C, respectively. The curve regression analysis revealed a strong correlation between altitude and nuclear DNA content in three Elymus species. In addition, the chromosomes of the St and Y genomes demonstrated higher polymorphism than that of the H genome. Larger genome size variations occurred in the mid-altitude populations (3,900–4,300 m) compared with other-altitude populations, suggesting a notable altitudinal pattern in genome size variation, which shaped genome evolution by altitude. This result supports our former hypothesis that genetic richness center at medium altitude is useful and valuable for species adaptation to highland environmental conditions, germplasm utilization, and conservation.

Pervasive Phylogenomic Incongruence Underlies Evolutionary Relationships in Eyebrights (Euphrasia, Orobanchaceae)

Disentangling the phylogenetic relationships of taxonomically complex plant groups is often mired by challenges associated with recent speciation, hybridization, complex mating systems, and polyploidy. Here, we perform a phylogenomic analysis of eyebrights (Euphrasia), a group renowned for taxonomic complexity, with the aim of documenting the extent of phylogenetic discordance at both deep and at shallow phylogenetic scales. We generate whole-genome sequencing data and integrate this with prior genomic data to perform a comprehensive analysis of nuclear genomic, nuclear ribosomal (nrDNA), and complete plastid genomes from 57 individuals representing 36 Euphrasia species. The species tree analysis of 3,454 conserved nuclear scaffolds (46 Mb) reveals that at shallow phylogenetic scales postglacial colonization of North Western Europe occurred in multiple waves from discrete source populations, with most species not being monophyletic, and instead combining genomic variants from across clades. At a deeper phylogenetic scale, the Euphrasia phylogeny is structured by geography and ploidy, and partially by taxonomy. Comparative analyses show Southern Hemisphere tetraploids include a distinct subgenome indicative of independent polyploidy events from Northern Hemisphere taxa. In contrast to the nuclear genome analyses, the plastid genome phylogeny reveals limited geographic structure, while the nrDNA phylogeny is informative of some geographic and taxonomic affinities but more thorough phylogenetic inference is impeded by the retention of ancestral polymorphisms in the polyploids. Overall our results reveal extensive phylogenetic discordance at both deeper and shallower nodes, with broad-scale geographic structure of genomic variation but a lack of definitive taxonomic signal. This suggests that Euphrasia species either have polytopic origins or are maintained by narrow genomic regions in the face of extensive homogenizing gene flow. Moreover, these results suggest genome skimming will not be an effective extended barcode to identify species in groups such as Euphrasia, or many other postglacial species groups.

Repeatome Analyses and Satellite DNA Chromosome Patterns in Deschampsia sukatschewii, D. cespitosa, and D. antarctica (Poaceae)

Subpolar and polar ecotypes of Deschampsia sukatschewii (Popl.) Roshev, D. cespitosa (L.) P. Beauv, and D. antarctica E. Desv. are well adapted to stressful environmental conditions, which make them useful model plants for genetic research and breeding. For the first time, the comparative repeatome analyses of subpolar and polar D. sukatschewii, D. cespitosa, and D. antarctica was performed using RepeatExplorer/TAREAN pipelines and FISH-based chromosomal mapping of the identified satellite DNA families (satDNAs). In the studied species, mobile genetic elements of class 1 made up the majority of their repetitive DNA; interspecific variations in the total amount of Ty3/Gypsy and Ty1/Copia retroelements, DNA transposons, ribosomal, and satellite DNA were revealed; 12–18 high confident and 7–9 low confident putative satDNAs were identified. According to BLAST, most D. sukatschewii satDNAs demonstrated sequence similarity with satDNAs of D. antarctica and D. cespitosa indicating their common origin. Chromosomal mapping of 45S rDNA, 5S rDNA, and satDNAs of D. sukatschewii allowed us to construct the species karyograms and detect new molecular chromosome markers important for Deschampsia species. Our findings confirmed that genomes of D. sukatschewii and D. cespitosa were more closely related compared to D. antarctica according to repeatome composition and patterns of satDNA chromosomal distribution.