Modelling tree population dynamics at the alpine and boreal tree-line ecotones in response to climate and land-use change

Summary Ecotones are sensitive to change because they contain high numbers of species living at the margin of their environmental tolerance. This is equally true of tree-lines, which are determined by attitudinal or latitudinal temperature gradients. In the current context of climate change, they are expected to undergo modifications in position, tree biomass and possibly species composition. Attitudinal and latitudinal tree-lines differ mainly in the steepness of the underlying temperature gradient: distances are larger at latitudinal tree-lines, which could have an impact on the ability of tree species to migrate in response to climate change. Aside from temperature, tree-lines are also affected on a more local level by pressure from human activities. These are also changing as a consequence of modifications in our societies and may interact with the effects of climate change. Forest dynamics models are often used for climate change simulations because of their mechanistic processes. The spatially-explicit model TreeMig was used as a base to develop a model specifically tuned for the northern European and Alpine tree-line ecotones. For the latter, a module for land-use change processes was also added. The temperature response parameters for the species in the model were first calibrated by means of tree-ring data from various species and sites at both tree-lines. This improved the growth response function in the model, but also lead to the conclusion that regeneration is probably more important than growth for controlling tree-line position and species' distributions. The second step was to implement the module for abandonment of agricultural land in the Alps, based on an existing spatial statistical model. The sensitivity of its most important variables was tested and the model's performance compared to other modelling approaches. The probability that agricultural land would be abandoned was strongly influenced by the distance from the nearest forest and the slope, bath of which are proxies for cultivation costs. When applied to a case study area, the resulting model, named TreeMig-LAb, gave the most realistic results. These were consistent with observed consequences of land-abandonment such as the expansion of the existing forest and closing up of gaps. This new model was then applied in two case study areas, one in the Swiss Alps and one in Finnish Lapland, under a variety of climate change scenarios. These were based on forecasts of temperature change over the next century by the IPCC and the HadCM3 climate model (ΔT: +1.3, +3.5 and +5.6 °C) and included a post-change stabilisation period of 300 years. The results showed radical disruptions at both tree-lines. With the most conservative climate change scenario, species' distributions simply shifted, but it took several centuries reach a new equilibrium. With the more extreme scenarios, some species disappeared from our study areas (e.g. Pinus cembra in the Alps) or dwindled to very low numbers, as they ran out of land into which they could migrate. The most striking result was the lag in the response of most species, independently from the climate change scenario or tree-line type considered. Finally, a statistical model of the effect of reindeer (Rangifer tarandus) browsing on the growth of Pinus sylvestris was developed, as a first step towards implementing human impacts at the boreal tree-line. The expected effect was an indirect one, as reindeer deplete the ground lichen cover, thought to protect the trees against adverse climate conditions. The model showed a small but significant effect of browsing, but as the link with the underlying climate variables was unclear and the model was not spatial, it was not usable as such. Developing the TreeMig-LAb model allowed to: a) establish a method for deriving species' parameters for the growth equation from tree-rings, b) highlight the importance of regeneration in determining tree-line position and species' distributions and c) improve the integration of social sciences into landscape modelling. Applying the model at the Alpine and northern European tree-lines under different climate change scenarios showed that with most forecasted levels of temperature increase, tree-lines would suffer major disruptions, with shifts in distributions and potential extinction of some tree-line species. However, these responses showed strong lags, so these effects would not become apparent before decades and could take centuries to stabilise. Résumé Les écotones son sensibles au changement en raison du nombre élevé d'espèces qui y vivent à la limite de leur tolérance environnementale. Ceci s'applique également aux limites des arbres définies par les gradients de température altitudinaux et latitudinaux. Dans le contexte actuel de changement climatique, on s'attend à ce qu'elles subissent des modifications de leur position, de la biomasse des arbres et éventuellement des essences qui les composent. Les limites altitudinales et latitudinales diffèrent essentiellement au niveau de la pente des gradients de température qui les sous-tendent les distance sont plus grandes pour les limites latitudinales, ce qui pourrait avoir un impact sur la capacité des espèces à migrer en réponse au changement climatique. En sus de la température, la limite des arbres est aussi influencée à un niveau plus local par les pressions dues aux activités humaines. Celles-ci sont aussi en mutation suite aux changements dans nos sociétés et peuvent interagir avec les effets du changement climatique. Les modèles de dynamique forestière sont souvent utilisés pour simuler les effets du changement climatique, car ils sont basés sur la modélisation de processus. Le modèle spatialement explicite TreeMig a été utilisé comme base pour développer un modèle spécialement adapté pour la limite des arbres en Europe du Nord et dans les Alpes. Pour cette dernière, un module servant à simuler des changements d'utilisation du sol a également été ajouté. Tout d'abord, les paramètres de la courbe de réponse à la température pour les espèces inclues dans le modèle ont été calibrées au moyen de données dendrochronologiques pour diverses espèces et divers sites des deux écotones. Ceci a permis d'améliorer la courbe de croissance du modèle, mais a également permis de conclure que la régénération est probablement plus déterminante que la croissance en ce qui concerne la position de la limite des arbres et la distribution des espèces. La seconde étape consistait à implémenter le module d'abandon du terrain agricole dans les Alpes, basé sur un modèle statistique spatial existant. La sensibilité des variables les plus importantes du modèle a été testée et la performance de ce dernier comparée à d'autres approches de modélisation. La probabilité qu'un terrain soit abandonné était fortement influencée par la distance à la forêt la plus proche et par la pente, qui sont tous deux des substituts pour les coûts liés à la mise en culture. Lors de l'application en situation réelle, le nouveau modèle, baptisé TreeMig-LAb, a donné les résultats les plus réalistes. Ceux-ci étaient comparables aux conséquences déjà observées de l'abandon de terrains agricoles, telles que l'expansion des forêts existantes et la fermeture des clairières. Ce nouveau modèle a ensuite été mis en application dans deux zones d'étude, l'une dans les Alpes suisses et l'autre en Laponie finlandaise, avec divers scénarios de changement climatique. Ces derniers étaient basés sur les prévisions de changement de température pour le siècle prochain établies par l'IPCC et le modèle climatique HadCM3 (ΔT: +1.3, +3.5 et +5.6 °C) et comprenaient une période de stabilisation post-changement climatique de 300 ans. Les résultats ont montré des perturbations majeures dans les deux types de limites de arbres. Avec le scénario de changement climatique le moins extrême, les distributions respectives des espèces ont subi un simple glissement, mais il a fallu plusieurs siècles pour qu'elles atteignent un nouvel équilibre. Avec les autres scénarios, certaines espèces ont disparu de la zone d'étude (p. ex. Pinus cembra dans les Alpes) ou ont vu leur population diminuer parce qu'il n'y avait plus assez de terrains disponibles dans lesquels elles puissent migrer. Le résultat le plus frappant a été le temps de latence dans la réponse de la plupart des espèces, indépendamment du scénario de changement climatique utilisé ou du type de limite des arbres. Finalement, un modèle statistique de l'effet de l'abroutissement par les rennes (Rangifer tarandus) sur la croissance de Pinus sylvestris a été développé, comme première étape en vue de l'implémentation des impacts humains sur la limite boréale des arbres. L'effet attendu était indirect, puisque les rennes réduisent la couverture de lichen sur le sol, dont on attend un effet protecteur contre les rigueurs climatiques. Le modèle a mis en évidence un effet modeste mais significatif, mais étant donné que le lien avec les variables climatiques sous jacentes était peu clair et que le modèle n'était pas appliqué dans l'espace, il n'était pas utilisable tel quel. Le développement du modèle TreeMig-LAb a permis : a) d'établir une méthode pour déduire les paramètres spécifiques de l'équation de croissance ä partir de données dendrochronologiques, b) de mettre en évidence l'importance de la régénération dans la position de la limite des arbres et la distribution des espèces et c) d'améliorer l'intégration des sciences sociales dans les modèles de paysage. L'application du modèle aux limites alpines et nord-européennes des arbres sous différents scénarios de changement climatique a montré qu'avec la plupart des niveaux d'augmentation de température prévus, la limite des arbres subirait des perturbations majeures, avec des glissements d'aires de répartition et l'extinction potentielle de certaines espèces. Cependant, ces réponses ont montré des temps de latence importants, si bien que ces effets ne seraient pas visibles avant des décennies et pourraient mettre plusieurs siècles à se stabiliser.

Population pharmacokinetic approach to evaluate the effect of CYP2D6, CYP3A, ABCB1, POR and NR1I2 genotypes on donepezil clearance.

AIMS: A large interindividual variability in plasma concentrations has been reported in patients treated with donepezil, the most frequently prescribed antidementia drug. We aimed to evaluate clinical and genetic factors influencing donepezil disposition in a patient population recruited from a naturalistic setting. METHODS: A population pharmacokinetic study was performed including data from 129 older patients treated with donepezil. The patients were genotyped for common polymorphisms in the metabolic enzymes CYP2D6 and CYP3A, in the electron transferring protein POR and the nuclear factor NR1I2 involved in CYP activity and expression, and in the drug transporter ABCB1. RESULTS: The average donepezil clearance was 7.3 l h(-1) with a 30% interindividual variability. Gender markedly influenced donepezil clearance (P < 0.01). Functional alleles of CYP2D6 were identified as unique significant genetic covariate for donepezil clearance (P < 0.01), with poor metabolizers and ultrarapid metabolizers demonstrating, respectively, a 32% slower and a 67% faster donepezil elimination compared with extensive metabolizers. CONCLUSION: The pharmacokinetic parameters of donepezil were well described by the developed population model. Functional alleles of CYP2D6 significantly contributed to the variability in donepezil disposition in the patient population and should be further investigated in the context of individual dose optimization to improve clinical outcome and tolerability of the treatment.

Effects of plant population and nitrogen fertilizer on yield and efficiency of maize-bean intercropping

Nitrogen supply and plant population are basic parameters for cereal-legume intercropping. In order to study plant population and nitrogen fertilizer effects on yield and yield efficiency of maize-bean intercropping, a field experiment was established. Three bean plant populations and three nitrogen levels were used. Maize dry matter accumulation decreased with increases in bean plant population. Competitive effect of intercrop beans on maize yields was high at higher plant populations, being decreased by nitrogen fertilizer; application of 50 kg ha-1 N was very efficient in increasing maize cob yield. Intercropping significantly decreased harvest index of beans in all plant population and nitrogen fertilizer situations. The efficiency of intercropping, compared to sole cropping, was evidenced by the values obtained for Land Equivalent Ratio (LER) for biomass, cob and pod yields that increased with increases in bean plant populations and nitrogen fertilizer levels.

Testing the competitive exclusion principle using various niche parameters in a native (Natrix maura) and an introduced (N. tessellata) colubrid

Despite the increase of animal and plant introductions worldwide and the strong augmentation of the reptile trade, few invasive snake populations have been studied. Dice snakes (Natrix tessellata) were introduced to the shores of Lake Geneva (Switzerland) in the early 1920s, and are now well established. This region of introduction was previously inhabited by Viperine snakes (N. maura). Ever since these two species have been under monitoring (which began in 1996) the Viperine snake population has shown drastic decline. We examine here the possibility of trophic competition by analysing diet composition, prey size and trophic niche overlap. Spatial distribution is also assessed in order to address the question of spatial competitive exclusion. We found very similar diets, and thus a high trophic niche overlap, indicating no partitioning of the trophic resource. No arguments in favour of spatial competitive exclusion were found. Our study suggests that trophic competition may occur between the two natricines and that it may give an explanation for the drastic decline of the Viperine snake in this area. Other pathways potentially playing a role in the exclusion of the Viperine snake are discussed.

Is use of nifurtimox for the treatment of Chagas disease compatible with breast feeding? A population pharmacokinetics analysis.

Introduction Women with Chagas disease receiving treatment with nifurtimox are discouraged from breast feeding. Many patients who would receive treatment with nifurtimox live in extreme poverty, have limited access to resources such as clean water and baby formula and may not have safe alternatives to breast milk. Aim We aimed to estimate, using limited available pharmacokinetics data, potential infant exposure to nifurtimox through breast milk. Methods Original nifurtimox plasma concentrations were obtained from published studies. Pharmacokinetic parameters were estimated using non-linear mixed-effect modelling with NONMEM V.VI. A total of 1000 nifurtimox plasma-concentration profiles were simulated and used to calculate the amount of drug that an infant would be exposed to, if breast fed 150 ml/kg/day. Results Breast milk concentrations on the basis of peak plasma levels (1361 ng/ml) and milk-plasma ratio were estimated. We calculated infant nifurtimox exposure of a breastfed infant of a mother treated with this drug to be below 10% of the maternal weight-adjusted dose, even if milk-plasma ratio were overestimated. Simulation led to similar estimates. Discussion Risk for significant infant exposure to nifurtimox through breast milk seems small and below the level of exposure of infants with Chagas disease receiving nifurtimox treatment. This potential degree of exposure may not justify discontinuation of breast feeding.

Basic density and pulp yield relationship with some chemical parameters in eucalyptus trees

The objective of this work was to evaluate the use of basic density and pulp yield correlations with some chemical parameters, in order to differentiate an homogeneous eucalyptus tree population, in terms of its potential for pulp production or some other technological applications. Basic density and kraft pulp yield were determined for 120 Eucalyptus globulus trees, and the values were plotted as frequency distributions. Homogenized samples from the first and fourth density quartiles and first and fourth yield quartiles were submitted to total phenols, total sugars and methoxyl group analysis. Syringyl/guaiacyl (S/G) and syringaldehyde/vanillin (S/V) ratios were determined on the kraft lignins from wood of the same quartiles. The results show the similarity between samples from high density and low yield quartiles, both with lower S/G (3.88-4.12) and S/V (3.99-4.09) ratios and higher total phenols (13.3-14.3 g gallic acid kg-1 ). Woods from the high yield quartile are statistically distinguished from all the others because of their higher S/G (5.15) and S/V (4.98) ratios and lower total phenols (8.7 g gallic acid kg-1 ). Methoxyl group and total sugars parameters are more adequate to distinguish wood samples with lower density.

Reference values and factors associated with renal resistive index in a family-based population study.

Increased renal resistive index (RRI) has been recently associated with target organ damage and cardiovascular or renal outcomes in patients with hypertension and diabetes mellitus. However, reference values in the general population and information on familial aggregation are largely lacking. We determined the distribution of RRI, associated factors, and heritability in a population-based study. Families of European ancestry were randomly selected in 3 Swiss cities. Anthropometric parameters and cardiovascular risk factors were assessed. A renal Doppler ultrasound was performed, and RRI was measured in 3 segmental arteries of both kidneys. We used multilevel linear regression analysis to explore the factors associated with RRI, adjusting for center and family relationships. Sex-specific reference values for RRI were generated according to age. Heritability was estimated by variance components using the ASSOC program (SAGE software). Four hundred women (mean age±SD, 44.9±16.7 years) and 326 men (42.1±16.8 years) with normal renal ultrasound had mean RRI of 0.64±0.05 and 0.62±0.05, respectively (P<0.001). In multivariable analyses, RRI was positively associated with female sex, age, systolic blood pressure, and body mass index. We observed an inverse correlation with diastolic blood pressure and heart rate. Age had a nonlinear association with RRI. We found no independent association of RRI with diabetes mellitus, hypertension treatment, smoking, cholesterol levels, or estimated glomerular filtration rate. The adjusted heritability estimate was 42±8% (P<0.001). In a population-based sample with normal renal ultrasound, RRI normal values depend on sex, age, blood pressure, heart rate, and body mass index. The significant heritability of RRI suggests that genes influence this phenotype.

Population pharmacokinetics in chronic myeloid leukaemia patients: observations under field-conditions

Introduction: Therapeutic drug monitoring (TDM) of imatinib has been increasingly proposed for chronic myeloid leukaemia (CML) patients, as several studies have found a correlation between trough concentrations (Cmin) >=1000ng/ml and improved response. The pharmacological monitoring project of EUTOS (European Treatment and Outcome Study) was launched to increase the availability of imatinib TDM, standardize labs, and validate proposed Cmin thresholds. Using the collected data, the objective of this analysis was to characterize imatinib Population pharmacokinetics (Pop-PK) in a large cohort of European patients, to quantify its variability and the influence of demographic factors and comedications, and to derive individual exposure variables suitable for further concentration-effect analyses.¦Methods: 4095 PK samples from 2478 adult patients were analyzed between 2006 and 2010 by LC-MS-MS and considered for Pop-PK analysis by NONMEM®. Model building used data from 973 patients with >=2 samples available (2590 samples). A sensitivity analysis was performed using all data. Available comedications (27%) were classified into inducers or inhibitors of P-glycoprotein, CYP3A4/5 and organic-cation-transporter-1 (hOCT-1).¦Results: A one-compartment model with linear elimination, zero-order absorption fitted the data best. Estimated Pop-PK parameters (interindividual variability, IIV %CV) for a 40-year old male patient were: clearance CL = 17.3 L/h (37.7%), volume V = 429L (51.1%), duration of absorption D1 = 3.2h. Outliers, reflecting potential compliance and time recording errors, were taken into account by estimating an IIV on the residual error (35.4%). Intra-individual residuals were 29.1% (proportional) plus ± 84.6 ng/mL (additive). Female patients had a 15.2% lower CL (14.6 L/h). A piece-wise linear effect of age estimated a CL of 18.7 L/h at 20 years, 17.3 L/h at 40 and 13.8 L/h at 60 years. These covariates explained 2% (CL) and 4.5% (V) of IIV variability. No effect of comedication was found. The sensitivity analysis expectedly estimated increased IIV, but similar fixed effect parameters.¦Conclusion: Imatinib PK was well described in a large cohort of CML patients under field conditions and results were concordant with previous studies. Patient characteristics explain only little IIV, confirming limited utility of prior dosage adjustment. As intra-variability is smaller than inter-patient variability, dose adjustment guided by TDM could however be beneficial in order to bring Cmin into a given therapeutic target.

Genetic improvement and population structure of the Nelore breed in Northern Brazil

The objective of this work was to evaluate the population structure and the genetic and phenotypic progress of Nelore cattle in Northern Brazil. Pedigree information concerning animals born between 1942 and 2006 were analyzed. Population structure was performed using the Endog program. Out of the 140,628 animals studied, 67.7, 14.52 and 3.18% had complete pedigree record of the first, second and third parental generation, respectively. Inbreeding and average relatedness coefficients were low: 0.2 and 0.13%, respectively. However, these parameters may have been underestimated, since information on pedigree was incomplete. The effective number of founders was 370 and the genetic contribution of 10, 50 and 448 most influent ancestors explained 13.2, 28 and 50% of the genetic variability in the population, respectively. The genetic variability for growth traits and population structure demonstrates high probability of increasing productivity through selective breeding. Moreover, management strategies to reduce the currently observed age at first calving and generation intervals are important for Nelore cattle genetic improvement.

Bayesian genetic parameters for body weight and survival of Nile tilapia farmed in Brazil

The objective of this study was to estimate genetic parameters for survival and weight of Nile tilapia (Oreochromis niloticus), farmed in cages and ponds in Brazil, and to predict genetic gain under different scenarios. Survival was recorded as a binary response (dead or alive), during harvest time in the 2008 grow-out period. Genetic parameters were estimated using a Bayesian mixed linear-threshold animal model via Gibbs sampling. The breeding population consisted of 2,912 individual fish, which were analyzed together with the pedigree of 5,394 fish. The heritabilities estimates, with 95% posterior credible intervals, for tagging weight, harvest weight and survival were 0.17 (0.09-0.27), 0.21 (0.12-0.32) and 0.32 (0.22-0.44), respectively. Credible intervals show a 95% probability that the true genetic correlations were in a favourable direction. The selection for weight has a positive impact on survival. Estimated genetic gain was high when selecting for harvest weight (5.07%), and indirect gain for tagging weight (2.17%) and survival (2.03%) were also considerable.

Detection and mapping of a lethal locus in a eucalyptus hybrid population

The objective of this work was to verify the existence of a lethal locus in a eucalyptus hybrid population, and to quantify the segregation distortion in the linkage group 3 of the Eucalyptus genome. A E. grandis x E. urophylla hybrid population, which segregates for rust resistance, was genotyped with 19 microsatellite markers belonging to linkage group 3 of the Eucalyptus genome. To quantify the segregation distortion, maximum likelihood (ML) models, specific to outbreeding populations, were used. These models consider the observed marker genotypes and the lethal locus viability as parameters. The ML solutions were obtained using the expectation‑maximization algorithm. A lethal locus in the linkage group 3 was verified and mapped, with high confidence, between the microssatellites EMBRA 189 e EMBRA 122. This lethal locus causes an intense gametic selection from the male side. Its map position is 25 cM from the locus which controls the rust resistance in this population.

Genetic parameters and estimated genetic gains in young rubber tree progenies

The objective of this work was to assess the genetic parameters and to estimate genetic gains in young rubber tree progenies. The experiments were carried out during three years, in a randomized block design, with six replicates and ten plants per plot, in three representative Hevea crop regions of the state of São Paulo, Brazil. Twenty-two progenies were evaluated, from three to five years old, for rubber yield and annual girth growth. Genetic gain was estimated with the multi-effect index (MEI). Selection by progenies means provided greater estimated genetic gain than selection based on individuals, since heritability values of progeny means were greater than the ones of individual heritability, for both evaluated variables, in all the assessment years. The selection of the three best progenies for rubber yield provided a selection gain of 1.28 g per plant. The genetic gains estimated with MEI using data from early assessments (from 3 to 5-year-old) were generally high for annual girth growth and rubber yield. The high genetic gains for annual girth growth in the first year of assessment indicate that progenies can be selected at the beginning of the breeding program. Population effective size was consistent with the three progenies selected, showing that they were not related and that the population genetic variability is ensured. Early selection with the genetic gains estimated by MEI can be made on rubber tree progenies.

Ambiguous nucleotide calls from population-based sequencing of HIV-1 are a marker for viral diversity and the age of infection.

Background. The time passed since the infection of a human immunodeficiency virus (HIV)-infected individual (the age of infection) is an important but often only poorly known quantity. We assessed whether the fraction of ambiguous nucleotides obtained from bulk sequencing as done for genotypic resistance testing can serve as a proxy of this parameter. Methods. We correlated the age of infection and the fraction of ambiguous nucleotides in partial pol sequences of HIV-1 sampled before initiation of antiretroviral therapy (ART). Three groups of Swiss HIV Cohort Study participants were analyzed, for whom the age of infection was estimated on the basis of Bayesian back calculation (n = 3,307), seroconversion (n = 366), or diagnoses of primary HIV infection (n = 130). In addition, we studied 124 patients for whom longitudinal genotypic resistance testing was performed while they were still ART-naive. Results. We found that the fraction of ambiguous nucleotides increased with the age of infection with a rate of .2% per year within the first 8 years but thereafter with a decreasing rate. We show that this pattern is consistent with population-genetic models for realistic parameters. Finally, we show that, in this highly representative population, a fraction of ambiguous nucleotides of >.5% provides strong evidence against a recent infection event < 1 year prior to sampling (negative predictive value, 98.7%). Conclusions. These findings show that the fraction of ambiguous nucleotides is a useful marker for the age of infection.

Seed biometric parameters in oil palm accessions from a Brazilian germplasm bank

The objective of this work was to evaluate the morphological diversity of oil palm seeds and to cluster the accessions according to their morphological characteristics. Forty-one accessions from the oil palm germplasm bank of Embrapa Amazônia Ocidental were evaluated - 18 of Elaeis oleifera and 23 of E. guineensis. The groups were formed based on morphological characteristics, by principal component analysis. In E. oleifera, four groups were formed, tied to their region of origin, but with significant morphological differences between accessions from the same population. For tenera-type E. guineensis seeds, three widely divergent groups were formed, especially as to external parameters, which differentiated them from the other ones. The parameter endocarp thickness stood out in intra- and inter-population differentiation. For dura-type E. guineensis, three groups were formed, with larger seeds and thicker endocarps, which differed from all the other ones. The variability observed for seed characteristics in the analyzed accessions allows the establishment of different groups, to define strategies for genetic improvement.

Linking landscape- and population ecology for large population management modelling: the case of Ibex (capra ibex) in Switzerland

La biologie de la conservation est communément associée à la protection de petites populations menacées d?extinction. Pourtant, il peut également être nécessaire de soumettre à gestion des populations surabondantes ou susceptibles d?une trop grande expansion, dans le but de prévenir les effets néfastes de la surpopulation. Du fait des différences tant quantitatives que qualitatives entre protection des petites populations et contrôle des grandes, il est nécessaire de disposer de modèles et de méthodes distinctes. L?objectif de ce travail a été de développer des modèles prédictifs de la dynamique des grandes populations, ainsi que des logiciels permettant de calculer les paramètres de ces modèles et de tester des scénarios de gestion. Le cas du Bouquetin des Alpes (Capra ibex ibex) - en forte expansion en Suisse depuis sa réintroduction au début du XXème siècle - servit d?exemple. Cette tâche fut accomplie en trois étapes : En premier lieu, un modèle de dynamique locale, spécifique au Bouquetin, fut développé : le modèle sous-jacent - structuré en classes d?âge et de sexe - est basé sur une matrice de Leslie à laquelle ont été ajoutées la densité-dépendance, la stochasticité environnementale et la chasse de régulation. Ce modèle fut implémenté dans un logiciel d?aide à la gestion - nommé SIM-Ibex - permettant la maintenance de données de recensements, l?estimation automatisée des paramètres, ainsi que l?ajustement et la simulation de stratégies de régulation. Mais la dynamique d?une population est influencée non seulement par des facteurs démographiques, mais aussi par la dispersion et la colonisation de nouveaux espaces. Il est donc nécessaire de pouvoir modéliser tant la qualité de l?habitat que les obstacles à la dispersion. Une collection de logiciels - nommée Biomapper - fut donc développée. Son module central est basé sur l?Analyse Factorielle de la Niche Ecologique (ENFA) dont le principe est de calculer des facteurs de marginalité et de spécialisation de la niche écologique à partir de prédicteurs environnementaux et de données d?observation de l?espèce. Tous les modules de Biomapper sont liés aux Systèmes d?Information Géographiques (SIG) ; ils couvrent toutes les opérations d?importation des données, préparation des prédicteurs, ENFA et calcul de la carte de qualité d?habitat, validation et traitement des résultats ; un module permet également de cartographier les barrières et les corridors de dispersion. Le domaine d?application de l?ENFA fut exploré par le biais d?une distribution d?espèce virtuelle. La comparaison à une méthode couramment utilisée pour construire des cartes de qualité d?habitat, le Modèle Linéaire Généralisé (GLM), montra qu?elle était particulièrement adaptée pour les espèces cryptiques ou en cours d?expansion. Les informations sur la démographie et le paysage furent finalement fusionnées en un modèle global. Une approche basée sur un automate cellulaire fut choisie, tant pour satisfaire aux contraintes du réalisme de la modélisation du paysage qu?à celles imposées par les grandes populations : la zone d?étude est modélisée par un pavage de cellules hexagonales, chacune caractérisée par des propriétés - une capacité de soutien et six taux d?imperméabilité quantifiant les échanges entre cellules adjacentes - et une variable, la densité de la population. Cette dernière varie en fonction de la reproduction et de la survie locale, ainsi que de la dispersion, sous l?influence de la densité-dépendance et de la stochasticité. Un logiciel - nommé HexaSpace - fut développé pour accomplir deux fonctions : 1° Calibrer l?automate sur la base de modèles de dynamique (par ex. calculés par SIM-Ibex) et d?une carte de qualité d?habitat (par ex. calculée par Biomapper). 2° Faire tourner des simulations. Il permet d?étudier l?expansion d?une espèce envahisseuse dans un paysage complexe composé de zones de qualité diverses et comportant des obstacles à la dispersion. Ce modèle fut appliqué à l?histoire de la réintroduction du Bouquetin dans les Alpes bernoises (Suisse). SIM-Ibex est actuellement utilisé par les gestionnaires de la faune et par les inspecteurs du gouvernement pour préparer et contrôler les plans de tir. Biomapper a été appliqué à plusieurs espèces (tant végétales qu?animales) à travers le Monde. De même, même si HexaSpace fut initialement conçu pour des espèces animales terrestres, il pourrait aisément être étndu à la propagation de plantes ou à la dispersion d?animaux volants. Ces logiciels étant conçus pour, à partir de données brutes, construire un modèle réaliste complexe, et du fait qu?ils sont dotés d?une interface d?utilisation intuitive, ils sont susceptibles de nombreuses applications en biologie de la conservation. En outre, ces approches peuvent également s?appliquer à des questions théoriques dans les domaines de l?écologie des populations et du paysage.<br/><br/>Conservation biology is commonly associated to small and endangered population protection. Nevertheless, large or potentially large populations may also need human management to prevent negative effects of overpopulation. As there are both qualitative and quantitative differences between small population protection and large population controlling, distinct methods and models are needed. The aim of this work was to develop theoretical models to predict large population dynamics, as well as computer tools to assess the parameters of these models and to test management scenarios. The alpine Ibex (Capra ibex ibex) - which experienced a spectacular increase since its reintroduction in Switzerland at the beginning of the 20th century - was used as paradigm species. This task was achieved in three steps: A local population dynamics model was first developed specifically for Ibex: the underlying age- and sex-structured model is based on a Leslie matrix approach with addition of density-dependence, environmental stochasticity and culling. This model was implemented into a management-support software - named SIM-Ibex - allowing census data maintenance, parameter automated assessment and culling strategies tuning and simulating. However population dynamics is driven not only by demographic factors, but also by dispersal and colonisation of new areas. Habitat suitability and obstacles modelling had therefore to be addressed. Thus, a software package - named Biomapper - was developed. Its central module is based on the Ecological Niche Factor Analysis (ENFA) whose principle is to compute niche marginality and specialisation factors from a set of environmental predictors and species presence data. All Biomapper modules are linked to Geographic Information Systems (GIS); they cover all operations of data importation, predictor preparation, ENFA and habitat suitability map computation, results validation and further processing; a module also allows mapping of dispersal barriers and corridors. ENFA application domain was then explored by means of a simulated species distribution. It was compared to a common habitat suitability assessing method, the Generalised Linear Model (GLM), and was proven better suited for spreading or cryptic species. Demography and landscape informations were finally merged into a global model. To cope with landscape realism and technical constraints of large population modelling, a cellular automaton approach was chosen: the study area is modelled by a lattice of hexagonal cells, each one characterised by a few fixed properties - a carrying capacity and six impermeability rates quantifying exchanges between adjacent cells - and one variable, population density. The later varies according to local reproduction/survival and dispersal dynamics, modified by density-dependence and stochasticity. A software - named HexaSpace - was developed, which achieves two functions: 1° Calibrating the automaton on the base of local population dynamics models (e.g., computed by SIM-Ibex) and a habitat suitability map (e.g. computed by Biomapper). 2° Running simulations. It allows studying the spreading of an invading species across a complex landscape made of variously suitable areas and dispersal barriers. This model was applied to the history of Ibex reintroduction in Bernese Alps (Switzerland). SIM-Ibex is now used by governmental wildlife managers to prepare and verify culling plans. Biomapper has been applied to several species (both plants and animals) all around the World. In the same way, whilst HexaSpace was originally designed for terrestrial animal species, it could be easily extended to model plant propagation or flying animals dispersal. As these softwares were designed to proceed from low-level data to build a complex realistic model and as they benefit from an intuitive user-interface, they may have many conservation applications. Moreover, theoretical questions in the fields of population and landscape ecology might also be addressed by these approaches.