Mechanisms of recent genome size variation in flowering plants

Conserved genome structure and phylogenetic insights for the heterogeneous subfamily of Convallarioideae (Asparagaceae) revealed from plastome data

BACKGROUND

Convallarioideae, a subfamily of Asparagaceae, encompasses a wide range of morphologically diverse lineages previously classified under different traditional families and holds significant economic value. Despite its importance, chloroplast genome data for Convallarioideae remain limited, hindering a comprehensive understanding of their genome structural evolution and phylogenetic relationships. This study aims to provide a detailed characterization of chloroplast genome features and to conduct robust phylogenetic analyses of this subfamily using an expanded dataset of chloroplast genomes.

RESULTS

The plastomes of the subfamily exhibit a typical circular quadripartite structure with conserved genomic organization and gene content. However, variations were observed in genome size, SSRs, and codon usage across the subfamily. Nine highly variable regions and positive selection genes were identified. Phylogenetic analyses based on complete plastid genomes resolved the non-monophyly of Polygonateae. Compared to Eriospermum mackenii, the chloroplast genomes of tribe Rusceae, tribe Dracaeneae, and the Polygonatum-Disporopsis lineage showed significant size reduction.

CONCLUSIONS

Chloroplast genomes across Convallarioideae exhibit remarkable structural conservation. The phylogenetic analyses revealed weakly resolved backbone relationships among core members of this subfamily. Indels in the LSC region and gene loss were identified as key drivers of structural divergence in plastome size. These results clarify the interplay between genomic architecture and phylogenetic discordance, advancing our understanding of genomic evolution within Convallarioideae.

Extensive genome evolution distinguishes maize within a stable tribe of grasses

Over the last 20 million years, the Andropogoneae tribe of grasses has evolved to dominate 17% of global land area. Domestication of these grasses in the last 10,000 years has yielded our most productive crops, including maize, sugarcane, and sorghum. The majority of Andropogoneae species, including maize, show a history of polyploidy - a condition that, while offering the evolutionary advantage of multiple gene copies, poses challenges to basic cellular processes, gene expression, and epigenetic regulation. Genomic studies of polyploidy have been limited by sparse sampling of taxa in groups with multiple polyploidy events. Here, we present 33 genome assemblies from 27 species, including chromosome-scale assemblies of maize relatives Zea and Tripsacum. In maize, the after-effects of polyploidy have been widely studied, showing reduced chromosome number, biased fractionation of duplicate genes, and transposable element (TE) expansions. While we observe these patterns within the genus Zea, 12 other polyploidy events deviate significantly. Those tetraploids and hexaploids retain elevated chromosome number, maintain nearly complete complements of duplicate genes, and have only stochastic TE amplifications. These genomes reveal variable outcomes of polyploidy, challenging simple predictions and providing a foundation for understanding its evolutionary implications in an ecologically and economically important clade.

New estimates and synthesis of chromosome numbers, ploidy levels and genome size variation in Allium sect. Codonoprasum: advancing our understanding of the unresolved diversification and evolution of this section

BACKGROUND

The genus Allium is known for its high chromosomal variability, but most chromosome counts are based on a few individuals and genome size (GS) reports are limited in certain taxonomic groups. This is evident in the Allium sect. Codonoprasum, a species-rich (> 150 species) and taxonomically complex section with weak morphological differences between taxa, the presence of polyploidy and frequent misidentification of taxa. Consequently, a significant proportion of older karyological reports may be unreliable and GS data are lacking for the majority of species within the section. This study, using chromosome counting and flow cytometry (FCM), provides the first comprehensive and detailed insight into variation in chromosome number, polyploid frequency and distribution, and GS in section members, marking a step towards understanding the unresolved diversification and evolution of this group.

RESULTS

We analysed 1578 individuals from 316 populations of 25 taxa and reported DNA ploidy levels and their GS, with calibration from chromosome counts in 22 taxa. Five taxa had multiple ploidy levels. First estimates of GS were obtained for 16 taxa. A comprehensive review of chromosome number and DNA-ploidy levels in 129 taxa of the section revealed that all taxa have x = 8, except A. rupestre with two polyploid series (x = 8, descending dysploidy x = 7), unique for this section. Diploid taxa dominated (72.1%), while di- & polyploid (12.4%) and exclusively polyploid (15.5%) taxa were less common. Ploidy diversity showed that diploid taxa dominated in the eastern Mediterranean and decreased towards the west and north, whereas only polyploid cytotypes of di- & polyploid taxa or exclusively polyploid taxa dominated in northern and northwestern Europe. A 4.1-fold variation in GS was observed across 33 taxa analysed so far (2C = 22.3-92.1 pg), mainly due to polyploidy, with GS downsizing observed in taxa with multiple ploidy levels. Intra-sectional GS variation suggests evolutionary relationships, and intraspecific GS variation within some taxa may indicate taxonomic heterogeneity and/or historical migration patterns.

CONCLUSIONS

Our study showed advantages of FCM as an effective tool for detecting ploidy levels and determining GS within the section. GS could be an additional character in understanding evolution and phylogenetic relationships within the section.

Exploring the complexity of genome size reduction in angiosperms

The genome sizes of angiosperms decreased significantly more than the genome sizes of their ancestors (pteridophytes and gymnosperms). Decreases in genome size involve a highly complex process, with remnants of the genome size reduction scattered across the genome and not directly linked to specific genomic structures. This is because the associated mechanisms operate on a much smaller scale than the mechanisms mediating increases in genome size. This review thoroughly summarizes the available literature regarding the molecular mechanisms underlying genome size reductions and introduces Utricularia gibba and Arabidopsis thaliana as model species for the examination of the effects of these molecular mechanisms. Additionally, we propose that phosphorus deficiency and drought stress are the major external factors contributing to decreases in genome size. Considering these factors affect almost all land plants, angiosperms likely gained the mechanisms for genome size reductions. These environmental factors may affect the retention rates of deletions, while also influencing the mutation rates of deletions via the functional diversification of the proteins facilitating double-strand break repair. The biased retention and mutation rates of deletions may have synergistic effects that enhance deletions in intergenic regions, introns, transposable elements, duplicates, and repeats, leading to a rapid decrease in genome size. We suggest that these selection pressures and associated molecular mechanisms may drive key changes in angiosperms during recurrent cycles of genome size decreases and increases.

Comparison of the Distribution Patterns of Microsatellites Across the Genomes of Reptiles

Microsatellites or simple sequence repeats (SSRs) are prevalent across various organisms' genomes. However, their distribution patterns and evolutionary dynamics in reptile genomes are rarely studied systematically. We herein conducted a comprehensive analysis of SSRs in the genomes of 36 reptile species. Our findings revealed that the total number of SSRs ranged from 1,840,965 to 7,664,452, accounting for 2.16%-8.19% of the genomes analyzed. The relative density ranged from 21,567.82 to 81,889.41 bp per megabase (Mbp). The abundance of different SSR categories followed the pattern of imperfect SSR (I-SSR) > perfect SSR (P-SSR) > compound SSR (C-SSR). A significant positive correlation was observed between the number of SSRs and genome size (p = 0.0034), whereas SSR frequency (p = 0.013) or density (p = 0.0099) showed a negative correlation with genome size. Furthermore, no correlation was found between SSR length and genome size. Mononucleotide repeats were the most common P-SSRs in crocodilians and turtles, whereas mononucleotides, trinucleotides, or tetranucleotides were the most common P-SSRs in snakes, lizards, and tuatara. P-SSRs of varying motif sizes showed nonrandom distribution across different genic regions, with AT-rich repeats being predominant. The genomic SSR content of the squamate lineage ranked the highest in abundance and variability, whereas crocodilians and turtles showed a slowly evolving and reduced microsatellite landscape. Gene ontology enrichment and Kyoto Encyclopedia of Genes and Genomes pathway analyses indicated that genes harboring P-SSRs in the coding DNA sequence regions were primarily involved in the regulation of transcription and translation processes. The SSR dataset generated in this study provides potential candidates for functional analysis and calls for broader-scale analyses across the evolutionary spectrum.

Comprehensive analysis of computational approaches in plant transcription factors binding regions discovery

Transcription factors (TFs) are regulatory proteins which bind to a specific DNA region known as the transcription factor binding regions (TFBRs) to regulate the rate of transcription process. The identification of TFBRs has been made possible by a number of experimental and computational techniques established during the past few years. The process of TFBR identification involves peak identification in the binding data, followed by the identification of motif characteristics. Using the same binding data attempts have been made to raise computational models to identify such binding regions which could save time and resources spent for binding experiments. These computational approaches depend a lot on what way they learn and how. These existing computational approaches are skewed heavily around human TFBRs discovery, while plants have drastically different genomic setup for regulation which these approaches have grossly ignored. Here, we provide a comprehensive study of the current state of the matters in plant specific TF discovery algorithms. While doing so, we encountered several software tools' issues rendering the tools not useable to researches. We fixed them and have also provided the corrected scripts for such tools. We expect this study to serve as a guide for better understanding of software tools' approaches for plant specific TFBRs discovery and the care to be taken while applying them, especially during cross-species applications. The corrected scripts of these software tools are made available at https://github.com/SCBB-LAB/Comparative-analysis-of-plant-TFBS-software.

Can DNA methylation shape climate response in trees?

Woody plants create the ecosystems they occupy and shape their biodiversity. Today's rapidly warming climate threatens these long-lived species by creating new environments in which existing populations become maladapted. Plants show enormous phenotypic diversity in response to environmental change, which can be caused by genotype or epigenetic mechanisms that influence the expression of the underlying DNA sequence. Whether epigenetics can affect ecologically important traits in trees is an important and controversial question. We explore the evidence that DNA methylation can affect gene expression, both directly and indirectly via its interaction with transposable elements (TEs), and subsequently shapes phenotypic variation in natural tree populations. Furthermore, we consider the potential of epigenetic approaches to assist in their conservation management strategies.

Herbgenomics meets Papaveraceae: a promising -omics perspective on medicinal plant research

Herbal medicines were widely used in ancient and modern societies as remedies for human ailments. Notably, the Papaveraceae family includes well-known species, such as Papaver somniferum and Chelidonium majus, which possess medicinal properties due to their latex content. Latex-bearing plants are a rich source of diverse bioactive compounds, with applications ranging from narcotics to analgesics and relaxants. With the advent of high-throughput technologies and advancements in sequencing tools, an opportunity exists to bridge the knowledge gap between the genetic information of herbs and the regulatory networks underlying their medicinal activities. This emerging discipline, known as herbgenomics, combines genomic information with other -omics studies to unravel the genetic foundations, including essential gene functions and secondary metabolite biosynthesis pathways. Furthermore, exploring the genomes of various medicinal plants enables the utilization of modern genetic manipulation techniques, such as Clustered Regularly-Interspaced Short Palindromic Repeats (CRISPR/Cas9) or RNA interference. This technological revolution has facilitated systematic studies of model herbs, targeted breeding of medicinal plants, the establishment of gene banks and the adoption of synthetic biology approaches. In this article, we provide a comprehensive overview of the recent advances in genomic, transcriptomic, proteomic and metabolomic research on species within the Papaveraceae family. Additionally, it briefly explores the potential applications and key opportunities offered by the -omics perspective in the pharmaceutical industry and the agrobiotechnology field.

First insight into the genomes of the Pulmonaria officinalis group (Boraginaceae) provided by repeatome analysis and comparative karyotyping

BACKGROUND

The genus Pulmonaria (Boraginaceae) represents a taxonomically complex group of species in which morphological similarity contrasts with striking karyological variation. The presence of different numbers of chromosomes in the diploid state suggests multiple hybridization/polyploidization events followed by chromosome rearrangements (dysploidy). Unfortunately, the phylogenetic relationships and evolution of the genome, have not yet been elucidated. Our study focused on the P. officinalis group, the most widespread species complex, which includes two morphologically similar species that differ in chromosome number, i.e. P. obscura (2n = 14) and P. officinalis (2n = 16). Ornamental cultivars, morphologically similar to P. officinalis (garden escapes), whose origin is unclear, were also studied. Here, we present a pilot study on genome size and repeatome dynamics of these closely related species in order to gain new information on their genome and chromosome structure.

RESULTS

Flow cytometry confirmed a significant difference in genome size between P. obscura and P. officinalis, corresponding to the number of chromosomes. Genome-wide repeatome analysis performed on genome skimming data showed that retrotransposons were the most abundant repeat type, with a higher proportion of Ty3/Gypsy elements, mainly represented by the Tekay lineage. Comparative analysis revealed no species-specific retrotransposons or striking differences in their copy number between the species. A new set of chromosome-specific cytogenetic markers, represented by satellite DNAs, showed that the chromosome structure in P. officinalis was more variable compared to that of P. obscura. Comparative karyotyping supported the hybrid origin of putative hybrids with 2n = 15 collected from a mixed population of both species and outlined the origin of ornamental garden escapes, presumably derived from the P. officinalis complex.

CONCLUSIONS

Large-scale genome size analysis and repeatome characterization of the two morphologically similar species of the P. officinalis group improved our knowledge of the genome dynamics and differences in the karyotype structure. A new set of chromosome-specific cytogenetic landmarks was identified and used to reveal the origin of putative hybrids and ornamental cultivars morphologically similar to P. officinalis.

Repeatome evolution across space and time: Unravelling repeats dynamics in the plant genus Erythrostemon Klotzsch (Leguminosae Juss)

Fluctuations in genomic repetitive fractions (repeatome) are known to impact several facets of evolution, such as ecological adaptation and speciation processes. Therefore, investigating the divergence of repetitive elements can provide insights into an important evolutionary force. However, it is not clear how the different repetitive element clades are impacted by the different factors such as ecological changes and/or phylogeny. To discuss this, we used the Neotropical legume genus Erythrostemon (Caesalpinioideae) as a model, given its ancient origin (~33 Mya), lineage-specific niche conservatism, macroecological heterogeneity, and disjunct distribution in Meso- and South American (MA and SA respectively) lineages. We performed a comparative repeatomic analysis of 18 Erythrostemon species to test the impact of environmental variables over repeats diversification. Overall, repeatome composition was diverse, with high abundances of satDNAs and Ty3/gypsy-Tekay transposable elements, predominantly in the MA and SA lineages respectively. However, unexpected repeatome profiles unrelated to the phylogeny/biogeography were found in a few MA (E. coccineus, E. pannosus and E. placidus) and SA (E. calycinus) species, related to reticulate evolution and incongruence between nuclear and plastid topology, suggesting ancient hybridizations. The plesiomorphic Tekay and satDNA pattern was altered in the MA-sensu stricto subclade with a striking genomic differentiation (expansion of satDNA and retraction of Tekay) associated with the colonization of a new environment in Central America around 20 Mya. Our data reveal that the current species-specific Tekay pool was the result of two bursts of amplification probably in the Miocene, with distinct patterns for the MA and SA repeatomes. This suggests a strong role of the Tekay elements as modulators of the genome-environment interaction in Erythrostemon, providing macroevolutionary insights about mechanisms of repeatome differentiation and plant diversification across space and time.

Transposable elements in Rosaceae: insights into genome evolution, expression dynamics, and syntenic gene regulation

Transposable elements (TEs) exert significant influence on plant genomic structure and gene expression. Here, we explored TE-related aspects across 14 Rosaceae genomes, investigating genomic distribution, transposition activity, expression patterns, and nearby differentially expressed genes (DEGs). Analyses unveiled distinct long terminal repeat retrotransposon (LTR-RT) evolutionary patterns, reflecting varied genome size changes among nine species over the past million years. In the past 2.5 million years, Rubus idaeus showed a transposition rate twice as fast as Fragaria vesca, while Pyrus bretschneideri displayed significantly faster transposition compared with Crataegus pinnatifida. Genes adjacent to recent TE insertions were linked to adversity resistance, while those near previous insertions were functionally enriched in morphogenesis, enzyme activity, and metabolic processes. Expression analysis revealed diverse responses of LTR-RTs to internal or external conditions. Furthermore, we identified 3695 pairs of syntenic DEGs proximal to TEs in Malus domestica cv. 'Gala' and M. domestica (GDDH13), suggesting TE insertions may contribute to varietal trait differences in these apple varieties. Our study across representative Rosaceae species underscores the pivotal role of TEs in plant genome evolution within this diverse family. It elucidates how these elements regulate syntenic DEGs on a genome-wide scale, offering insights into Rosaceae-specific genomic evolution.

PTFSpot: deep co-learning on transcription factors and their binding regions attains impeccable universality in plants

Unlike animals, variability in transcription factors (TFs) and their binding regions (TFBRs) across the plants species is a major problem that most of the existing TFBR finding software fail to tackle, rendering them hardly of any use. This limitation has resulted into underdevelopment of plant regulatory research and rampant use of Arabidopsis-like model species, generating misleading results. Here, we report a revolutionary transformers-based deep-learning approach, PTFSpot, which learns from TF structures and their binding regions' co-variability to bring a universal TF-DNA interaction model to detect TFBR with complete freedom from TF and species-specific models' limitations. During a series of extensive benchmarking studies over multiple experimentally validated data, it not only outperformed the existing software by >30% lead but also delivered consistently >90% accuracy even for those species and TF families that were never encountered during the model-building process. PTFSpot makes it possible now to accurately annotate TFBRs across any plant genome even in the total lack of any TF information, completely free from the bottlenecks of species and TF-specific models.

Ancient hybridization and repetitive element proliferation in the evolutionary history of the monocot genus Amomum (Zingiberaceae)

Genome size variation is a crucial aspect of plant evolution, influenced by a complex interplay of factors. Repetitive elements, which are fundamental components of genomic architecture, often play a role in genome expansion by selectively amplifying specific repeat motifs. This study focuses on Amomum, a genus in the ginger family (Zingiberaceae), known for its 4.4-fold variation in genome size. Using a robust methodology involving PhyloNet reconstruction, RepeatExplorer clustering, and repeat similarity-based phylogenetic network construction, we investigated the repeatome composition, analyzed repeat dynamics, and identified potential hybridization events within the genus. Our analysis confirmed the presence of four major infrageneric clades (A-D) within Amomum, with clades A-C exclusively comprising diploid species (2n = 48) and clade D encompassing both diploid and tetraploid species (2n = 48 and 96). We observed an increase in the repeat content within the genus, ranging from 84% to 89%, compared to outgroup species with 75% of the repeatome. The SIRE lineage of the Ty1-Copia repeat superfamily was prevalent in most analyzed ingroup genomes. We identified significant difference in repeatome structure between the basal Amomum clades (A, B, C) and the most diverged clade D. Our investigation revealed evidence of ancient hybridization events within Amomum, coinciding with a substantial proliferation of multiple repeat groups. This finding supports the hypothesis that ancient hybridization is a driving force in the genomic evolution of Amomum. Furthermore, we contextualize our findings within the broader context of genome size variations and repeatome dynamics observed across major monocot lineages. This study enhances our understanding of evolutionary processes within monocots by highlighting the crucial roles of repetitive elements in shaping genome size and suggesting the mechanisms that drive these changes.

Plastid genomes provide insights into the phylogeny and chloroplast evolution of the paper daisy tribe Gnaphalieae (Asteraceae)

Chloroplast genomes, as an essential source of phylogenetic information, are increasingly utilized in the evolutionary study of angiosperms. Gnaphalieae is a medium-sized tribe of the sunflower family of Asteraceae, with about 2,100 species in 178 genera distributed in temperate habitats worldwide. There has been considerable progress in our understanding of their phylogenetic evolution using both nuclear and chloroplast sequences, but no focus on chloroplast genomic data. In this study, we performed sequencing, assembly, and annotation of 16 representative chloroplast genomes from all the major lineages of Gnaphalieae. Our results showed that the plastomes exhibited a typical circular tetrad structure with similar genomic structure gene content. But there were differences in genome size, SSRs, and codon usage within the tribe. Phylogenetic analysis revealed Relhania clade is the earliest diverged lineages with the Lasiopogon clade and the Gnaphalium s.s. clade diverged subsequently. The core group includes FLAG clade sister to the HAP and Australasian group. Compared with the outgroup species, chloroplast genome size of the FLAG clade is much reduced whereas those of Australasian, HAP, Gnaphalium s.s., Lasiopogon and Relhania clades are relatively expanded. Insertions and deletions in the intergenic regions associated with repetitive sequence variations are supposed to be the main factor leading to length variations in the chloroplast genomes of Gnaphalieae. The comparative analyses of chloroplast genomes would provide useful implications into understanding the taxonomic and evolutionary history of Gnaphalieae.

Karyotype and LTR-RTs analysis provide insights into oak genomic evolution

BACKGROUND

Whole-genome duplication and long terminal repeat retrotransposons (LTR-RTs) amplification in organisms are essential factors that affect speciation, local adaptation, and diversification of organisms. Understanding the karyotype projection and LTR-RTs amplification could contribute to untangling evolutionary history. This study compared the karyotype and LTR-RTs evolution in the genomes of eight oaks, a dominant lineage in Northern Hemisphere forests.

RESULTS

Karyotype projections showed that chromosomal evolution was relatively conservative in oaks, especially on chromosomes 1 and 7. Modern oak chromosomes formed through multiple fusions, fissions, and rearrangements after an ancestral triplication event. Species-specific chromosomal rearrangements revealed fragments preserved through natural selection and adaptive evolution. A total of 441,449 full-length LTR-RTs were identified from eight oak genomes, and the number of LTR-RTs for oaks from section Cyclobalanopsis was larger than in other sections. Recent amplification of the species-specific LTR-RTs lineages resulted in significant variation in the abundance and composition of LTR-RTs among oaks. The LTR-RTs insertion suppresses gene expression, and the suppressed intensity in gene regions was larger than in promoter regions. Some centromere and rearrangement regions indicated high-density peaks of LTR/Copia and LTR/Gypsy. Different centromeric regional repeat units (32, 78, 79 bp) were detected on different Q. glauca chromosomes.

CONCLUSION

Chromosome fusions and arm exchanges contribute to the formation of oak karyotypes. The composition and abundance of LTR-RTs are affected by its recent amplification. LTR-RTs random retrotransposition suppresses gene expression and is enriched in centromere and chromosomal rearrangement regions. This study provides novel insights into the evolutionary history of oak karyotypes and the organization, amplification, and function of LTR-RTs.

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.

New whole-genome alignment tools are needed for tapping into plant diversity

Genome alignment is one of the most foundational methods for genome sequence studies. With rapid advances in sequencing and assembly technologies, these newly assembled genomes present challenges for alignment tools to meet the increased complexity and scale. Plant genome alignment is technologically challenging because of frequent whole-genome duplications (WGDs) as well as chromosome rearrangements and fractionation, high nucleotide diversity, widespread structural variation, and high transposable element (TE) activity causing large proportions of repeat elements. We summarize classical pairwise and multiple genome alignment (MGA) methods, and highlight techniques that are widely used or are being developed by the plant research community. We also outline the remaining challenges for precise genome alignment and the interpretation of alignment results in plants.

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.

The chromosome-level genome assembly of Fraxinus americana provides insights into the evolution of Oleaceae plants

White ash (Fraxinus americana linn.) originates from the southeastern United States. It is a tall and fast-growing tree species with strong salt-alkali resistance and cold tolerance, making it an important reforestation species and widely planted worldwide. Here, we completed the chromosome-level reference genome assembly of F. americana based on Illumina, PacBio, and Hi-C reads, with a genome size of 878.98 Mb, an N50 of 3.27 Mb, and a heterozygosity rate of 0.3 %. Based on de novo prediction, transcriptome prediction, and homology-based protein prediction, we obtained 39,538 genes. Approximately 843.21 Mb of the assembly genome was composed of 37,928 annotated protein-coding genes, with a gene function annotation rate of 95.93 %. 99.94 % of the overlap clusters (877.44 Mb) were anchored to 23 chromosomes. Synteny analysis of F. americana and other Oleaceae plants showed that F. americana underwent frequent chromosome rearrangements. The amplification of the Ale transposons effectively promoted the genome size of F. americana. Compared with other Oleaceae plants, the Glutathione S-transferase (GST) gene family in the F. americana genome has undergone significant expansion, which may help F. americana cope with adverse natural environments. Furthermore, we found that key enzyme-coding gene families related to lignin biosynthesis were expanded and highly expressed in F. americana leaves. These key genes drive lignin synthesis and benefit F. americana in fast-growing, as well as resisting biotic and abiotic stress. Overall, the F. americana genome assembly provides insights into the evolution of Oleaceae plants and provides abundant resources for breeding and germplasm conservation of white ash.

Combined analysis of transposable elements and structural variation in maize genomes reveals genome contraction outpaces expansion

Structural differences between genomes are a major source of genetic variation that contributes to phenotypic differences. Transposable elements, mobile genetic sequences capable of increasing their copy number and propagating themselves within genomes, can generate structural variation. However, their repetitive nature makes it difficult to characterize fine-scale differences in their presence at specific positions, limiting our understanding of their impact on genome variation. Domesticated maize is a particularly good system for exploring the impact of transposable element proliferation as over 70% of the genome is annotated as transposable elements. High-quality transposable element annotations were recently generated for de novo genome assemblies of 26 diverse inbred maize lines. We generated base-pair resolved pairwise alignments between the B73 maize reference genome and the remaining 25 inbred maize line assemblies. From this data, we classified transposable elements as either shared or polymorphic in a given pairwise comparison. Our analysis uncovered substantial structural variation between lines, representing both simple and complex connections between TEs and structural variants. Putative insertions in SNP depleted regions, which represent recently diverged identity by state blocks, suggest some TE families may still be active. However, our analysis reveals that within these recently diverged genomic regions, deletions of transposable elements likely account for more structural variation events and base pairs than insertions. These deletions are often large structural variants containing multiple transposable elements. Combined, our results highlight how transposable elements contribute to structural variation and demonstrate that deletion events are a major contributor to genomic differences.

A Flow Cytometry-Based Assessment of the Genomic Size and Ploidy Level of Wild Musa Species in India

The genome size variation is an important attribute in evolutionary and species characterization. Musa L. is regarded as one of the taxonomically complicated genera within the order Zingiberales, with more than 75 species from wild seeded to seedless cultivars that may be diploid, triploid or tetraploid. The knowledge of total nuclear DNA content in terms of genome size and ploidy level in wild species of Musa is absolutely important in evolutionary and genomic studies.

METHODS

In this paper, chromosome spreading was performed via protoplast isolation and a fast air-dry dropping method and flow cytometry were used with Raphanus sativus L. (Brassicaceae) as a standard for ploidy and genome size estimation.

RESULTS

The results showed that genome size (2C) varied amongst Musa species, based on the ratio of G1 peak positions. The lowest genome size (2C) was found in M. balbisiana var. andamanica (1.051 ± 0.060 pg) and the highest genome size (2C) was recorded for Musa ABB.cv. Meitei-hei (1.812 ± 0.108 pg) for the section Eumusa. Among the species belonging to the section Rhodochlamys, M. rosae had the lowest 2C content of 1.194 ± 0.033 pg whereas the highest nuclear DNA content (2C) was observed in M. velutina (1.488 ± 0.203 pg). Cytogenetic analysis revealed that the chromosome number of 14 wild Musa species was 2n = 22, while 1 species-Ensete glaucum-showed a chromosome number of 2n = 18 (diploid), and for 3 species, the chromosome number was 2n = 33 (triploids). An association study based on the Pearson correlation coefficient showed 2C nuclear DNA content was significant and positively correlated with ploidy level (R = 0.9) and chromosome number (R = 0.84).

CONCLUSIONS

The present study provides reliable information on the genome size and ploidy level of wild Musa species from the Indian region through flow cytometric analysis, which could be further utilized in taxonomic and crop improvement programs. For the first time, the nuclear DNA content of eight wild diploid and three triploid Indian species were estimated and reported. Genome size could be an effective indicator in identification of species and evolutionary studies in Musa with varying ploidy levels and morphological similarities.

Karyotype and genome size variation in Delphinium subg. Anthriscifolium (Ranunculaceae)

Five taxa of Delphiniumsubg.Anthriscifolium have been karyologically studied through chromosome counting, chromosomal measurement, and karyotype symmetry. Each taxon that we investigated has a basic chromosome number of x = 8, D.anthriscifoliumvar.savatieri, D.anthriscifoliumvar.majus, D.ecalcaratum, and D.callichromum were diploid with 2n = 16, while D.anthriscifoliumvar.anthriscifolium was tetraploid with 2n = 32. Monoploid chromosome sets of the investigated diploid taxa contained 1 metacentric chromosome, 3 submetacentric chromosomes, and 4 subtelocentric chromosomes. Higher interchromosomal asymmetry (CVCL) was present in D.ecalcaratum and D.callichromum than in other taxa. The highest levels of intrachromosomal asymmetry (MCA) and heterogeneity in centromere position (CVCI) were found in D.anthriscifoliumvar.majus. Diploid and tetraploid genome sizes varied by 3.02-3.92 pg and 6.04-6.60 pg, respectively. Karyotype and genome size of D.anthriscifoliumvar.savatieri, D.anthriscifoliumvar.majus, D.callichromum, and D.ecalcaratum were reported for the first time. Finally, based on cytological and morphological data, the classification of Delphiniumanthriscifolium was revised.

Chromosome-level reference genome assembly provides insights into the evolution of Pennisetum alopecuroides

Pennisetum alopecuroides is an important forage grass resource, which plays a vital role in ecological environment improvement. Therefore, the acquisition of P. alopecuroides genome resources is conducive to the study of the adaptability of Pennisetum species in ecological remediation and forage breeding development. Here we assembled a P. alopecuroides cv. 'Liqiu' genome at the chromosome level with a size of approximately 845.71 Mb, contig N50 of 84.83Mb, and genome integrity of 99.13% as assessed by CEGMA. A total of 833.41-Mb sequences were mounted on nine chromosomes by Hi-C technology. In total, 60.66% of the repetitive sequences and 34,312 genes were predicted. The genomic evolution analysis showed that P. alopecuroides cv. 'Liqiu' was isolated from Setaria 7.53-13.80 million years ago and from Cenchrus 5.33-8.99 million years ago, respectively. The whole-genome event analysis showed that P. alopecuroides cv. 'Liqiu' underwent two whole-genome duplication (WGD) events in the evolution process, and the duplication events occurred at a similar time to that of Oryza sativa and Setaria viridis. The completion of the genome sequencing of P. alopecuroides cv. 'Liqiu' provides data support for mining high-quality genetic resources of P. alopecuroides and provides a theoretical basis for the origin and evolutionary characteristics of Pennisetum.

Genome size and GC content of myxomycetes

More than 1272 myxomycetes species have been described, accounting for more than half of all Amoebozoa species. However, the genome size of only three myxomycetes species has been reported. Therefore, we used flow cytometry to present an extensive survey and a phylogeny-based analysis of genome size and GC content evolution in 144 myxomycetes species. The genome size of myxomycetes ranged from 18.7 Mb to 470.3 Mb, and the GC content ranged from 38.7% to 70.1%. Bright-spored clade showed larger genome sizes and more intra-order genome size variations than the dark-spored clade. GC content and genome size were positively correlated in both bright-spored and dark-spored clades, and spore size was positively correlated with genome size and GC content in the bright-spored clade. We provided the first genome size data set in Myxomycetes, and our results will provide helpful information for future Myxomycetes studies, such as genome sequencing.

Experimental estimates of germline mutation rate in eukaryotes: a phylogenetic meta-analysis

Mutation is the ultimate source of all genetic variation, and over the last 10 years the ready availability of whole-genome sequencing has permitted direct estimation of mutation rate for many non-model species across the tree of life. In this meta-analysis, we make a comprehensive search of the literature for mutation rate estimates in eukaryotes, identifying 140 mutation accumulation (MA) and parent-offspring (PO) sequencing studies covering 134 species. Based on these data, we revisit differences in the single-nucleotide mutation (SNM) rate between different phylogenetic lineages and update the known relationships between mutation rate and generation time, genome size, and nucleotide diversity-while accounting for phylogenetic nonindependence. We do not find a significant difference between MA and PO in estimated mutation rates, but we confirm that mammal and plant lineages have higher mutation rates than arthropods and that unicellular eukaryotes have the lowest mutation rates. We find that mutation rates are higher in species with longer generation times and larger genome sizes, even when accounting for phylogenetic relationships. Moreover, although nucleotide diversity is positively correlated with mutation rate, the gradient of the relationship is significantly less than one (on a logarithmic scale), consistent with higher mutation rates in populations with smaller effective size. For the 29 species for which data are available, we find that indel mutation rates are positively correlated with nucleotide mutation rates and that short deletions are generally more common than short insertions. Nevertheless, despite recent progress, no estimates of either SNM or indel mutation rates are available for the majority of deeply branching eukaryotic lineages-or even for most animal phyla. Even among charismatic megafauna, experimental mutation rate estimates remain unknown for amphibia and scarce for reptiles and fish.

Chromosomal-Scale Genome Assemblies of Two Coastal Plant Species, Scaevola taccada and S. hainanensis-Insight into Adaptation Outside of the Common Range

While most of the species in Goodeniaceae family, excluding the Scaevola genus, are endemic to Australasia, S. taccada and S. hainanensis have expanded their distribution range to the tropical coastlines of the Atlantic and Indian Oceans. S. taccada appears to be highly adapted to coastal sandy lands and cliffs, and it has become invasive in places. S. hainanensis is found mainly in salt marshes near mangrove forests, and is at risk of extinction. These two species provide a good system to investigate adaptive evolution outside the common distribution range of this taxonomic group. Here, we report their chromosomal-scale genome assemblies with the objective of probing their genomic mechanisms related to divergent adaptation after leaving Australasia. The scaffolds were assembled into eight chromosome-scale pseudomolecules, which covered 90.12% and 89.46% of the whole genome assembly for S. taccada and S. hainanensis, respectively. Interestingly, unlike many mangroves, neither species has undergone whole-genome duplication. We show that private genes, specifically copy-number expanded genes are essential for stress response, photosynthesis, and carbon fixation. The gene families that are expanded in S. hainanensis and contracted in S. taccada might have facilitated adaptation to high salinity in S. hainanensis. Moreover, the genes under positive selection in S. hainanensis have contributed to its response to stress and its tolerance of flooding and anoxic environments. In contrast, compared with S. hainanensis, the more drastic copy number expansion of FAR1 genes in S. taccada might have facilitated its adaptation to the stronger light radiation present in sandy coastal lands. In conclusion, our study of the chromosomal-scale genomes of S. taccada and S. hainanensis provides novel insights into their genomic evolution after leaving Australasia.

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).

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

Background

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

Methodology

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

Results

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

Conclusions

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

Retrotransposons: How the continuous evolutionary front shapes plant genomes for response to heat stress

Long terminal repeat retrotransposons (LTR retrotransposons) are the most abundant group of mobile genetic elements in eukaryotic genomes and are essential in organizing genomic architecture and phenotypic variations. The diverse families of retrotransposons are related to retroviruses. As retrotransposable elements are dispersed and ubiquitous, their "copy-out and paste-in" life cycle of replicative transposition leads to new genome insertions without the excision of the original element. The overall structure of retrotransposons and the domains responsible for the various phases of their replication is highly conserved in all eukaryotes. The two major superfamilies of LTR retrotransposons, Ty1/Copia and Ty3/Gypsy, are distinguished and dispersed across the chromosomes of higher plants. Members of these superfamilies can increase in copy number and are often activated by various biotic and abiotic stresses due to retrotransposition bursts. LTR retrotransposons are important drivers of species diversity and exhibit great variety in structure, size, and mechanisms of transposition, making them important putative actors in genome evolution. Additionally, LTR retrotransposons influence the gene expression patterns of adjacent genes by modulating potential small interfering RNA (siRNA) and RNA-directed DNA methylation (RdDM) pathways. Furthermore, comparative and evolutionary analysis of the most important crop genome sequences and advanced technologies have elucidated the epigenetics and structural and functional modifications driven by LTR retrotransposon during speciation. However, mechanistic insights into LTR retrotransposons remain obscure in plant development due to a lack of advancement in high throughput technologies. In this review, we focus on the key role of LTR retrotransposons response in plants during heat stress, the role of centromeric LTR retrotransposons, and the role of LTR retrotransposon markers in genome expression and evolution.

Transcriptome and Proteome Co-Profiling Offers an Understanding of Pre-Harvest Sprouting (PHS) Molecular Mechanisms in Wheat (Triticum aestivum)

While wheat (Triticum aestivum L.) is a widely grown and enjoyed crop, the diverse and complex global situation and climate are exacerbating the instability of its supply. In particular, pre-harvest sprouting (PHS) is one of the major abiotic stresses that frequently occurs due to irregular climate conditions, causing serious damage to wheat and its quality. In this study, transcriptomic analysis with RNA-seq and proteomic analysis with LC-MS/MS were performed in PHS-treated spikes from two wheat cultivars presenting PHS sensitivity and tolerance, respectively. A total of 13,154 differentially expressed genes (DEGs) and 706 differentially expressed proteins (DEPs) were identified in four comparison groups between the susceptible/tolerant cultivars. Gene function and correlation analysis were performed to determine the co-profiled genes and proteins affected by PHS treatment. In the functional annotation of each comparative group, similar functions were confirmed in each cultivar under PHS treatment; however, in Keumgang PHS+7 (K7) vs. Woori PHS+7 (W7), functional annotations presented clear differences in the "spliceosome" and "proteasome" pathways. In addition, our results indicate that alternative splicing and ubiquitin-proteasome support the regulation of germination and seed dormancy. This study provides an advanced understanding of the functions involved in transcription and translation related to PHS mechanisms, thus enabling specific proposals for the further analysis of germination and seed dormancy mechanisms and pathways in wheat.

Distinct composition and amplification dynamics of transposable elements in sacred lotus (Nelumbo nucifera Gaertn.)

Sacred lotus (Nelumbo nucifera Gaertn.) is a basal eudicot plant with a unique lifestyle, physiological features, and evolutionary characteristics. Here we report the unique profile of transposable elements (TEs) in the genome, using a manually curated repeat library. TEs account for 59% of the genome, and hAT (Ac/Ds) elements alone represent 8%, more than in any other known plant genome. About 18% of the lotus genome is comprised of Copia LTR retrotransposons, and over 25% of them are associated with non-canonical termini (non-TGCA). Such high abundance of non-canonical LTR retrotransposons has not been reported for any other organism. TEs are very abundant in genic regions, with retrotransposons enriched in introns and DNA transposons primarily in flanking regions of genes. The recent insertion of TEs in introns has led to significant intron size expansion, with a total of 200 Mb in the 28 455 genes. This is accompanied by declining TE activity in intergenic regions, suggesting distinct control efficacy of TE amplification in different genomic compartments. Despite the prevalence of TEs in genic regions, some genes are associated with fewer TEs, such as those involved in fruit ripening and stress responses. Other genes are enriched with TEs, and genes in epigenetic pathways are the most associated with TEs in introns, indicating a dynamic interaction between TEs and the host surveillance machinery. The dramatic differential abundance of TEs with genes involved in different biological processes as well as the variation of target preference of different TEs suggests the composition and activity of TEs influence the path of evolution.

Identification and characterization of transposable element AhMITE1 in the genomes of cultivated and two wild peanuts

Background

The cultivated peanut (Arachis hypogaea L., AABB) is an allotetraploid hybrid between two diploid peanuts, A. duranensis (AA genome) and A. ipaensis (BB genome). Miniature inverted-repeat transposable elements (MITEs), some of which are known as active nonautonomous DNA transposons with high copy numbers, play important roles in genome evolution and diversification. AhMITE1, a member of the MITE family of transposons, but information on the peanut genomes is still limited. Here, we analyzed AhMITE1, AuMITE1 and ApMITE1 in the cultivated (A. hypogaea) and two wild peanut (A. duranensis and A. ipaensis) genomes.

Results

The cultivated and the two wild peanut genomes harbored 142, 14 and 21 AhMITE1, AuMITE1 and ApMITE1 family members, respectively. These three family members exhibited highly conserved TIR sequences, and insertions preferentially occurred within 2 kb upstream and downstream of gene-coding and AT-rich regions. Phylogenetic and pairwise nucleotide diversity analysis showed that AhMITE1 and ApMITE1 family members have undergone one round of amplification bursts during the evolution of the peanut genome. PCR analyses were performed in 23 peanut varieties and demonstrated that AhMITE1 is an active transposon and that hybridization or chemical mutagenesis can promote the mobilization of AhMITE1.

Conclusions

AhMITE1, AuMITE1 and ApMITE1 family members were identified based on local BLAST search with MAK between the cultivated and the two wild peanut genomes. The phylogenetic, nucleotide diversity and variation copy numbers of AhMITE1, AuMITE1 and ApMITE1 members provides opportunities for investigating their roles during peanut evolution. These findings will contribute to knowledge on diversity of AhMITE1, provide information about the potential impact on the gene expression and promote the development of DNA markers in peanut.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08732-0.

Comparative Genomic and Phylogenetic Analysis of Chloroplast Genomes of Hawthorn (Crataegus spp.) in Southwest China

The hawthorns (Crataegus spp.) are widely distributed and famous for their edible and medicinal values. There are ∼18 species and seven varieties of hawthorn in China distributed throughout the country. We now report the chloroplast genome sequences from C. scabrifolia, C. chungtienensis and C. oresbia, from the southwest of China and compare them with the previously released six species in Crataegus and four species in Rosaceae. The chloroplast genome structure of Crataegus is typical and can be divided into four parts. The genome sizes are between 159,654 and 159,898bp. The three newly sequenced chloroplast genomes encode 132 genes, including 85 protein-coding genes, 37 tRNA genes, and eight rRNA genes. Comparative analysis of the chloroplast genomes revealed six divergent hotspot regions, including ndhA, rps16-trnQ-UUG, ndhF-rpl32, rps16-psbK, trnR-UCU-atpA and rpl32-trnL-UAG. According to the correlation and co-occurrence analysis of repeats with indels and SNPs, the relationship between them cannot be ignored. The phylogenetic tree constructed based on the complete chloroplast genome and intergenic region sequences indicated that C. scabrifolia has a different origin from C. chungtienensis and C. oresbia. We support the placement of C. hupehensis, C. cuneata, C. scabrifolia in C. subg. Crataegus and C. kansuensis, C. oresbia, C. kansuensis in C. subg. Sanguineae. In addition, based on the morphology, geographic distribution and phylogenetic relationships of C. chungtienensis and C. oresbia, we speculate that these two species may be the same species. In conclusion, this study has enriched the chloroplast genome resources of Crataegus and provided valuable information for the phylogeny and species identification of this genus.

Genome structure and evolutionary history of frankincense producing Boswellia sacra

Boswellia sacra Flueck (family Burseraceae) tree is wounded to produce frankincense. We report its de novo assembled genome (667.8 Mb) comprising 18,564 high-confidence protein-encoding genes. Comparing conserved single-copy genes across eudicots suggest >97% gene space assembly of B. sacra genome. Evolutionary history shows B. sacra gene-duplications derived from recent paralogous events and retained from ancient hexaploidy shared with other eudicots. The genome indicated a major expansion of Gypsy retroelements in last 2 million years. The B. sacra genetic diversity showed four clades intermixed with a primary genotype—dominating most resin-productive trees. Further, the stem transcriptome revealed that wounding concurrently activates phytohormones signaling, cell wall fortification, and resin terpenoid biosynthesis pathways leading to the synthesis of boswellic acid—a key chemotaxonomic marker of Boswellia. The sequence datasets reported here will serve as a foundation to investigate the genetic determinants of frankincense and other resin-producing species in Burseraceae.

•

Assembly and architecture of frankincense producing Boswellia sacra Flueck

•

Comparative genomics and evolutionary history of frankincense tree within orders

•

Transcriptome of stem part and gene expression patterns of wounding to the tree

•

Resin biosynthesis pathway and related CYP450 enzymes and gene families

Chromosome-level genome of Entada phaseoloides provides insights into genome evolution and triterpenoid saponins biosynthesis

As a medicinal herbal plant, Entada phaseoloides has high levels of secondary metabolites, particularly triterpenoid saponins, which are important resources for scientific research and medical applications. However, the lack of a reference genome for this genus has limited research on its evolution and utilization of its medicinal potential. In this study, we report a chromosome-scale genome assembly for E. phaseoloides using Illumina, Nanopore long reads, and high-throughput chromosome conformation capture technology. The assembled reference genome is 456.18 Mb (scaffold N50=30.9 Mb; contig N50=6.34 Mb) with 95.71 % of the sequences anchored onto 14 pseudochromosomes. E. phaseoloides was estimated to diverge from the Leguminosae lineage at approximately 72.0 million years ago. With the integration of transcriptomic and metabolomic data, gene expression patterns and metabolite profiling of E. phaseoloides were determined in different tissues. The pattern of gene expression and metabolic profile of the kernel were distinct from those of other tissues. Furthermore, the evolution of certain gene families involved in the biosynthesis of triterpenoid saponins and terpenes was analyzed and offer new insights into the formation of these two metabolites. Four CYP genes, one UGT gene and related transcription factors were identified as candidate genes contributing to regulation of triterpenoid saponins biosynthesis. As the first high-quality assembled reference genome in the genus Entada, it will not only provide new information for the evolutionary study of this genus and conservation biology of E. phaseoloides but also lay a foundation for the formation and utilization of secondary metabolites in medicinal plants.

Complete Chloroplast Genomes Provide Insights Into Evolution and Phylogeny of Campylotropis (Fabaceae)

The genus Campylotropis Bunge (Desmodieae, Papilionoideae) comprises about 37 species distributed in temperate and tropical Asia. Despite the great potential in soil conservation, horticulture, and medicine usage, little is known about the evolutionary history and phylogenetic relationships of Campylotropis due to insufficient genetic resources. Here, we sequenced and assembled 21 complete chloroplast genomes of Campylotropis species. In combination with the previously published chloroplast genomes of C. macrocarpa and closely related species, we conducted comparative genomics and phylogenomic analysis on these data. Comparative analysis of the genome size, structure, expansion and contraction of inverted repeat (IR) boundaries, number of genes, GC content, and pattern of simple sequence repeats (SSRs) revealed high similarities among the Campylotropis chloroplast genomes. The activities of long sequence repeats contributed to the variation in genome size and gene content in Campylotropis chloroplast genomes. The Campylotropis chloroplast genomes showed moderate sequence variation, and 13 highly variable regions were identified for species identification and further phylogenetic studies. We also reported one more case of matK pseudogene in the legume family. The phylogenetic analysis confirmed the monophyly of Campylotropis and the sister relationship between Lespedeza and Kummerowia, the latter two genera were then sister to Campylotropis. The intrageneric relationships of Campylotropis based on genomic scale data were firstly reported in this study. The two positively selected genes (atpF and rps19) and eight fast-evolving genes identified in this study may help us to understand the adaptation of Campylotropis species. Overall, this study enhances our understanding of the chloroplast genome evolution and phylogenetic relationships of Campylotropis.

Chromosome-level genome of Camellia lanceoleosa provides a valuable resource for understanding genome evolution and self-incompatibility

The section Oleifera (Theaceae) has attracted attention for the high levels of unsaturated fatty acids found in its seeds. Here, we report the chromosome-scale genome of the sect. Oleifera using diploid wild Camellia lanceoleosa with a final size of 3.00 Gb and an N50 scaffold size of 186.43 Mb. Repetitive sequences accounted for 80.63% and were distributed unevenly across the genome. Camellia lanceoleosa underwent a whole-genome duplication event approximately 65 million years ago (65 Mya), prior to the divergence of C. lanceoleosa and Camellia sinensis (approx. 6-7 Mya). Syntenic comparisons of these two species elucidated the genomic rearrangement, appearing to be driven in part by the activity of transposable elements. The expanded and positively selected genes in C. lanceoleosa were significantly enriched in oil biosynthesis, and the expansion of homomeric acetyl-coenzyme A carboxylase (ACCase) genes and the seed-biased expression of genes encoding heteromeric ACCase, diacylglycerol acyltransferase, glyceraldehyde-3-phosphate dehydrogenase and stearoyl-ACP desaturase could be of primary importance for the high oil and oleic acid content found in C. lanceoleosa. Theanine and catechins were present in the leaves of C. lanceoleosa. However, caffeine can not be dectected in the leaves but was abundant in the seeds and roots. The functional and transcriptional divergence of genes encoding SAM-dependent N-methyltransferases may be associated with caffeine accumulation and distribution. Gene expression profiles, structural composition and chromosomal location suggest that the late-acting self-incompatibility of C. lanceoleosa is likely to have favoured a novel mechanism co-occurring with gametophytic self-incompatibility. This study provides valuable resources for quantitative and qualitative improvements and genome assembly of polyploid plants in sect. Oleifera.

AnchorWave: Sensitive alignment of genomes with high sequence diversity, extensive structural polymorphism, and whole-genome duplication

Millions of species are currently being sequenced, and their genomes are being compared. Many of them have more complex genomes than model systems and raise novel challenges for genome alignment. Widely used local alignment strategies often produce limited or incongruous results when applied to genomes with dispersed repeats, long indels, and highly diverse sequences. Moreover, alignment using many-to-many or reciprocal best hit approaches conflicts with well-studied patterns between species with different rounds of whole-genome duplication. Here, we introduce Anchored Wavefront alignment (AnchorWave), which performs whole-genome duplication-informed collinear anchor identification between genomes and performs base pair-resolved global alignment for collinear blocks using a two-piece affine gap cost strategy. This strategy enables AnchorWave to precisely identify multikilobase indels generated by transposable element (TE) presence/absence variants (PAVs). When aligning two maize genomes, AnchorWave successfully recalled 87% of previously reported TE PAVs. By contrast, other genome alignment tools showed low power for TE PAV recall. AnchorWave precisely aligns up to three times more of the genome as position matches or indels than the closest competitive approach when comparing diverse genomes. Moreover, AnchorWave recalls transcription factor-binding sites at a rate of 1.05- to 74.85-fold higher than other tools with significantly lower false-positive alignments. AnchorWave complements available genome alignment tools by showing obvious improvement when applied to genomes with dispersed repeats, active TEs, high sequence diversity, and whole-genome duplication variation.

OUP accepted manuscript

Recognition of the important role of transposable elements (TEs) in eukaryotic genomes quickly led to a burgeoning literature modeling and estimating the effects of selection on TEs. Much of the empirical work on selection has focused on analyzing the site frequency spectrum (SFS) of TEs. But TE evolution differs from standard models in a number of ways that can impact the power and interpretation of the SFS. For example, rather than mutating under a clock-like model, transposition often occurs in bursts which can inflate particular frequency categories compared with expectations under a standard neutral model. If a TE burst has been recent, the excess of low-frequency polymorphisms can mimic the effect of purifying selection. Here, we investigate how transposition bursts affect the frequency distribution of TEs and the correlation between age and allele frequency. Using information on the TE age distribution, we propose an age-adjusted SFS to compare TEs and neutral polymorphisms to more effectively evaluate whether TEs are under selective constraints. We show that our approach can minimize instances of false inference of selective constraint, remains robust to simple demographic changes, and allows for a correct identification of even weak selection affecting TEs which experienced a transposition burst. The results presented here will help researchers working on TEs to more reliably identify the effects of selection on TEs without having to rely on the assumption of a constant transposition rate.

DARTS: An Algorithm for Domain-Associated Retrotransposon Search in Genome Assemblies

Retrotransposons comprise a substantial fraction of eukaryotic genomes, reaching the highest proportions in plants. Therefore, identification and annotation of retrotransposons is an important task in studying the regulation and evolution of plant genomes. The majority of computational tools for mining transposable elements (TEs) are designed for subsequent genome repeat masking, often leaving aside the element lineage classification and its protein domain composition. Additionally, studies focused on the diversity and evolution of a particular group of retrotransposons often require substantial customization efforts from researchers to adapt existing software to their needs. Here, we developed a computational pipeline to mine sequences of protein-coding retrotransposons based on the sequences of their conserved protein domains—DARTS (Domain-Associated Retrotransposon Search). Using the most abundant group of TEs in plants—long terminal repeat (LTR) retrotransposons (LTR-RTs)—we show that DARTS has radically higher sensitivity for LTR-RT identification compared to the widely accepted tool LTRharvest. DARTS can be easily customized for specific user needs. As a result, DARTS returns a set of structurally annotated nucleotide and amino acid sequences which can be readily used in subsequent comparative and phylogenetic analyses. DARTS may facilitate researchers interested in the discovery and detailed analysis of the diversity and evolution of retrotransposons, LTR-RTs, and other protein-coding TEs.

Whole-genome assembly and annotation of northern wild rice, Zizania palustris L., supports a whole-genome duplication in the Zizania genus

Zizania palustris L. (northern wild rice, NWR) is an aquatic grass native to North America that is notable for its nutritious grain. This is an important species with ecological, cultural and agricultural significance, specifically in the Great Lakes region of the USA. Using flow cytometry, we first estimated the NWR genome size to be 1.8 Gb. Using long- and short-range sequencing, Hi-C scaffolding and RNA-seq data from eight tissues, we generated an annotated whole-genome de novo assembly of NWR. The assembly was 1.29 Gb in length, highly repetitive (approx. 76.0%) and contained 46 421 putative protein-coding genes. The expansion of retrotransposons within the genome and a whole-genome duplication (WGD) after the Zizania-Oryza speciation event have both led to an increase in the genome size of NWR in comparison with Oryza sativa L. and Zizania latifolia. Both events depict a genome rapidly undergoing change over a short evolutionary time. Comparative analyses revealed the conservation of large syntenic blocks between NWR and O. sativa, which were used to identify putative seed-shattering genes. Estimates of divergence times revealed that the Zizania genus diverged from Oryza approximately 26-30 million years ago (26-30 MYA), whereas NWR and Z. latifolia diverged from one another approximately 6-8 MYA. Comparative genomics confirmed evidence of a WGD in the Zizania genus and provided support that the event occurred prior to the NWR-Z. latifolia speciation event. This genome assembly and annotation provides a valuable resource for comparative genomics in the Oryzeae tribe and provides an important resource for future conservation and breeding efforts of NWR.

Genome downsizing after polyploidy: mechanisms, rates and selection pressures

An analysis of over 10 000 plant genome sizes (GSs) indicates that most species have smaller genomes than expected given the incidence of polyploidy in their ancestries, suggesting selection for genome downsizing. However, comparing ancestral GS with the incidence of ancestral polyploidy suggests that the rate of DNA loss following polyploidy is likely to have been very low (4-70 Mb/million years, 4-482 bp/generation). This poses a problem. How might such small DNA losses be visible to selection, overcome the power of genetic drift and drive genome downsizing? Here we explore that problem, focussing on the role that double-strand break (DSB) repair pathways (non-homologous end joining and homologous recombination) may have played. We also explore two hypotheses that could explain how selection might favour genome downsizing following polyploidy: to reduce (i) nitrogen (N) and phosphate (P) costs associated with nucleic acid synthesis in the nucleus and the transcriptome and (ii) the impact of scaling effects of GS on cell size, which influences CO₂ uptake and water loss. We explore the hypothesis that losses of DNA must be fastest in early polyploid generations. Alternatively, if DNA loss is a more continuous process over evolutionary time, then we propose it is a byproduct of selection elsewhere, such as limiting the damaging activity of repetitive DNA. If so, then the impact of GS on photosynthesis, water use efficiency and/or nutrient costs at the nucleus level may be emergent properties, which have advantages, but not ones that could have been selected for over generational timescales.

Genome Size, Chromosome Number and Morphological Data Reveal Unexpected Infraspecific Variability in Festuca (Poaceae)

Polyploidy has played an important evolutionary role in the genus Festuca (Poaceae), and several ploidy levels (ranging from 2n = 2x = 14 to 2n = 12x = 84) have been detected to date. This study aimed to estimate the genome size and ploidy level of two subspecies belonging to the F. yvesii polyploid complex by flow cytometry and chromosome counting. The phenotypic variation of the cytotypes was also explored, based on herbarium material. The genome size of F. yvesii subsp. lagascae has been estimated for the first time. Nuclear 2C DNA content of F. yvesii subsp. summilusitana ranged from 21.44 to 31.91 pg, while that of F. yvesii subsp. lagascae was from 13.60 to 22.31 pg. We report the highest ploidy level detected for Festuca (2n = 14x = 98) and previously unknown cytotypes. A positive correlation between holoploid genome size and chromosome number counts shown herein was confirmed. The morphometric approach showed a slight trend towards an increase in the size of some organs consistent with the variation in the ploidy level. Differences in characters were usually significant only among the most extreme cytotypes of each subspecies, but, even in this case, the high overlapping ranges prevent their distinction.

Genome size variation and karyotype diversity in eight taxa of Sorbus sensu stricto (Rosaceae) from China

Eight taxa of Sorbus Linnaeus, 1753 sensu stricto (Rosaceae) from China have been studied karyologically through chromosome counting, chromosomal measurement and karyotype symmetry. Genome size was also estimated by flow cytometry. Six taxa, S. amabilis Cheng ex T.T.Yu et K.C.Kuan, 1963, S. hupehensis var. paucijuga (D.K. Zang et P.C. Huang, 1992) L.T. Lu, 2000, S. koehneana C.K. Schneider, 1906, S. pohuashanensis (Hance, 1875) Hedlund, 1901, S. scalaris Koehne, 1913 and S. wilsoniana C.K. Schneider, 1906 are diploids with 2n = 34, whereas two taxa, S. filipes Handel-Mazzetti,1933 and S. ovalis McAllister, 2005 are tetraploid with 2n = 68. In general, the chromosome size is mainly small, and karyotypes are symmetrical with predominance of metacentric chromosomes. Genome size variation of diploids and tetraploids is 1.401 pg -1.676 pg and 2.674 pg -2.684 pg, respectively. Chromosome numbers of S. amabilis and S. hupehensis var. paucijuga, and karyotype and genome size of eight taxa studied are reported for the first time. This study emphasised the reliability of flow cytometry in genome size determination to infer ploidy levels in Chinese native Sorbus species.