Diversity in genome size and GC content shows adaptive potential in orchids and is closely linked to partial endoreplication, plant life-history traits and climatic conditions

Rapid Targeted Assembly of the Proteome Reveals Evolutionary Variation of GC Content in Avian Lice

Nucleotide base composition plays an influential role in the molecular mechanisms involved in gene function, phenotype, and amino acid composition. GC content (proportion of guanine and cytosine in DNA sequences) shows a high level of variation within and among species. Many studies measure GC content in a small number of genes, which may not be representative of genome-wide GC variation. One challenge when assembling extensive genomic data sets for these studies is the significant amount of resources (monetary and computational) associated with data processing, and many bioinformatic tools have not been optimized for resource efficiency. Using a high-performance computing (HPC) cluster, we manipulated resources provided to the targeted gene assembly program, automated target restricted assembly method (aTRAM), to determine an optimum way to run the program to maximize resource use. Using our optimum assembly approach, we assembled and measured GC content of all of the protein-coding genes of a diverse group of parasitic feather lice. Of the 499 426 genes assembled across 57 species, feather lice were GC-poor (mean GC = 42.96%) with a significant amount of variation within and between species (GC range = 19.57%-73.33%). We found a significant correlation between GC content and standard deviation per taxon for overall GC and GC3, which could indicate selection for G and C nucleotides in some species. Phylogenetic signal of GC content was detected in both GC and GC3. This research provides a large-scale investigation of GC content in parasitic lice laying the foundation for understanding the basis of variation in base composition across species.

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

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

Phenotypic variation seems not to be associated with the genetic profile in Zygopetalum (Orchidaceae): a case study of a high-elevation rocky complex

BACKGROUND

Hybridization associated with polyploidy studies is rare in the tropics. The genus Zygopetalum (Orchidaceae) was investigated here as a case study of Neotropical plants. In the rocky highlands of the Ibitipoca State Park (ISP), southeast Brazil, individuals with intermediate colors and forms between the species Z. maculatum and Z. triste were commonly identified.

METHODS AND RESULTS

Chromosomal analysis and DNA quantity showed a uniform population. Regardless of the aspects related to the color and shape of floral structures, all individuals showed 2n = 96 chromosomes and an average of 14.05 pg of DNA. Irregularities in meiosis associated with chromosome number and C value suggest the occurrence of polyploidy. The genetic distance estimated using ISSR molecular markers revealed the existence of genetic variability not related to morphological clusters. Morphometric measurements of the flower pieces revealed that Z. maculatum shows higher variation than Z. triste although lacking a defined circumscription.

CONCLUSION

The observed variation can be explained by the polyploid and phenotypic plasticity resulting from the interaction of the genotypes with the heterogeneous environments observed in this habitat.

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.

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.

Genomic insights into the clonal reproductive Opuntia cochenillifera: mitochondrial and chloroplast genomes of the cochineal cactus for enhanced understanding of structural dynamics and evolutionary implications

Background

The cochineal cactus (Opuntia cochenillifera), notable for its substantial agricultural and industrial applications, predominantly undergoes clonal reproduction, which presents significant challenges in breeding and germplasm innovation. Recent developments in mitochondrial genome engineering offer promising avenues for introducing heritable mutations, potentially facilitating selective sexual reproduction through the creation of cytoplasmic male sterile genotypes. However, the lack of comprehensive mitochondrial genome information for Opuntia species hinders these efforts. Here, we intended to sequence and characterize its mitochondrial genome to maximize the potential of its genomes for evolutionary studies, molecular breeding, and molecular marker developments.

Results

We sequenced the total DNA of the O. cochenillifera using DNBSEQ and Nanopore platforms. The mitochondrial genome was then assembled using a hybrid assembly strategy using Unicycler software. We found that the mitochondrial genome of O. cochenillifera has a length of 1,156,235 bp, a GC content of 43.06%, and contains 54 unique protein-coding genes and 346 simple repeats. Comparative genomic analysis revealed 48 homologous fragments shared between mitochondrial and chloroplast genomes, with a total length of 47,935 bp. Additionally, the comparison of mitochondrial genomes from four Cactaceae species highlighted their dynamic nature and frequent mitogenomic reorganizations.

Conclusion

Our study provides a new perspective on the evolution of the organelle genome and its potential application in genetic breeding. These findings offer valuable insights into the mitochondrial genetics of Cactaceae, potentially facilitating future research and breeding programs aimed at enhancing the genetic diversity and adaptability of O. cochenillifera by leveraging its unique mitochondrial genome characteristics.

RNA barcode segments for SARS-CoV-2 identification from HCoVs and SARSr-CoV-2 lineages

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the pathogen responsible for coronavirus disease 2019 (COVID-19), continues to evolve, giving rise to more variants and global reinfections. Previous research has demonstrated that barcode segments can effectively and cost-efficiently identify specific species within closely related populations. In this study, we designed and tested RNA barcode segments based on genetic evolutionary relationships to facilitate the efficient and accurate identification of SARS-CoV-2 from extensive virus samples, including human coronaviruses (HCoVs) and SARSr-CoV-2 lineages. Nucleotide sequences sourced from NCBI and GISAID were meticulously selected and curated to construct training sets, encompassing 1733 complete genome sequences of HCoVs and SARSr-CoV-2 lineages. Through genetic-level species testing, we validated the accuracy and reliability of the barcode segments for identifying SARS-CoV-2. Subsequently, 75 main and subordinate species-specific barcode segments for SARS-CoV-2, located in ORF1ab, S, E, ORF7a, and N coding sequences, were intercepted and screened based on single-nucleotide polymorphism sites and weighted scores. Post-testing, these segments exhibited high recall rates (nearly 100%), specificity (almost 30% at the nucleotide level), and precision (100%) performance on identification. They were eventually visualized using one and two-dimensional combined barcodes and deposited in an online database (http://virusbarcodedatabase.top/). The successful integration of barcoding technology in SARS-CoV-2 identification provides valuable insights for future studies involving complete genome sequence polymorphism analysis. Moreover, this cost-effective and efficient identification approach also provides valuable reference for future research endeavors related to virus surveillance.

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.

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.

Repeat-based phylogenomics shed light on unclear relationships in the monocentric genus Juncus L. (Juncaceae)

The repetitive fraction (repeatome) of eukaryotic genomes is diverse and usually fast evolving, being an important tool for clarify plant systematics. The genus Juncus L. comprises 332 species, karyotypically recognized by having holocentric chromosomes. However, four species were recently described as monocentric, yet our understanding of their genome evolution is largely masked by unclear phylogenetic relationships. Here, we reassess the current Juncus systematics using low-coverage genome skimming data of 33 taxa to construct repeats, nuclear rDNA and plastome-based phylogenetic hypothesis. Furthermore, we characterize the repeatome and chromosomal distribution of Juncus-specific centromeric repeats/CENH3 protein to test the monocentricity reach in the genus. Repeat-base phylogenies revealed topologies congruent with the rDNA tree, but not with the plastome tree. The incongruence between nuclear and plastome chloroplast dataset suggest an ancient hybridization in the divergence of Juncotypus and Tenageia sections 40 Myr ago. The phylogenetic resolution at section level was better fitted with the rDNA/repeat-based approaches, with the recognition of two monophyletic sections (Stygiopsis and Tenageia). We found specific repeatome trends for the main lineages, such as the higher abundances of TEs in the Caespitosi and Iridifolii + Ozophyllum clades. CENH3 immunostaining confirmed the monocentricity of Juncus, which can be a generic synapomorphy for the genus. The heterogeneity of the repeatomes, with high phylogenetic informativeness, identified here may be correlated with their ancient origin (56 Mya) and reveals the potential of comparative genomic analyses for understanding plant systematics and evolution.

Complete mitogenome assembly of Selenicereus monacanthus revealed its molecular features, genome evolution, and phylogenetic implications

BACKGROUND

Mitochondria are the powerhouse of the cell and are critical for plant growth and development. Pitaya (Selenicereus or Hylocereus) is the most important economic crop in the family Cactaceae and is grown worldwide, however its mitogenome is unreported.

RESULTS

This study assembled the complete mitogenome of the red skin and flesh of pitaya (Selenicereus monacanthus). It is a full-length, 2,290,019 bp circular molecule encoding 59 unique genes that only occupy 2.17% of the entire length. In addition, 4,459 pairs of dispersed repeats (≥ 50 bp) were identified, accounting for 84.78% of the total length, and three repeats (394,588, 124,827, and 13,437 bp) mediating genomic recombination were identified by long read mapping and Sanger sequencing. RNA editing events were identified in all 32 protein-coding genes (PCGs), among which four sites (nad1-2, nad4L-2, atp9-copy3-223, and ccmFC-1309) were associated with the initiation or termination of PCGs. Seventy-eight homologous fragments of the chloroplast genome were identified in the mitogenome, the longest having 4,523 bp. In addition, evolutionary analyses suggest that S. monacanthus may have undergone multiple genomic reorganization events during evolution, with the loss of at least nine PCGs (rpl2, rpl10, rps2, rps3, rps10, rps11, rps14, rps19, and sdh3).

CONCLUSIONS

This study revealed the genetic basis of the S. monacanthus mitogenome, and provided a scientific basis for further research on phenotypic traits and germplasm resource development.

First Vanilla planifolia High-Density Genetic Linkage Map Provides Quantitative Trait Loci for Resistance to Fusarium oxysporum

Fusarium oxysporum f. sp. radicis-vanillae (Forv), the causal agent of root and stem rot disease, is the main pathogen affecting vanilla production. Sources of resistance have been reported in Vanilla planifolia G. Jackson ex Andrews, the main cultivated vanilla species. In this study, we developed the first high-density genetic map in this species with 1,804 genotyping-by-sequencing (GBS)-generated single nucleotide polymorphism (SNP) markers using 125 selfed progenies of the CR0040 traditional vanilla cultivar. Sixteen linkage groups (LG) were successfully constructed, with a mean of 113 SNPs and an average length of 207 cM per LG. The map had a high density with an average of 5.45 SNP every 10 cM and an average distance of 1.85 cM between adjacent markers. The first three LG were aligned against the first assembled chromosome of CR0040, and the other 13 LG were correctly associated with the other 13 assembled chromosomes. The population was challenged with the highly pathogenic Forv strain Fo072 using the root-dip inoculation method. Five traits were mapped, and 20 QTLs were associated with resistance to Fo072. Among the genes retrieved in the CR0040 physical regions associated with QTLs, genes potentially involved in biotic resistance mechanisms, coding for kinases, E3 ubiquitin ligases, pentatricopeptide repeat-containing proteins, and one leucine-rich repeat receptor underlying the qFo72_08.1 QTL have been highlighted. This study should provide useful resources for marker-assisted selection in V. planifolia.

Endoreplication-Why Are We Not Using Its Full Application Potential?

Endoreplication-a process that is common in plants and also accompanies changes in the development of animal organisms-has been seen from a new perspective in recent years. In the paper, we not only shed light on this view, but we would also like to promote an understanding of the application potential of this phenomenon in plant cultivation. Endoreplication is a pathway for cell development, slightly different from the classical somatic cell cycle, which ends with mitosis. Since many rounds of DNA synthesis take place within its course, endoreplication is a kind of evolutionary compensation for the relatively small amount of genetic material that plants possess. It allows for its multiplication and active use through transcription and translation. The presence of endoreplication in plants has many positive consequences. In this case, repeatedly produced copies of genes, through the corresponding transcripts, help the plant acquire the favorable properties for which proteins are responsible directly or indirectly. These include features that are desirable in terms of cultivation and marketing: a greater saturation of fruit and flower colors, a stronger aroma, a sweeter fruit taste, an accumulation of nutrients, an increased resistance to biotic and abiotic stress, superior tolerance to adverse environmental conditions, and faster organ growth (and consequently the faster growth of the whole plant and its biomass). The two last features are related to the nuclear-cytoplasmic ratio-the greater the content of DNA in the nucleus, the higher the volume of cytoplasm, and thus the larger the cell size. Endoreplication not only allows cells to reach larger sizes but also to save the materials used to build organelles, which are then passed on to daughter cells after division, thus ending the classic cell cycle. However, the content of genetic material in the cell nucleus determines the number of corresponding organelles. The article also draws attention to the potential practical applications of the phenomenon and the factors currently limiting its use.

Speciation across the Earth driven by global cooling in terrestrial orchids

Although climate change has been implicated as a major catalyst of diversification, its effects are thought to be inconsistent and much less pervasive than localized climate or the accumulation of species with time. Focused analyses of highly speciose clades are needed in order to disentangle the consequences of climate change, geography, and time. Here, we show that global cooling shapes the biodiversity of terrestrial orchids. Using a phylogeny of 1,475 species of Orchidoideae, the largest terrestrial orchid subfamily, we find that speciation rate is dependent on historic global cooling, not time, tropical distributions, elevation, variation in chromosome number, or other types of historic climate change. Relative to the gradual accumulation of species with time, models specifying speciation driven by historic global cooling are over 700 times more likely. Evidence ratios estimated for 212 other plant and animal groups reveal that terrestrial orchids represent one of the best-supported cases of temperature-spurred speciation yet reported. Employing >2.5 million georeferenced records, we find that global cooling drove contemporaneous diversification in each of the seven major orchid bioregions of the Earth. With current emphasis on understanding and predicting the immediate impacts of global warming, our study provides a clear case study of the long-term impacts of global climate change on biodiversity.

Comparative analyses and phylogenetic relationships of thirteen Pholidota species (Orchidaceae) inferred from complete chloroplast genomes

BACKGROUND

The orchid genus Pholidota Lindl. ex Hook. is economically important as some species has long been used in traditional medicine. However, the systematic status of the genus and intergeneric relationships inferred from previous molecular studies are unclear due to insufficient sampling and lack of informative sites. So far, only limited genomic information has been available. The taxonomy of Pholidota remains unresolved and somewhat controversial. In this study, the complete chloroplast (cp.) genomes of thirteen Pholidota species were sequenced and analyzed to gain insight into the phylogeny of Pholidota and mutation patterns in their cp. genomes.

RESULTS

All examined thirteen Pholidota cp. genomes exhibited typical quadripartite circular structures, with the size ranging from 158,786 to 159,781 bp. The annotation contained a total of 135 genes in each cp. genome, i.e., 89 protein-coding genes, 38 tRNA genes, and eight rRNA genes. The codon usage analysis indicated the preference of A/U-ending codons. Repeat sequence analysis identified 444 tandem repeats, 322 palindromic repeats and 189 dispersed repeats. A total of 525 SSRs, 13,834 SNPs and 8,630 InDels were detected. Six mutational hotspots were identified as potential molecular markers. These molecular markers and highly variable regions are expected to facilitate future genetic and genomic studies. Our phylogenetic analyses confirmed the polyphyletic status of the genus Pholidota, with species grouped into four main clades: Pholidota s.s. was resolved as the sister to a clade containing species of Coelogyne; the other two clades clustered together with species of Bulleyia and Panisea, respectively; species P. ventricosa was placed at the basal position, deviated from all other species.

CONCLUSION

This is the first study to comprehensively examine the genetic variations and systematically analyze the phylogeny and evolution of Pholidota based on plastid genomic data. These findings contribute to a better understanding of plastid genome evolution of Pholidota and provide new insights into the phylogeny of Pholidota and its closely related genera within the subtribe Coelogyninae. Our research has laid the foundation for future studies on the evolutionary mechanisms and classification of this economically and medicinally important genus.

Genome size variation and polyploidy prevalence in the genus Eragrostis are associated with the global dispersal in arid area

BACKGROUND

Biologists have long debated the drivers of the genome size evolution and variation ever since Darwin. Assumptions for the adaptive or maladaptive consequences of the associations between genome sizes and environmental factors have been proposed, but the significance of these hypotheses remains controversial. Eragrostis is a large genus in the grass family and is often used as crop or forage during the dry seasons. The wide range and complex ploidy levels make Eragrostis an excellent model for investigating how the genome size variation and evolution is associated with environmental factors and how these changes can ben interpreted.

METHODS

We reconstructed the Eragrostis phylogeny and estimated genome sizes through flow cytometric analyses. Phylogenetic comparative analyses were performed to explore how genome size variation and evolution is related to their climatic niches and geographical ranges. The genome size evolution and environmental factors were examined using different models to study the phylogenetic signal, mode and tempo throughout evolutionary history.

RESULTS

Our results support the monophyly of Eragrostis. The genome sizes in Eragrostis ranged from ~0.66 pg to ~3.80 pg. We found that a moderate phylogenetic conservatism existed in terms of the genome sizes but was absent from environmental factors. In addition, phylogeny-based associations revealed close correlations between genome sizes and precipitation-related variables, indicating that the genome size variation mainly caused by polyploidization may have evolved as an adaptation to various environments in the genus Eragrostis.

CONCLUSION

This is the first study to take a global perspective on the genome size variation and evolution in the genus Eragrostis. Our results suggest that the adaptation and conservatism are manifested in the genome size variation, allowing the arid species of Eragrostis to spread the xeric area throughout the world.

Genomic C-Value Variation Analysis in Jujube (Ziziphus jujuba Mill.) in the Middle Yellow River Basin

Chinese jujube (Ziziphus jujuba Mill.) originated in the Yellow River basin (YRB) of the Shanxi-Shaanxi region. The genomic C-value is a crucial indicator for plant breeding and germplasm evaluation. In this study, we used flow cytometry to determine the genomic C-values of jujube germplasms in the YRB of the Shanxi-Shaanxi region and evaluated their differences in different sub-regions. Of the 29 sub-regions, the highest and lowest variations were in Linxian and Xiaxian, respectively. The difference between jujube germplasms was highly significant (F = 14.89, p < 0.0001) in Linxian. Cluster analysis showed that both cluster 2 and 4 belonged to Linxian, which were clearly separated from other taxa but were cross-distributed in them. Linxian County is an important gene exchange center in the YRB of the Shanxi-Shaanxi region. Principal component analysis showed that cluster 1 had low genomic C-values and single-fruit weights and cluster 2 had high genomic C-values and vitamin C contents. The genomic C-value was correlated with single-fruit weight and vitamin C content. In addition, the genomic C-value was used to predict fruit agronomic traits, providing a reference for shortening the breeding cycle and genetic diversity-related studies of jujube germplasm.

Advances in Vertebrate (Cyto)Genomics Shed New Light on Fish Compositional Genome Evolution

Cytogenetic and compositional studies considered fish genomes rather poor in guanine-cytosine content (GC%) because of a putative "sharp increase in genic GC% during the evolution of higher vertebrates". However, the available genomic data have not been exploited to confirm this viewpoint. In contrast, further misunderstandings in GC%, mostly of fish genomes, originated from a misapprehension of the current flood of data. Utilizing public databases, we calculated the GC% in animal genomes of three different, technically well-established fractions: DNA (entire genome), cDNA (complementary DNA), and cds (exons). Our results across chordates help set borders of GC% values that are still incorrect in literature and show: (i) fish in their immense diversity possess comparably GC-rich (or even GC-richer) genomes as higher vertebrates, and fish exons are GC-enriched among vertebrates; (ii) animal genomes generally show a GC-enrichment from the DNA, over cDNA, to the cds level (i.e., not only the higher vertebrates); (iii) fish and invertebrates show a broad(er) inter-quartile range in GC%, while avian and mammalian genomes are more constrained in their GC%. These results indicate no sharp increase in the GC% of genes during the transition to higher vertebrates, as stated and numerously repeated before. We present our results in 2D and 3D space to explore the compositional genome landscape and prepared an online platform to explore the AT/GC compositional genome evolution.

The ecological adaptation of the unparalleled plastome character evolution in slipper orchids

Plastomes may have undergone adaptive evolution in the process of plant adaptation to diverse environments, whereby species may differ in plastome characters. Cypripedioideae successfully colonized distinct environments and could be an ideal group for studying the interspecific variation and adaptive evolution of plastomes. Comparative study of plastomes, ancestral state reconstruction, phylogenetic-based analysis, ecological niche modelling, and selective pressure analysis were conducted to reveal the evolutionary patterns of plastomes in Cypripedioideae and their relationship with environmental factors. The plastomes of the three evolved genera had reduced plastome size, increased GC content, and compacted gene content compared to the basal group. Variations in plastome size and GC content are proved to have clear relationships with climate regions. Furthermore, ecological niche modelling revealed that temperature and water factors are important climatic factors contributing to the distributional difference which is directly correlated with the climate regions. The temperature-sensitive genes ndh genes, infA, and rpl20 were found to be either lost/pseudogenized or under positive selection in the evolved groups. Unparalleled plastome character variations were discovered in slipper orchids. Our study indicates that variations in plastome characters have adaptive consequences and that temperature and water factors are important climatic factors that affect plastome evolution. This research highlights the expectation that plants can facilitate adaptation to different environmental conditions with the changes in plastome and has added critical insight for understanding the process of plastome evolution in plants.

The ecology of palm genomes: repeat-associated genome size expansion is constrained by aridity

Genome size varies 2400-fold across plants, influencing their evolution through changes in cell size and cell division rates which impact plants' environmental stress tolerance. Repetitive element expansion explains much genome size diversity, and the processes structuring repeat 'communities' are analogous to those structuring ecological communities. However, which environmental stressors influence repeat community dynamics has not yet been examined from an ecological perspective. We measured genome size and leveraged climatic data for 91% of genera within the ecologically diverse palm family (Arecaceae). We then generated genomic repeat profiles for 141 palm species, and analysed repeats using phylogenetically informed linear models to explore relationships between repeat dynamics and environmental factors. We show that palm genome size and repeat 'community' composition are best explained by aridity. Specifically, Ty3-gypsy and TIR elements were more abundant in palm species from wetter environments, which generally had larger genomes, suggesting amplification. By contrast, Ty1-copia and LINE elements were more abundant in drier environments. Our results suggest that water stress inhibits repeat expansion through selection on upper genome size limits. However, elements that may associate with stress-response genes (e.g. Ty1-copia) have amplified in arid-adapted palm species. Overall, we provide novel evidence of climate influencing the assembly of repeat 'communities'.

The relationship between transposable elements and ecological niches in the Greater Cape Floristic Region: A study on the genus Pteronia (Asteraceae)

Non-coding repetitive DNA (repeatome) is an active part of the nuclear genome, involved in its structure, evolution and function. It is dominated by transposable elements (TEs) and satellite DNA and is prone to the most rapid changes over time. The TEs activity presumably causes the global genome reorganization and may play an adaptive or regulatory role in response to environmental challenges. This assumption is applied here for the first time to plants from the Cape Floristic hotspot to determine whether changes in repetitive DNA are related to responses to a harsh, but extremely species-rich environment. The genus Pteronia (Asteraceae) serves as a suitable model group because it shows considerable variation in genome size at the diploid level and has high and nearly equal levels of endemism in the two main Cape biomes, Fynbos and Succulent Karoo. First, we constructed a phylogeny based on multiple low-copy genes that served as a phylogenetic framework for detecting quantitative and qualitative changes in the repeatome. Second, we performed a comparative analysis of the environments of two groups of Pteronia differing in their TEs bursts. Our results suggest that the environmental transition from the Succulent Karoo to the Fynbos is accompanied by TEs burst, which is likely also driving phylogenetic divergence. We thus hypothesize that analysis of rapidly evolving repeatome could serve as an important proxy for determining the molecular basis of lineage divergence in rapidly radiating groups.

In-depth analysis of genomes and functional genomics of orchid using cutting-edge high-throughput sequencing

High-throughput sequencing technology has been facilitated the development of new methodologies and approaches for studying the origin and evolution of plant genomes and subgenomes, population domestication, and functional genomics. Orchids have tens of thousands of members in nature. Many of them have promising application potential in the extension and conservation of the ecological chain, the horticultural use of ornamental blossoms, and the utilization of botanical medicines. However, a large-scale gene knockout mutant library and a sophisticated genetic transformation system are still lacking in the improvement of orchid germplasm resources. New gene editing tools, such as the favored CRISPR-Cas9 or some base editors, have not yet been widely applied in orchids. In addition to a large variety of orchid cultivars, the high-precision, high-throughput genome sequencing technology is also required for the mining of trait-related functional genes. Nowadays, the focus of orchid genomics research has been directed to the origin and classification of species, genome evolution and deletion, gene duplication and chromosomal polyploidy, and flower morphogenesis-related regulation. Here, the progressing achieved in orchid molecular biology and genomics over the past few decades have been discussed, including the evolution of genome size and polyploidization. The frequent incorporation of LTR retrotransposons play important role in the expansion and structural variation of the orchid genome. The large-scale gene duplication event of the nuclear genome generated plenty of recently tandem duplicated genes, which drove the evolution and functional divergency of new genes. The evolution and loss of the plastid genome, which mostly affected genes related to photosynthesis and autotrophy, demonstrated that orchids have experienced more separate transitions to heterotrophy than any other terrestrial plant. Moreover, large-scale resequencing provide useful SNP markers for constructing genetic maps, which will facilitate the breeding of novel orchid varieties. The significance of high-throughput sequencing and gene editing technologies in the identification and molecular breeding of the trait-related genes in orchids provides us with a representative trait-improving gene as well as some mechanisms worthy of further investigation. In addition, gene editing has promise for the improvement of orchid genetic transformation and the investigation of gene function. This knowledge may provide a scientific reference and theoretical basis for orchid genome studies.

A chromosome-level, haplotype-phased Vanilla planifolia genome highlights the challenge of partial endoreplication for accurate whole-genome assembly

Vanilla planifolia, the species cultivated to produce one of the world's most popular flavors, is highly prone to partial genome endoreplication, which leads to highly unbalanced DNA content in cells. We report here the first molecular evidence of partial endoreplication at the chromosome scale by the assembly and annotation of an accurate haplotype-phased genome of V. planifolia. Cytogenetic data demonstrated that the diploid genome size is 4.09 Gb, with 16 chromosome pairs, although aneuploid cells are frequently observed. Using PacBio HiFi and optical mapping, we assembled and phased a diploid genome of 3.4 Gb with a scaffold N50 of 1.2 Mb and 59 128 predicted protein-coding genes. The atypical k-mer frequencies and the uneven sequencing depth observed agreed with our expectation of unbalanced genome representation. Sixty-seven percent of the genes were scattered over only 30% of the genome, putatively linking gene-rich regions and the endoreplication phenomenon. By contrast, low-coverage regions (non-endoreplicated) were rich in repeated elements but also contained 33% of the annotated genes. Furthermore, this assembly showed distinct haplotype-specific sequencing depth variation patterns, suggesting complex molecular regulation of endoreplication along the chromosomes. This high-quality, anchored assembly represents 83% of the estimated V. planifolia genome. It provides a significant step toward the elucidation of this complex genome. To support post-genomics efforts, we developed the Vanilla Genome Hub, a user-friendly integrated web portal that enables centralized access to high-throughput genomic and other omics data and interoperable use of bioinformatics tools.

Stoichiometry versus ecology: the relationships between genome size and guanine-cytosine content, and tissue nitrogen and phosphorus in grassland herbs

BACKGROUND AND AIMS

Plant tissue nitrogen (N) and phosphorus (P) and genome traits, such as genome size and guanine-cytosine (GC) content, scale with growth or metabolic rates and are linked to plant ecological strategy spectra. Tissue NP stoichiometry and genome traits are reported to affect plant growth, metabolic rates or ecological strategies in contrasting ways, although the elemental costs for building and maintaining DNA are typically overlooked.

METHODS

We formulated stoichiometry- and ecology-based predictions on the relationship between genome size and GC content to tissue N, P and N : P and tested them on a set of 130 herbaceous species from a temperate grassland using ordinary, phylogenetic and quantile regression.

KEY RESULTS

Genome size was only negatively linked to plant N and N : P in species with very small genomes. We found no link between genome size and plant P. GC content was negatively linked to plant N and P but we found these significant links consistently in both GC-rich and GC-poor species. Finally, GC content correlated positively with plant N : P but only in species with GC-rich genomes.

CONCLUSIONS

Our results provide stronger support for the ecology-based predictions than the stoichiometry-based predictions, and for the links between GC content and plant N and P stoichiometry than for genome size. We argue that the theories of plant metabolic rates and ecological strategies (resource-acquisitive vs. conservative or ruderal vs. stress-tolerator spectra) better explain interspecific genome-NP stoichiometry relationships at the tissue level (although relatively weakly) than the stoichiometric theory based on the elemental costs for building and maintaining DNA.

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.

Genome Size Variation in Dianthus sylvestris Wulfen sensu lato (Caryophyllaceae)

Genome size (GS) is an important characteristic that may be helpful in delimitation of taxa, and multiple studies have shown correlations between intraspecific GS variation and morphological or environmental factors, as well as its geographical segregation. We estimated a relative GS (RGS) of 707 individuals from 162 populations of Dianthus sylvestris with a geographic focus on the Balkan Peninsula, but also including several populations from the European Alps. Dianthus sylvestris is morphologically variable species thriving in various habitats and six subspecies have been recognized from the Balkan Peninsula. Our RGS data backed-up with chromosome counts revealed that the majority of populations were diploid (2n = 30), but ten tetraploid populations have been recorded in D. sylvestris subsp. sylvestris from Istria (Croatia, Italy). Their monoploid RGS is significantly lower than that of the diploids, indicating genome downsizing. In addition, the tetraploids significantly differ from their diploid counterparts in an array of morphological and environmental characteristics. Within the diploid populations, the RGS is geographically and only partly taxonomically correlated, with the highest RGS inferred in the southern Balkan Peninsula and the Alps. We demonstrate greater RGS variation among the Balkan populations compared to the Alps, which is likely a result of more pronounced evolutionary differentiation within the Balkan Peninsula. In addition, a deep RGS divergence within the Alps likely points to persistence of the alpine populations in different Pleistocene refugia.

The specific DNA barcodes based on chloroplast genes for species identification of Theaceae plants

More than 600 species in over 40 genera have been identified in family Theaceae worldwide. The accurate identification of Theaceae plants can ensure the market economic order, and it plays a vital role in achieving the sustainable utilization of germplasm resources. DNA barcoding, one of the most potential species identification technologies at present, has advanced in the rapid, accurate and repetitive discrimination of species. In this study, matK + ndhF + ycf1 was observed as the optimal combined candidate gene sequence of DNA barcodes by analyzing genetic information of four single chloroplast DNA sequences, including matK, rbcL, ndhF and ycf1, as well as six combined gene sequences. Subsequently, the experiments were performed on phylogenetic analysis based on genetic distance to study the phylogenetic relationship of Theaceae plants and evaluate the species identification accuracy of matK + ndhF + ycf1. Lastly, the species-specific DNA barcodes were designed by searching the variable sites (one type of single nucleotide polymorphism sites) for the accurate identification of Camellia amplexicaulis, Franklinia alatamaha, Gordonia brandegeei and Stewartia micrantha. The previous methods of screening and testing candidate gene sequences were optimized, and innovation was made in the above methods. The process of making visual DNA barcodes was standardized. Besides, DNA barcoding technology increased the accuracy of species identification and DNA barcoding was analyzed in accordance with the theories of population genetics (e.g., neutral theory of molecular evolution). The results of the study will lay a basis for the identification and protection of Theaceae species and germplasm resources.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-022-01175-7.

Strategies of tolerance reflected in two North American maple genomes

The first chromosome‐scale assemblies for North American members of the Acer genus, sugar maple (Acer saccharum) and boxelder (Acer negundo), as well as transcriptomic evaluation of the abiotic stress response in A. saccharum are reported. This integrated study describes in‐depth aspects contributing to each species' approach to tolerance and applies current knowledge in many areas of plant genome biology with Acer physiology to help convey the genomic complexities underlying tolerance in broadleaf tree species.

Frequent ploidy changes in Salicaceae indicates widespread sharing of the salicoid whole genome duplication by the relatives of Populus L. and Salix L

BACKGROUNDS

Populus and Salix belong to Salicaceae and are used as models to investigate woody plant physiology. The variation of karyotype and nuclear DNA content can partly reflect the evolutionary history of the whole genome, and can provide critical information for understanding, predicting, and potentially ameliorating the woody plant traits. Therefore, it is essential to study the chromosome number (CN) and genome size in detail to provide information for revealing the evolutionary process of Salicaceae.

RESULTS

In this study, we report the somatic CNs of seventeen species from eight genera in Salicaceae. Of these, CNs for twelve species and for five genera are reported for the first time. Among the three subfamilies of Salicaceae, the available data indicate CN in Samydoideae is n = 21, 22, 42. The only two genera, Dianyuea and Scyphostegia, in Scyphostegioideae respectively have n = 9 and 18. In Salicoideae, Populus, Salix and five genera closely related to them (Bennettiodendron, Idesia, Carrierea, Poliothyrsis, Itoa) are based on relatively high CNs from n = 19, 20, 21, 22 to n = 95 in Salix. However, the other genera of Salicoideae are mainly based on relatively low CNs of n = 9, 10, 11. The genome sizes of 35 taxa belonging to 14 genera of Salicaceae were estimated. Of these, the genome sizes of 12 genera and all taxa except Populus euphratica are first reported. Except for Dianyuea, Idesia and Bennettiodendron, all examined species have relatively small genome sizes of less than 1 pg, although polyploidization exists.

CONCLUSIONS

The variation of CN and genome size across Salicaceae indicates frequent ploidy changes and a widespread sharing of the salicoid whole genome duplication (WGD) by the relatives of Populus and Salix. The shrinkage of genome size after WGD indicates massive loss of genomic components. The phylogenetic asymmetry in clade of Populus, Salix, and their close relatives suggests that there is a lag-time for the subsequent radiations after the salicoid WGD event. Our results provide useful data for studying the evolutionary events of Salicaceae.

Application-based guidelines for best practices in plant flow cytometry

Flow cytometry (FCM) is currently the most widely-used method to establish nuclear DNA content in plants. Since simple, 1-3-parameter, flow cytometers, which are sufficient for most plant applications, are commercially available at a reasonable price, the number of laboratories equipped with these instruments, and consequently new FCM users, has greatly increased over the last decade. This paper meets an urgent need for comprehensive recommendations for best practices in FCM for different plant science applications. We discuss advantages and limitations of establishing plant ploidy, genome size, DNA base composition, cell cycle activity, and level of endoreduplication. Applications of such measurements in plant systematics, ecology, molecular biology research, reproduction biology, tissue cultures, plant breeding, and seed sciences are described. Advice is included on how to obtain accurate and reliable results, as well as how to manage troubleshooting that may occur during sample preparation, cytometric measurements, and data handling. Each section is followed by best practice recommendations; tips as to what specific information should be provided in FCM papers are also provided.

The Welwitschia genome reveals a unique biology underpinning extreme longevity in deserts

The gymnosperm Welwitschia mirabilis belongs to the ancient, enigmatic gnetophyte lineage. It is a unique desert plant with extreme longevity and two ever-elongating leaves. We present a chromosome-level assembly of its genome (6.8 Gb/1 C) together with methylome and transcriptome data to explore its astonishing biology. We also present a refined, high-quality assembly of Gnetum montanum to enhance our understanding of gnetophyte genome evolution. The Welwitschia genome has been shaped by a lineage-specific ancient, whole genome duplication (~86 million years ago) and more recently (1-2 million years) by bursts of retrotransposon activity. High levels of cytosine methylation (particularly at CHH motifs) are associated with retrotransposons, whilst long-term deamination has resulted in an exceptionally GC-poor genome. Changes in copy number and/or expression of gene families and transcription factors (e.g. R2R3MYB, SAUR) controlling cell growth, differentiation and metabolism underpin the plant's longevity and tolerance to temperature, nutrient and water stress.

Repeat proliferation and partial endoreplication jointly shape the patterns of genome size evolution in orchids

Although the evolutionary drivers of genome size change are known, the general patterns and mechanisms of plant genome size evolution are yet to be established. Here we aim to assess the relative importance of proliferation of repetitive DNA, chromosomal variation (including polyploidy), and the type of endoreplication for genome size evolution of the Pleurothallidinae, the most species-rich orchid lineage. Phylogenetic relationships between 341 Pleurothallidinae representatives were refined using a target enrichment hybrid capture combined with high-throughput sequencing approach. Genome size and the type of endoreplication were assessed using flow cytometry supplemented with karyological analysis and low-coverage Illumina sequencing for repeatome analysis on a subset of samples. Data were analyzed using phylogeny-based models. Genome size diversity (0.2-5.1 Gbp) was mostly independent of profound chromosome count variation (2n = 12-90) but tightly linked with the overall content of repetitive DNA elements. Species with partial endoreplication (PE) had significantly greater genome sizes, and genomic repeat content was tightly correlated with the size of the non-endoreplicated part of the genome. In PE species, repetitive DNA is preferentially accumulated in the non-endoreplicated parts of their genomes. Our results demonstrate that proliferation of repetitive DNA elements and PE together shape the patterns of genome size diversity in orchids.

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.

Exploring adaptive landscapes across deep time: A case study using echinoid body size

Adaptive landscapes are a common way of conceptualizing the phenotypic evolution of lineages across deep time. Although multiple approaches exist to implement this concept into operational models of trait evolution, inferring adaptive landscapes from comparative datasets remains challenging. Here, I explore the macroevolutionary dynamics of echinoid body size using data from over 5000 specimens and a phylogenetic framework incorporating a dense fossil sampling and spanning approximately 270 million years. Furthermore, I implement a novel approach of exploring alternative parameterizations of adaptive landscapes that succeeds in finding simpler, yet better-fitting models. Echinoid body size has been constrained to evolve within a single adaptive optimum for much of the clade's history. However, most of the morphological disparity of echinoids was generated by multiple regime shifts that drove the repeated evolution of miniaturized and gigantic forms. Events of body size innovation occurred predominantly in the Late Cretaceous and were followed by a drastic slowdown following the Cretaceous-Paleogene mass extinction. The discovery of these patterns is contingent upon directly sampling fossil taxa. The macroevolution of echinoid body size is therefore characterized by a late increase in disparity (likely linked to an expansion of ecospace), generated by active processes driving lineages toward extreme morphologies.

Rapid radiation of angraecoids (Orchidaceae, Angraecinae) in tropical Africa characterised by multiple karyotypic shifts under major environmental instability

Angraecoid orchids present a remarkable diversity of chromosome numbers, which makes them a highly suitable system for exploring the impact of karyotypic changes on cladogenesis, diversification and morphological differentiation. We compiled an annotated cytotaxonomic checklist for 126 species of Angraecinae, which was utilised to reconstruct chromosomal evolution using a newly-produced, near-comprehensive phylogenetic tree that includes 245 angraecoid taxa. In tandem with this improved phylogenetic framework, using combined Bayesian, maximum likelihood and parsimony approaches on ITS-1 and five plastid markers, we propose a new cladistic nomenclature for the angraecoids, and we estimate a new timeframe for angraecoid radiation based on a secondary calibration, and calculate diversification rates using a Bayesian approach. Coincident divergence dates between clades with identical geographical distributions in the angraecoids and the pantropical orchid genus Bulbophyllum suggest that the same events may have intervened in the dispersal of these two epiphytic groups between Asia, continental Africa, Madagascar and the Neotropics. The major angraecoid lineages probably began to differentiate in the Middle Miocene, and most genera and species emerged respectively around the Late Miocene-Pliocene boundary and the Pleistocene. Ancestral state reconstruction using maximum likelihood estimation revealed an eventful karyotypic history dominated by descending dysploidy. Karyotypic shifts seem to have paralleled cladogenesis in continental tropical Africa, where approximately 90% of the species have descended from at least one inferred transition from n = 17-18 to n = 25 during the Middle Miocene Climatic Transition, followed by some clade-specific descending and ascending dysploidy from the Late Miocene to the Pleistocene. Conversely, detected polyploidy is restricted to a few species lineages mostly originating during the Pleistocene. No increases in net diversification could be related to chromosome number changes, and the apparent net diversification was found to be highest in Madagascar, where karyotypic stasis predominates. Finally, shifts in chromosome number appear to have paralleled the evolution of rostellum structure, leaflessness, and conspicuous changes in floral colour.

Correlation Analysis Reveals an Important Role of GC Content in Accumulation of Deletion Mutations in the Coding Region of Angiosperm Plastomes

Variation in GC content is assumed to correlate with various processes, including mutation biases, recombination, and environmental parameters. To date, most genomic studies exploring the evolution of GC content have focused on nuclear genomes, but relatively few have concentrated on organelle genomes. We explored the mechanisms maintaining the GC content in angiosperm plastomes, with a particular focus on the hypothesis of phylogenetic dependence and the correlation with deletion mutations. We measured three genetic traits, namely, GC content, A/T tracts, and G/C tracts, in the coding region of plastid genomes for 1382 angiosperm species representing 350 families and 64 orders, and tested the phylogenetic signal. Then, we performed correlation analyses and revealed the variation in evolutionary rate of selected traits using RRphylo. The plastid GC content in the coding region varied from 28.10% to 43.20% across angiosperms, with a few non-photosynthetic species showing highly reduced values, highlighting the significance of functional constraints. We found strong phylogenetic signal in A/T tracts, but weak ones in GC content and G/C tracts, indicating adaptive potential. GC content was positively and negatively correlated with G/C and A/T tracts, respectively, suggesting a trade-off between these two deletion events. GC content evolved at various rates across the phylogeny, with significant increases in monocots and Lamiids, and a decrease in Fabids, implying the effects of some other factors. We hypothesize that variation in plastid GC content might be a mixed strategy of species to optimize fitness in fluctuating climates, partly through influencing the trade-off between AT → GC and GC → AT mutations.

Integrative Study of Genotypic and Phenotypic Diversity in the Eurasian Orchid Genus Neotinea

Knowledge of population variation across species' ranges is a prerequisite for correctly assessing the overall variability of any group of organisms and provides an invaluable basis for unraveling evolutionary history, optimizing taxonomy and devising effective conservation strategies. Here, we examine the genus Neotinea, which represents a relatively recently delimited monophyletic genus of orchids, for which a detailed study of its overall variability was lacking. We applied a suite of biosystematic methods, consisting of flow cytometry, multivariate and geometric morphometrics, and analysis of genomic SNP data, to identify phylogenetic lineages within the genus, to delineate phenotypic variation relevant to these lineages, and to identify potential cryptic taxa within lineages. We found clear differentiation into four major lineages corresponding to the groups usually recognized within the genus: Neotinea maculata as a distinct and separate taxon, the Neotinea lactea group comprising two Mediterranean taxa N. lactea and Neotinea conica, the Neotinea ustulata group comprising two phenologically distinct varieties, and the rather complex Neotinea tridentata group comprising two major lineages and various minor lineages of unclear taxonomic value. N. conica constitutes both a monophyletic group within N. lactea and a distinct phenotype within the genus and merits its proposed subspecies-level recognition. By contrast, the spring and summer flowering forms of N. ustulata (var. ustulata and var. aestivalis) were confirmed to be distinct only morphologically, not phylogenetically. The most complex pattern emerged in the N. tridentata group, which splits into two main clades, one containing lineages from the Balkans and eastern Mediterranean and the other consisting of plants from Central Europe and the central Mediterranean. These individual lineages differ in genome size and show moderate degrees of morphological divergence. The tetraploid Neotinea commutata is closely related to the N. tridentata group, but our evidence points to an auto- rather than an allopolyploid origin. Our broad methodological approach proved effective in recognizing cryptic lineages among the orchids, and we propose the joint analysis of flow cytometric data on genome size and endopolyploidy as a useful and beneficial marker for delineating orchid species with partial endoreplication.

Substantial intraspecific genome size variation in golden-brown algae and its phenotypic consequences

BACKGROUND AND AIMS

While nuclear DNA content variation and its phenotypic consequences have been well described for animals, vascular plants and macroalgae, much less about this topic is known regarding unicellular algae and protists in general. The dearth of data is especially pronounced when it comes to intraspecific genome size variation. This study attempts to investigate the extent of intraspecific variability in genome size and its adaptive consequences in a microalgal species.

METHODS

Propidium iodide flow cytometry was used to estimate the absolute genome size of 131 strains (isolates) of the golden-brown alga Synura petersenii (Chrysophyceae, Stramenopiles), identified by identical internal transcribed spacer (ITS) rDNA barcodes. Cell size, growth rate and genomic GC content were further assessed on a sub-set of strains. Geographic location of 67 sampling sites across the Northern hemisphere was used to extract climatic database data and to evaluate the ecogeographical distribution of genome size diversity.

KEY RESULTS

Genome size ranged continuously from 0.97 to 2.02 pg of DNA across the investigated strains. The genome size was positively associated with cell size and negatively associated with growth rate. Bioclim variables were not correlated with genome size variation. No clear trends in the geographical distribution of strains of a particular genome size were detected, and strains of different genome size occasionally coexisted at the same locality. Genomic GC content was significantly associated only with genome size via a quadratic relationship.

CONCLUSIONS

Genome size variability in S. petersenii was probably triggered by an evolutionary mechanism operating via gradual changes in genome size accompanied by changes in genomic GC content, such as, for example, proliferation of transposable elements. The variation was reflected in cell size and relative growth rate, possibly with adaptive consequences.

Comparing DNA replication programs reveals large timing shifts at centromeres of endocycling cells in maize roots

Plant cells undergo two types of cell cycles–the mitotic cycle in which DNA replication is coupled to mitosis, and the endocycle in which DNA replication occurs in the absence of cell division. To investigate DNA replication programs in these two types of cell cycles, we pulse labeled intact root tips of maize (Zea mays) with 5-ethynyl-2’-deoxyuridine (EdU) and used flow sorting of nuclei to examine DNA replication timing (RT) during the transition from a mitotic cycle to an endocycle. Comparison of the sequence-based RT profiles showed that most regions of the maize genome replicate at the same time during S phase in mitotic and endocycling cells, despite the need to replicate twice as much DNA in the endocycle and the fact that endocycling is typically associated with cell differentiation. However, regions collectively corresponding to 2% of the genome displayed significant changes in timing between the two types of cell cycles. The majority of these regions are small with a median size of 135 kb, shift to a later RT in the endocycle, and are enriched for genes expressed in the root tip. We found larger regions that shifted RT in centromeres of seven of the ten maize chromosomes. These regions covered the majority of the previously defined functional centromere, which ranged between 1 and 2 Mb in size in the reference genome. They replicate mainly during mid S phase in mitotic cells but primarily in late S phase of the endocycle. In contrast, the immediately adjacent pericentromere sequences are primarily late replicating in both cell cycles. Analysis of CENH3 enrichment levels in 8C vs 2C nuclei suggested that there is only a partial replacement of CENH3 nucleosomes after endocycle replication is complete. The shift to later replication of centromeres and possible reduction in CENH3 enrichment after endocycle replication is consistent with a hypothesis that centromeres are inactivated when their function is no longer needed.

In traditional cell division, or mitosis, a cell’s genetic material is duplicated and then split between two daughter cells. In contrast, in some specialized cell types, the DNA is duplicated a second time without an intervening division step, resulting in cells that carry twice as much DNA. This phenomenon, which is called the endocycle, is common during plant development. At each step, DNA replication follows an ordered program in which highly compacted DNA is unraveled and replicated in sections at different times during the synthesis (S) phase. In plants, it is unclear whether traditional and endocycle programs are the same, especially since endocycling cells are typically in the process of differentiation. Using root tips of maize, we found that in comparison to replication in the mitotic cell cycle, there is a small portion of the genome whose replication in the endocycle is shifted in time, usually to later in S phase. Some of these regions are scattered around the genome and mostly coincide with active genes. However, the most prominent shifts occur in centromeres. The shift to later replication in centromeres is noteworthy because they orchestrate the process of separating duplicated chromosomes into daughter cells, a function that is not needed in the endocycle.

Early Diverging and Core Bromelioideae (Bromeliaceae) Reveal Contrasting Patterns of Genome Size Evolution and Polyploidy

The subfamily Bromelioideae is one of the most diverse groups among the neotropical Bromeliaceae. Previously, key innovations have been identified which account for the extraordinary radiation and species richness of this subfamily, especially in the so-called core Bromelioideae. However, in order to extend our understanding of the evolutionary mechanisms, the genomic mechanisms (e.g. polyploidy, dysploidy) that potentially underlie this accelerated speciation also need to be tested. Here, using PI and DAPI staining and flow cytometry we estimated genome size and GC content of 231 plants covering 30 genera and 165 species and combined it with published data. The evolutionary and ecological significance of all three genomic characters was tested within a previously generated dated phylogenetic framework using ancestral state reconstructions, comparative phylogenetic methods, and multiple regressions with climatic variables. The absolute genome size (2C) of Bromelioideae varied between 0.59 and 4.11 pg, and the GC content ranged between 36.73 and 41.43%. The monoploid genome sizes (Cx) differed significantly between core and early diverging lineages. The occurrence of dysploidy and polyploidy was, with few exceptions, limited to the phylogenetically isolated early diverging tank-less lineages. For Cx and GC content Ornstein-Uhlenbeck models outperformed the Brownian motion models suggesting adaptive potential linked to the temperature conditions. 2C-values revealed different rates of evolution in core and early diverging lineages also related to climatic conditions. Our results suggest that polyploidy is not associated with higher net diversification and fast radiation in core bromelioids. On the other hand, although coupled with higher extinction rates, dysploidy, polyploidy, and resulting genomic reorganizations might have played a role in the survival of the early diverging bromelioids in hot and arid environments.