Nuclear DNA Amounts in Chinese Bryophytes Estimated by Flow Cytometry: Variation Patterns and Biological Significances

There exists an obvious gap in our knowledge of the nuclear DNA amount of bryophytes, not only in terms of the low number of species represented, but also in systematic and geographic representation. In order to increase our knowledge of nuclear DNA amounts and variation patterns in bryophytes, and their potential phylogenetic significances and influences on phenotypes, we used flow cytometry to determine the DNA 1C values of 209 bryophyte accessions, which belong to 145 mosses and 18 liverworts collected from China, by using Physcomitrella patens as a standard. We quantified the differences in DNA 1C values among different orders and families and constructed a phylogenetic tree of 112 mosses with four gene sequences (nad5, rbcL, trnL-F, and 18S-ITS1-5.8S-ITS2-26S). DNA 1C values were mapped onto the phylogenetic tree to test a potential phylogenetic signal. We also evaluated the correlations of the DNA 1C value with the sizes of individuals, leaves, cells, and spores by using a phylogenetically controlled analysis. New estimates of nuclear DNA amounts were reported for 145 species. The DNA 1C values of 209 bryophyte accessions ranged from 0.422 pg to 0.860 pg, with an average value of 0.561 pg, and a 2.04-fold variation covered the extremes of all the accessions. Although the values are not significantly different (p = 0.355) between mosses (0.528 pg) and liverworts (0.542 pg), there are variations to varying extents between some families and orders. The DNA 1C value size exerts a positive effect on the sizes of plants, leaves, and cells, but a negative effect on spore size. A weak phylogenetic signal is detected across most moss species. Phylogenetic signals are comparatively strong for some lineages. Our findings show that bryophytes have very small and highly constrained nuclear DNA amounts. There are nucleotype effects of nuclear DNA amounts for bryophytes at the individual, organ, and cell levels. We speculate that smaller nuclear DNA amounts are advantageous for bryophytes in dry environments. Significant differences in the DNA 1C values among some moss families and orders, as well as phylogenetic signals for some lineages, imply that nuclear DNA amount evolution in mosses seems to be unidirectional.


Introduction
The amount of DNA in the whole chromosome complement irrespective of the ploidy level of the organism is referred to as the DNA C-value, while the amount of nuclear DNA in an unreplicated haploid chromosome complement is referred to as the genome size (DNA 1C value) [1,2]. Nuclear DNA amount is an important biodiversity character with fundamental biological significance [3,4], and remains vital for many diverse fields of biology.
Since 1976, lists of DNA amounts in different plant categories, complied for reference purposes, have been published [5], and the data have been pooled and released in a database since 1997. Nuclear DNA amount data for 12,273 species, comprising 10,770 angiosperms, 421 gymnosperms, 303 pteridophytes, 334 bryophytes (212 mosses, 109 liverworts, and 13 hornworts), and 445 algae, were included in the updated Version 4.0 of the Plant DNA C-value Database [2]. These lists and databases have been widely used for comparative studies [6][7][8][9]. However, the availability of nuclear DNA amount data varied widely among different plant categories [2]. Large gaps still exist in our knowledge of nuclear DNA amounts. Nuclear DNA amount data in many species and geographic regions have not yet been reported. Improving geographic and taxonomic representations were the two main targets in the following collection of nuclear DNA amount data [10,11].
DNA C-values in land plants (comprising bryophytes, lycophytes, monilophytes, gymnosperms, and angiosperms) vary ca. 1000-fold from approx. 0.11 to 127.4 pg [7]. The values highly vary among different angiosperm families; the average DNA 1C-vlaue differs 60-fold, from 0.62 pg (Salicaeae) to 26.58 pg (Alstroemeriaceae); nuclear DNA amounts are significantly lower in dicots than in monocots and are also significantly lower in nonperennials than in perennials [12]. However, the patterns of variation in nuclear DNA amounts among different bryophyte taxa are poorly understood.
Bryophytes are next to angiosperms in species number. Despite the numerous nuclear DNA amount data that exist in seed plants, the Plant DNA C-value Database revealed an obvious gap in our knowledge of the nuclear DNA amounts of bryophytes, not only in terms of the low number of species represented, but also in terms of systematic and geographic representation [2]. According to the updated DNA C-value Database (version 4.0), genome sizes (DNA 1C values) of 334 bryophyte species have been determined so far, mainly reported by Voglmayr (137 species) [13], Temsch et al. (77 species) [14], Bainard et al. (56 species) [10], Bainard (32 species) [15], Bainard and Villarreal (23 species) [16], and Greihuber et al. (5 species) [17]. Recently, the genome sizes of 33 moss species were reported by Bainard et al. [11]. Compared with those of angiosperms (approx. 4.2%), the nuclear DNA amount data of bryophyte species were relatively scarce (approx. 2.6% of ca. 12,800 bryophyte species, Goffinet and Shaw [18]). Concerning nuclear DNA amount data, taxonomic representation is also problematic. So far, no nuclear DNA amount estimate is available for more than 52% of bryophyte families. Additionally, although there is evidence of nuclear DNA amount data in some families, the proportion of their species with available nuclear DNA amount data was also rather low. For example, only 4 of the 447 species in the family Fissidentaceae, or less than 1% of the family's species, have been determined for their nuclear DNA amounts [2]. Despite the rich species of bryophytes in China, no nuclear DNA amount data were available for bryophyte accessions collected from China up until now.
Although previous studies revealed that genome size variation exhibits phylogenetic signals for some liverworts [10], mosses [11], and hornworts [16], more data were still needed to improve the systematic and geographic representation of the phylogenetic signal in bryophytes. Because of their dominant gametophytes, lack of vascular tissues, and poikilohydric strategy, bryophytes are unique among land plants [19]. Furthermore, bryophytes are characterized by their small size, high sensitivity to habitats and substrate specificity, and a short generation time, as well as their fast colonization-extinction rate [20,21]. However, there have been few studies on the relationship between bryophyte nuclear DNA amounts and morphological traits.
The objectives of the present work were to (1) increase our knowledge of nuclear DNA amounts and their variation patterns in bryophytes; (2) clarify whether nuclear DNA amounts exert potential nucleotype effects on the phenotypes of bryophytes, and their possible ecological significances; and (3) further confirm whether phylogenetic signals exist for nuclear DNA amounts in bryophytes.

General Aspects
The DNA 1C values of 209 bryophyte accessions were determined. These accessions included 145 mosses (belonging to 86 genera and 34 families) and 18 liverworts (belonging to 13 genera and 12 families). New DNA 1C values for 9 families, 64 genera, and 145 species were reported (Appendix A).

Relationships between DNA 1C Values and Morphological Traits
The DNA 1C values of the 112 mosses were positively related to plant sizes (p < 0.01), leaf length (p < 0.01) and width (p < 0.005), cell length (p < 0.1) and width (p < 0.5), but negatively related to spore sizes (p < 0.1). The above relationships also existed after taking into account phylogenetic non-independence ( Figure 3).

Phylogenetic Signal for DNA 1C Value among Moss Taxa
The phylogenetic tree included 112 moss species, accounting for 78.32% of the total moss species whose DNA 1C values were determined in the present study (Figures 4 and 5). Average DNA 1C values for these species ranged from 0.422 to 0.862 pg. The K-statistic for the whole tree with 112 moss species was 0.120, with a p-value of 0.203, indicating a weak phylogenetic signal for DNA 1C values across the tree. The phylogenetic signal was comparatively strong for some lineages, such as Dicranales (K-statistic = 0.488, p < 0.20) and its Clade A (K-statistic = 1.169, p < 0.05) (Figure 6), and Clade A of Hypnales (K-statistic = 2.429, p < 0.001) (Figure 7).

Discussion
The studies relevant to nuclear DNA amounts have been shown to be useful in elucidating taxonomic, evolutionary, and ecological problems in a number of plant taxa [9,12,[22][23][24][25][26]. Therefore, there is a continuing need to obtain more data on DNA 1C values for plant taxa.

Credibility of the Nuclear DNA Amounts We Determined
Our result of the moss accessions is mostly consistent with those of some previous reports. For example, a mean DNA 1C value of 0.52 pg for 209 moss accessions was recorded in the Bryophyte DNA C-value Database [2]. Voglmayr [13] reported a mean value of 0.509 pg for 137 moss species, Temsch et al. [27] reported a minimum of 0.39 pg and a maximum of 0.94 pg, and mean of 0.530 pg for 30 moss species, Greihuber et al. [17] reported a minimum of 0.44 pg, and maximum of 0.95 pg and a mean of 0.646 pg for five moss species, and Renzaglia, Rasch and Pike [28] reported an interspecific variation in DNA 1C values from 0.38 pg to 0.92 pg, with a mean of 0.544 pg for nine mosses by measuring their sperms using flow cytometry. Recently, Bainard et al. [11] [10]. The lower mean value of DNA 1C values for liverworts in this study was probably attributed to the insufficient liverwort samples in our work.
The DNA 1C value range is much broader in Voglmayr's paper, ranging from 0.174 to 2.16 pg, while it is much narrower in the present work, ranging from 0.422 pg to 0.860 pg. This is possibly due to different samples and determining methods. In fact, among the data of 289 accessions determined by Voglmayr [13] (determined by flow cytometer or by Feulgen densitometry), only 14 accessions each had a DNA 1C value larger than 1.0 pg and about 95% of the accessions had DNA 1C values smaller than 1.0 pg. In the DNA C-value database [2], among the 165 moss accessions determined by flow cytometry, only 7 accessions have a DNA 1C value larger than 1.0 pg, while about 96% of the accessions have a DNA 1C value smaller than 1.0 pg. Additionally, the data in the Kew database and Voglamyr's report were determined by flow cytometry or by Feulgen densitometry. According to our analyses of the data from the Kew database (version 4.0), the average 1C value (pg) determined by flow cytometry is 4.68 pg for 7626 angiosperm species, while that by Feulgen densitometry is 6.58 pg for 2931 angiosperm species. The case is the same for mosses. The average DNA 1C values determined by Feulgen densitometry (37 accessions) are 0.548 pg, while those determined by flow cytometry (165 accessions) are 0.509 pg. If the five largest data are excluded, the average DNA 1C values obtained by flow cytometry (160 accessions) are only 0.474 pg. Therefore, the data determined by flow cytometry are essentially smaller than those determined by Feulgen densitometry. The phenomenon that the data in the present work are generally smaller than those previously reported is possibly due to the fact that we used flow cytometry to determine the DNA 1C-values.

Variation Patterns of the Nuclear DNA Amounts in Bryophytes
The mean DNA 1C value of all the accessions was 0.53 pg, which was lower than that of 334 accessions (0.916 pg, DNA 1C value) in the Bryophyte DNA C-value Database [2]. Considering the high proportion of mosses in our collections, the result is reasonable and acceptable. Among the 209 accessions determined here, 188 are moss accessions; their DNA 1C values varied from 0.442 pg to 0.862 pg, with a mean value of 0.529 pg. The DNA 1C values of the mosses in the DNA C-value Database varied from 0.17 to 2.05 pg, with a mean of 0.519 pg [2].
Relatively weak interspecific variation in the DNA 1C value was detected for our accessions. In addition, 2.16-fold interspecific variation covers the extremes of the accessions, with an average value of 0.53 pg, a maximum of 0.862 pg in Rhodobryum giganteum, and a minimum of 0.398 pg in Macromitrium japonicum. This variation is much low compared with the 127.9-fold interspecific variation recorded in the Bryophyte DNA C-value Database [2]. Among the three groups of bryophytes in the database, liverworts have the largest DNA 1C value (a mean of 1.89 pg) and the largest interspecific variation (102 species, 97.43-fold), followed by mosses (a mean of 0.52 pg and an interspecific variation of 12.06-fold for 209 species), and hornworts have the smallest DNA 1C value (a mean of 0.249 pg) and interspecific variation (4.56-fold for 23 species). The 209 accessions did not contain hornworts and contained only 20 liverwort samples (accounting for less than 10% of the total accessions). Our result is mostly consistent with that of Voglmayr [13], who reported a two-fold interspecific variation in DNA 1C values from 0.3 to 0.6 pg for the majority of species. Although there existed a ca. 12-fold interspecific variation in DNA 1C values within mosses, 80% of the values were restricted to a range between 0.25 and 0.6 pg (2.4-fold variation) [13].
Among the ten moss orders we examined, the Bryales have the largest DNA 1C value (0.602 ± 0.022 pg), which is consistent with the speculation by Bainard et al. [11]. Considering significant differences in the DNA 1C value existed among some families, such as Leucobryaceae, Rhytidiaceae, Entodontaceae, Orthotrichaceae, Dicranaceae, Mniaceae, Bartramiaceae and Grimmiaceae (Table S1), the DNA 1C value evolution in mosses seems to be unidirectional.
Hookeria lucens was reported to have a relatively large DNA 1C value (1.61 pg) by Bainard et al. [11], but our investigation showed that Hookeria acutifolia Hook. & Grev. has a relatively small DNA 1C value (0.478 ± 0.017 pg). Therefore, more samples were needed to clarify the DNA 1C value throughout the order.

Nucleotype Effects and Possible Ecological Significance of Nuclear DNA Amounts in Bryophytes
Nuclear DNA content could affect the phenotype through the biophysical effects of its mass and volume, with the latter defined as nucleotype effects [29]. Nucleotype variation in nuclear DNA amount sets absolute limits on the minimum size and mass of cells. Such effects are additive in complex multicellular vascular plants, and the potential effects of the DNA 1C value can apply to cells, organs, and organisms, and act on many aspects of the life history of the plant [30]. In angiosperms, nuclear DNA amount positively correlated with the volume and weight of chromosomes [31][32][33], nuclear and cell sizes [34], epidermal cell size and leaf size of Lolium perenne L. [35], leaf width in the species of Nerine (Amaryllidaceae) [36], and plant height in Sencio in Australia [29]. However, the nucleotype effects of nuclear DNA content on phenotype have received little attention in bryophytes. Here, we detected significantly positive correlations of the DNA 1C value with plant size, leaf size, and cell size ( Figure 3), confirming the nucleotype effects of nuclear DNA content on the phenotype of bryophytes. The nucleotype effects of the DNA 1C value on the phenotype of bryophytes appeared at the individual, organ, and cell levels.
Early studies showed that the nuclear DNA amount negatively correlated with the duration of the mitotic and meiosis cycle [37,38], and minimum generation time [39]. Cutler et al. [40] suggested that smaller cells help plants to resist moisture stress because they maintain turgor with solute accumulation under lower water potential values compared to larger cells. Small cells often have small nuclear DNA amounts. According to Rejmánek [41], a low nuclear DNA content seems to be a result of selection for short minimum generation times in extreme cold environments. We speculate that a small nuclear DNA amount is also advantageous for bryophytes in dry environments. Bryophytes are generally sensitive to dry environments [21,42]. A small nuclear DNA amount allowed bryophytes to rapidly develop in a time-limited duration of favorable moisture availability in dry regions by a rapid mitotic cycle and a short duration of meiosis.
The spores of bryophytes are somewhat similar to the pollens of angiosperms as reproduction units. In angiosperms, significantly positive correlations of nuclear DNA content with pollen size were reported in Armeria maritima (Mill.) Willd. (Plumbaginaceae) [43], and in some cereal species [39]. However, we detected a negative relationship between spore size and DNA 1C value in 112 moss species (Figure 3). Löbel and Rydin [44] suggested that species with larger spores have a higher probability surviving in harsher habitats (e.g., dry habitats). Therefore, regardless of the fact that a small DNA 1C value allowed bryophytes to rapidly develop in a time-limited duration of favorable moisture availability in dry regions, another advantage for species with a smaller DNA 1C value and larger spores may be that larger spores have a higher probability to survive in a dry habitat. Proctor et al. [45], Baniaga et al. [26] and Bainard et al. [11] thought that desiccation tolerance might be an important selective pressure for plants to keep a relatively small DNA 1C value, which was consistent with our results. Nevertheless, more data will need to be collected before the relationship of spore size with the DNA 1C value in mosses can be clearly established.
To further rigorously clarify the relationships between nuclear DNA amounts and ecological adaptation for mosses, we must collect the global geographic distribution data of the mosses with known nuclear DNA amounts, and the corresponding climate data of these distribution points, then quantify the relationships among the climates, nuclear DNA amounts, and morphological traits. This will possibly allow us to better understand the ecological significances of the nuclear DNA amounts in mosses.

Phylogenetic Signals of Nuclear DNA Amounts in Bryophytes
According to Leitch et al. [2], the DNA 1C value is much smaller (a mean value of 0.52 pg) in mosses than in angiosperms (a mean of 5.13 pg), and the interspecific variation is much weaker in mosses (12.04-fold for most species) than in angiosperms (ca. 2000-fold). What is the reason that the DNA 1C values are so small and constant in mosses? DNA 1C value variation is likely to be a whole-organism phenomenon that can be studied at the developmental and ecological levels [46]. Renzaglia et al. [28] thought that the selection of biflagellated sperm may have favored a low nuclear DNA amount.
The variations in DNA amount have been found to be linked with phylogenetic signals across land plants [7], flowering plants [47,48] and a number of angiosperm taxa, such as Allium [49], Capsicum (Solanaceae) [50], Poaceae [51], and Bromelioideae of Bromeliaceae [52]. Bainard and Villarreal [16] reported a 20.46-fold interspecific variation in the DNA 1C-value from 0.27 to 20.46 pg for 67 hornwort species from 33 families using flow cytometry and detected a strong phylogenetic signal of DNA 1C-value across the liverwort phylogeny. Recently, Bainard et al. [11] detected a phylogenetic signal of DNA 1C values across the phylogeny of mosses based on the data they determined and those from previous studies. In our moss sampling, there existed only a two-fold interspecific variation in the DNA C-value from 0.422 to 0.862 pg and a weak phylogenetic signal across the phylogenetic tree, which we produced based on four gene regions available in the NCBI database. However, for some lineages of Dicranales and Hypnales, the phylogenetic signal was comparatively strong, and the variation in the DNA 1C value was roughly correlated with their phylogenetic relatedness (Figures 6 and 7). The above results were consistent with that of Bainard et al. [11]. The very small range of DNA 1C values across the 145 mosses we examined likely indicated that the DNA 1C values remain constrained in mosses and there has not been much divergence in the DNA 1C values over evolutionary history, which is consistent with the result of Baniaga et al. [26] that small nuclear genomes of Selaginella were associated with a low rate of DNA 1C value evolution. Additionally, many new DNA 1C value estimates reported for 145 bryophyte species are valuable for a better understanding of the phylogenetic signal of DNA 1C values across the phylogeny of the whole bryophyte group.

Materials
The plant materials are presented in Appendix A. Shuiliang Guo and Tong Cao identified voucher specimens, which were deposited at the bryophyte herbarium of Shanghai Normal University (SHTU). The taxonomy of species mainly follows the work of Jia and He [53].

Nuclei Isolation
The protocol for isolating nuclei was adapted from that of Johnston et al. [54]. Bryophyte tissues (5 to 15 moss fresh shoot tips and 1-2 cm 2 fresh liverwort tissue, ca. 10 mg air-dried tissue) were washed to remove soil, chemicals, and other organisms that might react with the chemicals and alter the results. The tissue was chopped at room temperature, with a razor blade in about 0.55 mL of isolation buffer to homogenize the tissues and release the nuclei. The composition of the isolation buffer (200 mL) contained 45 mM MgCl 2 , 30 mM sodium citrate, 20 mM MOPS and 0.1% (w/v) Triton X−100 (reminding deionized water, pH 7.0) [55]. The nucleus suspension was then filtered (with a 10 mL syringe and a 30 µm nylon mesh) to remove debris that might block the flow cell. The nuclei suspension was filtered into a 1.5 mL tube to centrifuge at 1600 r/min for 5 min. After removal of the supernatant liquid, the nuclei were stained with 150 µg·mL −1 propidium iodide in the presence of 0.5 µg·mL −1 RNase. The mixture was dyed at 4 • C under the dark for 20 min. From each species, three accessions were randomly selected for DNA amount measurement [56].

Nuclear DNA Amount Measurement
Physcomitrella patens (Hedw.) Bruch & Schimp. was used as a standard (0.53 pg/DNA 1C value) [57] because the material of the species was widely available, quick and easy to grow, suitable for FCM protocols, and with an appropriate genome size for bryophytes [11]. Additionally, Physcomitrella patens was the first non-seed plant to have its genome sequenced, with verified genome size stability within well-delimited species [58]. The species has very few secondary compounds, which will interfere with quantitative DNA staining [59]. Therefore, P. patens has been used as a standard in measuring the genome size of seed plants [28] and bryophytes [60]. We used the gametophytes of the species, which were cultured in a growth chamber.
A flow cytometer (FACSCalibur, BD Bioscience, Mountain View, CA, USA) was used for the nuclei suspension analysis. The laser-emission wavelength was adjusted to 488 nm. Each sample consisted of 300 µL of nuclei suspension, and analysis was conducted at a data rate of 100-150 nuclei per second. The histogram was analyzed by using the ModFit LT software to obtain the G0/G1 peak (namely the fluorescence value), and the variation coefficient of the G0/G1 peak (CV% = standard deviation/mean ×100) [61]. If the CV value is less than 5%, the results are acceptable; otherwise, they are not [62].
To determine the DNA 1C value, the relative position of the bryophyte 1C peak was compared to the position of the standard 1C peak. The 1C peak was observed for the bryophyte as the haploid (gametophytic) tissue was analyzed. The nuclear DNA amount of the sample could be calculated as follows: Sample 1C value (DNA pg) = Standard 1C value × sample 1C meak peak position of three replicates standard 1C mean peak position of three replicates (1) As is the case with most moss species, Physcomitrella patens is a typical endopolyploidy species [63] (Figure 8). If using the internal standard method, the 1C, 2C, and 4C peaks of the standard and the peaks of the sample are not easily separated and identified. Thus, we compared the genome size of Rhodobryum giganteum (Schwägr.) Paris determined by using P. patens as an internal standard with the value obtained by using P. patens as an external standard. We found that the genome size of R. giganteum obtained by using the external standard method was 0.862 ± 0.006, while that of the internal standard method was 0.860 ± 0.003 (with three replicates) (Figure 8), revealing that their difference was insignificant (p > 0.5). Therefore, the genome sizes were estimated in the present study by using P. patens as an external standard. Flow cytometry using an external standard was also used to measure the nuclear DNA amount in previous works [64][65][66][67][68][69][70][71]. To control the negative influences of the external standard method on the result as much as possible, the instrument settings were adjusted to control the "drift" in the peak location over time, and the histogram of the standard (P. patens) was obtained every time for comparison with that of the new species. Both the sample and the standard were measured with three replicates for each species.

Data Analysis
All the data were expressed as means ± standard errors with three replicates. One-way analysis of variance (ANOVA) was employed to test the differences among the taxa in their genome sizes using the procedures in the SPSS 22.0 statistical package (IBM Corp., Armonk, NY, USA). The least significant difference (LSD) method was employed.
To test whether the DNA 1C values are of evolutionary significance, we constructed a phylogenetic tree including 112 out of the 143 moss species whose genome sizes had been estimated. These species were selected because their four gene regions (nad5, rbcL, trnL-F; 18S-ITS1-5.8S-ITS2-26S) were available in the NCBI database (Table S2). We did not perform phylogenetic analyses of liverworts due to insufficient sampling.
Sequence chromatograms were compiled using Seqman II (DNASTAR Inc., Madison, WI, USA) and then automatically aligned in PhyDE 0.9971 [72]. Regions of partially incomplete data at the beginning and end of sequences were excluded from subsequent analyses. Gaps were treated as missing data. A total of 6106 base pairs, which included 3304 variable sites and 2477 parsimony-informative sites, were used to construct the phylogenetic tree.
We used MrModeltest v. 2.4 [73], which is incorporated in PAUP 4.0a168 [74], to select the best-fit nucleotide substitution model for each gene according to the corrected Akaike information criterion (AICc). The relevant parameters were set accordingly for each compartment. The phylogenetic tree was constructed using RAxML 8.2.10 [75]. The trees were visualized and annotated in TreeGraph 2 [76].
Using the ContMap function with default settings from the phytools package, DNA 1C values were mapped onto the phylogenetic tree. The K-statistic of the phylogenetic signal was calculated by using the phylosig function from the phytools package [77].
To detect the possible relationships between the DNA 1C value and morphological traits, data on plant size (at the individual level), leaf length and width (at the organ level), cell length and width, and spore diameter (at the cell level) were collected from relevant literature (Table S3). Plant size refers to the height of the main stem for acrocarpous mosses or the creeping stem for pleurocaropus mosses. Leaf length and width refer to those of branch leaves. Cell length and width refer to those of the median cells in branch leaves. The mean values of the above indices were used in relevant analyses.
The correlations of genome sizes with morphological traits were analyzed using ordinary least squares. A phylogenetically controlled analysis using the picante package [78] in R was also performed to fit a linear model to reveal the above relationships for the 112 species in the phylogenetic tree, which takes into account phylogenetic nonindependence between data points. Data of spore size were only available for 102 moss species from relevant literature (Table S3); thus, we analyzed the relationships of spore size with genome size within these moss species.
Supplementary Materials: The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/plants12071564/s1, Table S1: Significance levels of difference among families in genome sizes (each with more than three species) ; Table S2: List of the 112 moss  species including GenBank numbers for the four gene regions; Table S3: The data on plant size (at the individual level), leaf length and width (at the organ level), cell length and width, and spore diameter (at the cell level).

Data Availability Statement:
The data presented in this study are available from the corresponding author upon request. Note: Newly reported families, genera and species are marked with three asterisks, two asterisks, and one asterisk, respectively.