Electronic Supplementary Material
-
Materials and Methods
-
Results
-
References
-
Materials and Methods
Measuring Osteocyte Lacunae
The method of osteocyte lacuna size measuring is carried out using an automated approach. We utilize scanning electron microscopy to photograph the thin section slides to obtain the finest available resolution. Since lacunae are so small, a difference of a few pixel widths can make a large difference in calculating the area of an osteocyte lacuna. The threshold of the image is altered in the program ImageJ (Abramoff et al. 2004), so that the osteocyte lacunae are easily distinguishable from the rest of the bone matrix. The “Analyze Particles” feature is utilized to outline and automatically measure the lacunar areas in the photograph. The length and width calculated by this algorithm were subsequently used to calculate cell volume using the equation (4/3 x π x width axis radius2 x length axis radius) like in Organ et al. (2007).
Estimating Missing Values Using Phylogeny
The calculated volume values for each bone of each of the four studied taxa were then added to the most complete lacunar and genome size data set from Organ et al. (2009). As opposed to estimating the missing genome size values using the program BayesTraits (as used in Organ et al. (2007), Organ & Shedlock 2009, and Organ et al. (2009), a similar web based program called PhyloPars was employed. PhyloPars is extremely user friendly, able to easily handle missing data, and uses a maximum likelihood based method for inferring missing values (see Bruggeman et al. 2009 for more details). Additionally, the genome sizes used to check the accuracy of the estimated C-values from PhyloPars were obtained from the Animal Genome Size Database. To be consistent, only C-values obtained using the method of Feulgen Densitometry were averaged if more than one genome size was available for the given taxon.
The phylogenetic tree used in analysis was the extant taxon tree from Organ and Shedlock (2009). Each taxon in this study was added to this tree using Mesquite v. 2.73(Maddison & Maddison 2007) and the StratAdd package (Faure et al. 2006). Details of how this tree was assembled can be found in the electronic supplementary information of Organ et al. (2009).
In this sensitivity analysis, the “ln genome” for the 35 extant taxa was replaced with the minimum and maximum C-values reported in the Animal Genome Size Database. Subsequently, all of the analyses were rerun for each the minimum C-value and maximum C-value. It was shown that no matter which C-value (min, max or average) was used, the estimated genome size for the unknown value did not significantly differ (Table S3). Additionally, using the 3rd quartile of the cell size measurements did not make a significant difference in the estimated genome size as compared to using the mean (Table S4). It is important to note that Organ et al. (2011) has presented an improved multiple regression method for estimating the genome sizes of amphibians- that have inordinately large genomes that improves the accuracy of the estimation. We were not able to use this method due to the fact that the program we used, PhyloPars, does not support multiple regressions. When estimating genome size for a salamander or related amphibian, the multiple regression method takes into account the different relationship between osteocyte lacunae size and genome size, so it should be more accurate.
-
Results
Table S1: A summary of measured values from the 44 slides. The natural log of the mean lacuna volume presented here is the value that was used in reconstructing the genome size for that skeletal element in the program PhyloPars. The column “n lacunae” is how many lacunae were automatically measured per slide.
Taxon
|
Bone element
|
Ln Lacuna Volume (Mean)
|
Standard Error
|
n lacunae
|
Alligator sinensis
|
Tibia
|
3.908
|
0.749
|
15
|
Alligator sinensis
|
Fibula
|
3.776
|
0.666
|
13
|
Alligator sinensis
|
Rib
|
3.637
|
0.915
|
7
|
Alligator sinensis
|
Ulna
|
3.653
|
0.839
|
54
|
Alligator sinensis
|
Femur
|
4.953
|
0.706
|
10
|
Alligator sinensis
|
Thoracic
|
4.970
|
0.787
|
6
|
Alligator sinensis
|
Metatarsal
|
4.121
|
0.800
|
47
|
Alligator sinensis
|
Humerus
|
3.678
|
0.757
|
47
|
Alligator sinensis
|
Skull
|
4.588
|
0.560
|
8
|
Alligator sinensis
|
Caudal
|
2.539
|
1.646
|
11
|
Alligator sinensis
|
Radius
|
4.285
|
1.118
|
17
|
Ambystoma tigrinum
|
Tibia
|
8.634
|
0.867
|
7
|
Ambystoma tigrinum
|
Fibula
|
8.934
|
0.735
|
16
|
Ambystoma tigrinum
|
Rib
|
10.274
|
0.371
|
3
|
Ambystoma tigrinum
|
Ulna
|
8.584
|
0.615
|
8
|
Ambystoma tigrinum
|
Femur
|
9.172
|
0.572
|
6
|
Ambystoma tigrinum
|
Thoracic
|
9.378
|
0.584
|
5
|
Ambystoma tigrinum
|
Metatarsal
|
9.270
|
1.007
|
6
|
Ambystoma tigrinum
|
Humerus
|
8.192
|
0.740
|
10
|
Ambystoma tigrinum
|
Skull
|
8.979
|
0.699
|
7
|
Ambystoma tigrinum
|
Caudal
|
7.979
|
1.258
|
6
|
Ambystoma tigrinum
|
Radius
|
7.837
|
1.046
|
8
|
Marmota monax
|
Tibia
|
3.517
|
0.741
|
36
|
Marmota monax
|
Fibula
|
3.628
|
0.772
|
37
|
Marmota monax
|
Rib
|
3.655
|
0.851
|
42
|
Marmota monax
|
Ulna
|
4.076
|
1.200
|
46
|
Marmota monax
|
Femur
|
3.900
|
0.885
|
71
|
Marmota monax
|
Thoracic
|
4.320
|
1.028
|
46
|
Marmota monax
|
Metatarsal
|
4.997
|
1.230
|
40
|
Marmota monax
|
Humerus
|
4.874
|
1.651
|
60
|
Marmota monax
|
Skull
|
3.807
|
1.570
|
28
|
Marmota monax
|
Caudal
|
4.764
|
1.109
|
25
|
Marmota monax
|
Radius
|
4.338
|
1.045
|
54
|
Columba livia
|
Tibia
|
4.128
|
0.862
|
14
|
Columba livia
|
Fibula
|
3.447
|
0.501
|
6
|
Columba livia
|
Rib
|
4.070
|
0.921
|
8
|
Columba livia
|
Ulna
|
3.233
|
0.864
|
26
|
Columba livia
|
Femur
|
3.047
|
0.664
|
17
|
Columba livia
|
Thoracic
|
4.449
|
1.653
|
12
|
Columba livia
|
Metatarsal
|
3.011
|
0.723
|
8
|
Columba livia
|
Humerus
|
3.223
|
0.941
|
10
|
Columba livia
|
Skull
|
4.059
|
1.227
|
8
|
Columba livia
|
Caudal
|
3.950
|
1.621
|
6
|
Columba livia
|
Radius
|
3.104
|
0.836
|
8
|
|
Table S2: A summary of the genome sizes estimated from the program PhyloPars for each skeletal element. The estimated genome size in picograms (pg) is the exponential of the estimated ln genome value.
Taxon
|
Bone element
|
Est. ln genome
|
SE
|
Est. Genome in pg
|
Alligator sinensis
|
Tibia
|
0.372
|
0.126
|
1.451
|
Alligator sinensis
|
Fibula
|
0.317
|
0.126
|
1.373
|
Alligator sinensis
|
Rib
|
0.253
|
0.126
|
1.288
|
Alligator sinensis
|
Ulna
|
0.261
|
0.126
|
1.298
|
Alligator sinensis
|
Femur
|
0.813
|
0.126
|
2.255
|
Alligator sinensis
|
Thoracic
|
0.821
|
0.126
|
2.273
|
Alligator sinensis
|
Metatarsal
|
0.461
|
0.126
|
1.586
|
Alligator sinensis
|
Humerus
|
0.274
|
0.126
|
1.315
|
Alligator sinensis
|
Skull
|
0.660
|
0.126
|
1.935
|
Alligator sinensis
|
Caudal
|
-0.209
|
0.126
|
0.811
|
Alligator sinensis
|
Radius
|
0.533
|
0.126
|
1.704
|
Ambystoma tigrinum
|
Tibia
|
2.600
|
0.457
|
13.464
|
Ambystoma tigrinum
|
Fibula
|
2.720
|
0.457
|
15.180
|
Ambystoma tigrinum
|
Rib
|
3.292
|
0.457
|
26.897
|
Ambystoma tigrinum
|
Ulna
|
2.580
|
0.457
|
13.197
|
Ambystoma tigrinum
|
Femur
|
2.830
|
0.457
|
16.945
|
Ambystoma tigrinum
|
Thoracic
|
2.910
|
0.457
|
18.357
|
Ambystoma tigrinum
|
Metatarsal
|
2.870
|
0.457
|
17.637
|
Ambystoma tigrinum
|
Humerus
|
2.410
|
0.457
|
11.134
|
Ambystoma tigrinum
|
Skull
|
2.750
|
0.457
|
15.643
|
Ambystoma tigrinum
|
Caudal
|
2.320
|
0.457
|
10.176
|
Ambystoma tigrinum
|
Radius
|
2.260
|
0.457
|
9.583
|
Marmota monax
|
Tibia
|
0.440
|
0.212
|
1.553
|
Marmota monax
|
Fibula
|
0.487
|
0.212
|
1.627
|
Marmota monax
|
Rib
|
0.496
|
0.212
|
1.642
|
Marmota monax
|
Ulna
|
0.678
|
0.212
|
1.970
|
Marmota monax
|
Femur
|
0.780
|
0.212
|
2.181
|
Marmota monax
|
Thoracic
|
0.780
|
0.212
|
2.181
|
Marmota monax
|
Metatarsal
|
1.070
|
0.212
|
2.915
|
Marmota monax
|
Humerus
|
1.010
|
0.212
|
2.746
|
Marmota monax
|
Skull
|
0.563
|
0.212
|
1.756
|
Marmota monax
|
Caudal
|
0.966
|
0.212
|
2.627
|
Marmota monax
|
Radius
|
0.788
|
0.212
|
2.199
|
Columba livia
|
Tibia
|
0.411
|
0.240
|
1.508
|
Columba livia
|
Fibula
|
0.123
|
0.240
|
1.131
|
Columba livia
|
Rib
|
0.386
|
0.240
|
1.471
|
Columba livia
|
Ulna
|
0.030
|
0.240
|
1.030
|
Columba livia
|
Femur
|
-0.050
|
0.240
|
0.951
|
Columba livia
|
Thoracic
|
0.550
|
0.240
|
1.733
|
Columba livia
|
Metatarsal
|
-0.064
|
0.240
|
0.938
|
Columba livia
|
Humerus
|
0.025
|
0.240
|
1.025
|
Columba livia
|
Skull
|
0.382
|
0.240
|
1.465
|
Columba livia
|
Caudal
|
0.335
|
0.240
|
1.398
|
Columba livia
|
Radius
|
-0.025
|
0.240
|
0.975
|
Table S3: A summary of estimated values from the sensitivity analysis, where maximum, minimum and average C-values were used in three separate sets of analyses to determine if the value selected for C-value from the Animal Genome Size Database affected the estimated genome size for the one unknown value. The difference in mean estimated genome size is not significant with the standard error taken into account.
Taxon
|
Bone
|
Est. ln genome (Max C-val)
|
SE
|
Est. ln genome (Min C-val)
|
SE
|
Ambystoma tigrinum
|
Tibia
|
2.500
|
0.483
|
2.651
|
0.62
|
Ambystoma tigrinum
|
Caudal
|
2.274
|
0.483
|
2.334
|
0.62
|
Ambystoma tigrinum
|
Radius
|
2.225
|
0.483
|
2.266
|
0.62
|
Ambystoma tigrinum
|
Fibula
|
2.605
|
0.483
|
2.797
|
0.62
|
Ambystoma tigrinum
|
Rib
|
3.072
|
0.483
|
3.451
|
0.62
|
Ambystoma tigrinum
|
Ulna
|
2.483
|
0.483
|
2.483
|
0.62
|
Ambystoma tigrinum
|
Femur
|
2.689
|
0.483
|
2.914
|
0.62
|
Ambystoma tigrinum
|
Thoracic
|
2.762
|
0.483
|
3.017
|
0.62
|
Ambystoma tigrinum
|
Metatarsal
|
2.723
|
0.483
|
2.963
|
0.62
|
Ambystoma tigrinum
|
Humerus
|
2.347
|
0.483
|
2.436
|
0.62
|
Ambystoma tigrinum
|
Skull
|
2.622
|
0.483
|
2.822
|
0.62
|
Marmota monax
|
Tibia
|
0.651
|
0.224
|
0.191
|
0.288
|
Marmota monax
|
Caudal
|
1.083
|
0.224
|
0.795
|
0.288
|
Marmota monax
|
Ulna
|
0.846
|
0.224
|
0.464
|
0.288
|
Marmota monax
|
Fibula
|
0.690
|
0.224
|
0.244
|
0.288
|
Marmota monax
|
Rib
|
0.696
|
0.224
|
0.254
|
0.288
|
Marmota monax
|
Radius
|
0.937
|
0.224
|
0.590
|
0.288
|
Marmota monax
|
Femur
|
0.930
|
0.224
|
0.581
|
0.288
|
Marmota monax
|
Thoracic
|
0.930
|
0.224
|
0.581
|
0.288
|
Marmota monax
|
Metatarsal
|
1.167
|
0.224
|
0.912
|
0.288
|
Marmota monax
|
Humerus
|
1.122
|
0.224
|
0.849
|
0.288
|
Marmota monax
|
Skull
|
0.752
|
0.224
|
0.332
|
0.288
|
Columba livia
|
Tibia
|
0.446
|
0.254
|
0.493
|
0.325
|
Columba livia
|
Caudal
|
0.335
|
0.254
|
0.405
|
0.325
|
Columba livia
|
Radius
|
0.087
|
0.254
|
-0.010
|
0.325
|
Columba livia
|
Fibula
|
0.209
|
0.254
|
0.161
|
0.325
|
Columba livia
|
Rib
|
0.426
|
0.254
|
0.463
|
0.325
|
Columba livia
|
Ulna
|
0.133
|
0.254
|
0.054
|
0.325
|
Columba livia
|
Femur
|
0.070
|
0.254
|
-0.034
|
0.325
|
Columba livia
|
Thoracic
|
0.558
|
0.254
|
0.649
|
0.325
|
Columba livia
|
Metatarsal
|
0.056
|
0.254
|
-0.053
|
0.325
|
Columba livia
|
Humerus
|
0.129
|
0.254
|
0.049
|
0.325
|
Columba livia
|
Skull
|
0.422
|
0.254
|
0.459
|
0.325
|
Alligator sinensis
|
Tibia
|
0.549
|
0.135
|
0.275
|
0.171
|
Alligator sinensis
|
Caudal
|
0.075
|
0.135
|
-0.393
|
0.171
|
Alligator sinensis
|
Radius
|
0.685
|
0.135
|
0.460
|
0.171
|
Alligator sinensis
|
Fibula
|
0.504
|
0.135
|
0.212
|
0.171
|
Alligator sinensis
|
Rib
|
0.455
|
0.135
|
0.138
|
0.171
|
Alligator sinensis
|
Ulna
|
0.462
|
0.135
|
0.148
|
0.171
|
Alligator sinensis
|
Femur
|
0.915
|
0.135
|
0.782
|
0.171
|
Alligator sinensis
|
Thoracic
|
0.922
|
0.135
|
0.792
|
0.171
|
Alligator sinensis
|
Metatarsal
|
0.626
|
0.135
|
0.377
|
0.171
|
Alligator sinensis
|
Humerus
|
0.472
|
0.135
|
0.163
|
0.171
|
Alligator sinensis
|
Skull
|
0.790
|
0.135
|
0.607
|
0.171
|
Table S4: Due to the fact lacunae are 3-dimensional objects, the thin section preserve only one plane, and the most accurate representation of the size of that cell should be the largest 2-dimensional slice. Regression analysis of cell volume vs. genome size, however, does not demand an absolute, precise measurement of average cell volume for each taxon, but rather a consistent proxy for cell size in each taxon. This proxy, as long as it is consistently measured in all taxa, can then be regressed against genome size to obtain an equation to predict genome size from cell volume. Our study utilizes the size of all lacunae as a proxy for lacunar volume, because this method should be less prone to observer bias than arbitrarily averaging only a selected number of the largest osteocytes (as has been done in previous studies: Organ et al. 2007, 2009, 2011). However, to ease comparison with previous studies, we here present a supplementary analysis in which we use the average volume measurement of the largest 25% of osteocytes as a proxy for cell volume. The table below lists the average predicted genome size for each of the four taxa we studied, based on both the mean volume of all osteocyte lacunae and the mean volume of only the largest 25% of lacunae. As is shown, these estimates are statistically indistinguishable when standard deviation is taken into account, indicating that both methods of measuring cells return approximately the same estimated genome values.
Taxon
|
ln(Genome) est. 3rd Q
|
ln(Genome) est. mean
|
StDev
|
Alligator
|
0.49
|
0.26
|
0.13
|
Ambystoma
|
2.60
|
2.58
|
0.46
|
Columba
|
0.14
|
0.03
|
0.24
|
Marmota
|
1.12
|
0.79
|
0.21
|
Figure S1: An example of the intraspecific variation (95% CI) of reported C-values in the Animal Genome Size Database (www.genomesize.com). The variation in the database is most likely caused by different modes of measurement and associated error. Regardless of the source of this variation (measurement error or biological), the important result is that the spread of estimated genome size values (both the overall spread of estimates returned from all of the bones of a single taxon and the 95% CIs for the measurements of each individual bone) are still greater than the range of values presented in the database. In other words, estimated genome sizes using different bones are less precise than measured genome sizes. (Note: Not all species have multiple entries; the two plots below show our results for the pigeon and tiger salamander).
-
References
Abramoff, M.D., Magelhaes, P.J., Ram, S.J. 2004. Image processing with ImageJ.
Biophotonics International 11, 36-42.
Bruggeman, J., Heringa, J. & Brandt, B. W. 2009. PhyloPars: estimation of missing parameter values using phylogeny. Nucleic Acid Research 37, W179-W184.
Faure, E., Lony, E., Lovigny, R., Menegoz, A., Ting, Y. & Laurin, M. 2006 StratAdd module for Mesquite. (http://mesquiteproject.org/packages/stratigraphicTools/)
Maddison, W. P. & Maddison, D.R. 2010. Mesquite: a modular system for evolutionary analysis. Version 2.73 (http://mesquiteproject.org)
Organ, C.L., Canoville, A., Reisz, R.R. & Laurin, M. 2011. Paleogenomic data suggest
mammal-like genome size in the ancestral amniote and derived large genome size in
amphibians. Journal of Evolutionary Biology 24, 372-380.
Organ, C. L., Brusatte, S. L. & Stein, K. 2009. Sauropod dinosaurs evolved moderately sized genome unrelated to body size. Proceedings of the Royal Society B 276, 4303-4308.
Organ, C. L. & Shedlock, A. 2009. Palaeogenomics of pterosaurs and the evolution of small genome size in flying vertebrates. Biology Letters 5, 47-50.
Organ, C. L., Shedlock, A. M., Meade, A., Pagel, M. & Edwards, S. 2007. Origin of avian genome size and structure inc non-avian dinosaurs. Nature 446, 180-184.
|