Spatial statistical analysis and selection of genotypes in plant breeding

The objective of this study was to evaluate the efficiency of spatial statistical analysis in the selection of genotypes in a plant breeding program and, particularly, to demonstrate the benefits of the approach when experimental observations are not spatially independent. The basic material of this study was a yield trial of soybean lines, with five check varieties (of fixed effect) and 110 test lines (of random effects), in an augmented block design. The spatial analysis used a random field linear model (RFML), with a covariance function estimated from the residuals of the analysis considering independent errors. Results showed a residual autocorrelation of significant magnitude and extension (range), which allowed a better discrimination among genotypes (increase of the power of statistical tests, reduction in the standard errors of estimates and predictors, and a greater amplitude of predictor values) when the spatial analysis was applied. Furthermore, the spatial analysis led to a different ranking of the genetic materials, in comparison with the non-spatial analysis, and a selection less influenced by local variation effects was obtained.

Handling random errors and biases in methods used for short-term dietary assessment

Epidemiological studies have shown the effect of diet on the incidence of chronic diseases; however, proper planning, designing, and statistical modeling are necessary to obtain precise and accurate food consumption data. Evaluation methods used for short-term assessment of food consumption of a population, such as tracking of food intake over 24h or food diaries, can be affected by random errors or biases inherent to the method. Statistical modeling is used to handle random errors, whereas proper designing and sampling are essential for controlling biases. The present study aimed to analyze potential biases and random errors and determine how they affect the results. We also aimed to identify ways to prevent them and/or to use statistical approaches in epidemiological studies involving dietary assessments.

Population Dynamics Modeling of Arapaima gigas

Pirarucu (Arapaima gigas) has been of the most important natural fishing resources of the Amazon region. Due to its economic importance, and the necessity to preserve the species hand, field research concerning the habits and behavior of the pirarucu has been increasing for the last 20 years. The aim of this paper is to present a mathematical model for the pirarucu population dynamics considering the species peculiarities, particularly the male parental care over the offspring. The solution of the dynamical systems indicates three possible equilibrium points for the population. The first corresponds to extinction; the third corresponds to a stable population close to the environmental carrying capacity. The second corresponds to an unstable equilibrium located between extinction and full use of the carrying capacity. It is shown that lack of males’ parental care closes the gap between the point corresponding to the unstable equilibrium and the point of stable non-trivial equilibrium. If guarding failure reaches a critical point the two points coincide and the population tends irreversibly to extinction. If some event tends to destabilize the population equilibrium, as for instance inadequate parental care, the model responds in such a way as to restore the trajectory towards the stable equilibrium point avoiding the route to extinction. The parameters introduced to solve the system of equations are partially derived from limited but reliable field data collected at the Mamirauá Sustainable Development Reserve (MSDR) in the Brazilian Amazonian Region.

Uncertainties regarding dengue modeling in Rio de Janeiro, Brazil

Dengue fever is currently the most important arthropod-borne viral disease in Brazil. Mathematical modeling of disease dynamics is a very useful tool for the evaluation of control measures. To be used in decision-making, however, a mathematical model must be carefully parameterized and validated with epidemiological and entomological data. In this work, we developed a simple dengue model to answer three questions: (i) which parameters are worth pursuing in the field in order to develop a dengue transmission model for Brazilian cities; (ii) how vector density spatial heterogeneity influences control efforts; (iii) with a degree of uncertainty, what is the invasion potential of dengue virus type 4 (DEN-4) in Rio de Janeiro city. Our model consists of an expression for the basic reproductive number (R0) that incorporates vector density spatial heterogeneity. To deal with the uncertainty regarding parameter values, we parameterized the model using a priori probability density functions covering a range of plausible values for each parameter. Using the Latin Hypercube Sampling procedure, values for the parameters were generated. We conclude that, even in the presence of vector spatial heterogeneity, the two most important entomological parameters to be estimated in the field are the mortality rate and the extrinsic incubation period. The spatial heterogeneity of the vector population increases the risk of epidemics and makes the control strategies more complex. At last, we conclude that Rio de Janeiro is at risk of a DEN-4 invasion. Finally, we stress the point that epidemiologists, mathematicians, and entomologists need to interact more to find better approaches to the measuring and interpretation of the transmission dynamics of arthropod-borne diseases.

Trypanosoma cruzi strains isolated from human, vector, and animal reservoir in the same endemic region in Mexico and typed as T. cruzi I, discrete typing unit 1 exhibit considerable biological diversity

In this study, three strains of Trypanosoma cruzi were isolated at the same time and in the same endemic region in Mexico from a human patient with chronic chagasic cardiomyopathy (RyC-H); vector (Triatoma barberi) (RyC-V); and rodent reservoir (Peromyscus peromyscus) (RyC-R). The three strains were characterized by multilocus enzyme electrophoresis, random amplified polymorphic DNA, and by pathological profiles in experimental animals (biodemes). Based on the analysis of genetic markers the three parasite strains were typed as belonging to T. cruzi I major group, discrete typing unit 1. The pathological profile of RyC-H and RyC-V strains indicated medium virulence and low mortality and, accordingly, the strains should be considered as belonging to biodeme Type III. On the other hand, the parasites from RyC-R strain induced more severe inflammatory processes and high mortality (> 40%) and were considered as belonging to biodeme Type II. The relationship between genotypes and biological characteristics in T. cruzi strains is still debated and not clearly understood. An expert committee recommended in 1999 that Biodeme Type III would correspond to T. cruzi I group, whereas Biodeme Type II, to T. cruzi II group. Our findings suggest that, at least for Mexican isolates, this correlation does not stand and that biological characteristics such as pathogenicity and virulence could be determined by factors different from those identified in the genotypic characterization

Identification of a genetic marker associated with the resistance to Schistosoma mansoni infection using random amplified polymorphic DNA analysis

In schistosomiasis, the host/parasite interaction remains not completely understood. Many questions related to the susceptibility of snails to infection by respective trematode still remain unanswered. The control of schistosomiasis requires a good understanding of the host/parasite association. In this work, the susceptibility/resistance to Schistosoma mansoni infection within Biomphalaria alexandrina snails were studied starting one month post infection and continuing thereafter weekly up to 10 weeks after miracidia exposure. Genetic variations between susceptible and resistant strains to Schistosoma infection within B. alexandrina snails using random amplified polymorphic DNA analysis technique were also carried out. The results showed that 39.8% of the examined field snails were resistant, while 60.2% of these snails showed high infection rates.In the resistant genotype snails, OPA-02 primer produced a major low molecular weight marker 430 bp. Among the two snail strains there were interpopulational variations, while the individual specimens from the same snail strain, either susceptible or resistant, record semi-identical genetic bands. Also, the resistant character was ascendant in contrast to a decline in the susceptibility of snails from one generation to the next.

Genetic diversity of Leishmania infantum field populations from Brazil

Leishmania infantum (syn. Leishmania chagasi) is the etiological agent of visceral leishmaniasis (VL) in Brazil. The epidemiology of VL is poorly understood. Therefore, a more detailed molecular characterization at an intraspecific level is certainly needed. Herein, three independent molecular methods, multilocus microsatellite typing (MLMT), random amplification of polymorphic DNA (RAPD) and simple sequence repeats-polymerase chain reaction (SSR-PCR), were used to evaluate the genetic diversity of 53 L. infantum isolates from five different endemic areas in Brazil. Population structures were inferred by distance-based and Bayesian-based approaches. Eighteen very similar genotypes were detected by MLMT, most of them differed in only one locus and no correlation was found between MLMT profiles, geographical origin or the estimated population structure. However, complex profiles composed of 182 bands obtained by both RAPD and SSR-PCR assays gave different results. Unweighted pair group method with arithmetic mean trees built from these data revealed a high degree of homogeneity within isolates of L. infantum. Interestingly, despite this genetic homogeneity, most of the isolates clustered according to their geographical origin.

Fuzzy logic applied to the modeling of water dynamics in an Oxisol in northeastern Brazil

Modeling of water movement in non-saturated soil usually requires a large number of parameters and variables, such as initial soil water content, saturated water content and saturated hydraulic conductivity, which can be assessed relatively easily. Dimensional flow of water in the soil is usually modeled by a nonlinear partial differential equation, known as the Richards equation. Since this equation cannot be solved analytically in certain cases, one way to approach its solution is by numerical algorithms. The success of numerical models in describing the dynamics of water in the soil is closely related to the accuracy with which the water-physical parameters are determined. That has been a big challenge in the use of numerical models because these parameters are generally difficult to determine since they present great spatial variability in the soil. Therefore, it is necessary to develop and use methods that properly incorporate the uncertainties inherent to water displacement in soils. In this paper, a model based on fuzzy logic is used as an alternative to describe water flow in the vadose zone. This fuzzy model was developed to simulate the displacement of water in a non-vegetated crop soil during the period called the emergency phase. The principle of this model consists of a Mamdani fuzzy rule-based system in which the rules are based on the moisture content of adjacent soil layers. The performances of the results modeled by the fuzzy system were evaluated by the evolution of moisture profiles over time as compared to those obtained in the field. The results obtained through use of the fuzzy model provided satisfactory reproduction of soil moisture profiles.

MODELING OF SOIL LOAD-BEARING CAPACITY AS A FUNCTION OF SOIL MECHANICAL RESISTANCE TO PENETRATION

Estimation of soil load-bearing capacity from mathematical models that relate preconsolidation pressure (σp) to mechanical resistance to penetration (PR) and gravimetric soil water content (U) is important for defining strategies to prevent compaction of agricultural soils. Our objective was therefore to model the σp and compression index (CI) according to the PR (with an impact penetrometer in the field and a static penetrometer inserted at a constant rate in the laboratory) and U in a Rhodic Eutrudox. The experiment consisted of six treatments: no-tillage system (NT); NT with chiseling; and NT with additional compaction by combine traffic (passing 4, 8, 10, and 20 times). Soil bulk density, total porosity, PR (in field and laboratory measurements), U, σp, and CI values were determined in the 5.5-10.5 cm and 13.5-18.5 cm layers. Preconsolidation pressure (σp) and CI were modeled according to PR in different U. The σp increased and the CI decreased linearly with increases in the PR values. The correlations between σp and PR and PR and CI are influenced by U. From these correlations, the soil load-bearing capacity and compaction susceptibility can be estimated by PR readings evaluated in different U.

Phase diagram and complex patterns in the modeling of the bromate-oxalic acid-Ce-acetone oscillating reaction

Simulations have been carried out on the bromate - oxalic acid - Ce(IV) - acetone oscillating reaction, under flow conditions, using Field and Boyd's model (J. Phys. Chem. 1985, 89, 3707). Many different complex dynamic behaviors were found, including simple periodic oscillations, complex periodic oscillations, quasiperiodicity and chaos. Some of these complex oscillations can be understood as belonging to a Farey sequence. The many different behaviors were systematized in a phase diagram which shows that some regions of complex patterns were nested with one inside the other. The existence of almost all known dynamic behavior for this system allows the suggestion that it can be used as a model for some very complex phenomena that occur in biological systems.

MODELING A COMPUTER APPLICATION FOR MANAGEMENT OF MAINTENANCE ACTIVITIES OF UNPAVED ROADS

ABSTRACTThis study presents a contribution to the modeling of a computer application employing a method of serviceability performance for unpaved roads, aiming the management of maintenance/restoration activities of the primary surface layer. The proposed methodology consisted of field inspections during dry (April to September) and rainy (October to March) periods, during which objective evaluations were performed to survey of defects and their densities and degrees of severity. To aid the functional classification of analyzed road sections and the determination of the defect with major influence on the serviceability of these roads, the method of serviceability performance proposed by Silva (2009)was implemented in the Visual Basic for Applications (VBA) language in Microsoft Excel software. With the use of the computer application proposed it was possible to identify among the defects analyzed in field, through the index of serviceability of the sampling unit per defect type (ISUdef), which one had the greatest influence on determining the relative serviceability index per road section (IST). The results allow us to conclude that the computer application Road achieved satisfactory results, since the objective evaluation criteria applied to road sections denotes consistency regarding their serviceability.

MODELING AND SIMULATION OF RELIEF INFLUENCE ON EUCALYPTUS FORESTS: INTERACTION BETWEEN SOLAR IRRADIANCE AND PRODUCTIVITY

ABSTRACT This study aimed to verify the differences in radiation intensity as a function of distinct relief exposure surfaces and to quantify these effects on the leaf area index (LAI) and other variables expressing eucalyptus forest productivity for simulations in a process-based growth model. The study was carried out at two contrasting edaphoclimatic locations in the Rio Doce basin in Minas Gerais, Brazil. Two stands with 32-year-old plantations were used, allocating fixed plots in locations with northern and southern exposure surfaces. The meteorological data were obtained from two automated weather stations located near the study sites. Solar radiation was corrected for terrain inclination and exposure surfaces, as it is measured based on the plane, perpendicularly to the vertical location. The LAI values collected in the field were used. For the comparative simulations in productivity variation, the mechanistic 3PG model was used, considering the relief exposure surfaces. It was verified that during most of the year, the southern surfaces showed lower availability of incident solar radiation, resulting in up to 66% losses, compared to the same surface considered plane, probably related to its geographical location and higher declivity. Higher values were obtained for the plantings located on the northern surface for the variables LAI, volume and mean annual wood increase, with this tendency being repeated in the 3PG model simulations.

Modeling, optimizing and simulating robot calibration with accuracy improvement

This work describes techniques for modeling, optimizing and simulating calibration processes of robots using off-line programming. The identification of geometric parameters of the nominal kinematic model is optimized using techniques of numerical optimization of the mathematical model. The simulation of the actual robot and the measurement system is achieved by introducing random errors representing their physical behavior and its statistical repeatability. An evaluation of the corrected nominal kinematic model brings about a clear perception of the influence of distinct variables involved in the process for a suitable planning, and indicates a considerable accuracy improvement when the optimized model is compared to the non-optimized one.

Mathematical modeling of electro-rotation spectra of small particles in liquid solutions: Application to human erythrocyte aggregates

Electro-rotation can be used to determine the dielectric properties of cells, as well as to observe dynamic changes in both dielectric and morphological properties. Suspended biological cells and particles respond to alternating-field polarization by moving, deforming or rotating. While in linearly polarized alternating fields the particles are oriented along their axis of highest polarizability, in circularly polarized fields the axis of lowest polarizability aligns perpendicular to the plane of field rotation. Ellipsoidal models for cells are frequently applied, which include, beside sphere-shaped cells, also the limiting cases of rods and disks. Human erythrocyte cells, due to their particular shape, hardly resemble an ellipsoid. The additional effect of rouleaux formation with different numbers of aggregations suggests a model of circular cylinders of variable length. In the present study, the induced dipole moment of short cylinders was calculated and applied to rouleaux of human erythrocytes, which move freely in a suspending conductive medium under the effect of a rotating external field. Electro-rotation torque spectra are calculated for such aggregations of different length. Both the maximum rotation speeds and the peak frequencies of the torque are found to depend clearly on the size of the rouleaux. While the rotation speed grows with rouleaux length, the field frequency nup is lowest for the largest cell aggregations where the torque shows a maximum.