Dengue and the risk of urban yellow fever reintroduction in São Paulo State, Brazil

OBJECTIVE: To propose a mathematical method for the estimation of the Basic Reproduction Number, R0, of urban yellow fever in a dengue-infested area. METHODS: The method is based on the assumption that, as the same vector (Aedes aegypti) causes both infections, all the quantities related to the mosquito, estimated from the initial phase of dengue epidemic, could be applied to yellow fever dynamics. It is demonstrated that R0 for yellow fever is, on average, 43% lower than that for dengue. This difference is due to the longer dengue viremia and its shorter extrinsic incubation period. RESULTS: In this study the analysis was expanded to the epidemiological situation of dengue in São Paulo in the year 2001. The total number of dengue cases increased from 3,582 in 2000 to 51,348 in 2001. It was then calculated R0 for yellow fever for every city which have shown R0 of dengue greater than 1. It was also estimated the total number of unprotected people living in highly risky areas for urban yellow fever. CONCLUSIONS: Currently there is a great number of non-vaccinated people living in Aedes aegypti infested area in the state of São Paulo.

Pertussis booster vaccine for adolescents and young adults in São Paulo, Brazil

OBJECTIVE: To develop a model to assess different strategies of pertussis booster vaccination in the city of São Paulo. METHODS: A dynamic stationary age-dependent compartmental model with waning immunity was developed. The "Who Acquires Infection from Whom" matrix was used to modeling age-dependent transmission rates. There were tested different strategies including vaccine boosters to the current vaccination schedule and three of them were reported: (i) 35% coverage at age 12, or (ii) 70% coverage at age 12, and (iii) 35% coverage at age 12 and 70% coverage at age 20 at the same time. RESULTS: The strategy (i) achieved a 59% reduction of pertussis occurrence and a 53% reduction in infants while strategy (ii) produced 76% and 63% reduction and strategy (iii) 62% and 54%, respectively. CONCLUSION: Pertussis booster vaccination at age 12 proved to be the best strategy among those tested in this study as it achieves the highest overall reduction and the greatest impact among infants who are more susceptible to pertussis complications.

Bases para o estudo da genética de populações dos Hymenoptera em geral e dos Apinae sociais em particular

This paper deals with problems on population genetics in Hymenoptera and particularly in social Apidae. 1) The studies on populations of Hymenoptera were made according to the two basic types of reproduction: endogamy and panmixia. The populations of social Apinae have a mixed method of reproduction with higher percentage of panmixia and a lower of endogamy. This is shown by the following a) males can enter any hive in swarming time; b) males of Meliponini are expelled from hives which does not need them, and thus, are forced to look for some other place; c) Meliponini males were seen powdering themselves with pollen, thus becoming more acceptable in any other hive. The panmixia is not complete owing to the fact that the density of the breeding population as very low, even in the more frequent species as low as about 2 females and 160 males per reproductive area. We adopted as selection values (or survival indices) the expressions according to Brieger (1948,1950) which may be summarised as follows; a population: p2AA + ²pq Aa + q2aa became after selection: x p2AA + 2pq Aa + z q²aa. For alge-braics facilities Brieger divided the three selective values by y giving thus: x/y p2 AA + y/y 2 pq Aa + z/y q²aa. He called x/y of RA and z/y of Ra, that are survival or selective index, calculated in relation to the heterozygote. In our case all index were calculated in relation to the heterozygote, including the ones for haploid males; thus we have: RA surveval index of genotype AA Ra surveval index of genotype aa R'A surveval index of genotype A R'a surveval index of genotype a 1 surveval index of genotype Aa The index R'A ande R'a were equalized to RA and Ra, respectively, for facilities in the conclusions. 2) Panmitic populations of Hymenoptera, barring mutations, migrations and selection, should follow the Hardy-Weinberg law, thus all gens will be present in the population in the inicial frequency (see Graphifc 1). 3) Heterotic genes: If mutation for heterotic gene ( 1 > RA > Ra) occurs, an equilibrium will be reached in a population when: P = R A + Ra - 2R²a _____________ (9) 2(R A + Ra - R²A - R²a q = R A + Ra - 2R²A _____________ (10) 2(R A + Ra - R²A - R²a A heterotic gene in an hymenopteran population may be maintained without the aid of new mutation only if the survival index of the most viable mutant (RA) does not exced the limiting value given by the formula: R A = 1 + √1+Ra _________ 4 If RA has a value higher thah the one permitted by the formula, then only the more viable gene will remain present in the population (see Graphic 10). The only direct proof for heterotic genes in Hymenoptera was given by Mackensen and Roberts, who obtained offspring from Apis mellefera L. queens fertilized by their own sons. Such inbreeding resulted in a rapid loss of vigor the colony; inbred lines intercrossed gave a high hybrid vigor. Other fats correlated with the "heterosis" problem are; a) In a colony M. quadrifasciata Lep., which suffered severely from heat, the percentage of deths omong males was greater .than among females; b) Casteel and Phillips had shown that in their samples (Apis melifera L). the males had 7 times more abnormalities tian the workers (see Quadros IV to VIII); c) just after emerging the males have great variation, but the older ones show a variation equal to that of workers; d) The tongue lenght of males of Apis mellifera L., of Bombus rubicundus Smith (Quadro X), of Melipona marginata Lep. (Quadro XI), and of Melipona quadrifasciata Lep. Quadro IX, show greater variationthan that of workers of the respective species. If such variation were only caused by subviables genes a rapid increasse of homozigoty for the most viable alleles should be expected; then, these .wild populations, supposed to be in equilibrium, could .not show such variability among males. Thus we conclude that heterotic genes have a grat importance in these cases. 4) By means of mathematical models, we came to the conclusion tht isolating genes (Ra ^ Ra > 1), even in the case of mutations with more adaptability, have only the opor-tunity of survival when the population number is very low (thus the frequency of the gene in the breeding population will be large just after its appearence). A pair of such alleles can only remain present in a population when in border regions of two races or subspecies. For more details see Graphics 5 to 8. 5) Sex-limited genes affecting only females, are of great importance toHymenoptera, being subject to the same limits and formulas as diploid panmitic populations (see formulas 12 and 13). The following examples of these genes were given: a) caste-determining genes in the genus Melipona; b) genes permiting an easy response of females to differences in feeding in almost all social Hymenoptera; c) two genes, found in wild populations, one in Trigona (Plebéia) mosquito F. SMITH (quadro XII) and other in Melipona marginata marginata LEP. (Quadro XIII, colonies 76 and 56) showing sex-limited effects. Sex-limited genes affecting only males do not contribute to the plasticity or genie reserve in hymenopteran populations (see formula 14). 6) The factor time (life span) in Hymenoptera has a particular importance for heterotic genes. Supposing one year to be the time unit and a pair of heterotic genes with respective survival indice equal to RA = 0, 90 and Ra = 0,70 to be present; then if the life time of a population is either one or two years, only the more viable gene will remain present (see formula 11). If the species has a life time of three years, then both alleles will be maintained. Thus we conclude that in specis with long lif-time, the heterotic genes have more importance, and should be found more easily. 7) The colonies of social Hymenoptera behave as units in competition, thus in the studies of populations one must determine the survival index, of these units which may be subdivided in indice for egg-laying, for adaptive value of the queen, for working capacity of workers, etc. 8) A study of endogamic hymenopteran populations, reproduced by sister x brother mating (fig. 2), lead us to the following conclusions: a) without selection, a population, heterozygous for one pair of alleles, will consist after some generations (theoretically after an infinite number of generation) of females AA fecundated with males A and females aa fecundated with males a (see Quadro I). b) Even in endogamic population there is the theoretical possibility of the presence of heterotic genes, at equilibrium without the aid of new mutations (see Graphics 11 and 12), but the following! conditions must be satisfied: I - surveval index of both homozygotes (RA e Ra) should be below 0,75 (see Graphic 13); II - The most viable allele must riot exced the less viable one by more than is permited by the following formula (Pimentel Gomes 1950) (see Gra-fic 14) : 4 R5A + 8 Ra R4A - 4 Ra R³A (Ra - 1) R²A - - R²a (4 R²a + 4 Ra - 1) R A + 2 R³a < o Considering these two conditions, the existance of heterotic genes in endogamic populations of Hymenoptera \>ecames very improbable though not - impossible. 9) Genie mutation offects more hymenopteran than diploid populations. Thus we have for lethal genes in diploid populations: u = q2, and in Hymenoptera: u = s, being u the mutation ratio and s the frequency of the mutant in the male population. 10) Three factors, important to competition among species of Meliponini were analysed: flying capacity of workers, food gathering capacity of workers, egg-laying of the queen. In this connection we refer to the variability of the tongue lenght observed in colonies from several localites, to the method of transporting the pollen in the stomach, from some pots (Melliponi-ni storage alveolus) to others (e. g. in cases of pillage), and to the observation that the species with the most populous hives are almost always the most frequent ones also. 11) Several defensive ways used for Meliponini to avoid predation are cited, but special references are made upon the camouflage of both hive (fig. 5) and hive entrance (fig. 4) and on the mimetism (see list in page ). Also under the same heading we described the method of Lestrimelitta for pillage. 12) As mechanisms important for promoting genetic plasticity of hymenopteran species we cited: a) cytological variations and b) genie reserve. As to the former, duplications and numerical variations of chromosomes were studied. Diprion simile ATC was cited as example for polyploidy. Apis mellife-ra L. (n •= 16) also sugests polyploid origen since: a) The genus Melipona, which belongs to a" related tribe, presents in all species so far studied n = 9 chromosomes and b) there occurs formation of dyads in the firt spermatocyte division. It is su-gested that the origin of the sex-chromosome of Apis mellifera It. may be related to the possible origin of diplo-tetraploidy in this species. With regards to the genie reserve, several possible types of mutants were discussed. They were classified according to their survival indices; the heterotic and neutral mutants must be considered as more important for the genie reserve. 13) The mean radius from a mother to a daghter colony was estimated as 100 meters. Since the Meliponini hives swarm only once a year we may take 100 meters a year as the average dispersion of female Meliponini in ocordance to data obtained from Trigona (tetragonisca) jaty F. SMITH and Melipona marginata LEP., while other species may give different values. For males the flying distance was roughly estimated to be 10 times that for females. A review of the bibliography on Meliponini swarm was made (pg. 43 to 47) and new facts added. The population desity (breeding population) corresponds in may species of Meliponini to one male and one female per 10.000 square meters. Apparently the males are more frequent than the females, because there are sometimes many thousands, of males in a swarm; but for the genie frequency the individuals which have descendants are the ones computed. In the case of Apini and Meliponini, only one queen per hive and the males represented by. the spermatozoos in its spermateca are computed. In Meliponini only one male mate with the queen, while queens of Apis mellijera L. are fecundated by an average of about 1, 5 males. (Roberts, 1944). From the date cited, one clearly sees that, on the whole, populations of wild social bees (Meliponini) are so small that the Sewall Wright effect may become of great importance. In fact applying the Wright's formula: f = ( 1/aN♂ + 1/aN♀) (1 - 1/aN♂ + 1/aN♀) which measures the fixation and loss of genes per generation, we see that the fixation or loss of genes is of about 7% in the more frequent species, and rarer species about 11%. The variation in size, tergite color, background color, etc, of Melipona marginata Lep. is atributed to this genetic drift. A detail, important to the survival of Meliponini species, is the Constance of their breeding population. This Constance is due to the social organization, i. e., to the care given to the reproductive individuals (the queen with its sperm pack), to the way of swarming, to the food storage intended to control variations of feeding supply, etc. 14) Some species of the Meliponini are adapted to various ecological conditions and inhabit large geographical areas (e. g. T. (Tetragonisca jaty F. SMITH), and Trigona (Nanno-trigona testaceicornis LEP.) while others are limited to narrow regions with special ecological conditions (e. g. M. fuscata me-lanoventer SCHWARZ). Other species still, within the same geographical region, profit different ecological conditions, as do M. marginata LEP. and M. quadrifasciata LEP. The geographical distribution of Melipona quadrifasciata LEP. is different according to the subspecies: a) subsp anthidio-des LEP. (represented in Fig. 7 by black squares) inhabits a region fron the North of the S. Paulo State to Northeastern Brazil, ,b) subspecies quadrifasciata LEP., (marked in Fig. 7 with black triangles) accurs from the South of S. Paulo State to the middle of the State of Rio Grande do Sul (South Brazil). In the margined region between these two areas of distribution, hi-brid colonies were found (Fig. 7, white circles); they are shown with more details in fig. 8, while the zone of hybridization is roughly indicated in fig. 9 (gray zone). The subspecies quadrifasciata LEP., has 4 complete yellow bands on the abdominal tergites while anthidioides LEP. has interrupted ones. This character is determined by one or two genes and gives different adaptative properties to the subspecies. Figs. 10 shows certains meteorological isoclines which have aproximately the same configuration as the limits of the hybrid zone, suggesting different climatic adaptabilities for both genotypes. The exis-tance of a border zone between the areas of both subspecies, where were found a high frequency of hybrids, is explained as follows: being each subspecies adapted to a special climatic zone, we may suppose a poor adaptation of either one in the border region, which is also a region of intermediate climatic conditions. Thus, the hybrids, having a combination of the parent qualities, will be best adapted to the transition zone. Thus, the hybrids will become heterotic and an equilibrium will be reached with all genotypes present in the population in the border region.

Modeling AIDS vaccines: the cellular level

This paper discuses current strategies for the development of AIDS vaccines wich allow immunzation to disturb the natural course of HIV at different detailed stages of its life cycle. Mathematical models describing the main biological phenomena (i.e. virus and vaccine induced T4 cell growth; virus and vaccine induced activation latently infected T4 cells; incremental changes immune response as infection progress; antibody dependent enhancement and neutralization of infection) and allowing for different vaccination strategies serve as a backgroud for computer simulations. The mathematical models reproduce updated information on the behavior of immune cells, antibody concentrations and free viruses. The results point to some controversial outcomes of an AIDS vaccine such as an early increase in virus concentration among vaccinated when compared to nonvaccinated individuals.

The risk of acquiring the new influenza A(H1N1) for Brazilian travelers to Chile, Argentina and the USA

We estimate the risk of acquiring the new influenza A(H1N1) for Brazilian travelers to Chile, Argentina and the USA. This is done by a mathematical model that quantifies the intensity of transmission of the new virus in those countries and the probability that one individual has of acquiring the influenza depending on the date of arrival and time spent in the area. The maximum estimated risk reached 7.5 cases per 10,000 visitors to Chile, 17 cases per 10,000 travelers to Argentina and 23 cases per 10,000 travelers to the USA. The estimated number of imported cases until 27 July is 57 ± 9 from Chile, 136 ± 27 from the USA and 301 ± 21 from Argentina, which are in accord with the official figures. Estimating the number of imported cases was particularly important for the moment of the disease introduction into this country, but it will certainly be important again as a tool to calculate the number of future imported cases from northern countries in our next inter-epidemic season, were imported cases can constitute again the majority of the new influenza burden to the Brazilian health services.

New tools for epidemiology: a space odyssey

Geographical Information Systems (GIS) facilitate access to epidemiological data through visualization and may be consulted for the development of mathematical models and analysis by spatial statistics. Variables such as land-cover, land-use, elevations, surface temperatures, rainfall etc. emanating from earth-observing satellites, complement GIS as this information allows the analysis of disease distribution based on environmental characteristics. The strength of this approach issues from the specific environmental requirements of those causative infectious agents, which depend on intermediate hosts for their transmission. The distribution of these diseases is restricted, both by the environmental requirements of their intermediate hosts/vectors and by the ambient temperature inside these hosts, which effectively govern the speed of maturation of the parasite. This paper discusses the current capabilities with regard to satellite data collection in terms of resolution (spatial, temporal and spectral) of the sensor instruments on board drawing attention to the utility of computer-based models of the Earth for epidemiological research. Virtual globes, available from Google and other commercial firms, are superior to conventional maps as they do not only show geographical and man-made features, but also allow instant import of data-sets of specific interest, e.g. environmental parameters, demographic information etc., from the Internet.

A public health risk assessment for yellow fever vaccination: a model exemplified by an outbreak in the state of São Paulo, Brazil

We propose a method to analyse the 2009 outbreak in the region of Botucatu in the state of São Paulo (SP), Brazil, when 28 yellow fever (YF) cases were confirmed, including 11 deaths. At the time of the outbreak, the Secretary of Health of the State of São Paulo vaccinated one million people, causing the death of five individuals, an unprecedented number of YF vaccine-induced fatalities. We apply a mathematical model described previously to optimise the proportion of people who should be vaccinated to minimise the total number of deaths. The model was used to calculate the optimum proportion that should be vaccinated in the remaining, vaccine-free regions of SP, considering the risk of vaccine-induced fatalities and the risk of YF outbreaks in these regions.

Dielectric constant obtained from TDR and volumetric moisture of soils in southern Brazil

Soil moisture is the property which most greatly influences the soil dielectric constant, which is also influenced by soil mineralogy. The aim of this study was to determine mathematical models for soil moisture and the dielectric constant (Ka) for a Hapludalf, two clayey Hapludox and a very clayey Hapludox and test the reliability of universal models, such as those proposed by Topp and Ledieu and their co-workers in the 80's, and specific models to estimate soil moisture with a TDR. Soil samples were collected from the 0 to 0.30 m layer, sieved through a mesh of 0.002 m diameter and packed in PVC cylinders with a 0.1 m diameter and 0.3 m height. Seven samples of each soil class were saturated by capillarity and a probe composed of two rods was inserted in each one of them. Moisture readings began with the saturated soil and concluded when the soil was near permanent wilting point. In each step, the samples were weighed on a precision scale to calculate volumetric moisture. Linear and polynomial models were adjusted for each soil class and for all soils together between soil moisture and the dielectric constant. Accuracy of the models was evaluated by the coefficient of determination, the standard error of estimate and the 1:1 line. The models proposed by Topp and Ledieu and their co-workers were not adequate for estimating the moisture in the soil classes studied. The adjusted linear and polynomial models for the entire set of data of the four soil classes did not have sufficient accuracy for estimating soil moisture. The greater the soil clay and Fe oxide content, the greater the dielectric constant of the medium for a given volumetric moisture. The specific models, θ = 0.40283 - 0.04231 Ka + 0.00194 Ka² - 0.000022 Ka³ (Hapludox) θ = 0.01971 + 0.02902 Ka - 0.00086 Ka² + 0.000012 Ka³ (Hapludox -PF), θ = 0.01692 - 0.00507 Ka (Hapludalf) and θ = 0.08471 + 0.01145 Ka (Hapludox-CA), show greater accuracy and reliability for estimating soil moisture in the soil classes studied.

Use of the swat model for hydro-sedimentologic simulation in a small rural watershed

Mathematical models have great potential to support land use planning, with the goal of improving water and land quality. Before using a model, however, the model must demonstrate that it can correctly simulate the hydrological and erosive processes of a given site. The SWAT model (Soil and Water Assessment Tool) was developed in the United States to evaluate the effects of conservation agriculture on hydrological processes and water quality at the watershed scale. This model was initially proposed for use without calibration, which would eliminate the need for measured hydro-sedimentologic data. In this study, the SWAT model was evaluated in a small rural watershed (1.19 km²) located on the basalt slopes of the state of Rio Grande do Sul in southern Brazil, where farmers have been using cover crops associated with minimum tillage to control soil erosion. Values simulated by the model were compared with measured hydro-sedimentological data. Results for surface and total runoff on a daily basis were considered unsatisfactory (Nash-Sutcliffe efficiency coefficient - NSE < 0.5). However simulation results on monthly and annual scales were significantly better. With regard to the erosion process, the simulated sediment yields for all years of the study were unsatisfactory in comparison with the observed values on a daily and monthly basis (NSE values < -6), and overestimated the annual sediment yield by more than 100 %.

Correction of resistance to penetration by pedofunctions and a reference soil water content

The soil penetration resistance is an important indicator of soil compaction and is strongly influenced by soil water content. The objective of this study was to develop mathematical models to normalize soil penetration resistance (SPR), using a reference value of gravimetric soil water content (U). For this purpose, SPR was determined with an impact penetrometer, in an experiment on a Dystroferric Red Latossol (Rhodic Eutrudox), at six levels of soil compaction, induced by mechanical chiseling and additional compaction by the traffic of a harvester (four, eight, 10, and 20 passes); in addition to a control treatment under no-tillage, without chiseling or additional compaction. To broaden the range of U values, SPR was evaluated in different periods. Undisturbed soil cores were sampled to quantify the soil bulk density (BD). Pedotransfer functions were generated correlating the values of U and BD to the SPR values. By these functions, the SPR was adequately corrected for all U and BD data ranges. The method requires only SPR and U as input variables in the models. However, different pedofunctions are needed according to the soil layer evaluated. After adjusting the pedotransfer functions, the differences in the soil compaction levels among the treatments, previously masked by variations of U, became detectable.

Seleção de variáveis em QSAR

The process of building mathematical models in quantitative structure-activity relationship (QSAR) studies is generally limited by the size of the dataset used to select variables from. For huge datasets, the task of selecting a given number of variables that produces the best linear model can be enormous, if not unfeasible. In this case, some methods can be used to separate good parameter combinations from the bad ones. In this paper three methodologies are analyzed: systematic search, genetic algorithm and chemometric methods. These methods have been exposed and discussed through practical examples.

Moisture sorption isotherms and isosteric heat determination in Chilean papaya (Vasconcellea pubescens)

The moisture sorption isotherms of Chilean papaya were determined at 5, 20, and 45 ºC, over a relative humidity range of 10-95%. The GAB, BET, Oswin, Halsey, Henderson, Smith, Caurie and Iglesias-Chirife models were applied to the sorption experimental data. The goodness of fit of the mathematical models was statistically evaluated by means of the determination coefficient, mean relative percentage deviation, sum square error, root-mean-square error, and chi-square values. The GAB, Oswin and Halsey models were found to be the most suitable for the description of the sorption data. The sorption heats calculated using the Clausius-Clapeyron equation were 57.35 and 59.98 kJ·mol-1, for adsorption and desorption isotherms, respectively.

Modelagem da cinética de biodegradação de paclobutrazol em dois solos do semiárido do nordeste brasileiro

Mathematical models can help to prevent high levels of toxic substances in soil or fruits of plants treated with pesticides and indicate that such substances should be systematically monitored. The aim of this research was to study the kinetics of paclobutrazol biodegradation by soil native bacteria using mathematical models. Three models were used to assess the kinetics of paclobutrazol biodegradation obtained experimentally. Excellent fits were obtained using dual kinetic and logistic models. The use of glycerol as additional carbon source increased the biodegradation of PBZ and consequently decreased the time required for a given PBZ initial concentration be halved.

Numerical evaluation of the modulus of longitudinal elasticity in structural round timber elements of the Eucalyptus genus

Currently, the standards that deal with the determination of the properties of rigidity and strength for structural round timber elements do not take in consideration in their calculations and mathematical models the influence of the existing irregularities in the geometry of these elements. This study has as objective to determine the effective value of the modulus of longitudinal elasticity for structural round timber pieces of the Eucalyptus citriodora genus by a technique of optimization allied to the Inverse Analysis Method, to the Finite Element Method and the Least Square Method.

Software for calculation of reservoir active capacity with the sequent-peak algorithm

It is presented a software developed with Delphi programming language to compute the reservoir's annual regulated active storage, based on the sequent-peak algorithm. Mathematical models used for that purpose generally require extended hydrological series. Usually, the analysis of those series is performed with spreadsheets or graphical representations. Based on that, it was developed a software for calculation of reservoir active capacity. An example calculation is shown by 30-years (from 1977 to 2009) monthly mean flow historical data, from Corrente River, located at São Francisco River Basin, Brazil. As an additional tool, an interface was developed to manage water resources, helping to manipulate data and to point out information that it would be of interest to the user. Moreover, with that interface irrigation districts where water consumption is higher can be analyzed as a function of specific seasonal water demands situations. From a practical application, it is possible to conclude that the program provides the calculation originally proposed. It was designed to keep information organized and retrievable at any time, and to show simulation on seasonal water demands throughout the year, contributing with the elements of study concerning reservoir projects. This program, with its functionality, is an important tool for decision making in the water resources management.