Genome size variation in wild and cultivated maize along altitudinal gradients

Biotic and abiotic factors associated with genome size evolution in oaks

The evolutionary processes that underlie variation in plant genome size have been much debated. Abiotic factors are thought to have played an important role, with negative and positive correlations between genome size and seasonal or stressful climatic conditions being reported in several systems. In turn, variation in genome size may influence plant traits which affect interactions with other organisms, such as herbivores. The mechanisms underlying evolutionary linkages between plant genome size and biotic and abiotic factors nonetheless remain poorly understod. To address this gap, we conducted phylogenetically controlled analyses testing for associations between genome size, climatic variables, plant traits (defenses and nutrients), and herbivory across 29 oak (Quercus) species. Genome size is significantly associated with both temperature and precipitation seasonality, whereby oak species growing in climates with lower and less variable temperatures but more variable rainfall had larger genomes. In addition, we found a negative association between genome size and leaf nutrient concentration (found to be the main predictor of herbivory), which in turn led to an indirect effect on herbivory. A follow-up test suggested that the association between genome size and leaf nutrients influencing herbivory was mediated by variation in plant growth, whereby species with larger genomes have slower growth rates, which in turn are correlated with lower nutrients. Collectively, these findings reveal novel associations between plant genome size and biotic and abiotic factors that may influence life history evolution and ecological dynamics in this widespread tree genus.

An allometry perspective on crops

Understanding trait-trait coordination is essential for successful plant breeding and crop modeling. Notably, plant size drives variation in morphological, physiological, and performance-related traits, as described by allometric laws in ecology. Yet, as allometric relationships have been limitedly studied in crops, how they influence and possibly limit crop performance remains unknown. Here, we review how an allometry perspective on crops gains insights into the phenotypic evolution during crop domestication, the breeding of varieties adapted to novel conditions, and the prediction of crop yields. As allometry is an active field of research, modeling and manipulating crop allometric relationships can help to develop more resilient and productive agricultural systems to face future challenges.

Genome size and chromosome number are critical metrics for accurate genome assembly assessment in Eukaryota

The number of genome assemblies has rapidly increased in recent history, with NCBI databases reaching over 41,000 eukaryotic genome assemblies across about 2,300 species. Increases in read length and improvements in assembly algorithms have led to increased contiguity and larger genome assemblies. While this number of assemblies is impressive, only about a third of these assemblies have corresponding genome size estimations for their respective species on publicly available databases. In this paper, genome assemblies are assessed regarding their total size compared to their respective publicly available genome size estimations. These deviations in size are assessed related to genome size, kingdom, sequencing platform, and standard assembly metrics, such as N50 and BUSCO values. A large proportion of assemblies deviate from their estimated genome size by more than 10%, with increasing deviations in size with increased genome size, suggesting nonprotein coding and structural DNA may be to blame. Modest differences in performance of sequencing platforms are noted as well. While standard metrics of genome assessment are more likely to indicate an assembly approaching the estimated genome size, much of the variation in this deviation in size is not explained with these raw metrics. A new, proportional N50 metric is proposed, in which N50 values are made relative to the average chromosome size of each species. This new metric has a stronger relationship with complete genome assemblies and, due to its proportional nature, allows for a more direct comparison across assemblies for genomes with variation in sizes and architectures.

Adaptive gene loss in the common bean pan-genome during range expansion and domestication

The common bean (Phaseolus vulgaris L.) is a crucial legume crop and an ideal evolutionary model to study adaptive diversity in wild and domesticated populations. Here, we present a common bean pan-genome based on five high-quality genomes and whole-genome reads representing 339 genotypes. It reveals ~234 Mb of additional sequences containing 6,905 protein-coding genes missing from the reference, constituting 49% of all presence/absence variants (PAVs). More non-synonymous mutations are found in PAVs than core genes, probably reflecting the lower effective population size of PAVs and fitness advantages due to the purging effect of gene loss. Our results suggest pan-genome shrinkage occurred during wild range expansion. Selection signatures provide evidence that partial or complete gene loss was a key adaptive genetic change in common bean populations with major implications for plant adaptation. The pan-genome is a valuable resource for food legume research and breeding for climate change mitigation and sustainable agriculture.

Diploid and tetraploid cytotypes of the flagship Cape species Dicerothamnus rhinocerotis (Asteraceae): variation in distribution, ecological niche, morphology and genetics

BACKGROUND AND AIMS

The Greater Cape Floristic Region is one of the world's biodiversity hotspots and is considered poor in polyploids. To test this assumption, ploidy variation was investigated in a widespread Cape shrub, Dicerothamnus rhinocerotis (renosterbos, Asteraceae). The aim was to elucidate the cytotype distribution and population composition across the species range, and to assess differences in morphology, environmental niches and genetics.

METHODS

Ploidy level and genome size were determined via flow cytometry and cytotype assignment was confirmed by chromosome counting. Restriction site-associated DNA sequencing (RADseq) analyses were used to infer genetic relationships. Cytotype climatic and environmental niches were compared using a range of environmental layers and a soil model, while morphological differences were examined using multivariate methods.

KEY RESULTS

The survey of 171 populations and 2370 individuals showed that the species comprises diploid and tetraploid cytotypes, no intermediates and only 16.8 % of mixed populations. Mean 2C values were 1.80-2.06 pg for diploids and 3.48-3.80 pg for tetraploids, with very similar monoploid genome sizes. Intra-cytotype variation showed a significant positive correlation with altitude and longitude in both cytotypes and with latitude in diploids. Although niches of both cytotypes were highly equivalent and similar, their optima and breadth were shifted due to differences mainly in isothermality and available water capacity. Morphometric analyses showed significant differences in the leaves and corolla traits, the number of florets per capitulum, and cypsela dimensions between the two cytotypes. Genetic analyses revealed four groups, three of them including both cytotypes.

CONCLUSIONS

Dicerothamnus rhinocerotis includes two distinct cytotypes that are genetically similar. While tetraploids arise several times independently within different genetic groups, morphological and ecological differences are evident between cytotypes. Our results open up new avenues for questions regarding the importance of ploidy in the megadiverse Cape flora, and exemplify the need for population-based studies focused on ploidy variation.

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.

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.

Genome size variation and evolution during invasive range expansion in an introduced plant

Plants demonstrate exceptional variation in genome size across species, and their genome sizes can also vary dramatically across individuals and populations within species. This aspect of genetic variation can have consequences for traits and fitness, but few studies attributed genome size differentiation to ecological and evolutionary processes. Biological invasions present particularly useful natural laboratories to infer selective agents that might drive genome size shifts across environments and population histories. Here, we test hypotheses for the evolutionary causes of genome size variation across 14 invading populations of yellow starthistle, Centaurea solstitialis, in California, United States. We use a survey of genome sizes and trait variation to ask: (1) Is variation in genome size associated with developmental trait variation? (2) Are genome sizes smaller toward the leading edge of the expansion, consistent with selection for "colonizer" traits? Or alternatively, does genome size increase toward the leading edge of the expansion, consistent with predicted consequences of founder effects and drift? (3) Finally, are genome sizes smaller at higher elevations, consistent with selection for shorter development times? We found that 2C DNA content varied 1.21-fold among all samples, and was associated with flowering time variation, such that plants with larger genomes reproduced later, with lower lifetime capitula production. Genome sizes increased toward the leading edge of the invasion, but tended to decrease at higher elevations, consistent with genetic drift during range expansion but potentially strong selection for smaller genomes and faster development time at higher elevations. These results demonstrate how genome size variation can contribute to traits directly tied to reproductive success, and how selection and drift can shape that variation. We highlight the influence of genome size on dynamics underlying a rapid range expansion in a highly problematic invasive plant.

Quantifying the role of genome size and repeat content in adaptive variation and the architecture of flowering time in Amaranthus tuberculatus

Genome size variation, largely driven by repeat content, is poorly understood within and among populations, limiting our understanding of its significance for adaptation. Here we characterize intraspecific variation in genome size and repeat content across 186 individuals of Amaranthus tuberculatus, a ubiquitous native weed that shows flowering time adaptation to climate across its range and in response to agriculture. Sequence-based genome size estimates vary by up to 20% across individuals, consistent with the considerable variability in the abundance of transposable elements, unknown repeats, and rDNAs across individuals. The additive effect of this variation has important phenotypic consequences-individuals with more repeats, and thus larger genomes, show slower flowering times and growth rates. However, compared to newly-characterized gene copy number and polygenic nucleotide changes underlying variation in flowering time, we show that genome size is a marginal contributor. Differences in flowering time are reflected by genome size variation across sexes and marginally, habitats, while polygenic variation and a gene copy number variant within the ATP synthesis pathway show consistently stronger environmental clines than genome size. Repeat content nonetheless shows non-neutral distributions across the genome, and across latitudinal and environmental gradients, demonstrating the numerous governing processes that in turn influence quantitative genetic variation for phenotypes key to plant adaptation.

Genetic basis of genome size variation of wheat

Research on various species has revealed a connection between genome size variation and the physiological and ecological characteristics of the species, suggesting that it could be a crucial factor influencing a species' adaptability to different environments. Wheat, being one of the world's three primary grains, holds significance in this regard. Investigating the genome size of wheat and analyzing the genetic factors contributing to its variation could offer valuable insights for enhancing wheat agronomic traits. This project has developed a conservative site ratio calculation approach to determine the size of the wheat genome. Additionally, it employs flow cytometry and k-mer distribution analysis to validate this method. Furthermore, the researchers use re-sequencing data to investigate the impact of environmental selection pressure and transposon dynamics on the variation in the size of the wheat genome. The findings from this study demonstrate a strong relationship between the size of the wheat genome and several environmental factors. These results serve as a valuable reference for understanding the development of variation in the size of the hetero-hexaploid wheat genome. Moreover, they contribute to advancing fundamental research on the genetic mechanisms underlying wheat characteristics. Additionally, the study paves the way for exploring new research directions in wheat breeding, which holds promise for future advancements in this field.

Genome size variation within Crithmum maritimum: Clues on the colonization of insular environments

Angiosperms present an astonishing diversity of genome sizes that can vary intra- or interspecifically. The remarkable new cytogenomic data shed some light on our understanding of evolution, but few studies were performed with insular and mainland populations to test possible correlations with dispersal, speciation, and adaptations to insular environments. Here, patterns of cytogenomic diversity were assessed among geographic samples (ca. 114) of Crithmum maritimum (Apiaceae), collected across the Azores and Madeira archipelagos, as well as in adjacent continental areas of Portugal. Using flow cytometry, the results indicated a significant intraspecific genome size variation, spanning from reduced sizes in the insular populations to larger ones in the mainland populations. Moreover, there was a tendency for an increase in genome size along the mainland populations, associated with lower temperatures, higher precipitation, and lower precipitation seasonality. However, this gradient might be the result of historic phylogeographical events associated with previous dispersal and extinction of local populations. Overall, our findings provided evidence that smaller genome sizes might play a critical role in the colonization of islands, corroborating other studies that argue that organisms with smaller genomes use fewer resources, having a selective advantage under insular environments. Although further studies are needed to improve our understanding of the mechanisms underlying genome size evolution on islands, conservation strategies must be promoted to protect the rich cytogenomic diversity found among C. maritimum populations, which occur in coastal areas that are particularly threatened by human activity, pollution, invasive species, and climate changes.

Gradual genome size evolution and polyploidy in Allium from the Qinghai-Tibetan Plateau

BACKGROUND AND AIMS

Genome size is an important plant trait, with substantial interspecies variation. The mechanisms and selective pressures underlying genome size evolution are important topics in evolutionary biology. There is considerable diversity in Allium from the Qinghai-Tibetan Plateau, where genome size variation and related evolutionary mechanisms are poorly understood.

METHODS

We reconstructed the Allium phylogeny using DNA sequences from 71 species. We also estimated genome sizes of 62 species, and determined chromosome numbers in 65 species. We examined the phylogenetic signal associated with genome size variation, and tested how well the data fit different evolutionary models. Correlations between genome size variations and seed mass, altitude and 19 bioclimatic factors were determined.

KEY RESULTS

Allium genome sizes differed substantially between species and within diploids, triploids, tetraploids, hexaploids and octaploids. Size per monoploid genome (1Cx) tended to decrease with increasing ploidy levels. Allium polyploids tended to grow at a higher altitude than diploids. The phylogenetic tree was divided into three evolutionary branches. The genomes in Clade I were mostly close to the ancestral genome (18.781 pg) while those in Clades II and III tended to expand and contract, respectively. A weak phylogenetic signal was detected for Allium genome size. Furthermore, significant positive correlations were detected between genome size and seed mass, as well as between genome size and altitude. However, genome size was not correlated with 19 bioclimatic variables.

CONCLUSIONS

Allium genome size shows gradual evolution, followed by subsequent adaptive radiation. The three well-supported Allium clades are consistent with previous studies. The evolutionary patterns in different Allium clades revealed genome contraction, expansion and relative stasis. The Allium species in Clade II may follow adaptive radiation. The genome contraction in Clade III may be due to DNA loss after polyploidization. Allium genome size might be influenced by selective pressure due to the conditions on the Qinghai-Tibetan Plateau (low temperature, high UV irradiation and abundant phosphate in the soil).

Analysis of Nuclear DNA Content and Karyotype of Phaseolus vulgaris L

The common bean (Phaseolus vulgaris L.), whose annual production is 26 million tons worldwide, is one of the main sources of protein and is known as one of the most important food sources. In this study, the karyotype variations and the genome size of four common bean genotypes in Turkey were investigated to determine whether the geographic variables in these regions affected the genome size and the karyotype parameters. In addition, it is known that as that the cytological and chromosomal parameters change under the influence of the climatic conditions of each region, appropriate and stable cytological methods for each plant facilitate and enable the determination of the chromosomal structure and the identification of specific chromosomes in the genotypes of the relevant region. Correct and valuable information such as this enables breeders and researchers to determine the correct shape and actual size of chromosomes. The genome size of the genotypes was measured with a flow cytometer, and chromosome analyses were performed with the squash method. For each genotype, the karyotype parameters, such as the number of somatic chromosomes, the Mean Total Chromosome Length (MTCL), the Mean Centromere Index (MCI), and the Mean Arm Ratio (MAR), were measured. The results showed that the highest and the lowest amounts of DNA per nucleus (3.28 pg and 1.49 pg) were observed in the Bitlis and Elaziğ genotypes. In addition, all genotype chromosome numbers were counted to be 2n = 2x = 22. The Mean Total Chromosome Length varied from 15.65 µm in Elaziğ to 34.24 µm in the Bitlis genotype. The Mean Chromosome Length ranged between 1.42 µm and 3.11 µm in the Elaziğ and Bitlis genotypes. The Hakkari and Van genotypes consist of eleven metacentric chromosomes, while the Bitlis and Elaziğ genotypes consist of ten metacentric chromosomes and one sub-metacentric chromosome. However, the Mean Centromere Index and Arm Ratio differed considerably among the genotypes. The highest (46.88) and the lowest (43.18) values of the Mean Centromere Index were observed in the Hakkari and Elaziğ genotypes, respectively. On the other hand, the lowest (1.15) and the highest (1.36) values of the Mean Arm Ratio were obtained in the Bitlis and Elaziğ genotypes, respectively. Eventually, intraspecies variations in genome size and chromosomal parameters were observed, and it was determined that the changes in nuclear DNA content and different chromosomal parameters among the four Phaseolus genotypes from four different regions of Turkey indicate the effect of climate change in the regions on these parameters. Such information in these areas can be used as useful information for the improvement of this plant and breeding programs.

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.

Combining genome size and pollen morphology data to study species relationships in the genus Daucus (Apiaceae)

BACKGROUND

The genus Daucus (Apiaceae) comprises about 40 wild species and the cultivated carrot, a crop of great economic and nutritional importance. The rich genetic diversity of wild Daucus species makes them a valuable gene pool for carrot improvement breeding programs. Therefore, it is essential to have good knowledge of the genome structure and relationships among wild Daucus species. To broaden such knowledge, in this research, the nuclear DNA content for 14 Daucus accessions and four closely related species was estimated by flow cytometry and their pollen morphology was analyzed by light and scanning electron microscopy (SEM).

RESULTS

The flow cytometric analysis showed a 3.2-fold variation in the mean 2C values among Daucus taxa, ranging from 0.999 (D. carota subsp. sativus) to 3.228 pg (D. littoralis). Among the outgroup species, the mean 2C values were 1.775-2.882 pg. The pollen grains of Daucus were tricolporate, mainly prolate or perprolate (rarely) in shape, and mainly medium or small (rarely) in size (21.19-40.38 µm), whereas the outgroup species had tricolporate, perprolate-shaped, and medium-sized (26.01-49.86 µm) pollen grains. In the studied taxa, SEM analysis revealed that exine ornamentation was striate, rugulate, perforate, or the ornamentation pattern was mixed. At the time of shedding, all pollen grains were three-celled, as evidenced by DAPI staining. We also found high positive correlations between the length of the polar axis (P) and the length of the equatorial diameter (E) of pollen grains, as well as between P and P/E. However, when comparing cytogenetic information with palynological data, no significant correlations were observed.

CONCLUSIONS

This study complements the information on the nuclear DNA content in Daucus and provides comprehensive knowledge of the pollen morphology of its taxa. These findings may be important in elucidating the taxonomic relationships among Daucus species and can help in the correct identification of gene bank accessions. In a broader view, they could also be meaningful for the interpretation of evolutionary trends in the genus.

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.

Are chromosome number and genome size associated with habit and environmental niche variables? Insights from the Neotropical orchids

BACKGROUND AND AIMS

The entangled relationship of chromosome number and genome size with species distribution has been the subject of study for almost a century, but remains an open question due to previous ecological and phylogenetic knowledge constraints. To better address this subject, we used the clade Maxillariinae, a widely distributed and karyotypically known orchid group, as a model system to infer such relationships in a robust methodological framework.

METHODS

Based on the literature and new data, we gathered the chromosome number and genome size for 93 and 64 species, respectively. We built a phylogenetic hypothesis and assessed the best macroevolutionary model for both genomic traits. Additionally, we collected together ecological data (preferences for bioclimatic variables, elevation and habit) used as explanatory variables in multivariate phylogenetic models explaining genomic traits. Finally, the impact of polyploidy was estimated by running the analyses with and without polyploids in the sample.

KEY RESULTS

The association between genomic and ecological data varied depending on whether polyploids were considered or not. Without polyploids, chromosome number failed to present consistent associations with ecological variables. With polyploids, there was a tendency to waive epiphytism and colonize new habitats outside humid forests. The genome size showed association with ecological variables: without polyploids, genome increase was associated with flexible habits, with higher elevation and with drier summers; with polyploids, genome size increase was associated with colonizing drier environments.

CONCLUSIONS

The chromosome number and genome size variations, essential but neglected traits in the ecological niche, are shaped in the Maxillariinae by both neutral and adaptive evolution. Both genomic traits are partially correlated to bioclimatic variables and elevation, even when controlling for phylogenetic constraints. While polyploidy was associated with shifts in the environmental niche, the genome size emerges as a central trait in orchid evolution by the association between small genome size and epiphytism, a key innovation to Neotropical orchid diversification.

New Insights Into the Evolution of C4 Photosynthesis Offered by the Tarenaya Cluster of Cleomaceae

Cleomaceae is closely related to Brassicaceae and includes C3, C3-C4, and C4 species. Thus, this family represents an interesting system for studying the evolution of the carbon concentrating mechanism. However, inadequate genetic information on Cleomaceae limits their research applications. Here, we characterized 22 Cleomaceae accessions [3 genera (Cleoserrata, Gynandropsis, and Tarenaya) and 11 species] in terms of genome size; molecular phylogeny; as well as anatomical, biochemical, and photosynthetic traits. We clustered the species into seven groups based on genome size. Interestingly, despite clear differences in genome size (2C, ranging from 0.55 to 1.3 pg) in Tarenaya spp., this variation was not consistent with phylogenetic grouping based on the internal transcribed spacer (ITS) marker, suggesting the occurrence of multiple polyploidy events within this genus. Moreover, only G. gynandra, which possesses a large nuclear genome, exhibited the C4 metabolism. Among the C3-like species, we observed intra- and interspecific variation in nuclear genome size as well as in biochemical, physiological, and anatomical traits. Furthermore, the C3-like species had increased venation density and bundle sheath cell size, compared to C4 species, which likely predisposed the former lineages to C4 photosynthesis. Accordingly, our findings demonstrate the potential of Cleomaceae, mainly members of Tarenaya, in offering novel insights into the evolution of C4 photosynthesis.

The arches and spandrels of maize domestication, adaptation, and improvement

People living in the Balsas River basin in southwest México domesticated maize from the bushy grass teosinte. Nine thousand years later, in 2021, Ms. Deb Haaland - a member of the Pueblo of Laguna tribe of New Mexico - wore a dress adorned with a cornstalk when she was sworn in as the Secretary of Interior of the United States of America. This choice of garment highlights the importance of the coevolution of maize and the farmers who, through careful selection over thousands of years, domesticated maize and adapted the physiology and shoot architecture of maize to fit local environments and growth habits. Some traits such as tillering were directly selected on (arches), and others such as tassel size are the by-products (spandrels) of maize evolution. Here, we review current knowledge of the underlying cellular, developmental, physiological, and metabolic processes that were selected by farmers and breeders, which have positioned maize as a top global staple crop.

The extensive amplification of heterochromatin in Melipona bees revealed by high throughput genomic and chromosomal analysis

Satellite DNAs (satDNAs) and transposable elements (TEs) are among the main components of constitutive heterochromatin (c-heterochromatin) and are related to their functionality, dynamics, and evolution. A peculiar case regarding the quantity and distribution of c-heterochromatin is observed in the genus of bees, Melipona, with species having a low amount of heterochromatin and species with high amount occupying almost all chromosomes. By combining low-pass genome sequencing and chromosomal analysis, we characterized the satDNAs and TEs of Melipona quadrifasciata (low c-heterochromatin) and Melipona scutellaris (high low c-heterochromatin) to understand c-heterochromatin composition and evolution. We identified 15 satDNA families and 20 TEs for both species. Significant variations in the repeat landscapes were observed between the species. In M. quadrifasciata, the repetitive fraction corresponded to only 3.78% of the genome library studied, whereas in M. scutellaris, it represented 54.95%. Massive quantitative and qualitative changes contributed to the differential amplification of c-heterochromatin, mainly due to the amplification of exclusive repetitions in M. scutellaris, as the satDNA MscuSat01-195 and the TE LTR/Gypsy_1 that represent 38.20 and 14.4% of its genome, respectively. The amplification of these two repeats is evident at the chromosomal level, with observation of their occurrence on most c-heterochromatin. Moreover, we detected repeats shared between species, revealing that they experienced mainly quantitative variations and varied in the organization on chromosomes and evolutionary patterns. Together, our data allow the discussion of patterns of evolution of repetitive DNAs and c-heterochromatin that occurred in a short period of time, after separation of the Michmelia and Melipona subgenera.

Intragenomic Conflict between Knob Heterochromatin and B Chromosomes Is the Key to Understand Genome Size Variation along Altitudinal Clines in Maize

In maize, we studied the causes of genome size variation and their correlates with cultivation altitude that suggests the existence of adaptive clines. To discuss the biological role of the genome size variation, we focused on Bolivian maize landraces growing along a broad altitudinal range. These were analyzed together with previously studied populations from altitudinal clines of Northwestern Argentina (NWA). Bolivian populations exhibited numerical polymorphism for B chromosomes (Bs) (from 1 to 5), with frequencies varying from 16.6 to 81.8 and being positively correlated with cultivation altitude. The 2C values of individuals 0B (A-DNA) ranged between 4.73 and 7.71 pg, with 58.33% of variation. The heterochromatic knobs, detected by DAPI staining, were more numerous and larger in individuals 0B than in those with higher doses of Bs. Bolivian and NWA landraces exhibited the same pattern of A-DNA downsizing and fewer and smaller knobs with increasing cultivation altitude, suggesting a mechanistic link among heterochromatin, genome size and phenology. The negative association between the two types of supernumerary DNA (knob heterochromatin and Bs), mainly responsible for the genome size variation, may be considered as an example of intragenomic conflict. It could be postulated that the optimal nucleotype is the result of such conflict, where genome adjustment may lead to an appropriate length of the vegetative cycle for maize landraces growing across altitudinal clines.

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

BACKGROUND AND AIMS

Genome size varies considerably across the diversity of plant life. Although genome size is, by definition, affected by genetic presence/absence variants, which are ubiquitous in population sequencing studies, genome size is often treated as an intrinsic property of a species. Here, we studied intra- and interspecific genome size variation in taxonomically complex British eyebrights (Euphrasia, Orobanchaceae). Our aim is to document genome size diversity and investigate underlying evolutionary processes shaping variation between individuals, populations and species.

METHODS

We generated genome size data for 192 individuals of diploid and tetraploid Euphrasia and analysed genome size variation in relation to ploidy, taxonomy, population affiliation and geography. We further compared the genomic repeat content of 30 samples.

KEY RESULTS

We found considerable intraspecific genome size variation, and observed isolation-by-distance for genome size in outcrossing diploids. Tetraploid Euphrasia showed contrasting patterns, with genome size increasing with latitude in outcrossing Euphrasia arctica, but with little genome size variation in the highly selfing Euphrasia micrantha. Interspecific differences in genome size and the genomic proportions of repeat sequences were small.

CONCLUSIONS

We show the utility of treating genome size as the outcome of polygenic variation. Like other types of genetic variation, such as single nucleotide polymorphisms, genome size variation may be affected by ongoing hybridization and the extent of population subdivision. In addition to selection on associated traits, genome size is predicted to be affected indirectly by selection due to pleiotropy of the underlying presence/absence variants.

Increased seminal root number associated with domestication improves nitrogen and phosphorus acquisition in maize seedlings

BACKGROUND AND AIMS

Domesticated maize (Zea mays ssp. mays) generally forms between two and six seminal roots, while its wild ancestor, Mexican annual teosinte (Zea mays ssp. parviglumis), typically lacks seminal roots. Maize also produces larger seeds than teosinte, and it generally has higher growth rates as a seedling. Maize was originally domesticated in the tropical soils of southern Mexico, but it was later brought to the Mexican highlands before spreading to other parts of the continent, where it experienced different soil resource constraints. The aims of this study were to understand the impacts of increased seminal root number on seedling nitrogen and phosphorus acquisition and to model how differences in maize and teosinte phenotypes might have contributed to increased seminal root number in domesticated maize.

METHODS

Seedling root system architectural models of a teosinte accession and a maize landrace were constructed by parameterizing the functional-structural plant model OpenSimRoot using plants grown in mesocosms. Seedling growth was simulated in a low-phosphorus environment, multiple low-nitrogen environments, and at variable planting densities. Models were also constructed to combine individual components of the maize and teosinte phenotypes.

KEY RESULTS

Seminal roots contributed ~35 % of the nitrogen and phosphorus acquired by maize landrace seedlings in the first 25 d after planting. Increased seminal root number improved plant nitrogen acquisition under low-nitrogen environments with varying precipitation patterns, fertilization rates, soil textures and planting densities. Models suggested that the optimal number of seminal roots for nutrient acquisition in teosinte is constrained by its limited seed carbohydrate reserves.

CONCLUSIONS

Seminal roots can improve the acquisition of both nitrogen and phosphorus in maize seedlings, and the increase in seed size associated with maize domestication may have facilitated increased seminal root number.

Using carrot centromeric repeats to study karyotype relationships in the genus Daucus (Apiaceae)

Background

In the course of evolution, chromosomes undergo evolutionary changes; thus, karyotypes may differ considerably among groups of organisms, even within closely related taxa. The genus Daucus seems to be a promising model for exploring the dynamics of karyotype evolution. It comprises some 40 wild species and the cultivated carrot, a crop of great economic significance. However, Daucus species are very diverse morphologically and genetically, and despite extensive research, the taxonomic and phylogenetic relationships between them have still not been fully resolved. Although several molecular cytogenetic studies have been conducted to investigate the chromosomal structure and karyotype evolution of carrot and other Daucus species, detailed karyomorphological research has been limited to carrot and only a few wild species. Therefore, to better understand the karyotype relationships within Daucus, we (1) explored the chromosomal distribution of carrot centromeric repeats (CentDc) in 34 accessions of Daucus and related species by means of fluorescence in situ hybridization (FISH) and (2) performed detailed karyomorphological analysis in 16 of them.

Results

We determined the genomic organization of CentDc in 26 accessions of Daucus (belonging to both Daucus I and II subclades) and one accession of closely related species. The CentDc repeats were present in the centromeric regions of all chromosomes of 20 accessions (representing 11 taxa). In the other Daucus taxa, the number of chromosome pairs with CentDc signals varied depending on the species, yet their centromeric localization was conserved. In addition, precise chromosome measurements performed in 16 accessions showed the inter- and intraspecific karyological relationships among them.

Conclusions

The presence of the CentDc repeats in the genomes of taxa belonging to both Daucus subclades and one outgroup species indicated the ancestral status of the repeat. The results of our study provide useful information for further evolutionary, cytotaxonomic, and phylogenetic research on the genus Daucus and may contribute to a better understanding of the dynamic evolution of centromeric satellites in plants.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07853-2.

Genome size evolution is associated with climate seasonality and glucosinolates, but not life history, soil nutrients or range size, across a clade of mustards

BACKGROUND AND AIMS

We investigate patterns of evolution of genome size across a morphologically and ecologically diverse clade of Brassicaceae, in relation to ecological and life history traits. While numerous hypotheses have been put forward regarding autecological and environmental factors that could favour small vs. large genomes, a challenge in understanding genome size evolution in plants is that many hypothesized selective agents are intercorrelated.

METHODS

We contribute genome size estimates for 47 species of Streptanthus Nutt. and close relatives, and take advantage of many data collections for this group to assemble data on climate, life history, soil affinity and composition, geographic range and plant secondary chemistry to identify simultaneous correlates of variation in genome size in an evolutionary framework. We assess models of evolution across clades and use phylogenetically informed analyses as well as model selection and information criteria approaches to identify variables that can best explain genome size variation in this clade.

KEY RESULTS

We find differences in genome size and heterogeneity in its rate of evolution across subclades of Streptanthus and close relatives. We show that clade-wide genome size is positively associated with climate seasonality and glucosinolate compounds. Model selection and information criteria approaches identify a best model that includes temperature seasonality and fraction of aliphatic glucosinolates, suggesting a possible role for genome size in climatic adaptation or a role for biotic interactions in shaping the evolution of genome size. We find no evidence supporting hypotheses of life history, range size or soil nutrients as forces shaping genome size in this system.

CONCLUSIONS

Our findings suggest climate seasonality and biotic interactions as potential forces shaping the evolution of genome size and highlight the importance of evaluating multiple factors in the context of phylogeny to understand the effect of possible selective agents on genome size.

Linking genome size variation to population phenotypic variation within the rotifer, Brachionus asplanchnoidis

Eukaryotic organisms usually contain much more genomic DNA than expected from their biological complexity. In explaining this pattern, selection-based hypotheses suggest that genome size evolves through selection acting on correlated life history traits, implicitly assuming the existence of phenotypic effects of (extra) genomic DNA that are independent of its information content. Here, we present conclusive evidence of such phenotypic effects within a well-mixed natural population that shows heritable variation in genome size. We found that genome size is positively correlated with body size, egg size, and embryonic development time in a population of the monogonont rotifer Brachionus asplanchnoidis. The effect on embryonic development time was mediated partly by an indirect effect (via egg size), and a direct effect, the latter indicating an increased replication cost of the larger amounts of DNA during mitosis. Our results suggest that selection-based change of genome size can operate in this population, provided it is strong enough to overcome drift or mutational change of genome size.

Chromosomal-scale genome assembly of Eleutherococcus senticosus provides insights into chromosome evolution in Araliaceae

Siberian ginseng (Eleutherococcus senticosus, also known as ciwujia) belongs to the Araliaceae family, which contains more than 1,500 species in 41 genera with diverse chromosome numbers and genome sizes. General consensus posits that ancient whole-genome duplication events and rapid evolutionary radiation are the driving forces for this variation in genome properties. In an attempt to generate more genomic information for the Araliaceae family, we report a 1.30 Gb high-quality draft genome assembly (contig N50 = 309.43 kb) of E. senticosus via PacBio long reads and Hi-C chromatin interaction maps. We found that transposable elements accounted for 72.82% of the genome and a total of 36,372 protein-coding genes were predicted. Comparative analyses of the E. senticosus, Panax notoginseng and Daucus carota genomes revealed a burst expansion of Tekay chromoviral elements in Araliaceae after its divergence with Apiaceae. We also found that E. senticosus underwent a lineage-specific whole-genome duplication event Es-α and a whole-genome duplication event Araliaceae-β that was probably shared by all Araliaceae species. Even though the rediploidization of the E. senticosus genome is evident, pathway analyses show that these two whole-genome duplication events may have contributed to the adaptation of E. senticosus to a cold environment. Taken together, the high-quality genome assembly of E. senticosus provides a valuable genomic resource for future research on the evolution of Araliaceae.

A phylogenetic analysis of Momordica (Cucurbitaceae) in India based on karyo-morphology, nuclear DNA content and rDNA ITS1-5.8S-ITS2 sequences

The infrageneric delimitation of Momordica, a medicinally important genus of Cucurbitaceae, is ill-defined until date. Momordica chromosomes are extremely small and are difficult to stain and visualize because of the dense cytoplasmic background. We have conducted karyomorphometric analysis by EMA method in five Indian Momordica species, and the nuclear genome sizes were estimated by flow cytometry for the first time. The somatic chromosome numbers ranged from 2n = 18 to 56 in the species. We have resolved previously disputed chromosome numbers in M. cymbalaria and M. dioica as 2n = 18 (lowest) and 2n = 56, respectively. Chromosome counts in the other species were re-confirmed as 2n = 22 in M. charantia, 2n = 28 in M. cochinchinensis and 2n = 56 in M. subangulata. The largest genome size was recorded in M. cymbalaria (3.74 pg 2C^-1), while the smallest size (0.72 pg 2C^-1) was detected in M. charantia var. charantia. The nuclear genome sizes were analysed in comparison to chromosome numbers and total chromosome lengths of the species. Karyomorphometric indices showed comparable symmetric karyotypes in the species except in M. cymbalaria having tendency towards asymmetry. The UPGMA phenogram and principle component analysis based on nuclear DNA contents and karyomorphometric parameters demonstrated interspecies differences, intraspecific distinction within M. charantia varieties and highlighted distinction of M. cymbalaria. This study was further supported by the rDNA ITS sequence-based phylogenetic analysis which revealed the monophyletic origin of the Indian members of Momordica and clarified the intraspecies relationship among the studied members. As a whole, the study brought out new insights on species diversification within the genus Momordica in India and would benefit further studies on biosystematics and plant breeding programmes.

Rapid Genome Evolution and Adaptation of Thlaspi arvense Mediated by Recurrent RNA-Based and Tandem Gene Duplications

Retrotransposons are the most abundant group of transposable elements (TEs) in plants, providing an extraordinarily versatile source of genetic variation. Thlaspi arvense, a close relative of the model plant Arabidopsis thaliana with worldwide distribution, thrives from sea level to above 4,000 m elevation in the Qinghai-Tibet Plateau (QTP), China. Its strong adaptability renders it an ideal model system for studying plant adaptation in extreme environments. However, how the retrotransposons affect the T. arvense genome evolution and adaptation is largely unknown. We report a high-quality chromosome-scale genome assembly of T. arvense with a scaffold N50 of 59.10 Mb. Long terminal repeat retrotransposons (LTR-RTs) account for 56.94% of the genome assembly, and the Gypsy superfamily is the most abundant TEs. The amplification of LTR-RTs in the last six million years primarily contributed to the genome size expansion in T. arvense. We identified 351 retrogenes and 303 genes flanked by LTRs, respectively. A comparative analysis showed that orthogroups containing those retrogenes and genes flanked by LTRs have a higher percentage of significantly expanded orthogroups (SEOs), and these SEOs possess more recent tandem duplicated genes. All present results indicate that RNA-based gene duplication (retroduplication) accelerated the subsequent tandem duplication of homologous genes resulting in family expansions, and these expanded gene families were implicated in plant growth, development, and stress responses, which were one of the pivotal factors for T. arvense's adaptation to the harsh environment in the QTP regions. In conclusion, the high-quality assembly of the T. arvense genome provides insights into the retroduplication mediated mechanism of plant adaptation to extreme environments.

Genome size and repetitive sequences are driven by artificial selection on the length of the vegetative cycle in maize landraces from Northeastern Argentina

Abstract Variation in genome size and knob heterochromatin content was explored in relationship to altitudinal cline and length of the vegetative cycle in northern Argentina, USA and Mexico landraces. It was considering that the decrease in DNA and heterochromatin content could be an adaptation to a shorter growing season and the result of artificial selection by man. Guaraní landraces from Northeastern Argentina (NEA) show similar variation in genome size (3.81pg to 7.56pg) and knob heterochromatin content than maize growing across an altitudinal cline. The present analysis offers an overview of the genetic variability of NEA maize to explain why Guaraní landraces and those along an altitudinal cline share this similar variation. Karyotype and flow cytometry data were employed. The DNA content of Guaraní landraces which lacking B chromosomes, showed no significant relationship with knob heterochromatin, suggesting differences in the amount of interspersed DNA. A significant positive relationship was found between the length of the vegetative cycle and both number and percentage of knob heterochromatin. No significant correlation was found between genome size and vegetative cycle. All these results allow us to conclude that the variation in heterochromatin content among Guaraní maize is driven by the selection of farmers for flowering time.

Evolutionary Genomics of Structural Variation in Asian Rice (Oryza sativa) Domestication

Structural variants (SVs) are a largely unstudied feature of plant genome evolution, despite the fact that SVs contribute substantially to phenotypes. In this study, we discovered SVs across a population sample of 347 high-coverage, resequenced genomes of Asian rice (Oryza sativa) and its wild ancestor (O. rufipogon). In addition to this short-read data set, we also inferred SVs from whole-genome assemblies and long-read data. Comparisons among data sets revealed different features of genome variability. For example, genome alignment identified a large (∼4.3 Mb) inversion in indica rice varieties relative to japonica varieties, and long-read analyses suggest that ∼9% of genes from the outgroup (O. longistaminata) are hemizygous. We focused, however, on the resequencing sample to investigate the population genomics of SVs. Clustering analyses with SVs recapitulated the rice cultivar groups that were also inferred from SNPs. However, the site-frequency spectrum of each SV type-which included inversions, duplications, deletions, translocations, and mobile element insertions-was skewed toward lower frequency variants than synonymous SNPs, suggesting that SVs may be predominantly deleterious. Among transposable elements, SINE and mariner insertions were found at especially low frequency. We also used SVs to study domestication by contrasting between rice and O. rufipogon. Cultivated genomes contained ∼25% more derived SVs and mobile element insertions than O. rufipogon, indicating that SVs contribute to the cost of domestication in rice. Peaks of SV divergence were enriched for known domestication genes, but we also detected hundreds of genes gained and lost during domestication, some of which were enriched for traits of agronomic interest.

Adaptive introgression from maize has facilitated the establishment of teosinte as a noxious weed in Europe

Global trade has considerably accelerated biological invasions. The annual tropical teosintes, the closest wild relatives of maize, were recently reported as new agricultural weeds in two European countries, Spain and France. Their prompt settlement under climatic conditions differing drastically from that of their native range indicates rapid genetic evolution. We performed a phenotypic comparison of French and Mexican teosintes under European conditions and showed that only the former could complete their life cycle during maize cropping season. To test the hypothesis that crop-to-wild introgression triggered such rapid adaptation, we used single nucleotide polymorphisms to characterize patterns of genetic variation in French, Spanish, and Mexican teosintes as well as in maize germplasm. We showed that both Spanish and French teosintes originated from Zea mays ssp. mexicana race "Chalco," a weedy teosinte from the Mexican highlands. However, introduced teosintes differed markedly from their Mexican source by elevated levels of genetic introgression from the high latitude Dent maize grown in Europe. We identified a clear signature of divergent selection in a region of chromosome 8 introgressed from maize and encompassing ZCN8, a major flowering time gene associated with adaptation to high latitudes. Moreover, herbicide assays and sequencing revealed that French teosintes have acquired herbicide resistance via the introgression of a mutant herbicide-target gene (ACC1) present in herbicide-resistant maize cultivars. Altogether, our results demonstrate that adaptive crop-to-wild introgression has triggered both rapid adaptation to a new climatic niche and acquisition of herbicide resistance, thereby fostering the establishment of an emerging noxious weed.

Genetic Distinctiveness Highlights the Conservation Value of a Sicilian Manna Ash Germplasm Collection Assigned to Fraxinus angustifolia (Oleaceae)

The cosmopolitan genus Fraxinus comprises about 40 species occupying several habitats in the Northern Hemisphere. With some species hybridizing and sharing genetic variants, questions remain on the species assignment of germplasm within the genus Fraxinus despite numerous species-specific assessments. A multidisciplinary approach was employed to provide a definitive insight into the genetics of an endangered Fraxinus “manna ash” collection, located in a rich plant biodiversity hotspot of the Madonie Mountains (Sicily). Although the collection size was small, genetic diversity, assessed by chloroplast (cpSSR) and nuclear (nSSR) microsatellites (SSR—Simple Sequence Repeats), allowed identifying three different chloroplast haplotypes, with one (H5) dominant, and several polymorphic loci, able to discriminate most of the local accessions studied. Molecular data were linked to cytofluorimetric and phenotypic evaluations and, contrary to popular belief that manna ash is Fraxinus ornus L., the germplasm currently used for manna production belongs to Fraxinus angustifolia Vahl. Interestingly, joint analysis of our genetic panel with a large European dataset of Fraxinus spp. suggested the presence of a possible glacial refuge in Sicily, confirming its importance as biodiversity source. Our results will be helpful for the design of long-term conservation programs for genetic resources, such as in situ and ex situ conservation, seed collection and tree reintroduction.

Repetitive sequences and structural chromosome alterations promote intraspecific variations in Zea mays L. karyotype

LTR-retrotransposons, knobs and structural chromosome alterations contribute to shape the structure and organization of the Zea mays karyotype. Our initial nuclear DNA content data of Z. mays accessions revealed an intraspecific variation (2 C = 2.00 pg to 2 C = 6.10 pg), suggesting differences in their karyotypes. We aimed to compare the karyotypes of three Z. mays accessions in search of the differences and similarities among them. Karyotype divergences were demonstrated among the accessions, despite their common chromosome number (2n = 20) and ancestral origin. Cytogenomic analyses showed that repetitive sequences and structural chromosome alterations play a significant role in promoting intraspecific nuclear DNA content variation. In addition, heterozygous terminal deletion in chromosome 3 was pointed out as a cause of lower nuclear 2 C value. Besides this, translocation was also observed in the short arm of chromosome 1. Differently, higher 2 C value was associated with the more abundant distribution of LTR-retrotransposons from the family Grande in the karyotype. Moreover, heteromorphism involving the number and position of the 180-bp knob sequence was found among the accessions. Taken together, we provide insights on the pivotal role played by repetitive sequences and structural chromosome alterations in shaping the karyotype of Z. mays.

The relevance of gene flow with wild relatives in understanding the domestication process

The widespread use of genomic tools has allowed for a deeper understanding of the genetics and the evolutionary dynamics of domestication. Recent studies have suggested that multiple domestications and introgression are more common than previously thought. However, the ability to correctly infer the many aspects of domestication process depends on having an adequate representation of wild relatives. Cultivated maize (Zea mays ssp. mays) is one of the most important crops in the world, with a long and a relatively well-documented history of domestication. The current consensus points towards a single domestication event from teosinte Zea mays ssp. parviglumis from the Balsas Basin in Southwestern Mexico. However, the underlying diversity of teosintes from Z. mays ssp. parviglumis and Zea mays ssp. mexicana was not taken into account in early studies. We used 32 739 single nucleotide polymorphisms (SNPs) obtained from 29 teosinte populations and 43 maize landraces to explore the relationship between wild and cultivated members of Zea. We then inferred the levels of gene flow among teosinte populations and maize, the degree of population structure of Zea mays subspecies, and the potential domestication location of maize. We confirmed a strong geographic structure within Z. mays ssp. parviglumis and documented multiple gene flow events with other members of the genus, including an event between Z. mays ssp. mexicana and maize. Our results suggest that the likely ancestor of maize may have been domesticated in Jalisco or in the southern Pacific Coast and not in the Balsas Basin as previously thought. In this context, different populations of both teosinte subspecies have contributed to modern maize's gene pool. Our results point towards a long period of domestication marked by gene flow with wild relatives, confirming domestication as long and ongoing process.

Within-Population Genome Size Variation is Mediated by Multiple Genomic Elements That Segregate Independently during Meiosis

Within-species variation in genome size has been documented in many animals and plants. Despite its importance for understanding eukaryotic genome diversity, there is only sparse knowledge about how individual-level processes mediate genome size variation in populations. Here, we study a natural population of the rotifer Brachionus asplanchnoidis whose members differ up to 1.9-fold in diploid genome size, but were still able to interbreed and produce viable offspring. We show that genome size is highly heritable and can be artificially selected up or down, but not below a certain basal diploid genome size for this species. Analyses of segregation patterns in haploid males reveal that large genomic elements (several megabases in size) provide the substrate of genome size variation. These elements, and their segregation patterns, explain the generation of new genome size variants, the short-term evolutionary potential of genome size change in populations, and some seemingly paradoxical patterns, like an increase in genome size variation among highly inbred lines. Our study suggests that a conceptual model involving only two variables, 1) a basal genome size of the population, and 2) a vector containing information on additional elements that may increase genome size in this population (size, number, and meiotic segregation behavior), can effectively address most scenarios of short-term evolutionary change of genome size in a population.

Common gardens in teosintes reveal the establishment of a syndrome of adaptation to altitude

In plants, local adaptation across species range is frequent. Yet, much has to be discovered on its environmental drivers, the underlying functional traits and their molecular determinants. Genome scans are popular to uncover outlier loci potentially involved in the genetic architecture of local adaptation, however links between outliers and phenotypic variation are rarely addressed. Here we focused on adaptation of teosinte populations along two elevation gradients in Mexico that display continuous environmental changes at a short geographical scale. We used two common gardens, and phenotyped 18 traits in 1664 plants from 11 populations of annual teosintes. In parallel, we genotyped these plants for 38 microsatellite markers as well as for 171 outlier single nucleotide polymorphisms (SNPs) that displayed excess of allele differentiation between pairs of lowland and highland populations and/or correlation with environmental variables. Our results revealed that phenotypic differentiation at 10 out of the 18 traits was driven by local selection. Trait covariation along the elevation gradient indicated that adaptation to altitude results from the assembly of multiple co-adapted traits into a complex syndrome: as elevation increases, plants flower earlier, produce less tillers, display lower stomata density and carry larger, longer and heavier grains. The proportion of outlier SNPs associating with phenotypic variation, however, largely depended on whether we considered a neutral structure with 5 genetic groups (73.7%) or 11 populations (13.5%), indicating that population stratification greatly affected our results. Finally, chromosomal inversions were enriched for both SNPs whose allele frequencies shifted along elevation as well as phenotypically-associated SNPs. Altogether, our results are consistent with the establishment of an altitudinal syndrome promoted by local selective forces in teosinte populations in spite of detectable gene flow. Because elevation mimics climate change through space, SNPs that we found underlying phenotypic variation at adaptive traits may be relevant for future maize breeding.

The genome-wide dynamics of purging during selfing in maize

Self-fertilization (also known as selfing) is an important reproductive strategy in plants and a widely applied tool for plant genetics and plant breeding. Selfing can lead to inbreeding depression by uncovering recessive deleterious variants, unless these variants are purged by selection. Here we investigated the dynamics of purging in a set of eleven maize lines that were selfed for six generations. We show that heterozygous, putatively deleterious single nucleotide polymorphisms are preferentially lost from the genome during selfing. Deleterious single nucleotide polymorphisms were lost more rapidly in regions of high recombination, presumably because recombination increases the efficacy of selection by uncoupling linked variants. Overall, heterozygosity decreased more slowly than expected, by an estimated 35% to 40% per generation instead of the expected 50%, perhaps reflecting pervasive associative overdominance. Finally, three lines exhibited marked decreases in genome size due to the purging of transposable elements. Genome loss was more likely to occur for lineages that began with larger genomes with more transposable elements and chromosomal knobs. These three lines purged an average of 398 Mb from their genomes, an amount equivalent to three Arabidopsis thaliana genomes per lineage, in only a few generations.

Novel Insights into Plant Genome Evolution and Adaptation as Revealed through Transposable Elements and Non-Coding RNAs in Conifers

Plant genomes are punctuated by repeated bouts of proliferation of transposable elements (TEs), and these mobile bursts are followed by silencing and decay of most of the newly inserted elements. As such, plant genomes reflect TE-related genome expansion and shrinkage. In general, these genome activities involve two mechanisms: small RNA-mediated epigenetic repression and long-term mutational decay and deletion, that is, genome-purging. Furthermore, the spatial relationships between TE insertions and genes are an important force in shaping gene regulatory networks, their downstream metabolic and physiological outputs, and thus their phenotypes. Such cascading regulations finally set up a fitness differential among individuals. This brief review demonstrates factual evidence that unifies most updated conceptual frameworks covering genome size, architecture, epigenetic reprogramming, and gene expression. It aims to give an overview of the impact that TEs may have on genome and adaptive evolution and to provide novel insights into addressing possible causes and consequences of intimidating genome sizes (20⁻30 Gb) in a taxonomic group, conifers.

Phylogenetic trends and environmental correlates of nuclear genome size variation in Helianthus sunflowers

Flowering plants serve as a powerful model for studying the evolution of nuclear genome size (GS) given the tremendous GS variation that exists both within and across angiosperm lineages. Helianthus sunflowers consist of c. 50 species native to North America that occupy diverse habitats and vary in ploidy level. In the current study, we generated a comprehensive GS database for 49 Helianthus species using flow cytometric approaches. We examined variability across the genus and present a comparative phylogenetic analysis of GS evolution in diploid Helianthus species. Results demonstrated that different clades of diploid Helianthus species showed evolutionary patterns of GS contraction, expansion and relative stasis, with annual diploid species evolving smaller GS with the highest rate of evolution. Phylogenetic comparative analyses of diploids revealed significant negative associations of GS with temperature seasonality and cell production rate, indicating that the evolution of larger GS in Helianthus diploids may be more permissible in habitats with longer growing seasons where selection for more rapid growth may be relaxed. The Helianthus GS database presented here and corresponding analyses of environmental and phenotypic correlates will facilitate ongoing and future research on the ultimate drivers of GS evolution in this well-studied North American plant genus.

Genome Size Unaffected by Variation in Morphological Traits, Temperature, and Precipitation in Turnip

Genome size (GS) was proposed as proxy for gross phenotypic and environmental changes in plants. GS organismal complexity is an enigma in evolutionary biology. While studies pertaining to intraspecific GS variation are abundant, literatures reporting the adaptive significance of GS are largelymissing. During food shortage, Brassica rapa var. rapa (turnip) is used as food and fodder for sustaining the livelihood of residents in the Qinghai Tibetan Plateau (QTP), which is also known as “the roof of the world”. Thus, climatic extremities make this region a natural environment to test adaptive significance of GS variation in turnip landraces. Therefore, from the QTP and its adjacent regions (the Hengduanshan and the Himalayas), we investigated adaptive evolution of GS in turnip landraces. Tuber diameter of turnip landraces was found to be significantly correlated with most of the environmental factors. GS was also shown not to be associated with morphological traits, temperature, and precipitation. Moreover, principal component analyses based on the whole dataset trisected the landraces into three distinct populations based on landrace usage—Hengduanshan, QTP, and the Himalayas. Nonetheless, our cumulative dataset showed evidence of adaptation of turnip landrace to different environments throughnonassociated genomic and phenomic plasticity.

Modeling Interactions between Transposable Elements and the Plant Epigenetic Response: A Surprising Reliance on Element Retention

Transposable elements (TEs) compose the majority of angiosperm DNA. Plants counteract TE activity by silencing them epigenetically. One form of epigenetic silencing requires 21-22 nt small interfering RNAs that act to degrade TE mRNA and may also trigger DNA methylation. DNA methylation is reinforced by a second mechanism, the RNA-dependent DNA methylation (RdDM) pathway. RdDM relies on 24 nt small interfering RNAs and ultimately establishes TEs in a quiescent state. These host factors interact at a systems level, but there have been no system level analyses of their interactions. Here, we define a deterministic model that represents the propagation of active TEs, aspects of the host response and the accumulation of silenced TEs. We describe general properties of the model and also fit it to biological data in order to explore two questions. The first is why two overlapping pathways are maintained, given that both are likely energetically expensive. Under our model, RdDM silenced TEs effectively even when the initiation of silencing was weak. This relationship implies that only a small amount of RNAi is needed to initiate TE silencing, but reinforcement by RdDM is necessary to efficiently counter TE propagation. Second, we investigated the reliance of the host response on rates of TE deletion. The model predicted that low levels of deletion lead to few active TEs, suggesting that silencing is most efficient when methylated TEs are retained in the genome, thereby providing one explanation for the large size of plant genomes.

Coregulation of ribosomal RNA with hundreds of genes contributes to phenotypic variation

Ribosomal repeats occupy 5% of a plant genome, yet there has been little study of their diversity in the modern age of genomics. Ribosomal copy number and expression variation present an opportunity to tap a novel source of diversity. In the present study, we estimated the ribosomal DNA (rDNA) copy number and ribosomal RNA (rRNA) expression for a population of maize inbred lines and investigated the potential role of rDNA and rRNA dosage in regulating global gene expression. Extensive variation was found in both ribosomal DNA copy number and ribosomal RNA expression among maize inbred lines. However, rRNA abundance was not consistent with the copy number of the rDNA. We have not found that the rDNA gene dosage has a regulatory role in gene expression; however, thousands of genes are identified to be coregulated with rRNA expression, including genes participating in ribosome biogenesis and other functionally relevant pathways. We further investigated the potential roles of copy number and the expression level of rDNA on agronomic traits and found that both correlated with flowering time but through different regulatory mechanisms. This comprehensive analysis suggested that rRNA expression variation is a valuable source of functional diversity that affects gene expression variation and field-based phenotypic changes.

Demography and its effects on genomic variation in crop domestication

Over two thousand plant species have been modified morphologically through cultivation and human use. Here, we review three aspects of crop domestication that are currently undergoing marked revisions, due to analytical advancements and their application to whole genome resequencing (WGS) data. We begin by discussing the duration and demographic history of domestication. There has been debate as to whether domestication occurred quickly or slowly. The latter is tentatively supported both by fossil data and application of WGS data to sequentially Markovian coalescent methods that infer the history of effective population size. This history suggests the possibility of extended human impacts on domesticated lineages prior to their purposeful cultivation. We also make the point that demographic history matters, because it shapes patterns and levels of extant genetic diversity. We illustrate this point by discussing the evolutionary processes that contribute to the empirical observation that most crops examined to date have more putatively deleterious alleles than their wild relatives. These deleterious alleles may contribute to genetic load within crops and may be fitting targets for crop improvement. Finally, the same demographic factors are likely to shape the spectrum of structural variants (SVs) within crops. SVs are known to underlie many of the phenotypic changes associated with domestication and crop improvement, but we currently lack sufficient knowledge about the mechanisms that create SVs, their rates of origin, their population frequencies and their phenotypic effects.

Parallel altitudinal clines reveal trends in adaptive evolution of genome size in Zea mays

While the vast majority of genome size variation in plants is due to differences in repetitive sequence, we know little about how selection acts on repeat content in natural populations. Here we investigate parallel changes in intraspecific genome size and repeat content of domesticated maize (Zea mays) landraces and their wild relative teosinte across altitudinal gradients in Mesoamerica and South America. We combine genotyping, low coverage whole-genome sequence data, and flow cytometry to test for evidence of selection on genome size and individual repeat abundance. We find that population structure alone cannot explain the observed variation, implying that clinal patterns of genome size are maintained by natural selection. Our modeling additionally provides evidence of selection on individual heterochromatic knob repeats, likely due to their large individual contribution to genome size. To better understand the phenotypes driving selection on genome size, we conducted a growth chamber experiment using a population of highland teosinte exhibiting extensive variation in genome size. We find weak support for a positive correlation between genome size and cell size, but stronger support for a negative correlation between genome size and the rate of cell production. Reanalyzing published data of cell counts in maize shoot apical meristems, we then identify a negative correlation between cell production rate and flowering time. Together, our data suggest a model in which variation in genome size is driven by natural selection on flowering time across altitudinal clines, connecting intraspecific variation in repetitive sequence to important differences in adaptive phenotypes.

Genome size in plants can vary by orders of magnitude, but this variation has long been considered to be of little functional consequence. Studying three independent adaptations to high altitude in Zea mays, we find that genome size experiences parallel pressures from natural selection, causing a reduction in genome size with increasing altitude. Though reductions in overall repetitive content are responsible for the genome size change, we find that only those individual loci contributing most to the variation in genome size are individually targeted by selection. To identify the phenotype influenced by genome size, we study how variation in genome size within a single wild population impacts leaf growth and cell division. We find that genome size variation correlates negatively with the rate of cell division, suggesting that individuals with larger genomes require longer to complete a mitotic cycle. Finally, we reanalyze data from maize inbreds to show that faster cell division is correlated with earlier flowering, connecting observed variation in genome size to an important adaptive phenotype.