Comparative demography of the exotic Harmonia axyridis with other aphidophagous coccinellids reared on artificial diet

The objective of this work was to compare biological aspects and life table parameters of the coccinellids Harmonia axyridis, Cycloneda sanguineaand Hippodamia convergens. Insects were fed eggs of Anagasta kuehniella, and reared at 24.5±1ºC, 70±10% relative humidity, with a 12 hour photophase. Hippodamia convergenstook about 1.6 day to complete development, longer than H. axyridis, and 2.4 day longer than C. sanguinea.At immature stages, H. axyridisexhibited the highest survival percentage (49.2%), in comparison to the other coccinellids. For mean adult longevity, H. convergenswas deficient, in comparison with the other species. Mean period of pre oviposition was the longest in C. sanguinea; the longest oviposition time occurred for H. axyridis; and the post oviposition period was similar between the coccinellids. Considering the reproductive parameters, H. axyridisshowed the best performance in all aspects. For life table, the values of H. convergenswere higher than, although close, to those of H. axyridis. Nevertheless, the high net reproductive rate of H. axyridis showed this species potential to increase population size. The biological characteristics of the exotic H. axyridis favors its invasion and establishment in Brazil, corroborating results noticed in other countries.

Demogenetics in metapopulations

SUMMARY : Human-induced habitat fragmentation constitutes a major threat to biodiversity. Small and isolated populations suffer from increased stochasticity and from limited rescue effects. These two factors may be sufficient to cause local extinctions but fragmentation induces some genetic consequences that can also contribute significantly to extinction risks. Increased genetic drift reduces the effectiveness of selection against deleterious mutations, leading to their progressive accumulation. Drift also decreases both the standing genetic variation and the rate of fixation of beneficial mutations, limiting the evolutionary potential of isolated populations. Demography and genetics further interact and feed back on each other, progressively driving fragmented populations into "extinction vortices". The aim of the thesis was to better understand the processes occurring in fragmented populations. For this, I combined simulation studies and empirical data from three species that live in structured habitats. Chapter 1 and 2 investigate the demography of two shrew species in fragmented habitats. I showed that connectivity and habitat quality strongly affect the demography of the greater white-tooted shrew, although demographic stochasticity was extremely high. I also demonstrated that habitat fragmentation is one of the leading factors allowing the local coexistence of two competing shrew species. Chapter 3 and 4 focus on measuring connectivity in fragmented populations based on genetic data. In particular, I showed that genetic data can be used to detect the landscape elements impeding dispersal. In Chapter 5 that deals with the accumulation of deleterious mutations in fragmented populations, I demonstrated that mutation accumulation, as well a time to extinction, can be predicted from simple demographic and genetic measures. In the last two chapters, I monitored individual reproductive success in an isolated tree frogs population. These data allowed quantifying the effective population size, a measure closely linked to population evolutionary potential. To conclude, this thesis brings some new insights into the processes occurring in fragmented populations, and I hope it will contribute to the improvement of the management and conservation of fragmented populations.

Comparative biology of two congeneric stinkbugs, Chinavia impicticornis and C. ubica (Hemiptera: Pentatomidae)

The objective of this work was to compare the biology of Chinavia impicticornis and C. ubica (Hemiptera: Pentatomidae), two species of stinkbugs that occur as secondary pests in soybean. Life table procedures were used for evaluating nymphs, and fecundity tables for evaluating adults, in order to establish the demographics of the two species. The two species have similar demographic parameters, and the development of immature stages, from egg to adult, had similar duration periods of approximately 30 days. In both species, eggs and second-instar nymphs were the stages with higher mortality. Total egg production did not differ between females of both species. Fecundity and survival curves for adults showed similar trends in both species. However, C. ubica had greater potential to increase its populations, since its fecundity parameters were significantly higher than those of C. impicticornis. Moreover, the generational time and the time required to double the population size were shorter in C. ubica. Prolonged longevity, long oviposition period, high fecundity, and the ability to rapidly increase their populations indicate that both species can become potential pests in favorable environments.

RiskDiff: A web tool for the analysis of the difference due to risk and demographic factors for incidence or mortality data.

Background Analysing the observed differences for incidence or mortality of a particular disease between two different situations (such as time points, geographical areas, gender or other social characteristics) can be useful both for scientific or administrative purposes. From an epidemiological and public health point of view, it is of great interest to assess the effect of demographic factors in these observed differences in order to elucidate the effect of the risk of developing a disease or dying from it. The method proposed by Bashir and Estève, which splits the observed variation into three components: risk, population structure and population size is a common choice at practice. Results A web-based application, called RiskDiff has been implemented (available at http://rht.iconcologia.net/riskdiff.htm webcite), to perform this kind of statistical analyses, providing text and graphical summaries. Code from the implemented functions in R is also provided. An application to cancer mortality data from Catalonia is used for illustration. Conclusions Combining epidemiological with demographical factors is crucial for analysing incidence or mortality from a disease, especially if the population pyramids show substantial differences. The tool implemented may serve to promote and divulgate the use of this method to give advice for epidemiologic interpretation and decision making in public health.

Sex-chromosome differentiation and 'sex races' in the common frog (Rana temporaria).

Sex-chromosome differentiation was recently shown to vary among common frog populations in Fennoscandia, suggesting a trend of increased differentiation with latitude. By rearing families from two contrasted populations (respectively, from northern and southern Sweden), we show this disparity to stem from differences in sex-determination mechanisms rather than in XY-recombination patterns. Offspring from the northern population display equal sex ratios at metamorphosis, with phenotypic sexes that correlate strongly with paternal LG2 haplotypes (the sex chromosome); accordingly, Y haplotypes are markedly differentiated, with male-specific alleles and depressed diversity testifying to their smaller effective population size. In the southern population, by contrast, a majority of juveniles present ovaries at metamorphosis; only later in development do sex ratios return to equilibrium. Even at these later stages, phenotypic sexes correlate only mildly with paternal LG2 haplotypes; accordingly, there are no recognizable Y haplotypes. These distinct patterns of gonadal development fit the concept of 'sex races' proposed in the 1930s, with our two populations assigned to the 'differentiated' and 'semi-differentiated' races, respectively. Our results support the suggestion that 'sex races' differ in the genetic versus epigenetic components of sex determination. Analysing populations from the 'undifferentiated race' with high-density genetic maps should help to further test this hypothesis.

Normaalijakaumaan perustuva mutaatio-operaatio differentiaalievoluutioalgoritmiin

Evoluutioalgoritmit ovat viime vuosina osoittautuneet tehokkaiksi menetelmiksi globaalien optimointitehtävien ratkaisuun. Niiden vahvuutena on etenkin yleiskäyttöisyys ja kyky löytää globaali ratkaisu juuttumatta optimoitavan tavoitefunktion paikallisiin optimikohtiin. Tässä työssä on tavoitteena kehittää uusi, normaalijakaumaan perustuva mutaatio-operaatio differentiaalievoluutioalgoritmiin, joka on eräs uusimmista evoluutiopohjaisista optimointialgoritmeista. Menetelmän oletetaan vähentävän entisestään sekä populaation ennenaikaisen suppenemisen, että algoritmin tilojen juuttumisen riskiä ja se on teoreettisesti osoitettavissa suppenevaksi. Tämä ei päde alkuperäisen differentiaalievoluution tapauksessa, koska on voitu osoittaa, että sen tilanmuutokset voivat pienellä todennäköisyydellä juuttua. Työssä uuden menetelmän toimintaa tarkastellaan kokeellisesti käyttäen testiongelmina monirajoiteongelmia. Rajoitefunktioiden käsittelyyn käytetään Jouni Lampisen kehittämää, Pareto-optimaalisuuden periaatteeseen perustuvaa menetelmää. Samalla saadaan kerättyä lisää kokeellista näyttöä myös tämän menetelmän toiminnasta. Kaikki käytetyt testiongelmat kyettiin ratkaisemaan sekä alkuperäisellä differentiaalievoluutiolla, että uutta mutaatio-operaatiota käyttävällä versiolla. Uusi menetelmä osoittautui kuitenkin luotettavammaksi sellaisissa tapauksissa, joissa alkuperäisellä algoritmilla oli vaikeuksia. Lisäksi useimmat ongelmat kyettiin ratkaisemaan luotettavasti pienemmällä populaation koolla kuin alkuperäistä differentiaalievoluutiota käytettäessä. Uuden menetelmän käyttö myös mahdollistaa paremmin sellaisten kontrolliparametrien käytön, joilla hausta saadaan rotaatioinvariantti. Laskennallisesti uusi menetelmä on hieman alkuperäistä differentiaalievoluutiota raskaampi ja se tarvitsee yhden kontrolliparametrin enemmän. Uusille kontrolliparametreille määritettiin kuitenkin mahdollisimman yleiskäyttöiset arvot, joita käyttämällä on mahdollista ratkaista suuri joukko erilaisia ongelmia.

Populaation monimuotoisuuden mittaaminen liukulukukoodatuissa evoluutioalgoritmeissa

Diplomityössä esitetään menetelmä populaation monimuotoisuuden mittaamiseen liukulukukoodatuissa evoluutioalgoritmeissa, ja tarkastellaan kokeellisesti sen toimintaa. Evoluutioalgoritmit ovat populaatiopohjaisia menetelmiä, joilla pyritään ratkaisemaan optimointiongelmia. Evoluutioalgoritmeissa populaation monimuotoisuuden hallinta on välttämätöntä, jotta suoritettu haku olisi riittävän luotettavaa ja toisaalta riittävän nopeaa. Monimuotoisuuden mittaaminen on erityisen tarpeellista tutkittaessa evoluutioalgoritmien dynaamista käyttäytymistä. Työssä tarkastellaan haku- ja tavoitefunktioavaruuden monimuotoisuuden mittaamista. Toistaiseksi ei ole ollut olemassa täysin tyydyttäviä monimuotoisuuden mittareita, ja työn tavoitteena on kehittää yleiskäyttöinen menetelmä liukulukukoodattujen evoluutioalgoritmien suhteellisen ja absoluuttisen monimuotoisuuden mittaamiseen hakuavaruudessa. Kehitettyjen mittareiden toimintaa ja käyttökelpoisuutta tarkastellaan kokeellisesti ratkaisemalla optimointiongelmia differentiaalievoluutioalgoritmilla. Toteutettujen mittareiden toiminta perustuu keskihajontojen laskemiseen populaatiosta. Keskihajonnoille suoritetaan skaalaus, joko alkupopulaation tai nykyisen populaation suhteen, riippuen lasketaanko absoluuttista vai suhteellista monimuotoisuutta. Kokeellisessa tarkastelussa havaittiin kehitetyt mittarit toimiviksi ja käyttökelpoisiksi. Tavoitefunktion venyttäminen koordinaattiakseleiden suunnassa ei vaikuta mittarin toimintaan. Myöskään tavoitefunktion kiertäminen koordinaatistossa ei vaikuta mittareiden tuloksiin. Esitetyn menetelmän aikakompleksisuus riippuu lineaarisesti populaation koosta, ja mittarin toiminta on siten nopeaa suuriakin populaatioita käytettäessä. Suhteellinen monimuotoisuus antaa vertailukelpoisia tuloksia riippumatta parametrien lukumäärästä tai populaation koosta.

Online and Offline Participation at the Local Level : A Quantitative Analysis of the Catalan Municipalities

This paper presents a first analysis on local electronic participatory experiences in Catalonia. The analysis is based on a database constructed and collected by the authors. The paper carries out an explanatory analysis of local initiatives in eparticipationand off line participation taking into account political variables (usually not considered in this kind of analysis) but also classical socio-economic variables that characterise municipalities. Hence, we add a quantitative analysis to the numerous case studies on local e-participation experiences. We have chosen Catalonia because is one of the European regions with more initiatives and one that has enjoyed considerable local governmental support to citizen participation initiatives since the 80s. The paper offers a characterisation of these experiences and a first explanatory analysis, considering: i) the institutional context in which these experiences are embedded, ii) the characteristics of the citizen participation processes and mechanisms on-line, and iii) a set of explanatory variables composed by the population size, thepolitical adscription of the mayor, the electoral abstention rate, age, income and level ofeducation in the municipality. The model that we present is explanatory for the municipalities with more than 20,000 inhabitants but it is not for the fewer than 20,000inhabitants. Actually, the number of participatory activities developed by these last municipalities is very low. Among all the variables, population size becomes the mostinfluential variable. Political variables such as political party of the mayor and the localabstention rate have a certain influence but that have to be controlled by population size.

Online and Offline Participation at the Local level: A Quantitative Analysis of the Catalan Municipalities

This article presents an analysis on local participatory experiences in Catalonia,both online and in-person. The analysis is based on a database set up by theauthors. The article carries out an explanatory analysis of local participatoryinitiatives (on- and offline) taking into account political variables (not usually con-sidered in this kind of analysis) and also classical socio-economic variables thatcharacterize municipalities. Hence, we add a quantitative analysis to the numerouscase studies on local e-participation experiences. We have chosen Catalonia becauseit is one of the European regions with more initiatives and a considerable localgovernment support for citizen participation initiatives since the 1980s. Thearticle offers a characterization of these experiences and an explanatory analysis,considering: (i) the institutional context in which these experiences are embedded,(ii) the citizen participation processes and mechanisms online and (iii) a set ofexplanatory variables composed of the population size and the province to whichthe municipality belongs, the political tendency of the mayor, the electoral absten-tion rate, age, income, level of education, broadband connection and users of theInternet in the municipality. The model that we present is explanatory for munici-palities with more than 20,000 inhabitants but it is not for fewer than 20,000inhabitants. Actually, the majority of these latter municipalities have not developedany participatory activities. Among all the variables, population size is the mostinfluential variable and affects the influence of other variables, such as the politicalparty of the mayor, the local abstention rate and the province.

Deu anys de seguiment demogràfic i genètic d’Stachys maritima a Catalunya (2001-2010). Implicacions per a un pla de recuperació

Ten years of demographic and genetic monitoring of Stachys maritima in Catalonia (2001-2010). Implications for a recovery plan.- Stachys maritima is a species typical of the coastal dunes, with a wide distribution within the Mediterranean Basin. In spite of this, the species shows a clear regression. In Catalonia, it has been observed an important reduction of its populations since early 20th century, where it has disappeared from several localities in which it was relatively common (Tarragonès, Barcelonès). Herein we present the results of the demographic monitoring of the species during the last 10 years (2001-2010) in the known localities in Catalonia. Besides corroborating the disappearance (northern Sant Martí d’Empúries), the re-discovering (Llobregat Delta beach) and the detection of new populations (inner dunes of the Montgrí), a large year-to-year fluctuation of the monitored populations is stated; the possible reasons are discussed. In addition, the present work also includes the results of the allozyme diversity analysis of the new detected populations as well as the rediscoveries of the period 2004-2008, which were not included in a former study of genetic diversity carried out in 2002-2003. It is necessary to emphasize that the contribution of the new populations to the genetic diversity of Stachys maritima is very small, which can be attributed to their limited population size and /or to founder effects. Despite that the species is included in the Annex 2 (“En Perill d’Extinció”) within the Catàleg de Flora Amenaçada de Catalunya (Catalogue of Endangered Flora of Catalonia), and some “soft” conservation measures have been applied at local level (signposting of the beach accesses, environmental education, etc.) coupled with other more significant measures (e.g. translocation of individuals discovered in an artificial sandbank), it would be necessary the coordinated action and the scientific support of any initiative of conservation that could be carried out. The general frame to initiate actions of conservation should be the recovery plan of Stachys maritima, whose draft and application is mandatory in accordance to the Catàleg

Conservation Genetics of Exploited Finnish Salmonid Fishes

Genetic diversity is one of the levels of biodiversity that the World Conservation Union (IUCN) has recognized as being important to preserve. This is because genetic diversity is fundamental to the future evolution and to the adaptive flexibility of a species to respond to the inherently dynamic nature of the natural world. Therefore, the key to maintaining biodiversity and healthy ecosystems is to identify, monitor and maintain locally-adapted populations, along with their unique gene pools, upon which future adaptation depends. Thus, conservation genetics deals with the genetic factors that affect extinction risk and the genetic management regimes required to minimize the risk. The conservation of exploited species, such as salmonid fishes, is particularly challenging due to the conflicts between different interest groups. In this thesis, I conduct a series of conservation genetic studies on primarily Finnish populations of two salmonid fish species (European grayling, Thymallus thymallus, and lake-run brown trout, Salmo trutta) which are popular recreational game fishes in Finland. The general aim of these studies was to apply and develop population genetic approaches to assist conservation and sustainable harvest of these populations. The approaches applied included: i) the characterization of population genetic structure at national and local scales; ii) the identification of management units and the prioritization of populations for conservation based on evolutionary forces shaping indigenous gene pools; iii) the detection of population declines and the testing of the assumptions underlying these tests; and iv) the evaluation of the contribution of natural populations to a mixed stock fishery. Based on microsatellite analyses, clear genetic structuring of exploited Finnish grayling and brown trout populations was detected at both national and local scales. Finnish grayling were clustered into three genetically distinct groups, corresponding to northern, Baltic and south-eastern geographic areas of Finland. The genetic differentiation among and within population groups of grayling ranged from moderate to high levels. Such strong genetic structuring combined with low genetic diversity strongly indicates that genetic drift plays a major role in the evolution of grayling populations. Further analyses of European grayling covering the majority of the species’ distribution range indicated a strong global footprint of population decline. Using a coalescent approach the beginning of population reduction was dated back to 1 000-10 000 years ago (ca. 200-2 000 generations). Forward simulations demonstrated that the bottleneck footprints measured using the M ratio can persist within small populations much longer than previously anticipated in the face of low levels of gene flow. In contrast to the M ratio, two alternative methods for genetic bottleneck detection identified recent bottlenecks in six grayling populations that warrant future monitoring. Consistent with the predominant role of random genetic drift, the effective population size (N_e) estimates of all grayling populations were very low with the majority of N_e estimates below 50. Taken together, highly structured local populations, limited gene flow and the small N_e of grayling populations indicates that grayling populations are vulnerable to overexploitation and, hence, monitoring and careful management using the precautionary principles is required not only in Finland but throughout Europe. Population genetic analyses of lake-run brown trout populations in the Inari basin (northernmost Finland) revealed hierarchical population structure where individual populations were clustered into three population groups largely corresponding to different geographic regions of the basin. Similar to my earlier work with European grayling, the genetic differentiation among and within population groups of lake-run brown trout was relatively high. Such strong differentiation indicated that the power to determine the relative contribution of populations in mixed fisheries should be relatively high. Consistent with these expectations, high accuracy and precision in mixed stock analysis (MSA) simulations were observed. Application of MSA to indigenous fish caught in the Inari basin identified altogether twelve populations that contributed significantly to mixed stock fisheries with the Ivalojoki river system being the major contributor (70%) to the total catch. When the contribution of wild trout populations to the fisheries was evaluated regionally, geographically nearby populations were the main contributors to the local catches. MSA also revealed a clear separation between the lower and upper reaches of Ivalojoki river system – in contrast to lower reaches of the Ivalojoki river that contributed considerably to the catch, populations from the upper reaches of the Ivalojoki river system (>140 km from the river mouth) did not contribute significantly to the fishery. This could be related to the available habitat size but also associated with a resident type life history and increased cost of migration. The studies in my thesis highlight the importance of dense sampling and wide population coverage at the scale being studied and also demonstrate the importance of critical evaluation of the underlying assumptions of the population genetic models and methods used. These results have important implications for conservation and sustainable fisheries management of Finnish populations of European grayling and brown trout in the Inari basin.

Genetic structure and diversity of Copaifera langsdorffii Desf. in Cerrado fragments of the São Paulo State, Brazil

The loss of large areas of Cerrado (Brazilian savanna) in Brazil can lead to reduced biodiversity and to the extinction of species. Therefore, the present study aimed to investigate the genetic fragility of populations of Copaifera langsdorffii Desf exposed to different anthropic conditions in fragments of Cerrado in the state of São Paulo. The study was carried out in two Experimental Stations operated by the Forest Institute (Assis and Itirapina), in one fully protected conservation unit (Pedregulho) and in one private property (Brotas). Analyses were conducted using leaf samples from 353 adult specimens and eight pairs of microsatellite loci. The number of alleles per locus ranged from 13 to 15 in all populations, but the mean number of effective alleles was approximately half this value (7.2 to 9-1). Observed heterozygosity was significant and lower than the expected in all populations. Consequently, all populations deviated from Hardy-Weinberg expected frequencies. Fixation indexes were significant for all populations, with the Pedregulho population having the lowest value (0.189) and Itirapina having the highest (0.283). The analysis of spatial genetic structure detected family structures at distance classes of 20 to 65 m in the populations studied. No clones were detected in the populations. Estimates of effective population size were low, but the area occupied by each population studied was large enough for conservation, medium and long term. Recent reductions or bottlenecks were detected in all four populations. Mean Gst’ (genetic divergence) indicated that most of the variation was within populations. Cluster structure analysis based on the genotypes detected K= 4 clusters with distinct allele frequencies patterns. The genetic differentiation observed among populations is consistent with the hypothesis of genetic and geographic isolation. Therefore, it is essential to adopt conservation strategies that raise the gene flow between fragments.

Conservação genética de populações de Cryptocarya moschata Nees (Lauraceae) na Mata Atlântica do estado de São Paulo

Através da análise de 18 locos isoenzimáticos polimórficos, foram estimadas as freqüências alélicas referentes a 214 indivíduos adultos de quatro populações naturais de Cryptocarya moschata de duas regiões do estado de São Paulo. Com base nas heterozigosidades observadas e esperadas, foram obtidas estimativas das estatísticas F de Wright. Para fins de comparação, utilizou-se também o método da análise da variância para estimação dos parâmetros correspondentes F = FIT, qP= FST e f = FIS. Os dois métodos forneceram resultados bastante concordantes: IT = 0,142; ST = 0,140; IS = 0,002 e = 0,116; P = 0,123 e = ­0,008, indicando que os indivíduos dentro das populações devem ser panmíticos e que a diversidade entre populações é bastante alta, sendo similar à que se espera para famílias com estruturação de meios-irmãos. Calculando P com as populações tomadas duas a duas, notou-se tendência de P crescer com a distância geográfica o que sugere tendência de isolamento pela distância. O fluxo gênico foi estimado em 0,9 indivíduos por geração, o que corrobora a pronunciada diferenciação populacional encontrada. Devido ao valor negligível encontrado de IS, o tamanho efetivo de variância de cada população é equivalente ao número de indivíduos amostrados. As estratégias de manejo e conservação necessárias para a preservação da alta variabilidade genética intrapopulacional de C. moschata implicam na manutenção de populações com número grande de indivíduos. Além disso, para a preservação da espécie como um todo, a manutenção de muitas populações provavelmente é necessária.

The post-fire effects on the outcrossing rate of a Brazilian savannah shrub, Helicteres sacarolha A.St.-Hil

The effect of post-fire, plant density and number of flowers per plant on outcrossing rates was examined in a Brazilian savannah shrub, Helicteres sacarolha. Data on number of flowering plants and flowers per plant were collected during the flowering season in January and February of 1994. In October of 1994, a fire swept through the studied area after 30 years of fire absence. The burnt plants of H. sacarolha resprouted, producing flowers and fruits in 1995. Seeds from several plants in both years were collected. Allozyme markers were used to estimate the multilocus outcrossing rates for 1994 and 1995 progenies. After the fire, outcrossing rates increased. In 1995, plants flowered vigorously, increasing flower density and probably pollinator activity. Helicteres sacarolha seems to be fire resistant, like other plant species of the Brazilian savannah, but several plants tagged in 1994 were not found after the fire, and may have died. Thus, although genetic diversity and outcrossing rates increased following the fire in 1994, repeated events of fire may decrease drastically the population size of H. sacarolha, leading to a diminished genetic diversity and outcrossing rates.

Mating system in Myracrodruon urundeuva (Anarcadiaceae): implications for conservation genetics

The aims of this study were to investigate the mating system of a fragmented population of the dioecious tropical tree Myracrodruon urundeuva Allemão, using five microsatellite loci and the mixed mating and correlated mating models. The study was conducted in the Estação Ecológica de Paulo de Farias (436 ha), where the population occupies about 142 ha. The mating system was estimated using 514 open-pollinated offspring, collected from 30 seed-trees. Estimates of the multilocus outcrossing rate confirm that the species is dioecious (t m = 1.0). Low levels of mating among relatives were detected in the population (1 - t s = 0.020). The estimate of paternity correlation (r p(m)) indicated that offsprings were composed of mixtures of half-sibs and full-sibs, with the latter occurring at a low frequency (average of 0.148). The estimated coancestry coefficient within families (Θ = 0.147) was larger and the effective population size (Ne(v)) was lower (Ne(v) = 2.98) than expected in progenies from panmictic populations (Θ = 0.125, Ne(v) = 4, respectively). In terms of conservation, the results indicate that to retain an effective population size of 150, is necessary to collect seeds from at least 50 seed-trees.