Free Access
Issue
Ann. For. Sci.
Volume 67, Number 6, September 2010
Article Number 603
Number of page(s) 7
DOI https://doi.org/10.1051/forest/2010016
Published online 08 July 2010

© INRA, EDP Sciences, 2010

1. INTRODUCTION

The Spanish fir (Abies pinsapo Boiss.) is a conifer endemic to southern Spain and belongs to a Mediterranean basin group of firs, considered to be a relict species. In fact, it is an isolated fir because the nearest Abies populations are those of A. maroccana Trab. and A. tazaotana Villar in northern Morocco and A. alba Mill. in the Pyrenees.

The species is distributed in three isolated areas in “Sierra de Grazalema” Natural Park (N.P.) (Cadiz province) and “Sierra de las Nieves” N.P. and “Reales de Sierra Bermeja” Natural Reserve (N.R.) (Malaga province), declared protected zones. The most important stands are situated in “Sierra de las Nieves” N.P. with 2000 ha, followed by “Sierra de Grazalema” N.P. with 300 ha and “Reales de Sierra Bermeja” N.R. with 50 ha. In all the above areas A. pinsapo usually appears forming pure stands, although it is occasionally found with other tree species such as Quercus faginea Lam., Q. ilex L., Pinus halepensis Mill., P. pinaster Ait. and Juniperus oxycedrus L. (Arista and Talavera, 1997).

The Spanish fir forests represent one of the most unique ecosystems in the Iberian Peninsula due to their relict nature, and in order to ensure their conservation and protection, these sites have been declared “Biosphere Reserves” by UNESCO. However, there are several problems that affect the species as pests and diseases, low pollen production, climatic change and stand decay (Arista, 1995). Although there are some studies focused on the structure and dynamics of Spanish fir forests (Arista, 1995; Arista and Talavera, 1996; 1997), there appear to be no papers in the literature concerning the genetic structure of their populations. Specifically, most studies have been carried out to discriminate between the species of the Abies genus using different molecular markers as isozymes, mitochondrial DNA markers and chloroplast microsatellites (Pascual et al., 1993; Scaltsoyiannes et al., 1999; Terrab et al., 2007; Ziegenhagen et al., 2005).

A characteristic of gymnosperms that is useful for genetic studies is the haploid nature of megagametophyte. This tissue is derived from the female gametophyte that surrounds the embryo in a mature seed. Thus, analysis of the segregation pattern in parental trees may enable determination of the genetic control and linkage relationships without making crosses (El-Kassaby et al., 1987; Neale and Adams, 1981).

Genetic analysis of megagametophytes using isozymes has been carried out in conifers for the estimation of genetic diversity, heterozygosity, phylogeny and rate of gene flow in natural populations (Hamrick et al., 1992; Huang et al., 1994; Millar, 1983; Rogers, 1997; Wang and Liu, 1998; Wang and Nagasaka, 1997; Wheeler and Guries, 1982). Nevertheless, the role of seed storage proteins, that constitute an important component of the megagametophyte, has scarcely been studied. These proteins are, from a biological point of view, sources of amino acids for processes of synthesis occurring during germination. Nevertheless, the nature of these storage proteins is diverse, because different species can store these amino acids in four protein types (Gepts, 1990). In any case, their variations have no repercussion on any physiological function, which enables that these random changes to be maintained. Their main advantages as genetic markers are their high polymorphism, simple genetic control, environmental independence, and analysis efficiency.

Only few studies have been carried out on Gymnosperms megagametophyte storage proteins and in conifers they have mainly been based on biochemical characteristics of these proteins (Allona et al., 1994; Kurz et al., 1994; Piovesan et al., 1993). Likewise, our team used these proteins as markers to evaluate the genetic diversity of the stone pine (P. pinea L.), revealing a high level of variation in the studied populations which suggested that this technique could be useful to assess the genetic diversity in conifers (Alvarez et al., 2004).

In this study, we report for the first time, the inheritance of megagametophyte storage proteins in A.pinsapo, together with the use of these proteins as genetic marker for assessing the genetic diversity conserved in the extant stands of the species.

Table I

Description of the four populations of Abies pinsapo evaluated.

2. MATERIALS AND METHODS

2.1. Plant materials

Four populations were selected from the only three areas where this species is distributed. Samples were taken from a total of thirty trees, varying in number from six to nine per population (Tab. I). The sampling sites and the number of trees per population were selected according to the area covered by the species in each site and the heterogeneity of the stands. Further, the trees selected were a minimum of 50 m apart to avoid the sampling of closely related trees.

At least three cones were collected per tree given the low fertility expected for this species. Indeed, the number of viable seeds did not exceed 10% of the total collected per tree. For the analysis, ten megagametophytes were used per tree. According to Morris and Spieth (1978), this is an optimum number of female megagametophytes per tree for estimating polymorphism in conifers.

thumbnail Figure 1

SDS-PAGE of representative samples of the variation found in A. pinsapo. B1-D5 are the polymorphic bands found in the study. In square are shown the two allelism cases, C3 vs. C4 and D4 vs. D5.

2.2. Protein extraction and electrophoretic analysis

Megagametophytes were isolated by removing both the seed coat and the embryo. The lipids were removed with diethyl ether (50:1 v/w, 1 h, 4 °C), followed by a second extraction with acetone (50:1 v/w, 1 h, 4 °C). Albumin were extracted according the protocol proposed by Fonseca et al. (1997) for Q. suber L. and modified by Alvarez et al. (2004) for P. pinea. 500 μL of 10 mM Tris-HCl pH: 7.5 + 0.1% (w/v) dithiothreitol + 10 mM MgCl2+ 10 mM CaCl2 were added to each sample and incubated for 2 h at 4 °C. The samples were centrifuged at 14.000 g for 10 min and the supernatant containing the albumin fraction was transferred to a new tube. The samples were precipitated with 1 mL of cold-acetone, and the dried pellets were solubilised in buffer containing 625 mM Tris-HCl pH: 6.8, 2% (w/v) SDS, 10% (v/v) glycerol, 0.02% (w/v) bromophenol blue, and 2% (w/v) dithiothreitol in ratio 1:5 (w/v).

Proteins were fractionated in vertical SDS-PAGE slabs in a discontinuous Tris-HCl-SDS buffer system (pH: 6.8/8.8) at a 12% polyacrylamide concentration (w/v, C = 1.28%). The Tris-HCl/glycine buffer system of Laemmli (1970) was used. Electrophoresis was performed at a constant current of 30 mA/gel at 18 °C for 10 min after the tracking dye migrated off the gel. Gels were stained overnight with 12% (w/v) trichloroacetic acid solution containing 5% (v/v) ethanol and 0.05% (w/v) Coomassie Brilliant Blue R-250. De-staining was carried out with tap water.

2.3. Genetic analysis

To confirm Mendelian inheritance of polymorphic bands, the offspring of two heterozygous trees were analysed. Segregation of all polymorphic bands was tested for goodness of fit to the expected 1:1 Mendelian segregation ratio using chi-square (χ2) test. Similarly, linkage between bands was also tested. The linkage distance in cM was estimated by the recombination frequency (r = f/n, where f is the number of recombinant types observed in n megagametophytes).

The following genetic variability parameters were calculated in all the populations: number of alleles per locus (A), percentage of polymorphic bands (P), effective number of alleles per locus (He) and expected heterozygosity (He). The total genetic diversity over all populations (Ht) and the average genetic diversities within (Hs) and between (Dst) populations were also calculated according to Nei (1973). The relative magnitude of genetic differentiation between populations, Gst, was estimated as Dst/Ht. Popgene software version 1.32 (Yeh et al., 1997) was used for the estimation of the above parameters.

The genetic identity (I) values were estimated between populations (Nei, 1972). To visualize the genetic relationship among populations, a dendrogram was constructed based on Nei’s genetic distances (D = − lnI) by an unweighted pair-group method of cluster analysis using arithmetic averages (UPGMA) (Sneath and Sokal, 1973). The significance of the phenogram was tested by cophenetic correlation coefficient (Rohlf and Fisher, 1986).

3. RESULTS

3.1. Albumin differentiation

The sequential extraction process resulted in four possible protein fractions: albumins, globulins, prolamins and glutelins. Albumin was the fraction which showed the highest and the clearest resolution between bands with absence of ghost bands and other artefacts. Moreover, these bands were consistent, being the procedure higher repetitive.

Table II

Observed segregation of albumins in two trees. Chi-square test and goodness of fit to the 1:1 ratio.

According to the molecular weight marker used as reference, five zones were recorded in the gel (Fig. 1), classified as A (upper 66 KDa), B (66-45 KDa), C (45-30 KDa), D (30-14 KDa) and E (below 14 KDa). In total, up to 41 different bands were consistently scored, 14 of which were polymorphic (34%). No polymorphic bands were detected in zones A and E, whereas in the three others (zones B, C and D), five, four and five polymorphic bands were identified respectively (Fig. 1). The bands were named by the letter corresponding to the zone where they were situated and a corresponding number as a function of their mobility (Fig. 1).

The presence or absence of each polymorphic band was fit to a 1:1 ratio, and no significant deviation from the expected segregation was detected in any case (Tab. II). Of the 48 possible pairs of band combinations two allelic relationships were detected between C3 vs. C4 and D4 vs. D5 (Tab. II). In fact, in both cases each seed displayed either one or other of the two bands, but no simultaneous presence or absence was detected in any case (Fig. 1). Further, 38 cases of independence (fit a 1:1:1:1 segregation ratio) (Tab. III) and 8 cases of linkage were detected (Tab. IV).

Table III

Observed segregation using pairs of bands. Chi-square test, goodness of fit to the independence 1:1:1:1 ratio.

Table IV

Observed segregation using pairs of bands. Chi-square test, goodness of fit to the 1:1:1:1 ratio and linkage distance (d).

At population level, from all the studied samples, the bands that showed the highest frequencies were B3, C4 and D4, appearing in more than 78% of cases. The bands with the lowest frequencies, all situated in zone D (D3, D5 and D2, respectively), were present in less than 17% of the individuals (Tab. V). In general, bands displayed a wide distribution among the populations. The identified polymorphic bands were present in all populations, except band C4 that was monomorphic in the SG-1 population (Tab. V).

3.2. Variation within and between populations

All populations showed high polymorphism for the assessed bands. Populations SN-1, SN-2 and RSB-1 had the highest levels of polymorphic loci while this value was the lowest in population SG-1, where band C4 was monomorphic (Tab. VI). The mean expected heterozygosity (He) was 0.379, ranging from He = 0.412 for the SN-2 population and He = 0.322 in RSB-1 (Tab. VI).

In the four populations, approximately half of the genetic diversity was among trees within populations and the other half among seeds within trees (Tab. VII). The total genetic diversity (Ht) in the evaluated populations was 0.377, which was divided into two components, Hs and Dst, to assess the genetic diversity within and between populations. The total genetic differentiation among populations (Gst) was 5.8%, which indicates that 94.2% of the variation was within populations. Further, the highest differentiation among population pairs was found with RSB-1 and SG-1 (Gst = 7.0%) whereas the two populations from “Sierra de las Nieves” N.P. showed the lowest value (Gst = 1.2%) (Tab. VII).

Table V

Frequencies of each band in the 297 seeds analysed and in the 4 populations evaluated.

thumbnail Figure 2

Dendrogram of A. pinsapo populations based on Nei’s genetic distance. (Co-phenetic correlation = 0.852***).

3.3. Genetic distances

Genetic identity coefficients (I) were estimated for paired comparisons of all four populations based on the normalized identity of all polymorphic bands between each pair of populations. The mean value was I = 0.943, ranging from I = 0.926 for the pair SG-1/RSB-1 to I = 0.986 for SN-1/SN-2.

A dendrogram, based on Nei’s genetic distance matrix, was tested for the significance of the clustering method. The cophenetic correlation coefficient (r = 0.852***) indicated a good fit between the cluster and the data. Cluster analysis based on genetic distances, revealed that populations displayed a pattern of geographical variation among them (Fig. 2).

Table VI

Albumin genetic diversity based on the 14 polymorphic bands identified in the four Spanish fir populations.

4. DISCUSSION

Our data indicate that albumin was the most appropriate fraction to assess genetic diversity in A. pinsapo, due to the high polymorphism, consistent repeatability and simple genetic control. Nevertheless, this result contrasts with that obtained by other authors that indicate the globulin as the best fraction for genetic analysis in other conifers (Allona et al., 1994; Alvarez et al., 2004).

The segregation analysis of megagametophytes demonstrated the genetic control in polymorphic bands, indicating that they exhibit Mendelian inheritance. This result agrees with that obtained by Pascual et al. (1993) using isozymes in A. pinsapo. This being the first study of inheritance of these proteins, the two cases of allelism and the eight cases of linkage, constitute the beginning of the cataloguing of seed storage protein alleles in this species. As well as these two loci (allelic bands), ten more loci were detected, the majority being independent. This high percentage of independent loci provides more efficient and accurate estimates for some population analysis such as genetic distance and matting systems (Nei, 1975; Shaw and Allard, 1981).

Table VII

Differentiation of albumin diversity within and between 4 Spanish fir populations, and comparison between pairs of populations.

According to the results obtained for albumins, the number of polymorphic bands was 34.0%. The bands were very variable with a mean of polymorphic loci of 98.2%. This value was significantly higher than that found in the same species using isozymes (12.5%, Scaltsoyiannes et al. 1999) and similar to that detected in P. pinea (97.0%, Alvarez et al., 2004).

The mean genetic diversity (He = 0.379) was higher than the detected in the same species (He = 0.010) using isozymes (Scaltsoyiannes et al., 1999) but significantly lower than that found by Terrab et al. (2007) using chloroplast microsatellites (He = 0.833). Moreover, this value was higher than that reported for other species of the genus Abies as A. alba (0.149), A. cephalonica (He = 0.161) and A. nephrolepis (0.240) or found in other conifers (He = 0.173; Hamrick and Godt, 1990; Woo et al., 2008). However, most of the available population genetic data in the species are based on other markers, which may be the reason why the values are not comparable with our results.

In general, widely distributed Abies species have higher values of genetic diversity and lower levels of genetic differentiation than species with restricted distribution (Diebel and Feret, 1991; Fady and Conkle, 1993). Conversely, genetic differentiation detected among populations (Gst = 5.8%) was lower than the average value for gymnosperm (Gst = 18%) and woody species with an endemic distribution (Gst = 14.1%; Hamrick et al., 1992). The results agree with Terrab et al. (2007) who found a low differentiation among populations of the species (Fst = 0.04). This low value seems to be typical for out-crossing wind-pollinated species (Hamrik, 1989) and could indicate that gene flow through pollen among populations is fairly high and there are no strong barriers to gene exchange. In fact, the populations examined in the study are from a restricted area in southern Spain, separated by only a few kilometres. Indeed, the two populations from the same site (Sierra de las Nieves), were much more similar to each other than they were to the other populations, even though this area represents 85% of the whole area covered by this relict species.

In conclusion, albumins have been shown to be a good marker of the genetic diversity in A. pinsapo due to their high polymorphism, neutrality, consistent repeatability and simple genetic control. The preliminary analysis has highlighted a considerable genetic diversity in the species and low differentiation between populations. Moreover, the results suggest that increasing the number of locations could give interesting additional information on the spatial distribution of the genetic diversity in this species. Accordingly, the study of trees of different ages within locations would be very useful to determine the evolution between these populations and the increase or decrease in the genetic diversity as a consequence of the process of inbreeding. This knowledge could add complementary information to that currently available and thereby play a significant role in the formulation of appropriate conservation management strategies for the species.

Acknowledgments

The authors thank J. Lopez Quintanilla (Natural Park of “Sierra de las Nieves”, Malaga, Spain) for the supply of the seeds of Spanish fir.

References

  • Allona I., Collada C., Casado R. and Aragoncillo C., 1994. Electrophoretic analysis of seed storage proteins from gymnosperms. Electrophoresis 15: 1062–1067. [CrossRef] [PubMed] [Google Scholar]
  • Alvarez J.B., Toledo M.J., Abellanas B. and Martín L.M., 2004. Use of megagametophyte storage proteins as markers of the genetic diversity in stone pine (Pinus pinea L) in Andalusia, Spain. Genet. Resour. Crop Evol. 51: 621–627. [CrossRef] [Google Scholar]
  • Arista M., 1995. The structure and dynamics of an Abies pinsapo forest in southern Spain. For. Ecol. Manage. 74: 81–89. [CrossRef] [Google Scholar]
  • Arista M. and Talavera S., 1996. Density effect on the fruit-set, seed crop viability and seedling vigour of Abies pinsapo. Ann. Bot. 77: 187–189. [CrossRef] [Google Scholar]
  • Arista M. and Talavera S., 1997. Gender expression in Abies pinsapo Boiss., a Mediterranean fir. Ann. Bot. 79: 337–342. [CrossRef] [Google Scholar]
  • Diebel K.E. and Feret P.P., 1991. Isozyme variation within the Fraser fir (Abies fraseri (Pursh) Poir.) population on Mount Rogers, Virginia: Lack of microgeographic differentiation. Silvae Genet. 40: 79–85. [Google Scholar]
  • El-Kassaby Y.A., Meagher M.D., Parkinson J. and Portlock F.T., 1987. Allozyme inheritance, heterozygosity and outcrossing rate among Pinus monticola near Ladysmith, British Columbia. Heredity 58: 173–181. [CrossRef] [Google Scholar]
  • Fady B. and Conckle M.T., 1993. Allozyme and possible phylogenetic implications in Abies cephalonica Loudon and some related eastern Mediterranean firs. Silvae Genet. 42: 351–359. [Google Scholar]
  • Fonseca P.A., Ferreira R.B. and Texeira A.R., 1997. Seed proteins from Quercus suber. J. Agric. Food Chem. 45: 3443–3447. [CrossRef] [Google Scholar]
  • Gepts P., 1990. Genetic diversity of seed storage proteins in plants. In: Brown A.H.D., Clegg M.T., Kahler A.L., and Weir B.S. (Eds.), Plant population genetics, breeding and genetic resources. Sinauer Associates Inc Publishers, Suderland, Massachusetts, pp. 64–82. [Google Scholar]
  • Hamrick J.L., 1989. Isozymnes and the analysis of genetic structure in plant populations. In: Soltis D.E., and Soltis P.S. (Eds.), Isozymes in plant biology, Portland, Oregon, pp. 87–105. [Google Scholar]
  • Hamrick J.L. and Godt M.J.W., 1990. Allozyme diversity in plant species. In: Brown H.D., Clegg M.T., Khaler A.L., and Weir B.S. (Eds.), Plant Population Genetics, Breeding, and Genetic Resources. Sinauer Press, Sunderland, MA, pp. 43–63. [Google Scholar]
  • Hamrick J.L., Godt M.J.W. and Sherman-Broyles S.L., 1992. Factors influencing levels of genetic diversity in woody plant species. New For. 6: 95–124. [Google Scholar]
  • Huang Q.Q., Tomaru N., Wang L.H. and Ohba K., 1994. Genetic-control of isozyme variation in masson pine, Pinus massoniana Lamb. Silvae Genet. 43: 285–292. [Google Scholar]
  • Kurz M.L., Roberts D.R., Flinn B.S. and Vidaver W.E., 1994. A morphological and biochemical analysis of embryo maturation during yellow cypress seed development. Can. J. For. Res. 24: 431–441. [CrossRef] [Google Scholar]
  • Laemmli U.K., 1970. Cleavage of structural proteins during the assembly of the head of bacteriophage T4. Nature 227: 680–685. [CrossRef] [PubMed] [Google Scholar]
  • Millar C.I., 1983. A steep cline in Pinus muricata. Evolution 37: 311–319. [CrossRef] [PubMed] [Google Scholar]
  • Morris R.W. and Spieth P.T., 1978. Sampling strategies for using female gametophytes to estimate heterozygosity in conifers. Theor. Appl. Genet. 51: 217–222. [CrossRef] [PubMed] [Google Scholar]
  • Neale D. and Adams W.T., 1981. Inheritance of isozyme variants in seed tissues of balsam fir (Abies balsamea). Can. J. Bot. 59: 1285–1291. [CrossRef] [Google Scholar]
  • Nei M., 1972. Genetic distances between different populations. Am. Nat. 106: 283–292. [CrossRef] [Google Scholar]
  • Nei M., 1973. Analysis of gene diversity in subdivided populations. Proc. Natl. Acad. Sci. USA 70: 3321–3323. [Google Scholar]
  • Nei M., 1975. Molecular population genetics and evolution. Amsterdam, North-Holland. [Google Scholar]
  • Pascual L., Garcia F.J. and Perfectti F., 1993. Inheritance of isozyme variation in seed tissues of Abies pinsapo Boiss. Silvae Genet. 42: 335–340. [Google Scholar]
  • Piovesan G., Pelosi C., Schirone A. and Schirone B., 1993. Taxonomic evaluations of the genus Pinus (Pinaceae) based on electrophoretic data of salt soluble and insoluble seed storage proteins. Plant. Syst. Evol. 186: 57–68. [CrossRef] [Google Scholar]
  • Rogers D.L., 1997. Inheritance of allozymes from seed tissues of the hexaploid gymnosperm, Sequoia sempervirens (D. Don) Endl (Coast redwood). Heredity 78: 166–175. [CrossRef] [Google Scholar]
  • Rohlf F.J. and Fisher D.L., 1986. Test for hierarchical structure in random data sets. Syst. Zool. 17: 407–412. [CrossRef] [Google Scholar]
  • Scaltsoyiannes A., Tsaktsira M. and Drouzas A.D., 1999. Allozyme differentiation in the Mediterranean firs (Abies, Pinaceae). A first comparative study with phylogenetic implications. Plant Syst. Evol. 216: 289–307. [Google Scholar]
  • Shaw D.V. and Allard R.W., 1981. Analysis of mating system parameters and population structure in Douglas Fir using single-locus and multi-locus methods. In: Conkle M.T., and Berkeley C.A. (Eds.), Proc. Symp. Isozymes North Am. Forest Trees and Forest Insects, USDA Forest Serv. Gen. Techn. Rep. PSW-48, pp. 18–22. [Google Scholar]
  • Sneath P.H.A. and Sokal R.R., 1973. Numerical Taxonomy. San Francisco, Freedman. [Google Scholar]
  • Terrab A., Talavera S., Arista M., Paun O., Stuessy T.F. and Tremetsberger K., 2007. Genetic diversity at chloroplast microsatellites (cpSSR) and geographic structure in endangered West mediterranean firs (Abies spp., Pinaceae). Taxon 56: 409–416. [Google Scholar]
  • Wang C.T. and Liu T.P., 1998. Inheritance and linkage relationships of allozymes, and estimation of outcrossing rates in a seed orchard of Cunninghamia konishii Hay. Silvae Genet. 47: 33–37. [Google Scholar]
  • Wang Z.M. and Nagasaka K., 1997. Allozyme variation in natural populations of Picea glehnii in Hokkaido, Japan. Heredity 78: 470–475. [CrossRef] [Google Scholar]
  • Wheeler N.C. and Guries R.P., 1982. Population structure, genetic diversity, and morphological variation in Pinus contorta Dougl. Can. J. For. Res. 12: 595–606. [CrossRef] [Google Scholar]
  • Woo L.S., Hoon Y.B., Don H.S., Ho S.J. and Joo L.J., 2008. Ann. For. Sci. 65: 302–308. [Google Scholar]
  • Yeh F.C., Yang R.C., Boyle T.B.J., Ye Z.H., and Mao J.X., 1997. Popgene ver 1.32. The user-friendly software for population genetic analysis. Molecular Biology and Biotechnology Center, University of Alberta, Canada. [Google Scholar]
  • Ziegenhagen B., Fady B., Kuhlenkamp V. and Liepelt S., 2005. Differentiating groups of Abies species with a simple molecular marker. Silvae Genet 54: 123–126. [Google Scholar]

All Tables

Table I

Description of the four populations of Abies pinsapo evaluated.

Table II

Observed segregation of albumins in two trees. Chi-square test and goodness of fit to the 1:1 ratio.

Table III

Observed segregation using pairs of bands. Chi-square test, goodness of fit to the independence 1:1:1:1 ratio.

Table IV

Observed segregation using pairs of bands. Chi-square test, goodness of fit to the 1:1:1:1 ratio and linkage distance (d).

Table V

Frequencies of each band in the 297 seeds analysed and in the 4 populations evaluated.

Table VI

Albumin genetic diversity based on the 14 polymorphic bands identified in the four Spanish fir populations.

Table VII

Differentiation of albumin diversity within and between 4 Spanish fir populations, and comparison between pairs of populations.

All Figures

thumbnail Figure 1

SDS-PAGE of representative samples of the variation found in A. pinsapo. B1-D5 are the polymorphic bands found in the study. In square are shown the two allelism cases, C3 vs. C4 and D4 vs. D5.

In the text
thumbnail Figure 2

Dendrogram of A. pinsapo populations based on Nei’s genetic distance. (Co-phenetic correlation = 0.852***).

In the text