Human and climatic change impact modelling on the habitat suitability for the chimpanzee (Pan troglodytes ellioti) – Case study: The proposed Mount Cameroon National Park

Dissertation submitted in partial fulfilment of the requirements for the Degree of Master of Science in Geospatial Technologies

Predicting the distribution of forest tree species using topographic variables and vegetation index in eastern Acre, Brazil

Species distribution modeling has relevant implications for the studies of biodiversity, decision making about conservation and knowledge about ecological requirements of the species. The aim of this study was to evaluate if the use of forest inventories can improve the estimation of occurrence probability, identify the limits of the potential distribution and habitat preference of a group of timber tree species. The environmental predictor variables were: elevation, slope, aspect, normalized difference vegetation index (NDVI) and height above the nearest drainage (HAND). To estimate the distribution of species we used the maximum entropy method (Maxent). In comparison with a random distribution, using topographic variables and vegetation index as features, the Maxent method predicted with an average accuracy of 86% the geographical distribution of studied species. The altitude and NDVI were the most important variables. There were limitations to the interpolation of the models for non-sampled locations and that are outside of the elevation gradient associated with the occurrence data in approximately 7% of the basin area. Ceiba pentandra (samaúma), Castilla ulei (caucho) and Hura crepitans (assacu) is more likely to occur in nearby water course areas. Clarisia racemosa (guariúba), Amburana acreana (cerejeira), Aspidosperma macrocarpon (pereiro), Apuleia leiocarpa (cumaru cetim), Aspidosperma parvifolium (amarelão) and Astronium lecointei (aroeira) can also occur in upland forest and well drained soils. This modeling approach has potential for application on other tropical species still less studied, especially those that are under pressure from logging.

Brazilian Flavivirus phylogeny based on NS5

In this work, a comprehensive phylogenetic study based on 600 base pair nucleotide and on putative 200 amino acid sequences of NS5 was carried out in order to establish genetic relationships among 15 strains of 10 Brazilian flaviviruses: Bussuquara, Cacipacore, dengue type 1, 2 and 4, Iguape, Ilheus, Rocio, Saint Louis encephalitis (SLE), and yellow fever. Phylogenetic trees were created by neighbor-joining and maximum parsimony methods. These trees showed Brazilian flaviviruses grouped into three main branches: yellow fever branch, dengue branch subdivided in types 1, 2 and 4 branches, and Japanese encephalitis virus (JEV) complex branch including SLE virus strains, Cacipacore, Iguape, Rocio, Ilheus and Bussuquara. Viruses transmitted by Aedes mosquitoes, such as dengue and urban yellow fever, that are also the only Flavivirus causing hemorrhagic fevers in Brazil, were grouped in the same cluster. Encephalitis associated viruses, transmitted by Culex mosquitoes such as JEV complex branch including SLE virus strains, Cacipacore, Iguape, Rocio, Ilheus and Bussuquara were also grouped in the same clade.

Nava enhances ventilatory variability and diaphragmatic activity/tidal volume coupling

INTRODUCTION. Neurally Adjusted Ventilatory Assist (NAVA) is a new ventilatory mode in which ventilator settings are adjusted based on the electrical activity detected in the diaphragm (Eadi). This mode offers significant advantages in mechanical ventilation over standard pressure support (PS) modes, since ventilator input is determined directly from patient ventilatory demand. Therefore, it is expected that tidal volume (Vt) under NAVA would show better correlation with Eadi compared with PS, and exhibit greater variability due to the variability in the Eadi input to the ventilator. OBJECTIVES. To compare tidal volume variability in PS and NAVA ventilation modes, and its correlation with patient ventilatory demand (as characterized by maximum Eadi). METHODS. Acomparative study of patient-ventilator interaction was performed for 22 patients during standard PS with clinician determined ventilator settings; and NAVA, with NAVA gain set to ensure the same peak airway pressure as the total pressure obtained in PS. A 20 min continuous recording was performed in each ventilator mode. Respiratory rate, Vt, and Eadi were recorded. Tidal volume variance and Pearson correlation coefficient between Vt and Eadi were calculated for each patient. A periodogram was plotted for each ventilator mode and each patient, showing spectral power as a function of frequency to assess variability. RESULTS. Median, lower quartile and upper quartile values for Vt variance and Vt/Eadi correlation are shown in Table 1. The NAVA cohort exhibits substantially greater correlation and variance than the PS cohort. Power spectrums for Vt and Eadi are shown in Fig. 1 (PS and NAVA) for a typical patient. The enlarged section highlights how changes in Eadi are highly synchronized with NAVA ventilation, but less so for PS. CONCLUSIONS. There is greater variability in tidal volume and correlation between tidal volume and diaphragmatic electrical activity with NAVA compared to PS. These results are consistent with the improved patient-ventilator synchrony reported in the literature.

Variation and evolution of herkogamy in Exochaenium (Gentianaceae): implications for the evolution of distyly.

Backgrounds and Aims The spatial separation of stigmas and anthers (herkogamy) in flowering plants functions to reduce self-pollination and avoid interference between pollen dispersal and receipt. Little is known about the evolutionary relationships among the three main forms of herkogamy - approach, reverse and reciprocal herkogamy (distyly) - or about transitions to and from a non-herkogamous condition. This problem was examined in Exochaenium (Gentianaceae), a genus of African herbs that exhibits considerable variation in floral morphology, including the three forms of herkogamy. Methods Using maximum parsimony and maximum likelihood methods, the evolutionary history of herkogamic and non-herkogamic conditions was reconstructed from a molecular phylogeny of 15 species of Exochaenium and four outgroup taxa, based on three chloroplast regions, the nuclear ribosomal internal transcribed spacer (ITS1 and 2) and the 5·8S gene. Ancestral character states were determined and the reconstructions were used to evaluate competing models for the origin of reciprocal herkogamy. Key results Reciprocal herkogamy originated once in Exochaenium from an ancestor with approach herkogamy. Reverse herkogamy and the non-herkogamic condition homostyly were derived from heterostyly. Distylous species possessed pendent, slightly zygomorphic flowers, and the single transition to reverse herkogamy was associated with the hawkmoth pollination syndrome. Reductions in flower size characterized three of four independent transitions from reciprocal herkogamy to homostyly. Conclusions The results support Lloyd and Webb's model in which distyly originated from an ancestor with approach herkogamy. They also demonstrate the lability of sex organ deployment and implicate pollinators, or their absence, as playing an important role in driving transitions among herkogamic and non-herkogamic conditions.

The serine repeat antigen (SERA) gene family phylogeny in Plasmodium: the impact of GC content and reconciliation of gene and species trees.

Plasmodium falciparum is the parasite responsible for the most acute form of malaria in humans. Recently, the serine repeat antigen (SERA) in P. falciparum has attracted attention as a potential vaccine and drug target, and it has been shown to be a member of a large gene family. To clarify the relationships among the numerous P. falciparum SERAs and to identify orthologs to SERA5 and SERA6 in Plasmodium species affecting rodents, gene trees were inferred from nucleotide and amino acid sequence data for 33 putative SERA homologs in seven different species. (A distance method for nucleotide sequences that is specifically designed to accommodate differing GC content yielded results that were largely compatible with the amino acid tree. Standard-distance and maximum-likelihood methods for nucleotide sequences, on the other hand, yielded gene trees that differed in important respects.) To infer the pattern of duplication, speciation, and gene loss events in the SERA gene family history, the resulting gene trees were then "reconciled" with two competing Plasmodium species tree topologies that have been identified by previous phylogenetic studies. Parsimony of reconciliation was used as a criterion for selecting a gene tree/species tree pair and provided (1) support for one of the two species trees and for the core topology of the amino acid-derived gene tree, (2) a basis for critiquing fine detail in a poorly resolved region of the gene tree, (3) a set of predicted "missing genes" in some species, (4) clarification of the relationship among the P. falciparum SERA, and (5) some information about SERA5 and SERA6 orthologs in the rodent malaria parasites. Parsimony of reconciliation and a second criterion--implied mutational pattern at two key active sites in the SERA proteins-were also seen to be useful supplements to standard "bootstrap" analysis for inferred topologies.

Mapping the molecular characteristics of Brazilian human T-cell lymphotropic virus type 1 Env (gp46) and Pol amino acid sequences for vaccine design

This study was carried out to evaluate the molecular pattern of all available Brazilian human T-cell lymphotropic virus type 1 Env (n = 15) and Pol (n = 43) nucleotide sequences via epitope prediction, physico-chemical analysis, and protein potential sites identification, giving support to the Brazilian AIDS vaccine program. In 12 previously described peptides of the Env sequences we found 12 epitopes, while in 4 peptides of the Pol sequences we found 4 epitopes. The total variation on the amino acid composition was 9 and 17% for human leukocyte antigen (HLA) class I and class II Env epitopes, respectively. After analyzing the Pol sequences, results revealed a total amino acid variation of 0.75% for HLA-I and HLA-II epitopes. In 5 of the 12 Env epitopes the physico-chemical analysis demonstrated that the mutations magnified the antigenicity profile. The potential protein domain analysis of Env sequences showed the loss of a CK-2 phosphorylation site caused by D197N mutation in one epitope, and a N-glycosylation site caused by S246Y and V247I mutations in another epitope. Besides, the analysis of selection pressure have found 8 positive selected sites (w = 9.59) using the codon-based substitution models and maximum-likelihood methods. These studies underscore the importance of this Env region for the virus fitness, for the host immune response and, therefore, for the development of vaccine candidates.

Influence of environmental conditions on the distribution of Central Asian green toads with three ploidy levels.

We studied the distribution of Palearctic green toads (Bufo viridis subgroup), an anuran species group with three ploidy levels, inhabiting the Central Asian Amudarya River drainage. Various approaches (one-way, multivariate, components variance analyses and maximum entropy modelling) were used to estimate the effect of altitude, precipitation, temperature and land vegetation covers on the distribution of toads. It is usually assumed that polyploid species occur in regions with harsher climatic conditions (higher latitudes, elevations, etc.), but for the green toads complex, we revealed a more intricate situation. The diploid species (Bufo shaartusiensis and Bufo turanensis) inhabit the arid lowlands (from 44 to 789 m a.s.l.), while tetraploid Bufo pewzowi were recorded in mountainous regions (340-3492 m a.s.l.) with usually lower temperatures and higher precipitation rates than in the region inhabited by diploid species. The triploid species Bufo baturae was found in the Pamirs (Tajikistan) at the highest altitudes (2503-3859 m a.s.l.) under the harshest climatic conditions.

Climatic factors influencing triatomine occurrence in Central-West Brazil

We estimated the geographic distributions of triatomine species in Central-West Region of Brazil (CW) and analysed the climatic factors influencing their occurrence. A total of 3,396 records of 27 triatomine species were analysed. Using the maximum entropy method, ecological niche models were produced for eight species occurring in at least 20 municipalities based on 13 climatic variables and elevation. Triatoma sordida and Rhodnius neglectus were the species with the broadest geographic distributions in CW Brazil. The Cerrado areas in the state of Goiás were found to be more suitable for the occurrence of synanthropic triatomines than the Amazon forest areas in the northern part of the state of Mato Grosso. The variable that best explains the evaluated models is temperature seasonality. The results indicate that almost the entire region presents climatic conditions that are appropriate for at least one triatomine species. Therefore, it is recommended that entomological surveillance be reinforced in CW Brazil.

Predicting the geographic distribution of Lutzomyia longipalpis (Diptera: Psychodidae) and visceral leishmaniasis in the state of Mato Grosso do Sul, Brazil

To understand the geographic distribution of visceral leishmaniasis (VL) in the state of Mato Grosso do Sul (MS), Brazil, both the climatic niches of Lutzomyia longipalpis and VL cases were analysed. Distributional data were obtained from 55 of the 79 counties of MS between 2003-2012. Ecological niche models (ENM) of Lu. longipalpis and VL cases were produced using the maximum entropy algorithm based on eight climatic variables. Lu. longipalpis showed a wide distribution in MS. The highest climatic suitability for Lu. longipalpis was observed in southern MS. Temperature seasonality and annual mean precipitation were the variables that most influenced these models. Two areas of high climatic suitability for the occurrence of VL cases were predicted: one near Aquidauana and another encompassing several municipalities in the southeast region of MS. As expected, a large overlap between the models for Lu. longipalpis and VL cases was detected. Northern and northwestern areas of MS were suitable for the occurrence of cases, but did not show high climatic suitability for Lu. longipalpis . ENM of vectors and human cases provided a greater understanding of the geographic distribution of VL in MS, which can be applied to the development of future surveillance strategies.

On the Properties and Cosmology of Q-balls

The properties and cosmological importance of a class of non-topological solitons, Q-balls, are studied. Aspects of Q-ball solutions and Q-ball cosmology discussed in the literature are reviewed. Q-balls are particularly considered in the Minimal Supersymmetric Standard Model with supersymmetry broken by a hidden sector mechanism mediated by either gravity or gauge interactions. Q-ball profiles, charge-energy relations and evaporation rates for realistic Q-ball profiles are calculated for general polynomial potentials and for the gravity mediated scenario. In all of the cases, the evaporation rates are found to increase with decreasing charge. Q-ball collisions are studied by numerical means in the two supersymmetry breaking scenarios. It is noted that the collision processes can be divided into three types: fusion, charge transfer and elastic scattering. Cross-sections are calculated for the different types of processes in the different scenarios. The formation of Q-balls from the fragmentation of the Aflieck-Dine -condensate is studied by numerical and analytical means. The charge distribution is found to depend strongly on the initial energy-charge ratio of the condensate. The final state is typically noted to consist of Q- and anti-Q-balls in a state of maximum entropy. By studying the relaxation of excited Q-balls the rate at which excess energy can be emitted is calculated in the gravity mediated scenario. The Q-ball is also found to withstand excess energy well without significant charge loss. The possible cosmological consequences of these Q-ball properties are discussed.

Evidence for adaptive radiation from a phylogenetic study of plant defenses.

One signature of adaptive radiation is a high level of trait change early during the diversification process and a plateau toward the end of the radiation. Although the study of the tempo of evolution has historically been the domain of paleontologists, recently developed phylogenetic tools allow for the rigorous examination of trait evolution in a tremendous diversity of organisms. Enemy-driven adaptive radiation was a key prediction of Ehrlich and Raven's coevolutionary hypothesis [Ehrlich PR, Raven PH (1964) Evolution 18:586-608], yet has remained largely untested. Here we examine patterns of trait evolution in 51 North American milkweed species (Asclepias), using maximum likelihood methods. We study 7 traits of the milkweeds, ranging from seed size and foliar physiological traits to defense traits (cardenolides, latex, and trichomes) previously shown to impact herbivores, including the monarch butterfly. We compare the fit of simple random-walk models of trait evolution to models that incorporate stabilizing selection (Ornstein-Ulenbeck process), as well as time-varying rates of trait evolution. Early bursts of trait evolution were implicated for 2 traits, while stabilizing selection was implicated for several others. We further modeled the relationship between trait change and species diversification while allowing rates of trait evolution to vary during the radiation. Species-rich lineages underwent a proportionately greater decline in latex and cardenolides relative to species-poor lineages, and the rate of trait change was most rapid early in the radiation. An interpretation of this result is that reduced investment in defensive traits accelerated diversification, and disproportionately so, early in the adaptive radiation of milkweeds.

Impact of model complexity on cross-temporal transferability in Maxent species distribution models: An assessment using paleobotanical data

Maximum entropy modeling (Maxent) is a widely used algorithm for predicting species distributions across space and time. Properly assessing the uncertainty in such predictions is non-trivial and requires validation with independent datasets. Notably, model complexity (number of model parameters) remains a major concern in relation to overfitting and, hence, transferability of Maxent models. An emerging approach is to validate the cross-temporal transferability of model predictions using paleoecological data. In this study, we assess the effect of model complexity on the performance of Maxent projections across time using two European plant species (Alnus giutinosa (L.) Gaertn. and Corylus avellana L) with an extensive late Quaternary fossil record in Spain as a study case. We fit 110 models with different levels of complexity under present time and tested model performance using AUC (area under the receiver operating characteristic curve) and AlCc (corrected Akaike Information Criterion) through the standard procedure of randomly partitioning current occurrence data. We then compared these results to an independent validation by projecting the models to mid-Holocene (6000 years before present) climatic conditions in Spain to assess their ability to predict fossil pollen presence-absence and abundance. We find that calibrating Maxent models with default settings result in the generation of overly complex models. While model performance increased with model complexity when predicting current distributions, it was higher with intermediate complexity when predicting mid-Holocene distributions. Hence, models of intermediate complexity resulted in the best trade-off to predict species distributions across time. Reliable temporal model transferability is especially relevant for forecasting species distributions under future climate change. Consequently, species-specific model tuning should be used to find the best modeling settings to control for complexity, notably with paleoecological data to independently validate model projections. For cross-temporal projections of species distributions for which paleoecological data is not available, models of intermediate complexity should be selected.