Is There an Upper Limit to Genome Size?

Metabolomic Diversity and Defensive Phenolic Compounds in Cloud Forest Ferns

The few current metabolomic studies on ferns are mostly restricted to a single species or focused on specific compounds. We performed an untargeted metabolomic study on six of the most common fern species from the cloud forest, followed by a targeted analysis of 64 phenolic compounds, many of which have been associated with herbivore defense. The untargeted analysis revealed a total of 232 putative identified metabolites from 463 to 1427 signals per fern species, each with its proper chemical signature but not necessarily correlated to their phylogenetic relationship. The flavonoid, flavone, and flavonol biosynthesis were the most expressed pathways in all species except for Marattia laxa. Fern species also differed strongly in the concentrations of the 10 detected phenolic compounds. Our results show that ferns, including the most ancestral species, such as M. laxa, display a high metabolomic diversity comparable to seed plants. Each fern species held a different combination of defensive phenolic compounds. Further research is needed to explore the metabolic diversity, to identify the biochemical defenses of ferns, and, in particular, to detect the chemical compounds that act against their specific herbivorous insects in the cloud forest ecosystem.

Influence of Ty3/gypsy and Ty1/copia LTR-retrotransposons on the large genomes of Alstroemeriaceae: genome landscape of Bomarea edulis (Tussac) Herb

Repetitive elements are the main components of many plant genomes and play a crucial role in the variation of genome size and structure, ultimately impacting species diversification and adaptation. Alstroemeriaceae exhibits species with large genomes, not attributed to polyploidy. In this study, we analysed the repetitive fraction of the genome of Bomarea edulis through low-coverage sequencing and in silico characterization, and compared it to the repeats of Alstroemeria longistaminea, a species from a sister genus that has been previously characterized. LTR-retrotransposons were identified as the most abundant elements in the B. edulis genome (50.22%), with significant variations in abundance for specific lineages between the two species. The expansion of the B. edulis genome was likely due to three main lineages of LTR retrotransposons, Ty3/gypsy Tekay and Retand and Ty1/copia SIRE, all represented by truncated elements which were probably active in the past. Furthermore, the proportion of satDNA (~ 7%) was six times higher in B. edulis compared to A. longistaminea, with most families exhibiting a dispersed, uniform distribution in the genome. SatDNAs, thus, contributed to some extent to genome obesity. Despite diverging around 29 Mya, both species still share some satDNA families and retrotransposons. However, differences in repeat abundances and sequence variants led to genome differentiation despite their similar sizes and structure.

Germline ecology: Managed herds, tolerated flocks, and pest control

Multicopy sequences evolve adaptations for increasing their copy number within nuclei. The activities of multicopy sequences under constraints imposed by cellular and organismal selection result in a rich intranuclear ecology in germline cells. Mitochondrial and ribosomal DNA are managed as domestic herds subject to selective breeding by the genes of the single-copy genome. Transposable elements lead a peripatetic existence in which they must continually move to new sites to keep ahead of inactivating mutations at old sites and undergo exponential outbreaks when the production of new copies exceeds the rate of inactivation of old copies. Centromeres become populated by repeats that do little harm. Organisms with late sequestration of germ cells tend to evolve more "junk" in their genomes than organisms with early sequestration of germ cells.

Expansions and contractions of repetitive DNA elements reveal contrasting evolutionary responses to the polyploid genome shock hypothesis in Brachypodium model grasses

Brachypodium grass species have been selected as model plants for functional genomics of grass crops, and to elucidate the origins of allopolyploidy and perenniality in monocots, due to their small genome sizes and feasibility of cultivation. However, genome sizes differ greatly between diploid or polyploid Brachypodium lineages. We have used genome skimming sequencing data to uncover the composition, abundance, and phylogenetic value of repetitive elements in 44 representatives of the major Brachypodium lineages and cytotypes. We also aimed to test the possible mechanisms and consequences of the "polyploid genome shock hypothesis" (PGSH) under three different evolutionary scenarios of variation in repeats and genome sizes of Brachypodium allopolyploids. Our data indicated that the proportion of the genome covered by the repeatome in the Brachypodium species showed a 3.3-fold difference between the highest content of B. mexicanum-4x (67.97%) and the lowest of B. stacei-2x (20.77%), and that changes in the sizes of their genomes were a consequence of gains or losses in their repeat elements. LTR-Retand and Tekay retrotransposons were the most frequent repeat elements in the Brachypodium genomes, while Ogre retrotransposons were found exclusively in B. mexicanum. The repeatome phylogenetic network showed a high topological congruence with plastome and nuclear rDNA and transcriptome trees, differentiating the ancestral outcore lineages from the recently evolved core-perennial lineages. The 5S rDNA graph topologies had a strong match with the ploidy levels and nature of the subgenomes of the Brachypodium polyploids. The core-perennial B. sylvaticum presents a large repeatome and characteristics of a potential post-polyploid diploidized origin. Our study evidenced that expansions and contractions in the repeatome were responsible for the three contrasting responses to the PGSH. The exacerbated genome expansion of the ancestral allotetraploid B. mexicanum was a consequence of chromosome-wide proliferation of TEs and not of WGD, the additive repeatome pattern of young allotetraploid B. hybridum of stabilized post-WGD genome evolution, and the genomecontraction of recent core-perennials polyploids (B. pinnatum, B. phoenicoides) of repeat losses through recombination of these highly hybridizing lineages. Our analyses have contributed to unraveling the evolution of the repeatome and the genome size variation in model Brachypodium grasses.

A 160 Gbp fork fern genome shatters size record for eukaryotes

Vascular plants are exceptional among eukaryotes due to their outstanding genome size diversity which ranges ∼2,400-fold, including the largest genome so far recorded in the angiosperm Paris japonica (148.89 Gbp/1C). Despite available data showing that giant genomes are restricted across the Tree of Life, the biological limits to genome size expansion remain to be established. Here, we report the discovery of an even larger eukaryotic genome in Tmesipteris oblanceolata, a New Caledonian fork fern. At 160.45 Gbp/1C, this record-breaking genome challenges current understanding and opens new avenues to explore the evolutionary dynamics of genomic gigantism.

Genome engineering on size reduction and complexity simplification: A review

BACKGROUND

Genome simplification is an important topic in the field of life sciences that has attracted attention from its conception to the present day. It can help uncover the essential components of the genome and, in turn, shed light on the underlying operating principles of complex biological systems. This has made it a central focus of both basic and applied research in the life sciences. With the recent advancements in related technologies and our increasing knowledge of the genome, now is an opportune time to delve into this topic.

AIM OF REVIEW

Our review investigates the progress of genome simplification from two perspectives: genome size reduction and complexity simplification. In addition, we provide insights into the future development trends of genome simplification.

KEY SCIENTIFIC CONCEPTS OF REVIEW

Reducing genome size requires eliminating non-essential elements as much as possible. This process has been facilitated by advances in genome manipulation and synthesis techniques. However, we still need a better and clearer understanding of living systems to reduce genome complexity. As there is a lack of quantitative and clearly defined standards for this task, we have opted to approach the topic from various perspectives and present our findings accordingly.

An adventurous journey toward and away from fern apomixis: Insights from genome size and spore abortion patterns

PREMISE

Apomixis in ferns is relatively common and obligatory. Sterile hybrids may restore fertility via apomixis at a cost of long-term genetic stagnation. In this study, we outlined apomixis as a possible temporary phase leading to sexuality and analyzed factors relating to transitioning to and away from apomixis, such as unreduced and reduced spore formation in apomict and apo-sex hybrid ferns.

METHODS

We analyzed the genome size of 15 fern species or hybrids ("taxa") via flow cytometry. The number of reduced and unreduced gametophytes was established as a proxy for viable spore formation of either type. We also calculated the spore abortion ratio (sign of reduced spores) in several taxa, including the apo-sex hybrid Dryopteris × critica and its 16 apomictically formed offspring.

RESULTS

Four of 15 sampled taxa yielded offspring variable in genome size. Specifically, each variable taxon formed one viable reduced plant among 12-451 sampled gametophytes per taxon. Thus, haploid spore formation in the studied apomicts was very rare but possible. Spore abortion analyses indicated gradually decreasing abortion (haploid spore formation) over time. In Dryopteris × critica, abortion decreased from 93.8% to mean 89.5% in one generation.

CONCLUSIONS

Our results support apomixis as a transitionary phase toward sexuality. Newly formed apomicts hybridize with sexual relatives and continue to form haploid spores early on. Thus, they may get the genomic content necessary for regular meiosis and restore sexuality. If the missing relative goes extinct, the lineage gets locked into apomixis as may be the case with the Dryopteris affinis complex.

Genome composition and GC content influence loci distribution in reduced representation genomic studies

BACKGROUND

Genomic architecture is a key evolutionary trait for living organisms. Due to multiple complex adaptive and neutral forces which impose evolutionary pressures on genomes, there is a huge variability of genomic features. However, their variability and the extent to which genomic content determines the distribution of recovered loci in reduced representation sequencing studies is largely unexplored.

RESULTS

Here, by using 80 genome assemblies, we observed that whereas plants primarily increase their genome size by expanding their intergenic regions, animals expand both intergenic and intronic regions, although the expansion patterns differ between deuterostomes and protostomes. Loci mapping in introns, exons, and intergenic categories obtained by in silico digestion using 2b-enzymes are positively correlated with the percentage of these regions in the corresponding genomes, suggesting that loci distribution mostly mirrors genomic architecture of the selected taxon. However, exonic regions showed a significant enrichment of loci in all groups regardless of the used enzyme. Moreover, when using selective adaptors to obtain a secondarily reduced loci dataset, the percentage and distribution of retained loci also varied. Adaptors with G/C terminals recovered a lower percentage of selected loci, with a further enrichment of exonic regions, while adaptors with A/T terminals retained a higher percentage of loci and slightly selected more intronic regions than expected.

CONCLUSIONS

Our results highlight how genome composition, genome GC content, RAD enzyme choice and use of base-selective adaptors influence reduced genome representation techniques. This is important to acknowledge in population and conservation genomic studies, as it determines the abundance and distribution of loci.

DNA Damage, Genome Stability, and Adaptation: A Question of Chance or Necessity?

DNA damage causes the mutations that are the principal source of genetic variation. DNA damage detection and repair mechanisms therefore play a determining role in generating the genetic diversity on which natural selection acts. Speciation, it is commonly assumed, occurs at a rate set by the level of standing allelic diversity in a population. The process of speciation is driven by a combination of two evolutionary forces: genetic drift and ecological selection. Genetic drift takes place under the conditions of relaxed selection, and results in a balance between the rates of mutation and the rates of genetic substitution. These two processes, drift and selection, are necessarily mediated by a variety of mechanisms guaranteeing genome stability in any given species. One of the outstanding questions in evolutionary biology concerns the origin of the widely varying phylogenetic distribution of biodiversity across the Tree of Life and how the forces of drift and selection contribute to shaping that distribution. The following examines some of the molecular mechanisms underlying genome stability and the adaptive radiations that are associated with biodiversity and the widely varying species richness and evenness in the different eukaryotic lineages.

Characterization and analysis of multi-organ full-length transcriptomes in Sphaeropteris brunoniana and Alsophila latebrosa highlight secondary metabolism and chloroplast RNA editing pattern of tree ferns

BACKGROUND

Sphaeropteris brunoniana and Alsophila latebrosa are both old relict and rare tree ferns, which have experienced the constant changes of climate and environment. However, little is known about their high-quality genetic information and related research on environmental adaptation mechanisms of them. In this study, combined with PacBio and Illumina platforms, transcriptomic analysis was conducted on the roots, rachis, and pinna of S. brunoniana and A. latebrosa to identify genes and pathways involved in environmental adaptation. Additionally, based on the transcriptomic data of tree ferns, chloroplast genes were mined to analyze their gene expression levels and RNA editing events.

RESULTS

In the study, we obtained 11,625, 14,391 and 10,099 unigenes of S. brunoniana root, rachis, and pinna, respectively. Similarly, a total of 13,028, 11,431 and 12,144 unigenes were obtained of A. latebrosa root, rachis, and pinna, respectively. According to the enrichment results of differentially expressed genes, a large number of differentially expressed genes were enriched in photosynthesis and secondary metabolic pathways of S. brunoniana and A. latebrosa. Based on gene annotation results and phenylpropanoid synthesis pathways, two lignin synthesis pathways (H-lignin and G-lignin) were characterized of S. brunoniana. Among secondary metabolic pathways of A. latebrosa, three types of WRKY transcription factors were identified. Additionally, based on transcriptome data obtained in this study, reported transcriptome data, and laboratory available transcriptome data, positive selection sites were identified from 18 chloroplast protein-coding genes of four tree ferns. Among them, RNA editing was found in positive selection sites of four tree ferns. RNA editing affected the protein secondary structure of the rbcL gene. Furthermore, the expression level of chloroplast genes indicated high expression of genes related to the chloroplast photosynthetic system in all four species.

CONCLUSIONS

Overall, this work provides a comprehensive transcriptome resource of S. brunoniana and A. latebrosa, laying the foundation for future tree fern research.

Interploidy gene flow via a 'pentaploid bridge' and ploidy reduction in Cystopteris fragilis fern complex (Cystopteridaceae: Polypodiales)

KEY MESSAGE

Our results indicate the existence of interploidy gene flow in Cystopteris fragilis, resulting in sexual triploid and diploid gametophytes from pentaploid parents. Similar evolutionary dynamics might operate in other fern complexes and need further investigation. Polyploidization and hybridization are a key evolutionary processes in ferns. Here, we outline an interploidy gene flow pathway operating in the polyploid Cystopteris fragilis complex. The conditions necessary for the existence of this pathway were tested. A total of 365 C. fragilis individuals were collected covering representatives of all three predominant ploidy levels (tetraploid, pentaploid, and hexaploid), cultivated, had their ploidy level estimated by flow cytometry, and their spores collected. The spores, as well as gametophytes and sporophytes established from them, were analysed by flow cytometry. Spore abortion rate was also estimated. In tetraploids, we observed the formation of unreduced (tetraploid) spores (ca 2%). Collected pentaploid individuals indicate ongoing hybridization between ploidy levels. Pentaploids formed up to 52% viable spores, ca 79% of them reduced, i.e. diploid and triploid. Reduced spores formed viable gametophytes, and, in the case of triploids, filial hexaploid sporophytes, showing evidence of sexual reproduction. Some tetraploid sporophytes reproduce apomictically (based on uniform ploidy of their metagenesis up to filial sporophytes). Triploid and diploid gametophytes from pentaploid parents are able to mate among themselves, or with "normal" reduced gametophytes from the sexual tetraploid sporophytes (the dominant ploidy level in the sporophytes in this populations), to produce tetraploid, pentaploid, and hexaploid sporophytes, allowing for geneflow from the pentaploids to both the tetraploid and hexaploid populations. Similar evolutionary dynamics might operate in other fern complexes and need further investigation.

Genomic architecture constrains macromolecular allocation in dinoflagellates

Dinoflagellate genomes have a unique architecture that may constrain their physiological and biochemical responsiveness to environmental stressors. Here we quantified how nitrogen (N) starvation influenced macromolecular allocation and C:N:P of three photosynthetic marine dinoflagellates, representing different taxonomic classes and genome sizes. Dinoflagellates respond to nitrogen starvation by decreasing cellular nitrogen, protein and RNA content, but unlike many other eukaryotic phytoplankton examined RNA:protein is invariant. Additionally, 2 of the 3 species exhibit increases in cellular phosphorus and very little change in cellular carbon with N-starvation. As a consequence, N starvation induces moderate increases in C:N, but extreme decreases in N:P and C:P, relative to diatoms. Dinoflagellate DNA content relative to total C, N and P is much higher than similar sized diatoms, but similar to very small photosynthetic picoeukaryotes such as Ostreococcus. In aggregate these results indicate the accumulation of phosphate stores may be an important strategy employed by dinoflagellates to meet P requirements associated with the maintenance and replication of their large genomes.

MOMS: A pipeline for scaffolding using multi-optical maps

Here, we report a new multi-optical maps scaffolder (MOMS) aiming at utilizing complementary information among optical maps labelled by distinct enzymes. This pipeline was designed for data structure organization, scaffolding by path traversal, gap-filling and molecule reuse of optical maps. Our testing showed that this pipeline has uncapped enzyme tolerance in scaffolding. This means that there are no inbuilt limits as to the number of maps generated by different enzymes that can be utilized by MOMS. For the genome assembly of the human GM12878 cell line, MOMS significantly improved the contiguity and completeness with an up to 144-fold increase of scaffold N50 compared with initial assemblies. Benchmarking on the genomes of human and O. sativa showed that MOMS is more effective and robust compared with other optical-map-based scaffolders. We believe this pipeline will contribute to high-fidelity chromosome assembly and chromosome-level evolutionary analysis.

Relaxation of Natural Selection in the Evolution of the Giant Lungfish Genomes

Nonadaptive hypotheses on the evolution of eukaryotic genome size predict an expansion when the process of purifying selection becomes weak. Accordingly, species with huge genomes, such as lungfish, are expected to show a genome-wide relaxation signature of selection compared with other organisms. However, few studies have empirically tested this prediction using genomic data in a comparative framework. Here, we show that 1) the newly assembled transcriptome of the Australian lungfish, Neoceratodus forsteri, is characterized by an excess of pervasive transcription, or transcriptional leakage, possibly due to suboptimal transcriptional control, and 2) a significant relaxation signature in coding genes in lungfish species compared with other vertebrates. Based on these observations, we propose that the largest known animal genomes evolved in a nearly neutral scenario where genome expansion is less efficiently constrained.

Genetic basis of genome size variation of wheat

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

Genomic copy number variability at the genus, species and population levels impacts in situ ecological analyses of dinoflagellates and harmful algal blooms

The application of meta-barcoding, qPCR, and metagenomics to aquatic eukaryotic microbial communities requires knowledge of genomic copy number variability (CNV). CNV may be particularly relevant to functional genes, impacting dosage and expression, yet little is known of the scale and role of CNV in microbial eukaryotes. Here, we quantify CNV of rRNA and a gene involved in Paralytic Shellfish Toxin (PST) synthesis (sxtA4), in 51 strains of 4 Alexandrium (Dinophyceae) species. Genomes varied up to threefold within species and ~7-fold amongst species, with the largest (A. pacificum, 130 ± 1.3 pg cell-1 /~127 Gbp) in the largest size category of any eukaryote. Genomic copy numbers (GCN) of rRNA varied by 6 orders of magnitude amongst Alexandrium (102- 108 copies cell-1) and were significantly related to genome size. Within the population CNV of rRNA was 2 orders of magnitude (105 - 107 cell-1) in 15 isolates from one population, demonstrating that quantitative data based on rRNA genes needs considerable caution in interpretation, even if validated against locally isolated strains. Despite up to 30 years in laboratory culture, rRNA CNV and genome size variability were not correlated with time in culture. Cell volume was only weakly associated with rRNA GCN (20-22% variance explained across dinoflagellates, 4% in Gonyaulacales). GCN of sxtA4 varied from 0-102 copies cell-1, was significantly related to PSTs (ng cell-1), displaying a gene dosage effect modulating PST production. Our data indicate that in dinoflagellates, a major marine eukaryotic group, low-copy functional genes are more reliable and informative targets for quantification of ecological processes than unstable rRNA genes.

Tracing 100 million years of grass genome evolutionary plasticity

Grasses derive from a family of monocotyledonous plants that includes crops of major economic importance such as wheat, rice, sorghum and barley, sharing a common ancestor some 100 million years ago. The genomic attributes of plant adaptation remain obscure and the consequences of recurrent whole genome duplications (WGD) or polyploidization events, a major force in plant evolution, remain largely speculative. We conducted a comparative analysis of omics data from ten grass species to unveil structural (inversions, fusions, fissions, duplications, substitutions) and regulatory (expression and methylation) basis of genome plasticity, as possible attributes of plant long lasting evolution and adaptation. The present study demonstrates that diverged polyploid lineages sharing a common WGD event often present the same patterns of structural changes and evolutionary dynamics, but these patterns are difficult to generalize across independent WGD events as a result of non-WGD factors such as selection and domestication of crops. Polyploidy is unequivocally linked to the evolutionary success of grasses during the past 100 million years, although it remains difficult to attribute this success to particular genomic consequences of polyploidization, suggesting that polyploids harness the potential of genome duplication, at least partially, in lineage-specific ways. Overall, the present study clearly demonstrates that post-polyploidization reprogramming is more complex than traditionally reported in investigating single species and calls for a critical and comprehensive comparison across independently polyploidized lineages.

Molecular mechanisms of adaptive evolution in wild animals and plants

Wild animals and plants have developed a variety of adaptive traits driven by adaptive evolution, an important strategy for species survival and persistence. Uncovering the molecular mechanisms of adaptive evolution is the key to understanding species diversification, phenotypic convergence, and inter-species interaction. As the genome sequences of more and more non-model organisms are becoming available, the focus of studies on molecular mechanisms of adaptive evolution has shifted from the candidate gene method to genetic mapping based on genome-wide scanning. In this study, we reviewed the latest research advances in wild animals and plants, focusing on adaptive traits, convergent evolution, and coevolution. Firstly, we focused on the adaptive evolution of morphological, behavioral, and physiological traits. Secondly, we reviewed the phenotypic convergences of life history traits and responding to environmental pressures, and the underlying molecular convergence mechanisms. Thirdly, we summarized the advances of coevolution, including the four main types: mutualism, parasitism, predation and competition. Overall, these latest advances greatly increase our understanding of the underlying molecular mechanisms for diverse adaptive traits and species interaction, demonstrating that the development of evolutionary biology has been greatly accelerated by multi-omics technologies. Finally, we highlighted the emerging trends and future prospects around the above three aspects of adaptive evolution.

What Has Genomics Taught An Evolutionary Biologist?

Genomics, an interdisciplinary field of biology on the structure, function, and evolution of genomes, has revolutionized many subdisciplines of life sciences, including my field of evolutionary biology, by supplying huge data, bringing high-throughput technologies, and offering a new approach to biology. In this review, I describe what I have learned from genomics and highlight the fundamental knowledge and mechanistic insights gained. I focus on three broad topics that are central to evolutionary biology and beyond-variation, interaction, and selection-and use primarily my own research and study subjects as examples. In the next decade or two, I expect that the most important contributions of genomics to evolutionary biology will be to provide genome sequences of nearly all known species on Earth, facilitate high-throughput phenotyping of natural variants and systematically constructed mutants for mapping genotype-phenotype-fitness landscapes, and assist the determination of causality in evolutionary processes using experimental evolution.

The Plant DNA C-Values Database: A One-Stop Shop for Plant Genome Size Data

Genome size is a plant character with far-reaching implications, ranging from impacts on the financial and computing feasibility of sequencing and assembling genomes all the way to influencing the very ecology and evolution of species. The increasing recognition of the role of genome size in plant science has led to a rising demand for comprehensive and easily accessible sources of genome size data. The Plant DNA C-values database has established itself as a trusted and widely used central hub for users needing to access available plant genome size data, complemented with related cytogenetic (ploidy level) and karyological (chromosome number) information where available. Since its inception in 2001, the database has undergone six major updates to incorporate newly available genome size information, leading to the most recent release (Release 7.1), which comprises data for 12,273 species across all the major land plant and some algal lineages. Here we describe how to use the database efficiently, making use of its different query and filtering settings.

Chromosome size matters: genome evolution in the cyperid clade

BACKGROUND AND AIMS

While variation in genome size and chromosome numbers and their consequences are often investigated in plants, the biological relevance of variation in chromosome size remains poorly known. Here, we examine genome and mean chromosome size in the cyperid clade (families Cyperaceae, Juncaceae and Thurniaceae), which is the largest vascular plant lineage with predominantly holocentric chromosomes.

METHODS

We measured genome size in 436 species of cyperids using flow cytometry, and augment these data with previously published datasets. We then separately compared genome and mean chromosome sizes (2C/2n) amongst the major lineages of cyperids and analysed how these two genomic traits are associated with various environmental factors using phylogenetically informed methods.

KEY RESULTS

We show that cyperids have the smallest mean chromosome sizes recorded in seed plants, with a large divergence between the smallest and largest values. We found that cyperid species with smaller chromosomes have larger geographical distributions and that there is a strong inverse association between mean chromosome size and number across this lineage.

CONCLUSIONS

The distinct patterns in genome size and mean chromosome size across the cyperids might be explained by holokinetic drive. The numerous small chromosomes might function to increase genetic diversity in this lineage where crossovers are limited during meiosis.

Fern Spores-"Ready-to-Use" Standards for Plant Genome Size Estimation Using a Flow Cytometric Approach

Spores and pollen of plants were used as flow cytometric materials to efficiently infer genome sizes. Given this advantage, they hold great potential for various flow cytometric applications, particularly as plant genome size standards. To develop such novel standards, we investigated conditions of pretreatment (bead vortex), buffer, and reliable genome sizes of three fern spore collections-Cibotium taiwanense "Kuo4395", Sphaeropteris lepifera "Tang0001", and Alsophila metteniana "Lee s.n.". Additionally, up to 30 year-old spore collections were obtained from herbarium specimens and from samples stored at 4 °C; their spore nuclei were extracted, and the quality and quantity of these nucleus extractions through storage ages were examined. Nuclear extractions with a longer bead vortex duration or lower spore/bead ratio generally resulted in a higher recovered quantity but a lower quality or purity. For each spore standard, the protocol optimization was determined by their performance in bead vortex conditions, and a 1C genome size was further inferred by linear regression (C. taiwanense "Kuo4395" = 5.058 pg; S. lepifera "Tang0001" = 7.117 pg; and A. metteniana "Lee s.n." = 19.379 pg). Spore nucleus quality and quantity are significantly negatively correlated with storage ages. Nuclear extractions of 10-year-old refrigerated spores remained qualified as a genome size standard; however, none of the herbarium spore collections fit such criteria. Our study is the first to develop and apply dried and refrigerated spores for genome size standards. These standards are ready to use, easy to manipulate, and feature long-term storage in comparison with traditionally used standards of fresh leaves.

Intraspecies characterization of bacteria via evolutionary modeling of protein domains

The ability to detect and characterize bacteria within a biological sample is crucial for the monitoring of infections and epidemics, as well as for the study of human health and its relationship with commensal microorganisms. To this aim, a commonly used technique is the 16S rRNA gene targeted sequencing. PCR-amplified 16S sequences derived from the sample of interest are usually clustered into the so-called Operational Taxonomic Units (OTUs) based on pairwise similarities. Then, representative OTU sequences are compared with reference (human-made) databases to derive their phylogeny and taxonomic classification. Here, we propose a new reference-free approach to define the phylogenetic distance between bacteria based on protein domains, which are the evolving units of proteins. We extract the protein domain profiles of 3368 bacterial genomes and we use an ecological approach to model their Relative Species Abundance distribution. Based on the model parameters, we then derive a new measurement of phylogenetic distance. Finally, we show that such model-based distance is capable of detecting differences between bacteria in cases in which the 16S rRNA-based method fails, providing a possibly complementary approach , which is particularly promising for the analysis of bacterial populations measured by shotgun sequencing.

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

Genomic insights into genetic diploidization in the homosporous fern Adiantum nelumboides

Whole genome duplication has been recognized as a major process in speciation of land plants, especially in ferns. Whereas genome downsizing contributes greatly to the post-genome shock responses of polyploid flowering plants, diploidization of polyploid ferns diverges by maintaining most of the duplicated DNA and is thus expected to be dominated by genic processes. As a consequence, fern genomes provide excellent opportunities to study ecological speciation enforced by expansion of protein families via polyploidy. To test the key predictions of this hypothesis, we reported the de novo genome sequence of Adiantum nelumboides, a tetraploid homosporous fern. The obtained draft genome had a size of 6.27 Gb assembled into 11,767 scaffolds with the contig N50 of 1.37 Mb. Repetitive DNA sequences contributed with about 81.7%, a remarkably high proportion of the genome. With 69,568 the number of predicted protein-coding genes exceeded those reported in most other land plant genomes. Intragenomic synteny analyses recovered 443 blocks with the average block size of 1.29 Mb and the average gene content of 16 genes. The results are consistent with the hypothesis of high ancestral chromosome number, lack of substantial genome downsizing, and dominance of genic diploidization. As expected in the calciphilous plants, a notable number of detected genes were involved in calcium uptake and transport. In summary, the genome sequence of a tetraploid homosporous fern not only provides access to a genomic resource of a derived fern, but also supports the hypothesis of maintenance of high chromosome numbers and duplicated DNA in young polyploid ferns.

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.

Evolutionary Dynamics of the Repeatome Explains Contrasting Differences in Genome Sizes and Hybrid and Polyploid Origins of Grass Loliinae Lineages

The repeatome is composed of diverse families of repetitive DNA that keep signatures on the historical events that shaped the evolution of their hosting species. The cold seasonal Loliinae subtribe includes worldwide distributed taxa, some of which are the most important forage and lawn species (fescues and ray-grasses). The Loliinae are prone to hybridization and polyploidization. It has been observed a striking two-fold difference in genome size between the broad-leaved (BL) and fine-leaved (FL) Loliinae diploids and a general trend of genome reduction of some high polyploids. We have used genome skimming data to uncover the composition, abundance, and potential phylogenetic signal of repetitive elements across 47 representatives of the main Loliinae lineages. Independent and comparative analyses of repetitive sequences and of 5S rDNA loci were performed for all taxa under study and for four evolutionary Loliinae groups [Loliinae, Broad-leaved (BL), Fine-leaved (FL), and Schedonorus lineages]. Our data showed that the proportion of the genome covered by the repeatome in the Loliinae species was relatively high (average ∼ 51.8%), ranging from high percentages in some diploids (68.7%) to low percentages in some high-polyploids (30.7%), and that changes in their genome sizes were likely caused by gains or losses in their repeat elements. Ty3-gypsy Retand and Ty1-copia Angela retrotransposons were the most frequent repeat families in the Loliinae although the relatively more conservative Angela repeats presented the highest correlation of repeat content with genome size variation and the highest phylogenetic signal of the whole repeatome. By contrast, Athila retrotransposons presented evidence of recent proliferations almost exclusively in the Lolium clade. The repeatome evolutionary networks showed an overall topological congruence with the nuclear 35S rDNA phylogeny and a geographic-based structure for some lineages. The evolution of the Loliinae repeatome suggests a plausible scenario of recurrent allopolyploidizations followed by diploidizations that generated the large genome sizes of BL diploids as well as large genomic rearrangements in highly hybridogenous lineages that caused massive repeatome and genome contractions in the Schedonorus and Aulaxyper polyploids. Our study has contributed to disentangling the impact of the repeatome dynamics on the genome diversification and evolution of the Loliinae grasses.

Pervasive genome duplications across the plant tree of life and their links to major evolutionary innovations and transitions

Whole-genome duplication (WGD) has occurred repeatedly during plant evolution and diversification, providing genetic layers for evolving new functions and phenotypes. Advances in long-read sequencing technologies have enabled sequencing and assembly of over 1000 plant genomes spanning nearly 800 species, in which a large set of ancient WGDs has been uncovered. Here, we review the recently reported WGDs that occurred in major plant lineages and key evolutionary positions, and highlight their contributions to morphological innovation and adaptive evolution. Current gaps and challenges in integrating enormous volumes of sequenced plant genomes, accurately inferring WGDs, and developing web-based analysis tools are emphasized. Looking to the future, ambitious genome sequencing projects and global efforts may substantially recapitulate the plant tree of life based on broader sampling of phylogenetic diversity, reveal much of the timetable of ancient WGDs, and address the biological significance of WGDs in plant adaptation and radiation.

Genome size evolution in the beetle genus Diabrotica

Diabrocite corn rootworms are one of the most economically significant pests of maize in the United States and Europe and an emerging model for insect-plant interactions. Genome sizes of several species in the genus Diabrotica were estimated using flow cytometry along with that of Acalymma vittatum as an outgroup. Genome sizes ranged between 1.56 and 1.64 gigabase pairs and between 2.26 and 2.59 Gb, respectively, for the Diabrotica subgroups fucata and virgifera; the Acalymma vittatum genome size was around 1.65 Gb. This result indicated that a substantial increase in genome size occurred in the ancestor of the virgifera group. Further analysis of the fucata group and the virgifera group genome sequencing reads indicated that the genome size difference between the Diabrotica subgroups could be attributed to a higher content of transposable elements, mostly miniature inverted-transposable elements and gypsy-like long terminal repeat retroelements.

Full-length transcriptome analysis of Adiantum flabellulatum gametophyte

Ferns are important components of plant communities on earth, but their genomes are generally very large, with many redundant genes, making whole genome sequencing of ferns prohibitively expensive and time-consuming. This means there is a significant lack of fern reference genomes, making molecular biology research difficult. The gametophytes of ferns can survive independently, are responsible for sexual reproduction and the feeding of young sporophytes, and play an important role in the alternation of generations. For this study, we selected Adiantum flabellulatum as it has both ornamental and medicinal value and is also an indicator plant of acidic soil. The full-length transcriptome sequencing of its gametophytes was carried out using PacBio three-generation sequencing technology. A total of 354,228 transcripts were obtained, and 231,705 coding sequences (CDSs) were predicted, including 5,749 transcription factors (TFs), 2,214 transcription regulators (TRs) and 4,950 protein kinases (PKs). The transcripts annotated by non-redundant protein sequence database (NR), Kyoto encyclopedia of genes and genomes (KEGG), eukaryotic ortholog groups (KOG), Swissprot, protein family (Pfma), nucleotide sequence database (NT) and gene ontology (GO) were 251,501, 197,474, 193,630, 194,639, 195,956, 113,069 and 197,883, respectively. In addition, 138,995 simple sequence repeats (SSRs) and 111,793 long non-coding RNAs (lncRNAs) were obtained. We selected nine chlorophyll synthase genes for qRT-PCR, and the results showed that the full-length transcript sequences and the annotation information were reliable. This study can provide a reference gene set for subsequent gene expression quantification.

Genome size evolution of the extant lycophytes and ferns

Ferns and lycophytes have remarkably large genomes. However, little is known about how their genome size evolved in fern lineages. To explore the origins and evolution of chromosome numbers and genome size in ferns, we used flow cytometry to measure the genomes of 240 species (255 samples) of extant ferns and lycophytes comprising 27 families and 72 genera, of which 228 species (242 samples) represent new reports. We analyzed correlations among genome size, spore size, chromosomal features, phylogeny, and habitat type preference within a phylogenetic framework. We also applied ANOVA and multinomial logistic regression analysis to preference of habitat type and genome size. Using the phylogeny, we conducted ancestral character reconstruction for habitat types and tested whether genome size changes simultaneously with shifts in habitat preference. We found that 2C values had weak phylogenetic signal, whereas the base number of chromosomes (x) had a strong phylogenetic signal. Furthermore, our analyses revealed a positive correlation between genome size and chromosome traits, indicating that the base number of chromosomes (x), chromosome size, and polyploidization may be primary contributors to genome expansion in ferns and lycophytes. Genome sizes in different habitat types varied significantly and were significantly correlated with habitat types; specifically, multinomial logistic regression indicated that species with larger 2C values were more likely to be epiphytes. Terrestrial habitat is inferred to be ancestral for both extant ferns and lycophytes, whereas transitions to other habitat types occurred as the major clades emerged. Shifts in habitat types appear be followed by periods of genomic stability. Based on these results, we inferred that habitat type changes and multiple whole-genome duplications have contributed to the formation of large genomes of ferns and their allies during their evolutionary history.

Dependence of genome size and copy number of rRNA gene on cell volume in dinoflagellates

Dinoflagellates are an ecologically important group of protists in aquatic environment and have evolved many unusual and enigmatic genomic features such as immense genome sizes, high repeated genes, and a large portion of hydroxymethyluracil in DNA. Although previous studies have observed positive correlations between the large subunit (LSU) rRNA gene copy number and genome size of a variety of eukaryotic organisms (e.g. higher plants and animals), or between cell volume and LSU rRNA gene copy number, and/or between genome size and cell size, which suggests a possible co-evolution among these three features in different lineages of life, it remains an open question regarding the relationships among these three parameters in dinoflagellates. For the first time, we estimated the copy numbers of the LSU rRNA gene, the genome sizes, and cell volumes within a broad range of dinoflagellates (covering 15 species of 11 genera) using single-cell qPCR-based assay (determining LSU rRNA gene copy number), FlowCAM (cell volume measurement), and ultraviolet spectrophotometry (genome size estimation). The measured copy number of LSU rRNA gene ranged from 398 ± 184 (Prorocentrum minimum) to 152,078 ± 33,555 copies•cell^-1 (Alexandrium pacificum), while the genome size and the cell volume ranged from 5.6 ± 0.2 (Karlodinium veneficum) to 853 ± 19.9 pg•cell^-1 (Pseliodinium pirum), and from 1,070 ± 225 (Kar. veneficum) to 168,474 ± 124,180 μm³ (Ps. pirum), respectively. Together with the three parameters measured in literature, there are significant positive linear correlations between LSU rRNA gene copy numbers and genome sizes, cell volumes and LSU rRNA gene copy numbers, and between genome sizes and cell volumes via comparisons of multi-model regression analyses, suggesting a dependence of genome size and rRNA gene copy number on the cell volumes of dinoflagellates. Validation of the measurement methods was conducted via comparisons between reported data in the literature and that predicted using the linear equations we obtained, and between genome size measured by flow cytometry (FCM) and ultraviolet spectrophotometry (Nanodrop). These results provide insightful understandings of dinoflagellate evolution in terms of the relationships among genomes, gene copy number, and cell volume, and of rRNA gene-based studies in intra-populational and intra-individual genetic diversity, taxonomy, and diversity assessment in the environment of dinoflagellates. The results also provide a dataset useful for reads calibration in environmental metabarcoding studies of dinoflagellates and selection of candidate species for whole genome sequencing.

Exploring environmental selection on genome size in angiosperms

Angiosperms show a remarkable range in genome size (GS), yet most species have small genomes, despite the frequency of polyploidy and repeat amplification in the ancestries of most lineages. It has been suggested that larger genomes incur costs that have driven selection for GS reduction, although the nature of these costs and how they might impact selection remain unclear. We explore potential costs of increased GS encompassing impacts on minimum cell size with consequences for photosynthesis and water-use efficiency and effects of greater nitrogen and phosphorus demands of the nucleus leading to more severe trade-offs with photosynthesis. We suggest that nutrient-, water-, and/or CO₂-stressed conditions might favour species with smaller genomes, with implications for species' ecological and evolutionary dynamics.

Karyotype Reorganization in Wheat-Rye Hybrids Obtained via Unreduced Gametes: Is There a Limit to the Chromosome Number in Triticale?

To date, few data have been accumulated on the contribution of meiotic restitution to the formation of Triticum aestivum hybrid karyotypes. In this study, based on FISH and C-banding, karyotype reorganization was observed in three groups of F₅ wheat–rye hybrids 1R(1A) × R. Aberrations, including aneuploidy, telocentrics, and Robertsonian translocations, were detected in all groups. Some of the Group 1 plants and all of the Group 2 plants only had a 4R4R pair (in addition to 1R1R), which was either added or substituted for its homeolog in ABD subgenomes. In about 82% of meiocytes, 4R4R formed bivalents, which indicates its competitiveness. The rest of the Group 1 plants had 2R and 7R chromosomes in addition to 1R1R. Group 3 retained all their rye chromosomes, with a small aneuploidy on the wheat chromosomes. A feature of the meiosis in the Group 3 plants was asynchronous cell division and omission of the second division. Diploid gametes did not form because of the significant disturbances during gametogenesis. As a result, the frequency of occurrence of the formed dyads was negatively correlated (r = −0.73) with the seed sets. Thus, meiotic restitution in the 8n triticale does not contribute to fertility or increased ploidy in subsequent generations.

Comparative phylogenetics of repetitive elements in a diverse order of flowering plants (Brassicales)

Genome sizes of plants have long piqued the interest of researchers due to the vast differences among organisms. However, the mechanisms that drive size differences have yet to be fully understood. Two important contributing factors to genome size are expansions of repetitive elements, such as transposable elements (TEs), and whole-genome duplications (WGD). Although studies have found correlations between genome size and both TE abundance and polyploidy, these studies typically test for these patterns within a genus or species. The plant order Brassicales provides an excellent system to further test if genome size evolution patterns are consistent across larger time scales, as there are numerous WGDs. This order is also home to one of the smallest plant genomes, Arabidopsis thaliana-chosen as the model plant system for this reason-as well as to species with very large genomes. With new methods that allow for TE characterization from low-coverage genome shotgun data and 71 taxa across the Brassicales, we confirm the correlation between genome size and TE content, however, we are unable to reconstruct phylogenetic relationships and do not detect any shift in TE abundance associated with WGD.

Comparative Repeat Profiling of Two Closely Related Conifers (Larix decidua and Larix kaempferi) Reveals High Genome Similarity With Only Few Fast-Evolving Satellite DNAs

In eukaryotic genomes, cycles of repeat expansion and removal lead to large-scale genomic changes and propel organisms forward in evolution. However, in conifers, active repeat removal is thought to be limited, leading to expansions of their genomes, mostly exceeding 10 giga base pairs. As a result, conifer genomes are largely littered with fragmented and decayed repeats. Here, we aim to investigate how the repeat landscapes of two related conifers have diverged, given the conifers' accumulative genome evolution mode. For this, we applied low-coverage sequencing and read clustering to the genomes of European and Japanese larch, Larix decidua (Lamb.) Carrière and Larix kaempferi (Mill.), that arose from a common ancestor, but are now geographically isolated. We found that both Larix species harbored largely similar repeat landscapes, especially regarding the transposable element content. To pin down possible genomic changes, we focused on the repeat class with the fastest sequence turnover: satellite DNAs (satDNAs). Using comparative bioinformatics, Southern, and fluorescent in situ hybridization, we reveal the satDNAs' organizational patterns, their abundances, and chromosomal locations. Four out of the five identified satDNAs are widespread in the Larix genus, with two even present in the more distantly related Pseudotsuga and Abies genera. Unexpectedly, the EulaSat3 family was restricted to L. decidua and absent from L. kaempferi, indicating its evolutionarily young age. Taken together, our results exemplify how the accumulative genome evolution of conifers may limit the overall divergence of repeats after speciation, producing only few repeat-induced genomic novelties.

The maximum growth rate hypothesis is correct for eukaryotic photosynthetic organisms, but not cyanobacteria

The (maximum) growth rate (µ_max ) hypothesis predicts that cellular and tissue phosphorus (P) concentrations should increase with increasing growth rate, and RNA should also increase as most of the P is required to make ribosomes. Using published data, we show that though there is a strong positive relationship between the µ_max of all photosynthetic organisms and their P content (% dry weight), leading to a relatively constant P productivity, the relationship with RNA content is more complex. In eukaryotes there is a strong positive relationship between µ_max and RNA content expressed as % dry weight, and RNA constitutes a relatively constant 25% of total P. In prokaryotes the rRNA operon copy number is the important determinant of the amount of RNA present in the cell. The amount of phospholipid expressed as % dry weight increases with increasing µ_max in microalgae. The relative proportions of each of the five major P-containing constituents is remarkably constant, except that the proportion of RNA is greater and phospholipids smaller in prokaryotic than eukaryotic photosynthetic organisms. The effect of temperature differences between studies was minor. The evidence for and against P-containing constituents other than RNA being involved with ribosome synthesis and functioning is discussed.

Quantitative Approach to Fish Cytogenetics in the Context of Vertebrate Genome Evolution

Our novel Python-based tool EVANGELIST allows the visualization of GC and repeats percentages along chromosomes in sequenced genomes and has enabled us to perform quantitative large-scale analyses on the chromosome level in fish and other vertebrates. This is a different approach from the prevailing analyses, i.e., analyses of GC% in the coding sequences that make up not more than 2% in human. We identified GC content (GC%) elevations in microchromosomes in ancient fish lineages similar to avian microchromosomes and a large variability in the relationship between the chromosome size and their GC% across fish lineages. This raises the question as to what extent does the chromosome size drive GC% as posited by the currently accepted explanation based on the recombination rate. We ascribe the differences found across fishes to varying GC% of repetitive sequences. Generally, our results suggest that the GC% of repeats and proportion of repeats are independent of the chromosome size. This leaves an open space for another mechanism driving the GC evolution in vertebrates.

Effective double-digest RAD sequencing and genotyping despite large genome size

Obtaining informative data is the ambition of any genomic project, but in nonmodel species with very large genomes, pursuing such a goal requires surmounting a series of analytical challenges. Double-digest RAD sequencing is routinely used in nonmodel organisms and offers some control over the volume of data obtained. However, the volume of data recovered is not always an indication of the reliability of data sets, and quality checks are necessary to ensure that true and artefactual information is set apart. In the present study, we aim to fill the gap existing between the known applicability of RAD sequencing methods in plants with large genomes and the use of the retrieved loci for population genetic inference. By analysing two populations of Cypripedium calceolus, a nonmodel orchid species with a large genome size (1C ~ 31.6 Gbp), we provide a complete workflow from library preparation to bioinformatic filtering and inference of genetic diversity and differentiation. We show how filtering strategies to dismiss potentially misleading data need to be explored and adapted to data set-specific features. Moreover, we suggest that the occurrence of organellar sequences in libraries should not be neglected when planning the experiment and analysing the results. Finally, we explain how, in the absence of prior information about the genome of the species, seeking high standards of quality during library preparation and sequencing can provide an insurance against unpredicted technical or biological constraints.

The tight genome size of ants: diversity and evolution under ancestral state reconstruction and base composition

The mechanisms and processes driving change and variation in the genome size (GS) are not well known, and only a small set of ant species has been studied. Ants are an ecologically successful insect group present in most distinct ecosystems worldwide. Considering their wide distribution and ecological plasticity in different environmental contexts, we aimed to expand GS estimation within Formicidae to examine distribution patterns and variation in GS and base composition and to reconstruct the ancestral state of this character in an attempt to elucidate the generalized pattern of genomic expansions. Genome size estimates were generated for 99 ant species, including new GS estimates for 91 species of ants, and the mean GS of Formicidae was found to be 0.38 pg. The AT/GC ratio was 62.40/37.60. The phylogenetic reconstruction suggested an ancestral GS of 0.38 pg according to the Bayesian inference/Markov chain Monte Carlo method and 0.37 pg according to maximum likelihood and parsimony methods; significant differences in GS were observed between the subfamilies sampled. Our results suggest that the evolution of GS in Formicidae occurred through loss and accumulation of non-coding regions, mainly transposable elements, and occasionally by whole genome duplication. However, further studies are needed to verify whether these changes in DNA content are related to colonization processes, as suggested at the intraspecific level.

Genome Size Evolution and Dynamics in Iris, with Special Focus on the Section Oncocyclus

Insights into genome size dynamics and its evolutionary impact remain limited by the lack of data for many plant groups. One of these is the genus Iris, of which only 53 out of c. 260 species have available genome sizes. In this study, we estimated the C-values for 41 species and subspecies of Iris mainly from the Eastern Mediterranean region. We constructed a phylogenetic framework to shed light on the distribution of genome sizes across subgenera and sections of Iris. Finally, we tested evolutionary models to explore the mode and tempo of genome size evolution during the radiation of section Oncocyclus. Iris as a whole displayed a great variety of C-values; however, they were unequally distributed across the subgenera and sections, suggesting that lineage-specific patterns of genome size diversification have taken place within the genus. The evolutionary model that best fitted our data was the speciational model, as changes in genome size appeared to be mainly associated with speciation events. These results suggest that genome size dynamics may have contributed to the radiation of Oncocyclus irises. In addition, our phylogenetic analysis provided evidence that supports the segregation of the Lebanese population currently attributed to Iris persica as a distinct species.

DNA methylation enables transposable element-driven genome expansion

Multicellular eukaryotic genomes show enormous differences in size. A substantial part of this variation is due to the presence of transposable elements (TEs). They contribute significantly to a cell's mass of DNA and have the potential to become involved in host gene control. We argue that the suppression of their activities by methylation of the C-phosphate-G (CpG) dinucleotide in DNA is essential for their long-term accommodation in the host genome and, therefore, to its expansion. An inevitable consequence of cytosine methylation is an increase in C-to-T transition mutations via deamination, which causes CpG loss. Cytosine deamination is often needed for TEs to take on regulatory functions in the host genome. Our study of the whole-genome sequences of 53 organisms showed a positive correlation between the size of a genome and the percentage of TEs it contains, as well as a negative correlation between size and the CpG observed/expected (O/E) ratio in both TEs and the host DNA. TEs are seldom found at promoters and transcription start sites, but they are found more at enhancers, particularly after they have accumulated C-to-T and other mutations. Therefore, the methylation of TE DNA allows for genome expansion and also leads to new opportunities for gene control by TE-based regulatory sites.

Environmental pressures on stomatal size may drive plant genome size evolution: evidence from a natural experiment with Cape geophytes

BACKGROUND AND AIMS

The idea that genome (size) evolution in eukaryotes could be driven by environmental factors is still vigorously debated. In extant plants, genome size correlates positively with stomatal size, leading to the idea that conditions enabling the existence of large stomata in fossil plants also supported growth of their genome size. We test this inductive assumption in drought-adapted, prostrate-leaved Cape (South Africa) geophytes where, compared with their upright-leaved geophytic ancestors, stomata develop in a favourably humid microclimate formed underneath their leaves.

METHODS

Stomatal parameters (leaf cuticle imprints) and genome size (flow cytometry) were measured in 16 closely related geophytic species pairs from seven plant families. In each pair, representing a different genus, we contrasted a prostrate-leaved species with its upright-leaved phylogenetic relative, the latter whose stomata are exposed to the ambient arid climate.

KEY RESULTS

Except for one, all prostrate-leaves species had larger stomata, and in 13 of 16 pairs they also had larger genomes than their upright-leaved relatives. Stomatal density and theoretical maximum conductance were less in prostrate-leaved species with small guard cells (<1 pL) but showed no systematic difference in species pairs with larger guard cells (>1 pL). Giant stomata were observed in the prostrate-leaved Satyrium bicorne (89-137 µm long), despite its relatively small genome (2C = 9 Gbp).

CONCLUSIONS

Our results imply that climate, through selection on stomatal size, might be able to drive genome size evolution in plants. The data support the idea that plants from 'greenhouse' geological periods with large stomata might have generally had larger genome sizes when compared with extant plants, though this might not have been solely due to higher atmospheric CO2 in these periods but could also have been due to humid conditions prevailing at fossil deposit sites.

3 tera-basepairs as a fundamental limit for robust DNA replication

In order to maintain functional robustness and species integrity, organisms must ensure high fidelity of the genome duplication process. This is particularly true during early development, where cell division is often occurring both rapidly and coherently. By studying the extreme limits of suppressing DNA replication failure due to double fork stall errors, we uncover a fundamental constant that describes a trade-off between genome size and architectural complexity of the developing organism. This constant has the approximate value N U ≈ 3 × 1012 basepairs, and depends only on two highly conserved molecular properties of DNA biology. We show that our theory is successful in interpreting a diverse range of data across the Eukaryota.

Cross-Species Annotation of Expressed Genes and Detection of Different Functional Gene Modules Between 10 Cold- and 10 Hot-Propertied Chinese Herbal Medicines

According to the traditional Chinese medicine (TCM) system, Chinese herbal medicines (HMs) can be divided into four categories: hot, warm, cold, and cool. A cool nature usually is categorized as a cold nature, and a warm nature is classified as a hot nature. However, the detectable characteristics of the gene expression profile associated with the cold and hot properties have not been studied. To address this question, a strategy for the cross-species annotation of conserved genes was established in the present study by using transcriptome data of 20 HMs with cold and hot properties. Functional enrichment analysis was performed on group-specific expressed genes inferred from the functional genome of the reference species (i.e., Arabidopsis). Results showed that metabolic pathways relevant to chrysoeriol, luteolin, paniculatin, and wogonin were enriched for cold-specific genes, and pathways of inositol, heptadecane, lauric acid, octanoic acid, hexadecanoic acid, and pentadecanoic acid were enriched for hot-specific genes. Six functional modules were identified in the HMs with the cold property: nucleotide biosynthetic process, peptidy-L-cysteine S-palmitoylation, lipid modification, base-excision repair, dipeptide transport, and response to endoplasmic reticulum stress. For the hot HMs, another six functional modules were identified: embryonic meristem development, embryonic pattern specification, axis specification, regulation of RNA polymerase II transcriptional preinitiation complex assembly, mitochondrial RNA modification, and cell redox homeostasis. The research provided a new insight into HMs’ cold and hot properties from the perspective of the gene expression profile of plants.