Evolution of genome size and chromosome number in the carnivorous plant genus Genlisea (Lentibulariaceae), with a new estimate of the minimum genome size in angiosperms

Genome Size Evolution Differs Between Drosophila Subgenera with Striking Differences in Male and Female Genome Size in Sophophora

Genome size varies across the tree of life, with no clear correlation to organismal complexity or coding sequence, but with differences in non-coding regions. Phylogenetic methods have recently been incorporated to further disentangle this enigma, yet most of these studies have focused on widely diverged species. Few have compared patterns of genome size change in closely related species with known structural differences in the genome. As a consequence, the relationship between genome size and differences in chromosome number or inter-sexual differences attributed to XY systems are largely unstudied. We hypothesize that structural differences associated with chromosome number and X-Y chromosome differentiation, should result in differing rates and patterns of genome size change. In this study, we utilize the subgenera within the Drosophila to ask if patterns and rates of genome size change differ between closely related species with differences in chromosome numbers and states of the XY system. Genome sizes for males and females of 152 species are used to answer these questions (with 92 newly added or updated estimates). While we find no relationship between chromosome number and genome size or chromosome number and inter-sexual differences in genome size, we find evidence for differing patterns of genome size change between the subgenera, and increasing rates of change throughout time. Estimated shifts in rates of change in sex differences in genome size occur more often in Sophophora and correspond to known neo-sex events.

The genome size of clusterbean (Cyamopsis tetragonoloba) is significantly smaller compared to its wild relatives as estimated by flow cytometry

Clusterbean (C. tetragonoloba) is an important, leguminous vegetable and industrial crop with vast genetic diversity but meager genetic, cytological and genomic information. In the present study, an optimized procedure of flow cytometry was used to estimate the genome size of three clusterbean species, represented by C. tetragonoloba (cv. RGC-936) and two wild relatives (C. serreta and C. senegalensis). For accurate estimation of genomic content, singlet G₀/G₁ populations of multiple tissues such as leaves, hypocotyl, and matured seeds were determined and used along with three different plant species viz. Pisum sativum (as primary), Oryza sativa, and Glycine max (secondary), as external and internal reference standards. Seed tissue of the test sample and G. max provided the best estimate of nuclear DNA content in comparison to other sample tissues and reference standards. The genome size of C. tetragonoloba was detemined at 580.9±0.02Mbp (1C), while that of C. serreta and C. senegalensis was estimated at 979.6±0.02Mbp (1C) and 943.4±0.03Mbp (1C), respectively. Thus, the wild relatives harbor, nearly double the genome content of the cultivated cluster bean. Findings of this study will enrich genomic database of the legume family and can serve as the starting point for clusterbean evolutionary and genomics studies.

Genome survey and SSR analysis of Apocynum venetum

Apocynum venetum is an eco-economic plant that exhibits high stress resistance. In the present paper, we carried out a whole-genome survey of A. venetum in order to provide a foundation for its whole-genome sequencing. High-throughput sequencing technology (Illumina NovaSep) was first used to measure the genome size of A. venetum, and bioinformatics methods were employed for the evaluation of the genome size, heterozygosity ratio, repeated sequences, and GC content in order to provide a foundation for subsequent whole-genome sequencing. The sequencing analysis results indicated that the preliminary estimated genome size of A. venetum was 254.40 Mbp, and its heterozygosity ratio and percentage of repeated sequences were 0.63 and 40.87%, respectively, indicating that it has a complex genome. We used k-mer = 41 to carry out a preliminary assembly and obtained contig N50, which was 3841 bp with a total length of 223949699 bp. We carried out further assembly to obtain scaffold N50, which was 6196 bp with a total length of 227322054 bp. We performed simple sequence repeat (SSR) molecular marker prediction based on the A. venetum genome data and identified a total of 101918 SSRs. The differences between the different types of nucleotide repeats were large, with mononucleotide repeats being most numerous and hexanucleotide repeats being least numerous. We recommend the use of the ‘2+3’ (Illumina+PacBio) sequencing combination to supplement the Hi-C technique and resequencing technique in future whole-genome research in A. venetum.

PROBABILISTIC MODELING AND ANALYSIS OF DNA FRAGMENTATION

Deoxyribonucleic Acid (DNA) sequencing has become indispensable to the modern biological and medicine sciences. DNA fragmentation is a preliminary step in a dominant technique called shotgun sequencing that provides a time and cost effective strategy for the DNA sequencing. In this paper, we propose a probabilistic model for the random DNA fragmentation and derive an average number of fragments with the suitable length along with the probability of covering the entire DNA strand through the de novo assembly or the referenced-based mapping assembly. We formulate the coverage problem in terms of the probability of bond breaking between nucleotides and the number of DNA molecules participating in the fragmentation process, and provide insights into the optimal DNA fragmentation. We obtain the lower bound for the minimum number of suitable fragments required to reconstruct the DNA strand with the specified reliability. We evaluate the derived results with our DNA Fragmentation Tool which demonstrate, the validity of these results based on our model. Finally, we update our model with respect to the fragments’ size distribution of real data.

Ancient Polyploidy and Genome Evolution in Palms

Mechanisms of genome evolution are fundamental to our understanding of adaptation and the generation and maintenance of biodiversity, yet genome dynamics are still poorly characterized in many clades. Strong correlations between variation in genomic attributes and species diversity across the plant tree of life suggest that polyploidy or other mechanisms of genome size change confer selective advantages due to the introduction of genomic novelty. Palms (order Arecales, family Arecaceae) are diverse, widespread, and dominant in tropical ecosystems, yet little is known about genome evolution in this ecologically and economically important clade. Here, we take a phylogenetic comparative approach to investigate palm genome dynamics using genomic and transcriptomic data in combination with a recent, densely sampled, phylogenetic tree. We find conclusive evidence of a paleopolyploid event shared by the ancestor of palms but not with the sister clade, Dasypogonales. We find evidence of incremental chromosome number change in the palms as opposed to one of recurrent polyploidy. We find strong phylogenetic signal in chromosome number, but no signal in genome size, and further no correlation between the two when correcting for phylogenetic relationships. Palms thus add to a growing number of diverse, ecologically successful clades with evidence of whole-genome duplication, sister to a species-poor clade with no evidence of such an event. Disentangling the causes of genome size variation in palms moves us closer to understanding the genomic conditions facilitating adaptive radiation and ecological dominance in an evolutionarily successful, emblematic tropical clade.

Plant genome sequences: past, present, future

The green plants (Viridiplantae) are an essential kingdom of life, responsible via photosynthesis for the majority of global primary production, and directly utilized by humankind for nutrition, animal feed, fuel, clothing, medicine and other purposes. There are an estimated 391 000 species of land plants, in addition to 8000 species of green algae. Their genomes are unusually diverse compared to those of other kingdoms, ranging in size from ∼10 Mb to over 100 Gb. Knowledge of plant genomes initially lagged behind those of other kingdoms but has greatly increased with the development of new technologies for DNA sequencing; bioinformatic analysis, rather than data production, is increasingly the bottleneck to further knowledge. Recent proposals are now contemplating the sequencing, assembly and annotation of the genomes of all of the world's plant species; meanwhile, low coverage sequencing to measure diversity across collections and wild populations has already become commonplace for many species, especially those utilized as crops.

A critical comparison of technologies for a plant genome sequencing project

BACKGROUND

A high-quality genome sequence of any model organism is an essential starting point for genetic and other studies. Older clone-based methods are slow and expensive, whereas faster, cheaper short-read-only assemblies can be incomplete and highly fragmented, which minimizes their usefulness. The last few years have seen the introduction of many new technologies for genome assembly. These new technologies and associated new algorithms are typically benchmarked on microbial genomes or, if they scale appropriately, on larger (e.g., human) genomes. However, plant genomes can be much more repetitive and larger than the human genome, and plant biochemistry often makes obtaining high-quality DNA that is free from contaminants difficult. Reflecting their challenging nature, we observe that plant genome assembly statistics are typically poorer than for vertebrates.

RESULTS

Here, we compare Illumina short read, Pacific Biosciences long read, 10x Genomics linked reads, Dovetail Hi-C, and BioNano Genomics optical maps, singly and combined, in producing high-quality long-range genome assemblies of the potato species Solanum verrucosum. We benchmark the assemblies for completeness and accuracy, as well as DNA compute requirements and sequencing costs.

CONCLUSIONS

The field of genome sequencing and assembly is reaching maturity, and the differences we observe between assemblies are surprisingly small. We expect that our results will be helpful to other genome projects, and that these datasets will be used in benchmarking by assembly algorithm developers.

Plastome-Wide Rearrangements and Gene Losses in Carnivorous Droseraceae

The plastid genomes of four related carnivorous plants (Drosera regia, Drosera erythrorhiza, Aldrovanda vesiculosa, and Dionaea muscipula) were sequenced to examine changes potentially induced by the transition to carnivory. The plastid genomes of the Droseraceae show multiple rearrangements, gene losses, and large expansions or contractions of the inverted repeat. All the ndh genes are lost or nonfunctional, as well as in some of the species, clpP1, ycf1, ycf2 and some tRNA genes. Uniquely, among land plants, the trnK gene has no intron. Carnivory in the Droseraceae coincides with changes in plastid gene content similar to those induced by parasitism and mycoheterotrophy, suggesting parallel changes in chloroplast function due to the similar switch from autotrophy to (mixo-) heterotrophy. A molecular phylogeny of the taxa based on all shared plastid genes indicates that the "snap-traps" of Aldrovanda and Dionaea have a common origin.

Progress in the study of genome size evolution in Asteraceae: analysis of the last update

The Genome Size in Asteraceae Database (GSAD, http://www.asteraceaegenomesize.com) has been recently updated, with data from papers published or in press until July 2018. This constitutes the third release of GSAD, currently containing 4350 data entries for 1496 species, which represent a growth of 22.52% in the number of species with available genome size data compared with the previous release, and a growth of 57.72% in terms of entries. Approximately 6% of Asteraceae species are covered in terms of known genome sizes. The number of source papers included in this release (198) means a 48.87% increase with respect to release 2.0. The significant data increase was exploited to study the genome size evolution in the family from a phylogenetic perspective. Our results suggest that the role of chromosome number in genome size diversity within Asteraceae is basically associated to polyploidy, while dysploidy would only cause minor variation in the DNA amount along the family. Among diploid taxa, we found that the evolution of genome size shows a strong phylogenetic signal. However, this trait does not seem to evolve evenly across the phylogeny, but there could be significant scale and clade-dependent patterns. Our analyses indicate that the phylogenetic signal is stronger at low taxonomic levels, with certain tribes standing out as hotspots of autocorrelation between genome size and phylogeny. Finally, we also observe meaningful associations among nuclear DNA content on Asteraceae species and other phenotypical and ecological traits (i.e. plant habit and invasion ability). Overall, this study emphasizes the need to continue generating and analysing genome size data in order to puzzle out the evolution of this parameter and its many biological correlates.

On the Road to Breeding 4.0: Unraveling the Good, the Bad, and the Boring of Crop Quantitative Genomics

Understanding the quantitative genetics of crops has been and will continue to be central to maintaining and improving global food security. We outline four stages that plant breeding either has already achieved or will probably soon achieve. Top-of-the-line breeding programs are currently in Breeding 3.0, where inexpensive, genome-wide data coupled with powerful algorithms allow us to start breeding on predicted instead of measured phenotypes. We focus on three major questions that must be answered to move from current Breeding 3.0 practices to Breeding 4.0: ( a) How do we adapt crops to better fit agricultural environments? ( b) What is the nature of the diversity upon which breeding can act? ( c) How do we deal with deleterious variants? Answering these questions and then translating them to actual gains for farmers will be a significant part of achieving global food security in the twenty-first century.

A Mechanism for Genome Size Reduction Following Genomic Rearrangements

The factors behind genome size evolution have been of great interest, considering that eukaryotic genomes vary in size by more than three orders of magnitude. Using a model of two wild peanut relatives, Arachis duranensis and Arachis ipaensis, in which one genome experienced large rearrangements, we find that the main determinant in genome size reduction is a set of inversions that occurred in A. duranensis, and subsequent net sequence removal in the inverted regions. We observe a general pattern in which sequence is lost more rapidly at newly distal (telomeric) regions than it is gained at newly proximal (pericentromeric) regions - resulting in net sequence loss in the inverted regions. The major driver of this process is recombination, determined by the chromosomal location. Any type of genomic rearrangement that exposes proximal regions to higher recombination rates can cause genome size reduction by this mechanism. In comparisons between A. duranensis and A. ipaensis, we find that the inversions all occurred in A. duranensis. Sequence loss in those regions was primarily due to removal of transposable elements. Illegitimate recombination is likely the major mechanism responsible for the sequence removal, rather than unequal intrastrand recombination. We also measure the relative rate of genome size reduction in these two Arachis diploids. We also test our model in other plant species and find that it applies in all cases examined, suggesting our model is widely applicable.

Genome Size Diversity and Its Impact on the Evolution of Land Plants

Genome size is a biodiversity trait that shows staggering diversity across eukaryotes, varying over 64,000-fold. Of all major taxonomic groups, land plants stand out due to their staggering genome size diversity, ranging ca. 2400-fold. As our understanding of the implications and significance of this remarkable genome size diversity in land plants grows, it is becoming increasingly evident that this trait plays not only an important role in shaping the evolution of plant genomes, but also in influencing plant community assemblages at the ecosystem level. Recent advances and improvements in novel sequencing technologies, as well as analytical tools, make it possible to gain critical insights into the genomic and epigenetic mechanisms underpinning genome size changes. In this review we provide an overview of our current understanding of genome size diversity across the different land plant groups, its implications on the biology of the genome and what future directions need to be addressed to fill key knowledge gaps.

Comparative genomic analysis of Genlisea (corkscrew plants-Lentibulariaceae) chloroplast genomes reveals an increasing loss of the ndh genes

In the carnivorous plant family Lentibulariaceae, all three genome compartments (nuclear, chloroplast, and mitochondria) have some of the highest rates of nucleotide substitutions across angiosperms. While the genera Genlisea and Utricularia have the smallest known flowering plant nuclear genomes, the chloroplast genomes (cpDNA) are mostly structurally conserved except for deletion and/or pseudogenization of the NAD(P)H-dehydrogenase complex (ndh) genes known to be involved in stress conditions of low light or CO2 concentrations. In order to determine how the cpDNA are changing, and to better understand the evolutionary history within the Genlisea genus, we sequenced, assembled and analyzed complete cpDNA from six species (G. aurea, G. filiformis, G. pygmaea, G. repens, G. tuberosa and G. violacea) together with the publicly available G. margaretae cpDNA. In general, the cpDNA structure among the analyzed Genlisea species is highly similar. However, we found that the plastidial ndh genes underwent a progressive process of degradation similar to the other terrestrial Lentibulariaceae cpDNA analyzed to date, but in contrast to the aquatic species. Contrary to current thinking that the terrestrial environment is a more stressful environment and thus requiring the ndh genes, we provide evidence that in the Lentibulariaceae the terrestrial forms have progressive loss while the aquatic forms have the eleven plastidial ndh genes intact. Therefore, the Lentibulariaceae system provides an important opportunity to understand the evolutionary forces that govern the transition to an aquatic environment and may provide insight into how plants manage water stress at a genome scale.

Development of microsatellite markers for the carnivorous plant Genlisea aurea (Lentibulariaceae) using genomics data of NGS

Genlisea aurea A.St.-Hil. is a carnivorous plant endemic species to Brazil in the Lentibulariaceae family. Very few studies have addressed the genetic structure and conservation status of G. aurea and the Lentibulariaceae. Microsatellites markers are advantageous tools that can be employed to predict the vulnerability of Lentibulariaceae species. Therefore, the development of molecular markers focusing the population analyses of Genlisea for future genetic studies and conservation actions are essential. Thus, we developed simple sequence repeats (SSRs) based on in silico analyses of G. aurea draft genome assembly. We characterized 40 individuals from several populations and identified 12 loci that were polymorphic, with heterozygosity between 0.123 and 0.650. We demonstrated that the G. aurea SSR markers work cross-species in Genlisea filiformis, G. repens, G. tuberosa and G. violacea. These markers will be important for future population, phylogeographic and conservation studies in G. aurea and other Genlisea species.

Angiosperm Plant Desiccation Tolerance: Hints from Transcriptomics and Genome Sequencing

Desiccation tolerance (DT) in angiosperms is present in the small group of resurrection plants and in seeds. DT requires the presence of protective proteins, specific carbohydrates, restructuring of membrane lipids, and regulatory mechanisms directing a dedicated gene expression program. Many components are common to resurrection plants and seeds; however, some are specific for resurrection plants. Understanding how each component contributes to DT is challenging. Recent transcriptome analyses and genome sequencing indicate that increased expression is essential of genes encoding protective components, recently evolved, species-specific genes and non-protein-coding RNAs. Modification and reshuffling of existing cis-regulatory promoter elements seems to play a role in the rewiring of regulatory networks required for increased expression of DT-related genes in resurrection species.

Inferring the evolutionary mechanism of the chloroplast genome size by comparing whole-chloroplast genome sequences in seed plants

The chloroplast genome originated from photosynthetic organisms and has retained the core genes that mainly encode components of photosynthesis. However, the causes of variations in chloroplast genome size in seed plants have only been thoroughly analyzed within small subsets of spermatophytes. In this study, we conducted the first comparative analysis on a large scale to examine the relationship between sequence characteristics and genome size in 272 seed plants based on cross-species and phylogenetic signal analysis. Our results showed that inverted repeat regions, large or small single copies, intergenic regions, and gene number can be attributed to the variations in chloroplast genome size among closely related species. However, chloroplast gene length underwent evolution affecting chloroplast genome size in seed plants irrespective of whether phylogenetic information was incorporated. Among chloroplast genes, atpA, accD and ycf1 account for 13% of the variation in genome size, and the average Ka/Ks values of homologous pairs of the three genes are larger than 1. The relationship between chloroplast genome size and gene length might be affected by selection during the evolution of spermatophytes. The variation in chloroplast genome size may influence energy generation and ecological strategy in seed plants.

Evolution of genome size and genomic GC content in carnivorous holokinetics (Droseraceae)

BACKGROUND AND AIMS

Studies in the carnivorous family Lentibulariaceae in the last years resulted in the discovery of the smallest plant genomes and an unusual pattern of genomic GC content evolution. However, scarcity of genomic data in other carnivorous clades still prevents a generalization of the observed patterns. Here the aim was to fill this gap by mapping genome evolution in the second largest carnivorous family, Droseraceae, where this evolution may be affected by chromosomal holokinetism in Drosera METHODS: The genome size and genomic GC content of 71 Droseraceae species were measured by flow cytometry. A dated phylogeny was constructed, and the evolution of both genomic parameters and their relationship to species climatic niches were tested using phylogeny-based statistics.

KEY RESULTS

The 2C genome size of Droseraceae varied between 488 and 10 927 Mbp, and the GC content ranged between 37·1 and 44·7 %. The genome sizes and genomic GC content of carnivorous and holocentric species did not differ from those of their non-carnivorous and monocentric relatives. The genomic GC content positively correlated with genome size and annual temperature fluctuations. The genome size and chromosome numbers were inversely correlated in the Australian clade of Drosera CONCLUSIONS: Our results indicate that neither carnivory (nutrient scarcity) nor the holokinetism have a prominent effect on size and DNA base composition of Droseraceae genomes. However, the holokinetic drive seems to affect karyotype evolution in one of the major clades of Drosera Our survey confirmed that the evolution of GC content is tightly connected with the evolution of genome size and also with environmental conditions.

The Chloroplast Genome of Utricularia reniformis Sheds Light on the Evolution of the ndh Gene Complex of Terrestrial Carnivorous Plants from the Lentibulariaceae Family

Lentibulariaceae is the richest family of carnivorous plants spanning three genera including Pinguicula, Genlisea, and Utricularia. Utricularia is globally distributed, and, unlike Pinguicula and Genlisea, has both aquatic and terrestrial forms. In this study we present the analysis of the chloroplast (cp) genome of the terrestrial Utricularia reniformis. U. reniformis has a standard cp genome of 139,725bp, encoding a gene repertoire similar to essentially all photosynthetic organisms. However, an exclusive combination of losses and pseudogenization of the plastid NAD(P)H-dehydrogenase (ndh) gene complex were observed. Comparisons among aquatic and terrestrial forms of Pinguicula, Genlisea, and Utricularia indicate that, whereas the aquatic forms retained functional copies of the eleven ndh genes, these have been lost or truncated in terrestrial forms, suggesting that the ndh function may be dispensable in terrestrial Lentibulariaceae. Phylogenetic scenarios of the ndh gene loss and recovery among Pinguicula, Genlisea, and Utricularia to the ancestral Lentibulariaceae cladeare proposed. Interestingly, RNAseq analysis evidenced that U. reniformis cp genes are transcribed, including the truncated ndh genes, suggesting that these are not completely inactivated. In addition, potential novel RNA-editing sites were identified in at least six U. reniformis cp genes, while none were identified in the truncated ndh genes. Moreover, phylogenomic analyses support that Lentibulariaceae is monophyletic, belonging to the higher core Lamiales clade, corroborating the hypothesis that the first Utricularia lineage emerged in terrestrial habitats and then evolved to epiphytic and aquatic forms. Furthermore, several truncated cp genes were found interspersed with U. reniformis mitochondrial and nuclear genome scaffolds, indicating that as observed in other smaller plant genomes, such as Arabidopsis thaliana, and the related and carnivorous Genlisea nigrocaulis and G. hispidula, the endosymbiotic gene transfer may also shape the U. reniformis genome in a similar fashion. Overall the comparative analysis of the U. reniformis cp genome provides new insight into the ndh genes and cp genome evolution of carnivorous plants from Lentibulariaceae family.

The effects of fresh and rapid desiccated tissue on estimates of Ophiopogoneae genome size

Fresh plant material is usually used for genome size estimation by flow cytometry (FCM). Lack of fresh material is cited as one of the main reasons for the dearth of studies on plants from remote locations. Genome sizes in fresh versus desiccated tissue of 16 Ophiopogoneae species and five model plant species were estimated. Our results indicated that desiccated tissue was suitable for genome size estimation; this method enables broader geographic sampling of plants when fresh tissue collection is not feasible. To be useful, after dessication the Ophiopogoneae sample should be green without brown or yellow markings; it should be stored in deep freezer at -80 °C, and the storage time should be no more than 6 months.

The Small Nuclear Genomes of Selaginella Are Associated with a Low Rate of Genome Size Evolution

The haploid nuclear genome size (1C DNA) of vascular land plants varies over several orders of magnitude. Much of this observed diversity in genome size is due to the proliferation and deletion of transposable elements. To date, all vascular land plant lineages with extremely small nuclear genomes represent recently derived states, having ancestors with much larger genome sizes. The Selaginellaceae represent an ancient lineage with extremely small genomes. It is unclear how small nuclear genomes evolved in Selaginella We compared the rates of nuclear genome size evolution in Selaginella and major vascular plant clades in a comparative phylogenetic framework. For the analyses, we collected 29 new flow cytometry estimates of haploid genome size in Selaginella to augment publicly available data. Selaginella possess some of the smallest known haploid nuclear genome sizes, as well as the lowest rate of genome size evolution observed across all vascular land plants included in our analyses. Additionally, our analyses provide strong support for a history of haploid nuclear genome size stasis in Selaginella Our results indicate that Selaginella, similar to other early diverging lineages of vascular land plants, has relatively low rates of genome size evolution. Further, our analyses highlight that a rapid transition to a small genome size is only one route to an extremely small genome.

Genome size and ploidy influence angiosperm species' biomass under nitrogen and phosphorus limitation

Angiosperm genome sizes (GS) range c. 2400-fold, and as nucleic acids are amongst the most phosphorus- (P) and nitrogen (N)-demanding cellular biomolecules, we test the hypothesis that a key influence on plant biomass and species composition is the interaction between N and P availability and plant GS. We analysed the impact of different nutrient regimes on above-ground biomass of angiosperm species with different GS, ploidy level and Grime's C-S-R (competitive, stress-tolerant, ruderal) plant strategies growing at the Park Grass Experiment (Rothamsted, UK), established in 1856. The biomass-weighted mean GS of species growing on plots with the addition of both N and P fertilizer were significantly higher than that of plants growing on control plots and plots with either N or P. The plants on these N + P plots are dominated by polyploids with large GS and a competitive plant strategy. The results are consistent with our hypothesis that large genomes are costly to build and maintain under N and P limitation. Hence GS and ploidy are significant traits affecting biomass growth under different nutrient regimes, influencing plant community composition and ecosystem dynamics. We propose that GS is a critical factor needed in models that bridge the knowledge gap between biodiversity and ecosystem functioning.

Chromosome identification for the carnivorous plant Genlisea margaretae

Genlisea margaretae, subgenus Genlisea, section Recurvatae (184 Mbp/1C), belongs to a plant genus with a 25-fold genome size difference and an extreme genome plasticity. Its 19 chromosome pairs could be distinguished individually by an approach combining optimized probe pooling and consecutive rounds of multicolor fluorescence in situ hybridization (mcFISH) with bacterial artificial chromosomes (BACs) selected for repeat-free inserts. Fifty-one BACs were assigned to 18 chromosome pairs. They provide a tool for future assignment of genomic sequence contigs to distinct chromosomes as well as for identification of homeologous chromosome regions in other species of the carnivorous Lentibulariaceae family, and potentially of chromosome rearrangements, in cases where more than one BAC per chromosome pair was identified.

Venus flytrap carnivorous lifestyle builds on herbivore defense strategies

Although the concept of botanical carnivory has been known since Darwin's time, the molecular mechanisms that allow animal feeding remain unknown, primarily due to a complete lack of genomic information. Here, we show that the transcriptomic landscape of the Dionaea trap is dramatically shifted toward signal transduction and nutrient transport upon insect feeding, with touch hormone signaling and protein secretion prevailing. At the same time, a massive induction of general defense responses is accompanied by the repression of cell death-related genes/processes. We hypothesize that the carnivory syndrome of Dionaea evolved by exaptation of ancient defense pathways, replacing cell death with nutrient acquisition.

Cenh3: An Emerging Player in Haploid Induction Technology

True-breeding lines are required for the development and production of crop varieties. In a classical breeding approach these lines are obtained through inbreeding, and often 7-9 generations of inbreeding is performed to achieve the desired level of homozygosity, over a period of several years. In contrast, the chromosomes of haploids can be doubled to produce true-breeding lines in a single generation. Over the last century, scientists have developed a variety of techniques to induce haploids and doubled haploids, though these techniques apply only to particular crop varieties. Ravi and Chan (2010) discovered that haploids could be obtained in Arabidopsis through the manipulation of the centromere-specific histone 3 variant, CENH3. Their approach, which involved extensive modifications to a transgenic CENH3, held promise of being translated to crop species, and has been successfully employed in maize (see Kelliher et al., 2016). Refinements of this technology have since been developed which indicate that non-transgenic modifications to CENH3 will also induce haploids. The complementation of a cenh3 null by CENH3 from closely related plant species can result in plants that are fertile but haploid-inducing on crossing by CENH3 wt plants- suggesting that introgression of alien CENH3 may produce non-transgenic haploid inducers. Similarly, a remarkably wide variety of point mutations in CENH3, inducible by chemical agents, have recently been shown to result in haploid induction on crossing by wild-type CENH3 plants. These CENH3-variant plants grow normally, are fully fertile on self-pollination, and may be present in existing mutagenized collections.

Centromere and telomere sequence alterations reflect the rapid genome evolution within the carnivorous plant genus Genlisea

Linear chromosomes of eukaryotic organisms invariably possess centromeres and telomeres to ensure proper chromosome segregation during nuclear divisions and to protect the chromosome ends from deterioration and fusion, respectively. While centromeric sequences may differ between species, with arrays of tandemly repeated sequences and retrotransposons being the most abundant sequence types in plant centromeres, telomeric sequences are usually highly conserved among plants and other organisms. The genome size of the carnivorous genus Genlisea (Lentibulariaceae) is highly variable. Here we study evolutionary sequence plasticity of these chromosomal domains at an intrageneric level. We show that Genlisea nigrocaulis (1C = 86 Mbp; 2n = 40) and G. hispidula (1C = 1550 Mbp; 2n = 40) differ as to their DNA composition at centromeres and telomeres. G. nigrocaulis and its close relative G. pygmaea revealed mainly 161 bp tandem repeats, while G. hispidula and its close relative G. subglabra displayed a combination of four retroelements at centromeric positions. G. nigrocaulis and G. pygmaea chromosome ends are characterized by the Arabidopsis-type telomeric repeats (TTTAGGG); G. hispidula and G. subglabra instead revealed two intermingled sequence variants (TTCAGG and TTTCAGG). These differences in centromeric and, surprisingly, also in telomeric DNA sequences, uncovered between groups with on average a > 9-fold genome size difference, emphasize the fast genome evolution within this genus. Such intrageneric evolutionary alteration of telomeric repeats with cytosine in the guanine-rich strand, not yet known for plants, might impact the epigenetic telomere chromatin modification.

In Depth Characterization of Repetitive DNA in 23 Plant Genomes Reveals Sources of Genome Size Variation in the Legume Tribe Fabeae

The differential accumulation and elimination of repetitive DNA are key drivers of genome size variation in flowering plants, yet there have been few studies which have analysed how different types of repeats in related species contribute to genome size evolution within a phylogenetic context. This question is addressed here by conducting large-scale comparative analysis of repeats in 23 species from four genera of the monophyletic legume tribe Fabeae, representing a 7.6-fold variation in genome size. Phylogenetic analysis and genome size reconstruction revealed that this diversity arose from genome size expansions and contractions in different lineages during the evolution of Fabeae. Employing a combination of low-pass genome sequencing with novel bioinformatic approaches resulted in identification and quantification of repeats making up 55–83% of the investigated genomes. In turn, this enabled an analysis of how each major repeat type contributed to the genome size variation encountered. Differential accumulation of repetitive DNA was found to account for 85% of the genome size differences between the species, and most (57%) of this variation was found to be driven by a single lineage of Ty3/gypsy LTR-retrotransposons, the Ogre elements. Although the amounts of several other lineages of LTR-retrotransposons and the total amount of satellite DNA were also positively correlated with genome size, their contributions to genome size variation were much smaller (up to 6%). Repeat analysis within a phylogenetic framework also revealed profound differences in the extent of sequence conservation between different repeat types across Fabeae. In addition to these findings, the study has provided a proof of concept for the approach combining recent developments in sequencing and bioinformatics to perform comparative analyses of repetitive DNAs in a large number of non-model species without the need to assemble their genomes.

Fastest predators in the plant kingdom: functional morphology and biomechanics of suction traps found in the largest genus of carnivorous plants

Understanding the physics of plant movements, which describe the interplay between plant architecture, movement speed and actuation principles, is essential for the comprehension of important processes like plant morphogenesis. Recent investigations especially on rapid plant movements at the interface of biology, physics and engineering sciences highlight how such fast motions can be achieved without the presence of muscles, nerves and technical hinge analogies. The suction traps (bladders) of carnivorous bladderworts (Utricularia spp., Lentibulariaceae, Lamiales) are considered as some of the most elaborate moving structures in the plant kingdom. A complex interplay of morphological and physiological adaptations allows the traps to pump water out of their body and to store elastic energy in the deformed bladder walls. Mechanical stimulation by prey entails opening of the otherwise watertight trapdoor, followed by trap wall relaxation, sucking in of water and prey, and consecutive trapdoor closure. Suction can also occur spontaneously in non-stimulated traps. We review the current state of knowledge about the suction trap mechanism with a focus on architectonically homogeneous traps of aquatic bladderwort species from section Utricularia (the so-called 'Utricularia vulgaris trap type'). The functional morphology and biomechanics of the traps are described in detail. We discuss open questions and propose promising aspects for future studies on these sophisticated ultra-fast trapping devices.

Comparative Genome Analysis Reveals Divergent Genome Size Evolution in a Carnivorous Plant Genus

The C-value paradox remains incompletely resolved after >40 yr and is exemplified by 2,350-fold variation in genome sizes of flowering plants. The carnivorous Lentibulariaceae genus Genlisea, displaying a 25-fold range of genome sizes, is a promising subject to study mechanisms and consequences of evolutionary genome size variation. Applying genomic, phylogenetic, and cytogenetic approaches, we uncovered bidirectional genome size evolution within the genus Genlisea. The Genlisea nigrocaulis Steyerm. genome (86 Mbp) has probably shrunk by retroelement silencing and deletion-biased double-strand break (DSB) repair, from an ancestral size of 400 to 800 Mbp to become one of the smallest among flowering plants. The G. hispidula Stapf genome has expanded by whole-genome duplication (WGD) and retrotransposition to 1550 Mbp. Genlisea hispidula became allotetraploid after the split from the G. nigrocaulis clade ∼29 Ma. Genlisea pygmaea A. St.-Hil. (179 Mbp), a close relative of G. nigrocaulis, proved to be a recent (auto)tetraploid. Our analyses suggest a common ancestor of the genus Genlisea with an intermediate 1C value (400-800 Mbp) and subsequent rapid genome size evolution in opposite directions. Many abundant repeats of the larger genome are absent in the smaller, casting doubt on their functionality for the organism, while recurrent WGD seems to safeguard against the loss of essential elements in the face of genome shrinkage. We cannot identify any consistent differences in habitat or life strategy that correlate with genome size changes, raising the possibility that these changes may be selectively neutral.

The blue butterfly Polyommatus (Plebicula) atlanticus (Lepidoptera, Lycaenidae) holds the record of the highest number of chromosomes in the non-polyploid eukaryotic organisms

The blue butterfly species Polyommatus (Plebicula) atlanticus (Elwes, 1906) (Lepidoptera, Lycaenidae) is known to have a very high haploid number of chromosomes (n= circa 223). However, this approximate count made by Hugo de Lesse 45 years ago was based on analysis of a single meiotic I metaphase plate, not confirmed by study of diploid chromosome set and not documented by microphotographs. Here I demonstrate that (1) Polyommatus atlanticus is a diploid (non-polyploid) species, (2) its meiotic I chromosome complement includes at least 224-226 countable chromosome bodies, and (3) all (or nearly all) chromosome elements in meiotic I karyotype are represented by bivalents. I also provide the first data on the diploid karyotype and estimate the diploid chromosome number as 2n=ca448-452. Thus, Polyommatus atlanticus is confirmed to possess the highest chromosome number among all the non-polyploid eukaryotic organisms.

Chromatin organization and cytological features of carnivorous Genlisea species with large genome size differences

The monophyletic carnivorous genus Genlisea (Lentibulariaceae) is characterized by a bi-directional genome size evolution resulting in a 25-fold difference in nuclear DNA content. This is one of the largest ranges found within a genus so far and makes Genlisea an interesting subject to study mechanisms of genome and karyotype evolution. Genlisea nigrocaulis, with 86 Mbp one of the smallest plant genomes, and the 18-fold larger genome of G. hispidula (1,550 Mbp) possess identical chromosome numbers (2n = 40) but differ considerably in chromatin organization, nuclear and cell size. Interphase nuclei of G. nigrocaulis and of related species with small genomes, G. aurea (133 Mbp, 2n ≈ 104) and G. pygmaea (179 Mbp, 2n = 80), are hallmarked by intensely DAPI-stained chromocenters, carrying typical heterochromatin-associated methylation marks (5-methylcytosine, H3K9me2), while in G. hispidula and surprisingly also in the small genome of G. margaretae (184 Mbp, 2n = 38) the heterochromatin marks are more evenly distributed. Probes of tandem repetitive sequences together with rDNA allow the unequivocal discrimination of 13 out of 20 chromosome pairs of G. hispidula. One of the repetitive sequences labeled half of the chromosome set almost homogenously supporting an allopolyploid status of G. hispidula and its close relative G. subglabra (1,622 Mbp, 2n = 40). In G. nigrocaulis 11 chromosome pairs could be individualized using a combination of rDNA and unique genomic probes. The presented data provide a basis for future studies of karyotype evolution within the genus Genlisea.

Quantitative testing of the methodology for genome size estimation in plants using flow cytometry: a case study of the Primulina genus

Flow cytometry (FCM) is a commonly used method for estimating genome size in many organisms. The use of FCM in plants is influenced by endogenous fluorescence inhibitors and may cause an inaccurate estimation of genome size; thus, falsifying the relationship between genome size and phenotypic traits/ecological performance. Quantitative optimization of FCM methodology minimizes such errors, yet there are few studies detailing this methodology. We selected the genus Primulina, one of the most representative and diverse genera of the Old World Gesneriaceae, to evaluate the methodology effect on determining genome size. Our results showed that buffer choice significantly affected genome size estimation in six out of the eight species examined and altered the 2C-value (DNA content) by as much as 21.4%. The staining duration and propidium iodide (PI) concentration slightly affected the 2C-value. Our experiments showed better histogram quality when the samples were stained for 40 min at a PI concentration of 100 μg ml(-1). The quality of the estimates was not improved by 1-day incubation in the dark at 4°C or by centrifugation. Thus, our study determined an optimum protocol for genome size measurement in Primulina: LB01 buffer supplemented with 100 μg ml(-1) PI and stained for 40 min. This protocol also demonstrated a high universality in other Gesneriaceae genera. We report the genome size of nine Gesneriaceae species for the first time. The results showed substantial genome size variation both within and among the species, with the 2C-value ranging between 1.62 and 2.71 pg. Our study highlights the necessity of optimizing the FCM methodology prior to obtaining reliable genome size estimates in a given taxon.

Progress, challenges and the future of crop genomes

The availability of plant reference genomes has ushered in a new era of crop genomics. More than 100 plant genomes have been sequenced since 2000, 63% of which are crop species. These genome sequences provide insight into architecture, evolution and novel aspects of crop genomes such as the retention of key agronomic traits after whole genome duplication events. Some crops have very large, polyploid, repeat-rich genomes, which require innovative strategies for sequencing, assembly and analysis. Even low quality reference genomes have the potential to improve crop germplasm through genome-wide molecular markers, which decrease expensive phenotyping and breeding cycles. The next stage of plant genomics will require draft genome refinement, building resources for crop wild relatives, resequencing broad diversity panels, and plant ENCODE projects to better understand the complexities of these highly diverse genomes.