- Download PDF-
SEEFOR 6 (1): 39-51
Article ID: 43
Original scientific paper
Leaf Morphology Variation of Populus nigra L. in Natural Populations along the Rivers in Croatia and Bosnia and Herzegovina
Davorin Kajba 1*, Dalibor Ballian 2, Marilena Idžojtić 1, Igor Poljak 1
1 University of Zagreb, Faculty of Forestry, Department of Forest Genetic, Dendrology and Botany, Svetošimunska 25, HR-10000 Zagreb, Croatia
2 University of Sarajevo, Faculty of Forestry, Department of Silviculture and Urban Greening, Zagrebačka 20, BA-71000 Sarajevo, Bosnia and Herzegovina
KAJBA D, BALLIAN D, IDŽOJTIĆ M, POLJAK I 2015 Leaf Morphology Variation of Populus nigra L. in Natural Populations along the Rivers in Croatia and Bosnia and Herzegovina. South-east Eur for 6 (1): 39-51. DOI: http://dx.doi.org/10.15177/seefor.15-06
Received: 10 Jan 2015; Accepted: 13 Mar 2015; Published online: 29 Mar 2015
Background and Purpose: The aim was to determine the morphological differences between the hairy type of European black poplar (Populus nigra subsp. caudina) and the typical type from the riparian forests populations as well as between the river systems. Hairy black poplar spreads in a mosaic pattern across the Submediterranean climatic type along the River Neretva and the typical European black poplar is growing on alluvial soils along large rivers in the territory of Croatia and Bosnia and Herzegovina.
Material and Methods: Samples for leaf morphometric analysis were collected in 17 natural populations of European black poplar along six rivers in Croatia and Bosnia and Herzegovina.
Results: Discriminant analyses have determined that in the differentiation of population groups largely contribute some characters such as the distance between the leaf widest part and the leaf base (DBW) and the petiole length (PL). The differences between populations and analysed groups, as well as the differences between populations belonging to a particular river system, were confirmed for all studied characteristics.
Conclusions: Significant differences have been determined between the typical and the hairy type of European black poplar in the studied morphological traits and these dissimilarities are in accordance with the climatic differences in respective habitats of continental riparian forests and the Submediterranean type of climate. Populations sampled in the lower course of the River Neretva, which correspond to the hairy type of the European black poplar, have smaller leaves and a greater angle between the first lower lateral vein and the midrib.
Keywords: leaf morphological traits, European black poplar, Populus nigra subsp. caudina, genetic variability
The conservation of genetic resources of some forest tree species, particularly the endangered ones, should have high priority in any large-scale programme of forest ecosystem management. The European black poplar (Populus nigra L.) is a species growing on alluvial soils along large rivers in the temperate zones in Europe, Asia and Northern Africa. In the territory of the Republic of Croatia, the European black poplar grows in riparian floodplain forests along rivers Mura, Sava, Drava and the Danube. Its presence in the total forest area of Croatia is small, yet their universal role is significant. Because of continuous human activities in the European black poplar’s natural habitats, its range has been reduced to individual, partly separated areas only, although for the European conditions the habitats are very well preserved riparian forests . The overall length of rivers and creeks in Bosnia and Herzegovina is around 11000 km, but the habitats of European black poplar are limited only to the valleys of bigger rivers, such as Drina, Bosna, Una, Sana and Neretva with the total length of 1200 km . In Herzegovina, along the River Neretva, a special subspecies of black poplar can be found, known as the hairy poplar (P. nigra subsp. caudina), which is characterised by distinctly hairy shoots and leaves . Hairiness retains on one-year and two-year-old shoots and partially on older shoots as well. Apart from the leaves and shoots, the petioles of inflorescence and fruits are also permanently hairy. So far, these poplars have been recorded in two more southern Balkan areas, in southern Macedonia and in Thessaly in central Greece, while some populations from Albania and Turkey have also attributed to this taxon , , . The studied population of the hairy type of black poplar is found in the Submediterranean type of climate. This climate is characterised by mild winters and dry summers. The value of gene conservation in these populations may be crucial for the whole species with regard to the trends of climate changes , , , . The survival of the European black poplar, as an autochthonous species, has been endangered also by the cultivation of superior hybrids of Euroamerican poplars (P. × canadensis) and Eastern Cottonwood (P. deltoides) clones within the area of its natural range. Some of these imported clones were female and their flowering was synchronized with the European black poplar. Crossing between these various poplar species is possible and the question is whether in the European black poplar juvenile populations there is the introgression of the Eastern Cottonwood genes or the Euroamerican hybrid genes , , . The genetic resources of the European black poplar are also threaten by a non-planned land use along the rivers, the regulation of river flows, gravel pits, dumps, and the lack of legislative regulations.
In this study the intention is, apart from hairiness, to find out which dissimilarities can be found in other morphological traits of the leaves of the typical and hairy black poplar, as well as between the populations and the rivers’ systems. It also examines the importance of the hairy type of the European black poplar concerning conservation and the breeding of poplars with regard to the increasing effect of climate changes.
MATERIALS AND METHOD
The samples for leaf morphometric analysis were collected in 17 natural populations of the European black poplar in the area of Croatia and Bosnia and Herzegovina (Figure 1). In Croatia, the investigations were conducted among five populations along the basin of the River Drava (Žirovnjak, Tršćana, Bobrovac, Šućurica and Topolje) and three populations along the River Sava (Jarun, Zaprešić 1, Zaprešić 2) from 5 to 15 years old. In Bosnia and Herzegovina the samples were collected in nine populations from four rivers: one population was along the River Drina (Tegare), two populations were collected along the River Vrbas and the River Bosna (Zlavast, Banja Luka, Bilješevo, Doboj) and four populations were sampled along the River Neretva (Ovčari, Počitelj, Čapljina, Metković), all being at the age up to mature trees (30 years). To simplify the processing of data, the studied populations were divided into four separate groups according to the ecogeographical zones. Groups CRO-01 and CRO-02 belong to Croatian rivers (Drava and Sava). Group B&H-01 is com-posed of six populations from the inner part of Bosnia and Herzegovina and group B&H-02 consists of three populations that were sampled in the lower course of the River Neretva in the Submediterranean region. Leaf samples were collected from short fertile shoots and in total the analysis included 514 trees and 6513 leaves, from which 32565 data have been measured and statistically processed (Table 1).
FIGURE 1. The location of the analyzed Populus nigra populations (after Vanden Broeck 2003)
TABLE 1. Sampled populations
The fully developed, undamaged European black poplar leaves have been investigated. Within one tree the leaf dimorphism and the seasonal heterophylly in poplars was also found , . The leaves were taken only from short shoots only in the central part of tree crown, which proved to give a more reliable description of the studied morphological properties than the leaves from the long shoots , , .
The measurements of the following morpho-logical parameters have been taken: leaf blade length (LBL), leaf blade width (LBW), petiole length (PL), the angle between the first lower lateral vein and the midrib (α) and the distance between the leaf widest part and the leaf base (DBW), as shown in Figure 2.
FIGURE 2. The investigated leaf morphological traits (LBL - leaf blade length, LBW - leaf blade width, DBW - distance between the leaf widest part and the leaf base, PL - petiole length, α - angle between the first lower lateral vein and the midrib).
Measured morphological traits were summarized using descriptive statistical parameters, with standard algorithms using descriptive statistical analysis. Data are presented following univariate statistical parameters: arithmetical mean (M), standard deviation (SD) and the coefficient of variation (CV). Statistically significant differences between the investigated objects were determined using univariate nonparametric Kruskal-Wallis analysis of variance tests (Kruskal-Wallis ANOVA).
The analysis was carried out on the entire sample in a manner that examined differences between all populations, among the groups CRO-01, CRO-02, B&H-01 and B&H-02, between the groups CRO-01 and CRO-02, groups B&H-01 and B&H-02 and between the populations on the basis of belonging to a particular river system.
According to McGarigal et al., multi-variate statistical methods were used (the principal components analysis and the discri-minant analysis) to show the similarities or the differences between the studied populations and groups on the basis of the measured morphological characteristics of leaves. In both analyses the arithmetic means of each population for each parameter were used. Five principal components in total were derived by principal component analysis. The first two of them were selected by Kaiser criterion (eigenvalue greater than 1) for further analysis of relations between populations . When it comes to the interpretation of the results, the first two principal components are shown in the graph to determine the trend of the relationship between the populations. To determine the features that best distinguish predefined groups, a discriminant analysis was carried out. These statistical analyses were performed using the statistical program STATISTICA 8.0 .
The results of the descriptive statistical analysis in Table 2 are shown according to the populations, and in Table 3 are listed according to the investigated groups. The results show that the Croatian population from the group CRO-01 has the highest average values for the traits such as leaf blade length (LBL), leaf blade width (LBW), the distance between the leaf widest part and the leaf base (DBW) and the petiole length (PL). On the contrary, the highest average values for the trait angle between the first lower lateral vein and the midrib (α) were peculiar to populations from the group CRO-2. The two groups of studied populations from Bosnia and Herzegovina had larger and wider leaf blades (LBL, LBW) and longer petioles (PL) that characterized populations from the inner part of Bosnia and Herzegovina (group B&H-01). Populations which were sampled in the lower course of the River Neretva (group of B&H-02), except that they had smaller leaves, also had the higher values for the angle between the first lower lateral vein and the midrib (α). Also it can be seen that the highest degree of variability shows the trait of petiole length (PL) as it is shown in Tables 2 and 3. Higher variability for all measured characters for Croatian river populations was determined for the group CRO-02 (the River Sava) and in Bosnia and Herzegovina for the group B&H-02 (lower part of the River Neretva), except for the trait angle between the first lower lateral vein and the midrib (α), which shows a slightly greater degree of variability.
TABLE 2. Descriptive statistics of the measured morphological traits at the population level
TABLE 3. Descriptive statistics of the measured morphological traits at the group level
The differences between populations and analysed groups, as well as the differences between the populations belonging to a particular river system, were confirmed for all studied characteristics. Table 4 shows the results of non-parametric univariate analysis of variance (Kruskal-Wallis ANOVA). Analysed populations and groups, as well as the populations belonging to a particular river system, differ from each other with statistical significance level of 0.0001 for all the measured characters of leaves.
TABLE 4. Kruskal-Wallis ANOVA for the measured morphological traits
The results of the principal component analysis show that the cumulative variability captured for the first two PC-axes is approximately 90%, with the largest share of the first PC-axis of approximately 60% (Table 5). The remaining three axes, whose eigenvalue is less than 1 in total variability, contribute to only about 10% and can be ignored. The contribution of each of the measured morphological characters to the each of the calculated PC-axis is shown in Table 6. Figure 3 shows the relation between studied populations based on the first and second principal components. From the Table 6 we can assume that all of the original variables are positively correlated with the first principal component, regarding that from the right side of the diagram populations with larger leaves are grouped and on the left side populations with smaller leaves. The distribution of the studied populations along the second PC-axis largely contributes to a positive correlation for the trait angle between the first lower lateral vein and the midrib (α), and a negative correlation for the character distance between the leaf widest part and the leaf base (DBW).
TABLE 5. Eigenvalues, the percent of variance and cumulative variance
TABLE 6. Factor-variable correlations (factor loadings)
FIGURE 3. Scatterplot of the PC analysis for 17 studied populations of Populus nigra. Particular values for the first PC are on the x axis, and for the second PC on the y axis.
The discriminant analysis was used to de-termine which leaf morphological characters best separate the studied groups. The results of the discriminant analysis suggest that the variability between the experimental groups is significant (Wilk’s λ = 0.00032; F (15, 25) = 29.336; p<0.00001). For five variables and four groups of canonical analysis, three discriminant functions were obtained. From the means of the canonical variables it is evident that the discriminant function 1 in the fullest extent particularly distinguishes group B&H-01 from other studied groups (Figure 4, Table 7). In the Table 8 it is evident that in this differentiation largely contribute characters such as the distance between the leaf widest part and the leaf base (DBW) and petiole length (PL). Since this table contains the eigenvalue and the cumulative proportion of the explained variation for each discriminant function, it can be determined that the discriminant function 1 contains 68.2% of the explained variability. The discriminant function 2, which together with discriminant function 1 explains 95.9% of the variability, clearly distinguished group of B&H-02 from the other three studied populations which were sampled on the River Sava and which belong to the group CRO-02 (Figure 4, Table 7). Variables α (the angle between the first lower lateral vein and the midrib) and LBL (leaf blade length) largely contribute in the differentiation between these two groups, as it is shown in Table 8. The third discriminant function, which explains only 4.1% of the total variability, best distinguished the two groups of Croatian population belonging to the basin of the River Drava and the River Sava (Table 7).
FIGURE 4. Scatterplot of the canonical scores for four studied groups of Populus nigra. Particular values for the first discriminant function are on the x axis, and for the second function on the y axis
TABLE 7. Means of canonical variables
TABLE 8. Standardised coefficients for canonical variables
Interpopulation variability and the variation between the river systems were statistically significant for all the investigated morphological traits that correspond to the previous research in the juvenile populations of the European black poplar along the River Drava and the River Sava, , . Some of the morphological parameters, such as petiole length (PL), show a statistically important difference between the investigated populations, and have also been determined in previous investigations . Interpopulation variability in the Danube basin was in most of morphological characters between 10 and 20% . The highest contribution of variation between populations had a number of veins on the left and the right blade side. The highest contribution of variation within the populations (more than 50%) had a distance between the base of the leaf blade and the widest part . This character corresponds to the results of the interpopulation differentiation in our study. The mentioned characters is under a genetic control and is less influenced by the environmental changes, so it could serve to estimate the intrapopulation and interpopulation variability, as well as to show the introgression of the genes of other poplar species in local population of the European black poplar. The morphological variation of seven leaf characters from 84 trees of eastern cottonwood (Populus deltoides Bartr. ex. Marsh. var. deltoides) from nine populations in Ontario was also studied . Leaf characters showed significant variation among populations (P≤0.003) and the number of veins showed a significant correlation with the longitude (r=0.81, P<0.009) and the latitude (r=-0.69, P<0.05). None of the other leaf characters were related to latitude and longitude.
The genetic diversity of ten black poplar populations along the main river systems across Slovenia and Croatia (rivers Soča, Sava, Drava and Mura) was revealed by using six microsatellite loci. The results of this study indicate that the gene pool of the remaining Populus nigra populations in this territory is still well preserved. In natural populations, the introgression of the genes of Populus deltoides was very low. One of the populations along the River Drava (Suhopolje, near locality Šućurica, population Drava 4 in Table 1 and 2) was distinct from the others and clustered alone as a separate group.
The assessment of the population genetic structure of the Populus nigra L. along the rivers revealed a significant spatial genetic structure which was probably caused by a limited gene flow. This helps to explain the reduced diversity in the juveniles, , .
The results of the research of the European black poplar leaf morphological variability in the natural stands along the Sava and Drava show that in the studied progeny the phenotypes resemble each other and mostly prevail among the European black poplar, , with no introgression with P. deltoides or poplar hybrids. These results confirm the researches of Van den Broeck et al. , carried out in Belgium. By studying the poplar flowering phenology in clonal archives they noted that the male poplar clone (P. nigra var. italica) and 22 clones of P. × canadensis flower earlier than the studied clones of local European black poplar (P. nigra). Based on this finding, they concluded that there is no great danger of the introgression of the American black poplar genes into the European black poplar genom, but also that the introgression of genes should not be excluded completely, having in mind that considerable variations in the poplar flowering phenology have been noticed.
Leaf morphometric characters have sig-nificantly contributed in the intra- and interpopulation variability of Populus nigra populations in the basin of the River Danube. The discriminative model with all characters achieved less than 60.38% of correct allocation due to a high variability of examined characters within populations. The results of morphological analysis do not entirely support the results of the molecular analysis, but they generally support the similarity among neighbouring populations and, according to this, the results provide a proper selection and the use of forest reproductive material.
Hairy subspecies of European black poplar (Populus nigra subsp. caudina) spreads in a mosaic pattern across a warm and dry area of Herzegovina along the River Neretva, while along the rivers in the region of central and mountainous Bosnian belt only individual trees are mixed with the typical black poplar. It can be assumed that the area of its current distribution range is what has been left of its former and probably much larger population. The habitat of the hairy European black poplar in Herzegovina is the consequence of ecogeographical conditions and of rather historical factors and their effect on the development of flora in that area. This belief leads to the idea that the current distribution of the hairy black poplar is a result of a recent spread of the black poplar from the North (from the River Drava and the River Sava, Groups CRO-01 and CRO-02) toward the South. Individual trees of the hairy black poplar in northern and central Bosnia grow at higher altitudes, in the colder and more humid mountainous region . The isolated Submediterranean populations (Group B&H-02), growing along the River Neretva, can be found at 15 m a.s.l. The average temperature during vegetation period ranges between 15.9°C and 16.9°C for populations from inner Bosnia and Herzegovina (Group B&H-01) and is 20.7°C for populations from the lower course of the River Neretva (Group B&H-02). Ecological and climate conditions along the River Neretva are completely different with characteristically long, warm and dry summers (with the average precipitation of 365 mm in vegetation period). Moreover, these conditions are characterised by a distinctive fluctuation in the water level, as well as the saturation of the pedological layer with lime-stones. We can assume that the hairy black poplar found along the River Neretva is a remnant of an earlier continuous land area, spreading from Greece to Herzegovina via Macedonia and Albania. The hairy type of black poplar (P. nigra subsp. caudina) growing along the River Neretva differs considerably from the black poplar in the riparian populations along the European rivers. This rare hairy type of P. nigra is considered to be a xeromorphic form in the Submediterranean region .
In a preliminary attempt to assess the taxo-nomic status of P. nigra subsp. caudina ten trees from each of seven populations were studied for nuclear and chloroplast markers that have shown differences between species, subspecies, and single clones in Populus, , . One of the populations was a pure hairy type, three were hybrids and the others were typical spp. nigra. The results show no indication that the hairy type is a separate taxonomic unit because all the polymorphisms are the same as the ones of the European black poplars from Bosnia and Herzegovina, whether they are hairy or not. These results do not support the hypothesis of a possible tertiary relict but rather they support the classification of the hairy type as a xeromorphic ecotype of the European black poplar. For the characterised strong trichome development on petiols of P. nigra subsp. caudina, candidate genes were analysed in order to understand the basis of its distinctness as well as the genetic variation , .
There are significant differences between the typical and hairy types of European black poplar in the studied morphological traits and these differences are in accordance with the climate differences in respective habitats of continental riparian forests and Submediterranean type of climate.
- VANDEN BROECK A 2003 EUFORGEN Technical Guidelines for genetic conservation and use of European black poplar (Populus nigra). IPGRI, Rome, Italy, 6 p. URL: http://www.euforgen.org/fileadmin/bioversity/publications/pdfs/EUFORGEN/927_Technical_guidelines_for_genetic_conservation_and_use_for_black_poplar__Populus_nigra_.pdf (2 December 2014)
- KAJBA D, VRATARIĆ P 2000 Conservation of European black poplar (Populus nigra) genetic resources in Croatia. In: Isebrands JG, Richardson J (eds) 21st Session of the IPC-2000: poplar and willow culture: meeting the needs of society and the environment, Vancouver, Canada, 24-28 September 2000. US Department of Agriculture, Forest Service, North Central Research Station, Gen Tech Rep NC-215, St Paul MN, USA, pp 87. URL: http://www.nrs.fs.fed.us/pubs/gtr/gtr_nc215.pdf (2 December 2014)
- BALLIAN D 2004 The status of black and white poplars (Populus nigra, P. alba L.) in Bosnia and Herzegovina. In: Koskela J, de Vries SMG, Kajba D, von Wühlisch G (eds) Populus nigra Network, Report of the seventh (25-27 October 2001, Osijek, Croatia) and eight (22-24 May 2003 Treppeln, Germany) meetings. IPGRI, Rome, Italy, pp 17-20
- JANJIĆ N 1984 Prilog poznavanju svojti crne topole ( nigra L.) u Bosni i Hercegovini. Glasnik Zemaljskog muzeja BiH u Sarajevu 23: 255-278
- DZEKOV S 1960 Domašnata crna topola ( nigra) vo NR Makedonija. Godišen Zbornik na Zemjodelsko-šumarskiot fakultet, 44 p
- PAPAIOANNOU J 1963 The Geographical Distribution of nigra var. pubescens Parl. and P. thevestina Dode and Their Hybrids. Food and Agriculture Organization of the United Nations, Rome, Italy, pp 1-15
- BUGALA W 1967 Systematyka euroazjatyckich topoli z grupy nigra L. Rocznik Arboretum Kornickie 12: 47-220
- PETERS RL 1990 Effects of global warming on forests. For Ecol Manag 35 (1-2): 13-33. DOI: http://dx.doi.org/10.1016/0378-1127(90)90229-5
- ERWIN TL 1991 An evolutionary basis for conservation strategies. Science 253 (5021): 750-752. DOI: http://dx.doi.org/10.1126/science.253.5021.750
- ERIKSSON G, NAMKOONG G, ROBERDS JH 1993 Dynamic gene conservation for uncertain futures. For Ecol Manag 62 (1-4): 15-37. DOI: http://dx.doi.org/10.1016/0378-1127(93)90039-p
- AITKEN SN, YEAMAN S, HOLLIDAY JA, WANG T, CURTIS-MCLANE S 2008 Adaptation, migration or extirpation: climate change outcomes for tree populations. Evol Appl 1: 95-111. DOI: http://dx.doi.org/10.1111/j.1752-4571.2007.00013.x
- CAGELLI L, LEFÉVRE F 1995 The conservation of Populus nigra and gene flow with cultivated poplars in Europe. Forest Genetics 2: 135-144
- LEFÈVRE F, KAJBA D, HEINZE B, ROTACH P, DE VRIES SMG, TUROK J 2001 Black poplar: a model for gene resource conservation in forest ecosystems. Forest Chron 77 (2): 239-244. DOI: http://dx.doi.org/10.5558/tfc77239-2
- LEFÈVRE F, BARSOUM N, HEINZE B, KAJBA D, ROTACH P, DE VRIES SMG, TUROK 2001 EUFORGEN Technical Bulletin: In situ Conservation of Populus nigra. IPGRI, Rome, Italy, 58 p
- KAJBA D, BALLIAN D, IDŽOJTIĆ M, BOGDAN S 2004 The differences among hairy and typical European black poplars and the possible role of the hairy type in relation to climatic changes. For Ecol Manag 197 (1-3): 279-284. DOI: http://dx.doi.org/10.1016/j.foreco.2004.05.024
- BISSOFFI S, CAGELLI L 1992 Leaf shape as a tool for the discrimination among poplar clones. 19th Session FAO/IPC, Zaragoza, Spain, FO:CIP:NR 92/2
- ZSUFFA L 1974 The genetics of Populus nigra Annales Forestales 6 (2): 29-53
- ECKENWALDER JE 1996 Systematics and evolution of Populus. In: Stettler RF, Bradshaw Jr HD, Heilman PE, Hinckley TM (eds) Biology of Populus and its Implications for Management and Conservation. National Research Council of Canada, NCR Research Press, Ottawa, Canada, pp 7-32
- KRSTINIĆ A, TRINAJSTIĆ I, KAJBA D, SAMARDŽIĆ J 1997 Morphological variability of the leaves of black poplar (Populus nigra) in natural stands along the Sava river (Croatia). In: Turok J, Lèfevre F, de Vries S, Alba N, Heinze B, van Slycken J (eds) Populus nigra Network, Report of the fourth meeting, Geraardsbergen, Belgium, 3-5 October 1997. IPGRI, Rome, Italy, pp 71-77
- KAJBA D, IDŽOJTIĆ M, BOGDAN S 1999 Discriminant analysis of leaf morphological characters of the European Black Poplar (Populus nigra) in natural population in Croatia. In: Turok J, Lefévre F, de Vries S, Alba N, Heinze B, Voloyanchuk R, Lipman E (eds) Populus nigra Network, Report of the fifth meeting, Kyiv, Ukraine, 5-9 May 1999. IPGRI, Rome, Italy, pp 73-76
- MAKSIMOVIĆ Z, ŠIJAČIĆ-NIKOLIĆ M 2013 Morphometric characteristics of black poplar (Populus nigra) leaves in the area of Great War Island. Bulletin of the Faculty of Forestry 108: 93-108. DOI: http://dx.doi.org/10.2298/gsf1308093m
- SOKAL RR, ROHLF FJ 1989 Biometry. Freeman and Co., San Francisco, CA, USA, 859 p
- MCGARIGAL K, CUSHMAN S, STATFFORD S 2000 Multivariate statistics for wildlife and ecology research. Springer-Verlag, New York, NY, USA, 283 p. DOI: http://dx.doi.org/10.1007/978-1-4612-1288-1
- KAISER HF 1960 The application of electronic computers to factor analysis. Educ Psychol Meast 20: 141-151. DOI: http://dx.doi.org/10.1177/001316446002000116
- STATSOFT INC 2001 STATISTICA (data analysis software system), version 8.0. URL: http://www.statsoft.com (4 December 2014)
- KAJBA D, ROMANIĆ B 2002 Morphological leaf variability of the European black poplar (Populus nigra) in natural populations in the Drava River Basin in Croatia. In: van Damm B, Bordács S (eds) Genetic diversity in river populations of European Black Poplar – implications for riparian eco-system management, Szeksárd, Hungary, 16-20 May 2001. Csiszár Nyomda, Budapest, Hungary, pp 221-227
- ŠILER B, SKORIĆ M, MIŠIĆ D, KOVAČEVIĆ B, JELIĆ M, PATENKOVIĆ A, KURBALIJA NOVIČIĆ Z 2014 Variability of European Black Poplar (Populus nigra) in the Danube Basin. Vojvodinašume, Novi Sad, 128 p.
- RAJORA OP, ZSUFFA L, DANCIK BP 1991 Allozyme and Leaf Morphological Variation of Eastern Cottonwood at the Northern Limits of Its Range in Ontario. Forest Sci 37 (2): 688-702
- BOŽIĆ G, HEINZE B, KAJBA D, KRYSTUFEK V, VANDEN BROECK A 2010 Genetic diversity in Populus nigra populations along rivers in Eastern Austria, Slovenia and Croatia as revealed by microsatellite markers. In: Spano P, Cherubini L (eds) Fifth International Poplar Symposium, Poplars and Willows: from research models to multipurpose trees for a bio-based society, Orvieto, Italy, 20-25 September 2010, pp 22
- POSPIŠKOVÁ M, ŠÁLKOVÁ I 2006 Population structure and parentage analysis of black poplar along the Morava River. Can J Forest Res 36 (5): 1067-1076. DOI: http://dx.doi.org/10.1139/x06-003
- SMULDERS MJM, COTTRELL JE, LEFEVRE F, VAN DER SCHOOT J, ARENS P, VOSMAN B, TABBENER HE, GRASSI F, et al. 2008 Structure of the genetic diversity in black poplar (Populus nigra) populations across European river systems: Consequences for conservation and restoration. For Ecol Manag 255 (5-6): 1388-1399. DOI: http://dx.doi.org/10.1016/j.foreco.2007.10.063
- RATHMACHER G, NIGGEMANN M, KÖHNEN M, ZIEGENHAGEN B, BIALOZYT R 2010 Short-distance gene flow in Populus nigra accounts for small-scale spatial genetic structures: implications for in situ conservation measures. Conserv Genet 11 (4): 1327-1338. DOI: http://dx.doi.org/10.1007/s10592-009-9961-6
- BRUS R, GALIEN U, BOŽIČ G, JARNI K 2010 Morphological study of the leaves of two European black poplar (Populus nigra) populations in Slovenia. Period Biol 112 (3): 317-325
- VANDEN BROECK A, COX K, VAN SLYCKEN J, HALFMAERTEN D 2000 Genetic pollution and mating systems in an artificial stand of black poplar (Populus nigra). In: Isebrands JG, Richardson J (eds) 21st Session of the IPC-2000: poplar and willow culture: meeting the needs of society and the environment, Vancouver, Canada, 24-28 September 2000. US Department of Agriculture, Forest Service, North Central Research Station, Gen Tech Rep NC-215, St Paul, MN, USA, pp 186
- HEINZE B 1998 PCR-based chloroplast DNA assays for the identification of native Populus nigra and introduced poplar hybrids in Europe. Forest Genetics 5 (1): 31-38
- LEXER C, FAY MF, JOSEPH JA, NICA MS, HEINZE B 2005 Barrier to gene flow between two ecologically divergent Populus species, alba (white poplar) and P. tremula (European aspen): the role of ecology and life history in gene introgression. Mol Ecol 14 (4): 1045-1057. DOI: http://dx.doi.org/10.1111/j.1365-294x.2005.02469.x
- KAJBA D, BALLIAN D, HEINZE B, Idžojtić M, Bogdan S 2004 Populus nigra ssp. caudina and ist importance for forest tree improvement and conservation of poplar genetic resources. In: El-Lakany MH (ed) 22nd Session of the International Poplar Commission, The contribution of Poplars and Willows to Suistanable Forestry and Rural Development, Santiago, Chile, 29 November-2 December 2004. FAO, Rome, Italy, pp 29
- KAJBA D, BALLIAN D, HEINZE B, Idžojtić M, Bogdan S, Katičić-Trupčević I 2005 Importance of Conservation of the Hairy Type of Populus nigra in the Changing Climatic and Environmental Conditions. In: Čolić Franekić J, Ugarković Đ (eds) Second Congress of Croatian Geneticists, Supetar, Croatia, 24-27 September 2005. Croatian Genetic Society, Zagreb, Croatia, pp 90
- HASE I, HEINZE B, ARTHOFER W, KAJBA D, BALLIAN D 2007 Candidate genes for hairness in poplar. Plant and Animal Genomes XV Conference, San Diego, USA, 13-17 January 2007
© 2015 by the Croatian Forest Research Institute. This is an Open Access paper distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0).