The quest for suitable reference standards in genome size research

Reference standards for flow cytometric estimation of absolute nuclear DNA content in plants

The estimation of nuclear DNA content has been by far the most popular application of flow cytometry in plants. Because flow cytometry measures relative fluorescence intensities of nuclei stained by a DNA fluorochrome, ploidy determination, and estimation of the nuclear DNA content in absolute units both require comparison to a reference standard of known DNA content. This implies that the quality of the results obtained depends on the standard selection and use. Internal standardization, when the nuclei of an unknown sample and the reference standard are isolated, stained, and measured simultaneously, is mandatory for precise measurements. As DNA peaks representing G1 /G0 nuclei of the sample and standard appear on the same histogram of fluorescence intensity, the quotient of their position on the fluorescence intensity axis provides the quotient of DNA amounts. For the estimation of DNA amounts in absolute units, a number of well-established standards are now available to cover the range of known plant genome sizes. Since there are different standards in use, the standard and the genome size assigned to it has always to be reported. When none of the established standards fits, the introduction of a new standard species is needed. For this purpose, the regression line approach or simultaneous analysis of the candidate standard with several established standards should be prioritized. Moreover, the newly selected standard organism has to fulfill a number of requirements: it should be easy to identify and maintain, taxonomically unambiguous, globally available, with known genome size stability, lacking problematic metabolites, suitable for isolation of sufficient amounts of nuclei, and enabling measurements with low coefficients of variation of DNA peaks, hence suitable for the preparation of high quality samples.

Measuring the Invisible: The Sequences Causal of Genome Size Differences in Eyebrights (Euphrasia) Revealed by k-mers

Genome size variation within plant taxa is due to presence/absence variation, which may affect low-copy sequences or genomic repeats of various frequency classes. However, identifying the sequences underpinning genome size variation is challenging because genome assemblies commonly contain collapsed representations of repetitive sequences and because genome skimming studies by design miss low-copy number sequences. Here, we take a novel approach based on k-mers, short sub-sequences of equal length k, generated from whole-genome sequencing data of diploid eyebrights (Euphrasia), a group of plants that have considerable genome size variation within a ploidy level. We compare k-mer inventories within and between closely related species, and quantify the contribution of different copy number classes to genome size differences. We further match high-copy number k-mers to specific repeat types as retrieved from the RepeatExplorer2 pipeline. We find genome size differences of up to 230Mbp, equivalent to more than 20% genome size variation. The largest contributions to these differences come from rDNA sequences, a 145-nt genomic satellite and a repeat associated with an Angela transposable element. We also find size differences in the low-copy number class (copy number ≤ 10×) of up to 27 Mbp, possibly indicating differences in gene space between our samples. We demonstrate that it is possible to pinpoint the sequences causing genome size variation within species without the use of a reference genome. Such sequences can serve as targets for future cytogenetic studies. We also show that studies of genome size variation should go beyond repeats if they aim to characterise the full range of genomic variants. To allow future work with other taxonomic groups, we share our k-mer analysis pipeline, which is straightforward to run, relying largely on standard GNU command line tools.

Application-based guidelines for best practices in plant flow cytometry

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

Biogeography and genome size evolution of the oldest extant vascular plant genus, Equisetum (Equisetaceae)

BACKGROUND AND AIMS

Extant plant groups with a long fossil history are key elements in understanding vascular plant evolution. Horsetails (Equisetum, Equisetaceae) have a nearly continuous fossil record dating back to the Carboniferous, but their phylogenetic and biogeographic patterns are still poorly understood. We use here the most extensive phylogenetic analysis to date as a framework to evaluate their age, biogeography and genome size evolution.

METHODS

DNA sequences of four plastid loci were used to estimate divergence times and investigate the biogeographic history of all extant species of Equisetum. Flow cytometry was used to study genome size evolution against the framework of phylogenetic relationships in Equisetum.

KEY RESULTS

On a well-supported phylogenetic tree including all extant Equisetum species, a molecular clock calibrated with multiple fossils places the node at which the outgroup and Equisetum diverged at 343 Mya (Early Carboniferous), with the first major split among extant species occurring 170 Mya (Middle Jurassic). These dates are older than those reported in some other recent molecular clock studies but are largely in agreement with a timeline established by fossil appearance in the geological record. Representatives of evergreen subgenus Hippochaete have much larger genome sizes than those of deciduous subgenus Equisetum, despite their shared conserved chromosome number. Subgenus Paramochaete has an intermediate genome size and maintains the same number of chromosomes.

CONCLUSIONS

The first divergences among extant members of the genus coincided with the break-up of Pangaea and the resulting more humid, warmer climate. Subsequent tectonic activity most likely involved vicariance events that led to species divergences combined with some more recent, long-distance dispersal events. We hypothesize that differences in genome size between subgenera may be related to the number of sperm flagellae.

On genetic diversity in caraway: Genotyping of a large germplasm collection

Caraway (Carum carvi) is a widespread and frequently used spice and medicinal plant with a long history of cultivation. However, due to ongoing climatic changes, the cultivation is becoming increasingly risky. To secure caraway cultivation in future, timely breeding efforts to develop adapted material are necessary. Analysis of genetic diversity can accompany this process, for instance, by revealing untapped gene pools. Here, we analyzed 137 accessions using genotyping by sequencing (GBS). Hence, we can report a broad overview of population structure and genetic diversity of caraway. Population structure was determined using a principal coordinate analysis, a Bayesian clustering analysis, phylogenetic trees and a neighbor network based on 13,155 SNPs. Genotypic data indicate a clear separation of accessions into two subpopulations, which correlates with the flowering type (annual vs. biennial). Four winter-annual accessions were closer related to biennial accessions. In an analysis of molecular variance, genetic variation between the two subpopulations was 7.84%. In addition, we estimated the genome size for 35 accessions by flow cytometry. An average genome size of 4.282 pg/2C (± 0.0096 S.E.) was estimated. Therefore, we suggest a significantly smaller genome size than stated in literature.

Two cryptic species of California mustard within Caulanthus lasiophyllus

PREMISE

Cryptic species are evolutionarily distinct lineages lacking distinguishing morphological traits. Hidden diversity may be lurking in widespread species whose distributions cross phylogeographic barriers. This study investigates molecular and morphological variation in the widely distributed Caulanthus lasiophyllus (Brassicaceae) in comparison to its closest relatives.

METHODS

Fifty-two individuals of C. lasiophyllus from across the species' range were sequenced for the nuclear ribosomal internal transcribed spacer region (ITS) and the chloroplast trnL-F region. A subset of these samples were examined for the chloroplast ndhF gene. All 52 individuals were scored for 13 morphological traits, as well as monthly and annual climate conditions at the collection locality. Morphological and molecular results are compared with the closest relatives-C. anceps and C. flavescens-in the "Guillenia Clade." To test for polyploidy, genome size estimates were made for four populations.

RESULTS

Caulanthus lasiophyllus consists of two distinct lineages separated by eight ITS differences-eight times more variation than what distinguishes C. anceps and C. flavescens. Fewer variable sites were detected in trnL-F and ndhF regions, yet these data are consistent with the ITS results. The two lineages of C. lasiophyllus are geographically and climatically distinct; yet morphologically overlapping. Their genome sizes are not consistently different.

CONCLUSIONS

Two cryptic species within C. lasiophyllus are distinguished at the molecular, geographic, and climatic scales. They have similar genome sizes and are morphologically broadly overlapping, but an ephemeral basal leaf character may help distinguish the species.

Genome Size of 17 Species From Caelifera (Orthoptera) and Determination of Internal Standards With Very Large Genome Size in Insecta

Comparative studies of insect genome size show that Orthoptera is a unique group of Insecta with a significantly enlarged genome. To determine a suitable internal standard for a large genome and to compare the effects of different internal standards on estimates of genome size, we used four internal standards to estimate nuclear DNA content in nine insect species with large genomes. The results showed that the combination of two internal standards, Locusta migratoria (♂1C = 6.20 pg, ♀1C = 6.60 pg) and Periplaneta americana♂ (1C = 3.41 pg), was suitable for estimating large genome of Caelifera by flow cytometry. Using these two internal standards, we estimated the genome sizes of 17 species of Caelifera (12 genera in Acrididae, 2 genera in Pamphagidae, 1 genus in Pyrgomorphidae) using flow cytometry. Genomes ranged from 6.57 pg (Shirakiacris shirakii) to 18.64 pg (Bryodemella holdereri), the largest described in insects to date. These species showed significant genomic dimorphism based on sex: females had a 0.56 pg larger genome than males on average, which might be due to the sex chromosome determinism mechanism of X0(♂)/XX(♀). To test the results obtained by flow cytometry, we used k-mers of Illumina sequencing data to gauge the C-value of Calliptamus abbreviatus and Haplotropis brunneriana. The results of the two methods are slightly different. Genomes were estimated to be about 0.28 and 0.26 pg smaller, respectively, than the flow cytometry values. Furthermore, we also reconstructed the evolutionary relationships of these taxa and discuss the genome size evolution in a phylogenetic framework.

Nuclear Genome Size in Contrast to Sex Chromosome Number Variability in the Human Bed Bug, Cimex lectularius (Heteroptera: Cimicidae)

The human bed bug Cimex lectularius is one of the most prevalent human ectoparasites in temperate climate zones. The cytogenetic features of this resilient pest include holokinetic chromosomes, special chromosome behavior in meiosis, and numerical variation of chromosomes, where the diploid number ranges from 26 + X₁ X₂ Y to 26 + X_1-20 Y. It is desirable to assess the nuclear DNA content of various cytotypes for a further detailed study of the C. lectularius genome. Detailed knowledge of the DNA content of this parasite could also clarify the origin of additional chromosomes. The average nuclear genome size C. lectularius with 2n = 26 + X₁ X₂ Y is 2C = 1.94 pg for males and 1.95 pg for females. There is a significant correlation between genome size and the number of chromosomes, but in some specimens with additional chromosomes, nuclear genome size decreases or remains average. Several species used as the internal reference standard were tested for further investigations of genome size in C. lectularius, and the plant Solanum pseudocaspicum turned out to be the most suitable. © 2019 International Society for Advancement of Cytometry.

Nuclear DNA content of Pongamia pinnata L. and genome size stability of in vitro-regenerated plantlets

Pongamia pinnata L. is a multipurpose versatile legume that is well known as a prospective feedstock biodiesel species. However, to date, there has been little genomic research aimed at the exploitation of the biotechnological potential of this species. Genetic characterization of any plant is a challenging task when there is no information about the genome size and organization of the species. Therefore, the genome size of P. pinnata was estimated by flow cytometry with respect to two standards (Zea mays and Pisum sativum), and compared with that of in vitro-raised plants (nodal segment, in vitro-rooted plantlets and acclimatized in vitro plants) to study the potential effect of somaclonal variation on genome size. This method can be used to support the establishment of true-to-type plants to encourage afforestation programs. Modified propidium iodide/hypotonic citrate buffer was used for isolation of the intact nuclei. The 2C DNA value of this species was estimated to be 2.51 ± 0.01 pg. Statistically, there was no significant difference in the DNA content of the in vitro-grown plants and mother plant at α = 0.05. As a result of the low genome size of P. pinnata, a species that has adapted itself to a wide range of edaphic and ecological condition, we can now proceed for its next generation sequencing and genomic diversity studies.

Recent updates and developments to plant genome size databases

Two plant genome size databases have been recently updated and/or extended: the Plant DNA C-values database (http://data.kew.org/cvalues), and GSAD, the Genome Size in Asteraceae database (http://www.asteraceaegenomesize.com). While the first provides information on nuclear DNA contents across land plants and some algal groups, the second is focused on one of the largest and most economically important angiosperm families, Asteraceae. Genome size data have numerous applications: they can be used in comparative studies on genome evolution, or as a tool to appraise the cost of whole-genome sequencing programs. The growing interest in genome size and increasing rate of data accumulation has necessitated the continued update of these databases. Currently, the Plant DNA C-values database (Release 6.0, Dec. 2012) contains data for 8510 species, while GSAD has 1219 species (Release 2.0, June 2013), representing increases of 17 and 51%, respectively, in the number of species with genome size data, compared with previous releases. Here we provide overviews of the most recent releases of each database, and outline new features of GSAD. The latter include (i) a tool to visually compare genome size data between species, (ii) the option to export data and (iii) a webpage containing information about flow cytometry protocols.

Development of flow cytometric protocol for nuclear DNA content estimation and determination of chromosome number in Pongamia pinnata L., a valuable biodiesel plant

The potentiality of Pongamia pinnata L. as a sustainable source of feedstock for the biodiesel industry is dependent on an extensive knowledge of the genome structure of the plant. Flow cytometry, with propidium iodide (PI) as the DNA stain, was used to estimate the nuclear DNA content of P. pinnata, with respect to Zea mays 'CE-777' as standard. The internal and pseudo-internal standardization was followed on account of the inhibitory effect of secondary compounds on PI intercalation. The antioxidants (PVP-40 and β-mercaptoethanol) were added to the nuclear isolation buffer for the reduction of inhibitory effect of P. pinnata cytosol. Nuclear DNA content estimation was done for P. pinnata leaves from different altitudes (37-117 m height from sea level) of Assam. Flow cytometry analysis indicated that the nuclear DNA content of P. pinnata is 2.66 pg with predicted 1C value of 1,300 Mb using Z. mays as standard. Coefficient of variation in flow cytometric analysis was within the limit of 5 % indicating that the results were reliable. Somatic chromosome numbers were counted from root-tip cells and was found to be 2n = 22 corresponding to the diploid level (x = 11). A decreasing trend in the nuclear DNA content was observed for the species of different altitudes.

Genome size variation among and within Camellia species by using flow cytometric analysis

BACKGROUND

The genus Camellia, belonging to the family Theaceae, is economically important group in flowering plants. Frequent interspecific hybridization together with polyploidization has made them become taxonomically "difficult taxa". The DNA content is often used to measure genome size variation and has largely advanced our understanding of plant evolution and genome variation. The goals of this study were to investigate patterns of interspecific and intraspecific variation of DNA contents and further explore genome size evolution in a phylogenetic context of the genus.

METHODOLOGY/PRINCIPAL FINDINGS

The DNA amount in the genus was determined by using propidium iodide flow cytometry analysis for a total of 139 individual plants representing almost all sections of the two subgenera, Camellia and Thea. An improved WPB buffer was proven to be suitable for the Camellia species, which was able to counteract the negative effects of secondary metabolite and generated high-quality results with low coefficient of variation values (CV) <5%. Our results showed trivial effects on different tissues of flowers, leaves and buds as well as cytosolic compounds on the estimation of DNA amount. The DNA content of C. sinensis var. assamica was estimated to be 1C = 3.01 pg by flow cytometric analysis, which is equal to a genome size of about 2940 Mb.

CONCLUSION

Intraspecific and interspecific variations were observed in the genus Camellia, and as expected, the latter was larger than the former. Our study suggests a directional trend of increasing genome size in the genus Camellia probably owing to the frequent polyploidization events.

Genome size in Anthurium evaluated in the context of karyotypes and phenotypes

BACKGROUND AND AIMS

Anthurium is an important horticultural crop from the family Araceae, order Alismatales, a lineage considered to have diverged from other monocots prior to the cereals. Genome size and its distribution in Anthurium were investigated to gain a basic understanding of genome organization in this large genus and to forge a firm foundation for advancement of molecular approaches for the study of Anthurium. Currently, genome size estimates have been reported for only two Anthurium samples.

METHODOLOGY

Bulk nuclear DNA content estimates were obtained by flow cell cytometry using leaf tissue collected from Anthurium species of different subgeneric groups and from commercial cultivars. The most current and well-supported topology of subgeneric, sectional relationships was applied to present genome size estimates in the context of reported chromosome counts, karyotypes, putative phylogenetic relationships, observed phenotypes and pedigree.

PRINCIPAL RESULTS

Genome size estimates based on bulk nuclear DNA content for 77 accessions representing 34 species and 9 cultivars were obtained, including initial estimates for 33 Anthurium species, and both the smallest (Anthurium obtusum; Tetraspermium) and largest (Anthurium roseospadix; Calomystrium) Anthurium genome sizes reported to date. Genome size did not distinguish any subgeneric section, but ranged 5-fold (4.42-20.83 pg/2 C) despite consistent 2N= 30 chromosome counts. Intraspecies genome size variation >20 % is reported for Anthurium ravenii, A. watermaliense and A. gracile.

CONCLUSIONS

Genome size estimates for Anthurium species spanning 13 recognized subgeneric sections indicate that genome size does not generally correlate with chromosome count or phylogenetic relationships. Mechanisms of genome expansion and contraction, including amplification and reduction of repetitive elements, polyploidy, chromosome reorganization/loss, may be involved in genome evolution in Anthurium as in other species. The new information on Anthurium genome sizes provides a platform for molecular studies supporting further research on genome evolution as well as cultivar development.

C-value reassessment of plant standards: an image cytometry approach

Image cytometry (ICM) has been used to measure DNA 2C-values by evaluating the optical density of Feulgen-stained nuclei. This optical measurement is carried out using three basic tools: microscopy, digital video camera, and image analysis software. Because ICM has been applied to plants, some authors have remarked that studies should be performed before this technique can be accepted as an accurate method for determination of plant genome size. Based on this, the 2C-value of eight plants, which are widely used as standards in DNA quantifications, was reassessed in a cascade-like manner, from A. thaliana through R. sativus, S. lycopersicum, Glycine max, Z. mays, P. sativum, V. faba, to A. cepa. The mean 2C-values of all plants were statistically compared to the values reported by other authors using flow cytometry and/or ICM. These analyses demonstrated that ICM is an accurate and reliable method for 2C-value measurement, representing an attractive alternative to flow cytometry. Statistical comparison of the results also indicated Glycine max 'Polanka' as the most adequate primary standard. However, distinct authors have been advised that 2C DNA content of the reference standard should be close to that of the sample. As three further approaches also revisited the 2C-value of these eight plants, we have thus proposed a mean 2C-value for each eight species.

Applications of flow cytometry in plant pathology for genome size determination, detection and physiological status

Flow cytometers are probably the most multipurpose laboratory devices available. They can analyse a vast and very diverse range of cell parameters. This technique has left its mark on cancer, human immunodeficiency virus and immunology research, and is indispensable in routine clinical diagnostics. Flow cytometry (FCM) is also a well-known tool for the detection and physiological status assessment of microorganisms in drinking water, marine environments, food and fermentation processes. However, flow cytometers are seldom used in plant pathology, despite FCM's major advantages as both a detection method and a research tool. Potential uses of FCM include the characterization of genome sizes of fungal and oomycete populations, multiplexed pathogen detection and the monitoring of the viability, culturability and gene expression of plant pathogens, and many others. This review provides an overview of the history, advantages and disadvantages of FCM, and focuses on the current applications and future possibilities of FCM in plant pathology.

Organisation of the plant genome in chromosomes

The plant genome is organized into chromosomes that provide the structure for the genetic linkage groups and allow faithful replication, transcription and transmission of the hereditary information. Genome sizes in plants are remarkably diverse, with a 2350-fold range from 63 to 149,000 Mb, divided into n=2 to n= approximately 600 chromosomes. Despite this huge range, structural features of chromosomes like centromeres, telomeres and chromatin packaging are well-conserved. The smallest genomes consist of mostly coding and regulatory DNA sequences present in low copy, along with highly repeated rDNA (rRNA genes and intergenic spacers), centromeric and telomeric repetitive DNA and some transposable elements. The larger genomes have similar numbers of genes, with abundant tandemly repeated sequence motifs, and transposable elements alone represent more than half the DNA present. Chromosomes evolve by fission, fusion, duplication and insertion events, allowing evolution of chromosome size and chromosome number. A combination of sequence analysis, genetic mapping and molecular cytogenetic methods with comparative analysis, all only becoming widely available in the 21st century, is elucidating the exact nature of the chromosome evolution events at all timescales, from the base of the plant kingdom, to intraspecific or hybridization events associated with recent plant breeding. As well as being of fundamental interest, understanding and exploiting evolutionary mechanisms in plant genomes is likely to be a key to crop development for food production.

Distribution and habitat segregation on different spatial scales among diploid, tetraploid and hexaploid cytotypes of Senecio carniolicus (Asteraceae) in the Eastern Alps

BACKGROUND AND AIMS

The spatial distribution of cytotypes can provide valuable insights into evolutionary patterns of polyploid complexes. In a previous study the macro-scale distribution of the three main cytotypes in Senecio carniolicus (Asteraceae) within the Eastern Alps was characterized. Employing a roughly 12-fold extended sampling, the present study focuses on unravelling patterns of cytotype distribution on the meso- and microscale and on correlating those with ecological properties of the growing sites.

METHODS

DAPI flow cytometry of dried samples was used to determine DNA ploidy level in 5033 individuals from 100 populations spread over the entire Eastern Alpine distribution area of S. carniolicus. Descriptors of microhabitats as well as spatial data were recorded in the field, and analysed with a mixed-effects ANOVA.

KEY RESULTS

Extensive variation in DNA ploidy levels (2x, 3x, 4x, 5x, 6x, 7x, 8x, 9x) was detected. Of the main cytotypes, diploids and hexaploids were widespread and had strongly overlapping distributions resulting in the frequent occurrence of cytotype mixtures (half of the investigated populations), whereas tetraploids were disjunctly distributed and occurred in the south-west and the east of the species' distribution area. In spite of the frequent co-occurrence of cytotypes, only 1 % of the samples belonged to secondary cytotypes (3x, 5x, 7x, 8x, 9x). Diploids, tetraploids and hexaploids were altitudinally segregated, but with broad overlap. Similarly, highly significant differences in vegetation and rock cover as well as microhabitat exposure were found between the main cytotypes.

CONCLUSIONS

Senecio carniolicus shows a remarkable diversity of cytotypes. The distribution of the three main cytotypes (2x, 4x, 6x) has been shaped by Pleistocene glaciations to different extents. Whereas tetraploids are nearly entirely restricted to refugia, hexaploids colonized areas that were extensively glaciated. Diploid and hexaploid individuals often co-occur in mixed populations, where they are spatially and ecologically segregated at both the meso-scale (altitudinal differentiation, exposure of the growing site) and the micro-scale (cover of vegetation and bare rock). With regard to the ecological parameters investigated, the tetraploid cytotype occupies an intermediate position. The rareness of secondary cytotypes suggests the presence of strong pre- or post-zygotic mating barriers.